有限元法基础 教学课件 作者 赵维涛 陈孝珍4 杆件系统有限元法.ppt

《有限元法基础 教学课件 作者 赵维涛 陈孝珍4 杆件系统有限元法.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第四章杆件系统有限元法第四章杆件系统有限元法§4.1平面及空间桁架结构有限元§4.2平面及空间刚架结构有限元1.单元位移模式如图所示的杆单元，横截面面积为，长度为，材料弹性模量为，设单元有两个节点，由于该单元只承受轴向载荷，故节点位移只有轴向位移，节点位移向量为节点力向量4.1.1桁架单元刚度矩阵由于单元上两个节点的位移都只有一个自由度，因此可设单元的位移模式为坐标x的一次函数，即式中：，为待定常数，可根据单元的节点位移确定，即：在节点时，，；在节点时，，，代入上式，得上式可写为形函数表示为4.1.1桁架单元

2、刚度矩阵2.单元应变单元只有轴向应变，将位移代入，得式中：——单元应变矩阵，如下3.单元应力将应变代入物理方程，则4.1.1桁架单元刚度矩阵式中：——单元应力矩阵，为4.单元刚度矩阵单元的单元刚度矩阵为将单元应变矩阵代入上式，可得局部坐标系下桁架单元刚度矩阵为4.1.1桁架单元刚度矩阵由于平面桁架单元中每个节点有两个自由度，因此将上式进一步扩展，可得到在局部平面坐标系xoy下桁架单元刚度矩阵为同理，由于空间桁架单元中每个节点有三个自由度，因此桁架单元在局部空间坐标系下的单元刚度矩阵为4.1.1桁架单元刚度矩阵

3、5.等效节点力若载荷直接作用在节点上则可作为节点力处理，对于其它载荷情况下的基本原理和相应的公式形式同平面问题，这里不再赘述。4.1.1桁架单元刚度矩阵在以上的推导过程中可以发现，局部坐标是根据单元的几何形状选取的，但不同的单元一般具有不同的局部坐标系，因而不能进行不同单元刚度矩阵间的混合运算，即得不到整个计算模型的有限元计算格式。解决这个问题的唯一途径就是建立整体坐标与局部坐标间的坐标变换，并通过这种交换来得到单元刚度矩阵在整体坐标系下的显式。平面桁架单元如图所示。局部坐标系xoy与整体坐标系的夹角为。节点

4、在局部坐标系中的位移与该节点在整体坐标系中的位移间的关系如下4.1.2桁架单元转换矩阵同理，节点在局部坐标系中的位移与节点在整体坐标系中的位移间的关系为4.1.2桁架单元转换矩阵则在整体坐标系中的节点位移用局部坐标系中的节点位移表示为式中，称为平面桁架单元的转换矩阵，如下4.1.2桁架单元转换矩阵在整体坐标系和局部坐标系中节点力之间的转换关系经推导后，有如下关系式经过同样的推导过程，可得空间桁架单元的转换矩阵为4.1.2桁架单元转换矩阵可以证明，转换矩阵的逆矩阵等于它的转置矩阵，即在局部坐标系下单元的平衡方程

5、表示为在整体坐标系下单元的平衡方程为将式转换矩阵代入式，得即得4.1.3整体坐标系下的单元刚度矩阵则桁架单元在整体坐标系下的单元刚度矩阵为将式局部坐标系下的单元刚度矩阵和转换矩阵代入上式，得平面桁架单元在整体坐标系下的单元刚度矩阵为4.1.3整体坐标系下的单元刚度矩阵同理，得到整体坐标系下空间桁架单元刚度矩阵为式中：、和——局部坐标轴在整体坐标系下的方向余弦。4.1.3整体坐标系下的单元刚度矩阵在组装桁架结构总体刚度矩阵时，首先要计算每一个单元在整体坐标系下的刚度矩阵，这里需要强调的是：单元刚度矩阵一定是在整

6、体坐标系下的，而不是局部坐标系下的；然后进行总体刚度的组装，桁架结构总体刚度矩阵的组装过程与第3章中的平面问题（如三角形单元）是一致的。4.1.4总体刚度矩阵用有限元法求解图示结构的桁架内力，设各杆的EA为常数。4.1.5算例解：（1）单元和节点编码如图所示，图中箭头的指向为局部坐标系的正向。4.1.5算例4.1.5算例（2）计算各单元的单元刚度矩阵。单元①，；单元②，；单元③，；单元④，；单元⑤，；单元⑥，，则；（3）形成整体刚度矩阵4.1.5算例（4）形成整体载荷列阵式中：、、、——节点1、2所对应的支座

7、反力。（5）建立整体平衡方程，求解位移分量4.1.5算例由节点1和节点2的约束条件，即，划去上式的1、2、3、4行与1、2、3、4列，得解得节点位移为4.1.5算例（6）求解内力单元在局部坐标系下的内力为，则单元1的内力为同理得4.1.5算例由此可得，单元①为零杆；单元②受拉，轴力为0.56P；单元③受压，轴力为-0.44P；单元④受压，轴力为-0.44P；单元⑤受拉，轴力为0.63P；单元⑥受压，轴力为0.79P。1.无轴向变形的平面梁单元1）单元位移模式如图所示的梁单元为一个无轴向变形的等截面直杆，共有两

8、个节点，节点位移包括挠度和转角，节点力包括剪力和弯矩。4.2.1平面刚架单元单元的节点位移向量表示为节点力向量表示为设单元的挠度的位移模式取为将单元上两个节点的坐标和相应的位移代入上式，可得单元的位移模式为4.2.1平面刚架单元形函数表示为2）单元应变由材料力学可知，梁在发生弯曲变形而引起梁的轴向变形产生的应变称为梁的弯曲应变，可由下式计算4.2.1平面刚架单元将位移代如，可得式中：——单元应变矩阵