矩形截面梁运动控制方程的建立

《矩形截面梁运动控制方程的建立》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、矩形截面梁运动控制方程的建立在建立控制微分方程的时候，本文主要应用了弹性力学中的能量变分原理和Reissner原理。（1）能量变分原理弹性梁的非线性振动问题在数学上要通过积分控制方程来精确求解是相当困难的，且大多数情况下不能求得精确解。因此，人们就致力于用能量变分法或其他方法来近似求解，而避免了控制微分方程的困难。能量变分法有两种用处，一是可以用能量泛函通过变分近似求解薄板的各种边值问题，而无需建立板的控制微分方程。二是薄板的控制微分方程及边界条件可以通过那难量法来建立，特别是疑难的边界条件，如自由边的边界条件就是克希霍夫用能量法建立的。能

2、量变分法是近似解法中最有成效的方法之一，而且，它是半解析法和有限元等数值计算方法的理论基础。对于能量变分法，它的本质其实就是把求解弹性力学的基本方程的定解问题，变化成求泛函数的极大极小值（或驻值）问题。而在求解问题的近似解的时候，泛函数的极大极小值（或驻值）问题又进而变为函数的极大极小值（或驻值）问题。最后可以把问题归纳总结为求解线性代数方程组的问题。（2）Reissner原理如果弹性体处于运动状态，则根据达朗贝尔原理，在个质点上加上惯性力后，就可以将运动问题当作静力问)题来处理。2.1梁的理论及其应用建立并研究如下图所示的矩形截面梁，假设

3、截面高度h和截面宽度b的尺寸都远远小于梁的长度l，同时假设梁的上表面承受横向均布荷载qx。图1梁的模型图则我们可以写出梁在直角坐标系下的位移场为0uu(x,t)z(x,t)ww(x,t)（2-1）式中0u(x,t)——梁中线沿x方向的位移。w(x,t)——梁中线沿z方向的位移。u——梁沿方向的位移。w——梁z方向的位移。w(x,t)(x,t)(x,t)——梁的横向扭转角，且对浅梁有x。2．2梁的运动控制方程的推导2.2.1基本方程的推导(1)建立梁的几何方程分析过程中考虑梁的几何非线性因素，则根据经典非线性弹性理论，

4、梁内任一点在某一瞬时的应变可由格林应变张量来描述，其在笛卡儿坐标系中的分量形式为：1(uuuu)（2-2）iji,jj,ik,ik,j2则有00zkxxx（2-3）0xzxz其中0012uwx,x,x20wxz,x00而扭率k，对浅梁有kw。x,xx,xx(2)建立弹性体（梁）的本构关系根据经典梁理论，其在笛卡儿坐标系下的广义内力可定义为h2Ndzxhx2h2Mzdz（2-4）xhx2h2Qdzxhxz2式中N——梁的轴向力；xM——梁所承受的弯矩；xQ——梁所承受的横向

5、剪力。x则弹性体（梁）的本构关系可表示为如下形式0NA00（2-5）M0Dk及QC（2-xxz6）式中A——梁的薄膜刚度；D——梁的弯曲刚度；CGh,G是梁的横向剪切刚度。000且,kx,x2．2．2Reissner变分原理及其应用(1)弹性体（梁）运动（平衡）微分方程采用Reissner变分原理建立弹性体（梁）的运动平衡微分方程，以及相应边界条件，首先引入Reissner函数B()dVfudVpudS（2-7）ijijijiiiiVVA式中

6、B()弹性体的余能密度，f为沿坐标i方向上每单位体积内的体积力，p为ijii沿坐标i方向上每单位面积所承受的表面力，V为弹性体所占的空间，S为弹性体表面上面力已被给定的那部分面积。若设梁高为h，梁的质量密度为，梁在上表面承受分布荷载q(x,t)。将以0上（2-1）、(2-3)、(2-4)、(2-5)各式代入（2-7）,并对其沿高度积分后，（2-7）式右边第一项变为T0N1bl0QxxzMkT11NA0N1QxGQxdx2M0DM2（

7、2-8）又根据达朗贝尔（DalembertJbe.R）原理，将惯性力u当作分布力，而忽略质itt,量体力，则(2-7)式第二项变为300hbhuuwwdx（2-9）2l0,tt,tt12,tt又令N、N和M是已经给定的沿梁的周边所施加的外力和外力矩，则（2-xxzx7）式右边第三项变为03blq(x,t)wdxS(NuM)ds（2-10）根据Reissner变分原理，有0，即123)0将（2-8）、（2-9）、（2-10）式代入上式，并结合上文的推导，可解得用应力表示

8、的梁的运动（平衡）微分方程0Nhux,x0,ttQNwhwq（2-11）x,xx,x,x0,tt3h0MQx,xx,tt12及相应的边界条件0NN或u