ansys非线性分析指南.pdf

《ansys非线性分析指南.pdf》由会员分享，可在线阅读，更多相关《ansys非线性分析指南.pdf（11页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、ANSYS非线形分析指南 弹塑性分析 弹塑性分析 弹塑性分析 在这一册中,我们将详细地介绍由于塑性变性引起的非线性问题-弹塑性分析,我们的介绍人为以下几个方面:什么是塑性 塑性理论简介 ANSYS 程序中所用的性选项 怎样使用塑性 塑性分析练习题 什么是塑性什么是塑性 塑性是一种在某种给定载荷下,材料产生永久变形的材料特性,对大多的工程材料来说,当其应力低于比例极限时,应力一应变关系是线性的 另外,大多数材料在其应力低于屈服点时表现为弹性行为也 就 是说当 移 走 载 荷 时其应变也完全消失 由于屈服点和比例极限相差很小因此在 ANSYS 程序中假定它们相同在应力一应变的曲线中 低于屈服点的叫

2、作弹性部分 超过屈服点的叫作塑性部分 也叫作应变强化部分塑性分析中考虑了塑性区域的材料特性 路径相关性路径相关性 即然塑性是不可恢复的 那么这种问题的就与加载历史有关 这类非线性问题叫作与路径相关的或非保守的非线性 路径相关性是指对一种给定的边界条件可能有多个正确的解内部的应力应变分布存在为了得到真正正确的结果我们必须按照系统真正经历的加载过程加载 率相关性率相关性 塑性应变的大小可能是加载速度快慢的函数 如果塑性应变的大小与时间有关 这种塑性叫作率无关性塑性相反与应变率有关的性叫作率相关的塑性 大多的材料都有某种程度上的率相关性 但在大多数静 力分 析所经历的应变率范围两者的应力应变曲线差别

3、不大所以在一般的分析中我们变为是与率无关的 工程应力工程应力应变与真实的应力应变与真实的应力应变应变 塑性材料的数据一般以拉伸的应力应变曲线形式给出材料数据可能是工程应力PA0与工程应变ll0也可能是真实应力P/A与真实应变nLll(0)大应变的塑性分析一般采用真实的应力 应变数据而小应变分析一般采用工程的应力应变数据 什么时候激活塑性什么时候激活塑性 当材料中的应力超过屈服点时塑性被激活也就是说有塑性应变发生而屈服应力本身可能是下列某个参数的函数 温度 应变率 以前的应变历史 侧限压力 其它参数 塑性理论介绍 塑性理论介绍 在这一章中我们将依次介绍塑性的三个主要方面 屈服准则 流动准则 强化

4、准则 屈服准则屈服准则 对单向受拉试件 我们可以通过简单的比较轴向应力与材料的屈服应力来决定是否有塑性变形发生然而对于一般的应力状态是否到达屈服点并不是明显的 屈服准则是一个可以用来与单轴测试的屈服应力相比较的应力状态的标量表示 因此知道了应力状态和屈服准则程序就能确定是否有塑性应变产生 第1页 ANSYS非线形分析指南 弹塑性分析 屈服准则的值有时候也叫作等效应力一个通用的屈服准则是 Von Mises 屈服准则当等效应力超过材料的屈服应力时将会发生塑性变形 可以在主应力空间中画出 Mises 屈服准则见 图 31 在 3D 中屈服面是一个以为轴的圆柱面123=在 2D 中屈服面是一个椭圆在

5、屈服面内部的任何应力状态都是弹性的屈服面外部的任何应力状态都会引起屈服注意静水压应力状态123=不会导致屈服屈服与静水压应力无关而只与偏差应力有关因此1180=230=123180=的应力状态比的应力状态接近屈 服Mises 屈服准则是一种除了土壤和脆性材料外典型使用的屈服准则在土壤和脆性材料中屈服应力是与静水压应力侧限压力有关的侧限压力越高发生屈服所需要的剪应力越大 流动准则流动准则 流动准则描述了发生屈服时塑性应变的方向也就是说流动准则定义了单个塑性应变分量xplypl 等随着屈服是怎样发展的 一般来说 流动方程是塑性应变在垂直于屈服面的方向发展的屈服准则中推导出来的这种流动准则叫作相关流

6、动准则如果不用其它的流动准 则从其它不同的函数推导出来则叫作不相关的流动准则 强化准则强化准则 强化准则描述了初始屈服准则随着塑性应变的增加是怎样发展的 一般来说屈服面的变化是以前应变历史的函数在 ANSYS 程序中使用了两种强化准则 等向强化等向强化是指屈服面以材料中所作塑性功的大小为基础在尺寸上扩张对 Mises 屈服准则来说屈服面在所有方向均匀扩张见图 3-2 第2页 ANSYS非线形分析指南 弹塑性分析 图 3-2 等向强化时的屈服面变化图 由于等向强化在受压方向的屈服应力等于受拉过程中所达到的最高应力 随动强化随动强化假定屈服面的大小保持不变而仅在屈服的方向上移 动当某个方向的屈服应

7、力升高时其相反方向的屈服应力应该降低见图 3-3 图 3-3 随动强化时的屈服面变化图 在随动强化中 由于拉伸方向屈服应力的增加导致压缩方向屈服应力的降低 所以在对应的两个屈服应力之间总存一个的差值2y初始各向同性的材料在屈服后将不再是向同性的 塑性选项 塑性选项 ANSYS 程序提供了多种塑性材料选项 在此主要介绍四种典型的材料选项可以通过激活一个数据表来选择这些选项 经典双线性随动强化 BKIN 双线性等向强化 BISO 多线性随动强化 MKIN 多线性等向强化 MISO 经典的双线性随动强化经典的双线性随动强化BKINBKIN使用一个双线性来表示应力应变曲线所以有两个斜率弹 第3页 AN

8、SYS非线形分析指南 弹塑性分析 性斜率和塑性斜率由于随动强化的 Vonmises 屈服准 则被使用所以包含有鲍辛格效应此选项适用于遵守 Von Mises 屈服准则初始为各向同性材料的小应变问题这包括大多数的金属 需要输入的常数是屈服应力和切向斜率yTE可以定义高达六条不同温度下的曲线 注意 使用 MP 命令来定义弹性模量 弹性模量也可以是与温度相关的 切向斜率 Et 不可以是负数也不能大于弹性模量 在使用经典的双线性随动强化时可以分下面三步来定义材料特性 1 定义弹性模量 2 激活双线性随动强化选项 3 使用数据表来定义非线性特性 双线性等向强化双线性等向强化BIS0BIS0也是使用双线性

9、来表示应力应变曲线在此选项中等向强化的 Von Mises 屈服准则被使用这个选项一般用于初始各向同性材料的大应变问题需要输入的常数与 BKIN 选项相同 多线性随动强化多线性随动强化MKINMKIN使用多线性来表示应力应变曲线模拟随动强化效应这个选项使用 Von Mises 屈服准则对使用双线性选项BKIN不 能足够表示应力应变曲线的小应变分析是有用的 需要的输入包括最多五个应力应变数据点 用 数 据 表 输 入可以定义五条不同温度下的曲线 在使用多线性随动强化时可以使用与 BKIN 相同的步骤来定义材料特性所不同的是在数据表中输入的常数不同下面是一个用命令流定义多线性随动强化的标准输入 M

10、PTEMP1070 MPDATAEX330ES25ES TBMK2N3 TBTEMPSTRA2N TBDATA0.010.050.1 TBTEMP10 TBDATA300003700038000 TBTEMP70 TBDATA2250003100033000 多线性等向强化多线性等向强化MISO使用多线性来表示使用 Von Mises 屈服 准则的等向强化的应力应变曲线它适用于比例加载的情况和大应变分析 需要输入最多 100 个应力应变曲线最多可以定义 20 条不同温度下的曲线 其材料特性的定义步骤如下 1 定义弹性模量 2 定义 MISO 数据表 3 为输入的应力应变数据指定温度值 4 输入

11、应力应变数据 5 画材料的应力应变曲线 与 MKIN 数据表不同的是MISO 的数据表对不同的温度可以有不同的应变值因此每条温度曲线有它自己的输入表 怎 样 使 用 塑 性 怎 样 使 用 塑 性 在这一章中我们将介绍在程序中怎样使用塑性重点介绍以下几个方面 可 用 的 ANSYS 输 入 ANSYS 输 出 量 使 用 塑 性 的 一 些 原 则 第4页 ANSYS非线形分析指南 弹塑性分析 加 强 收 敛 性 的 方 法 查 看 塑 性 分 析 的 结 果 ANSYS 输 入ANSYS 输 入 当使用 TB 命令选择塑性选项和输入所需常数时应该考虑到 常数应该是塑性选项所期望的形式 例如我

12、们总是需要应力和总的应变而不是应力与塑性应变 如果还在进行大应变分析应力应变曲线数据应该是真实应力真实应 变 对双线性选项BKINBISO输入常数和可以按下述方法来决定yTE如果材料没有明显的屈服应力y通常以产生 0.2%的塑性应变所对应的应力作为屈服应力而可以通过在分析中所预期的应变范围内来拟合实验曲线得到TE 其它有用的载荷步选项其它有用的载荷步选项 使用的子步数使用的时间步长既然塑性是一种与路径相关的非线性因此需要使用许多载荷增量来加载 激活自动时间步长 如果在分析所经历的应变范围内应力应变曲线是光滑的使用预测器选项这能够极大的降低塑性分析中的总体迭代数 输出量 输出量 在塑性分析中对每

13、个节点都可以输出下列量 EPPL塑性应变分量,等等 xplyplEPEQ累加的等效塑性应变 SEPL根据输入的应力应变曲线估算出的对于 EPEQ 的等效应 力 HPRES静水压应力 PSV塑性状态变量 PLWK单位体积内累加的塑性功 上面所列节点的塑性输出量实际上是离节点最近的那个积分点的值 如果一个单元的所有积分点都是弹性的 EPEQ0那么节点的弹性应变和应力从积分点外插得到如果任一积分点是塑性的EPEQ0那么节点的弹性应变和应力实际上是积分点的值这是程序的缺省情况但可 以人为的改变它 程序使用中的一些基本原则程序使用中的一些基本原则 下面的这些原则应该有助于可执行一个精确的塑性分析 1 所

14、需要的塑性材料常数必须能够足以描述所经历的应力或应变范围内的材料特性 2 缓慢加载应该保证在一个时间步内最大的塑性应变增量小于 5%一 般 来说如果 Fy 是系统刚开始屈服时的载荷那么在塑性范围内的载荷增量应近似为 0.05*Fy 对用面力或集中力加载的情况 Fy 对用位移加载的情况 3 当模拟类似梁或壳的几何体时必须有足够的网格密度为了能够足够的模拟弯曲反应在厚度方向必须至少有二个单元 4 除非那个区域的单元足够大应该避免应力奇异由于建模而导致的应力奇异有 单点加载或单点约束 凹角 模型之间采用单点连接 单点耦合或接触条件 5 如果模型的大部分区域都保持在弹性区内那么可以采用下列方法来降低计

15、算时间 在弹性区内仅仅使用线性材料特性 不 使 用 TB 命 令 在线性部分使用子结构 加强收敛性的方法加强收敛性的方法 如果不收敛是由于数值计算导致的可以采用下述方法来加强问题的收敛性 1使用小的时间步长 第5页 ANSYS非线形分析指南 弹塑性分析 2 如果自适应下降因子是关闭的打开它相反如果它是打开的 且割线刚度正在被连续地使用那么关闭它 3使用线性搜索特别是当大变形或大应变被激活时 4预测器选项有助于加速缓慢收敛的问题但也可能使其它的问题变得不稳定 5可以将缺省的牛顿拉普森选项转换成修正的MODI或初始刚度INIT牛顿拉普森选项这两个选项比全牛顿拉普森选项更稳定 需要更的迭代但这两个选

16、项仅在小挠度和小应变塑性分析中有效 查 看 结 果 查 看 结 果 1 感兴趣的输出项例如应力变形支反力等等对加载历史的响应应该是光滑的一个不光滑的曲线可能表明使用了太大的时间步长或太粗的网 格 2 每个时间步长内的塑性应变增量应该小于 5这个值在输出文件中以Max plastic Strain Step输出也可以使用 POST26 来显示这个值Main Menu:Time Hist Postpro Define Variables 3 塑性应变等值线应该是光滑的通过任一单元的梯度不应该太大 4 画出某点的应力应变图应力是指输出量 SEQVMises 等 效 应 力总应变由累加的塑性应变 EP

17、EQ 和弹性应变得来 塑性分析实例塑性分析实例GUI 方法GUI 方法 在这个实例分析中我们将进行一个圆盘在周期载荷作用下的塑性分析 问题描述问题描述 一个周边简支的圆盘在其中心受到一个冲杆的周期作用由于冲杆被假定是刚性的因此在建模时不考虑冲杆而将圆盘上和冲杆接触的结点的 Y 方向上的位移耦合起来 由于模型和载荷都是轴对称的 因此用轴对称模型来进行计算 求解通过四个载荷步实现 问题详细说明问题详细说明 材料性质 EX=70000 (杨氏模量 NUXY=0.325泊松比 塑性时的应力应变关系如下 应变 应力 0.0007857 55 0.00575 112 0.02925 172 0.1 241

18、 加载历史 时间 载荷 0 0 1 -6000 2 750 3 -6000 问题描述图问题描述图 第6页 ANSYS非线形分析指南 弹塑性分析 6.513010F 步骤一步骤一建立计算所需要的模型建立计算所需要的模型 在这一步中建立计算分析所需要的模型包括定义单元类型划分网格给定边界条件并将数据库文件保存为exercise2.db 在此对这一步的过程不作详细叙述 步骤二步骤二恢复数据库文件恢复数据库文件exercise.dbexercise.db Utility MenuFileResume from 步骤三步骤三定义材料性质 定义材料性质 1选 择 菜 单 路 径 Main MenuPrep

19、rocessorMatersal Props-Constant-Isotropic.Isotropic Matersal Properties(各向同性材料性质对话框出现 2单击 OK 来指定材料号为 1另一个 I sotropic Material Properties 对话框出现 3对杨氏模量EX键入 EXX 4对泊松比NUXY键入 0.325 5单击 OK 步骤四步骤四定义和填充多线性随动强化数据表定义和填充多线性随动强化数据表MKINMKIN 1选择菜单路径 Main MenuPreprocessorMaterial PropsDataTablesDefine/Activate.Def

20、ine/Activate Data Table(激活数据表对话框出现 2在关于 type of data table(数据表类型的卷动框中卷动到Multi kinem MKIN且选中它 3在 material refersuce number(材料参考号中健入 1 4对 number of temperatures(温度数键入 1单击 OK 5选 择 菜 单 路 径Main MenuPreprocessorMaterial PropsData TablesEdit Active.Data Table MKIN 对话框出现 6在Strain一行中从第二列起分别输入 STN1,STN2,STN3,

21、STN4 7在Curve 1一行中从第二列起分别输入 STS1,STS2,STS3,STS4 8选择 FileApply&Quit 9选择菜单路径 Main MenuPreprosessorMaterial PorpsData TablesGraph.Graph Data Tables(图形表示数据表对话框出现 10单击 OK 接受绘制 MKIN 表的缺省一个 MKIN 表的标绘图出现在 ANSYS 图形窗口中 步骤五步骤五进入求解器 进入求解器 选择菜单路径 Main MenuSolution 步骤六步骤六定义分析类型和选项 定义分析类型和选项 第7页 ANSYS非线形分析指南 弹塑性分析

22、1选择菜单路径 Main MenuSolution-Analysis Type-New Analysis.2单击Static来选中它然后单击 OK 步骤七步骤七打开预测器打开预测器设置输出控制设置输出控制 1选择菜单路径 Main menusolution-Load Set Opts-NonlinearPredictor 2将 predictor 的状态设置为ON 3选择菜单路径 Main MenuSolution-Load Step Options-Output Ctrls DB/Results File.Coutrols for Database and Results File Writ

23、ing(对数据库和结果文件写入的控制对话框出现 4单击Every substep且选中它 步骤八步骤八设置载荷步选项 设置载荷步选项 1选 择 菜 单 路 径Main MenuSolution-Load Step Options-Time/Frequenc time&Substep Time&Substep Option(时间和子步数选项对话框出现 2对 time at end of Load Step(载荷步终止时间键入 1e-6 3对 Number of substeps(子步数键入 1 步骤九步骤九对第一个载荷步加载 对第一个载荷步加载 在结点 3 的 Y 方向施加一大小为 0 的集中力

24、载荷 步骤十步骤十将第一个载荷步写入载荷步文件将第一个载荷步写入载荷步文件 1选择菜单路径 Main MenuSolution-Write Ls File出现对话框 2在LSNUM的输入框中键入 1 步骤十一步骤十一对第二个载荷步加载对第二个载荷步加载并写入载荷步文件并写入载荷步文件 1选择菜单路径Main MenuSolution-Load Step Options-Time/Frequenctime&Substep Time&Substep Option(时间和时间步选项对话框出现 2对 time at end of Load Step(载荷步终止时间键入 1 3对 Number of s

25、ubsteps(子步数键入 10 4单击 automatic time stepping option自动时间步长选项使之为 ON然后单击 OK 5在结点 3 的 Y 方向施加一大小为-6000 的集中力载荷 6选择菜单路径 Main MenuSolution-Write Ls File出现对话框 7在LSNUM的输入框中键入 2 步骤十二步骤十二对第三个载荷步加载对第三个载荷步加载并写入载荷步文件并写入载荷步文件 1选择菜单路径Main MenuSolution-Load Step Options-Time/Frequenctime&Substep Time&Substep Option(时

26、间和时间步选项对话框出现 2对 time at end of Load Step(载荷步终止时间键入 2 3在结点 3 的 Y 方向施加一大小为 750 的集中力载荷 4选择菜单路径 Main MenuSolution-Write Ls File出现对话框 5在LSNUM的输入框中键入 3 步骤十三步骤十三对第四个载荷步加载对第四个载荷步加载并写入载荷步文件并写入载荷步文件 1选择菜单路径 Main MenuSolution-Load Step Options-Time/Frequenctime&Substep Time&Substep Option(时间和 时间步选项对话框出现 2对 tim

27、e at end of Load Step(载荷步终止时间键入 3 3在结点 3 的 Y 方向施加一大小为-6000 的集中力载荷 4选择菜单路径 Main MenuSolution-Write Ls File出现对话框 5在LSNUM的输入框中键入 4 步骤十三步骤十三求解问题 求解问题 1选择菜单路径 Main MenuSolution-Solve-From Ls Files对话框出现 2对LSMIN键入 1对LSMAX键入 4 3单击对话框中的 OK 开始求解 第8页 ANSYS非线形分析指南 弹塑性分析 步骤十四步骤十四进行后处理进行后处理 在这一步中可以进行所想要的后处理在此不进行详

28、述 非线性静态实例分析非线性静态实例分析命令流方式命令流方式 你可以用下面显示的 ANSYS 命令替代 GUI 选择进行上面这个例题的塑性分析 fini/cle/title,circular plate loaded by a circular punch-kinematic hardening rpl=65 rpu=5 h=6.5 exx=70000 sts1=55 !yield stress stn1=sts1/exx sts2=112 stn2=0.00575 sts3=172 stn3=0.02925 sts4=241 stn4=0.1 nex=15 !No.of elements a

29、long the radius net=2 !No.of elements in the plates thickness /prep7 et,1,42,1!define node for convenient postprocessing n,1,rpl,0 n,2,0,0 n,3,0,h/2 !define geometry k,1,-(h/2)k,2,rpu,-(h/2)k,3,rpl,-(h/2)kgen,2,1,3,1,h,3 nex1=nint(0.8*net)nex2=nex-nex1 l,1,2 l,4,5 l,2,3 l,5,6 a,1,2,5,4 a,2,3,6,5 les

30、ize,1,nex1 lesize,2,nex1 lesize,3,nex2,2.5 lesize,4,nex2,2.5 esize,(h/net)amesh,all nummrg,node nodes 第9页 ANSYS非线形分析指南 弹塑性分析 nsel,s,loc,x,0,rpu nsel,r,loc,y,(h/2)cp,1,uy,all nsel,all esel,all fini /solu nsel,s,loc,x,rpl nsel,r,loc,y,0 d,all,uy nsel,s,loc,x,0 d,all,ux nsel,all outres,all,all fini /pr

31、ep7 mp,ex,1,exx mp,nuxy,1,0.325 tb,mkin,1 tbtemp,strain tbdata,stn1,stn2,stn3,stn4 tbtemp,tbdata,sts1,sts2,sts3,sts4 fini /solu pred,on outres,all,all nsubst,1 time,1e-6 f,3,fy,0 lswrite autots,on nsubst,10 time,1 f,3,fy,-6000 lswrite time,2 f,3,fy,750 lswrite time,3 f,3,fy,-6000 lswrite lssolve,1,4 fini /post1 set,2/dscal,1,1 pldisp,2 第10页 ANSYS非线形分析指南 弹塑性分析 fini /post26 nsol,2,2,u,y,uy2 rforce,3,1,f,y,ry1 add,2,2,uy2,-1/grid,1/axlab,x,deflection mm/axlab,y,force n xval,2 plvar,3 prvar,2,3 fini 第11页