通用后处理和时间历程后处理的区别
发布网友
发布时间:2022-04-29 15:29
我来回答
共2个回答
热心网友
时间:2023-11-06 14:47
APDL
APDL换行与续行-
APDL规定每行72个字符
如果要写表达式A=C1+C2 (C1与C2都为表达式
可以用
B=C1
A=B+C2
将一行拆成两行来做但是如果不是表达式,而是输入一个命令参数过多的话,可以用续行命令RMORE,格式如下:
RMORE, R7, R8, R9, R10, R11, R12
这个命令每次也只能输入6个参数,如果多于6个,可以重复使用RMORE就可以输入13-18,19-24等等。另外,于上面续行相应的是换行,一行命令太短可以使用多个命令共一行
$”,没有双引号。这样就可以将一行变成多行使。:) 换行符是“
ANSYS常见术语/命令对照表
Utility Menu 实用菜单
SAVE_DB 存储数据库 RESUME_DB 恢复数据库
Select Entity 选择实体 Comp/Assembly 组元/集合
Plot/Replot 画图/重新画图 Pan,Zoom,Rotate… 平移,缩放,旋转…
WorkPlane(WP) 工作平面 Coordinate System(CS) 坐标系
Macro 宏 Preference… 优先设置…
Preprocessor 前处理 General Postproc 通用后处理
TimeHist Postproc 时间历程后处理 APDL ANSYS参数化设计语言
Line Fillet 在两条线的过渡生成线 Arbitrary 任意形状
Cylinder 圆柱体 Prism 棱柱体 Cone 圆锥形
Sphere 球形 Polygon 多边形 Stress 应力
Strain 应变 Displacement 位移 DOF 自由度
Von Mises(Stress) 平均应力 Contour 等高线(图)
Deformed/Undeformed shape 变形后/未变形的形状 Results Summary 结果摘要
Radiation Matrix 辐射矩阵 Modeling 建模
Meshing 划分网格 Attribute 属性 LS (Load Step) 载荷步
ansys的常用命令介绍
对ANSYS学习也有一个来月的时间了,可是还是什么都不会!郁闷!整理了一些ANSYS常用的命令;但深知自己的水平,还不敢保证完全正确;给大家一些参考,望指正:
1. A,P1,P2,…,P17,P18(以点定义面)
2. AADD,NA1,NA2,…NA8,NA9(面相加)
3. AATT,MAT,REAL,TYPE,ESYS,SECN(指定面的单元属性)
【注】ESYS为坐标系统号、SECN为截面类型号。
4. *ABBR,Abbr,String(定义一个缩略词)
5. ABBRES,Lab,Fname,Ext(从文件中读取缩略词)
6. ABBSAVE,Lab,Fname,Ext(将当前定义的缩略词写入文件)
7. ABS,IR,IA,--,--,Name,--,--,FACTA(取绝对值)
8. ACCAT,NA1,NA2(连接面)
9. ACEL,ACEX,ACEY,ACEZ(定义结构的线性加速度)
10. ACLEAR,NA1,NA2,NINC(清除面单元网格)
11. ADAMS,NMODES,KSTRESS,KSHELL
12. ADAPT, NSOLN, STARGT, TTARGT, FACMN, FACMX, KYKPS, KYMAC
13. ADD,IR, IA, IB, IC, Name, --,-- , FACTA, FACTB, FACTC(变量加运算)
14. ADELE,NA1,NA2,NINC,KSWP(删除面)
【注】KSWP =0删除面但保留面上关键点、1删除面及面上关键点。
15. ADRAG,NL1,NL2,…,NL6,NLP1,NLP2,…,NLP6(将既有线沿一定路径拖拉成面)
16. AESIZE,ANUM,SIZE(指定面上划分单元大小)
17. AFILLT,NA1,NA1,RAD(两面之间生成倒角面)
18. AFSURF,SAREA,TLINE(在既有面单元上生成重叠的表面单元)
19. *AFUN, Lab(指定参数表达式中角度单位)
20. AGEN, ITIME, NA1, NA2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE(复制面)
21. AGLUE,NA1,NA2,…,NA8,NA9(面间相互粘接)
22. AINA,NA1,NA2,…,NA8,NA9(被选面的交集)
23. AINP,NA1,NA2,…,NA8,NA9(面集两两相交)
24. AINV,NA,NV(面体相交)
25. AL,L1,L2,…,L9,L10(以线定义面)
26. ALIST,NA1,NA2,NINC,Lab(列表显示面的信息)
【注】Lab=HPT时,显示面上硬点信息,默认为空。
27. ALLSEL,LabT,Entity(选择所有实体)
【注】LabT=ALL(指定实体及其所有下层实体)、BELOW(指定实体及其下一层实体);
Entity=ALL、VOLU、AREA、LINE、KP、ELEM、NODE。
28. AMESH,NA1,NA2,NINC(划分面生成面单元)
AMESH,AREA,KP1,KP2,KP3,KP4(通过点划分面单元)
29. /AN3D,Kywrd,KEY(三维注释)
30. ANCNTR,NFRAM,DELAY,NCYCL(在POST1中生成结构变形梯度线的动画)
31. ANCUT,NFRAM,DELAY,NCYCL,QOFF,KTOP,TOPOFF,NODE1,NODE2,NODE3(在POST1中生成等势切面云图动画)
32. ANDATA,DELAY,NCYCL,RSLTDAT,MIN,MAX,INCR,FRCLST,AUTOCNTRKY(生成某一范围内的结果数据的顺序梯度线动画)
33. ANDSCL,NFRAM,DELAY,NCYCL(在POST1中生成结构变形的动画)
34. ANFLOW,NFRAM,DELAY,NCYCL,TIME,SPACING,SIZE,LENGTH(生成粒子流或带电粒子运动的动画)
35. /ANGLE,WN,THETA,Axis,KINCR(绕指定轴旋转视图)
36. ANHARM,NFRAM,DELAY,NCYCL(生成谐波分析的动画)
37. ANIM,NCYCL,KCYCL,DELAY(动画显示图形序列)
38. ANISOS,NFRAM,DELAY,NCYCL(在POST1中生成等势面云图动画)
39. ANMODE,NFRAM,DELAY,NCYCL,KACCEL(在POST1中生成结构变形模态的动画)
40. /ANNOT,Lab,VAL1,VAL2(激活图形显示注释)
【注】Lab=OFF、ON、DELE、SAVE、SCALE、XORIG、YORIG、SNAP、STAT、DEFA、REFR、TMODE。
41. ANORM,ANUM,NOEFLIP(重新定义面的法线方向)
【注】NOEFLIP=0、1。
42. ANTIME,NFRAM,DELAY,NCYCL,AUTOCNTRKY,RSLTDAT,MIN,MAX(在指定时间段内生成动画)
详细见网站:www.101315.com
有限元基本概念和原理
有限元分析(FEA,Finite Element Analysis)的基本概念是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成,对每一单元假定一个合适的 (较简单的)近似解,然后推导求解这个域总的满足条件(如结构的平衡条件),从而得到问题的解。这个解不是准确解,而是近似解,因为实际问题被较简单的问 题所代替。由于大多数实际问题难以得到准确解,而有限元不仅计算精度高,而且能适应各种复杂形状,因而成为行之有效的工程分析手段。
有限元是那些集合在一起能够表示实际连续域的离散单元。有限元的概念早在几个世纪前就已产生并得到了应用,例如用多边形(有限个直线单元)*近圆来求得圆 的周长,但作为一种方法而被提出,则是最近的事。有限元法最初被称为矩阵近似方法,应用于航空器的结构强度计算,并由于其方便性、实用性和有效性而引起从 事力学研究的科学家的浓厚兴趣。经过短短数十年的努力,随着计算机技术的快速发展和普及,有限元方法迅速从结构工程强度分析计算扩展到几乎所有的科学技术 领域,成为一种丰富多彩、应用广泛并且实用高效的数值分析方法。
有限元方法与其他求解边值问题近似方法的根本区别在于它的近似性仅限于相对小的子域中。20世纪60年代初首次提出结构力学计算有限元概念的克拉夫(Clough)教授形象地将其描绘为:“有限元法=Rayleigh Ritz法+分片函数”,即有限元法是Rayleigh Ritz法的一种局部化情况。不同于求解(往往是困难的)满足整个定义域边界条件的允许函数的Rayleigh Ritz法,有限元法将函数定义在简单几何形状(如二维问题中的三角形或任意四边形)的单元域上(分片函数),且不考虑整个定义域的复杂边界条件,这是有限元法优于其他近似方法的原因之一。
对于不同物理性质和数学模型的问题,有限元求解法的基本步骤是相同的,只是具体公式推导和运算求解不同。有限元求解问题的基本步骤通常为:
第一步:问题及求解域定义:根据实际问题近似确定求解域的物理性质和几何区域。
第二步:求解域离散化:将求解域近似为具有不同有限大小和形状且彼此相连的有限个单元组成的离散域,习惯上称为有限元网络划分。显然单元越小(网络越细)则离散域的近似程度越好,计算结果也越精确,但计算量及误差都将增大,因此求解域的离散化是有限元法的核心技术之一。
第三步:确定状态变量及控制方法:一个具体的物理问题通常可以用一组包含问题状态变量边界条件的微分方程式表示,为适合有限元求解,通常将微分方程化为等价的泛函形式。
第四步:单元推导:对单元构造一个适合的近似解,即推导有限单元的列式,其中包括选择合理的单元坐标系,建立单元试函数,以某种方法给出单元各状态变量的离散关系,从而形成单元矩阵(结构力学中称刚度阵或柔度阵)。
为保证问题求解的收敛性,单元推导有许多原则要遵循。 对工程应用而言,重要的是应注意每一种单元的解题性能与约束。例如,单元形状应以规则为好,畸形时不仅精度低,而且有缺秩的危险,将导致无法求解。
第五步:总装求解:将单元总装形成离散域的总矩阵方程(联合方程组),反映对近似求解域的离散域的要求,即单元函数的连续性要满足一定的连续条件。总装是在相邻单元结点进行,状态变量及其导数(可能的话)连续性建立在结点处。
第六步:联立方程组求解和结果解释:有限元法最终导致联立方程组。联立方程组的求解可用直接法、选代法和随机法。求解结果是单元结点处状态变量的近似值。对于计算结果的质量,将通过与设计准则提供的允许值比较来评价并确定是否需要重复计算。
简言之,有限元分析可分成三个阶段,前处理、处理和后处理。前处理是建立有限元模型,完成单元网格划分;后处理则是采集处理分析结果,使用户能简便提取信息,了解计算结果。
对称边界条件
有些结构由於具有某些对称性,我们可以在对称面上施加适当的对称边界条件,这样只要建立部分的模型,既省时又省力。
这里讲的对称,不只是几何形状的对称,还包括边界条件、外力施加、材料性质的对称。如果仅有形状对称,其他条件只要有一样不对称,就不能用对称模型求解。
举一个简单的例子。图一是一个含圆孔的平板,边界条件如图所示。我们知道图一的模型上下、左右对称,因此可以简化成图二。
为什麽图二的边界条件可以设为如此?为回答这个问题,我们先将图一分为上下两部分。假设A点与B点是上下两部分相对称的点(图三),则A点与B点的「y」 方向位移方向相反,大小相等。想像A与B点越来越靠近对称面,直到几乎重合在一起(图四),此时A、B两点仍可视为对称於对称面,两点的y方向位移仍为大 小相等,方向相反。然而此时两点几乎重合,因此A、B两点的位移应该相同。为了满足上面的条件,唯一的可能就是在对称面上y方向位移等於零。图二中的另一 边x方向位移等於零,可用相同方法解释。
事实上,图二的两个对称边界,还有一个条件,那就是剪应力xy为零。这个道理的解释方法与位移类似。参考图三,A、B两点由於对称,剪应力也要一模一样; 当A、B两点越来越靠近对称面(图四),两点的剪应力就必须满足牛顿第三运动定律(作用力等於反作用力,方向相反),此时唯一的可能就是剪应力为零。因此 在对称面上,剪应力为零
总之,在对称面上,垂直於对称面的位移以及作用於对称面上的剪应力皆为
在破坏力学的应用
破坏力学是固体力学的一个分支,这门学科主要在探讨四个部分:
1. 含裂纹结构的受到外力时的应力分布;
2. 含裂纹结构受到多大的外力,裂纹会成长;
3. 结构中,裂纹一旦成长,会往那个方向成长;
4. 各种工程结构抵抗裂纹成长的能力,这部分通常由实验决定。
首先先介绍破坏力学的基本概念。图一为含裂纹结构的破坏模式,共分为第一型(mode I or opening mode)、第二型(mode II or sliding mode)以及第三型(mode III or tearing mode)。
裂纹尖端为奇异点(应力正比於根号 r分之一(stress ~ 1/r^0.5),r为结构中某一点与裂纹尖端的距离),去探讨裂纹尖端(或附近)的应力有多大是没有意义的。有鉴於此,我们引入了一个参数:应力强度因子(stress intensity factor)。应力强度因子共有三个,一般写做KI、KII以及KIII,分别对应到三种不同的破坏模式。
在线弹性破坏力学(linear elastic fracture mechanics,简称LEFM)中,外力的大小与应力强度因子成正比关系。此外,应力强度因子也和几何参数(例如裂纹长度、外力与裂纹的距离等)有 关。对於同一结构,应力强度因子可视为含裂纹结构所受外力大小的一个指标;换言之,应力强度因子越大,结构越危险。当应力强度因子超过「破坏韧性 (fracture toughness),记为Kc」,裂纹开始成长。破坏韧性通常由制作标准试片由实验求得,不同的材料有不同的破坏韧性,所以破坏韧性可视为材料性质。
至於裂缝会往那个方向成长?有好几个理论可以预测,例如能量释放率理论(energy release rate)、最大周向应力理论(maximum circumferential stress)、J积分理论(J-integral)、应变能密度理论(strain energy density theory)。以上几个理论,除了最大周向应力之外,基本上都和应变能有关,这个部分较为深入,不详述了。
接着我们要探讨有限元素法在破坏力学的应用。前面讲过,探讨裂缝尖端附近的应力有多大是没有意义的。因此我们的重点在於,如何运用有限元素法求得裂纹尖端的应力强度因子。最常见的就是利用「四分之一节点(quarter point)」元素来模拟裂纹尖端,如图二。严格来说,四分之一节点元素不是「三角形元素」,它是由四边形元素退化而成,退化方式如图三。一般我们利用节点上的位移,求得裂纹尖端的开口位移(crack tip opening displacement)来反推得到应力强度因子。MARC及ANSYS皆用四分之一节点元素来算应力强度因子。
以上所说的四分之一节点元素,其最原始的统御方程式和弹性力学并没有两样,只是将其中几个节点变了位置,因此形状函数也跟着改变,最後得到一个结果:这种元素内部的应力呈现根号r分之一的奇异性,因此可用来模拟裂纹。
另外有些软体,例如NASTRAN,则利用另一种特殊元素,共有18个节点,如图四。这种元素最原始的统御方程式,和一般的弹性力学稍有不同,它一开始就假设应力正比於根号 r分之一。
不管如何,裂纹尖端的元素与一般的元素不同,裂纹尖端网格的大小,由算出来的应力强度因子准不准来决定。如何知道准不准?可建立一个简单模型,有理论解的 题目,即可比较,例如图四。以平面问题为例,若采用四分之一节点元素,元素边长大约在裂纹长度的百分之六左右;若采用NASTRAN的特殊元素,其长度更 只要裂纹长度的百分之十至二十即可。当然以上的准则并非万用,不同的题目最好还是多试几个网格粗细,再做决定。
方向
元素有方向性。有些元素如果方向搞错,则跑不出来;有些虽然跑得出来,但结果却有问题。 以平面元素(包括平面应变、平面应力以及轴对称)为例,四个节点的编号必须是「逆时钟方向」。如果是顺时钟方向,则在有限元素的定义中,这个元素的面积小於零,就跑不出来了。
板(或壳)元素也有方向性。板元素的方向由右手定则决定,也就是元素的正向方向,由元素编号顺序的方向决定,如此也定义了板元素的上表面或下表面。当板受 到bending时,上表面和下表面分别会受到张应力以及压应力(当然也可能颠倒,视bending的方向而定)。图一为两个相邻的板元素受到 bending作用,元素1的方向和元素2的方向不同,因此它们的上、下表面也不一致(参考图一)。在受到图一所示的bending时,元素1的上表面为 张应力,下表面为压应力;元素2的上表面则为压应力,下表面为张应力。有些有限元素软体(如MARC)在计算板元素的「上表面」节点应力时,是根据两个元 素的「上表面」高斯点上应力经外差到节点上再平均求得,但注意图一中两个元素的上下表面定义颠倒,所以平均後的上表面节点应力就很不准了。
再解释清楚一点,我们先将两个元素分开,同时画上厚度,比较清楚,如图二,图一中的A点在图二则可分为A1至A4四点。如果单独看元素1,A1点应该受到 张应力,A2为压应力;单独看元素2,A3为张应力,A4为压应力。当有限元素法在计算A点「上表面」应力时,是将元素1以及元素2的「上表面」应力平均 (即将A1及A4的应力平均),本来应该上表面是张应力的,被这麽一搞,应力下降,误差就来了。
一般有限元素法软体,在做automesh时,元素的方向与surface(或area)的方向一致,因此同一个surface(或area)做出来的网 格,方向都相同。但如果两个surface(或area)相邻,就要特别注意方向有没有一致。如果是「手动」建立网格,更要特别注意建立出来的元素的方 向,以免造成困扰。
形状函数
形状函数在有限元素法中是非常重要的一个概念,它定义了元素内部位移的分布。以一维线性元素为例(图一) 此种元素有两个节点,分别以节点1及2表示。节点1的座标为xi = -1,节点2则为xi = -2。元素中有几个节点,就有几个形状函数。因此我们有两个形状函数。形状函数有个特点:考虑第n个形状函数,若代入第n个节点的座标,其函数值为1;若 代入其他节点的座标,则函数值为0。由於图一为两节点元素,所以形状函数为线性,其表示式为:其他种类的元素的形状函数,一般有限元素法书籍(或商用有限 元素法软体的使用手册)皆有提到,这里不多写了。
元素的特性主要由形状函数所主宰。例如三节点之三角形元素,沿着元素内部任一方向,位移皆为线性分布。为什麽我们知道它是线性分布?只要查查该元素的形状 函数立刻就明白了。由於应变为位移的对空间的一次微分,所以在三节点三角形元素内部,应变是保持不变的,因此该元素又称为等应变元素(constant strain triangle)。应变既然不变,应力也不变。同样道理,四节点金字塔3D元素,应力及应变也是保持不变的。而四边形四
热心网友
时间:2023-11-06 14:47
ANSYS两个后处理器:
★通用后处理器POST1:查看整个模型在各个时间点上的结果
★时间历程后处理器POST26:查看整个模型上的某一点结果随时间的变化曲线。后处理可在求解完后直接进入,也可在重新进入ANSYS后读入文件进入后处理。
热心网友
时间:2023-11-06 14:47
APDL
APDL换行与续行-
APDL规定每行72个字符
如果要写表达式A=C1+C2 (C1与C2都为表达式
可以用
B=C1
A=B+C2
将一行拆成两行来做但是如果不是表达式,而是输入一个命令参数过多的话,可以用续行命令RMORE,格式如下:
RMORE, R7, R8, R9, R10, R11, R12
这个命令每次也只能输入6个参数,如果多于6个,可以重复使用RMORE就可以输入13-18,19-24等等。另外,于上面续行相应的是换行,一行命令太短可以使用多个命令共一行
$”,没有双引号。这样就可以将一行变成多行使。:) 换行符是“
ANSYS常见术语/命令对照表
Utility Menu 实用菜单
SAVE_DB 存储数据库 RESUME_DB 恢复数据库
Select Entity 选择实体 Comp/Assembly 组元/集合
Plot/Replot 画图/重新画图 Pan,Zoom,Rotate… 平移,缩放,旋转…
WorkPlane(WP) 工作平面 Coordinate System(CS) 坐标系
Macro 宏 Preference… 优先设置…
Preprocessor 前处理 General Postproc 通用后处理
TimeHist Postproc 时间历程后处理 APDL ANSYS参数化设计语言
Line Fillet 在两条线的过渡生成线 Arbitrary 任意形状
Cylinder 圆柱体 Prism 棱柱体 Cone 圆锥形
Sphere 球形 Polygon 多边形 Stress 应力
Strain 应变 Displacement 位移 DOF 自由度
Von Mises(Stress) 平均应力 Contour 等高线(图)
Deformed/Undeformed shape 变形后/未变形的形状 Results Summary 结果摘要
Radiation Matrix 辐射矩阵 Modeling 建模
Meshing 划分网格 Attribute 属性 LS (Load Step) 载荷步
ansys的常用命令介绍
对ANSYS学习也有一个来月的时间了,可是还是什么都不会!郁闷!整理了一些ANSYS常用的命令;但深知自己的水平,还不敢保证完全正确;给大家一些参考,望指正:
1. A,P1,P2,…,P17,P18(以点定义面)
2. AADD,NA1,NA2,…NA8,NA9(面相加)
3. AATT,MAT,REAL,TYPE,ESYS,SECN(指定面的单元属性)
【注】ESYS为坐标系统号、SECN为截面类型号。
4. *ABBR,Abbr,String(定义一个缩略词)
5. ABBRES,Lab,Fname,Ext(从文件中读取缩略词)
6. ABBSAVE,Lab,Fname,Ext(将当前定义的缩略词写入文件)
7. ABS,IR,IA,--,--,Name,--,--,FACTA(取绝对值)
8. ACCAT,NA1,NA2(连接面)
9. ACEL,ACEX,ACEY,ACEZ(定义结构的线性加速度)
10. ACLEAR,NA1,NA2,NINC(清除面单元网格)
11. ADAMS,NMODES,KSTRESS,KSHELL
12. ADAPT, NSOLN, STARGT, TTARGT, FACMN, FACMX, KYKPS, KYMAC
13. ADD,IR, IA, IB, IC, Name, --,-- , FACTA, FACTB, FACTC(变量加运算)
14. ADELE,NA1,NA2,NINC,KSWP(删除面)
【注】KSWP =0删除面但保留面上关键点、1删除面及面上关键点。
15. ADRAG,NL1,NL2,…,NL6,NLP1,NLP2,…,NLP6(将既有线沿一定路径拖拉成面)
16. AESIZE,ANUM,SIZE(指定面上划分单元大小)
17. AFILLT,NA1,NA1,RAD(两面之间生成倒角面)
18. AFSURF,SAREA,TLINE(在既有面单元上生成重叠的表面单元)
19. *AFUN, Lab(指定参数表达式中角度单位)
20. AGEN, ITIME, NA1, NA2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE(复制面)
21. AGLUE,NA1,NA2,…,NA8,NA9(面间相互粘接)
22. AINA,NA1,NA2,…,NA8,NA9(被选面的交集)
23. AINP,NA1,NA2,…,NA8,NA9(面集两两相交)
24. AINV,NA,NV(面体相交)
25. AL,L1,L2,…,L9,L10(以线定义面)
26. ALIST,NA1,NA2,NINC,Lab(列表显示面的信息)
【注】Lab=HPT时,显示面上硬点信息,默认为空。
27. ALLSEL,LabT,Entity(选择所有实体)
【注】LabT=ALL(指定实体及其所有下层实体)、BELOW(指定实体及其下一层实体);
Entity=ALL、VOLU、AREA、LINE、KP、ELEM、NODE。
28. AMESH,NA1,NA2,NINC(划分面生成面单元)
AMESH,AREA,KP1,KP2,KP3,KP4(通过点划分面单元)
29. /AN3D,Kywrd,KEY(三维注释)
30. ANCNTR,NFRAM,DELAY,NCYCL(在POST1中生成结构变形梯度线的动画)
31. ANCUT,NFRAM,DELAY,NCYCL,QOFF,KTOP,TOPOFF,NODE1,NODE2,NODE3(在POST1中生成等势切面云图动画)
32. ANDATA,DELAY,NCYCL,RSLTDAT,MIN,MAX,INCR,FRCLST,AUTOCNTRKY(生成某一范围内的结果数据的顺序梯度线动画)
33. ANDSCL,NFRAM,DELAY,NCYCL(在POST1中生成结构变形的动画)
34. ANFLOW,NFRAM,DELAY,NCYCL,TIME,SPACING,SIZE,LENGTH(生成粒子流或带电粒子运动的动画)
35. /ANGLE,WN,THETA,Axis,KINCR(绕指定轴旋转视图)
36. ANHARM,NFRAM,DELAY,NCYCL(生成谐波分析的动画)
37. ANIM,NCYCL,KCYCL,DELAY(动画显示图形序列)
38. ANISOS,NFRAM,DELAY,NCYCL(在POST1中生成等势面云图动画)
39. ANMODE,NFRAM,DELAY,NCYCL,KACCEL(在POST1中生成结构变形模态的动画)
40. /ANNOT,Lab,VAL1,VAL2(激活图形显示注释)
【注】Lab=OFF、ON、DELE、SAVE、SCALE、XORIG、YORIG、SNAP、STAT、DEFA、REFR、TMODE。
41. ANORM,ANUM,NOEFLIP(重新定义面的法线方向)
【注】NOEFLIP=0、1。
42. ANTIME,NFRAM,DELAY,NCYCL,AUTOCNTRKY,RSLTDAT,MIN,MAX(在指定时间段内生成动画)
详细见网站:www.101315.com
有限元基本概念和原理
有限元分析(FEA,Finite Element Analysis)的基本概念是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成,对每一单元假定一个合适的 (较简单的)近似解,然后推导求解这个域总的满足条件(如结构的平衡条件),从而得到问题的解。这个解不是准确解,而是近似解,因为实际问题被较简单的问 题所代替。由于大多数实际问题难以得到准确解,而有限元不仅计算精度高,而且能适应各种复杂形状,因而成为行之有效的工程分析手段。
有限元是那些集合在一起能够表示实际连续域的离散单元。有限元的概念早在几个世纪前就已产生并得到了应用,例如用多边形(有限个直线单元)*近圆来求得圆 的周长,但作为一种方法而被提出,则是最近的事。有限元法最初被称为矩阵近似方法,应用于航空器的结构强度计算,并由于其方便性、实用性和有效性而引起从 事力学研究的科学家的浓厚兴趣。经过短短数十年的努力,随着计算机技术的快速发展和普及,有限元方法迅速从结构工程强度分析计算扩展到几乎所有的科学技术 领域,成为一种丰富多彩、应用广泛并且实用高效的数值分析方法。
有限元方法与其他求解边值问题近似方法的根本区别在于它的近似性仅限于相对小的子域中。20世纪60年代初首次提出结构力学计算有限元概念的克拉夫(Clough)教授形象地将其描绘为:“有限元法=Rayleigh Ritz法+分片函数”,即有限元法是Rayleigh Ritz法的一种局部化情况。不同于求解(往往是困难的)满足整个定义域边界条件的允许函数的Rayleigh Ritz法,有限元法将函数定义在简单几何形状(如二维问题中的三角形或任意四边形)的单元域上(分片函数),且不考虑整个定义域的复杂边界条件,这是有限元法优于其他近似方法的原因之一。
对于不同物理性质和数学模型的问题,有限元求解法的基本步骤是相同的,只是具体公式推导和运算求解不同。有限元求解问题的基本步骤通常为:
第一步:问题及求解域定义:根据实际问题近似确定求解域的物理性质和几何区域。
第二步:求解域离散化:将求解域近似为具有不同有限大小和形状且彼此相连的有限个单元组成的离散域,习惯上称为有限元网络划分。显然单元越小(网络越细)则离散域的近似程度越好,计算结果也越精确,但计算量及误差都将增大,因此求解域的离散化是有限元法的核心技术之一。
第三步:确定状态变量及控制方法:一个具体的物理问题通常可以用一组包含问题状态变量边界条件的微分方程式表示,为适合有限元求解,通常将微分方程化为等价的泛函形式。
第四步:单元推导:对单元构造一个适合的近似解,即推导有限单元的列式,其中包括选择合理的单元坐标系,建立单元试函数,以某种方法给出单元各状态变量的离散关系,从而形成单元矩阵(结构力学中称刚度阵或柔度阵)。
为保证问题求解的收敛性,单元推导有许多原则要遵循。 对工程应用而言,重要的是应注意每一种单元的解题性能与约束。例如,单元形状应以规则为好,畸形时不仅精度低,而且有缺秩的危险,将导致无法求解。
第五步:总装求解:将单元总装形成离散域的总矩阵方程(联合方程组),反映对近似求解域的离散域的要求,即单元函数的连续性要满足一定的连续条件。总装是在相邻单元结点进行,状态变量及其导数(可能的话)连续性建立在结点处。
第六步:联立方程组求解和结果解释:有限元法最终导致联立方程组。联立方程组的求解可用直接法、选代法和随机法。求解结果是单元结点处状态变量的近似值。对于计算结果的质量,将通过与设计准则提供的允许值比较来评价并确定是否需要重复计算。
简言之,有限元分析可分成三个阶段,前处理、处理和后处理。前处理是建立有限元模型,完成单元网格划分;后处理则是采集处理分析结果,使用户能简便提取信息,了解计算结果。
对称边界条件
有些结构由於具有某些对称性,我们可以在对称面上施加适当的对称边界条件,这样只要建立部分的模型,既省时又省力。
这里讲的对称,不只是几何形状的对称,还包括边界条件、外力施加、材料性质的对称。如果仅有形状对称,其他条件只要有一样不对称,就不能用对称模型求解。
举一个简单的例子。图一是一个含圆孔的平板,边界条件如图所示。我们知道图一的模型上下、左右对称,因此可以简化成图二。
为什麽图二的边界条件可以设为如此?为回答这个问题,我们先将图一分为上下两部分。假设A点与B点是上下两部分相对称的点(图三),则A点与B点的「y」 方向位移方向相反,大小相等。想像A与B点越来越靠近对称面,直到几乎重合在一起(图四),此时A、B两点仍可视为对称於对称面,两点的y方向位移仍为大 小相等,方向相反。然而此时两点几乎重合,因此A、B两点的位移应该相同。为了满足上面的条件,唯一的可能就是在对称面上y方向位移等於零。图二中的另一 边x方向位移等於零,可用相同方法解释。
事实上,图二的两个对称边界,还有一个条件,那就是剪应力xy为零。这个道理的解释方法与位移类似。参考图三,A、B两点由於对称,剪应力也要一模一样; 当A、B两点越来越靠近对称面(图四),两点的剪应力就必须满足牛顿第三运动定律(作用力等於反作用力,方向相反),此时唯一的可能就是剪应力为零。因此 在对称面上,剪应力为零
总之,在对称面上,垂直於对称面的位移以及作用於对称面上的剪应力皆为
在破坏力学的应用
破坏力学是固体力学的一个分支,这门学科主要在探讨四个部分:
1. 含裂纹结构的受到外力时的应力分布;
2. 含裂纹结构受到多大的外力,裂纹会成长;
3. 结构中,裂纹一旦成长,会往那个方向成长;
4. 各种工程结构抵抗裂纹成长的能力,这部分通常由实验决定。
首先先介绍破坏力学的基本概念。图一为含裂纹结构的破坏模式,共分为第一型(mode I or opening mode)、第二型(mode II or sliding mode)以及第三型(mode III or tearing mode)。
裂纹尖端为奇异点(应力正比於根号 r分之一(stress ~ 1/r^0.5),r为结构中某一点与裂纹尖端的距离),去探讨裂纹尖端(或附近)的应力有多大是没有意义的。有鉴於此,我们引入了一个参数:应力强度因子(stress intensity factor)。应力强度因子共有三个,一般写做KI、KII以及KIII,分别对应到三种不同的破坏模式。
在线弹性破坏力学(linear elastic fracture mechanics,简称LEFM)中,外力的大小与应力强度因子成正比关系。此外,应力强度因子也和几何参数(例如裂纹长度、外力与裂纹的距离等)有 关。对於同一结构,应力强度因子可视为含裂纹结构所受外力大小的一个指标;换言之,应力强度因子越大,结构越危险。当应力强度因子超过「破坏韧性 (fracture toughness),记为Kc」,裂纹开始成长。破坏韧性通常由制作标准试片由实验求得,不同的材料有不同的破坏韧性,所以破坏韧性可视为材料性质。
至於裂缝会往那个方向成长?有好几个理论可以预测,例如能量释放率理论(energy release rate)、最大周向应力理论(maximum circumferential stress)、J积分理论(J-integral)、应变能密度理论(strain energy density theory)。以上几个理论,除了最大周向应力之外,基本上都和应变能有关,这个部分较为深入,不详述了。
接着我们要探讨有限元素法在破坏力学的应用。前面讲过,探讨裂缝尖端附近的应力有多大是没有意义的。因此我们的重点在於,如何运用有限元素法求得裂纹尖端的应力强度因子。最常见的就是利用「四分之一节点(quarter point)」元素来模拟裂纹尖端,如图二。严格来说,四分之一节点元素不是「三角形元素」,它是由四边形元素退化而成,退化方式如图三。一般我们利用节点上的位移,求得裂纹尖端的开口位移(crack tip opening displacement)来反推得到应力强度因子。MARC及ANSYS皆用四分之一节点元素来算应力强度因子。
以上所说的四分之一节点元素,其最原始的统御方程式和弹性力学并没有两样,只是将其中几个节点变了位置,因此形状函数也跟着改变,最後得到一个结果:这种元素内部的应力呈现根号r分之一的奇异性,因此可用来模拟裂纹。
另外有些软体,例如NASTRAN,则利用另一种特殊元素,共有18个节点,如图四。这种元素最原始的统御方程式,和一般的弹性力学稍有不同,它一开始就假设应力正比於根号 r分之一。
不管如何,裂纹尖端的元素与一般的元素不同,裂纹尖端网格的大小,由算出来的应力强度因子准不准来决定。如何知道准不准?可建立一个简单模型,有理论解的 题目,即可比较,例如图四。以平面问题为例,若采用四分之一节点元素,元素边长大约在裂纹长度的百分之六左右;若采用NASTRAN的特殊元素,其长度更 只要裂纹长度的百分之十至二十即可。当然以上的准则并非万用,不同的题目最好还是多试几个网格粗细,再做决定。
方向
元素有方向性。有些元素如果方向搞错,则跑不出来;有些虽然跑得出来,但结果却有问题。 以平面元素(包括平面应变、平面应力以及轴对称)为例,四个节点的编号必须是「逆时钟方向」。如果是顺时钟方向,则在有限元素的定义中,这个元素的面积小於零,就跑不出来了。
板(或壳)元素也有方向性。板元素的方向由右手定则决定,也就是元素的正向方向,由元素编号顺序的方向决定,如此也定义了板元素的上表面或下表面。当板受 到bending时,上表面和下表面分别会受到张应力以及压应力(当然也可能颠倒,视bending的方向而定)。图一为两个相邻的板元素受到 bending作用,元素1的方向和元素2的方向不同,因此它们的上、下表面也不一致(参考图一)。在受到图一所示的bending时,元素1的上表面为 张应力,下表面为压应力;元素2的上表面则为压应力,下表面为张应力。有些有限元素软体(如MARC)在计算板元素的「上表面」节点应力时,是根据两个元 素的「上表面」高斯点上应力经外差到节点上再平均求得,但注意图一中两个元素的上下表面定义颠倒,所以平均後的上表面节点应力就很不准了。
再解释清楚一点,我们先将两个元素分开,同时画上厚度,比较清楚,如图二,图一中的A点在图二则可分为A1至A4四点。如果单独看元素1,A1点应该受到 张应力,A2为压应力;单独看元素2,A3为张应力,A4为压应力。当有限元素法在计算A点「上表面」应力时,是将元素1以及元素2的「上表面」应力平均 (即将A1及A4的应力平均),本来应该上表面是张应力的,被这麽一搞,应力下降,误差就来了。
一般有限元素法软体,在做automesh时,元素的方向与surface(或area)的方向一致,因此同一个surface(或area)做出来的网 格,方向都相同。但如果两个surface(或area)相邻,就要特别注意方向有没有一致。如果是「手动」建立网格,更要特别注意建立出来的元素的方 向,以免造成困扰。
形状函数
形状函数在有限元素法中是非常重要的一个概念,它定义了元素内部位移的分布。以一维线性元素为例(图一) 此种元素有两个节点,分别以节点1及2表示。节点1的座标为xi = -1,节点2则为xi = -2。元素中有几个节点,就有几个形状函数。因此我们有两个形状函数。形状函数有个特点:考虑第n个形状函数,若代入第n个节点的座标,其函数值为1;若 代入其他节点的座标,则函数值为0。由於图一为两节点元素,所以形状函数为线性,其表示式为:其他种类的元素的形状函数,一般有限元素法书籍(或商用有限 元素法软体的使用手册)皆有提到,这里不多写了。
元素的特性主要由形状函数所主宰。例如三节点之三角形元素,沿着元素内部任一方向,位移皆为线性分布。为什麽我们知道它是线性分布?只要查查该元素的形状 函数立刻就明白了。由於应变为位移的对空间的一次微分,所以在三节点三角形元素内部,应变是保持不变的,因此该元素又称为等应变元素(constant strain triangle)。应变既然不变,应力也不变。同样道理,四节点金字塔3D元素,应力及应变也是保持不变的。而四边形四
热心网友
时间:2023-11-06 14:47
ANSYS两个后处理器:
★通用后处理器POST1:查看整个模型在各个时间点上的结果
★时间历程后处理器POST26:查看整个模型上的某一点结果随时间的变化曲线。后处理可在求解完后直接进入,也可在重新进入ANSYS后读入文件进入后处理。
热心网友
时间:2023-10-15 13:49
APDL
APDL换行与续行-
APDL规定每行72个字符
如果要写表达式A=C1+C2 (C1与C2都为表达式
可以用
B=C1
A=B+C2
将一行拆成两行来做但是如果不是表达式,而是输入一个命令参数过多的话,可以用续行命令RMORE,格式如下:
RMORE, R7, R8, R9, R10, R11, R12
这个命令每次也只能输入6个参数,如果多于6个,可以重复使用RMORE就可以输入13-18,19-24等等。另外,于上面续行相应的是换行,一行命令太短可以使用多个命令共一行
$”,没有双引号。这样就可以将一行变成多行使。:) 换行符是“
ANSYS常见术语/命令对照表
Utility Menu 实用菜单
SAVE_DB 存储数据库 RESUME_DB 恢复数据库
Select Entity 选择实体 Comp/Assembly 组元/集合
Plot/Replot 画图/重新画图 Pan,Zoom,Rotate… 平移,缩放,旋转…
WorkPlane(WP) 工作平面 Coordinate System(CS) 坐标系
Macro 宏 Preference… 优先设置…
Preprocessor 前处理 General Postproc 通用后处理
TimeHist Postproc 时间历程后处理 APDL ANSYS参数化设计语言
Line Fillet 在两条线的过渡生成线 Arbitrary 任意形状
Cylinder 圆柱体 Prism 棱柱体 Cone 圆锥形
Sphere 球形 Polygon 多边形 Stress 应力
Strain 应变 Displacement 位移 DOF 自由度
Von Mises(Stress) 平均应力 Contour 等高线(图)
Deformed/Undeformed shape 变形后/未变形的形状 Results Summary 结果摘要
Radiation Matrix 辐射矩阵 Modeling 建模
Meshing 划分网格 Attribute 属性 LS (Load Step) 载荷步
ansys的常用命令介绍
对ANSYS学习也有一个来月的时间了,可是还是什么都不会!郁闷!整理了一些ANSYS常用的命令;但深知自己的水平,还不敢保证完全正确;给大家一些参考,望指正:
1. A,P1,P2,…,P17,P18(以点定义面)
2. AADD,NA1,NA2,…NA8,NA9(面相加)
3. AATT,MAT,REAL,TYPE,ESYS,SECN(指定面的单元属性)
【注】ESYS为坐标系统号、SECN为截面类型号。
4. *ABBR,Abbr,String(定义一个缩略词)
5. ABBRES,Lab,Fname,Ext(从文件中读取缩略词)
6. ABBSAVE,Lab,Fname,Ext(将当前定义的缩略词写入文件)
7. ABS,IR,IA,--,--,Name,--,--,FACTA(取绝对值)
8. ACCAT,NA1,NA2(连接面)
9. ACEL,ACEX,ACEY,ACEZ(定义结构的线性加速度)
10. ACLEAR,NA1,NA2,NINC(清除面单元网格)
11. ADAMS,NMODES,KSTRESS,KSHELL
12. ADAPT, NSOLN, STARGT, TTARGT, FACMN, FACMX, KYKPS, KYMAC
13. ADD,IR, IA, IB, IC, Name, --,-- , FACTA, FACTB, FACTC(变量加运算)
14. ADELE,NA1,NA2,NINC,KSWP(删除面)
【注】KSWP =0删除面但保留面上关键点、1删除面及面上关键点。
15. ADRAG,NL1,NL2,…,NL6,NLP1,NLP2,…,NLP6(将既有线沿一定路径拖拉成面)
16. AESIZE,ANUM,SIZE(指定面上划分单元大小)
17. AFILLT,NA1,NA1,RAD(两面之间生成倒角面)
18. AFSURF,SAREA,TLINE(在既有面单元上生成重叠的表面单元)
19. *AFUN, Lab(指定参数表达式中角度单位)
20. AGEN, ITIME, NA1, NA2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE(复制面)
21. AGLUE,NA1,NA2,…,NA8,NA9(面间相互粘接)
22. AINA,NA1,NA2,…,NA8,NA9(被选面的交集)
23. AINP,NA1,NA2,…,NA8,NA9(面集两两相交)
24. AINV,NA,NV(面体相交)
25. AL,L1,L2,…,L9,L10(以线定义面)
26. ALIST,NA1,NA2,NINC,Lab(列表显示面的信息)
【注】Lab=HPT时,显示面上硬点信息,默认为空。
27. ALLSEL,LabT,Entity(选择所有实体)
【注】LabT=ALL(指定实体及其所有下层实体)、BELOW(指定实体及其下一层实体);
Entity=ALL、VOLU、AREA、LINE、KP、ELEM、NODE。
28. AMESH,NA1,NA2,NINC(划分面生成面单元)
AMESH,AREA,KP1,KP2,KP3,KP4(通过点划分面单元)
29. /AN3D,Kywrd,KEY(三维注释)
30. ANCNTR,NFRAM,DELAY,NCYCL(在POST1中生成结构变形梯度线的动画)
31. ANCUT,NFRAM,DELAY,NCYCL,QOFF,KTOP,TOPOFF,NODE1,NODE2,NODE3(在POST1中生成等势切面云图动画)
32. ANDATA,DELAY,NCYCL,RSLTDAT,MIN,MAX,INCR,FRCLST,AUTOCNTRKY(生成某一范围内的结果数据的顺序梯度线动画)
33. ANDSCL,NFRAM,DELAY,NCYCL(在POST1中生成结构变形的动画)
34. ANFLOW,NFRAM,DELAY,NCYCL,TIME,SPACING,SIZE,LENGTH(生成粒子流或带电粒子运动的动画)
35. /ANGLE,WN,THETA,Axis,KINCR(绕指定轴旋转视图)
36. ANHARM,NFRAM,DELAY,NCYCL(生成谐波分析的动画)
37. ANIM,NCYCL,KCYCL,DELAY(动画显示图形序列)
38. ANISOS,NFRAM,DELAY,NCYCL(在POST1中生成等势面云图动画)
39. ANMODE,NFRAM,DELAY,NCYCL,KACCEL(在POST1中生成结构变形模态的动画)
40. /ANNOT,Lab,VAL1,VAL2(激活图形显示注释)
【注】Lab=OFF、ON、DELE、SAVE、SCALE、XORIG、YORIG、SNAP、STAT、DEFA、REFR、TMODE。
41. ANORM,ANUM,NOEFLIP(重新定义面的法线方向)
【注】NOEFLIP=0、1。
42. ANTIME,NFRAM,DELAY,NCYCL,AUTOCNTRKY,RSLTDAT,MIN,MAX(在指定时间段内生成动画)
详细见网站:www.101315.com
有限元基本概念和原理
有限元分析(FEA,Finite Element Analysis)的基本概念是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成,对每一单元假定一个合适的 (较简单的)近似解,然后推导求解这个域总的满足条件(如结构的平衡条件),从而得到问题的解。这个解不是准确解,而是近似解,因为实际问题被较简单的问 题所代替。由于大多数实际问题难以得到准确解,而有限元不仅计算精度高,而且能适应各种复杂形状,因而成为行之有效的工程分析手段。
有限元是那些集合在一起能够表示实际连续域的离散单元。有限元的概念早在几个世纪前就已产生并得到了应用,例如用多边形(有限个直线单元)*近圆来求得圆 的周长,但作为一种方法而被提出,则是最近的事。有限元法最初被称为矩阵近似方法,应用于航空器的结构强度计算,并由于其方便性、实用性和有效性而引起从 事力学研究的科学家的浓厚兴趣。经过短短数十年的努力,随着计算机技术的快速发展和普及,有限元方法迅速从结构工程强度分析计算扩展到几乎所有的科学技术 领域,成为一种丰富多彩、应用广泛并且实用高效的数值分析方法。
有限元方法与其他求解边值问题近似方法的根本区别在于它的近似性仅限于相对小的子域中。20世纪60年代初首次提出结构力学计算有限元概念的克拉夫(Clough)教授形象地将其描绘为:“有限元法=Rayleigh Ritz法+分片函数”,即有限元法是Rayleigh Ritz法的一种局部化情况。不同于求解(往往是困难的)满足整个定义域边界条件的允许函数的Rayleigh Ritz法,有限元法将函数定义在简单几何形状(如二维问题中的三角形或任意四边形)的单元域上(分片函数),且不考虑整个定义域的复杂边界条件,这是有限元法优于其他近似方法的原因之一。
对于不同物理性质和数学模型的问题,有限元求解法的基本步骤是相同的,只是具体公式推导和运算求解不同。有限元求解问题的基本步骤通常为:
第一步:问题及求解域定义:根据实际问题近似确定求解域的物理性质和几何区域。
第二步:求解域离散化:将求解域近似为具有不同有限大小和形状且彼此相连的有限个单元组成的离散域,习惯上称为有限元网络划分。显然单元越小(网络越细)则离散域的近似程度越好,计算结果也越精确,但计算量及误差都将增大,因此求解域的离散化是有限元法的核心技术之一。
第三步:确定状态变量及控制方法:一个具体的物理问题通常可以用一组包含问题状态变量边界条件的微分方程式表示,为适合有限元求解,通常将微分方程化为等价的泛函形式。
第四步:单元推导:对单元构造一个适合的近似解,即推导有限单元的列式,其中包括选择合理的单元坐标系,建立单元试函数,以某种方法给出单元各状态变量的离散关系,从而形成单元矩阵(结构力学中称刚度阵或柔度阵)。
为保证问题求解的收敛性,单元推导有许多原则要遵循。 对工程应用而言,重要的是应注意每一种单元的解题性能与约束。例如,单元形状应以规则为好,畸形时不仅精度低,而且有缺秩的危险,将导致无法求解。
第五步:总装求解:将单元总装形成离散域的总矩阵方程(联合方程组),反映对近似求解域的离散域的要求,即单元函数的连续性要满足一定的连续条件。总装是在相邻单元结点进行,状态变量及其导数(可能的话)连续性建立在结点处。
第六步:联立方程组求解和结果解释:有限元法最终导致联立方程组。联立方程组的求解可用直接法、选代法和随机法。求解结果是单元结点处状态变量的近似值。对于计算结果的质量,将通过与设计准则提供的允许值比较来评价并确定是否需要重复计算。
简言之,有限元分析可分成三个阶段,前处理、处理和后处理。前处理是建立有限元模型,完成单元网格划分;后处理则是采集处理分析结果,使用户能简便提取信息,了解计算结果。
对称边界条件
有些结构由於具有某些对称性,我们可以在对称面上施加适当的对称边界条件,这样只要建立部分的模型,既省时又省力。
这里讲的对称,不只是几何形状的对称,还包括边界条件、外力施加、材料性质的对称。如果仅有形状对称,其他条件只要有一样不对称,就不能用对称模型求解。
举一个简单的例子。图一是一个含圆孔的平板,边界条件如图所示。我们知道图一的模型上下、左右对称,因此可以简化成图二。
为什麽图二的边界条件可以设为如此?为回答这个问题,我们先将图一分为上下两部分。假设A点与B点是上下两部分相对称的点(图三),则A点与B点的「y」 方向位移方向相反,大小相等。想像A与B点越来越靠近对称面,直到几乎重合在一起(图四),此时A、B两点仍可视为对称於对称面,两点的y方向位移仍为大 小相等,方向相反。然而此时两点几乎重合,因此A、B两点的位移应该相同。为了满足上面的条件,唯一的可能就是在对称面上y方向位移等於零。图二中的另一 边x方向位移等於零,可用相同方法解释。
事实上,图二的两个对称边界,还有一个条件,那就是剪应力xy为零。这个道理的解释方法与位移类似。参考图三,A、B两点由於对称,剪应力也要一模一样; 当A、B两点越来越靠近对称面(图四),两点的剪应力就必须满足牛顿第三运动定律(作用力等於反作用力,方向相反),此时唯一的可能就是剪应力为零。因此 在对称面上,剪应力为零
总之,在对称面上,垂直於对称面的位移以及作用於对称面上的剪应力皆为
在破坏力学的应用
破坏力学是固体力学的一个分支,这门学科主要在探讨四个部分:
1. 含裂纹结构的受到外力时的应力分布;
2. 含裂纹结构受到多大的外力,裂纹会成长;
3. 结构中,裂纹一旦成长,会往那个方向成长;
4. 各种工程结构抵抗裂纹成长的能力,这部分通常由实验决定。
首先先介绍破坏力学的基本概念。图一为含裂纹结构的破坏模式,共分为第一型(mode I or opening mode)、第二型(mode II or sliding mode)以及第三型(mode III or tearing mode)。
裂纹尖端为奇异点(应力正比於根号 r分之一(stress ~ 1/r^0.5),r为结构中某一点与裂纹尖端的距离),去探讨裂纹尖端(或附近)的应力有多大是没有意义的。有鉴於此,我们引入了一个参数:应力强度因子(stress intensity factor)。应力强度因子共有三个,一般写做KI、KII以及KIII,分别对应到三种不同的破坏模式。
在线弹性破坏力学(linear elastic fracture mechanics,简称LEFM)中,外力的大小与应力强度因子成正比关系。此外,应力强度因子也和几何参数(例如裂纹长度、外力与裂纹的距离等)有 关。对於同一结构,应力强度因子可视为含裂纹结构所受外力大小的一个指标;换言之,应力强度因子越大,结构越危险。当应力强度因子超过「破坏韧性 (fracture toughness),记为Kc」,裂纹开始成长。破坏韧性通常由制作标准试片由实验求得,不同的材料有不同的破坏韧性,所以破坏韧性可视为材料性质。
至於裂缝会往那个方向成长?有好几个理论可以预测,例如能量释放率理论(energy release rate)、最大周向应力理论(maximum circumferential stress)、J积分理论(J-integral)、应变能密度理论(strain energy density theory)。以上几个理论,除了最大周向应力之外,基本上都和应变能有关,这个部分较为深入,不详述了。
接着我们要探讨有限元素法在破坏力学的应用。前面讲过,探讨裂缝尖端附近的应力有多大是没有意义的。因此我们的重点在於,如何运用有限元素法求得裂纹尖端的应力强度因子。最常见的就是利用「四分之一节点(quarter point)」元素来模拟裂纹尖端,如图二。严格来说,四分之一节点元素不是「三角形元素」,它是由四边形元素退化而成,退化方式如图三。一般我们利用节点上的位移,求得裂纹尖端的开口位移(crack tip opening displacement)来反推得到应力强度因子。MARC及ANSYS皆用四分之一节点元素来算应力强度因子。
以上所说的四分之一节点元素,其最原始的统御方程式和弹性力学并没有两样,只是将其中几个节点变了位置,因此形状函数也跟着改变,最後得到一个结果:这种元素内部的应力呈现根号r分之一的奇异性,因此可用来模拟裂纹。
另外有些软体,例如NASTRAN,则利用另一种特殊元素,共有18个节点,如图四。这种元素最原始的统御方程式,和一般的弹性力学稍有不同,它一开始就假设应力正比於根号 r分之一。
不管如何,裂纹尖端的元素与一般的元素不同,裂纹尖端网格的大小,由算出来的应力强度因子准不准来决定。如何知道准不准?可建立一个简单模型,有理论解的 题目,即可比较,例如图四。以平面问题为例,若采用四分之一节点元素,元素边长大约在裂纹长度的百分之六左右;若采用NASTRAN的特殊元素,其长度更 只要裂纹长度的百分之十至二十即可。当然以上的准则并非万用,不同的题目最好还是多试几个网格粗细,再做决定。
方向
元素有方向性。有些元素如果方向搞错,则跑不出来;有些虽然跑得出来,但结果却有问题。 以平面元素(包括平面应变、平面应力以及轴对称)为例,四个节点的编号必须是「逆时钟方向」。如果是顺时钟方向,则在有限元素的定义中,这个元素的面积小於零,就跑不出来了。
板(或壳)元素也有方向性。板元素的方向由右手定则决定,也就是元素的正向方向,由元素编号顺序的方向决定,如此也定义了板元素的上表面或下表面。当板受 到bending时,上表面和下表面分别会受到张应力以及压应力(当然也可能颠倒,视bending的方向而定)。图一为两个相邻的板元素受到 bending作用,元素1的方向和元素2的方向不同,因此它们的上、下表面也不一致(参考图一)。在受到图一所示的bending时,元素1的上表面为 张应力,下表面为压应力;元素2的上表面则为压应力,下表面为张应力。有些有限元素软体(如MARC)在计算板元素的「上表面」节点应力时,是根据两个元 素的「上表面」高斯点上应力经外差到节点上再平均求得,但注意图一中两个元素的上下表面定义颠倒,所以平均後的上表面节点应力就很不准了。
再解释清楚一点,我们先将两个元素分开,同时画上厚度,比较清楚,如图二,图一中的A点在图二则可分为A1至A4四点。如果单独看元素1,A1点应该受到 张应力,A2为压应力;单独看元素2,A3为张应力,A4为压应力。当有限元素法在计算A点「上表面」应力时,是将元素1以及元素2的「上表面」应力平均 (即将A1及A4的应力平均),本来应该上表面是张应力的,被这麽一搞,应力下降,误差就来了。
一般有限元素法软体,在做automesh时,元素的方向与surface(或area)的方向一致,因此同一个surface(或area)做出来的网 格,方向都相同。但如果两个surface(或area)相邻,就要特别注意方向有没有一致。如果是「手动」建立网格,更要特别注意建立出来的元素的方 向,以免造成困扰。
形状函数
形状函数在有限元素法中是非常重要的一个概念,它定义了元素内部位移的分布。以一维线性元素为例(图一) 此种元素有两个节点,分别以节点1及2表示。节点1的座标为xi = -1,节点2则为xi = -2。元素中有几个节点,就有几个形状函数。因此我们有两个形状函数。形状函数有个特点:考虑第n个形状函数,若代入第n个节点的座标,其函数值为1;若 代入其他节点的座标,则函数值为0。由於图一为两节点元素,所以形状函数为线性,其表示式为:其他种类的元素的形状函数,一般有限元素法书籍(或商用有限 元素法软体的使用手册)皆有提到,这里不多写了。
元素的特性主要由形状函数所主宰。例如三节点之三角形元素,沿着元素内部任一方向,位移皆为线性分布。为什麽我们知道它是线性分布?只要查查该元素的形状 函数立刻就明白了。由於应变为位移的对空间的一次微分,所以在三节点三角形元素内部,应变是保持不变的,因此该元素又称为等应变元素(constant strain triangle)。应变既然不变,应力也不变。同样道理,四节点金字塔3D元素,应力及应变也是保持不变的。而四边形四
热心网友
时间:2023-10-15 13:49
ANSYS两个后处理器:
★通用后处理器POST1:查看整个模型在各个时间点上的结果
★时间历程后处理器POST26:查看整个模型上的某一点结果随时间的变化曲线。后处理可在求解完后直接进入,也可在重新进入ANSYS后读入文件进入后处理。
热心网友
时间:2023-10-15 13:49
APDL
APDL换行与续行-
APDL规定每行72个字符
如果要写表达式A=C1+C2 (C1与C2都为表达式
可以用
B=C1
A=B+C2
将一行拆成两行来做但是如果不是表达式,而是输入一个命令参数过多的话,可以用续行命令RMORE,格式如下:
RMORE, R7, R8, R9, R10, R11, R12
这个命令每次也只能输入6个参数,如果多于6个,可以重复使用RMORE就可以输入13-18,19-24等等。另外,于上面续行相应的是换行,一行命令太短可以使用多个命令共一行
$”,没有双引号。这样就可以将一行变成多行使。:) 换行符是“
ANSYS常见术语/命令对照表
Utility Menu 实用菜单
SAVE_DB 存储数据库 RESUME_DB 恢复数据库
Select Entity 选择实体 Comp/Assembly 组元/集合
Plot/Replot 画图/重新画图 Pan,Zoom,Rotate… 平移,缩放,旋转…
WorkPlane(WP) 工作平面 Coordinate System(CS) 坐标系
Macro 宏 Preference… 优先设置…
Preprocessor 前处理 General Postproc 通用后处理
TimeHist Postproc 时间历程后处理 APDL ANSYS参数化设计语言
Line Fillet 在两条线的过渡生成线 Arbitrary 任意形状
Cylinder 圆柱体 Prism 棱柱体 Cone 圆锥形
Sphere 球形 Polygon 多边形 Stress 应力
Strain 应变 Displacement 位移 DOF 自由度
Von Mises(Stress) 平均应力 Contour 等高线(图)
Deformed/Undeformed shape 变形后/未变形的形状 Results Summary 结果摘要
Radiation Matrix 辐射矩阵 Modeling 建模
Meshing 划分网格 Attribute 属性 LS (Load Step) 载荷步
ansys的常用命令介绍
对ANSYS学习也有一个来月的时间了,可是还是什么都不会!郁闷!整理了一些ANSYS常用的命令;但深知自己的水平,还不敢保证完全正确;给大家一些参考,望指正:
1. A,P1,P2,…,P17,P18(以点定义面)
2. AADD,NA1,NA2,…NA8,NA9(面相加)
3. AATT,MAT,REAL,TYPE,ESYS,SECN(指定面的单元属性)
【注】ESYS为坐标系统号、SECN为截面类型号。
4. *ABBR,Abbr,String(定义一个缩略词)
5. ABBRES,Lab,Fname,Ext(从文件中读取缩略词)
6. ABBSAVE,Lab,Fname,Ext(将当前定义的缩略词写入文件)
7. ABS,IR,IA,--,--,Name,--,--,FACTA(取绝对值)
8. ACCAT,NA1,NA2(连接面)
9. ACEL,ACEX,ACEY,ACEZ(定义结构的线性加速度)
10. ACLEAR,NA1,NA2,NINC(清除面单元网格)
11. ADAMS,NMODES,KSTRESS,KSHELL
12. ADAPT, NSOLN, STARGT, TTARGT, FACMN, FACMX, KYKPS, KYMAC
13. ADD,IR, IA, IB, IC, Name, --,-- , FACTA, FACTB, FACTC(变量加运算)
14. ADELE,NA1,NA2,NINC,KSWP(删除面)
【注】KSWP =0删除面但保留面上关键点、1删除面及面上关键点。
15. ADRAG,NL1,NL2,…,NL6,NLP1,NLP2,…,NLP6(将既有线沿一定路径拖拉成面)
16. AESIZE,ANUM,SIZE(指定面上划分单元大小)
17. AFILLT,NA1,NA1,RAD(两面之间生成倒角面)
18. AFSURF,SAREA,TLINE(在既有面单元上生成重叠的表面单元)
19. *AFUN, Lab(指定参数表达式中角度单位)
20. AGEN, ITIME, NA1, NA2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE(复制面)
21. AGLUE,NA1,NA2,…,NA8,NA9(面间相互粘接)
22. AINA,NA1,NA2,…,NA8,NA9(被选面的交集)
23. AINP,NA1,NA2,…,NA8,NA9(面集两两相交)
24. AINV,NA,NV(面体相交)
25. AL,L1,L2,…,L9,L10(以线定义面)
26. ALIST,NA1,NA2,NINC,Lab(列表显示面的信息)
【注】Lab=HPT时,显示面上硬点信息,默认为空。
27. ALLSEL,LabT,Entity(选择所有实体)
【注】LabT=ALL(指定实体及其所有下层实体)、BELOW(指定实体及其下一层实体);
Entity=ALL、VOLU、AREA、LINE、KP、ELEM、NODE。
28. AMESH,NA1,NA2,NINC(划分面生成面单元)
AMESH,AREA,KP1,KP2,KP3,KP4(通过点划分面单元)
29. /AN3D,Kywrd,KEY(三维注释)
30. ANCNTR,NFRAM,DELAY,NCYCL(在POST1中生成结构变形梯度线的动画)
31. ANCUT,NFRAM,DELAY,NCYCL,QOFF,KTOP,TOPOFF,NODE1,NODE2,NODE3(在POST1中生成等势切面云图动画)
32. ANDATA,DELAY,NCYCL,RSLTDAT,MIN,MAX,INCR,FRCLST,AUTOCNTRKY(生成某一范围内的结果数据的顺序梯度线动画)
33. ANDSCL,NFRAM,DELAY,NCYCL(在POST1中生成结构变形的动画)
34. ANFLOW,NFRAM,DELAY,NCYCL,TIME,SPACING,SIZE,LENGTH(生成粒子流或带电粒子运动的动画)
35. /ANGLE,WN,THETA,Axis,KINCR(绕指定轴旋转视图)
36. ANHARM,NFRAM,DELAY,NCYCL(生成谐波分析的动画)
37. ANIM,NCYCL,KCYCL,DELAY(动画显示图形序列)
38. ANISOS,NFRAM,DELAY,NCYCL(在POST1中生成等势面云图动画)
39. ANMODE,NFRAM,DELAY,NCYCL,KACCEL(在POST1中生成结构变形模态的动画)
40. /ANNOT,Lab,VAL1,VAL2(激活图形显示注释)
【注】Lab=OFF、ON、DELE、SAVE、SCALE、XORIG、YORIG、SNAP、STAT、DEFA、REFR、TMODE。
41. ANORM,ANUM,NOEFLIP(重新定义面的法线方向)
【注】NOEFLIP=0、1。
42. ANTIME,NFRAM,DELAY,NCYCL,AUTOCNTRKY,RSLTDAT,MIN,MAX(在指定时间段内生成动画)
详细见网站:www.101315.com
有限元基本概念和原理
有限元分析(FEA,Finite Element Analysis)的基本概念是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成,对每一单元假定一个合适的 (较简单的)近似解,然后推导求解这个域总的满足条件(如结构的平衡条件),从而得到问题的解。这个解不是准确解,而是近似解,因为实际问题被较简单的问 题所代替。由于大多数实际问题难以得到准确解,而有限元不仅计算精度高,而且能适应各种复杂形状,因而成为行之有效的工程分析手段。
有限元是那些集合在一起能够表示实际连续域的离散单元。有限元的概念早在几个世纪前就已产生并得到了应用,例如用多边形(有限个直线单元)*近圆来求得圆 的周长,但作为一种方法而被提出,则是最近的事。有限元法最初被称为矩阵近似方法,应用于航空器的结构强度计算,并由于其方便性、实用性和有效性而引起从 事力学研究的科学家的浓厚兴趣。经过短短数十年的努力,随着计算机技术的快速发展和普及,有限元方法迅速从结构工程强度分析计算扩展到几乎所有的科学技术 领域,成为一种丰富多彩、应用广泛并且实用高效的数值分析方法。
有限元方法与其他求解边值问题近似方法的根本区别在于它的近似性仅限于相对小的子域中。20世纪60年代初首次提出结构力学计算有限元概念的克拉夫(Clough)教授形象地将其描绘为:“有限元法=Rayleigh Ritz法+分片函数”,即有限元法是Rayleigh Ritz法的一种局部化情况。不同于求解(往往是困难的)满足整个定义域边界条件的允许函数的Rayleigh Ritz法,有限元法将函数定义在简单几何形状(如二维问题中的三角形或任意四边形)的单元域上(分片函数),且不考虑整个定义域的复杂边界条件,这是有限元法优于其他近似方法的原因之一。
对于不同物理性质和数学模型的问题,有限元求解法的基本步骤是相同的,只是具体公式推导和运算求解不同。有限元求解问题的基本步骤通常为:
第一步:问题及求解域定义:根据实际问题近似确定求解域的物理性质和几何区域。
第二步:求解域离散化:将求解域近似为具有不同有限大小和形状且彼此相连的有限个单元组成的离散域,习惯上称为有限元网络划分。显然单元越小(网络越细)则离散域的近似程度越好,计算结果也越精确,但计算量及误差都将增大,因此求解域的离散化是有限元法的核心技术之一。
第三步:确定状态变量及控制方法:一个具体的物理问题通常可以用一组包含问题状态变量边界条件的微分方程式表示,为适合有限元求解,通常将微分方程化为等价的泛函形式。
第四步:单元推导:对单元构造一个适合的近似解,即推导有限单元的列式,其中包括选择合理的单元坐标系,建立单元试函数,以某种方法给出单元各状态变量的离散关系,从而形成单元矩阵(结构力学中称刚度阵或柔度阵)。
为保证问题求解的收敛性,单元推导有许多原则要遵循。 对工程应用而言,重要的是应注意每一种单元的解题性能与约束。例如,单元形状应以规则为好,畸形时不仅精度低,而且有缺秩的危险,将导致无法求解。
第五步:总装求解:将单元总装形成离散域的总矩阵方程(联合方程组),反映对近似求解域的离散域的要求,即单元函数的连续性要满足一定的连续条件。总装是在相邻单元结点进行,状态变量及其导数(可能的话)连续性建立在结点处。
第六步:联立方程组求解和结果解释:有限元法最终导致联立方程组。联立方程组的求解可用直接法、选代法和随机法。求解结果是单元结点处状态变量的近似值。对于计算结果的质量,将通过与设计准则提供的允许值比较来评价并确定是否需要重复计算。
简言之,有限元分析可分成三个阶段,前处理、处理和后处理。前处理是建立有限元模型,完成单元网格划分;后处理则是采集处理分析结果,使用户能简便提取信息,了解计算结果。
对称边界条件
有些结构由於具有某些对称性,我们可以在对称面上施加适当的对称边界条件,这样只要建立部分的模型,既省时又省力。
这里讲的对称,不只是几何形状的对称,还包括边界条件、外力施加、材料性质的对称。如果仅有形状对称,其他条件只要有一样不对称,就不能用对称模型求解。
举一个简单的例子。图一是一个含圆孔的平板,边界条件如图所示。我们知道图一的模型上下、左右对称,因此可以简化成图二。
为什麽图二的边界条件可以设为如此?为回答这个问题,我们先将图一分为上下两部分。假设A点与B点是上下两部分相对称的点(图三),则A点与B点的「y」 方向位移方向相反,大小相等。想像A与B点越来越靠近对称面,直到几乎重合在一起(图四),此时A、B两点仍可视为对称於对称面,两点的y方向位移仍为大 小相等,方向相反。然而此时两点几乎重合,因此A、B两点的位移应该相同。为了满足上面的条件,唯一的可能就是在对称面上y方向位移等於零。图二中的另一 边x方向位移等於零,可用相同方法解释。
事实上,图二的两个对称边界,还有一个条件,那就是剪应力xy为零。这个道理的解释方法与位移类似。参考图三,A、B两点由於对称,剪应力也要一模一样; 当A、B两点越来越靠近对称面(图四),两点的剪应力就必须满足牛顿第三运动定律(作用力等於反作用力,方向相反),此时唯一的可能就是剪应力为零。因此 在对称面上,剪应力为零
总之,在对称面上,垂直於对称面的位移以及作用於对称面上的剪应力皆为
在破坏力学的应用
破坏力学是固体力学的一个分支,这门学科主要在探讨四个部分:
1. 含裂纹结构的受到外力时的应力分布;
2. 含裂纹结构受到多大的外力,裂纹会成长;
3. 结构中,裂纹一旦成长,会往那个方向成长;
4. 各种工程结构抵抗裂纹成长的能力,这部分通常由实验决定。
首先先介绍破坏力学的基本概念。图一为含裂纹结构的破坏模式,共分为第一型(mode I or opening mode)、第二型(mode II or sliding mode)以及第三型(mode III or tearing mode)。
裂纹尖端为奇异点(应力正比於根号 r分之一(stress ~ 1/r^0.5),r为结构中某一点与裂纹尖端的距离),去探讨裂纹尖端(或附近)的应力有多大是没有意义的。有鉴於此,我们引入了一个参数:应力强度因子(stress intensity factor)。应力强度因子共有三个,一般写做KI、KII以及KIII,分别对应到三种不同的破坏模式。
在线弹性破坏力学(linear elastic fracture mechanics,简称LEFM)中,外力的大小与应力强度因子成正比关系。此外,应力强度因子也和几何参数(例如裂纹长度、外力与裂纹的距离等)有 关。对於同一结构,应力强度因子可视为含裂纹结构所受外力大小的一个指标;换言之,应力强度因子越大,结构越危险。当应力强度因子超过「破坏韧性 (fracture toughness),记为Kc」,裂纹开始成长。破坏韧性通常由制作标准试片由实验求得,不同的材料有不同的破坏韧性,所以破坏韧性可视为材料性质。
至於裂缝会往那个方向成长?有好几个理论可以预测,例如能量释放率理论(energy release rate)、最大周向应力理论(maximum circumferential stress)、J积分理论(J-integral)、应变能密度理论(strain energy density theory)。以上几个理论,除了最大周向应力之外,基本上都和应变能有关,这个部分较为深入,不详述了。
接着我们要探讨有限元素法在破坏力学的应用。前面讲过,探讨裂缝尖端附近的应力有多大是没有意义的。因此我们的重点在於,如何运用有限元素法求得裂纹尖端的应力强度因子。最常见的就是利用「四分之一节点(quarter point)」元素来模拟裂纹尖端,如图二。严格来说,四分之一节点元素不是「三角形元素」,它是由四边形元素退化而成,退化方式如图三。一般我们利用节点上的位移,求得裂纹尖端的开口位移(crack tip opening displacement)来反推得到应力强度因子。MARC及ANSYS皆用四分之一节点元素来算应力强度因子。
以上所说的四分之一节点元素,其最原始的统御方程式和弹性力学并没有两样,只是将其中几个节点变了位置,因此形状函数也跟着改变,最後得到一个结果:这种元素内部的应力呈现根号r分之一的奇异性,因此可用来模拟裂纹。
另外有些软体,例如NASTRAN,则利用另一种特殊元素,共有18个节点,如图四。这种元素最原始的统御方程式,和一般的弹性力学稍有不同,它一开始就假设应力正比於根号 r分之一。
不管如何,裂纹尖端的元素与一般的元素不同,裂纹尖端网格的大小,由算出来的应力强度因子准不准来决定。如何知道准不准?可建立一个简单模型,有理论解的 题目,即可比较,例如图四。以平面问题为例,若采用四分之一节点元素,元素边长大约在裂纹长度的百分之六左右;若采用NASTRAN的特殊元素,其长度更 只要裂纹长度的百分之十至二十即可。当然以上的准则并非万用,不同的题目最好还是多试几个网格粗细,再做决定。
方向
元素有方向性。有些元素如果方向搞错,则跑不出来;有些虽然跑得出来,但结果却有问题。 以平面元素(包括平面应变、平面应力以及轴对称)为例,四个节点的编号必须是「逆时钟方向」。如果是顺时钟方向,则在有限元素的定义中,这个元素的面积小於零,就跑不出来了。
板(或壳)元素也有方向性。板元素的方向由右手定则决定,也就是元素的正向方向,由元素编号顺序的方向决定,如此也定义了板元素的上表面或下表面。当板受 到bending时,上表面和下表面分别会受到张应力以及压应力(当然也可能颠倒,视bending的方向而定)。图一为两个相邻的板元素受到 bending作用,元素1的方向和元素2的方向不同,因此它们的上、下表面也不一致(参考图一)。在受到图一所示的bending时,元素1的上表面为 张应力,下表面为压应力;元素2的上表面则为压应力,下表面为张应力。有些有限元素软体(如MARC)在计算板元素的「上表面」节点应力时,是根据两个元 素的「上表面」高斯点上应力经外差到节点上再平均求得,但注意图一中两个元素的上下表面定义颠倒,所以平均後的上表面节点应力就很不准了。
再解释清楚一点,我们先将两个元素分开,同时画上厚度,比较清楚,如图二,图一中的A点在图二则可分为A1至A4四点。如果单独看元素1,A1点应该受到 张应力,A2为压应力;单独看元素2,A3为张应力,A4为压应力。当有限元素法在计算A点「上表面」应力时,是将元素1以及元素2的「上表面」应力平均 (即将A1及A4的应力平均),本来应该上表面是张应力的,被这麽一搞,应力下降,误差就来了。
一般有限元素法软体,在做automesh时,元素的方向与surface(或area)的方向一致,因此同一个surface(或area)做出来的网 格,方向都相同。但如果两个surface(或area)相邻,就要特别注意方向有没有一致。如果是「手动」建立网格,更要特别注意建立出来的元素的方 向,以免造成困扰。
形状函数
形状函数在有限元素法中是非常重要的一个概念,它定义了元素内部位移的分布。以一维线性元素为例(图一) 此种元素有两个节点,分别以节点1及2表示。节点1的座标为xi = -1,节点2则为xi = -2。元素中有几个节点,就有几个形状函数。因此我们有两个形状函数。形状函数有个特点:考虑第n个形状函数,若代入第n个节点的座标,其函数值为1;若 代入其他节点的座标,则函数值为0。由於图一为两节点元素,所以形状函数为线性,其表示式为:其他种类的元素的形状函数,一般有限元素法书籍(或商用有限 元素法软体的使用手册)皆有提到,这里不多写了。
元素的特性主要由形状函数所主宰。例如三节点之三角形元素,沿着元素内部任一方向,位移皆为线性分布。为什麽我们知道它是线性分布?只要查查该元素的形状 函数立刻就明白了。由於应变为位移的对空间的一次微分,所以在三节点三角形元素内部,应变是保持不变的,因此该元素又称为等应变元素(constant strain triangle)。应变既然不变,应力也不变。同样道理,四节点金字塔3D元素,应力及应变也是保持不变的。而四边形四
热心网友
时间:2023-10-15 13:49
ANSYS两个后处理器:
★通用后处理器POST1:查看整个模型在各个时间点上的结果
★时间历程后处理器POST26:查看整个模型上的某一点结果随时间的变化曲线。后处理可在求解完后直接进入,也可在重新进入ANSYS后读入文件进入后处理。