1、a1 1 1-1 1-1 利用动态静力法进行动力分析利用动态静力法进行动力分析 一、思路一、思路动静法：动静法：第一章第一章 单自由度单自由度机械系统的动力学分析机械系统的动力学分析根据达朗贝尔原理将惯性力计入静力平衡根据达朗贝尔原理将惯性力计入静力平衡 方程，来求解未知力（如原动件上施加的力、方程，来求解未知力（如原动件上施加的力、约束反力等）。约束反力等）。用静力平衡方程解决动力学问题用静力平衡方程解决动力学问题基本方程为：基本方程为：JMmaFa2 2二、典型实例二、典型实例例例1 1：已知：已知：求：角加速度求：角加速度解：利用动静法拆开机构解：利用动静法拆开机构 轮轮1 1：有反力：
2、有反力R R，惯性力矩惯性力矩 ,M,M1 1 轮轮2 2：有反力：有反力R R，惯性力矩，惯性力矩 ,M,M2 2则有方程：则有方程：得得:12!212,z zJ JM M111J22J1212121212(/)(/)MMzzJJzz结论：结论：1 1、加惯性力（力矩）、加惯性力（力矩）核心核心2 2、约束反力、约束反力 纽带纽带3 3、一个构件列一个受力平衡方程、一个构件列一个受力平衡方程基础基础01111 JRrM02222 JMRra3 3例例2 2：已知从动件推程方程已知从动件推程方程:求：凸轮角加速度求：凸轮角加速度解：忽略摩擦时反力解：忽略摩擦时反力R R，沿法线方向，沿法线方向
3、 凸轮：有反力凸轮：有反力R R ，惯性力矩，惯性力矩，M M1 1 推杆：有反力推杆：有反力R R，惯性力矩，惯性力矩，F F2 2则有方程：则有方程：得：得：1212,AShJmMF 111022sin()0cos0AMJRrSRFm S121212(/)(/)AMF hJm h00/vhtgrSrS结论：结论：例例1 1的角加速度是用传动比的角加速度是用传动比例例2 2的角加速度是用推杆位移方程的角加速度是用推杆位移方程a4 4例例3 3：已知：已知：求：求：建立运动方程建立运动方程解：设解：设杆杆1 1转角转角 杆杆3 3位移位移则有方程：则有方程：1331,(),AlJmMF驱111
4、133 3sin00AMRlJRFm s2331111(cossin)RFm l 2222131311131111sinsincossin0AMF lm lm lJ 1 3sa5 5 1-2 1-2 利用等效力学模型法进行动力学分析利用等效力学模型法进行动力学分析 一、等效力学模型概念一、等效力学模型概念 1 1、思路、思路动能定理：动能定理：EW合外力所做功的增量合外力所做功的增量=系统动能的增量系统动能的增量 质点：质点：)21(2mvddsFa6 6222231122223 31111()2222AssMdFdsdJJm vm v1 12 2、实例：、实例：已知如图，构建动力学方程已知如
5、图，构建动力学方程2222323211222111111()()()()2sAsvvvMFdtdJJmm21112VVMdtdJ等效力矩等效力矩 Mv等效转动惯量等效转动惯量JvM2s3m等效力学模型等效力学模型a7 7力矩与转速同力矩与转速同向取正，反向向取正，反向取负取负1.1.等效力矩等效力矩2.2.等效转动惯量等效转动惯量3.3.等效质量等效质量4.4.等效力等效力以上可以看出，这些等效参数仅与传动比有关，而与真实以上可以看出，这些等效参数仅与传动比有关，而与真实 速度无关。速度无关。1(cos)niiViiiivMMF221()()nsiiViisivJmJ221()()nsiiVi
6、siivmmJvv0()()cos)niiViiiivFMFvv为力与速为力与速度夹角度夹角二、等效参数二、等效参数1ia8 81.1.瞬心法瞬心法2.2.解析法解析法3.3.特例特例 齿轮传动，凸轮传动等齿轮传动，凸轮传动等2421ABBPll1331APvl223121()cos(sin)sflll求传动比方法：求传动比方法：a9 9根据动能定理根据动能定理 有：有：1.1.微分形式微分形式2.2.积分形式积分形式WE 21()2VVM ddJ211222VVVdJdMJdtd 22vvvdJMJd 的函数的函数 的函数的函数212VVVdmFm ssds同理：002211()22VVVo
7、MdJJ002211()22sVVvosF s dsm vm v同理：三、方程形式三、方程形式dJdMvv)21(2a1010例例1 1.已知已知求：角加速度求：角加速度解：解：以构件以构件1 1为等效件为等效件12!212,zzJ JMM1222121211,()VVMMMJJJ21112VVVdJMJd选微分形式：2121VMMM2111121222()/()zzMMJJzz四、典型实例四、典型实例a1111例例2.2.已知从动件的推程方程已知从动件的推程方程 求：凸轮的角加速度求：凸轮的角加速度（略杆的重力）（略杆的重力）解：选凸轮为等效件解：选凸轮为等效件1212,AShJmMF 21
8、2()VVAvMMFvJJm21212212(/)()()(/)VVAAMF hhhMJMFJmJm hhShvS，a1212例例3.3.已知已知 求：求：建立系统运动方程建立系统运动方程（略（略mm2 2，mm2 2g g）解：选解：选1 1为等效件为等效件1331,(),AlJmMF驱313123131()()VVAvMMFvJJm33331cossinsinvSSlSll cossin)sin(sin23223131lmlmJlFMA a1313例例4.4.已知：已知：，略重力及质量略重力及质量求：求：1 1）启动力矩）启动力矩M M1 1最小值；最小值；2 2）如启动）如启动3 3秒后
9、秒后n n1 1=600rpm=600rpm，求，求M M1 1。解：解：1)1)选中心轮选中心轮1 1为等效件为等效件1324212320,40,0.18,0.38,0.22kgm15HHzzzzJJJJMN m11()HvHvMMMJC241141 3113HHz ziiz z 1115/35N mHHMM iH1Hia1414222212311()()()0.8kg mHVHJJJJJ11111115/321.76N m0.8HVVHMJMMMM1121600rpm20 rad/s20/3rad/sn若不忽略齿轮若不忽略齿轮2 2，3 3的质量？的质量？2 2）a.a.若匀速转动若匀速转
10、动MM1 1=？b b.若去掉若去掉M1，多长时间停车？多长时间停车？Ha1515五、运动方程的求解五、运动方程的求解 1.=1.=常数常数VJVM11/VVVVMJMJ或022011()22VVMdJJVMVM22vvvdJMJd()()VVVVddMJdtJdtM00()VVdttJM3 3）为角速度的函数：为角速度的函数：1 1）为常数（用微分形式）：为常数（用微分形式）：2 2）为转角的函数：为转角的函数：a16162.2.不为常数不为常数VJ1 1）=常数常数VM00220011()()22VVVM dJJ 积分形式：02200011()()()22VVVMJJ ()()ddfdtd
11、tf2 2）：利用积分方程：利用积分方程()VVMM3 3）：利用微分方程：利用微分方程()VVMM2()1()2VVVdJMJda1717例例1.1.已知：已知：求：求：1 1）由静止启动）由静止启动5 5秒时蜗杆秒时蜗杆1 1的角速度；的角速度；2 2）若）若 ，其它条件不变，其它条件不变，求蜗杆求蜗杆1 1的角速度。的角速度。123421234142()40,20,30,9,12,8,9(kg m)15N m150 N mzzzzJJJJMM右旋，（驱），（阻）11152M解：解：1 1）314141421515010N mVz zMMMz z22224123411()()()9.06k
12、g mVJJJJJa18181/1.10VVMJ555.5rad/st2 2）分析分析41411102VMMMVJ 不变15111001102VVVVVJddMJdtddtJdtM101159.06ln(102)|3.3356rad/s21 151ta1919例例2.2.已知：弹簧压缩产生的力矩已知：弹簧压缩产生的力矩 求：断电后角速度为求：断电后角速度为0 0时杆的转角时杆的转角,MabJ220011()22abdJJ利用积分形式得：利用积分形式得：2220111222abJJ0=020,aba2020例例3.3.已知：从动件推程方程已知：从动件推程方程 求：凸轮运动参数的变化规律求：凸轮运
13、动参数的变化规律解：选凸轮为等效件解：选凸轮为等效件2212122,AShJmMF2121VvMMF222vSSkk212(2)VAJJmk2221212212(2)(2)22AMk FJmkmk1M2F1 1a2121练习练习：已知：已知：13132()(),()AMFlJmm驱，阻，略求：运动方程求：运动方程分析：选分析：选1 1为等效件为等效件31313211223131321123131122111122231131311325111cos()()cos1,coscos2sin1()4cos2cosVVAAVVVAvlMMFMFvlJJmJmsdJlslMJdvlMFJmm l3131
14、3211223131321123131122111122231131311325111cos()()cos1,coscos2sin1()4cos2cosVVAAVVVAvlMMFMFvlJJmJmsdJlslMJdvlMFJmm l31313211223131321123131122111122231131311325111cos()()cos1,coscos2sin1()4cos2cosVVAAVVVAvlMMFMFvlJJmJmsdJlslMJdvlMFJmm l31313211223131321123131122111122231131311325111cos()()cos1,cosco
15、s2sin1()4cos2cosVVAAVVVAvlMMFMFvlJJmJmsdJlslMJdvlMFJmm la2222 1-3 1-3 利用拉格朗日法进行动力学分析利用拉格朗日法进行动力学分析 一、分析力学的基础知识一、分析力学的基础知识 1.1.分析力学分析力学a23232.2.约束及分类、约束方程约束及分类、约束方程 约束：约束：分类：分类：双面约束（刚杆的约束）双面约束（刚杆的约束）单面约束（绳子的约束）单面约束（绳子的约束）完整约束（几何约束）完整约束（几何约束）非完整约束（运动约束）非完整约束（运动约束）稳定约束（定常约束）稳定约束（定常约束）非稳定约束（非定常约束）非稳定约束（
16、非定常约束）对位置进行限制的约束对位置进行限制的约束-对速度、加速度进行限制对速度、加速度进行限制对构件的位置或运动进行限制对构件的位置或运动进行限制根据约束对限制的不同情况：根据约束对限制的不同情况：-不随时间变化而变化不随时间变化而变化-随时间变化而变化随时间变化而变化-用等式方程表示的约束用等式方程表示的约束-用不等式方程表示的约束用不等式方程表示的约束u 约束方程：约束方程：将约束条件用数学形式表示出来的方程将约束条件用数学形式表示出来的方程a24245.5.理想约束：理想约束：6.6.广义坐标广义坐标 ：10niiiRr这里的广义坐标是杆这里的广义坐标是杆1 1转角还是转角还是B B
17、点直角坐标点直角坐标,为什么？为什么？在任意虚位移上系统约束反力所作元功在任意虚位移上系统约束反力所作元功之和为零（略摩擦）之和为零（略摩擦）用以确定机构位置的一组独立参数用以确定机构位置的一组独立参数qa25257.7.自由度：自由度：8.8.广义速度：广义坐标广义速度：广义坐标q q 对时间对时间t t 的一阶导数的一阶导数 12,.,nq qq12,.,nq qq广义坐标的独立变分数目广义坐标的独立变分数目自由度数自由度数在完整系统中，广义坐标数在完整系统中，广义坐标数=独立变分数独立变分数=自由度数自由度数a2626例：例：如图平面机械手如图平面机械手广义坐标：广义坐标：112233,
18、qqqm m点坐标：点坐标：112233coscoscosxlll112233sinsinsinyllliiiiirrrqtq偏导数中广义坐标是相偏导数中广义坐标是相互独立的，均为时间互独立的，均为时间t t的函数的函数)sinsinsin(333222111lllx333222111coscoscoslllyjqjq a27279.9.广义加速度广义加速度 ：10.10.虚位移原理：虚位移原理：q 22211nniiiijkjkjkrrrrqq qtqqq 证明证明稳定理想约束系统处于平衡的充分必稳定理想约束系统处于平衡的充分必要条件是在任一给定虚位移上，主动要条件是在任一给定虚位移上，主动
19、力所做元功之和为零。力所做元功之和为零。jnjjiqqr 1a2828虚位移原理的表达形式：虚位移原理的表达形式：0iiWFr形式形式1 1：形式形式2 2：形式形式3 3：广义坐标表达式广义坐标表达式01jNjjqQW0)(iiziiyiixzFyFxFjQ广义力广义力a2929例例1 1：已知如图，已知如图，1112112231122coscoscoscosssxlxllxll1111122211221233112212sin,0sin,sinsin,sinssxxlqqxxllqqxxllqq 广义力：广义力：111121111222222(sin)(sin)(sin)(sin)(sin
20、)ssQm glm glFlQm glFl1112112231122coscoscoscosssxlxllxll1112112231122coscoscoscosssxlxllxll1111122211221233112212sin,0sin,sinsin,sinssxxlqqxxllqqxxllqq 1111122211221233112212sin,0sin,sinsin,sinssxxlqqxxllqqxxllqq 111121111222222(sin)(sin)(sin)(sin)(sin)ssQm glm glFlQm glFl1111122211221233112212sin,0
21、sin,sinsin,sinssxxlqqxxllqqxxllqq 1111122211221233112212sin,0sin,sinsin,sinssxxlqqxxllqqxxllqq 1111122211221233112212sin,0sin,sinsin,sinssxxlqqxxllqqxxllqq 2cos求广义力。求广义力。,2121SSllll,21（垂直向下）Fgmgm解：解：2211,qq取a3030问题：如有力矩问题：如有力矩M M是否影响广义力？是否影响广义力？广义力应用的是虚位移广义力应用的是虚位移原理，所以有影响。原理，所以有影响。()0iiiiF rMa3131例
22、例2 2：已知已知求：广义力求：广义力1123410N m9N m20,40,20,80HMMzzzz，解解：自由度数：自由度数=广义坐标数广义坐标数取取11q11niHiiHiiiQMMMqq11111/HHiHHHHitQMMMMMiMqt241141 311()9HHz ziiz z 1110(9)9N m9Q jQjq1q1QHi1a3232例例3 3：解：解：设设A A点虚位移点虚位移BCBC杆虚位移杆虚位移cosBCaCECE杆位移杆位移cos(902)cosCE2Ma tgF已知六杆机构中的力已知六杆机构中的力F F，求平，求平衡时的驱动力矩衡时的驱动力矩M M。虚位移原理应用虚
23、位移原理应用用于解决静力学问题用于解决静力学问题则：则：0EFMa3333例例4 4：已知已知求：平衡时，求：平衡时，112342N m20,40,20,80Mzzzz解解：分析：分析0iiQ q1124,qq111018N mHHHQMMM 24440=16N mHHQMMM4,?HMM取取因广义坐标为独立参数，不互相影响因广义坐标为独立参数，不互相影响轮轮4 4不动，轮不动，轮1 1有虚位移，得：有虚位移，得：轮轮1 1不动，轮不动，轮4 4有虚位移，得：有虚位移，得：1/98/9a3434惯性力为惯性力为 ，11.11.动力学普遍方程（第一类拉格朗日方程）：动力学普遍方程（第一类拉格朗日
24、方程）：imiFiima1()0niiiiiFmar动力学普遍方程：具有理想约束的质点系运动时，动力学普遍方程：具有理想约束的质点系运动时，在任一瞬时作用在质点系上的所在任一瞬时作用在质点系上的所 有有主动力和惯性力主动力和惯性力在在任意虚位移任意虚位移 上所做的上所做的元功之和等于零元功之和等于零。若系统具有理想约束，并由若系统具有理想约束，并由n n个质点组成，个质点组成，任一质点为任一质点为 ，主动力为主动力为 ，根据虚位移原理在任一瞬时有：根据虚位移原理在任一瞬时有：0)(iixiixxamF0)(iiyiiyyamF0)(iiziizzamFa3535例：例：121,HHr r m
25、mM用功率表示功用功率表示功111111212()()0HAHHHHMJJm rrrr 又又1212112111212111111,HHHHHzrriizrrrrrrr 212()3HAmrrJ211 112Jm r已知标准齿轮标准安装，系统在已知标准齿轮标准安装，系统在水平面内运动，水平面内运动，求：运动与受力关系求：运动与受力关系分析：分析：212113()()32HHHMrrmmH a363612.12.第二类拉格朗日方程：第二类拉格朗日方程：设理想、完整约束系统由设理想、完整约束系统由n n个质点组成，个质点组成，12(,.,)(1,2,.,)iiNrr q qqin上式变分得上式变分
26、得（变分运算如同微分运算，进行微分运算后，变分运算如同微分运算，进行微分运算后，将微分符号改为变分符号将微分符号改为变分符号）1Niijjjrrqq矢径对时间求导矢径对时间求导1Niiiijjjrrrvqtq有有N N个自由度，个自由度，系统中任一质点的矢径可表示为：系统中任一质点的矢径可表示为：a3737将上式对将上式对 求偏导有：求偏导有：iiijjjrvrqqqjq 将上式对将上式对t t求全导数：求全导数：()iiijjjrrvddtqqq将第一类拉氏方程打开，有：将第一类拉氏方程打开，有：11110,()nnnniiiiiiijjjiiijjrF rmarF rQqQFq惯性力所做元
27、功之和：惯性力所做元功之和：11111nnNNniiiiiijiijiijjijjrrm arm aqm aqqq1Niiiijjjrrrvqtqa3838()()iiijiiiiiijjrd mvqrrdmamvqdtdtq和和带入有：带入有：2211()()()22()()iiiiiiijiiiiiijjjjrd mvmvmvqrrddmamvqdtdtqdtqq引入系统动能：引入系统动能：2112niiiTmv得得a3939由于广义坐标的变分都是独立的，因此上式中必有：由于广义坐标的变分都是独立的，因此上式中必有：说明：说明：1 1、拉氏方程是一个由、拉氏方程是一个由N N个方程组成的二
28、阶方程个方程组成的二阶方程 组，其特点是不含约束反力。组，其特点是不含约束反力。2 2、拉氏方程是以能量的角度研究问题，、拉氏方程是以能量的角度研究问题，因此避免了加速度的分析。因此避免了加速度的分析。3 3、方程表明了动能变化和主动力之间的关系。、方程表明了动能变化和主动力之间的关系。拉氏方程：拉氏方程：a4040例例1 1：已知标准齿轮标准安装，系已知标准齿轮标准安装，系统在水平面内运动统在水平面内运动求：运动与受力关系求：运动与受力关系121,HHr r m mM分析：系统具有一个自由度分析：系统具有一个自由度又又HHQM212113()()32HHHMrrmm二、利用拉式方程进行动力学
29、分析二、利用拉式方程进行动力学分析取取HHq211212212121 JvmJTBHA22211222112)(21)(2121HHHArrrJrrmJB Ba4141例例2.2.已知：已知：求：用拉格朗日方程动力学求：用拉格朗日方程动力学方程方程1331,(),AlJmMF驱解：解：系统一个自由度，取系统一个自由度，取系统动能：系统动能：22113 31122ATJm v33111(cos)/sinvsd ldtl 则则2231111sincosTm l 11q2112231sin21lmJTAa4242从虚功率角度求广义力从虚功率角度求广义力1 122.NQqQ q此机构的虚功率：此机构的
30、虚功率：11331311(sin)NMFsMFl由拉氏方程：由拉氏方程：得得:222213113111131(sin)sincossinAJm lm lMFl 也可由虚功来求也可由虚功来求Q Q1 1a4343课堂练习课堂练习已知：已知：132421322421H20,60,0.1kg m,0.4kg m,0.2kg m,12.75N mHzzzzJJJJJM，求？(略2、3质量及重力影响)用：等效法、拉氏一法、拉氏二法mNMH16a44441 1、等效法：、等效法：选选H H为等效件为等效件1241141 3111 98HHHz ziiz z 423343111 32HHHziiiz 等效力
31、矩：等效力矩：1112.75 8 1686N mVHHMMM 2222231123()()()0.1 80.5 20.28.6kg mVHHHJJJJJ因为因为 为常数，选微分方程为常数，选微分方程VJ2/10rad/sVVHHVVMJMJa45452 2、动力学普遍方程（拉氏一法）：、动力学普遍方程（拉氏一法）：给定：给定：111112322()0HHHHHMJJJMJ 11228822HHHH HH a46463 3、拉氏二法：、拉氏二法：取：取：22211232111()222JHTJJJJ128,2HH 221(64 0.14 0.50.2)4.32HHT HHHHqq，a4747广义
32、力的求法一般有两种方法广义力的求法一般有两种方法:1 1、按广义力定义求解、按广义力定义求解2 2、采用虚功方法进行求解、采用虚功方法进行求解由于采用虚功方法进行求解相比较而言容易由于采用虚功方法进行求解相比较而言容易一些，因此本课程中涉及到广义力的求解都一些，因此本课程中涉及到广义力的求解都是采用虚功方法。是采用虚功方法。a4848例例：如图示机构，求平衡时：如图示机构，求平衡时机构自由度数为机构自由度数为3 3，构件，构件1 1、4 4、7 7运动定义为广义坐标，即运动定义为广义坐标，即371,F F M112437,qqqS平衡时，在平衡位置的虚功为零，平衡时，在平衡位置的虚功为零，又广
33、义力为零又广义力为零123.0.0.0QQQ可以求出三个未知数可以求出三个未知数分析：分析：a4949例例：五杆机构，取：五杆机构，取1124,qqS构件构件1 1由由 控制，控制，构件构件4 4由由 控制，控制，件件2 2、3 3由由 共同控制。共同控制。第二章第二章 两自由度机构动力学分析两自由度机构动力学分析 2-1 2-1 两自由度机构的运动分析两自由度机构的运动分析1.1.构件上某点速度：构件上某点速度：上式也可以表示为：上式也可以表示为：分析：分析：)(11q02q)(42sq01q21,qqir称为类线速称为类线速度（矢量）度（矢量）a5050的物理意义：的物理意义：当当 时，时
34、，的大小、方向即为的大小、方向即为 的大小方向的大小方向量纲由广义坐标决定量纲由广义坐标决定iv2.2.构件角速度构件角速度注意：角速度在平面机构中为标量，在空间机构中为矢量注意：角速度在平面机构中为标量，在空间机构中为矢量21,iiuu:1 iu1,012qq1 iu)(ii),(21qqii如研究杆如研究杆2 2、杆、杆3 3：不是传动比不是传动比21,iiii第第i i个件对广义坐标个件对广义坐标1 1，2 2的的类角速度类角速度（标量）（标量）的物理意义？的物理意义？21,iiiia5151212122222121212122222121221221qqiiqiqiJqquuququm
35、iiiiisisisisisi 2-2 2-2 利用拉格朗日方程建立两自由度机构利用拉格朗日方程建立两自由度机构 的动力学方程的动力学方程1 1221 122iiiiiivu qu qi qi q 拉格朗日方程：拉格朗日方程：()jjjdTTQdtqq一、惯性系数一、惯性系数1q 1q 2q 2q 求求1 1个构件动能：个构件动能：a5252对于对于 ：与与 和和 均相关件的质量和转动惯量才能计入，均相关件的质量和转动惯量才能计入，则系统动能：则系统动能：说明：说明：11J1q12J1q2q对于对于 ：件件i 的运动必须与的运动必须与 有关，有关，即与即与 相关件的质量和转动惯量才能计入，相关
36、件的质量和转动惯量才能计入，1q惯性系数惯性系数11J为正；为正；12J可为正、为负、为零。可为正、为负、为零。a5353例例1 1：已知：已知：求：求：12121212,ssssm mJJl l ll111222,JJJ分析：广义坐标可以设为：分析：广义坐标可以设为：1122,qq则：则：,111q211101qq,21212qq,21212qq则：则：a5454222111112111222(2cos)sssssJmlJm lll lqJ21222 12 22 21 22221222(cos)ssssssJm uuJ i im ll lqJ,21sl2212221cos2qllllss22
37、222222iJumssa55551324132420,40,0.1,0.25zzzzJJJJ111222,JJJ例例2 2：已知差动轮系已知差动轮系轮轮2 2、3 3质量略，质量略，H转动惯量略。转动惯量略。求：求：分析：广义坐标可以设为：分析：广义坐标可以设为：112,Hqq11121,0ii则：则：a5656a5757二、计算动能二、计算动能用惯性系数表示的动能：用惯性系数表示的动能：11 11221TJ qJ qq22112212121211111122JJJTqqq qqqqq a5858则拉氏方程为：则拉氏方程为：22112212121211111122JJJTqqq qqqqq
38、两个自由度的两个自由度的拉氏方程拉氏方程a5959例例3 3：已知：已知：求：建立运动方程求：建立运动方程分析：选广义坐标分析：选广义坐标:1122,qq则：则：u 求类线速度求类线速度:F，a6060常数常数a6161求广义力：求广义力：方程：方程：此为二阶非线性微分方程，用数值解法求解。此为二阶非线性微分方程，用数值解法求解。a62621,H 例例4 4：已知差动轮系中：已知差动轮系中：，各轮质量略。，各轮质量略。112,Hqq11212,2,3,qq q分析：取广义坐标：分析：取广义坐标：则：则：2qH1H求：求：a6363方法方法1 1：方法方法2 2：同理同理求：求：2221,ii2
39、123,212221ii,02q 令即即H H不动，则：不动，则：,01q 令即即1 1轮不动，则：轮不动，则：2212iiH2223i求：求：3231,ii)81(1a646411 1122121 12222J qJ qQJ qJ qQ计算广义力：计算广义力：动力学方程：动力学方程：差动轮系动力学方程，可以直接应用此结论式。差动轮系动力学方程，可以直接应用此结论式。22HHiJa656512212,AsJm JM F12,qqr例例5 5：已知：已知：重力略，建立运动方程。重力略，建立运动方程。分析：分析：选广义坐标：选广义坐标：则：则：a666611 1212222222J qm r rQ
40、J qm rQ计算广义力：计算广义力：动力学方程：动力学方程：a6767iiidTTQdtqq0iidLLdtqqiiidLLQdtqq例例第三章第三章 多自由度机构的动力学分析多自由度机构的动力学分析 3-1 拉格朗日方程拉格朗日方程a68681 12211 12233222133332()()()()m xkxxk xm xk xxkxxm xFk xxa6969 3-2 多自由度机构的动力分析多自由度机构的动力分析一、运动关系一、运动关系1 1、某构件运动与一个广义坐标相关、某构件运动与一个广义坐标相关2 2、某构件运动与几个广义坐标相关、某构件运动与几个广义坐标相关3 3、各构件在广义
41、坐标下的表示、各构件在广义坐标下的表示4 4、构件速度、角速度表示、构件速度、角速度表示5 5、构件质心的坐标、速度表示、构件质心的坐标、速度表示()1 1()22()33.isisisisvuququq类线速度类线速度1 12233.isiiii qi qi q类角速度类角速度a7070二、系统动能二、系统动能221122iiisisiTmvJ21122xxxyxzxiiisxiyizixyyyyzyixzyzzzziijjjmvjjjjjj4222211 1222333444121211313141423232424343411112222iiTJ qJ qJ qJ qJ q qJ q q
42、J q qJ q qJ q qJ q q a7171以平面以平面4 4自由度为例自由度为例(表格形式表格形式)：11 1J q 122J q 133J q 144J q 21 1J q 222J q 233J q 244J q 31 1J q 322J q 333J q 344J q 41 1J q 422J q 433J q 444J q iiidTTQdtqqa72721NkjjjJ q22111()22Nkjjjkkjkjjkkj kJJJqqqqq1Nkkjkjjj kJq qq 111()NNkMkLMLMLML MLMkMk L kJJJq qqqq a7373 的下标的含义：与的
43、下标的含义：与i i、j j广义坐标同时有关的构件广义坐标同时有关的构件的等效质量或惯量。的等效质量或惯量。ijJ23232323()ii xi xi yi yisiiJm uuuuJii如如以以3 3自由度第自由度第3 3个方程为例：个方程为例：31 13223332223132331122123132333333313212313212123213111()()222()J qJ qJ qJJJJJqqqqqqqqJJq qq qqqJJJq qQqqq 1NkjjjJ q22111()22Nkjjjkkjkjjkkj kJJJqqqqq1Nkkjkjjj kJq qq 111()NNkM
44、kLMLMLML MLMkM k L kJJJq qqqq a7474空间任一运动的刚体空间任一运动的刚体证明证明21122xxxyxzxiiiisxiyizixyyyyzyixzyzzzziijjjTmvjjjjjja757522222222222112211()()(2)22112211()()02211()()(22iCiCiiCCiiCiiCCiiiCiiCiiiCixyzyxzzxTm vm vvm vvvvm vvv vvvmvvmv vvmvmvmrrmvmrrrrr 222222)2()11(222)221122yxy x yxz x zyz y zCxxxyyyzzzxyxy
45、yzyzxzxzTCrr rr rr rmvJJJJJJmvJ 如果质心速度如果质心速度为零，刚体动为零，刚体动量也为零量也为零根据转动惯根据转动惯量计算公式量计算公式a76762211()22iisisiTmvJ系统动能：系统动能：三、系统势能三、系统势能势能只与位置有关，即仅与广义坐标本身有关，因势能只与位置有关，即仅与广义坐标本身有关，因此在系统运动明确之后，势能也可求得，一般在拉此在系统运动明确之后，势能也可求得，一般在拉格朗日方程中用格朗日方程中用“U U”表示。表示。四、广义力四、广义力广义力一般用虚位移原理求得，如果系统仅有广义力一般用虚位移原理求得，如果系统仅有有势有势力力做功
46、，引入拉氏函数广义力为零（如一些震动系做功，引入拉氏函数广义力为零（如一些震动系统）。统）。引入拉氏函数后广义力不包括有势力引入拉氏函数后广义力不包括有势力常见势能有常见势能有哪些？哪些？a7777例例1 1：如图，已知各转动惯量、力矩如图，已知各转动惯量、力矩其余略，求动力学方程其余略，求动力学方程124520zzzz分析：系统自由度为：分析：系统自由度为：112336,qqq12311221231,3,62,4,5,qqqHq qHq q q3660zz3 72 743F 设：设：a7878112336,QM QMQM113263,广义力用虚功广义力用虚功原理求解原理求解动能均为角速度动能
47、均为角速度(广义速度广义速度)的函数，的函数，121122433,24H21233125312123,168Ha7979注：注：轮系中，一般类角速度是轮系中，一般类角速度是定值。所以有惯性系数为定值。定值。所以有惯性系数为定值。11 1122133121 1222233231 13223333J qJ qJ qQJ qJ qJ qQJ qJ qJ qQ1NkjjjJ q22111()22Nkjjjkkjkjjkkj kJJJqqqqq1Nkkjkjjj kJq qq 111()NNkMkLMLMLML MLMkM k L kJJJq qqqq a8080例例2 2：如图，杆长已知，质心如图，杆
48、长已知，质心位置已知，各杆受力矩、转动位置已知，各杆受力矩、转动惯量已知。建立系统动力学方惯量已知。建立系统动力学方程。程。分析分析1 1：系统为平面：系统为平面N N自由度自由度开链机构，广义力为重力、开链机构，广义力为重力、外力矩和手爪部外力。外力矩和手爪部外力。分析分析2 2：动能函数为质心速度、：动能函数为质心速度、角速度函数，势能为广义坐角速度函数，势能为广义坐标函数。标函数。问题问题1 1：广义力如何求？：广义力如何求？问题问题2 2：T T或或L L函数的表达？函数的表达？思考：动能、势能思考：动能、势能的广义力表达式的广义力表达式a81811NkjjjJ q22111()22N
49、kjjjkkjkjjkkj kJJJqqqqq1Nkkjkjjj kJq qq 111()NNkMkLMLMLML MLMkM k L kJJJq qqqq 12121212()ii xi xi yi yisiiJm u uuuJii221111112nsskkJJmlm l各杆转动部分仅与各各杆转动部分仅与各自的广义坐标有关。自的广义坐标有关。21122112211 2211 223(sin)(sin)coscos(sinsincoscos)ssnkkmllllm llll2 12121 2123cos()cos()nskkm llm lla8282广义力广义力:1Q 1M121sinsm
50、gl112sinnkkm gl11sinxF l11cosyF l通式通式:211nnkjkjjjjjjjj kJJ qqQqa8383例例3 3：如图已知：如图已知：其余略，求动力学方程。其余略，求动力学方程。135720zzzz246860zzzz213570.27kg mJJJJ224680.81kg mJJJJ121488N m,9N mHHMMMMMa8484分析分析1 1：系统自由度数：系统自由度数？分析分析2 2：各构件与广义坐标关系：各构件与广义坐标关系112132,HHqqq1122312231;,5;2,3,4,6,7,8,q Hq Hqq qq q11 1122121 1
