ansys高级非线性分析四蠕变（课堂PPT）

蠕变,四章,S,隐式和显式蠕变本章综述,前一章探讨了率无关塑性,本章将讨论ANSYS中蠕变行为的分析方法。尽管从材料的观点看蠕变和粘塑性是相同的,但本构模型的使用不同。因此,率相关的塑性主题被分为两部分,这部分属于蠕变。本章将介绍ANSYS中可用的大量的隐式和显式蠕变法则。主要讨论金属的蠕变。然而,讨论的各种观点也适用于塑料或混凝土等其它材料的蠕变。ANSYS有隐式和显式两种蠕变分析过程。首先讨论蠕变的一般知识，然后是进行隐式或显式蠕变分析的细节。,S,隐式和显式蠕变本章综述,本章包括下列主题：A.蠕变现象背景B.术语定义C.一般蠕变方程D.隐式蠕变过程E.显式蠕变过程F.ANSYS蠕变模型的求解过程G.隐式蠕变与显式蠕变的比较,S,隐式和显式蠕变A.蠕变背景,晶体材料中,如金属,蠕变机理与空隙的扩散流动和位错运动有关。空隙是点缺陷,倾向于形成与所施加应力方向垂直(而不是平行)的晶界。空隙由高集中区向低集中区运动。在低应力状态下发生扩散流动，但通常需要高温条件。晶粒的位错是线缺陷.位错运动(攀升、滑动、偏移)在高应力状态下被激活,尽管在中温时也可能发生位错运动。有时晶界滑移被认为是一种独立的导致蠕变变形的机理。,S,隐式和显式蠕变.蠕变背景,尽管对材料科学的详细论述超越了本书的范围，也足以说明前面提及的物理机理导致蠕变。蠕变变形与应力、应变、时间和温度的相关性一般用与下式相似的形式来模拟：函数f1-4与选择的蠕变法则有关。通常由不同应变速率和温度条件下的各种拉伸实验确定相关蠕变常数。假设各向同性行为,vonMises方程用于计算有效应力,在蠕变应变率方程中使用等效应变(与率无关塑性相似)。,S,隐式和显式蠕变.蠕变背景,当计算弹性、塑性和蠕变应变时,ANSYS使用附加的应变分量:计算塑性应变(流动法则)类似于前一章的描述。根据蠕变应变率方程计算蠕变应变，它的详细形式将在后面讨论。弹性、蠕变和塑性应变都基于(当前的)应力状态，但它们是独立计算的(彼此互不依赖)。注意在利用隐式蠕变与显式蠕变进行计算时有差别。,S,隐式和显式蠕变.蠕变背景,与塑性类似,蠕变是一种基于偏量行为的不可逆(非弹性)应变。在蠕变流动条件下,假设材料是不可压缩的。另一方面,与率无关塑性不同，蠕变在发生非弹性应变时没有屈服面。因此,蠕变不需要高应力值来发生更多的蠕变应变。假设在所有非零应力值时都会发生蠕变应变。前面提到,从材料的角度看,蠕变和粘塑性是相同的。工程应用中,通常蠕变用于描述低应变速率的热激活过程。率无关塑性和隐式蠕变应变以弱耦合方式处理。相反,ANSYS中粘塑性本构模型用于描述高应变速率的应用(例如,冲击载荷)。非弹性应变以强耦合方式处理。,S,隐式和显式蠕变.蠕变分析实例,焊球蠕变分析的例子(热循环)。,BretZahn,ChipPAC授权的ANSYS模型(,单元185(B-Bar),双曲正弦隐式蠕变模型,S,隐式和显式蠕变B.术语的定义,蠕变的三个阶段:在常载荷条件下,蠕变的单轴应变与时间的关系如下图所示。第一阶段,应变率随时间减小，该阶段在较短时间内完成。第二阶段具有常应变率。第三阶段,应变率迅速增加直到破坏(断裂)。,S,隐式和显式蠕变.术语的定义,蠕变的三个阶段(续)蠕变应变率可能是应力、应变、温度、和/或时间的函数。对于工程分析,通常最关心蠕变的第一和第二阶段，第三阶段蠕变通常与开始破坏有关(颈缩、损坏),且时间较短，因此,ANSYS中不模拟第三阶段。第一阶段蠕变的应变率通常远大于第二阶段蠕变。然而,应变率在第一阶段逐渐降低而在第二阶段几乎为常值(对于前面提及的常应力、恒温下单轴试验情况而言)，而且第一阶段蠕变时间比第二阶段短。,S,隐式和显式蠕变.术语的定义,蠕变在施加恒应力条件下,蠕变应变增加。应力松弛在施加恒应变条件下,应力降低。,S,隐式和显式蠕变.术语的定义,时间强化假设蠕变应变率仅与从蠕变过程开始的时间有关，也就是说,该曲线上/下移动，当应力从s1到s2变化时,计算A到B点的不同蠕变速率。应变强化假设蠕变速率仅与材料现有的应变有关，也就是说,曲线左/右移动.当应力从s1到s2变化时,计算A到B点的不同蠕变应变率。,S,隐式和显式蠕变.术语的定义,显式蠕变显式蠕变是指应用Euler向前法进行蠕变应变演化的计算。在每个时间步使用的蠕变应变率与该时间步开始时的速率一致，并假设在整个时间步Dt内为常量，因此，需要很小的时间步以使减小误差。对于有塑性的显式蠕变,首先进行塑性修正,然后进行蠕变修正。两种修正在不同应力值处进行，因此精度较差。,S,隐式和显式蠕变.术语的定义,隐式蠕变隐式蠕变应用了Euler向后积分法求解蠕变应变。该方法在数值上无条件稳定，这意味着不必象显式蠕变方法那样,使用小的时间步,所以总体上会更快。对于隐式蠕变加上率无关塑性,塑性修正和蠕变修正同时进行,而不是分别进行。因此,隐式蠕变一般比显式蠕变更精确,但它仍与时间步大小有关，必须使用足够小的时间步来精确捕捉路径相关行为。基于上述原因,隐式蠕变是ANSYS推荐的方法(高效、精确)。后面将详细讨论两种蠕变过程。,S,隐式和显式蠕变C.一般蠕变方程,如前所述,蠕变方程通常是一种速率形式,类似于下式:然而,被分析的材料类型决定了具体蠕变方程的选择。现在讨论一些一般的特性,具体模型包括在隐式蠕变和显式蠕变两部分中。单元手册第2.5章节也包括了隐式和显式蠕变方程。,S,隐式和显式蠕变.一般蠕变方程,温度相关性蠕变效应被热激活,它的温度相关性通常由Arrhenius法则来表示:式中Q为激活能,R为普适气体常数,T为绝对温度。应力相关性蠕变应变通常也与应力有关,尤其是位错蠕变。Norton法则为:对上述幂定律的常见修正如下:,S,隐式和显式蠕变.一般蠕变方程,通常第一阶段蠕变显示时间或应变强化。时间强化包含在一个时间相关项中:应变强化包含在一个应变相关项中:由可用的材料数据来决定使用哪一项(应变或时间强化)。第二阶段不具有时间或应变强化，第二阶段的蠕变应变率通常是常数。,S,隐式和显式蠕变.一般蠕变方程,ANSYS中可用的蠕变法则汇总如下:如前所述,只要可能,推荐采用隐式蠕变。,隐式蠕变,D节,S,隐式和显式蠕变D.隐式蠕变过程,本节讨论进行隐式蠕变分析的过程。如前所述,因为隐式蠕变比显式蠕变更稳定、精确和高效(更快),所以是首选的方法。这是由于该公式中使用了隐式积分法(Euler向后法)。有些情况下采用的蠕变法则或单元类型需要使用显式蠕变方法。,S,隐式和显式蠕变.支持的单元类型,隐式蠕变材料支持的单元类型:“核心”单元:PLANE42,SOLID45,PLANE82,SOLID92和SOLID95“18x”系列单元:LINK180,SHELL181,PLANE182,PLANE183,SOLID185,SOLID186,SOLID187,BEAM188和BEAM189隐式蠕变分析中,推荐选择“18x”系列单元。因为在“18x”系列单元中有大量的单元技术可用(参考第二章),这些单元提供了更强的灵活性和能力。这些公式包括B-Bar,URI,增强应变及混合U-P。“18x”系列单元比核心单元支持更多的本构模型,包括超弹性。,S,隐式和显式蠕变.支持的塑性模型,前已述及,蠕变与率无关塑性解耦。隐式蠕变允许与下述率无关塑性模型相结合(见第三章):有蠕变的BISO,MISO和NLISO(等向强化蠕变)有蠕变的BKIN(随动强化蠕变)蠕变与HILL(各向异性蠕变)蠕变与BISO,MISO和有HILL的NLISO(等向强化的各向异性蠕变)蠕变与有HILL的BKIN(随动强化的各向异性蠕变),S,隐式和显式蠕变.定义隐式蠕变模型,可以用命令或通过GUI定义隐式蠕变模型(在下面的幻灯片中讨论)。对于隐式蠕变,以两种方式定义温度相关性:由TBDATA(或材料GUI)定义温度相关常数很多蠕变方程包括前面提及的Arrhenius方程由用户确定使用任一个或两个方法来包含温度相关性。,S,隐式和显式蠕变.定义隐式蠕变模型,通过命令定义隐式蠕变时,使用TB,CREEP,具体蠕变模型由TBOPT定义。TB,CREEP,mat,ntemp,npts,TBOPTTBTEMP定义与温度相关的常数TBDATA定义实际常数对下面的例子,TBOPT=2是指采用时间-强化蠕变方程.使用TBTEMP命令指定温度相关常数,与该方程相关的四个常数由TBDATA命令指定为自变量。TB,CREEP,1,1,4,2TBTEMP,100TBDATA,1,C1,C2,C3,C4,S,隐式和显式蠕变.定义隐式蠕变模型,所有隐式蠕变模型可以在如下的材料GUI中选择:StructuralNonlinearInelasticRateDependentCreep确保首先定义必需的线弹性材料属性(EX和PRXY)。,S,隐式和显式蠕变.定义隐式蠕变模型,当选择合适的隐式蠕变模型后,出现一个独立对话框显示所需的输入项。下面例子中,已经定义了第一阶段蠕变方程,提示用户输入四个蠕变常数。也可以输入温度相关的常数。,S,隐式与显式蠕变.可用的隐式蠕变模型,下表显示可用的隐式蠕变法则，方程在下面幻灯片中说明。,S,隐式和显式蠕变.可用的隐式蠕变模型,1)应变强化TBOPT=1第一阶段蠕变2)时间强化TBOPT=2第一阶段蠕变3)广义指数TBOPT=3第一阶段蠕变4)广义GrahamTBOPT=4第一阶段蠕变,S,隐式和显式蠕变.可用的隐式蠕变模型,5)广义BlackburnTBOPT=5第一阶段蠕变6)修正的时间强化TBOPT=6第一阶段蠕变7)修正的应变强化TBOPT=7第一阶段蠕变,S,隐式和显式蠕变.可用的隐式蠕变模型,8)广义GarofaloTBOPT=8第二阶段蠕变9)指数形式TBOPT=9第二阶段蠕变10)NortonTBOPT=10第二阶段蠕变,S,隐式和显式蠕变.可用的隐式蠕变模型,11)时间强化TBOPT=11第一阶段+第二阶段12)有理多项式TBOPT=12第一阶段+第二阶段13)广义时间强化TBOPT=13第一阶段蠕变,S,隐式和显式蠕变.练习,请参考附加练习题:练习6:应力松弛,显式蠕变,E节,S,隐式和显式蠕变E.显式蠕变过程,本节讨论进行显式蠕变分析的过程。前已述及由于隐式蠕变方法比显式蠕变更有效和精确而成为首选方法。显式蠕变使用需要非常小的时间步的Euler向前法,因此会有很多次迭代。与隐式蠕变不同,塑性计算不是同时进行的。首先进行塑性分析,然后是蠕变计算(叠加)，该时间步的塑性应变等不进行重新调整。只要可能,都应使用隐式蠕变，然而,有些情况下采用的蠕变法则或单元类型需要使用显式蠕变。,S,隐式和显式蠕变.显式蠕变过程,显式蠕变材料支持的单元类型:“核心”单元:PLANE42,SOLID45,PLANE82,SOLID92和SOLID95其它单元:LINK1,PLANE2,LINK8,PIPE20,BEAM23,BEAM24,SHELL43,SHELL51,PIPE60,SOLID62和SOLID65注意18x系列单元不支持显式蠕变。显式蠕变支持的塑性基本模型:任何使用中的单元类型所支持的塑性模型都可以与显式蠕变相结合(如BISO,MISO,BKIN,KINH/MKIN,DP)。前已述及,这是塑性和蠕变的非同时性模拟(首先进行塑性计算,然后是蠕变计算)。,S,隐式和显式蠕变.显式蠕变过程,可以使用命令或通过GUI定义一个显式蠕变模型(在下面的幻灯片中讨论)。显式蠕变没有与温度有关的常数。由蠕变方程说明温度相关性。显式蠕变常数定义并输入为C1,C2等等,其中C1为第一个常数,C6为第六个常数等。不必定义所有的常数，需要使用的常数数目与选择的蠕变方程有关。,S,隐式和显式蠕变.显式蠕变过程,当通过命令定义显式蠕变时,使用TB,CREEP命令,TBOPT=0(或保留空白)。TB,CREEP,mat,ntemp,nptsTBDATA定义实际的常数。通常用C6常数指定第一阶段蠕变(选择C6=0到15)，若C10或T+Toffset0,则不计算第一阶段蠕变。用C12常数选择第二阶段蠕变(C12=0或1)，第一阶段蠕变方程C6=9-11,13-15已经包括第二阶段蠕变效应,所以第二阶段蠕变C12常数被忽略，若C70或T+Toffset0,则不计算第二阶段蠕变。用C66常数指定辐射感应蠕变，C6=0至11时使用该方程，若C550且C610,或T+Toffset0,则不计算第二阶段蠕变。,S,隐式和显式蠕变.定义显式蠕变模型,如下的材料GUI中可以选择所有的显式蠕变模型:StructuralNonlinearInelasticRateDependentCreep确保首先定义必需的线弹性材料属性(EX和PRXY)。,S,隐式和显式蠕变.定义显式蠕变模型,当选择合适的显式蠕变模型后,出现一个独立的对话框显示需要的输入项。下面例子中,已经定义了一个蠕变方程,提示用户输入各种蠕变常数。,S,隐式和显式蠕变.可用的显式蠕变模型,下表为可用的显式蠕变法则，其方程将在下面的幻灯片说明。注意对于具体材料,很多蠕变法则有固有的单位(例如C6=11)。请查阅单元手册的第2.5章节以保证使用正确的单位系统。,S,隐式和显式蠕变.可用的显式蠕变模型,1)应变强化C6=0第一阶段蠕变2)时间强化C6=1第一阶段蠕变3)广义指数C6=2第一阶段蠕变4)钢的蠕变法则C6=9,10,11,14,or15第一和第二阶段蠕变参阅单元手册第2.5章中的定义ec的选项。,S,隐式和显式蠕变.可用的显式蠕变模型,5)幂函数蠕变法则C6=12第一阶段蠕变6)真幂函数C6=13第一和第二阶段蠕变7)指数形式C12=0第二阶段蠕变8)NortonC12=1第二阶段蠕变,蠕变求解过程,F节,S,隐式和显式蠕变F.求解蠕变问题,前面讨论了ANSYS中隐式和显式蠕变的一些区别。可用的蠕变法则取决于蠕变方法蠕变常数的输入不同支持的单元类型不同隐式蠕变是首选方法本节将回顾蠕变材料模型的求解过程:求解选项后处理适当的时候将提示隐式和显式蠕变求解过程的区别。,S,隐式和显式蠕变.蠕变的求解选项,含蠕变材料的模型求解与其他非线性问题类似,但求解蠕变问题时有一些特殊的考虑事项。蠕变是大应变还是小应变取决于该问题。若不能确定,最好打开大变形效应。与其它率无关材料的静态非线性分析不同,蠕变分析中“时间”有重要意义。确保结束时间对于模型和所关心的时间范围是合适的。注意分析不是必需为瞬态分析，惯性效应(TIMINT)打开或关闭取决于该问题。然而,一般来说,由于时间范围较长,蠕变分析不考虑惯性效应(ANTYPE,STATIC或TIMINT,OFF)。,S,隐式和显式蠕变.蠕变的求解选项,指定蠕变分析的求解类型MainMenuSolution-AnalysisType-SolnControlSolutionControls-BasicTab-AnalysisOptions需要时指定大位移求解（NLGEOM,ON）推荐使用缺省的求解控制设置(SOLCONTROL)，缺省时,求解控制打开。,S,隐式和显式蠕变.蠕变的求解选项,在载荷步末使用合适的时间(TIME)。SolutionControls-BasicTab-TimeControl与其它静态分析不同,因为蠕变是率相关的,在蠕变分析中TIME有其含义。分析可以是静态或瞬态的(不包括或包括惯性效应)，但TIME应该是真实单位。指定一个合适的初始时间步长(DELTIM)来捕捉非线性效应，确保最小和最大时间步长也合理。,右面例子中,在求解控制对话框中,指定结束时间为1300.初始时间步是长1.0,最小为0.5,最大为10。(仅供示范),S,隐式和显式蠕变.切换隐式蠕变效应,仅对于隐式蠕变,RATE(“包括应变率效应”)命令在分析中用于打开或关闭蠕变效应。SolutionControls-NonlinearTab-CreepOption这对建立初始条件有用。这种情况下,应该设定一个非常小的结束时间值(如1e-8),关闭速率效应(RATE,OFF)，正常求解，然后用RATE,ON命令打开蠕变效应,并指定真实结束时间。RATE命令仅适用于下述情况:18x系列单元的隐式蠕变(vonMises势)18x或核心单元的各向异性隐式蠕变(Hill势),S,隐式和显式蠕变.蠕变应变比率,因为蠕变是一种路径相关现象,确保充分地捕捉响应是重要的。ANSYS采用的一个措施是蠕变比率Cs,Cs定义为:式中,Decr为等效蠕变应变增量,eet为修正的等效弹性应变(详见理论手册第4.2/4.3章)。,S,隐式和显式蠕变.蠕变应变比率,需要时使用后退控制(CUTCONTROL)指定等效蠕变应变比率的最大值。SolutionControls-NonlinearTab-CutbackControlCUTCONTROL,CRPLIMIT,crvalue,1将为隐式蠕变施加蠕变比率的最大值-crvalue。缺省时不指定隐式蠕变极限控制。CUTCONTROL,CRPLIMIT,crvalue,0将为显式蠕变施加蠕变比率的最大值-crvalue。缺省时显式蠕变极限为10%。在一个时间步长中,若ANSYS计算的一个蠕变应变比率大于crvalue,则求解步自动二分,直到满足蠕变极限或达到最小时间步长。,S,隐式和显式蠕变.蠕变应变比率,如前所述,显式蠕变有一个稳定性极限，这与蠕变极限0.25(25%)一致,所以指定的极限不能大于该值，使用时,缺省值0.10(10%)比较合适。另一方面,隐式蠕变是绝对稳定的,所以没有稳定性极限，因而,缺省时不强加极限。然而,这不意味着隐式蠕变是“绝对精确。”一般推荐的蠕变比率极限为1到10(100-1000%),并确保指定一个足够小的初始时间步长及最小、最大时间步长。蠕变比率被打印在如下所示的输出文件/窗口:*LOADSTEP1SUBSTEP94COMPLETED.CUMITER=94*TIME=940.