05平面问题高精度单元.pptVIP

节点映射条件 坐标变换的形函数 (1) 坐标的映射和变换 节点坐标 坐标变换的形函数 坐标变换的形函数是基准坐标的函数，采用坐标变换的形函数实现物理坐标和基准坐标的变换满足几何条件，即由母单元畸变而得到的各实际单元之间无缝隙、无重叠的连接。 实际单元的位移模式形函数＝坐标变换的形函数，可保证位移在各单元之间连续； 形函数中有常数项和线性项，反映了单元的刚体位移和常量应变，满足了相对完备性，即解的收敛性的必要条件。 等参单元：坐标变换的形函数插值阶次 = 位移模式形函数插值阶次 超参单元：坐标变换的形函数插值阶次 位移模式形函数插值阶次 亚参单元：坐标变换的形函数插值阶次 位移模式形函数插值阶次 坐标变换形函数 位移模式形函数 Jacobi矩阵 Jacobi矩阵的伴随行列式 (2) 形函数对整体坐标的导数 其中： (3) 面积和体积积分变换 d?、d?在物理坐标系中 三维问题，xyz坐标系中d?、d?和d?围成的微六面体体积 （1）几何矩阵和应力矩阵 3 二维等参元的刚度矩阵 平面应力问题 （2）单元刚度矩阵 h——单元厚度 平面应力问题 由于形函数对整体坐标x 和y 的导数中包含雅可比矩阵的逆矩阵，需采用数值积分来求得等参数单元的刚度矩阵。 （3）荷载列阵 当单元在任意一点受集中力Q作用时 当单元受体力P作用时 当单元在某一边界上受面力S作用时 4 八节点曲边四边形单元 母单元 八结点四边形单元 （1）母单元的形函数 （3）位移模式 （2）坐标变换 （4）几何矩阵 （5）单元刚度矩阵 （i，j＝1，2，3，4，5，6，7，8） 5 数值积分 计算单元刚度矩阵和荷载列阵时，要计算 基本思想 函数的积分可近似用若干个点的函数值进行加权相加 n为积分点数，Ak为积分权系数，?k为积分点位置。 如何确定Ak（积分权系数）、?k（积分点位置）？ 一维高斯（GAUSS）积分公式 积分点n 积分点?k 积分权系数Ak 1 0.0000000000 2.000000 2 ?0.577350269 13 ?0.774596669 0.000000000 0.555555556 0.888888888 4 ?0.8611363116 ?0.3399810436 0.347854845 0.652145455 用“高斯点”表示积分点，对称于积分区域的中心成对出现。 对称的高斯点处有相同的权系数。 高斯积分的阶数越大，计算量越大。 n阶高斯积分公式对于2n－1次多项式都精确满足。 积分计算工作量小精度较高的数值积分方法——高斯积分。 二维高斯积分公式 三维高斯积分公式 对不同的积分变量顺序使用一维的高斯积分公式可得到多维情形的高斯积分公式。 4－8节点等参元刚度矩阵数值积分公式 平面4节点四边形等参元的刚度矩阵的计算 试采用4点Gauss积分计算单元的刚度矩阵。材料的弹性模量为E=30×106MPa，泊松比为μ=0.3，厚度为t=0.1m。 坐标的映射函数 雅可比矩阵 1/8 (0.5+0.5?) (2.5+0.5?) 雅可比矩阵的行列式 位移函数关于坐标系的偏导数变换 应变分量关于两个坐标系的计算表达式 (19+4?-?)/64 -J12 -J12 -J12 平面应力单元的弹性系数矩阵 选择4点Gauss积分 单元的刚度矩阵为 本 章 小 结 介绍了平面高精度单元及其构成方法： 矩形单元 高阶单元的构造方法：拉格朗日单元，埃尔米特单元 平面等参元：物理坐标系和基准坐标系的三种映射关系、数值积分 实践课题2 在matlab平台上编写4节点四边形等参元的刚度矩阵计算程序，要求输入单元的节点坐标及材料参数，输出单元刚度矩阵。 k = 1.0e+012 * 1.4619 0.2211 -1.2150 0.0616 -0.3716 -0.4236 0.1248 0.1409 0.2211 1.1491 0.1440 -0.0419 -0.4236 -0.2582 0.0585 -0.8489 -1.2150 0.1440 2.1968 -0.8933 0.5212 -0.0151 -1.5030 0.7645 0.0616 -0.0419 -0.8933 1.8399 0.0673 -0.5250 0.7645 -1.2729 -0.3716 -0.4236 0.5212 0.0673 1.1645 0.2763 -1.3141 0.0800 -0.4236 -0.2582