有限元 2-弹性力学平面问题(2.4矩形单元,2.5六节点三角形单元).ppt

第2章 弹性力学平面问题有限单元法 2.1 三角形单元 2.2 三角形单元中几个问题的讨论 2.3 平面问题有限元程序设计? 2.4 矩形单元 2.5 六结点三角形单元 2.6 四结点四边形单元 2.7 八结点曲线四边形等参元? 2.8 几个问题的补充 图示深梁，设划分成具有八个矩形单元的网格，从中取出一典形单元e，单元有４个结点(i,j,m,p)。 一、局部坐标 设单元的边界平行于结构整体坐标x, y，边长分别为2a, 2b， 结点编号从左下角开始逆时钟编成1，2, 3, 4。 单元局部坐标的原点取在单元的形心(x0, y0)上，并采用无因次坐标ξ,η。ξ,η轴分别平行于x, y轴， 二、结点位移列阵和结点力列阵 每个结点2个位移分量，共８个位移分量，设结点位移和结点力列阵分别为： 位移函数是单值连续的。在平行于x轴的直线上，位移分量是x的线性函数；在平行于y轴的直线上，它是y的线性函数，故2-4-4称为坐标x, y的双线性函数。 将各结点坐标代入位移模式 以矩阵形式表示： 如果引进参数: ξ0=ξiξ, η0=ηiη(i=1, 2, 3, 4), (ξi, ηi)是矩形单元4个结点的局部坐标。结点i(ξi, ηi)的坐标值分别是 (-1,-1), (1,-1),(1,1), (-1,-1)。代入上式，则可将上式简记成： 矩形单元的形函数，具有与三角形单元类似的特点，即： 在结点上的形函数值： 四、几何矩阵[B]，应力矩阵[S] 在三角形单元中，我们已经推知： 应变： 应力： 矩形单元自然具有相同关系，只是[B]， [S]的内容有所区别。 仿照三角形单元： 式中: 单元应力矩阵： 式中： 五、单元刚度矩阵 将上述公式代入单刚的一般表达式 第2章 弹性力学平面问题有限单元法 2.1 三角形单元 2.2 三角形单元中几个问题的讨论 2.3 平面问题有限元程序设计? 2.4 矩形单元 2.5 六结点三角形单元 2.6 四结点四边形单元 2.7 八结点曲线四边形等参元? 2.8 几个问题的补充 2.5 六结点三角形单元 一、面积坐标 1．面积坐标的定义 在第2节中，我们讨论过三角形单元形函数的几何意义，如果从三角形单元e内任意点p向三个顶点i,j,m 引连线，使其分割成三个小三角形便可得到这些小三角形面积与形函数的关系。 反过来说，只要给定了三角形面积 Ai , Aj , Am，p点即可确定。因为 Ai, Aj ,Am 中只有两个是独立的(Ai＋Aj＋Am=A)，这与P点坐标(x, y)对应。因此，我们可以用三角形面积来定义三角形内的任意一点P的位置，并称它为面积坐标。P点的面积坐标定义为以下三个面积比: 由此定义,必有：Li +Lj +Lm =1。由于面积坐标Li , Lj , Lm 只限于用来确定三角形单元上(内)点的位置（只有当L(i,j,m)≤1才有定义）,并且只是在单元分析中使用，所以L(i,j,m)属于局部坐标。 直角坐标与坐标原点、平移转动等有关，可因人而异，故称它为人为坐标，而面积坐标具有不变性，故称其为自然坐标。） 2．面积坐标Li和整体坐标x,y之间的关系 三角形内部点和边界点上的面积坐标值可以用右图中的一些等值线表示（面积坐标的等值线）。从图中可看出：在平行于jm边的任何一条直线上各点都具有相同的Li 值(可从几何作图上看出)其值等于该直线到对应jm边界的距离与结点i到jm边界距离的比值。其它亦然。 3．面积坐标与直角坐标的关系 将式(2-5-1)写成矩阵形式得到用x、y表示的面积坐标: 其反变换式为: 数学意义:三角形ijm内某点的坐标(x,y),可以用三角形顶点(3个结点)的坐标(xi , yi ), (xj , yj ), (xm , ym )进行函数插值而获得，其插值多项式是三角形的面积坐标Li , Lj , Lm， 也就是三角形的形函数Ni , Nj , Nm 。 等参单元 “等参元”的另一种定义是：（参考书：湖南科技社《有限元法概论》 ） 设位移变量u的近似式: 亚参元概念： 有时为了减少结点坐标数据的输入量，用于位移描述的结点数和用于坐标变换的结点数可以不一样。如图示四边形单元，可将位移取为8个结点,即位移函数为： 而坐标变换式则取为: 我们称这种坐标结点少，而位移结点多的单元为亚参元。 二、六结点三角形单元 在三角形单元三条边界的中点上各增设一个结点，使三角形单元共有12个位移自由度。 １、结点位移和节点力列阵： ２、单元位移函数和形函数 设为２次多项式： 为了计算简单，采用面积坐标Li, Lj