摘要:设一个平面梁单元有两个节点,如图39.5-2所示。在局部坐标系内,平面梁单元定义有6个自由度: 图39.5一平面梁单元 略去轴向位移,平面梁单元有如下4个自由度 对于平面梁单元,其弯曲变形的位移场 设为下式: 以及 代人节点位移和节点坐标, (其中,I为梁单元的
设一个平面梁单元有两个节点,如图39.5-2所示。在局部坐标系内,平面梁单元定义有6个自由度:
![](/uploads/allimg/140928/10-14092Q34533J7.png)
图39.5一平面梁单元
略去轴向位移,平面梁单元有如下4个自由度
![](/uploads/allimg/140928/10-14092Q34RE04.png)
对于平面梁单元,其弯曲变形的位移场
设为下式:
![](/uploads/allimg/140928/10-14092Q34T2395.png)
![](/uploads/allimg/140928/10-14092Q34T2395.png)
以及
![](/uploads/allimg/140928/10-14092Q34939563.png)
![](/uploads/allimg/140928/10-14092Q34939563.png)
代人节点位移和节点坐标,
![](/uploads/allimg/140928/10-14092Q35004323.png)
![](/uploads/allimg/140928/10-14092Q3501G51.png)
(其中,I为梁单元的长度),得到
![](/uploads/allimg/140928/10-14092Q3525R07.png)
得到用形函数矩阵表达的单元内任一点的位移是
![](/uploads/allimg/140928/10-14092Q35344600.png)
式中,N(x)为平面梁单元的形函数:
为节点位移向量,
。对于
![](/uploads/allimg/140928/10-14092Q3564N42.png)
![](/uploads/allimg/140928/10-14092Q35004323.png)
![](/uploads/allimg/140928/10-14092Q3501G51.png)
(其中,I为梁单元的长度),得到
![](/uploads/allimg/140928/10-14092Q3525R07.png)
得到用形函数矩阵表达的单元内任一点的位移是
![](/uploads/allimg/140928/10-14092Q35344600.png)
式中,N(x)为平面梁单元的形函数:
![](/uploads/allimg/140928/10-14092Q40430958.png)
![](/uploads/allimg/140928/10-14092Q35621W3.png)
![](/uploads/allimg/140928/10-14092Q3564N42.png)
式中N(x)的具体表达式是
![](/uploads/allimg/140928/10-14092Q35I2418.png)
根据最小势能原理导出单元刚度矩阵。弯曲梁的应变能是
![](/uploads/allimg/140928/10-14092Q359413A.png)
式中
![](/uploads/allimg/140928/10-14092Q400151b.png)
其中
![](/uploads/allimg/140928/10-14092Q40054L9.png)
代人粱单元应变能公式,同时假设El对于该单元而言是常童,得单元应变能
![](/uploads/allimg/140928/10-14092Q4051I11.png)
节点位移向量
不是x的函数,上式可以写成
![](/uploads/allimg/140928/10-14092Q40Q2353.png)
![](/uploads/allimg/140928/10-14092Q40545X1.png)
![](/uploads/allimg/140928/10-14092Q40Q2353.png)
应变能的一般形式可以表达成
![](/uploads/allimg/140928/10-14092Q40U2136.png)
式中,
为平面梁单元的单元刚度矩阵,即
![](/uploads/allimg/140928/10-14092Q41031600.png)
![](/uploads/allimg/140928/10-14092Q40Z5Q9.png)
考虑到B是x的函数,上式积分后得局部坐标系下的平面梁单元的单元刚度矩阵
![](/uploads/allimg/140928/10-14092Q4093cA.png)
前面给出的平面单元刚度矩阵是局部坐标系下的表达式,其坐标方向是由单元方向确定的。在各自的局部坐标系下,各个不同方向的梁单元都具有统一形式的单元刚度矩阵。 在组集整体刚度矩阵时,必须建立一个统一的整体坐标系,将所有单元上的节点力、节点位移和单元刚度矩阵都进行坐标变换,变成整体坐标系下的表达式之后,再组集成整体 刚 度矩阵。局部坐标系向整体坐标系的转换关系如下:
设
,分别表示局部坐标系
下的单元节点力(包括等效节点力)、节点位移和单元刚度矩阵,
分别表示整体坐标系0xyz下的单元节点力、节点位移和单元刚度矩阵,T是两种坐
设
![](/uploads/allimg/140928/10-14092Q4145OS.png)
![](/uploads/allimg/140928/10-14092Q4151M08.png)
![](/uploads/allimg/140928/10-14092Q41545c7.png)
标系之间的转换矩阵。两种坐标系下的节点载荷、节点位移和单元刚度矩阵的变换关系为
![](/uploads/allimg/140928/10-14092Q41609549.png)
其中坐标变换矩阵为
![](/uploads/allimg/140928/10-14092Q416312B.png)
式中,θ是
![](/uploads/allimg/140928/10-14092Q41S5518.png)
(责任编辑:laugh521521)
上一篇:杆梁问题的有限元-杆单元
下一篇:杆梁问题的有限元-空间梁单元
文章分享: |
|