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

矩形截面梁运动控制方程的建立 在建立控制微分方程的时候 本文主要应用了弹性力学中的能量变分原理和 Reissner 原理 1 能量变分原理 弹性梁的非线性振动问题在数学上要通过积分控制方程来精确求解是相当困 难的 且大多数情况下不能求得精确解 因此 人们就致力于用能量变分法或其 他方法来近似求解 而避免了控制微分方程的困难 能量变分法有两种用处 一 是可以用能量泛函通过变分近似求解薄板的各种边值问题 而无需建立板的控制 微分方程 二是薄板的控制微分方程及边界条件可以通过那难量法来建立 特别 是疑难的边界条件 如自由边的边界条件就是克希霍夫用能量法建立的 能量变 分法是近似解法中最有成效的方法之一 而且 它是半解析法和有限元等数值计 算方法的理论基础 对于能量变分法 它的本质其实就是把求解弹性力学的基本 方程的定解问题 变化成求泛函数的极大极小值 或驻值 问题 而在求解问题 的近似解的时候 泛函数的极大极小值 或驻值 问题又进而变为函数的极大极 小值 或驻值 问题 最后可以把问题归纳总结为求解线性代数方程组的问题 2 Reissner 原理 如果弹性体处于运动状态 则根据达朗贝尔原理 在个质点上加上惯性力 后 就可以将运动问题当作静力问 题来处理 2 1 梁的理论及其应用梁的理论及其应用 建立并研究如下图所示的矩形截面梁 假设截面高度 h 和截面宽度 b 的尺寸都 远远小于梁的长度l 同时假设梁的上表面承受横向均布荷载 xq 图 1 梁的模型图 则我们可以写出梁在直角坐标系下的位移场为 0 txww txztxuu 2 1 式中 0 txu 梁中线沿x方向的位移 txw 梁中线沿z方向的位移 u 梁沿方向的位移 w 梁z方向的位移 tx 梁的横向扭转角 且对浅梁有 x txw tx 2 2 梁的运动控制方程的推导梁的运动控制方程的推导 2 2 1 基本方程的推导 1 建立梁的几何方程 分析过程中考虑梁的几何非线性因素 则根据经典非线性弹性理论 梁内任 一点在某一瞬时的应变可由格林应变张量 来描述 其在笛卡儿坐标系中的分量形 式为 2 1 jkikijjiij uuuu 2 2 则有 0 00 xzxz xxx zk 2 3 其中 xxz xxx w wu 0 2 0 0 2 1 而扭率 xx k 0 对浅梁有 xxx wk 0 2 建立弹性体 梁 的本构关系 根据经典梁理论 其在笛卡儿坐标系下的广义内力可定义为 2 2 2 2 2 2 h h h h h h dzQ zdzM dzN xzx xx xx 2 4 式中 x N 梁的轴向力 x M 梁所承受的弯矩 x Q 梁所承受的横向剪力 则弹性体 梁 的本构关系可表示为如下形式 0 0 0 0 k D A M N 2 5 及 xzx CQ 2 6 式中 A 梁的薄膜刚度 D 梁的弯曲刚度 GhC G 是梁的横向剪切刚度 且 xx k 000 2 2 2 Reissner 变分原理及其应用 1 弹性体 梁 运动 平衡 微分方程 采用 Reissner 变分原理建立弹性体 梁 的运动平衡微分方程 以及相应边界 条件 首先引入 Reissner 函数 dVB V ijijij dSupdVuf VA iiii 2 7 式中 ij B 弹性体的余能密度 i f为沿坐标i方向上每单位体积内的体积力 i p为 沿坐标i方向上每单位面积所承受的表面力 V 为弹性体所占的空间 S为弹性体 表面上面力已被给定的那部分面积 若设梁高为h 梁的质量密度为 0 梁在上表面承受分布荷载 txq 将以 上 2 1 2 3 2 4 2 5 各式代入 2 7 并对其沿高度积分后 2 7 式 右边第一项变为 l xzx T Q k M N b 0 0 1 dxGQQ M N D A M N xx T 2 1 0 02 1 1 2 8 又根据达朗贝尔 Dalembert J be R 原理 将惯性力 i tt u 当作分布力 而忽略质 量体力 则 2 7 式第二项变为 dx h wwuuhb l tttttt 2 3 00 0 12 2 9 又令 x N xz N和 x M是已经给定的沿梁的周边所施加的外力和外力矩 则 2 7 式右边第三项变为 dsMuNwdxtxqb lS 0 3 2 10 根据 Reissner 变分原理 有0 即 0 123 将 2 8 2 9 2 10 式代入上式 并结合上文的推导 可解得用 应力表示的梁的运动 平衡 微分方程 ttxxx tt x xxxx ttxx h QM qhwwNQ huN 3 0 0 0 0 12 2 11 及相应的边界条件 或 或 或 xx xzxxx xx MM wwNwNQ uuNN 0 2 12 上式中 0 u w和 分别为梁中面边界上的已知位移和转角 将 2 11 中的第三式对x求导 则可求出 xx Q 代入 2 11 中的第二式 可得出梁的横向运动控制方程 qhw h wNM ttttx x xxxxx 0 3 0 12 2 13 又将几何方程 2 3 中 xz 的表达式代入本构关系 2 6 中得 xx wCQ 2 14 由 2 11 中的第三式有 ttxxx h MQ 3 0 12 2 15 将 2 15 式代入 2 14 式得 ttxxsx h M C Tw 3 0 12 1 2 16 式 2 16 中 s T为跟踪常数 0 s T表示不考虑横向剪切变形 1 s T表示横向 剪切变形被考虑 将应力表达式 2 5 代入方程 2 11 可得用位移表示的梁的运动微分方程 形式 ttxxsx ttxtt xxxxxxxxxxxxx ttxxxxx C h C D Tw D h w D h D q ww D A wu D A wu D A w u A h wwu 3 0 3 0 0 2 0 0 0 0 0 12 12 2 3 2 17 2 3 方程的无量纲参数化方程的无量纲参数化 下面引入无量刚量 l x h u u 0 h w w l h 0 2 Ah l t 则 0 u l h u x 2 0 u l h u xx w l h wx 2 w l h wxx 1 hl xx 3 w l h w xxx 3 w l h w xxx 4 w l h wxxxx 0 4 2 0 u l Ah u tt 0 4 2 w l Ah wtt 0 3 l A tt 0 4 l A ttx 将以上各式代入 2 17 式并整理得到 2 uwwu 0 2 12 2 3 wq wwwuwuw 2 18 3 2 12 S Tw 其中 D Ah 2 l h 2 Ghl D Gl A D ql q 3 0 EhA 1 12 2 3 Eh D 1 2 E G 现对 2 18 式进行必要的