ANSYS复习指南.docx
- 文档编号:29414821
- 上传时间:2023-07-23
- 格式:DOCX
- 页数:31
- 大小:163.32KB
ANSYS复习指南.docx
《ANSYS复习指南.docx》由会员分享,可在线阅读,更多相关《ANSYS复习指南.docx(31页珍藏版)》请在冰豆网上搜索。
ANSYS复习指南
有限元复习指南
欧拉方程
归纳法可论证:
则欧拉方程转化为常系数线性方程:
耦合场分析
耦合场分析的定义
耦合场分析是指在有限元分析的过程中考虑了两种或者多种工程学科物理场的交叉作用和相互影响耦合。
例如,压电分析考虑了结构和电场的相互作用,它主要解决由于所施加的位移载荷引起的电压分布问题,反之亦然;其他的耦合场分析还有热-应力耦合分析、热-电耦合分析、流体-结构耦合分析、磁-热耦合分析和磁-结构耦合分析等等。
耦合场分析的过程取决于所需解决的问题是由哪些场的耦合作用,但是耦合场的分析最终可归结为两种不同的方法“序贯耦合方法和直接耦合方法”:
序贯耦合解法是按照顺序进行两次或更多次的相关场分析,它是通过把第一次场分析的结果作为第二次场分析的载荷来实现两种场的耦合的,例如序贯热-应力耦合分析是将热分析得到的节点温度作为体力载荷施加在后序的应力分析中来实现耦合的;
直接耦合解法利用包含所有必须自由度的耦合单元类型,仅仅通过一次求解就能得出耦合场分析结果。
在这种情形下耦合是通过计算包含所有必须项的单元矩阵或单元载荷向量来实现的,例如利用单元SOLID5、PLANE13或SOLID98可直接进行压电分析。
不存在高度非线性:
序贯耦合方法
存在高度非线性:
直接耦合方法
典型ANSYS顺序耦合分析应用包括:
✓热应力
✓感应加热
✓感应搅拌
✓稳态流体-结构耦合
✓磁-结构耦合
✓静电-结构耦合
✓电流传导-静磁
什么是物理环境?
单元类型及KEYOPT设定;
实常数;
材料属性;
单元坐标系
求解分析选项;
载荷步选项;
约束方程;
耦合节点集;
施加的边界条件和载荷;
GUI过滤设置;
分析标题(/TITLE);
进行顺序耦合场分析可使用间接法或物理环境法。
在直接耦合场分析中,只需用耦合场单元进行一次分析。
表4-1具有耦合场分析能力的单元
单元名称
描述
SOLID5
耦合场六面体
PLANE13
耦合场四边形
FLUID29
声学四边形
FLUID30
声学六面体
CONTAC48
2-D点对面接触
CONTAC49
3-D点对面接触
CONTA171
2-D面对面接触
CONTA172
2-D面对面接触
CONTA173
3-D面对面接触
CONTA174
3-D面对面接触
SOLID62
3-D磁-结构单元
FLUID116
热-流体管道单元
PLANE67
热电四边形单元
LINK68
热电线单元
SOLID69
热电六面体单元
SOLID98
耦合场四面体单元
CIRCU124
通用电路单元
SHELL157
热电壳单元
TRANS126
机-电换能器
什么是直接耦合场分析
耦合场单元包含所有必要的自由度,通过计算适当的单元矩阵(矩阵耦合)或是单元载荷矢量(载荷矢量耦合)来实现场的耦合在用矩阵耦合方法计算的线性问题中,通过一次迭代即可完成耦合场相互作用的计算,而载荷矢量耦合方法在完成一次耦合响应中至少需要二次迭代。
对于非线性问题,矩阵方法和载荷矢量耦合方法均需迭代。
×记忆要点:
+耦合中的自由度方向(UX,UY,等)是节点坐标系中的方向.
+求解器只保留耦合中的第一个自由度,并把它作为主自由度,而不保留其余自由度.
+施加在耦合节点上的载荷(在耦合自由度方向)求和后作用在主节点上.
+耦合自由度上的约束只能施加在主节点上.
约束方程一般应用于:
×连接不同的网格
×连接不同类型的单元
×建立刚性区
×过盈装配
APDL语言
APDL是ANSYS参数化设计语言的的缩写,它是一种允许使用参数并能完成一系列任务的强大的程序语言
一些命名规则:
×参数名不超过8个字符,并以字母开头。
×参数名中只能出现字母,数字和下划线。
×避免以下划线开头,这在ANSYS中另有它用。
×参数名不分大小写,如“RAD”和“Rad”是一样的。
所有的参数都以大写形式存储。
避免使用ANSYS标识,如STAT,DEFA,和ALL
记住:
你可以用参数化表达式和/或if-then-else分支来修改参数。
有关宏的注意事项:
×建立短小且简单的宏。
×当创建宏的时候,可以在ANSYS命令窗口中剪切并粘贴命令来检测这些命令的顺序是否正确。
×运用注释来描述意图或期望的执行结果。
×在自己的目录中保存你自己的宏。
×在ANSYS_MACROLIB环境变量目录中保存那些任何人都可以进入的宏
识别无效分析结果
应该知道所分析的对象的一些基本的行为:
•重力方向总是竖直向下的
•离心力总是沿径向向外的
•物体受热一般要膨胀
•没有一种材料能抵抗1,000,000psi的应力
•轴对称的物体几乎没有为零的环向应力
•弯曲载荷造成的应力使一侧受压,另一侧受拉
如果只有一个载荷施加在结构上,检验结果比较容易.如果有多个载荷,可单独施加一个或几个载荷分别检验,然后施加所有载荷检验分析结果
定义误差估计的基础
ANSYS通用后处理包含网格离散误差估计.
误差估计是依据沿单元内边界的应力或热流的不连续性,是平均与未平均节点应力间的差值
误差估计在同时符合以下情况有效:
•线性静力结构分析及线性稳态热分析
•大多数2-D或3-D实体或壳单元
不符合以上条件的分析,或者使用PowerGraphics时,ANSYS自动关闭误差估计
单元死活
×不与任何单元连接的节点会发生“漂移”。
在某些情况下,可能想约束不被激活的自由度,以减少要求解的方程数目或避免出现病态条件。
×单元被激活时,如果想维持其一定的形状,那么约束不被激活的自由度是很重要的。
但在重新激活单元时一定要删除这些人为约束。
×注意,约束方程(CE或CEINTF)不能用于不被激活的自由度。
对于死亡单元载荷矢量中的单元载荷是自动置零的,而对激活单元则自动恢复,质量矩阵中的质量也被置为零(加速度载荷对死亡单元无影响)
集中力不能从不被激活的自由度上自动删除(比如那些不与任何活的单元连接的节点)。
务必删除死节点上的集中载荷
有限元法是古典变分法与分片插值法相结合的产物。
它不是在全域范围内选取基函数,而是先将全域分成单元,在单元范围内用低次多项式分片插值,再将它们组合起来,形成全域内的函数,用以逼近问题的真解。
H方法和p方法的含义及区别
h-方法可用于任何类型的结构分析,而p-方法只能用于线性结构静力分析。
根据所求的问题,h-方法通常需要比p-方法更密的网格。
p-方法在应用较粗糙的网格时,提供了求得适当精度的一种很好的途径。
h-方法和p-方法是两种处理有限元单元的方法,他们之间的区别主要反映在提高计算精度的办法上,H方法是用简单的单元,但减小单元尺寸,即增加单元数(细化网格)来实现,而P方法则采用复杂的单元,如增加个体单元的节点,形成非线性单元,但保持单元尺寸不变(即不用细化)来提高计算精度。
P-method是通过提高形函数阶次来提高计算精度的,H-method则通过减小单元尺寸来提高精度。
两种方法都可以达到根据计算误差调整单元阶次或网格尺寸,实现自适应分析的目的。
目前,H-method相对成熟一些。
什么是泛函?
泛函数是以函数作为变元的函数(通俗地说,泛函数就是函数的函数)。
泛函数概念的产生与变分学问题的研究发展有密切关系。
有限元分析的发展阶段
•失效分析阶段
•设计验证分析阶段
•保证产品正常工作分析阶段
•提高产品性能分析阶段
•设计创新分析阶段
•产品寿命分析阶段
有限元分析的几大功能
•结构分析
•热分析
•电磁分析
•流体分析(CFD)
•耦合场分析-多物理场
有限元分析的基本步骤
1.创建有限元模型
–创建或读入几何模型.
–定义材料属性.
–划分单元(节点及单元).
2.施加载荷进行求解
–施加载荷及载荷选项.
–求解.
3.查看结果
–查看分析结果.
–检验结果.(分析是否正确)
–
–
ANSYS在分析过程中需要读写文件.
文件格式为jobname.ext,其中jobname是设定的工作文件名,ext是由ANSYS定义的扩展名,用于区分文件的用途和类型.
默认的工作文件名是file.
分析完成后,您必须保存如下文件:
log文件(.log),数据库文件(.db),结果文件(.rst,.rth等),载荷步文件(.s01,.s02,...),输出文件(.out),物理环境文件(.ph1,.ph2,...).
电磁场理论
B满足的边界条件
H满足的边界条件
磁通连续性)B的法向分量连续
H的切向分量连续,只有在跨越理想完纯导体或超导体的边界时,K≠0,它才会不连续
磁场四个区,.磁化曲线的非线性
起始磁导率区:
在外磁场的作用下,磁畴壁做可逆移动;
不可逆磁化区:
这是畴壁做不可逆跳跃式移动的过程,磁化强度随着磁场增大急剧增大;
旋转磁化区:
畴壁位移基本完毕,只有靠磁畴内磁矩的转动来进一步增大磁化强度;
接近饱和区:
磁畴内磁矩的可逆转动造成,此时增大磁场,磁化强度增加很小。
工程上常用的磁性材料通常分软磁材料和硬磁材料两类。
软磁材料:
磁性能的主要特点磁导率高,矫顽力低。
对某些软磁材料来说,还要求磁化损耗尽可能低;
硬磁材科:
磁性能的主要特点是矫顽力高,剩磁感应强度和磁能积也是其最重要的参数;
磁场使用单元场合
磁场分析假设和限制
Ø单元面积必须大于零;
Ø在做轴对称分析时,在全局坐标系下的Y轴必须是对称轴(轴向),而且模型必须在右半平面(X>0);
Ø当单元只有AZ自由度(KEYOPT
(1)=0),在做谐波或瞬态分析时,如果用BFE,,JS命令施加了电流密度载荷,那么该单元表征一个绞线区域,没有施加电流密度载荷时,则表征一个包含涡流效应的实心导体。
Ø谐波分析时不允许有永磁体
Ø在做静磁场(magnetostatic)分析时,不能有VOLT,AZ自由度
Ø对于谐波和瞬态分析,不支持速度和电路耦合效应
Ø当用电压源(KEYOPT
(1)=2)或电路(KEYOPT
(1)=3/4)驱动时:
◆必须用MKS单位制
◆磁导率和电导率是定常和各向同性的
◆线圈区域所有的CURR和EMF自由度必须耦合在一起
Ø在做电路耦合瞬态分析时,若用TINTP命令指定反向欧拉算法的参数时,使用默认值THETA=1.0。
Ø如果同一材料同时给定B-H曲线和ur,那么使用ur
Ø可以指定各向异性相对磁导率,MURX,MURY,MURZ
Ø可以同时用ur和B-H曲线指定各方向磁导率,若ur=0,则使用B-H曲线
Ø一般采用Sparse或Frontalsolver
Ø模型很大是,可以采用JCG或PCGsolver
Ø如果是电压激励或包含速度效应时,只能用sparse,frontal,JCG,ICCGsolver
Ø电路激励时,只能用Sparse或Frontalsolver
径向磁轴承转子涡流分析
原因:
对于径向磁轴承,当转子高速旋转时,内部的交变磁场会感应出涡流,由于集肤效应,磁场趋于转子表面的浅层区域内并造成磁饱和,从而影响气隙磁密分布,减小轴承的电磁力,并最终影响系统的性能。
目的:
通过有限元分析研究转子中涡流、磁密随转速的变化,并计算轴承的电磁力,最后计算磁轴承的电流刚度,从而将涡流的影响包括进来。
分析方法:
Ø静态分析:
转子损耗,静态电磁力
Ø时谐分析:
无实际意义,仅作为例子说明方法
Ø瞬态分析:
转子损耗,动态电磁力
避免出现矢量位函数A对空间坐标的偏导数。
实现这一目的有以下三种技术:
Ø边界积分法
Ø重新剖分法
Ø独立坐标系法
转子边界和定子边界耦合方法
虚节点法,插值法
什么是模态分析?
•模态分析是用来确定结构的振动特性的一种技术:
–自然频率
–振型
–振型参与系数(即在特定方向上某个振型在多大程度上
参与了振动)
模态分析是所有动力学分析类型的最基础的内容
模态分析的好处:
•使结构设计避免共振或以特定频率进行振动(例如扬声器);
•使工程师可以认识到结构对于不同类型的动力载荷是如何响应的;
•有助于在其它动力分析中估算求解控制参数(如时间步长)。
建议:
由于结构的振动特性决定结构对于各种动力载荷的响应情
况,所以在准备进行其它动力分析之前首先要进行模态分
析。
•在ANSYS中有以下几种提取模态的方法:
–BlockLanczos法
–子空间法
–PowerDynamics法
–缩减法
–不对称法
–阻尼法
模态分析四个步骤
•建模
•选择分析类型和选项
•施加边界条件并求解
•观察结果
什么是有预应力的模态分析?
为什么要做有预应力的模态分析?
•具有预应力结构的模态分析;
•同样的结构在不同的应力状态下表现出不同的动力特性。
–例如,一根琴弦随着拉力的增加,它的振动频率也随之增大。
–涡轮叶片旋转时,由于离心力引起的预应力的作用,它的自然频率逐渐具有增大的趋势。
为了恰当地设计这些结构,必须要做具有预应力和无预应力的模型的模态分析
有预应力的模态分析
•建模
•在静态分析中给模型施加预应力
•做具有预应力的模态分析
什么是循环对称结构的模态分析?
•利用循环对称的模态分析;
•可以只模拟结构的一个扇形区,然后观察整个结构的振型。
–节省了建模时间-不需要模拟整个结构。
–节省了计算时间和硬盘空间-只需要较少的单元和自由度。
应用:
可用于任何具有循环对称的结构:
如涡轮、叶轮
循环对称结构的模态分析步骤
•基本扇区的建模
•确定循环对称平面
•复制一个基本扇区
•在两个扇区上施加边界条件
•指定分析类型和选项
•用CYCSOL命令求解
•将求解结果扩展到3600,对结果进行评价
为什么节径范围很重要?
•由于只模拟了一个基本扇区,所以ANSYS需要知道将要提取哪些振型。
是提取对某一给定节径的所有振型还是提取所给节径范围内的前几阶振型?
•结构的低阶振型通常是前几节径的前几阶振型;
•通常,只需对前面少数几条节径提取少数几阶振型。
模态分析可分为实验模态分析与计算模态分析两种方法
一般来说,稳态分析中网格上结点温度比实际温度要低。
也就是说,如果加密网格,温度将增加,但加密到一定程度,结果将不显著增加(i.e.,结果收敛)。
热分析
如何使热传递分析包括非线性?
传热问题可以是非线性的,因为材料特性受温度或因边界条件是非线性的。
通常情况下,性能材料非线性与温度相关的是轻度的,该属性不随温度的变化迅速。
然而,当潜热效应计算在内,分析可能是严重的非线性的。
何时需要定义比热和密度?
•瞬态问题,这些数值用于形成比热矩阵(该矩阵表示瞬态分析中需要的热能存储效果).
•稳态分析中包括有热质量传递效果(i.e.,模型中有流动导体介质).
单元类型
控制单元(恒温器),谐波单元(线性)热流单元,热平面效果单元
热载荷和边界条件注意事项
–在ANSYS中,没有施加载荷的边界作为完全绝热处理。
–对称边界条件的施加是使边界绝热得到的。
–如果模型某一区域的温度已知,就可以固定为该数值。
–响应热流率只在固定温度自由度时使用
稳态热分析中关于材料特性的总体说明
–对于稳态分析,热材料特性必须输入热传导率“k”-KXX,和可选的KYY,KZZ。
–如果用户不定义,KYY和KZZ缺省等于KXX。
–密度(DENS)和比热(C)或热焓(ENTH)在没有质量传递的稳态热分析中不需要。
–随温度变化的材料导热系数k,使得热分析为非线性。
–与温度有关的换热系数也被处理为材料特性。
迭代(快速求解)选项
快速求解选项可以用于除了以下情况的任意线性,非线性,稳态或瞬态热分析:
–可能不能用于旧的平面效果单元,SURF19和SURF22。
–可能不能用于辐射计算。
–不推荐用于带相变的热传递问题。
–温度偏移是可选的,但在下列情况下必须使用:
–有辐射效果,并且使用了°F或°C。
–使用了随温度变化的热生成率(MASS71)。
要使用载荷步方法:
1.指定分析控制,指定第一载荷步的载荷和选项,改变标题。
2.写载荷步文件。
3.改变为下一个载荷步中的载荷和选项,改变标题。
4.写下一个载荷步文件。
5.对于剩余的载荷步重复步骤3和4。
6.求解载荷步文件。
7.后处理
辐射是一种非线性现象,因此需要进行平衡迭代来得到收敛解。
非线性线单元LINK31用于计算两节点间或多对节点间的简单辐射热传递,节点的位置是任意的,可作为其它单元的节点。
LINK31需要定义如下数据:
材料属性:
EMIS辐射率(可以随温度变化)
实常数:
AREA(Ai)(有效辐射面积)
FORMF(Fij)(形状系数)
SBCONST(Stefan-Boltzman常数
使用AUX12计算面与面之间的辐射可分为以下三个步骤:
(1)定义辐射面
(2)生成辐射矩阵
(3)将辐射矩阵用于热分析
由于辐射表面不支持对称条件,包含辐射的模型就无法利用几何上的对称性,因此必须建立完整的分析模型。
•无论是隐藏的方法还是非隐藏的方法,通常辐射表面的网格越细,越规则,形状系数计算精度越高。
但是,对于隐藏法而言,如果要得到相同精度的形状系数,其对网格的要求比非隐藏法更高。
如果网格太差,即使将射线的数量增加到其最大值,也无法获得所需的求解精度。
•对于轴对称情况,NDIV设为20,可以得到足够精度的形状系数。
单元在拓展到3D时应有合理的形状(长细比应在合理的范围内)。
•用于生成2D辐射矩阵的LINK32单元,并不直接支持轴对称选项。
因此,对于轴对称模型,确认在运行热分析以前删除或不选择此单元。
•理论上讲,对于封闭系统,由任意一个辐射表面到所有其它辐射面的形状系数的和为1;对于开放系统则应小于1。
可以通过执行MPRINT,1命令将形状系数如下打印出来“***FORMFACTORS***TOTAL=Value”,由此可检查每一个辐射面形状系数的计算是否正确。
如果超过1则肯定错误。
尤其在两辐射面间有障碍时,不留意地使用了非隐藏方法计算,就会出现这种情况。
•只有所有的辐射面之间可以完全地看到对方时,才能使用非隐藏方法。
否则形状系数的计算是错误的,热分析的结果不正确甚至不会收敛。
•隐藏方法需要明显更长的计算时间,所以只有辐射面间有障碍存在或无法分组计算时才选用。
•对于有些情况可以对辐射面分组,各组之间在辐射传热上是完全独立的。
由于在一个组中的辐射面之间没有阻碍,可以用非隐藏方法计算形状系数,分别写入辐射矩阵文件。
这样可以节省大量CPU时间。
要对辐射面进行分组,在写矩阵之前选择的需要的辐射面组。
•对于隐藏方法,增大射线数量会提高形状系数的计算精度。
只可以描述系统的“焓”变了多少,而不能测出系统的“焓”具体是多少。
焓变>0系统经历了吸热反应。
焓变<0系统经历了放热反应。
进行相变分析时,使用:
–打开时间积分的瞬态分析。
–时间步初始数值较小,时间步也很小。
–自动时间步。
–低阶单元类型(PLANE55或SOLID70)。
–如果选择的高阶单元,打开对角比热矩阵选项。
如果载荷在这个载荷步是恒定的,需要设为阶越选项;如果载荷值随时间线性变化,则要设定为渐变选项:
引起结构非线性的原因有很多,它们可分成以下三种主要类型:
1.几何非线性
大应变,大位移,大旋转
2.材料非线性
塑性,超弹性,粘弹性,蠕变
3.状态改变非线性
接触,单元死活
非线性分析得到的结果
不能使用叠加原理!
结构响应与路径有关,也就是说加载的顺序可能是重要的。
结构响应与施加的载荷可能不成比例。
一些典型的非线性分析的应用包括:
非线性屈曲失稳分析
金属成形研究
碰撞与冲击分析
制造过程分析(装配、部件接触等)
材料非线性分析(弹性材料、聚合物)
ANSYS缺省的收敛判据是力/力矩和位移/旋转增量。
对于力/力矩缺省的容限是0.5%,对于位移/旋转增量的容限是5%。
经验表明这些容限对于大多数问题具有足够的精确度。
缺省的设置对于广泛的工程问题既不“太紧”也不“太松”。
为得到平方的收敛速度,切向刚度矩阵需要是全一致的。
切向刚度矩阵[KT]由四部分组成:
[KT]=[Kinc]+[Ku]+[K]-[Ka]
这里[Kinc]=主切向刚度矩阵
[Ku]=初始位移矩阵
[K]=初始应力矩阵
[Ka]=初始载荷矩阵
对于非线性分析,有三种方程求解器供选择:
Sparse(直接求解器,求解控制缺省设置)
波前求解器(直接求解器)
PCG(迭代求解器)
线性搜索
线性搜索是一个非常强大的改进收敛工具。
激活后它并不会降低求解的稳定性,而且在许多情况下,它可以改进收敛缓慢的求解。
但是它会需要增加一些CPU时间来计算线性搜索参数。
线性搜索是一种非常有效的改进求解振荡问题收敛情况的方法。
如果看到输出窗口中的MAXDOFINC在正负值间振荡,激活线性搜索!
如果非线性响应平稳(而且时间步长适当小),预测可以加速收敛。
如果非线性响应不平稳,或分析中涉及大旋转,预测可能导致发散!
对大旋转分析不要使用预测。
自适应下降不能与线性搜索合用!
这两个选项相互排斥。
下面列出了完成非线性分析所需的典型步骤:
1.指定分析类型
2.指定几何非线性打开或关闭
3.为载荷步指定“时间”
4.用NSUBST或DELTIM设定子步数
5.施加载荷与边界条件
6.指定输出控制与监视值
7.保存数据库
8.求解载荷步
如果收敛顺利,自动时间步长控制将增加载荷增量。
如果收敛困难,自动时间步长控制将二分或减少载荷增量。
变形体几何形态的改变将明显影响物体的载荷-位移(如刚度)特性。
几何非线性并不只是指大位移,而且还包括几何状态改变所引起的任何结构响应的变化。
它包括大应变、大位移和大旋转。
大应变分析设定应变不再是无限小的,而是有限的或相当大的
当应变超过一定百分比及不能忽视几何形状的改变时,可认为是大应变。
大应变理论考虑了形状的改变(例如厚度,面积等等)及任何大旋转。
注意载荷方向
×在大应变分析中预测网格扭曲,划分适当的网格。
参考预测网格扭曲的指南。
×避免过分约束边界处的变形
×避免使用带中间节点的单元
×使用适当的单元类型和积分准则以解决网格自锁问题。
(在单元选择一章中有相关的更详细信息。
)
×时间步大小应控制在每个子步中的最大旋转度数小于5或10度。
×对大旋转分析不要使用预测。
×梁单元和壳单元使用足够的网格密度;没有一个单元可承受超过30度的弯曲度。
×如果自动
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- ANSYS 复习 指南