有限元应用第4章等参元.pptx

第4章 等参数单元本章主要内容§4.1 等参数单元的意义及来由§4.2 平面问题四边形等参单元§4.3 三维问题六面体等参单元§4.4 高斯数值积分法§4.5 高斯数值积分法应用广西大学土木建筑工程学院§4.1 等参数单元的意义及来由实际问题常常需要使用一些几何形状不太规整的单元来逼近原问题。直接研究这些不规整单元的表达式比较困难（在整体坐标系下构造单元形函数、刚度矩阵及等效节点载荷列阵都十分冗繁）。事实上，形状不规整的单元和形状规整的单元（正方形单元、正六面体单元）可以建立一种映射关系，使得物理坐标系中的整体坐标和自然坐标系中的局部坐标一一对应。利用形状规整的单元的插值函数可以描述不规则单元的力学性能，于是就产生了等参元。广西大学土木建筑工程学院§4.1 等参数单元的意义及来由坐标映射四边形正方形子单元母单元单元局部（自然坐标）整体坐标（物理坐标）映射坐标映射几何坐标：位移插值：广西大学土木建筑工程学院§4.1 等参数单元的意义及来由等参数单元：采用同样的节点、形函数描述单元的几何坐标和位移的单元，也称为等参元。等参元的插值函数用自然坐标给出。如果坐标变换节点数多于位移插值的节点数，称为超参变换。反之，如果坐标变换节点数少于位移插值的节点数，则称为亚参变换。等参单元的提出为有限元法成为现代工程领域最有效的数值分析方法迈出了极为重要的一步。利用等参元可以构造高精度单元。等参单元定义的给出广西大学土木建筑工程学院§4.2 平面问题四边形等参单元映射1. 几何坐标变换：构造插值函数:节点几何坐标条件：广西大学土木建筑工程学院§4.2 平面问题四边形等参单元1. 几何坐标变换：单元形函数：即：广西大学土木建筑工程学院§4.2 平面问题四边形等参单元1. 几何坐标变换：子单元结点坐标子单元几何坐标单元形函数矩阵几何坐标的变换：子单元的几何坐标为母单元局部坐标的函数。广西大学土木建筑工程学院§4.2 平面问题四边形等参单元2. 位移插值函数：位移变换位移插值函数：节点位移条件：广西大学土木建筑工程学院§4.2 平面问题四边形等参单元2. 位移插值函数：同理可得：广西大学土木建筑工程学院§4.2 平面问题四边形等参单元2. 位移插值函数：形函数：单元形函数矩阵单元位移单元结点位移向量简写成：广西大学土木建筑工程学院§4.2 平面问题四边形等参单元等参元：单元的几何坐标、位移都采用同样的节点和相同的形函数通过插值的方式表示。形函数用自然坐标给出。坐标：位移：形函数：广西大学土木建筑工程学院§4.2 平面问题四边形等参单元3. 单元应变：单元应变矩阵：问题的关键：是局部坐标的直接函数，可直接对自然坐标求导。是整体坐标(x,y)的复合函数，可利用复合函数对x,y求导。广西大学土木建筑工程学院§4.2 平面问题四边形等参单元3. 单元应变：形函数：自然坐标的偏导数：几何坐标：几何坐标的偏导数：雅可比矩阵：广西大学土木建筑工程学院§4.2 平面问题四边形等参单元3. 单元应变：利用复合函数求偏导数：伴随矩阵：逆矩阵：雅可比矩阵：广西大学土木建筑工程学院§4.2 平面问题四边形等参单元3. 单元应变：逆矩阵：都是自然坐标的函数。代入可求得节点的应变矩阵：则单元任意点的应变：广西大学土木建筑工程学院§4.2 平面问题四边形等参单元4. 单元应力： 四边形等参元的应力： 单元应力矩阵： 对于平面问题，节点的应力矩阵：广西大学土木建筑工程学院§4.2 平面问题四边形等参单元5. 单元刚度矩阵：Ke可按节点分块 :刚度子块 :(r, s = 1~4) 令 :则 :广西大学土木建筑工程学院§4.2 平面问题四边形等参单元5. 单元刚度矩阵：令 :则刚度子块 :(r, s = 1~4) 广西大学土木建筑工程学院§4.2 平面问题四边形等参单元6. 单元等效节点荷载向量： 1. 集中力的等效节点力: 由于计算复杂，把有集中力处设置为节点，转化为直接节点荷载。 2. 体力的等效节点力:展开后单元等效到节点上的力：（i=1~4）广西大学土木建筑工程学院§4.2 平面问题四边形等参单元6. 单元等效节点荷载向量： 3. 面力的等效节点力:单元等效到节点上的力：如果是边界为直边，形函数在边界上仍为线性函数，如果是均布荷载作用，与三节点常应变单元处理相同，只需把合力平均分配到此边界的两个节点上去即可。广西大学土木建筑工程学院§4.2 平面问题四边形等参单元6. 单元等效节点荷载向量： 3. 面力的等效节点力:如果是边界为曲边，任意分布正应力σ和剪应力τ ，由于表面力的运算是曲线积分，通常把表面的正应力和切应力沿坐标轴分解为 qx ,qy广西大学土木建筑工程学院§4.2 平面问题四边形等参单元6. 单元等效节点荷载向量： 3. 面力的等效节点力:( 12, 3