中文收敛要求.pptVIP

* * * * Properties of the Shape Functions Property 1. Linear independence Property 2. Delta function properties Property 3. Partitions of unity property Property 4. Linear field reproduction 检查单元行为的一个办法是看其单刚的模态 Example: Q4单元( E = 100 Pa and n = 0.3) Q4单元变形行为检查 例：用 CST 和LST单元求解悬臂梁弯曲问题 分别用3节点CST和6节点LST平面应力单元离散上面的悬臂梁并进行有限元计算，比较端点挠度的计算精度。 Mesh generation: CST Mesh generation: LST sxx stress contours CST LST 两种特殊的单元 1. 无限元(Infinite element) 很多物理问题中用到了无限元，例如: - 航空工程: 飞行器翅膀在空气中运动 - 海洋学(Oceanography): 海水在海岛周围的绕射(diffraction) - 土木工程: 大地提供的支撑载荷 考虑如下问题 P infinite x r 有限元方法： ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 对于静态分析: 可行，单计算成本较高。 对于动态分析: 在边界上会出现发射波，与实际不符。 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? Infinite element a b 更好的方法: 使用无边界的单元，从而减少单元数目和解决边界反射问题。 基本思想: 同时使用两类形函数: - 标准形函数 [ N ] - 递减形函数 [ Nd ] 或递增形函数 [ Ng ] 于是，用映射单元的描述方法: - 把 [ N ] 用于坐标插值 ，把[ Nd ] 用于场量(如位移)插值 - 把 [ N ] 用于场量插值，把 [ Ng ] 用于坐标插值 用一维问题来说明第二种方法 ? ? ? x x=-1 x=0 x=1 ? ? ? 1 2 3 x=x1 x=x2 a a x r 1 2 3 1 解 用标准的形函数进行插值 通常 U3 = 0 (因为它位于无限远处) 或者 U3 = 常数。 列式推导与常规的有限元基本相同，唯一的不同之处在于使用了递增的形函数进行几何映射。 ? ? ? 1 2 3 x=x1 x=x2 a a x r 于是 问题： 如果极点 (r=0) 不属于单元，则不会有问题 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? Infinite element a b 2. 奇异单元(Singular element) 为了反映裂纹尖端处 1/√r 奇异性, 使用 “quarter-node element” ? ? ? ? focused mesh at crack tip ? ? 在裂纹尖端附近 (Mode I fracture) 现用一维等参元进行示例 ? ? 1 2 ? 3 x x=-1 x=0 x=1 ? ? 1 2 ? 3 x x=0 x=L/4 x=L 几何映射 把形函数和 x1 = 0, x2 = L, x3 = L/4 代入几何映射中 位移及应变场为： 解出 但是 于是在x趋近于0的过程中正比于 1/√x。 示例 Exact solution: u(x)=√x u(x) u(x) x x ? 50 equal size 3-node isoparametric elements no quarter-point (singular) element u(x) u(x) x x ? 50 equal size 3-node isoparametric elements with quarter-point (singular) element at x=0 Comparison (N = 5, 20 and 50) 移动最小二乘拟合(Moving Least Squares，MLS) 是x(一维)或者x,y(二维)的单项式： 1D (1, x, x2, x3…) 2D (1, x, y, x2, xy, y2 ….) Note aj are functions of x 定义加权平方误差 其中 W(x-xi) 为权重函数，随x-xi的增加而减小。 q=x/h 选取系数a使得加