ANSYS高级非线性培训手册第08章

1、Training ManualDrucker-Prager 与混凝土与混凝土Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual 本章包括一些土木本构模型本章包括一些土木本构模型, 即即 ANSYS 中可用中可用的的 Drucker-Prager 塑性和混凝土。塑性和混凝土。 Drucker-Prager 用于颗粒状材料用于颗粒状材料, 如土壤、岩石、混如土壤、岩石、混凝土。凝土。 混

2、凝土模型用于表示脆性材料的特性混凝土模型用于表示脆性材料的特性, 包括岩石和某包括岩石和某些陶瓷材料。介绍了断裂和压碎选项。些陶瓷材料。介绍了断裂和压碎选项。Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual 本章讨论下列内容本章讨论下列内容:A. Drucker-Prager 塑性塑性B. 混凝土模型混凝土模型Advanced Structural Nonlinearities 1

3、1.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual Drucker-Prager (DP) 塑性应用于颗粒状塑性应用于颗粒状(摩擦摩擦)材材料料, 如土壤、岩石和混凝土。如土壤、岩石和混凝土。 与金属塑性不同与金属塑性不同, 对于对于DP, 屈服面是与压力有关屈服面是与压力有关的的 von Mises面面: 式中式中 s se 是修正的等效应力是修正的等效应力, s sm 是静水压力是静水压力, b b是材料常数。是材料常数。 在主应力空间画出的屈服面在

4、主应力空间画出的屈服面 是一个圆锥。是一个圆锥。 21213sMsTmebsss1s2s3s1 s2 s3Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual 一些重要的注意点一些重要的注意点: 压缩时压缩时, 静水压力的增加导致屈服强度的增加。静水压力的增加导致屈服强度的增加。 因为体积应变与静水压力有关因为体积应变与静水压力有关, 所以能考虑由于屈服所以能考虑由于屈服引起的材料的体

5、积膨胀。引起的材料的体积膨胀。 假设没有硬化假设没有硬化, 因此材料行为是弹性因此材料行为是弹性-理想塑性。理想塑性。Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual 可将屈服准则写成如下形式可将屈服准则写成如下形式. 材料参数材料参数 b b 和和 s sy被定义为被定义为 式中式中 f f 是内摩擦角是内摩擦角, , c 为粘滞力。为粘滞力。 DP 模型需要输入粘滞力模型需要输

6、入粘滞力(剪切屈服应力剪切屈服应力) “c”, 单位为单位为应力的单位。应力的单位。 还需输入内摩擦角还需输入内摩擦角“f f” ，单位为度。单位为度。 yTmsMsFsbs21213ffbsin33sin2 ffssin33cos6cyAdvanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual 注意压缩时的屈服应力大于拉伸时的屈服应力。注意压缩时的屈服应力大于拉伸时的屈服应力。若有单轴拉伸

7、若有单轴拉伸 s st 和压缩和压缩 s s c 屈服应力作为原始数屈服应力作为原始数据据, 可由下列式子将它们转换为材料参数可由下列式子将它们转换为材料参数 f f 和和 c:tctcytctcsssssssssb323ffsbbfcos6sin333233sin1ycAdvanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual 方程的简单变换说明在主应力空间原点和拉伸屈服之间方程的简单变换

8、说明在主应力空间原点和拉伸屈服之间的距离等于的距离等于(c)cot(f f) 。s1s2s3s1 s2 s3 fcotc fsffsffsbscotsin33cos6sin33sin233ccmmymAdvanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual 除了前面提及的两个参数除了前面提及的两个参数 f f 和和 c, 还有另外一个还有另外一个称为剪胀角称为剪胀角 f ff 的参数的参数

9、, 需要为需要为 DP 模型输入。模型输入。 剪胀角剪胀角 f ff 控制将要发生的体积膨胀的数量。控制将要发生的体积膨胀的数量。 颗粒在材料剪切时相互颗粒在材料剪切时相互“隆起隆起”是致密颗粒状材料的是致密颗粒状材料的一个例子。一个例子。 图示它的一种方法是在子午平面上画出屈服面。图示它的一种方法是在子午平面上画出屈服面。“p” 是静水压力是静水压力, “q” 是修正的等效应力。是修正的等效应力。pqf f f ffAdvanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced St

10、ructural Nonlinearities 11.0Training Manual 在如下的子午平面在如下的子午平面, f ff 表示塑性流动的方向表示塑性流动的方向(剪胀角剪胀角)，另一方面另一方面, f f 表示屈服面外法线的方向表示屈服面外法线的方向(内摩擦角内摩擦角)。 若若f ff=f, f, 则流动法则称为则流动法则称为关联的关联的，结果发生明显的体，结果发生明显的体积膨胀。积膨胀。 若若f ff Preprocessor Material Props Material Models Structural Nonlinear Inelastic Non-metal Plasti

11、city应输入所有的常数应输入所有的常数(即粘滞力不能为即粘滞力不能为0)，注意还需输入弹性材料属性，注意还需输入弹性材料属性(杨氏模量杨氏模量EX)， 本材料模型不考虑温度相关性。本材料模型不考虑温度相关性。TB,DP,1,0 TBDATA,1,cohesionTBDATA,2,fricangleTBDATA,3,flowangleAdvanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual

12、 Drucker-Prager 是率无关塑性模型，对求解选项是率无关塑性模型，对求解选项, 与其它与其它率无关塑性模型的考虑事项相同。率无关塑性模型的考虑事项相同。 需要时需要时, 指定非线性几何效应指定非线性几何效应(NLGEOM,ON) 指定适当的子步数以捕捉路径相关响应。指定适当的子步数以捕捉路径相关响应。 后处理考虑事项后处理考虑事项: 若材料屈服若材料屈服, 则等效塑性应变则等效塑性应变(NL,EPEQ) 为非零。为非零。 等效应力参数等效应力参数 s spl (NL,SEPL) 是在当前静水应力水平下的是在当前静水应力水平下的von Mises 等效应力等效应力: 注意对等效应变注

13、意对等效应变(EPPL,EQV), ANSYS 采取不可压缩非弹性应变采取不可压缩非弹性应变(n n=0.5)，然而然而, 若若f ff 0, 这是不真实的这是不真实的(屈服时发生体积膨胀屈服时发生体积膨胀)。当当f ff 0时时, 考虑非弹性应变的下列情况考虑非弹性应变的下列情况, 其中其中 e eeqv 应是非零的应是非零的:myplebsss330 . 0 , ,pleqvplzplyplxAAAeeeeAdvanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Struct

14、ural Nonlinearities 11.0Training Manual ANSYS 中混凝土材料模型用于模拟脆性材料中混凝土材料模型用于模拟脆性材料, 如混凝土如混凝土, 岩石和陶瓷。岩石和陶瓷。 包括断裂和压碎破坏模式。包括断裂和压碎破坏模式。 破坏前破坏前, 假设行为是线弹性的，然而假设行为是线弹性的，然而, 塑性和塑性和/或蠕变或蠕变可以与混凝土结合可以与混凝土结合, 以提供破坏前的非线性行为。以提供破坏前的非线性行为。 该本构模型适用于低拉伸强度、高压缩承载能力。该本构模型适用于低拉伸强度、高压缩承载能力。 由实常数指定沿三个单元坐标方向的由实常数指定沿三个单元坐标方向的 “分

15、布的分布的”加筋加筋, 或者由或者由LINK 或或 COMBIN 单元分别添加离散的加筋。单元分别添加离散的加筋。Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual 混凝土材料有以下特性混凝土材料有以下特性: 在单元的每个积分点上进行材料计算。在单元的每个积分点上进行材料计算。 混凝土模型在破坏点前呈现线弹性行为，在破坏混凝土模型在破坏点前呈现线弹性行为，在破坏载荷载荷(s sc o

16、r s st)下下, 发生发生压碎或断裂压碎或断裂，并且在该点，并且在该点材料完全破坏。材料完全破坏。 压碎情况压碎情况(压缩压缩)下下, 材料完全破坏。材料完全破坏。 允许在每个积分点的三个正交方向上断裂允许在每个积分点的三个正交方向上断裂(拉伸拉伸)，断裂在一个或几个方向上发生。对于发生断裂的断裂在一个或几个方向上发生。对于发生断裂的方向方向, 拉伸强度实质上变为零拉伸强度实质上变为零, 尽管断裂结束时尽管断裂结束时, 垂直于裂缝的压应力可以转移，在没有发生断裂垂直于裂缝的压应力可以转移，在没有发生断裂的方向上材料性能保持不变。的方向上材料性能保持不变。 剪切转移系数剪切转移系数 b bt

17、 和和 b bc 定义了在引起沿断裂面滑定义了在引起沿断裂面滑移的载荷作用下的剪切强度减缩系数。移的载荷作用下的剪切强度减缩系数。e es sftfcAdvanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual 混凝土材料能与其它非线性组合混凝土材料能与其它非线性组合: 混凝土可以包括塑性和蠕变。通常混凝土可以包括塑性和蠕变。通常, 多线性弹性或多线性弹性或 Drucker-Prager 塑性

18、用于混凝土。注意塑性屈服面必须位塑性用于混凝土。注意塑性屈服面必须位于混凝土破坏面于混凝土破坏面内部内部，否则不会发生屈服。，否则不会发生屈服。 右图为在主应力空间画出右图为在主应力空间画出 的混凝土破坏面。因此的混凝土破坏面。因此, 任何其它非线性材料行任何其它非线性材料行 为为(即塑性即塑性) 的屈服面必须的屈服面必须 位于混凝土破坏面内部。位于混凝土破坏面内部。 否则否则, 材料将完全破坏而材料将完全破坏而 从未屈服。从未屈服。 在断裂在断裂/压碎检查之前压碎检查之前, 进行由于塑性的应力调整。进行由于塑性的应力调整。Advanced Structural Nonlinearities

19、11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual 对于材料破坏对于材料破坏, 必须考虑四个区域必须考虑四个区域: 0 s s1 s s2 s s3(压缩压缩 压缩压缩 压缩压缩) s s1 0 s s2 s s3(拉伸拉伸 压缩压缩 压缩压缩) s s1 s s2 0 s s3(拉伸拉伸 - 拉伸拉伸 - 压缩压缩) s s1 s s2 s s3 0 (拉伸拉伸 - 拉伸拉伸 - 拉伸拉伸) 对于三维应力状态对于三维应力状态, 破坏面是主应力和下面讨

20、论破坏面是主应力和下面讨论的五个输入参数的函数，对上面所述的四个范围的五个输入参数的函数，对上面所述的四个范围, 破坏面破坏面都各不相同都各不相同。 有关方程的详情请查阅有关方程的详情请查阅 ANSYS理论手册理论手册，第，第 4.7 节。节。Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual 需要的混凝土材料常数如下需要的混凝土材料常数如下: 极限抗拉强度极限抗拉强度, ft 极限

21、抗压强度极限抗压强度, fc 极限双轴抗压强度极限双轴抗压强度, fcb 周围静水应力状态周围静水应力状态, s sah 双轴压缩与静水应力状态的叠加状态的极限抗压强度双轴压缩与静水应力状态的叠加状态的极限抗压强度, f1 单轴压缩与静水应力状态的叠加状态的极限抗压强度单轴压缩与静水应力状态的叠加状态的极限抗压强度, f2 首先需要两个常数首先需要两个常数ft 和和 fc. 其它的缺省为其它的缺省为 该近似值仅对低静水应力分量情况该近似值仅对低静水应力分量情况, 或或 有效有效. 否否则则, 用户应提供上述所有值。用户应提供上述所有值。chf3sccccbffffff725. 145. 12

22、. 121Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual 通过实常数通过实常数, 也可把加筋指定为体积比也可把加筋指定为体积比(VR)。 加筋是加筋是“分布的分布的”，若需要模拟离散加筋，若需要模拟离散加筋, 建议采用建议采用 LINK 或或 COMBIN单元。单元。 需要的输入包括钢筋材料号、体积比及方向角需要的输入包括钢筋材料号、体积比及方向角q q 和和 f f。 加筋也可

23、以包括塑性和蠕变，一般地加筋也可以包括塑性和蠕变，一般地, 对加筋采用一对加筋采用一种常用的随动强化法则。种常用的随动强化法则。 可指定多达可指定多达 3 组加筋。组加筋。 加筋体积比的总和不能超过加筋体积比的总和不能超过1.0。Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual 只有只有 SOLID65，8节点六面体单元支持混凝土。节点六面体单元支持混凝土。 Main Menu P

24、reprocessor Element Type Add/Edit/Delete 缺省时缺省时(KEYOPT(1), SOLID65 还包括特殊位移形状。还包括特殊位移形状。 后面讨论断裂选项的应力松弛。后面讨论断裂选项的应力松弛。Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual SOLID65 实常数涉及加筋实常数涉及加筋(若若存在存在): Main Menu Real Cons

25、tants Add/Edit/Delete 取向角取向角 q q 和和 f f 如下所示如下所示:Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual 通过材料通过材料 GUI 或或TB命令输入混凝土材料参数命令输入混凝土材料参数: Main Menu Preprocessor Material Props Material Models Structural Nonlinear In

26、elastic Non-metal PlasticitycahcbctctTfffff21sbbAdvanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual 可以指定多达可以指定多达 6 组温度相关的常数。组温度相关的常数。 前面讨论了常数前面讨论了常数 1-8，常数，常数 9 在后面讨论，常数在后面讨论，常数 3 或或 4 的的“-1”值分别去除断裂或压碎行为，常数值分别去除断裂或压碎行为，

27、常数 5-8 缺省值已在缺省值已在前面幻灯片讨论。前面幻灯片讨论。(左边所示为命令输入。左边所示为命令输入。)ConstantSymbolMeaning1btShear transfer coefficients for an open crack. (defaults to 1e-6)2bcShear transfer coefficients for a closed crack. (defaults to 1e-6)3ftUniaxial tensile cracking stress.4fcUniaxial crushing stress (positive).5fcbBiaxial

28、crushing stress (positive).6sahAmbient hydrostatic stress state for use with constants 7 and 8. (default is 0.0)7f1Biaxial crushing stress (positive) under the ambient hydrostatic stress state (constant 6).8f2Uniaxial crushing stress (positive) under the ambient hydrostatic stress state (constant 6)

29、.9TcStiffness multiplier for cracked tensile condition, used if KEYOPT(7) = 1 (defaults to 0.6).TB,CONC,1,1,9, TBTEMP,0TBDATA,1,ShrCf-OpTBDATA,2,ShrCf-ClTBDATA,3,UnTensStTBDATA,4,UnCompStTBDATA,5,BiCompStTBDATA,6,HydroPrsTBDATA,7,BiCompStTBDATA,8,UnTensStTBDATA,9,TenCrFacAdvanced Structural Nonlinea

30、rities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual 若存在加筋若存在加筋, 它们的方向通过它们的方向通过 GUI 显示显示: Utility Menu PlotCtrls Device Options Vector mode ON Utility Menu PlotCtrls Style Size and Shape Display of element shapes based on real constant descriptions

31、 ON Utility Menu Plot Elements或通过命令或通过命令: /DEV,VECTOR,1 /ESHAPE,1 EPLOT钢筋方向为红色所示。钢筋方向为红色所示。Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual 求解后求解后, 可以画出断裂可以画出断裂: Main Menu General Postproc Plot Results Concrete Plot

32、Crack/Crush 或通过命令或通过命令: PLCRACK也能得到其它项也能得到其它项, 如状态如状态(未破坏的未破坏的, 压碎压碎, 开口裂缝开口裂缝, 闭合裂缝闭合裂缝), 断裂方向角断裂方向角 及加筋解也可以获得。及加筋解也可以获得。右图中右图中, 注意断裂方向注意断裂方向 和平面在每个积分和平面在每个积分 点上画出。点上画出。Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manu

33、al混凝土分析的提示混凝土分析的提示: 确保混凝土常数对于特定的应用是有效的。例如确保混凝土常数对于特定的应用是有效的。例如, 若有大的静水应力分量若有大的静水应力分量(典型的大混凝土结构典型的大混凝土结构, 如如大坝大坝), fcb, f1和和 f2 的缺省的计算值是的缺省的计算值是不不合适的。合适的。 类似地类似地, 零抗拉强度会引起收敛困难且在物理上零抗拉强度会引起收敛困难且在物理上不真实。不真实。 若混凝土结构在加载区域完全破坏若混凝土结构在加载区域完全破坏, 因为刚度为因为刚度为零零, 所以不收敛，这说明几何不稳定性所以不收敛，这说明几何不稳定性(见第见第 9 章章)。 为确保正确的

34、载荷转移为确保正确的载荷转移, 若材料断裂或压碎若材料断裂或压碎, 则则需要使用大量的子步数。需要使用大量的子步数。Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Advanced Structural Nonlinearities 11.0Training Manual混凝土分析提示混凝土分析提示: 对于断裂问题对于断裂问题, 使用使用KEYOPT(7)=1 有助于收敛，这是断有助于收敛，这是断裂后的应力松弛。裂后的应力松弛。 破坏后破坏后, 材料刚度突然下降为零材料刚度突然下降为零(左图左图)。 输入常数输入常数 Tc 作为混凝土材料特性的第九个常数，这是一个系数作为混凝土材料特性的第九个常数，这是一个系数(缺省值缺省值=0.6), 它作为应力松弛的一个乘子。它作为应力松弛的一