现代设计方法杆有限元分析

现代设计方法杆有限元分析第一页，共三十页，2022年，8月28日例：一个受载的悬臂梁和用三角形单元离散化的模型第二页，共三十页，2022年，8月28日注意：1)节点是有限元法的重要概念，有限元模型中，相邻单元的作用通过节点传递，而单元边界不传递力，这是离散结构与实际结构的重大差别；

2）节点力与节点载荷的差别单元：即原始结构离散后，满足一定几何特性和物理特性的最小结构域节点：单元与单元间的连接点。节点力：单元与单元间通过节点的相互作用力。节点载荷：作用于节点上的外载。节点自由度(DOFs)

：用于描述一个物理场的响应特性。第三页，共三十页，2022年，8月28日分离但节点重叠的单元A和B之间没有信息传递（需进行节点合并处理）信息是通过单元之间的公共节点传递的。具有公共节点的单元之间存在信息传递第四页，共三十页，2022年，8月28日非法结构离散不同材料节点不合法第五页，共三十页，2022年，8月28日典型单元类型

单元类型单元形状节点数节点自由度杆单元22梁单元23平面单元32平面四边形42轴对称问题32板壳单元43四面体单元43第六页，共三十页，2022年，8月28日2.单元分析

主要内容：由节点位移求内部任意点的位移，由节点位移求单元应变，应力和节点力。3.整体分析

(1)由节点平衡方程，建立以整体刚度矩阵[K]为系数的，整体节点位移{d}和外载{R}的关系式——整体平衡方程。

(2)考虑几何边界条件，修改总体刚度矩阵，求解全部未知位移分量。第七页，共三十页，2022年，8月28日第八页，共三十页，2022年，8月28日a)受拉阶梯杆示意图①u1u2u3E(1),A(1)②F3123E(2),A(2)l(1)l(2)R1x二.有限元求解基本原理(一维问题)引例:用有限元法求图1所示受拉阶梯杆的位移和应力。已知杆截面面积A(1)=2×10-4m2，A(2)=1×10-4m2，各段杆长l(1)=l(2)=0.1m；材料弹性模量E(1)=E(2)=2×107Pa，作用于杆端的拉力F3=10N。第九页，共三十页，2022年，8月28日1.单元划分根据材料力学的平面假设，等截面受拉杆的同一截面可认为具有相同的位移和应力，即位移只与截面的轴向坐标(x)有关，所以可将阶梯杆看作由两个“一维单元”组成，同一个单元内截面面积及材料特性不变。最简单的情况是，每一个单元有两个节点，他们分别位于单元两端。相邻两单元靠公共节点联结。受拉阶梯杆就简化为由两个一维单元和三个节点构成的有限单元模型。图中①和②是单元号，1,2,3是节点号。取节点位移作为基本未知量，应力由求得的节点位移算出。第十页，共三十页，2022年，8月28日c)单元图b)有限元模型①u1E(1),A(1)Node1Node2l(1)R1u2u2u3②F3Node3E(2),A(2)l(2)Node2uiAeEeNodeiNodejlePiujPjx图5-6②①123第十一页，共三十页，2022年，8月28日2.确定单元插值函数(形函数)有限元法将整个求解域离散为一系列仅靠公共节点联结的单元，而每一个单元本身却视为光滑连续体。单元内任一点的场变量(如位移)可由本单元的节点值根据场变量在单元中的假定分布规律(插值函数)插值求得。本例中，每单元有两个节点，采用线性插值。图c是一典型单元图，两节点分别为i和j，节点场变量值分别记为ui和uj

。设单元中坐标为x处的场变量为u(x)。第十二页，共三十页，2022年，8月28日单元的位移场为u(x)，由两个端点的位移来进行线形插值确定，设u(x)为：

(1.a)单元节点条件：

(1.b)将节点条件(1.b)带入(1.a)，可以求得a0和a1：

(1.c)第十三页，共三十页，2022年，8月28日

则

其中N(x)叫做形状函数矩阵(shapefunctionmatrix),为qe叫做节点位移列阵(nodaldisplacementvector),即(2)形函数矩阵的分量数目应与单元自由度数目相等第十四页，共三十页，2022年，8月28日3.单元方程(单元节点位移与节点力的关系)由等截面杆变形与拉力的关系(虎克定律)得到

(3)式中，

Pi和Pj分别为作用于单元e的节点i和节点j的节点力。第十五页，共三十页，2022年，8月28日式(3)写成矩阵形式为：

(4)或简记为：[k]eqe=Pe

(5)[k]e常称为单元刚度矩阵(stiffnessmatrixofelement)，简称单元刚阵:

P

e={Pi

Pj}T称为单元节点力列阵(nodalforcevector)。式(5)称为单元方程。第十六页，共三十页，2022年，8月28日到目前为止，单元方程(4)或(5)尚不能求解，因为节点力列阵Pe尚属未知。Pe的分量Pi和Pj为相邻单元作用于单元e的节点i和j的力，即属于单元之间的作用力。只有将具有公共节点的单元“组集”在一起才能确定上述节点力和节点外载荷之间的关系。4.单元组集建立总体方程组为获得总体方程组，必须先将单元方程按照局部自由度(ui和uj)和总体自由度(u1、u2和u3)的对应关系进行扩展。第十七页，共三十页，2022年，8月28日①u1E(1),A(1)Node1Node2l(1)R1u2u2u3②F3Node3E(2),A(2)l(2)Node2iijj单元1i=1;j=2单元局部坐标全局坐标单元2i=2;j=3第十八页，共三十页，2022年，8月28日(6)式中，各项上角码表示单元序号；下角码表示自由度总体序号。

具体来说，单元1的扩展方程为：第十九页，共三十页，2022年，8月28日(7)由于相邻两单元公共节点上的基本场变量（位移）相同，所以可将扩展后的各单元方程相加。单元2的扩展方程为：第二十页，共三十页，2022年，8月28日(8)将式(6)和式(7)相加得：第二十一页，共三十页，2022年，8月28日(9)组集后的结果简记为：Kq=P式中，K称为总体特性矩阵(常称为总体刚度矩阵和总刚阵)，P称为总体节点载荷列阵。需指出的是，对单元的一个公共节点而言，除了有相邻单元作用于该节点的力之外，还可能有做用于该节点的外载荷。若一节点上无外载荷作用(如本例中节点2)，则说明各相邻单元作用于该节点的力是平衡的，即该节点的节点合力为零。上述组集过程可记为：有限元模型单元总数第二十二页，共三十页，2022年，8月28日若某节点上有外载荷作用(如节点3)，则各单元作用于该节点的内力和(即方程(8)中第3式左端项的负值)与该节点的外载荷(F3)相平衡，即：

(10)即，列阵F各分量的含义是作用于相应自由度(节点位移)上的节点外载荷。将相应数据代入式(8)得：

(11)第二十三页，共三十页，2022年，8月28日

上式即为本题的总体线性代数方程组，但不能获得唯一解，因为上式中的矩阵是奇异的。这种奇异性不是因数据巧合造成的，而是有其必然性。原因在于总体方程组式(8)只考虑了力平衡条件，而只根据力平衡不能唯一地确定系统的位移，因为系统在有任意刚性位移的情况下仍可处于力平衡状态。为获得各节点位移的唯一解，必须消除可能产生的刚体位移，即必须计入位移边界条件。第二十四页，共三十页，2022年，8月28日本题的位移边界条件为u1=0,那么，式(11)中只剩下两个待求的自由度u2和u3。也就是说，可从式(11)中消去一个方程。譬如，舍去第一个方程并将u1=0代入后得：

(12)解得：

u2=2.5×10-4m；u3=7.5×10-4m。

q=

[u1u2u3]T=[02.5×10-47.5×10-4]Tm.这与材料力学求得的结果相同。5.计入边界条件，解方程组第二十五页，共三十页，2022年，8月28日应变的表达由几何方程得知，1D单元中任一点的应变

(13)其中

(14)B(x)称为单元应变矩阵，或称为几何函数矩阵(strain-displacementmatrix).6.计算单元应变和应力第二十六页，共三十页，2022年，8月28日(15)其中(16)S(x)叫做应力矩阵

(stress-displacementmatrix).应力的表达第二十七页，共三十页，2022年，8月28日对于单元1对于单元2第二十八页，共三十页，2022年，8月28日7.求支反力具体对单元①，有

(17)其中R1为节点1的外力，即为支反力，P2为单元①的节点2所受的力，将u1和u2的值带入式(17)，有

第二十九页，共三十页，2022年，8月28日作业用有限元法求图示受拉阶梯杆的位移和应力。已知杆截面面积A(1)=4×10-4m2，A(2)=2×10-4m2，，A(3)