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：MainMenuSolutionSol'nControl:BasicTab或MainMenuSolutionUnabridgedMenuTime/FrequencTimeandSubstps)。如果取了太大的时间步，则程序将二分时间步，并重新求解。其他类型的非线性行为可以与塑性同时产生。实际上，大位移和大应变几何非线性经常伴随有塑性材料响应，如果用户预计在结构中存在大变形，则必须在分析中用NLGEOM命令激活这些效应(GUI:MainMenuSolutionSol'nControl:BasicTaborMainMenuSolutionUnabridgedMenuAnalysis)。对于大应变分析，材料应力—应变特性必须按真实应力和对数应变输入。在这一节中，我们将依次介绍塑性理论的三个主要方面：屈服准则流动准则强化准则4.2.1.1屈服准则对单向受拉试件，我们可以通过简单的比较轴向应力与材料的屈服应力来决定是否有塑性变形发生，然而，对于一般的应力状态，是否到达屈服点并不是明显的。屈服准则是一个可以用来与单轴测试的屈服应力相比较的应力状态的标量表示。因此，知道了应力状态和屈服准则，程序就能确定是否有塑性应变产生。在多轴应力状态下，屈服准则可以用下式来表示：yef其中e为等效应力，y为屈服应力。当等效应力超过材料的屈服应力时，将会发生塑性变形。1．VonMises屈服准则VonMises屈服准则是一个比较通用的屈服准则，尤其适用于金属材料。对于VonMises屈服准则，其等效应力为：23123222121e其中1，2，3为三个主应力。可以在主应力空间中画出Mises屈服准则，见图4-2。图4-2主应力空间中的Mises屈服面在3D主应力空间中，Mises屈服面是一个以321为轴的圆柱面，在2D中，屈服面是一个椭圆，在屈服面内部的任何应力状态，都是弹性的，屈服面外部的任何应力状态都会引起屈服。注意：静水压应力状态（321）不会导致屈服：屈服与静水压应力无关，而只与偏差应力有关，因此，1801，032的应力状态比180321的应力状态接近屈服。Mises屈服准则是一种除了土壤和脆性材料外典型使用的屈服准则，在土壤和脆性材料中，屈服应力是与静水压应力（侧限压力）有关的，侧限压力越高，发生屈服所需要的剪应力越大。2．Hill屈服准则以上介绍的VonMises屈服准则是一个各向同性的屈服准则，而Hill屈服准则是各向异性的，可以考虑材料的弹性参数的各向异性和屈服强度的各向异性。它是VonMises屈服准则的延伸。Hill屈服准则的等效应力被表达为：222222222xzyzxyzxzyyxeMLNGFH其中六个材料常数H,F,G,N,L,M可以由试验结果取得，如下各式所示：21xxRHG，21yyRHF，21zzRGF223xyRN，223yzRL，223xzRM以上各式右端的六个常书为材料的屈服应力比率，如下：0yxxxxR，0yyyyyR，0yzzzzR03yxyxyR，03yyzyzR，03yxzxzR其中yij为各个方向的屈服应力，0为参考屈服应力。Hill屈服面在主应力空间中如图4-3。图4-3主应力空间中的Hill屈服面在3D主应力空间中，Hill屈服面是一个以321为轴的椭圆柱面，在2D中，屈服面也是一个椭圆，在屈服面内部的任何应力状态，都是弹性的，屈服面外部的任何应力状态都会引起屈服。3．广义Hill屈服准则广义Hill屈服准则是对Hill屈服准则的进一步延伸，广义希尔准则不仅考虑了在材料的三个正交方向屈服强度的不同，而且还考虑了拉伸状态和压缩状态下屈服强度的不同。广义Hill屈服准则的等效应力可以表示为：213131LMTTe其中：665544332313232212131211000000000000000000000000MMMMMMMMMMMMM，TLLLL0,0,0,,,321jjjjKM,j=1to6)(2133221112MMMM)(2133221113MMMM)(2133221123MMMM)(jjjjjMLj=1to3以上各式中，j和j分别是j方向的拉伸和压缩屈服强度，在此，压缩屈服应力被作为正值处理。对剪切屈服，jj。令111M，则有：xxK单轴屈服强度值的选取需使得以下两式成立0zzzzyyyyxxxx0)(2331133222211233222211MMMMMMMMM广义Hill屈服准则使用等向强化准则，因此在应力空间中，屈服面是一个经过移动的椭圆柱面，其大小随塑性应变而胀缩，见图4-4。图4-4主应力空间中的广义Hill屈服面4．Drucker-Prager(DP)屈服准则Drucker-prager屈服准则是对Mohr-Coulomb准则的近似，它修正了VonMises屈服准则，即在VonMises表达式中包含一个附加项。其屈服面并不随着材料的逐渐屈服而改变，因此没有强化准则,塑性行为被假定为理想弹塑性，然而其屈服强度随着侧限压力(静水应力)的增加而相应增加，另外，这种材料考虑了由于屈服而引起的体积膨胀，但不考虑温度变化的影响。此材料选项适用于混凝土、岩石和土壤等颗粒状材料。对DP材料，其屈服准则表达式为：yTmeSMS21213其中：200000020000002000000100000010000001M，m＝)(31zyx，S＝偏差应力材料常数和屈服强度y的表达式如下：)sin3(3sin2，)sin3(3cos6Cy其中：为内摩擦角，C为粘滞力。对DP材料，当材料参数y,给定后，其屈服面为一圆锥面，此圆锥面是六角形的摩尔－库仑屈服面的外接锥面，如图4-5所示：图4-5Drucker-Prager屈服面4.2.1.2流动准则流动准则描述了发生屈服时，塑性应变的方向，也就是说，流动准则定义了单个塑性应变分量（xpl，ypl等）随着屈服是怎样发展的。流动准则由以下方程给出：Qdpl其中：是塑性乘子（决定了塑性应变量），Q为塑性势，是应力的函数（决定了塑性应变方向）。一般来说，流动方程是塑性应变在垂直于屈服面的方向发展的屈服准则中推导出来的，即Q等于屈服函数，这种流动准则叫作关联流动准则，如果使用其它的流动准则（从其它不同的函数推导出来），则叫作不关联的流动准则。4.2.1.3强化准则强化准则描述了初始屈服准则随着塑性应变的增加是怎样发展的。一般来说，屈服面的变化是以前应变历史的函数，在ANSYS程序中，使用了三种强化准则:等向强化是指屈服面以材料中所作塑性功的大小为基础在尺寸上扩张。对Mises屈服准则来说，屈服面在所有方向均匀扩张。见图4-6。图4-6等向强化时的屈服面变化图由于等向强化，在受压方向的屈服应力等于受拉过程中所达到的最高应力。随动强化假定屈服面的大小保持不变而仅在屈服的方向上移动，当某个方向的屈服应力升高时，其相反方向的屈服应力应该降低。见图4-7。图4-7随动强化时的屈服面变化图在随动强化中，由于拉伸方向屈服应力的增加导致压缩方向屈服应力的降低，所以在对应的两个屈服应力之间总存一个y2的差值，初始各向同性的材料在屈服后将不再是各向同性的。混合强化是等向强化和随动强化的结合，屈服面不仅在大小上扩张，而且还在屈服的方向上移动。见图4-8。图4-8混合强化时的屈服面变化图4.2.2塑性材料选项一些选项可用于描述塑性行为。用户也可以应用Ansys可编程特性在程序中加入其他特性。参见《ANSYSGuidetoUserProgrammableFeatures》。1、双线性随动强化(BKIN)选项该选项假设总应力范围等于屈服应力的二倍，以包括包辛格效应(见图4-10)。对于服从VonMisses屈服准则的一般小应变分析，建议用这一选项。不应该应用于大应变分析。用户可以用BKIN和HILL选项组合来模拟各向异性随动强化塑性。应力-应变-温度数据的示例如下。图4-9说明了这种材料选项的典型显示[TBPLOT]。MPTEMP,1,0,500!DefinetemperaturesforYoung'smodulusMP,EX,1,12E6,-8E3!C0andC1termsforYoung'smodulusInitialyieldsurfaceSubsequentyieldsurface12TB,BKIN,1,2!ActivateadatatableTBTEMP,0.0!Temperature=0.0TBDATA,1,44E3,1.2E6!Yield=44,000;Tangentmodulus=1.2E6TBTEMP,500!Temperature=500TBDATA,1,29.33E3,0.8E6!Yield=29,330;Tangentmodulus=0.8E6TBLIST,BKIN,1!Listthedatatable/XRANGE,0,0.01!X-axisofTBPLOTtoextendfromvarepsilon=0to0.01TBPLOT,BKIN,1!Displaythedatatable上面的命令MPTEMP,MP,TB,TBTEMP,TBDATA,TBLIST,/XRANGE和TBPLOT参见《ANSYSCommandsRefer