2013机械振动数值分析大作业.doc

研究生大作业共页第1页二O一二~二O一三学年第2学期《机械振动数值分析》课程作业日期:年月日课程编号:、质量矩阵以及初始条件分别为编制基于雅可比迭代法的特征值问题计算子程序(程序要求:该子程序的输入变量为刚度矩阵、质量矩阵和收敛准则,输出为固有频率和对应的固有振型);利用所编制的程序计算该系统的全部固有频率和固有振型,将求解出的固有频率与精确解进行对比,并绘制出各固有频率随迭代次数变化的曲线;编制基于中心差分法的瞬态响应分析子程序,计算第一个自由度在10s内的自由振动位移响应(程序要求:该子程序的输入变量为刚度矩阵、质量矩阵和计算时间,输出为位移响应)。(1到12号节点固定),。全部杆截面均为圆截面,截面直径为。各杆之间为铰接。钢的杨氏模量取为210GPa,,密度为7800kg/m3。编制基于逆迭代法的特征值问题计算子程序(程序要求:该子程序的输入变量为刚度矩阵、质量矩阵、收敛准则及所需要求解出的特征值个数,输出为固有频率和对应的固有振型);共页第2页编制基于Lanczos法的特征值问题计算子程序(程序要求:该子程序的输入变量为刚度矩阵、质量矩阵、收敛准则及所需要求解出的特征值个数,输出为固有频率和对应的固有振型。编制出的子空间迭代法程序必须调用雅可比迭代法进行特征值问题的计算);利用所编制的逆迭代法子程序和Lanczos法子程序完成图1所示模型的前10阶固有频率和固有振型的计算,绘制出固有振型图;假设初始时刻在桁架中心点处施加三角形脉冲激励,激励幅值为100N,脉冲宽度为50ms,考虑比例阻尼模型,取,,编制基于Newmark方法的瞬态响应分析程序,考虑两种不同的时间步长,分别计算载荷作用点处的位移响应,并进行比较,说明时间步长对计算