ansys高级非线性分析四蠕变
蠕变分析【精选文档】
4。
4 蠕变分析4.4.1 蠕变理论4.4.1。
1 定义蠕变是率相关材料非线性,即在常荷载作用下,材料连续变形的特性。
相反如果位移固定,反力或应力将随时间而变小,这种特性有时也称为应力松驰,见图4—18a .图4-18 应力松弛和蠕变蠕变的三个阶段如图4-18b所示.在初始蠕变阶段,应变率随时间而减小,这个阶段一般发生在一个相当短的时期。
在第二期蠕变阶段,有一个常应变率,所以应变以常速率发展,在第三期蠕变阶段,应变率迅速增加直到材料失效.由于第三期蠕变阶段所经历的时间很短,材料将失效,所以通常情况下,我们感兴趣的是初始蠕变和第二期蠕变。
ANSYS程序中的蠕变行为用来模拟初始蠕变和第二期蠕变。
蠕变系数可以是应力、应变、温度、时间或其它变量的函数。
在高温应力分析中(如核反应堆等),蠕变分析非常重要。
例如,假设在核反应堆中施加了预荷载,以保证与相邻部件保持接触而不松开。
在高温下过了一段时间后,预荷载将降低(应力松驰),可能使接触部件松开。
对于一些材料如预应力砼,蠕变也可能十分重要。
最重要的是要记住,蠕变是永久变形。
4。
4。
1。
2 理论介绍蠕变方程:我们通过一个方程来模拟蠕变行为,此方程描述了在实验中观测到的主要特征(特别是在一维的拉伸实验中)。
这个方程以蠕应变率的方式表示出来,其形式如下:上式中,A、B、C、D是从实验中得到的材料常数,常数本身也可能是应力,应变,时间或温度的函数,这种形式的方程被称为状态方程。
上式中,当常数D为负值时,蠕应变率随时间下降,材料处于初始蠕变阶段,当D为0时,蠕应变率为常值,材料处于第二期蠕变阶段。
对于2-D或3-D应力状态,使用VON Mises方程计算蠕应变率方程中所使用的标量等效应力和等效应变。
对蠕变方程积分时,我们使用经过修改的总应变,其表达式为:经过修改的等效总应变为:其等效应力由下式算出:其中:G=剪切模量=等效蠕应变增量由程序给出的某一种公式进行计算,一般为正值,如果在数据表中,则使用的是衰减的蠕应变率而不是常蠕变率,但这个选项一般不被推荐,因为在初始蠕变所产生的应力为主的情况下,它可能会严重的低估蠕变值.如果,程序使用修正的等效蠕应变增量来代替蠕应变增量。
ANSYS本构关系
2021/10/10
5
(4)多线性等向强化(MISO) 使用多线性来表示使用Von Mises屈服准则
的等向强化的应力-应变曲线,它适用于比例加 载的情况和大应变分析。 需要输入最多100个应力-应变曲线,最多可以 定义20条不同温度下的曲线。 其材料特性的定义步骤如下: 1.定义弹性模量 2.定义MISO数据表 3.为输入的应力-应变数据指定温度值 4.输入应力-应变数据 5.画材料的应力-应变曲线 与MKIN 数据表不同的是,MISO的数据表对不同 的温度可以有不同的应变值,因此,每条温度曲 线有它自己的输入表。
率,弹性斜率和塑性斜率,由于随动强化的Von mises屈 服准则被使用,所以包含有鲍辛格效应,此选项适用于遵 守Von Mises 屈服准则,初始为各向同性材料的小应变问 题,这包括大多数的金属。
需要输入的常数是屈服应力和切向斜率,可以定义高 达六条不同温度下的曲线。 注意: 1.使用MP命令来定义弹性模量 2.弹性模量也可以是与温度相关的 3.切向斜率Et不可以是负数,也不能大于弹性模量 在使用经典的双线性随动强化时,可以分下面三步来定义 材料特性。 1.定义弹性模量 2.激活双线性随动强化选项 3.使用数据表来定义非线性特性
2021/10/10
3
(2)双线性等向强化(BIS0)
也是使用双线性来表示应力-应变曲线,在此选 项中,等向强化的Von Mises 屈服准则被使用, 这个选项一般用于初始各向同性材料的大应变问 题。需要输入的常数与BKIN选项相同。 举例如下: MP,EX,1,200e9 MP,NUXY,1,0.25 MP,GXY,1,150e9 TB,BISO,1
变关系,而开裂和压碎后采用其给出的本构关系。但屈服准则尚可另
Ansys Workbench蠕变分析
Workbench蠕变分析
1.问题描述:一方形实体,尺寸为5mX1mX1m,在200MPa压力,温度分别为20℃、100℃、500℃下蠕变分析。
边界条件如下图所示。
2.新建static structure模块,如图所示,分别对应20℃、100℃、500℃。
3.双击Engineering Data,定义材料蠕变特性,本示例使用Time Hardening模型(关于蠕变的模型选择请自行补脑)。
4.关闭Engineering Data,进入到Mechnical Systems,划分网格,施加边界条件。
5.对于A5、B5、C5温度分别为22℃、100℃、500℃。
6.对A5模块定义载荷步,分两次加载(示例仅计算9秒的蠕变).载荷步详细见下图。
7.同样的,对B5、C5做相同的载荷步设置。
并进行求解。
完成后如下图所示。
上图图中显示的是22℃下的蠕变。
在200MPa载荷下,变形为2.5E-3m(第一个载荷步的结果,时间点为1秒),保持载荷不变,持续9秒,变形达到0.034m(第二个载荷步,时间点为10s)。
100℃、500℃的蠕变如下图所示
100℃
500℃
PS:本例重在分享Workbench中进行蠕变分析的过程,请忽略其结果的合理性。
By问道真人
2015/11/27。
ansys几何非线性+塑性+接触+蠕变
除去蠕变,这个模型的结果可靠性是不错的。
作了一系列接触问题,通过试验验证符合的很好。
模型解释:(1)一个弹性结构受压(接触)变形,到发生塑性变形。
(2)拿开压缩板,结构回弹,但不会回到原始位置。
(3)这时计算蠕变,释放掉应力。
(4)再压弹性结构到开始压缩位置。
比较这四步的接触力。
结果:第二,三步当然没有接触力,(若没有应力释放,第一、第四步接触力应一样,)有了应力释放,第四步接触力比第一步减小。
这个模型中的蠕变没用太好。
用的是隐式6号蠕变方程,蠕变是时间和应力的函数,参数是乱定的(应力释放太快)。
想请教有关蠕变方面的资料,尤其是材料蠕变方程选用及参数方面的资料。
/prep7!------------CuSn8----------ET,1,182,,,3mp,ex,1,115e9mp,prxy,1,0.3r,1,0.3TB,BKIN,1TBDA TA,1,470E6,0tm=100*SET,C1,1.5625E-14 !ASSIGN VALUE*SET,C2,1.5 !ASSIGN V ALUE*SET,C3, !ASSIGN V ALUE*SET,C4,0 !ASSIGN V ALUETB,CREEP,1,,,6 !ACTIV ATE DA TA TABLETBDA TA,1,C1,C2,C3,C4 !DEFINE DATA FOR TABLE!-----------contact-----------------ET,9,169ET,10,171R,9,,,0.1,0.1,,!RMORE,,,1.0E20,0.0,1.0,!RMORE,0.0,0,1.0,0,0,0.5!RMORE,,,1.0,0.0MP,MU,9,0.0!----------------geometryk,,2k,,2,0.2k,,,0.2k,,-0.2k,,-0.2,1.2k,,,1k,,2,1.2k,,1,1k,,1.25,1k,,2,1L,8,9,k,,1.5,1.2k,,1.75,1.45L, 1, 2 L, 1, 4 L, 4, 5 L, 5, 11larc,7,12,11,0.25 larc,11,12,7,0.25L, 7, 10 L, 10, 9 L, 8, 6 L, 6, 3 L, 3, 2LFILLT,11,10,0.3, ,!*LFILLT,4,5,0.5, ,!*LFILLT,11,12,0.3, ,!*LFILLT,4,3,0.5, ,FLST,2,16,4 FITEM,2,12 FITEM,2,15 FITEM,2,11 FITEM,2,13FITEM,2,10FITEM,2,1FITEM,2,9FITEM,2,8FITEM,2,7FITEM,2,6FITEM,2,5FITEM,2,14FITEM,2,4FITEM,2,16FITEM,2,3FITEM,2,2AL,P51Xrect,1,3,1.45+0.001,1.5type,1mat,1esize,0.05amesh,all!---------contact------------ allstype,10mat,9real,9lsel,s,,,6,7nsll,s,1esln,s,0esurf,alltype,9mat,9real,9lsel,s,,,17nsll,s,1esln,s,0esurf,all!------boundarylsel,s,,,3nsll,,1d,all,uxd,all,uylsel,s,,,19nsll,,1cp,11,uy,allcplgen,11,ux*get,nmin,node,,num,min d,nmin,uxksel,s,,,10nslk*get,ndis,node,,num,minfini/soluantype,staticnlgeom,onautots,onallssaverate,offtime,1e-8d,nmin,uy,-0.3nsub,20outres,all,allsolve*get,rf1,node,nmin,rf,fy *get,dis1,node,ndis,u,ytime,2e-8d,nmin,uy,0.0nsub,20outres,all,allsolve*get,rf2,node,nmin,rf,fy*get,dis2,node,ndis,u,y!BFUNIF,TEMP,90rate,onTIME,tm!NSUBST,10OUTPR,BASIC,10 ! PRINT BASIC SOLUTION FOR EVERY 10TH SUBSTEP OUTRES,ESOL,1 ! STORE ELEMENT SOLUTION FOR EVERY SUBSTEP SOLVE*get,rf3,node,nmin,rf,fy*get,dis3,node,ndis,u,yrate,offtime,tm+1e-8d,nmin,uy,-0.3nsub,20outres,all,allsolve*get,rf4,node,nmin,rf,fy*get,dis4,node,ndis,u,y/EOFtime,11d,nmin,uy,-0.0nsub,20outres,all,allsolve*get,rf11,node,nmin,rf,fy*get,dis11,node,ndis,u,y/eoffini/post1*get,rf2,node,nmin,rf,fy fini/eof。
ansys功能及其简介
ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件。
由世界上最大的有限元分析软件公司之一的美国ANSYS开发,它能与多数CAD软件接口,实现数据的共享和交换,如Pro/Engineer, NASTRAN, Alogor, I-DEAS, AutoCAD等,是现代产品设计中的高级CAE 工具之一。
CAE的技术种类CAE的技术种类有很多,其中包括有限元法(FEM,即Finite Element Method),边界元法(BEM,即Boundary Element Method),有限差法(FDM,即Finite Difference Element Method)等。
每一种方法各有其应用的领域,而其中有限元法应用的领域越来越广,现已应用于结构力学、结构动力学、热力学、流体力学、电路学、电磁学等。
ANSYS有限元软件包是一个多用途的有限元法计算机设计程序,可以用来求解结构、流体、电力、电磁场及碰撞等问题。
因此它可应用于以下工业领域:航空航天、汽车工业、生物医学、桥梁、建筑、电子产品、重型机械、微机电系统、运动器械等。
软件主要包括三个部分:前处理模块,分析计算模块和后处理模块。
前处理模块提供了一个强大的实体建模及网格划分工具,用户可以方便地构造有限元模型;分析计算模块包括结构分析(可进行线性分析、非线性分析和高度非线性分析)、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析,可模拟多种物理介质的相互作用,具有灵敏度分析及优化分析能力;后处理模块可将计算结果以彩色等值线显示、梯度显示、矢量显示、粒子流迹显示、立体切片显示、透明及半透明显示(可看到结构内部)等图形方式显示出来,也可将计算结果以图表、曲线形式显示或输出。
软件提供了100种以上的单元类型,用来模拟工程中的各种结构和材料。
该软件有多种不同版本,可以运行在从个人机到大型机的多种计算机设备上,如PC,SGI,HP,SUN,DEC,IBM,CRAY等。
ANSYS讲义非线性分析
t1
t2
时间 t
XJTU
自动时间步(续)
• 自动时间步算法是 非线性求解控制 中包含的多种算法的一种。
(在以后的非线性求解控制中有进一步的讨论。) • 基于前一步的求解历史与问题的本质,自动时间步算法或者增加
或者减小子步的时间步大小。
XJTU
5) 输出文件的信息
在非线性求解过程中,输出窗口显示许多关于收敛的信息。输出 窗口包括:
子步
时间 ”相关联。
“时间
两个载荷步的求解 ”
XJTU
在非线性求解中的 “ 时间 ”
• 每个载荷步与子步都与 “ 时间 ”相关联。 子步 也叫时间步。
• 在率相关分析(蠕变,粘塑性)与瞬态分析中,“ 时间 ”代表真实 的时间。
• 对于率无关的静态分析,“ 时间 ” 表示加载次序。在静态分析中, “ 时间 ” 可设置为任何适当的值。
最终结果偏离平衡。
u 位移
XJTU
1) Newton-Raphson 法
ANSYS 使用Newton-Raphson平衡迭代法 克服了增量
求解的问题。 在每个载荷增量步结束时,平衡迭代驱 使解回到平衡状态。
载荷
F
4 3 2
1
u 位移
一个载荷增量中全 Newton-Raphson 迭代 求解。(四个迭代步如 图所示)
XJTU
非线性分析的应用(续)
宽翼悬臂梁的侧边扭转失 稳
一个由于几何非线性造 成的结构稳定性问题
XJTU
非线性分析的应用(续)
橡胶底密封 一个包含几何非线 性(大应变与大变 形),材料非线性 (橡胶),及状态 非线性(接触的例 子。
XJTU
非线性分析的应用(续)
Ansys软件简介及求解步骤
内容:ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件。
由世界上最大的有限元分析软件公司之一的美国ANSYS开发,它能与多数CAD软件接口,实现数据的共享和交换,如Pro/Engineer, NASTRAN, Alogor, I-DEAS, AutoCAD等,是现代产品设计中的高级CAD工具之一。
一、软件功能简介软件主要包括三个部分:前处理模块,分析计算模块和后处理模块。
前处理模块提供了一个强大的实体建模及网格划分工具,用户可以方便地构造有限元模型;分析计算模块包括结构分析(可进行线性分析、非线性分析和高度非线性分析)、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析,可模拟多种物理介质的相互作用,具有灵敏度分析及优化分析能力;后处理模块可将计算结果以彩色等值线显示、梯度显示、矢量显示、粒子流迹显示、立体切片显示、透明及半透明显示(可看到结构内部)等图形方式显示出来,也可将计算结果以图表、曲线形式显示或输出。
软件提供了100种以上的单元类型,用来模拟工程中的各种结构和材料。
该软件有多种不同版本,可以运行在从个人机到大型机的多种计算机设备上,如PC,SGI,HP,SUN,DEC,IBM,CRAY等。
目前版本为ANSYS5.7版,其微机版本要求的操作系统为Windows 95/98或Windows NT,也可运行于UNIX系统下。
微机版的基本硬件要求为:显示分辨率为1024×768,显示内存为2M以上,硬盘大于350M,推荐使用17英寸显示器。
启动ANSYS,进入欢迎画面以后,程序停留在开始平台。
从开始平台(主菜单)可以进入各处理模块:PREP7(通用前处理模块),SOLUTION(求解模块),POST1(通用后处理模块),POST26(时间历程后处理模块)。
ANSYS用户手册的全部内容都可以联机查阅。
用户的指令可以通过鼠标点击菜单项选取和执行,也可以在命令输入窗口通过键盘输入。
基于ANSYS的钢筋混凝土结构非线性有限元分析
2、应力-应变曲线:描述了混凝土和钢筋的在往复荷载作用下的变形和能量吸收能力,显示 了结构的塑性变形和损伤演化过程。
参考内容
引言
钢筋混凝土结构在建筑工程中具有重要地位,其非线性行为对结构性能影响 显著。因此,进行钢筋混凝土结构的非线性有限元分析对于预测结构响应、优化 结构设计具有实际意义。本次演示将根据输入的关键词和内容,建立钢筋混凝土 结构非线性有限元分析模型,并详细描述分析过程、结果及结论。
基于ANSYS的钢筋混凝土结构 非线性有限元分析
基本内容
引言:
钢筋混凝土结构是一种广泛应用于建筑工程的重要材料,其非线性力学行为 对结构设计的安全性和稳定性具有重要影响。为了精确模拟钢筋混凝土结构的真 实行为,需要借助先进的数值计算方法,如非线性有限元分析。ANSYS作为一种 广泛使用的有限元分析软件,为钢筋混凝土结构的非线性分析提供了强大的支持。
对于钢筋混凝土,其非线性行为主要来自两个方面:混凝土的本构关系和钢 筋与混凝土之间的相互作用。在非线性有限元分析中,需要建立合适的模型来描 述这些行为。例如,可以采用各向异性本构模型来描述钢筋混凝土的力学行为, 该模型可以捕捉到材料在不同主应力方向上的不同响应。
二、ANSYS中混凝土本构关系研 究
在进行荷载试验时,通过施加不同大小和方向的荷载,检测结构的变形和破 坏过程。采用静力荷载试验和动力荷载试验两种方式,分别模拟实际结构在不同 荷载条件下的响应。在试验过程中,记录各阶段的位移、应变和荷载数据。
在进行有限元分析时,采用ANSYS软件对试验数据进行模拟分析。首先进行 模态分解,了解结构的基本振动特性。随后进行屈曲分析,预测结构的失稳趋势。 通过调整模型参数和网格划分,对比分析不同方案下的有限元计算结果,为结构 的优化设计提供依据。
ansys高级非线性分析四蠕变
ANSYS高级非线性分析四蠕变简介蠕变是物质受力作用下的长时间形变现象。
在工程领域中,蠕变是一个重要的现象,因为它可能导致结构疲劳、失效以及安全问题。
在工程设计中,了解和预测材料蠕变行为是非常重要的,特别是在高温环境下。
在这方面,ANSYS提供了强大的高级非线性分析工具,可以模拟和预测材料的蠕变行为。
ANSYS高级非线性分析ANSYS是一款流行的有限元分析软件,广泛应用于工程领域。
在高级非线性分析中,ANSYS可以模拟材料的非线性行为,包括蠕变。
蠕变分析是一种长时间的形变分析,可以用来评估材料和结构在高温条件下的稳定性。
ANSYS的高级非线性分析功能可以通过几个步骤来实现。
首先,需要定义材料的力学性质,包括弹性模量、屈服强度和蠕变参数等。
其次,需要定义结构的几何形状和边界条件。
然后,可以设置加载条件,包括时间和温度的变化。
最后,可以进行蠕变分析,并输出结果。
四蠕变四蠕变是一种特殊的蠕变现象,出现在某些金属材料中。
四蠕变是指材料在连续加载下显示出四个明显的蠕变阶段,其形变速率逐渐增加,并最终稳定在一个较高的速率上。
这种蠕变现象对于结构设计和材料选择具有重要意义。
ANSYS可以模拟和分析四蠕变现象。
在进行四蠕变分析时,需要输入材料的蠕变参数。
这些参数可以通过试验或经验获得。
通过分析材料的力学性质和加载条件,可以获得材料的蠕变行为,并预测材料在长时间形变下的稳定性。
ANSYS高级非线性分析四蠕变的应用ANSYS高级非线性分析四蠕变在工程设计中具有广泛的应用。
以下列举了几个常见的应用领域:1. 高温结构设计在高温条件下,材料的蠕变行为是一个重要的考虑因素。
使用ANSYS高级非线性分析四蠕变,可以模拟和预测高温结构的蠕变行为,从而在结构设计中选择合适的材料和几何形状。
2. 轴承和齿轮设计轴承和齿轮是机械系统中常见的零件,需要承受高强度和高载荷。
使用ANSYS高级非线性分析四蠕变,可以评估轴承和齿轮材料的蠕变行为,从而选择合适的材料和结构参数。
ansys非线性瞬态结构分析重要命令
非线性各向同性硬化模型(TB,NLISO)选项基于Voce硬化规律或power硬化规律。该模型的优势在于材料行为由函数确定,而函数由TBDATA命令定义的四个材料常数确定。你可以通过拟合材料拉伸应力-应变曲线来得到这四个常数。不同于MISO,不需要担心如何恰当选定应力-应变点来输入。但是该选项只是适用于如下图所示的拉伸曲线。该选项适合大应变分析。可以综合Chaboche, creep, viscoplastic, and Hill anisotropy等选项来反正复杂材料行为。
Swelling Material Model
User-Defined Material Model
2 Plasticity塑性
大多数工程材料在达到所谓的弹性比例极限前都表现出线性的应力-应变关系。超出该极限,应力-应变关系变为非线性,但也不会变成完全没有弹性。塑性以不可恢复的变形为特点,当应力超过屈服极限材料即表现塑性。一般弹性极限与屈服极限差别很小,ANSYS中一般将这两点当成一点。塑性是一个不可恢复、与路径相关的现象。换句话说,载荷施加顺序及塑性响应顺序都影响最终结果。如果分析中会产生塑性形变,最好以较小的载荷步和时间步求解,以便模型会更遵循载荷-响应曲线。
双线性随动硬化模型(TB,BKIN)假设总应力范围等于屈服强度的两倍,以便包括包辛格效应。建议该选项使用于遵循von Mises屈服准则的一般小形变情况。不建议做大变形应用。BKIN选项可以综合蠕变和希尔各向异性选项来仿真更复杂的材料行为。
Multilinear Kinematic Hardening Material Model多线性随动硬化模型
ANSYS结构非线性分析指南(一至三章)
ANSYS结构⾮线性分析指南(⼀⾄三章)ANSYS结构⾮线性分析指南(⼀到三章)屈服准则概念:1.理想弹性材料物体发⽣弹性变形时,应⼒与应变完全成线性关系,并可假定它从弹性变形过渡到塑性变形是突然的。
2.理想塑性材料(⼜称全塑性材料)材料发⽣塑性变形时不产⽣硬化的材料,这种材料在进⼊塑性状态之后,应⼒不再增加,也即在中性载荷时即可连续产⽣塑性变形。
3.弹塑性材料在研究材料塑性变形时,需要考虑塑性变形之前的弹性变形的材料这⾥可分两种情况:Ⅰ.理想弹塑性材料在塑性变形时,需要考虑塑性变形之前的弹性变形,⽽不考虑硬化的材料,也即材料进⼊塑性状态后,应⼒不再增加可连续产⽣塑性变形。
Ⅱ.弹塑性硬化材料在塑性变形时,既要考虑塑性变形之前的弹性变形,⼜要考虑加⼯硬化的材料,这种材料在进⼊塑性状态后,如应⼒保持不变,则不能进⼀步变形。
只有在应⼒不断增加,也即在加载条件下才能连续产⽣塑性变形。
4.刚塑性材料在研究塑性变形时不考虑塑性变形之前的弹性变形。
这⼜可分两种情况:Ⅰ.理想刚塑性材料在研究塑性变形时,既不考虑弹性变形,⼜不考虑变形过程中的加⼯硬化的材料。
Ⅱ.刚塑性硬化材料在研究塑性变形时,不考虑塑性变形之前的弹性变形,但需要考虑变形过程中的加⼯硬化材料。
屈服准则的条件:1.受⼒物体内质点处于单向应⼒状态时,只要单向应⼒⼤到材料的屈服点时,则该质点开始由弹性状态进⼊塑性状态,即处于屈服。
2.受⼒物体内质点处于多向应⼒状态时,必须同时考虑所有的应⼒分量。
在⼀定的变形条件(变形温度、变形速度等)下,只有当各应⼒分量之间符合⼀定关系时,质点才开始进⼊塑性状态,这种关系称为屈服准则,也称塑性条件。
它是描述受⼒物体中不同应⼒状态下的质点进⼊塑性状态并使塑性变形继续进⾏所必须遵守的⼒学条件,这种⼒学条件⼀般可表⽰为f(σi j)=C⼜称为屈服函数,式中C是与材料性质有关⽽与应⼒状态⽆关的常数,可通过试验求得。
屈服准则是求解塑性成形问题必要的补充⽅程。
非线性有限元ANSYS(liuheng)
弹塑性力学
对象;任意变形体 特征;变形(屈服、
非线性) 任意形状的体 变量: 1)材料弹塑物性描述 2)变形方面描述 3)力的平衡描述 方程(微体dxdydz): 1)物理本构方程 (屈服、非线性) 2)几何变形方程 3)力的平衡方程 三大变量—三大方程 求解:解析半解析法
非变形体(刚体)
变形体
非线性有限元ANSYS讲义——刘恒
变形体及其受力情况的描述
基本变量:
u
==》 ui ij ij
(位移) (应变) (应力)
(如研究xyz三个方向,对应张量描述)
基本方程: 1)材料方面
三大类变量
2)几何方面
3)力平衡方面 求解方法:
三大类方程
1)经典解析方法
2)半解析方法
3)传统数值求解方法
难易程度
很难
较难
简单
求解精度
高
较高
低
方程最后形式
常微分方程
差分方程 积分方程->线性方程组 线性方程组
方法规范性 不规范,技巧要求高 较规范 只要试函数确定,后续过程很规范
方法通用性
不好
较好
较好
很好
解题范围 简单问题(非常有限) 较复杂问题
较大
大
其它
具一定物理背景(最小残差) 最小势能
第二章
非线性有限元概述
Fa - Fnr = 0
收敛是平衡的度量。
非线性有限元ANSYS讲义——刘恒
收 敛(续)
Newton-Raphson 迭代过程如下所示。基于 u0 时的结 构构形,计算出切向刚度KT,基于F 计算出的位移增量 是u ,结构构形更新为 u1。
载荷
ANSYS结构非线性分析指南_第四章
第四章材料非线性分析4.1 材料非线性概述许多与材料有关的参数可以使结构刚度在分析期间改变。
塑性、非线性弹性、超弹性材料、混凝土材料的非线性应力—应变关系,可以使结构刚度在不同载荷水平下(以及在不同温度下)改变。
蠕变、粘塑性和粘弹性可以引起与时间、率、温度和应力相关的非线性。
膨胀可以引起作为温度、时间、中子流水平(或其他类似量)函数的应变。
ANSYS程序应可以考虑多种材料非线性特性:1.率不相关塑性指材料中产生的不可恢复的即时应变。
2.率相关塑性也可称之为粘塑性,材料的塑性应变大小将是加载速度与时间的函数。
3.材料的蠕变行为也是率相关的,产生随时间变化的不可恢复应变,但蠕变的时间尺度要比率相关塑性大的多。
4.非线性弹性允许材料的非线性应力应变关系,但应变是可以恢复的。
5.超弹性材料应力应变关系由一个应变能密度势函数定义,用于模拟橡胶、泡沫类材料,变形是可以恢复的。
6.粘弹性是一种率相关的材料特性,这种材料应变中包含了弹性应变和粘性应变。
7.混凝土材料具有模拟断裂和压碎的能力。
8.膨胀是指材料在中子流作用下的体积扩大效应。
4.2 塑性分析4.2.1 塑性理论简介许多常用的工程材料,在应力水平低于比例极限时,应力—应变关系为线性的。
超过这一极限后,应力—应变关系变成非线性,但却不一定是非弹性的。
以不可恢复的应变为特征的塑性,则在应力超过屈服点后开始出现。
由于屈服极限与比例极限相差很小,ANSYS程序在塑性分析中,假设这二个点相同,见图4-1。
图4-1 弹塑性应力-应变曲线塑性是一种非保守的(不可逆的),与路径相关的现象。
换句话说,荷载施加的顺序,以及什么时候发生塑性响应,影响最终求解结果。
如果用户预计在分析中会出现塑性响应,则应把荷载处理成一系列的小增量荷载步或时间步,以使模型尽可能附合荷载—响应路径。
最大塑性应变是在输出(Jobname.OUT)文件的子步信息中打印的。
在一个子步中,如果执行了大量的平衡迭代,或得到大于15%的塑性应变增量,则塑性将激活自动时间步选项[AUTOTS ](GUI :Main Menu>Solution> Sol'n Control:Basic Tab 或 MainMenu>Solution>Unabridged Menu> Time /Frequenc>Time and Substps)。
ansys蠕变应变率命令流
ansys蠕变应变率命令流1.引言1.1 概述概述:蠕变应变是材料在长时间作用下产生的变形现象,它在工程领域具有重要的作用。
蠕变应变率是描述蠕变现象发生速率的参数,它反映了材料在一定应力下产生蠕变变形的能力。
ANSYS是一种常用的工程仿真软件,它提供了许多功能强大的命令流供用户使用。
其中,蠕变应变率命令流被广泛应用于材料蠕变性能的研究和工程实践中。
本文将对ANSYS蠕变应变率命令流进行详细介绍和使用方法的讲解。
通过学习和掌握这一命令流,工程师和研究人员可以更加准确地预测材料在长时间持续加载下的变形情况,为工程设计提供科学依据。
在介绍ANSYS蠕变应变率命令流之前,我们将先简要概述蠕变现象、蠕变应变和蠕变应变率的基本概念。
然后,我们将详细讲解蠕变应变率命令流的使用方法,包括命令的输入格式、参数的设置以及结果的分析等内容。
最后,我们将通过实例的展示和分析,进一步说明蠕变应变率命令流在实际工程中的应用价值。
通过本文的学习,读者将能够全面了解ANSYS蠕变应变率命令流的原理和功能,掌握蠕变应变率命令流的使用方法,并能够灵活应用于工程设计和材料研究中。
同时,本文也将对该命令流存在的一些局限性进行探讨,并对未来的研究方向进行展望。
1.2 文章结构本文主要介绍了ANSYS蠕变应变率命令流的使用方法。
文章共分为引言、正文和结论三个部分。
引言部分首先概述了本文的主题和目的。
其次,介绍了ANSYS蠕变应变率命令流在工程领域的重要性和应用背景。
最后,明确了本文的目的和意义,即通过详细介绍ANSYS蠕变应变率命令流的使用方法,为读者提供一个全面了解和掌握该技术的指南。
正文部分主要包括两个小节,分别是ANSYS蠕变应变率命令流的介绍和使用方法。
- 在2.1小节中,将详细介绍ANSYS蠕变应变率命令流的背景和原理。
首先,解释了蠕变现象的定义和原因。
接着,介绍了应变率是描述蠕变行为的重要参数。
随后,详细阐述了ANSYS蠕变应变率命令流在工程仿真中的作用和应用范围。
ansys非线性瞬态结构分析重要命令
3.1 瞬态分析的三种方法
瞬态分析有三种方法:full, mode-superposition ,andreduced。对于涉及非线性(plasticity, large deflections, large strain, and so on)的情况一般使用全积分方法。全积分也是最费时的方法。
自动时间步
AUTOTS,ON
DELTIM, DTIME, DTMIN, DTMAX, Carry
如果使用自动时间步,当Carry=OFF,以DTIME为起始时间步长,最小时间步不小于DTMIN,最大不大于DTMAX,当Carry=ON,以上一载荷步的最后子时间步长为起始时间步长。
另一个效果相同的命令组合:
OUTRES,Item,Freq,Cname
Item:NSOL,节点结果;ESOL,单元结果;ALL,所有。
Freq:n,每第n个子步;-n,均分成n段;NONE,一个也不存;ALL,每一子步;LAST,最后一子步;%array%按数组提供的时刻来存储。
例1:
NSUBST,6
OUTRES,ERASE设置到默认值,对于静态和瞬态分析,默认的是输出每一载荷步的最后子步的所有结果;谐态分析是每一子步。
线搜索选项
LNSRCH, Key
线搜索选项(LNSRCH)。该选项可代替自适应下降选项。如果线搜索选项是打开的,程序将自动关闭自适应选项。
非线性分析收敛标准
CNVTOL, Lab, VALUE, TOLER, NORM, MINREF
设置分析终结标准
NCNV, KSTOP, DLIM, ITLIM, ETLIM, CPLIM
Lab:
SPARSE-Sparse direct equation solver
ansys高级非线性分析四蠕变
第一阶段 第二阶段
断裂 第三阶段
t
September 30, 2001 Inventory #001491 4-9
Advanced Structural Nonlinearities 6.0
隐式和显式蠕变
... 术语的定义
Training Manual
• 蠕变的三个阶段(续)
– 蠕变应变率可能是应力、应变、温度、和/或时间的函数。
• 本章将介绍ANSYS中可用的大量的隐式和显式蠕变法则。
– 主要讨论金属的蠕变。 然而, 讨论的各种观点也适用于塑料或混凝土 等其它材料的蠕变。
– ANSYS有隐式和显式两种蠕变分析过程。 首先讨论蠕变的一般知识, 然后是进行隐式或显式蠕变分析的细节。
September 30, 2001 Inventory #001491 4-2
September 30, 2001 Inventory #001491 4-17
Advanced Structural Nonlinearities 6.0
隐式和显式蠕变
... 一般蠕变方程
Training Manual
• ANSYS 中可用的蠕变法则汇总如下:
Creep Equation Description Strain Hardening Time Hardening Generalized Exponential Generalized Graham Generalized Blackburn Modified Time Hardening Modified Strain Hardening Generalized Garofalo (Hyperbolic sine) Exponential Form Norton Time Hardening Rational Polynomial Generalized Time Hardening User Creep Annealed 304 Stainless Steel Annealed 316 Stainless Steel Annealed 2.25 Cr - 1 Mo Low Alloy Steel Power Function Creep Law Sterling Power Function Creep Law Annealed 316 Stainless Steel 20% Cold Worked 316 SS (Irradiation-Induced)
基于ANSYS的蠕变屈曲工程评定方法与实例
图1 蠕变屈曲变形-时间的特征基于ANSYS的蠕变屈曲工程评定方法与实例唐艳芳(上海时鉴工程技术有限公司, 上海 201203)[摘 要] 随着石油化工和核电工业的快速发展,越来越多的工业装置在高温下运行。
当设备或部件因压力或重力等载荷产生压缩应力时,就有可能发生蠕变屈曲。
本文介绍了ASME III-NH规范中关于蠕变屈曲的相关概念、评定方法、评定原理,同时用一个工程实例说明基于通用有限元软件进行蠕变屈曲的校核过程和实施步骤,为工程中防止蠕变屈曲失效提供了参考方案。
[关键词] 压力容器;高温;蠕变屈曲;ANSYS;实例作者简介:唐艳芳(1985—),女,江西南昌人,硕士研究生,中级工程师,从事压力容器分析及设计工作。
压力容器广泛应用于石化化工、核电等行业,如外压容器、夹套容器的内层、大型储罐球罐、核电站球形安全壳等结构,大多为薄膜结构,在承受压缩载荷的情况下极易发生屈曲失效。
随着压力容器的大型化、轻量化以及操作条件的复杂化,高温下的蠕变屈曲问题日益凸显。
本文基于ASME III-NH 分卷[1],阐述了蠕变屈曲的相关概念和原理、蠕变屈曲的类型、评定方法及评定原理,结合NB 分卷[2]介绍了与时间无关屈曲和与时间有关屈曲两种设计方法的不同之处。
同时,以某工程项目中的支承圆筒为例,展示了基于通用有限元软件ANSYS Workbench 实施蠕变屈曲计算、分析与评定的过程,为工程中蠕变屈曲的校核提供了较全面的参考。
1 屈曲的类型结构的屈曲或失稳一般需要考虑两种类型,一种是弹性或弹塑性屈曲,另一种是蠕变屈曲。
弹性或弹塑性屈曲可能在寿命期内的任意时间瞬间发生,只取决于几何结构和在时间上的短时材料响应,一般与时间无关。
蠕变屈曲可能是因为初始缺陷随时间得到增强从而导致几何失去稳定,这种失效与时间有关,甚至当载荷恒定时也有可能发生。
由于蠕变应变随时间进行累积,发生蠕变屈曲时的载荷要小于发生弹性和弹塑性屈曲的载荷。
ansys的非线性命令解析
引用小健哥的ANSYS 非线性分析命令解析ANSYS应用基于问题物理特性的自动求解控制方法,把各种非线性分析控制参数设置到合适的值。
如果用户对这些设置不满意,还可以手工设置。
下列命令的缺省设置已进行了优化处理:AUTOTS PRED MONITORDELTIM NROPT NEQITNSUBST TINTP SSTIFCNVTOL CUTCONTROL KBCLNSRCH OPNCONTROL EQSLVARCLEN CDWRITE LSWRITE这些命令及其设置在将在后面讨论。
参见《ANSYS Commands Reference》。
如果用户选择自己的设置而不是ANSYS的缺省设置,或希望用以前版本的ANSYS的输入列表,则可用/ SOLU 模块的SOLCONTROL ,OFF命令,或在/ BATCH 命令后用/ CONFIG ,NLCONTROL,OFF命令。
参见SOLCONTROL 命令的详细描述。
ANSYS对下面的分析激活自动求解控制单场的非线性或瞬态结构以及固体力学分析,在求解自由度为UX、UY、UZ、ROTX、ROTY、ROTZ 的结合时;单场的非线性或瞬态热分析,在求解自由度为TEMP时;注意--本章后面讨论的求解控制对话框,不能对热分析做设置。
用户必须应用标准的ANSYS 求解命令或GUI来设置。
2.2 非线性静态分析步骤尽管非线性分析比线性分析变得更加复杂,但处理基本相同。
只是在非线形分析的过程中,添加了需要的非线形特性。
非线性静态分析是静态分析的一种特殊形式。
如同任何静态分析,处理流程主要由以下主要步骤组成:建模;设置求解控制;设置附加求解控制;加载;求解;考察结果。
2.2.1 建模这一步对线性和非线性分析基本上是一样的,尽管非线性分析在这一步中可能包括特殊的单元或非线性材料性质,参考§4《材料非线性分析》,和§6.1《单元非线性》。
如果模型中包含大应变效应,应力─应变数据必须依据真实应力和真实(或对数)应变表示。
蠕变分析实例
图2
坐标轴设置对话框
ห้องสมุดไป่ตู้3
时间为 1000 小时的轴向应力结果显示
ANSYS 显示窗口将显示螺栓的有限元图,如图 1。
图1 五、施加载荷
螺栓有限元模型图
1) 施 加位 移 约束 。选 择 Preprocessor → Loads → Define Loads → Apply → Structural→Displacement→On Nodes 命令, 出现 Apply U, ROT on Nodes 拾取菜单,单击 Pick All 和 OK,在 Lab2 DOFs to be contrained 复选框 中选择 All DOF, 取 VALUE Displacement value 为 0。 2) 施加温度。选择 Preprocessor→Loads→Define Loads→Apply→Structural→Temperature→Uniform Temp 命令,取 Uniform Temperature 为 900。 六、求解计算 1) 定义分析类型。分析类型为 Static。 2)选择 Solution→Load Step Opts→Solution Ctrl 命令,出现 Nonlinear Solution Control 对话框,使 Solution Control 状态从 ON 变为 OFF,在 Pressure load stiffness 下拉菜单中选择 Program Chosen。 3) 定义求解时间步。选择 Solution→Load Step Opts→Time/Frequenc→Time and Substps 命令,取 Time at end of load step 为 3600000,取 Number of substeps 为 100,并选中 Stepped, 其余采用默认设置。 4) 求解输出控制。选择 Solution → Load Step Opts → Output Ctrls → Solu Printout 命令,在 Item for printout control 下拉菜单中选择 Basic quantities,并选中 Every Nth substp, 取 N 值为 36000,Component name 为 All entities。 5) 写入数据库和结果文件控制。选择 Solution → Load Step Opts → Output Ctrls→DB/Result Files 命令,在 Item to be controlled 下拉菜单中选 择 Element solution, 并选中 Every Nth substp, 取 N 值为 1,Component name 为 All entities。 6) 选择 Solution→Solve→Current LS 命令,单击 OK, ANSYS 将开始求解计 算,求解结束时,出现 Note 对话框。 七、查看求解结果 1) 轴向应力和时间变化关系曲线显示 ·定义时间-历程变量。选择 TimeHist Postpro→Define Variables 命令,单 击 Add 按钮,选中 by seq no.,单击 OK 按钮,出现单元拾取菜单,在输入栏
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
September 30, 2001 Inventory #001491 4-6
隐式和显式蠕变
... 蠕变背景
Training Manual
Advanced Structural Nonlinearities 6.0
• 与塑性类似, 蠕变是一种基于偏量行为的不可逆(非弹性)应变。 在 蠕变流动条件下, 假设材料是不可压缩的。
September 30, 2001 Inventory #001491 4-17
隐式和显式蠕变
... 一般蠕变方程
• ANSYS 中可用的蠕变法则汇总如下:
Creep Equation Description Strain Hardening Time Hardening Generalized Exponential Generalized Graham Generalized Blackburn Modified Time Hardening Modified Strain Hardening Generalized Garofalo (Hyperbolic sine) Exponential Form Norton Time Hardening Rational Polynomial Generalized Time Hardening User Creep Annealed 304 Stainless Steel Annealed 316 Stainless Steel Annealed 2.25 Cr - 1 Mo Low Alloy Steel Power Function Creep Law Sterling Power Function Creep Law Annealed 316 Stainless Steel 20% Cold Worked 316 SS (Irradiation-Induced) Type Primary Primary Primary Primary Primary Primary Primary Secondary Secondary Secondary Both Both Primary Both Both Both Primary Both Both Both Explicit Implicit C6/C12 value TBOPT value C6=0 1 C6=1 2 C6=2 3 4 5 6 7 8 C12=0 9 C12=1 10 11 C6=15 12 13 C6=100 100 C6=9 C6=10 C6=11 C6=12 C6=13 C6=14 C66=5 -
Training Manual
Advanced Structural Nonlinearities 6.0
– 隐式蠕变应用了Euler向后积分法求解蠕变应变。 该方法在数值上无条 件稳定,这意味着不必象显式蠕变方法那样, 使用小的时间步, 所以总 体上会更快。
cr f t Dt , t Dt , T t Dt ,
D : el
el pl cr
计算塑性应变(流动法则)类似于前一章的描述。 根据蠕变应变率方 程计算蠕变应变,它的详细形式将在后面讨论。
• 弹性、蠕变和塑性应变都基于(当前的)应力状态,但它们是独立计 算的(彼此互不依赖)。
– 注意在利用隐式蠕变与显式蠕变进行计算时有差别。
cr e
应力相关性
Q RT
式中Q为激活能, R为普适气体常数, T为绝对温度。
– 蠕变应变通常也与应力有关, 尤其是位错蠕变。Norton 法则为:
cr n
对上述幂定律的常见修正如下:
cr eC
September 30, 2001 Inventory #001491 4-16
Training Manual
Advanced Structural Nonlinearities 6.0
– 第一阶段, 应变率随时间减小,该阶段在较短时间内完成。第二阶段 具有常应变率。第三阶段, 应变率迅速增加直到破坏(断裂)。
断裂
第一阶段
第二阶段
第三阶段
t
September 30, 2001 Inventory #001491 4-9
• 本章将介绍ANSYS中可用的大量的隐式和显式蠕变法则。
– 主要讨论金属的蠕变。 然而, 讨论的各种观点也适用于塑料或混凝土 等其它材料的蠕变。 – ANSYS有隐式和显式两种蠕变分析过程。 首先讨论蠕变的一般知识, 然后是进行隐式或显式蠕变分析的细节。
September 30, 2001 Inventory #001491 4-2
– 第一阶段蠕变的应变率通常远大于第二阶段蠕变。然而, 应变率在第 一阶段逐渐降低而在第二阶段几乎为常值 (对于前面提及的常应力、恒 温下单轴试验情况而言), 而且第一阶段蠕变时间比第二阶段短。
September 30, 2001 Inventory #001491 4-10
隐式和显式蠕变
... 术语的定义
– 对于隐式蠕变加上率无关塑性, 塑性修正和蠕变修正同时进行, 而不是 分别进行。 因此, 隐式蠕变一般比显式蠕变更精确, 但它仍与时间步大 小有关,必须使用足够小的时间步来精确捕捉路径相关行为。
• 基于上述原因, 隐式蠕变是 ANSYS 推荐的方法(高效、精确)。 后 面将详细讨论两种蠕变过程。
September 30, 2001 Inventory #001491 4-14
Bret Zahn, ChipPAC 授权的ANSYS 模型 (/)
September 30, 2001 Inventory #001491 4-8
隐式和显式蠕变
B. 术语的定义
• 蠕变的三个阶段:
– 在常载荷条件下, 蠕变的单轴应变与时间的关系如下图所示。
... 术语的定义
• 时间强化
– 假设蠕变应变率仅与从蠕变过程开始 的时间有关, 也就是说, 该曲线 上/下移动,当应力从1 到 2变化时, 计算A 到 B点的不同蠕变速率。
Training Manual
Advanced Structural Nonlinearities 6.0
1
A
2
• 应变强化
cr f1 f 2 f 3 t f 4 T
• 函数f1-4 与选择的蠕变法则有关。
– 通常由不同应变速率和温度条件下的各种拉伸实验确定相关蠕变常数。
• 假设各向同性行为, von Mises 方程用于计算有效应力, 在蠕变应 变率方程中使用等效应变(与率无关塑性相似)。
• 另一方面, 与率无关塑性不同,蠕变在发生非弹性应变时没有屈服 面。
– 因此, 蠕变不需要高应力值来发生更多的蠕变应变。 假设在所有非零 应力值时都会发生蠕变应变。
• 前面提到, 从材料的角度看, 蠕变和粘塑性是相同的。
– 工程应用中, 通常蠕变用于描述低应变速率的热激活过程。 率无关塑 性和隐式蠕变应变以弱耦合方式处理。 – 相反, ANSYS 中粘塑性本构模型用于描述高应变速率的应用(例如, 冲 击载荷)。 非弹性应变以强耦合方式处理。
隐式和显式蠕变
... 术语的定义
• 蠕变的三个阶段(续)
– 蠕变应变率可能是应力、应变、温度、和/或时间的函数。
Training Manual
Advanced Structural Nonlinearities 6.0
– 对于工程分析, 通常最关心蠕变的第一和第二阶段, 第三阶段蠕变通 常与开始破坏有关(颈缩、损坏), 且时间较短, 因此, ANSYS 中不模 拟第三阶段。
September 30, 2001 Inventory #001491 4-7
隐式和显式蠕变
... 蠕变分析实例
焊球蠕变分析的例子(热循环)。
Training Manual
Advanced Structural Nonlinearities 6.0
单元185 (B-Bar), 双曲 正弦隐式蠕变模型
四章
蠕变g Manual
Advanced Structural Nonlinearities 6.0
• 前一章探讨了率无关塑性, 本章将讨论 ANSYS 中蠕变行为的分析 方法。
– 尽管从材料的观点看蠕变 和粘塑性是相同的, 但本构模型的使用不同 。 因此, 率相关的塑性主题被分为两部分, 这部分属于蠕变。
September 30, 2001 Inventory #001491 4-3
隐式和显式蠕变
A. 蠕变背景
Training Manual
Advanced Structural Nonlinearities 6.0
• 晶体材料中, 如金属, 蠕变机理与空隙的扩散流动和位错运动有关。
– 空隙是点缺陷, 倾向于形成与所施加应力方向垂直(而不是平行)的晶界 。空隙由高集中区向低集中区运动。 在低应力状态下发生扩散流动, 但通常需要高温条件。 – 晶粒的位错是线缺陷. 位错运动(攀升、滑动、偏移)在高应力状态下被 激活, 尽管在中温时也可能发生位错运动。 – 有时晶界滑移被认为是一种独立的导致蠕变变形的机理。
隐式和显式蠕变
... 一般蠕变方程
通常第一阶段蠕变显示时间或应变强化。
– 时间强化包含在一个时间相关项中:
Training Manual
Advanced Structural Nonlinearities 6.0
cr t m
– 应变强化包含在一个应变相关项中:
cr n
– 由可用的材料数据来决定使用哪一项(应变或时间强化)。 – 第二阶段不具有时间或应变强化,第二阶段的蠕变应变率通常是常数。
September 30, 2001 Inventory #001491 4-4
隐式和显式蠕变