哈工大结构动力学大作业2012春
哈尔滨工业大学 结构力学II 第二套张金生 结构动力学-9
1 T 2 最大动能为 Ti max X i m X i i 2 1 T 最大势能为 U i max X i k X i 2 Ti max U i max 由能量守恒，有
mN
y N (t )
X T k X i i i2 X T mX i i
2
5.04 1 m 4 9.08 归一化 X 1.8 k 2.24 11.22
1
4.667 m 8.334 归一化 k 10.334 4.99 m 8.98 归一化 k 11.19
k
k
1.取自重引起的位移 y1 3mg / k y2 y1 2mg / k 5mg / k
mg
mg
y1
y3 y1 mg / k 6mg / k
3 X 1 5 精确解: 2 6 1 0.198k / m T X 1 k X 1 2 1 0.445 k / m 1 T X 1 mX 1 14k 0.2k / m 1 0.447 k / m 70m
X n D n 1 n 一直算到 X 与 X 成比例为止. X n 为基本振型. 这时下式成立 1 n 1 n n X DX X 2
1
X
基本频率由下式计算
令பைடு நூலகம்
X 2 mX D m ---动力矩阵 1 X DX ---标准特征值问题 2
j
1 2 j
( a b 2 ) j
例.求图示体系的正交阻尼矩阵 c 和阻尼比 3 . 已知: 1 2 0.05
m m
k k
3 2 1
哈工大结构动力学大作业2012春
结构动力学大作业对于如下结构，是研究质量块的质量变化和在简支梁上位置的变化对整个系统模态的影响。
1以上为一个简支梁结构。
集中质量块放于梁上，质量块距简支梁的左端点距离为L.将该简支梁简化为欧拉伯努利梁，并离散为N 个单元。
每个单元有两个节点，四个自由度。
单元的节点位移可表示为：]1122,,,e v v δθθ⎡=⎣则单元内一点的挠度可计作：带入边界条件：1332210)(x a x a x a a x v +++=01)0(a v x v ===3322102)(L a L a L a a v L x v +++===110d d a xv x ===θ2321232d d L a L a a xv Lx ++===θ10v a =[]1234N N N N N =建立了单元位移模式后，其动能势能均可用节点位移表示。
单元的动能为：00111()222l l T TT ke e e e e y E dx q N Ndxq q mq t ρρ∂===∂⎰⎰ 其中m 为单元质量阵，并有：lT m N Ndx ρ=⎰带入公式后积分可得：2222156225413224133541315622420133224l l l l l l l m l l ll ll ρ-⎡⎤⎢⎥-⎢⎥=⎢⎥-⎢⎥---⎣⎦单元势能可表示为2220011()()222T l lT Te pe e e e q y E EI dx EI N N dxq q Kq x ∂''''===∂⎰⎰其中K 为单元刚度矩阵，并有()lT K EI N N dx ''''=⎰2232212612664621261266264l l l l l l EI k l l l l lll -⎡⎤⎢⎥-⎢⎥=⎢⎥---⎢⎥-⎣⎦以上为单元类型矩阵，通过定义全局位移矩阵，可以得到系统刚度矩阵和系统质量矩11θ=a )2(1)(3211222θθ+--=Lv v L a )(1)(22122133θθ++-=Lv v L a 1232133222231)(θ⎪⎪⎭⎫ ⎝⎛+-+⎪⎪⎭⎫ ⎝⎛+-=L x L x x v L x L x x v 22232332223θ⎪⎪⎭⎫ ⎝⎛-+⎪⎪⎭⎫ ⎝⎛-+L x L x v L x L x 24231211)()()()()(θθx N v x N x N v x N x v +++=阵。
哈工程机械动力学大作业
机械动力学大作业含弹性摆杆的铰链四杆机构动力学仿真学号：院系名称：机电工程学院专业：机械工程学生姓名：本次进行设计和分析的对象为平面铰链四杆机构，在Adams的环境下，通过对四杆机构进行建模以及运动仿真，绘制出摆杆的相关曲线图。
为了形成有效的对比，先建立含有刚性摆杆的四杆机构，进行运动仿真，绘制出摆杆的相关曲线。
再建立含有柔性摆杆的铰链四杆机构，所有参数设置均和刚性摆杆一样。
考虑到弹性摇杆可能发生较大的形变，不利于观测，绘制摇杆运动曲线时选择摇杆的质心作为参考点。
在Adams中主要有三种方法创建柔性构件，第一种是将刚性构件离散化后采用柔性梁连接；第二种是直接将刚体替换为柔性体；第三种是运用有限元分析的方法建立柔性构件。
本次建模，主要采用前两种方法建立柔性摆杆。
运用有限元建立柔性构件，等以后再进行深入研究。
同时两种方法建立的柔性杆可以形成对比。
通过本次设计，主要学习了Adams 软件建模以及运动仿真、图形处理、刚柔混合建模的操作方法，对自己也是一个很大锻炼和提升。
设计的为平面曲柄摇杆机构。
相关参数如：曲柄长L=200mm，宽W=60mm,高D=30mm；连杆长L=427mm,宽W=30mm,高D=20mm；摇杆长L=403mm，宽W=40mm，高D=20mm；机架长L=600mm，宽W=40mm，高D=20mm；曲柄角速度为40deg/sec。
经过验证，最短杆长度加上最长杆长度小于中间两根杆的长度之和，满足曲柄存在的条件，且最长杆为机架，故为曲柄摇杆机构。
一、建模过程1、建立四个标记点，这四个点依次连接就可以确定一个铰链四杆机构。
2、建立四根杆的模型3、在杆件之间添加转动副4、选择最长杆为机架并固定5、给曲柄添加驱动，使曲柄角速度为40deg/sec。
6、使模型的显示方式为实体显示7、进行运动仿真，主要通过时间和步长来控制仿真运动的快慢。
8、载入动画，进行图像处理，绘制曲线图。
9、将刚性摇杆换成柔性的摇杆并添加转动副和驱动10、进行运动仿真11、载入动画，进行图像处理，绘制曲线。
最新哈工大春季钢结构基本原理及设计试题参考答案
一、填空题：(每空1分，共22分)1、《建筑钢结构焊接技术规程》（JGJ－2002）推荐使用碳当量(或C E) 来衡量低合金钢的焊接性能。
2、硫、磷、氮、氧和氢均为有害元素，其中磷和氮易引起钢材的低温冷脆。
3、影响结构疲劳寿命的最主要因素是构造状态、循环荷载和循环次数。
4、钢材的机械性能指标为屈服强度、抗拉强度、伸长率、冷弯性能、Z向收缩率和冲击韧性。
5、主要焊接结构不能采用Q235 A 级钢，因含碳量不作交货条件，无法保证含碳量。
6、将一块Q235B级钢板与Q345B级钢板通过焊接进行连接时，宜选择E43型焊条。
7、«钢结构设计规范»（GB50017－2003）规定：质量属于三级的焊缝，焊缝的抗拉设计强度等于母材抗拉设计强度的0.85 倍。
8、单轴对称的T形截面轴心受压构件，绕对称轴失稳时易发生弯扭失稳，绕非对称轴时易发生弯曲失稳。
9、轴心受压构件失稳时可能绕两主轴屈曲，设计时应遵循等稳定原则，如进行梯形钢屋架设计时屋架端斜杆应采用不等边角钢长肢相连的截面型式。
10、在轴心受压构件中，确定箱形截面板件满足局部稳定的宽（高）厚比限值的原则是构件应力达到屈服前其板件不发生局部屈曲（或局部屈曲临界应力不低于屈服应力，或不先于屈服），确定工字形截面确定板件宽（高）厚比限值的原则是构件整体屈曲前其板件不发生局部屈曲（或局部屈曲临界应力不低于整体屈曲临界应力或等稳定或不先于整体失稳）。
11、荷载作用点的位置对梁的整体稳定有影响，相对于荷载作用于工字形截面简支梁受拉翼缘，当荷载作用于梁的受压翼缘时，其梁的整体稳定性将降低。
12、某梯形钢屋架，下弦支座处应设置刚性系杆，该系杆需要按受压杆设计。
13、 某工字形组合截面简支梁，若腹板的高厚比为100，应设置 横向 加劲肋，若腹板高厚比为210，应设置 纵向 加劲肋。
二、单项选择题(每题2分，共14分)1、 最易产生脆性破坏的应力状态是 B 。
哈尔滨工业大学01-03年春秋结构力学试卷及答案精品文档29页
哈工大 2001 年春季学期结构力学试卷(请考生注意：本试卷共5页)一．是非题（将判断结果填入括弧：以O 表示正确，X 表示错误）（本大题分4小题，共11分）1 . (本小题 3分)图示结构中DE 杆的轴力F NDE =F P /3。
（ ）.2 . (本小题 4分)用力法解超静定结构时，只能采用多余约束力作为基本未知量。
（ ） 3 . (本小题 2分)力矩分配中的传递系数等于传递弯矩与分配弯矩之比，它与外因无关。
（） 4 . (本小题 2分)用位移法解超静定结构时，基本结构超静定次数一定比原结构高。
（ ）二．选择题（将选中答案的字母填入括弧内）（本大题分5小题，共21分）1 （本小题6分）图示结构EI=常数，截面A 右侧的弯矩为：（ ） A ．2/M ； B ．M ； C ．0； D. )2/(EI M 。
2. （本小题4分）图示桁架下弦承载，下面画出的杆件内力影响线，此杆件是：（ ） Ａ.ch ; B.ci; C.dj; D .cj .3. (本小题 4分)图a 结构的最后弯矩图为：A. 图b;B. 图c;C. 图d;D.都不对。
（ ）( a) (b) (c) (d)4. (本小题 4分)用图乘法求位移的必要条件之一是：A.单位荷载下的弯矩图为一直线；2=1/4 M /8 M /2B.结构可分为等截面直杆段；C.所有杆件EI 为常数且相同；D.结构必须是静定的。
( ) 5. （本小题3分）图示梁A 点的竖向位移为（向下为正）：( )Ａ．F P l 3/(24EI ); B . F P l 3/(!6EI ); C . 5F P l 3/(96EI ); D. 5F P l 3/(48EI ).三（本大题 5分）对图示体系进行几何组成分析。
四（本大题 9分）图示结构B 支座下沉4 mm ，各杆EI=2.0×105 kN ·m 2，用力法计算并作M 图。
五（本大题 11分） 用力矩分配法计算图示结构，并作M 图。
结构动力学大作业
结构动力学作业姓名：学号：目录1.力插值法 (1)1.1分段常数插值法 (1)1.2分段线性插值法 (4)2.加速度插值法 (7)2.1常加速度法 (7)2.2线加速度法 (9)附录 (12)分段常数插值法源程序 (12)分段线性插值法源程序 (12)常加速度法源程序 (13)线加速度法源程序 (13)1.力插值法力插值法对结构的外荷载进行插值，分为分段常数插值法和分段线性插值法，这两种方法均适用于线性结构的动力反应计算。
1.1分段常数插值法图1-1为一个单自由度无阻尼系统，结构的刚度为k ，质量为m ，位移为y (t )，施加的外力为P (t )。
图1-2为矩形脉冲荷载的示意图，图中t d 表示作用的时间，P 0表示脉冲荷载的大小。
图1-1 单自由度无阻尼系统示意图图1-2 矩形脉冲荷载示意图对于一个满足静止初始条件的无阻尼单自由度体系来说，当施加一个t d 时间的矩形脉冲荷载，此时结构在t d 时间内的位移反应可以用杜哈梅积分得到：0()sin ()2 (1cos )(1cos ) (0)tst st d P y t t d m ty t y t t Tωττωπω=-=-=-≤≤⎰(1-1)如果结构本身有初始的位移和速度，那么叠加上结构自由振动的部分，结构的位移反应为：002()cos sin (1cos) (0)st d yty t y t t y t t Tπωωω=++-≤≤ (1-2)图1-3分段常数插值法微段示意图对于施加于结构任意大小的力，将其划分为Δt 的微段，每一段的荷载都为一个常数(每段相当于一个矩形的脉冲荷载)，如图1-3所示，则将每一段的位移和速度写成增量的形式为：1cos t sin t (1cos t)i ii i yP y y kωωωω+=∆+∆+-∆ (1-3)i+1/sin t cos t sin t i ii yP yy kωωωωω=-∆+∆+∆ (1-4)程序流程图如下图1-4分段常数插值法流程图根据流程图可以编写相应的算法，利用MATLAB 进行编程，程序源代码见附录。
哈工大研究生动力学试卷08年秋(A)-final
结构动力学试卷
（2008秋，研究生学位课程）
姓名： 学号： 得分：
第一部分：非主观性试题（共20分）
回答下面问题。
（共20分）
（1） 何谓结构动力特性？（4分）
（2） 简述结构线性振动与非线性振动的区别；（4分） （3） 简述线加速度法计算精度受哪些因素影响；（ 6分） （4） 试确定下面图（a ）,（b ）体系的动力自由度；(6分)
第二部分：主观性试题（共80分）
一、(15分)试列出图示体系的运动方程，图中杆件CD 为匀质刚性杆，质量分布集度为
m ，杆件AB 为无重刚杆，均布弹性支撑的刚度系数为 a k
k a
, 弹簧BC 刚度系数为k ， 阻尼器阻尼系数为c ，如图所示。
（a ） 平面刚架
平面桁架
题一图
二、(20分)试求题二图所示体系：（1）自振周期；（2）阻尼比；（3）稳态振幅。
质点质量为m , 阻尼系数为C ，各杆为无重杆，各杆尺寸与刚度如图所示,ω=2θ。
二、 (20分)试求题三图所示体系的稳态振幅。
0sin t
θ
()sin p t p t
θ=题三图 ()sin p t p t
θ=题三图
四、（15分）对题四图所示结构，试：（1）列出振型方程，并写出边界条件；（2）求稳态振幅。
不计杆BC 段质量。
五、（10分）试求题五图所示结构各构件的一致刚度矩阵与一致质量矩阵，并集成结
构的刚度与质量矩阵。
题五图
,
m
P P 0=(1,0,2)
,0,(1,0,3)(1,
0,4)(1,0,5)。
2012年哈工大考研结构力学试题一
填空∙路基的干湿类型划分为4类，即（干燥，中湿，潮湿，过湿）。
∙路基路面整体结构应有足够的承载能力，体现在（强度）和（刚度）两方面。
∙路基的最小填土高度，一般应保证路基处于（干燥）或（中湿）状态。
∙为避免挖方边坡零星土块下落堆积，保护边沟不致阻塞，可在挖方坡脚处设置（碎落台）。
∙在重力式挡土墙中，为了减少圬工砌体因硬化收缩和温度变化作用而产生的裂缝，须设置（伸缩缝）。
∙衡量土基压实程度的指标是（干容重）。
∙普通水泥混凝土路面的抗滑标准以（构造深度）为指标表示。
∙对于柔性路面，当采用两层沥青混凝土面层时，为增加层与层间的结合，应在中间设置（粘层）。
∙公路自然区划中二级区划得划分主要以（潮湿系数）为依据。
∙挡土墙按（极限状态设计的分项系数法）进行设计。
∙为防止挡土墙（不均匀沉降）引起强身开裂应设置沉降缝。
∙在排水纵坡陡于10%水头高差大于1米的陡坡地段，可以设置（跌水）和（急流槽）排除。
∙（跌水）和（急水槽）是地面排水沟渠的两种特殊方式，通常设在都坡处。
∙常用的坡面防护设施有（植物防护）和（工程防护）。
∙路基的典型横断面型式可分为（路堑）（路堤）和（填挖结合）等三种类型。
∙土基回弹模量可以采用（查表法）（现场实测法）室内试验法，换算法等方法获得。
∙《公路沥青路面设计规范》规定，路面设计以（双轮组单轴轴载100KN）为标准轴载，以（BZZ-100）表示。
∙石灰土中，石灰质量应符合（Ⅲ）级以上标准。
∙一班沥青混合料具有较高的（抗压）强度。
∙水泥混凝土路面以（抗弯拉强度）作为设计控制指标，用（劈裂实验）试验方法确定。
∙路基的填筑方法可分为（分层平铺），（竖向填筑）。
∙河滩路堤在水位变化时，除了受外力及自重外，还要受到（浮力），（动水压力）作用。
∙《公路沥青路面设计规范》规定，路面设计应采用（双圆垂直均布荷载）作用下的（多层弹性层状体系）理论，以（设计弯沉值）为路面整体刚度的设计指标。
∙沥青混合料的沥青最佳用量，通常以（马歇尔）试验来确定。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
结构动力学大作业对于如下结构，是研究质量块的质量变化和在简支梁上位置的变化对整个系统模态的影响。
1以上为一个简支梁结构。
集中质量块放于梁上，质量块距简支梁的左端点距离为L.将该简支梁简化为欧拉伯努利梁，并离散为N 个单元。
每个单元有两个节点，四个自由度。
单元的节点位移可表示为：]1122,,,e v v δθθ⎡=⎣则单元内一点的挠度可计作：带入边界条件：1332210)(x a x a x a a x v +++=01)0(a v x v ===3322102)(L a L a L a a v L x v +++===110d d a x vx ===θ2321232d d L a L a a xv Lx ++===θ10v a =[]1234N N N N N =建立了单元位移模式后，其动能势能均可用节点位移表示。
单元的动能为：00111()222l l T T Tke e e e e y E dx qN Ndxq q mq t ρρ∂===∂⎰⎰ 其中m 为单元质量阵，并有：lT m N Ndx ρ=⎰带入公式后积分可得：2222156225413224133541315622420133224l l l l l l l m l l l l l l ρ-⎡⎤⎢⎥-⎢⎥=⎢⎥-⎢⎥---⎣⎦单元势能可表示为2220011()()222Tl lTT e pe e e e q y E EI dx EI N N dxq q Kq x ∂''''===∂⎰⎰ 其中K 为单元刚度矩阵，并有()lT K EI N N dx ''''=⎰2232212612664621261266264l l l l l l EI k l l l l l l l -⎡⎤⎢⎥-⎢⎥=⎢⎥---⎢⎥-⎣⎦以上为单元类型矩阵，通过定义全局位移矩阵，可以得到系统刚度矩阵和系统质量矩11θ=a )2(1)(3211222θθ+--=Lv v L a )(1)(22122133θθ++-=Lv v L a 1232133222231)(θ⎪⎪⎭⎫⎝⎛+-+⎪⎪⎭⎫ ⎝⎛+-=L x L x x v L x L x x v 22232332223θ⎪⎪⎭⎫ ⎝⎛-+⎪⎪⎭⎫ ⎝⎛-+L x L x v L x L x 24231211)()()()()(θθx N v x N x N v x N x v +++=阵。
当集中质量块加到简支梁上时，可以认为系统的刚度矩阵不变，而质量矩阵在 有基础上有所增加。
当集中质量块位于第i 个单元内，距单元左节点的位移为0x ，则其动能可表示为：000()()2T Te e m V N x N x δδ=则质量块附加到总体刚度矩阵的子矩阵尾为：000()()T m N x N x例如当质量为0m 的质量块位于单元的右节点上时，其质量阵可表示为00000000000100000m m ⎡⎤⎢⎥⎢⎥'=⎢⎥⎢⎥⎣⎦根据以上理论，编写matlab 程序。
（相关参数：密度7860KG/M3， 长度L=1m ，截面尺寸0.02m*0.02m ）以下为matlab 源代码： 主程序clear; clc;Beam_InputData541; %输入相关参数，划分了40个单 %元k=zeros(No_nel*No_dof,No_nel*No_dof); %初始化单元刚度阵和质量阵 m=zeros(No_nel*No_dof,No_nel*No_dof); kk=zeros(Sys_dof,Sys_dof);mm=zeros(Sys_dof,Sys_dof); %初始化总体质量阵和刚度阵index=zeros(No_nel*No_dof,1); foriel=1:No_el nd(1)=iel; nd(2)=iel+1; leng=0.025;[k,m]=BeamElement11(prop,leng); %调用子函数形成单元刚度矩阵 %质量矩阵index=femEldof(nd,No_nel,No_dof);kk=femAssemble1(kk,k,index);mm=femAssemble1(mm,m,index); %调用子函数形成总体刚度矩阵 end %和质量矩阵Mzhiliangkuai=1/5*prop(3)*prop(6)*1; %将集中质量块加进总体质量矩Ms=zeros(82,82); %之中Ms(41,41)=Mzhiliangkuai; %此为质量块为于梁中央时的处mm=Ms+mm; 理方法kk(1,:)=[]; %引入边界条件，划去相应的kk(80,:)=[]; %行和列kk(:,1)=[];kk(:,80)=[];mm(1,:)=[];mm(80,:)=[];mm(:,1)=[];mm(:,80)=[];[V,D]=eig(kk,mm); %求解特征值和特征向量%或使用雅克比迭代求解特征值与特征向量，以下为代码%function[aa,v]=jac(a)%此程序用jacobi方法求实对称矩阵的全部特征值和特征向量%输入x：n x n矩阵%----------------------------------------------------%% n=length(a);% aa=zeros(1,n);% u=zeros(n,n);% l=0;% v是n阶单位矩阵% v=eye(n,n);% while (1)% l=l+1;% fm=0.0;% for i=1:n% for j=1:n% if i==j% continue;% else% d=abs(a(i,j));% if d>fm% fm=d;% p=i;% q=j;% end% end% end% end% if fm<0.0001% break;% end% x=2*a(p,q)*(sign(a(p,p)-a(q,q)));% y=abs(a(p,p)-a(q,q));% if y~=0% cn=sqrt((1+y/sqrt(x^2+y^2))/2.0);% else% cn=1/sqrt(2);% end% if y~=0% sn=x/sqrt(x^2+y^2)/(2*cn);% elseif y==0 & a(p,q)>0% sn=1/sqrt(2);% elseif y==0 & a(p,q)<0% sn=-1/sqrt(2);% end% %更新a% a_t=a;% for i=1:n% if i~=p & i~=q% a(i,p)=a_t(i,p)*cn+a_t(i,q)*sn;% a(i,q)=-a_t(i,p)*sn+a_t(i,q)*cn;% a(p,i)=a(i,p);% a(q,i)=a(i,q);% end% end% a(p,p)=a_t(p,p)*cn*cn+2*a_t(p,q)*cn*sn+a_t(q,q)*sn*sn; % a(q,q)=a_t(p,p)*sn*sn-2*a_t(p,q)*sn*cn+a_t(q,q)*cn*cn; % a(p,q)=0.0;% a(q,p)=a(p,q);% %更新v% v_t=v;% for i=1:n% v(i,p)=v_t(i,p)*cn+v_t(i,q)*sn;% v(i,q)=-v_t(i,p)*sn+v_t(i,q)*cn;% end% end% for i=1:n% aa(1,i)=a(i,i);% end%end%for i=1:10;% w(i,1)=sqrt(D(i,i))/(2*pi);%end%E=prop(1); Iz=prop(8); rho=prop(3)*prop(6); L=1;%i=(1:10)';%omega2=i.*i*pi^2*sqrt(E*Iz/(rho*L^4));%omega3=omega2/(2*pi);%for i=1:39;% x1(i)=i*0.025;% y1(i)=V(2*i,3);%end%plot(x1,y1);各子函数：Beam_InputData541No_el=40; number of elementsNo_nel=2; No_dof=2;No_node=(No_nel-1)*No_el+1;Sys_dof=No_node*No_dof;prop(1)=2.1e11; % elastic modulusprop(3)=7860; prop(4)=0.02;prop(5)=0.02; prop(6)=0.02*0.02;prop(7)=0.02*0.02^3/12; prop(8)=0.02*0.02^3/12;BeamElement11function [k,m]=BeamElement11(prop,leng);%,Opt_mass)E=prop(1);u=prop(2);rho=prop(3);A=prop(6); Iz=prop(8);G=E/(2*(1+u));c=E*Iz/(leng^3);k0=[12 6*leng -12 6*leng;6*leng 4*leng^2 -6*leng 2*leng^2;-12 -6*leng 12 -6*leng;6*leng 2*leng^2 -6*leng 4*leng^2];k=c*k0;mass=rho*A*leng;m0=[156 22*leng 54 -13*leng;22*leng 4*leng^2 13*leng -3*leng^2;54 13*leng 156 -22*leng;-13*leng -3*leng^2 -22*leng 4*leng^2];m=mass/420*m0;femAssemble1function [kk]=femAssemble1(kk,k,index)eldof = length(index);for i=1:eldofii=index(i);for j=1:eldofjj=index(j);kk(ii,jj)=kk(ii,jj)+k(i,j);endendfunction [index]=femEldof(nd,No_nel,No_dof)k=0;for i=1:No_nelstart = (nd(i)-1)*No_dof;for j=1:No_dofk=k+1;index(k)=start+j;endend首先我分析了无集中质量块时梁的模态特性，提取了前十阶固有频率和前四阶振型。
然后分析了当集中质量块分别为简支梁质量1/40,1/20.1/10,1/5,1/2,1 倍和2倍时这八种工况。