用拉格朗日方程建立系统数学模型.doc

第二章 用拉格朗日方程建立系统的数学模型2.1概述拉格朗日方程属于能量法，推导中使用标量，直接对整个系统建模特点：列式简洁、考虑全面、建模容易、过程规范适合于线性系统也适合于非线性系统，适合于保守系统，也适合于非保守系统。2.2拉格朗日方程1 哈密尔顿原理系统总动能（21）系统总势能（22）非保守力的虚功（23）哈密尔顿原理的数学描述： （24）2 拉格朗日方程：拉格朗日方程的表达式： （25）（推导：）将系统总动能、总势能和非保守力的虚功的表达式代入哈密尔顿原理式中（变分驻值原理），有（26）利用分步积分 （27）并注意到端点不变分（端点变分为零） （28）故 （29）从而有（210）由变分学原理的基本引理：（设 n维向量函数M(t)，在区间内处处连续，在内具有二阶连续导数，在处为零，并对任意选取的n维向量函数，有则在整个区间内，有 ）我们可以得到：（211）即（212）对非保守系统，阻尼力是一种典型的非保守力，如果采用线性粘性阻尼模型，则阻尼力与广义速度成正比，在这种情况下，可引入瑞利耗散（耗能）函数D， （213）阻尼力产生的广义非保守力为： （214）对于仅受有势力和线性阻尼力作用的系统，其拉格朗日方程为： （215）如果系统上还作用了除有势力和阻尼力以外的非保守力，如结构受到的外激励力（对应的广义非保守力可通过非保守力的虚功求得，仍记为），则系统的拉格朗日方程为： （216）2.3 拉格朗日方程在振动系统建模中应用在某些结构振动问题中，取分离体、确定各分离体的受力情况，然后利用牛顿第二定律建立方程的方法不一定可用，或者很不方便，这时，采用拉格朗日方程来建立振动方程就很方便。LuO1 集中参数模型中应用【例】质量为M的长直杆上有一个集中质量m可在杆上滑动。杆绕固定点摆动，建立其自由振动方程。势能( 以O点为势能零点)动能选广义坐标为，且，代入拉格朗日方程得到：以上是对离散系统应用拉格朗日方程建立振动方程，如果利用拉格朗日方程建立连续系统的方程，则它是一种同时将系统离散化、变量分离并达到系统降阶的途径。2 连续参数模型中应用与假设模态法联合使用对一维连续系统，假设位移为：（217）则系统具有N个自由度，N个广义坐标为，不一定是系统的真实模态，可以是假设的一种变形模态。只要满足以下条件：（1） 是位移形函数，反映某种可能的位移形状（2） 构成一组线性无关向量（3） 连续导数阶次满足势能中所要求的阶次（4） 满足位移边界条件（不一定满足力边界条件）p(x,t)Lxu(x,t)2.1 杆的纵向振动轴向位移为 (214) (215)将代得到： (218)(219)其中 (220)分布轴力p(x，t) 在广义坐标上的虚功 （221）广义力 （222）代入拉格朗日方程得：(223)或xyp(x,t)u(x,t)L（224）（2008-3-26）2.2 梁的横向振动横向位移函数 （222）动能 （225）势能 （226）， （227）分布外力做的功：（228） （229）代入拉格朗日方程： （230）或矩阵方程： （231）注意假设模态法与有限元素法的区别：这里的是对整个结构的假设模态（相当于整个结构变形的形函数），不是单元的位移形函数，对复杂结构，确定精度（品质）较高的假设模态是比较困难的。3. 粘性阻尼系统中阻尼的处理假设结构中具有分布粘性阻尼力 （232）广义力 （233） （234）代入拉格朗日方程得到 （235）上式中为其他的广义非保守力2.4 坐标约束与拉格朗日乘子通常对一个N 维结构系统，采用拉格朗日方程建立振动方程时，广义坐标是线性独立的，但是实际问题中，有时希望采用一套不是独立的坐标来建立方程，可能更加方便。记系统不独立的坐标为则被约束坐标数C=M-N（236）对广义坐标，有C个约束方程：（237）如果令每一个坐标取变分，则：（238）（239）上式说明这M个不独立，而是由上述C个方程联系起来。在哈密尔顿原理式中，将坐标数由N扩展到M，即得到：（240）注意，由于此时的不独立，不能直接由变分学基本原理，得出方括号内的项等于零的结论。对上面的约束方程引入拉格朗日乘子（或称为拉格朗日乘子函数），得到：（241）代入哈密尔顿原理方程式中，（242）我们可以选择C个，使C个相应的方括号表达式为零，那么其余N=M-C个独立的对应的方括号内的项必为零。从而得到带约束的拉格朗日方程（修正的拉格朗日方程）为： (243)联立上两个方程，就可确定M+C个未知数【应用实例】Lxu(x, t)求两端固定杆的轴向自由振动微分方程。【解】令， （244）即假设模态为 （245）约束边界条件：（246）第一个条件由形函数满足，第二个条件实际为：（247）这就是约束方程。根据2.3节的公式（2-20），可以求出轴向振动的杆的质量矩阵和刚度矩阵为： (248)现在用修正的拉格朗日方程来建立方程:本例只有一个约束方程，故只需一个拉格朗日乘子，即在拉格朗日方程中引入和项，且（249）代入修正的拉格朗日方程中，并联立约束方程得到：（250）可由此解出。2.5 受约束结构的振动此处的约束是指结构的附加惯性约束或附加弹性约束，即一个结构由于添加了质量或弹簧而对结构表现为一种约束，而不是指通常的坐标约束。一般说来，给结构添加一个弹簧，弹簧将对结构的运动表现为一种弹性约束而使系统的固有频率增加，相反，添加一个质量，也表现为对系统的惯性约束，但使其固有频率降低。看一个例子：对于一个未受约束的一维结构，受分布力和分布力矩作用，现求其强迫振动。假定其主模态、固有频率已知，则其任一点处的挠度可以表示为：（251）代入拉格朗日方程可得到广义坐标满足的方程： (252)为对应主模态的广义质量。如果在处还作用了集中力和集中力矩，则相应的广义力虚功由下式确定： （253）则广义力为： （254）所以，运动方程为： (255)方程 (253)就是本节分析受约束结构振动的基本方程。当结构上处添加一个刚度系数为的弹簧时 （256）当结构上处添加一个集中质量时 （257）【注意】在讨论受约束结构时，均假定未受约束结构的模态参数是已知的，相当于在结构有修改时，求解修改后结构的振动问题。Kka【例】如图所示的一个简支梁，受弹簧约束，求其运动方程对线弹簧， （258）对扭簧， （259）代入上基本方程： （260）受弹性约束后结构的主振动仍然是简谐的，所以：（261）代入上方程得到：（262）由上式中N个的系数行列式为零，就得到受约束结构的频率方程。a【例】如图所示，在简支梁的处添加一个质量，求运动方程。（263）与上例推导相同，可得到：（264）从而： （265）当仅取未受约束简支梁的第一阶振型时，上方程简化为：（266） （267）由此式可以看到，如果添加的质量放置在该阶振型的节点处，由于，故此时。同样，如果添加的弹簧放置在该处，同样不会对固有频率产生影响。因此，要通过附加质量或附件弹簧（刚度）的方式来对结构进行修改，要避开原结构固有振型的节点。当然，实际中质量和刚度元件