空间及轴对称问题有限元的分析.pptVIP

* * * * 注意矩阵黑体。 * 说明平面应力、平面应变的特点，为什么各有三个量为零。指出应力用应变表示是从应变用应力表示的虎刻定律得到的。指出平面应变的物理方程可由置换得到。 哈尔滨工业大学 王焕定 哈尔滨工业大学 王焕定 空间及轴对称问题有限元 概述 空间问题(四面体、六面体类) 轴对称问题 轴对称问题非轴对称荷载 概 述 三个方向尺寸属于同一数量级，所受荷载或形体复杂，不可能像上一章那样简化成平面问题处理，这时必须按空间问题求解。 与平面分析不同，空间有限元分析有如下两个困难： 1）对空间物体进行离散化时不像平面问题那样直观，人工进行离散时很容易产生错误； 2）未知量的数量剧增。 建立网格自动生成前处理程序 采用高阶单元来提高单元精度 平面图形绕面内一轴旋转所产生的空间物体，称为轴对称物体，是一类特殊的空间问题。 空间问题 1 常应变四面体单元形函数 与平面三角形单元相对应，四面体单元内任一点可用“体积坐标”来表示。 各子四面体体积 与三角形单元一样，体积坐标为Ti =Vi /V ,三个是独立的，它有“本1，它0，总和1”的性质。 P 1 2 3 四面体总体积 (右旋体积正) 1 2 3 4 P 2 3 4 P 1 2 4 P 1 3 4 P 剩下来的工作基本和三角形常应变单元类似。作业：自学单元列式内容。 空间问题 2 十结点（二次）四面体单元形函数 类似于平面六结点二次三角形单元，采用试凑法建立结点的形函数。 T1 T2 T3 O 1 2 3 4 5 6 7 8 N1=a×785×234=a×(T1-1/2)×T1 为使N1满足本点为1，可得a=2， 代回后得 N1 =T1 (2T1-1) 余者类似，也可按如下通式得到： 式中p为形函数阶次，分子为不通过i点的平面方程左端项，分母中括号内为i点体积坐标。 请大家自行验证！ 空间问题 3 形成四面体的对角线划分方法 先划分成六面体再分为四面体 1 2 4 3 5 6 8 7 1 4 6 7 1 2 4 6 1 4 3 7 4 8 7 6 1）六面体划分为5个四面体 A5型 1467间连6根对角线 1 5 6 7 空间问题 3 形成四面体的对角线划分方法 1）六面体划分为5个四面体 1 2 4 3 5 6 8 7 1 2 3 5 2 5 6 8 2 4 3 8 3 5 8 7 2 3 5 8 B5型 2358间连6根对角线 相邻六面体必须一个为A5另一个为B5 共同点 相对面对角线 相互空间交叉 空间问题 3 形成四面体的对角线划分方法 2）先划为五面体再划分为6个四面体 1 2 4 3 5 6 8 7 1 2 4 3 5 6 4 3 5 6 8 7 连47、76、63 6874、5673、4763 连23、25、63 2351、3562、3642 A6型 以折面3564分 空间问题 3 形成四面体的对角线划分方法 2）先划为五面体再划分为6个四面体 1 2 4 3 5 6 8 7 连35、52、63 3562、5673、2351 连47、46、63 3764、6874、3642 A6型 以折面2376分 2 4 3 6 8 7 1 2 3 5 6 7 两种A6划分 结果完全相同 空间问题 3 形成四面体的对角线划分方法 2）先划为五面体再划分为6个四面体 1 2 4 3 5 6 8 7 连23、35、45 2453、4753、2351 连45、46、67 4562、5674、6874 B6型 以折面2475分 2 4 5 6 8 7 1 2 4 3 5 7 空间问题 3 形成四面体的对角线划分方法 2）先划为五面体再划分为6个四面体 1 2 4 3 5 6 8 7 连47、76、54 4753、5674、6874 连32、25、54 2351、4352、4562 B6型 以折面3465分 1 2 4 3 5 6 4 3 5 6 8 7 两种B6划分 结果也完全相同 作业：P.95给出了由六面体8个角点点号，按式(4.1.25)求A6和A5型四面体结点号的方法。请考虑B6和B5型的计算公式。 空间问题 4 六面体类单元的形函数 1）八结点单元 1 2 3 4 5 6 7 8 ξ η ζ 类似平面问题矩形线性单元，由试凑法可建立形函数如下： 2）二十结点单元 和平面问题一样，基于试凑法，可以根据上述八结点低阶单元形函数构造各顶点形函数。 ξ η ζ 1 2 3 4 5 6 7 8 9 10 11 12 14 17 20 作业：32结点三次单元 空间问题 5 五面体类单元的形函数 1）试凑法建立六结点形函数 用于与六面体单元联合，解决边界形状不规则物体的分析。 课堂练