文档介绍：有限元分析及应用
第1页，共65页，编辑于2022年，星期六
*
:
重力方向总是竖直向下的
离心力总是沿径向向外的
没有一种材料能抵抗 1,000,000 psi 的应力
轴对称的物体几乎没有为
察看能量误差：Plot Results > Element Solu > Error Estimation > Energy error (ENER).
能量误差
第12页，共65页，编辑于2022年，星期六
*
应力上下限可以确定由于网格离散误差对模型的应力最大值的影响.
显示或列出的应力上下限包括:
估计的上限 - SMXB
估计的下限 - SMNB
应力上下限限并不是估计实际的最高或最小应力。它定义了一个确信范围。如果没有其他的确凿的验证，就不能认为实际的最大应力低于 SMXB.
例如:SMX=32750是节点解的实际值
SMXB=33200是估计的上限
X stress SMAX ~ 32,750 psi
SMXB ~ 33,200 psi
(difference ~ 450 psi ~ %)
应力上下限
第13页，共65页，编辑于2022年，星期六
局部的细化
采用plane42单元网格局部细化与未细化
能量百分比误差
局部细化
未细化
Displacement
DMX=-03
SEPC=
DMX=-03
应力偏差
Element Solution(SDSG)
SMN=
SMX=
SDSG
SMN=
SMX=
能量误差估计
Element Solution(SERR)
SMN=-03
SMX=
SERR
SMN=
SMX=
应力上下限
Nodal solution(SEQV)
SMN=
SMNB=
SMX=4579
SMXB=4623
SEQV
SMN=
SMNB=
SMX=4421
SMXB=4999
第14页，共65页，编辑于2022年，星期六
*
P方法及p单元的应用
P方法应用控制：
P方法用于线弹性结构分析—实体和壳体。
P单元由以下5种单元：
2-D Quadrilateral (Plane145)
2-D Triangle (Plane146)
3-D Brick (Solid 147)
3-D Tetrehedron (Solid 148)
3-D Shell (Solid 150)
规定 % 局部应力差，使用p方法计算的最大X方向应力约为 34,700 psi
(比普通h方法高出大约 5% )
P 单元的位移形函数
u=a1+a2x+a3y+a4x2+a5xy+a6y2
v=a7+a8x+a9y+ a10x2+a11xy+a12y2
P方法的优点：
如果使用 p-方法进行结构分析，可以依靠p单元自动调整单元多项式阶数（2-8），达到收敛到设定的精度. 对这种方法的相信程度，与使用经验有关.
第15页，共65页，编辑于2022年，星期六

GUI:Main Menu > Preference > P-Method
定义一个P单元，P方法被激活。

建模过程与H-单元分析相同，单元类型必须用P单元
（a）指定P单元水平
定义局部P-水平等级
定义P单元时用KeyOpt选项定义
定义整体p-水平等级
命令： PPRANGE , START, MAX
GUI: Main Menu > Solution > P-Method > Set P Range
(b)定义几何模型应用实体建模
(c) 用P单元分网。自适应网格对P方法是无效的
、求解
应用实体模型加载，而不是有限元模型
求解：推荐采用条件共轭梯度法（PCG），但PCG对于壳体P单元无效
察看结果
P方法进行静力分析的步骤
第16页，共65页，编辑于2022年，星期六
举例：
E=30e6 lb/in2
V=
Thick=
在节点（0，5，0）处的收敛标准设为1%
第17页，共65页，编辑于2022年，星期六
高级网格划分技术
延伸网格划分
映射网格划分
层状网格划分
第18页，共65页，编辑于2022年，星期六
延伸网格划分 & 举例
将一个二