有限元法技术报告.ppt

由§7—9节的内容，单元上作用集中力的结果： P x1 1 4 3 2 P y 1 P y 4 P x 4 P y3 P x 3 P x 2 P y 2 F x F y F 设单元在Ae内作用体积力，其体力密度为 ，在边界 上作用边界面力，面力密度为 。 Ae内微分体上的体力： 可视为集中力； 上微分段的面力： 可视为集中力； 结点荷载 和 ，由下式给出： 如果只考虑单元的自重，设单元总自重为W，将自重视为作用在重心的集中力： ξ η W (0,0) 重心处的形函数值由式（8-2）得： 5．解答的收敛性 由式（8—1）、(8—2)可知，单元位移为双线性函数： 显然，上述位移模式已经包含了单元的刚体位移状态和常应变状态，及满足收敛性条件（1）（2）。 在单元边界上，x=常数或y=常数，位移是线性变化的，此线性函数可由边界两端位移完全确定。故只要相邻两个单元在公共结点处保持位移相等，则在公共边界上每一点的位移仍保持相等，即满足收敛条件（3）。 §8—3 四结点四边形等参元 1．坐标映射 坐标映射式： 其中 ( , )是母元的形状函数，由式(8—3)所定义，即： 母单元的形函数 与坐标映射式中的 相同， 故称这种单元为等参数单元。 证明坐标映射式成立： 1）在四结点处（如结点1）： 在母单元的结点1，坐标 ，代入 ，由形函数的性质（1）得 ，代入变换式得： ，即坐 标映射式能实现结点1两套坐标 与 之间的一一对应，对其它三结点也可得出同样结论。 1）在四条边上（如12边）： 在母元的12边上： 代入映射式得： 消去参数 得： 上式为过等参元结点1和12边中点的直线，即等参元的12边。其它三条边也具有同样的结果。故映射式成立。 1、母单元与等参元的映射具有如下性质： 母元中的坐标线对应于等参元的直线。 由坐标映射式证明2）知：母元中的某一坐标线 ＝常数对应等参元中的线性函数： 即直线参数方程。即母元中的坐标线对应于等参元中的一条直线。 四结点正方形母元对应于四个结点可以任意布置的直边四边形等参元。 映射式能保证相邻等参元彼此相容。 2．等参元的整体坐标和局部坐标 直角坐标系XY是等参元的整体坐标系。而曲线坐标系 是等参元的局部坐标系，坐标原点位于各单元的形心， 、 的取值范围为（-1，-1），仅在单元内有意义。 整体坐标系XY与单元局部坐标系 的比较： 整体坐标系XY作用于整个结构的所有单元，将在整体分析中使用。局部坐标系 仅适用于单个单元，将在单元分析中使用。 以局部坐标 表示单元，具有统一的形式，且形式简单。以整体坐标XY表示单元，不具有统一的形式，且形式复杂。 由以上比较可知，在单元分析中，采用局部坐标系 较方便，在单元分析中一切运算都在局部坐标系中进行。 3．两类坐标系中导数的变换关系（用于建立位移与应变的关系） 据复合函数的求导法则： 表成矩阵： 其中 称为雅可比(Jacobi)矩阵 其中： 是 的逆阵。 对于四结点四边形等参元， 、 和 的表达式可由坐标映射式求出，式（8-25）~式（8-29）。 4．位移模式 ( ， )即为母单元的形函数，与坐标变换式（8-18）比较，两者的 ( ， )一致，故称等参数单元。 5．单元刚度矩阵 由几何方程求应变： 将位移表达式代入： 其中： 由式（8-22）得[Bi]表达式： 求单元刚度矩阵 单元刚度短阵的分块形式： 现求局部坐标系( ， )中微分面积dA的表示式： 所围成的微分面积ABCD ，即由 造成的 、 。 ，即由 造成的 、 。 微元面积： 故： 由于被积函数 很复杂，需利用高斯求积法求数值积分。 6．等效结点荷载 （1）体积力移置：体积力密度： （2）面力移置：设在边界线上 作用边界面力，面力密度为： 由式（8-16）的第二项： i=1,2,3,4