《1 前言》
1 前言
多数岩土工程都处于弹塑性状态, 因而岩土塑性在岩土工程的设计中至关重要。早在1773年Coulomb提出了土体破坏条件, 其后推广为Mohr-Coulomb条件。1857年Rankine研究了半无限体的极限平衡, 提出了滑移面概念。1903年Kotter建立了滑移线方法。Fellenius (1929) 提出了极限平衡法。以后 Terzaghi、Sokolovskii又将其发展形成了较完善的岩土滑移线场方法与极限平衡法。1975年, W. F. Chen在极限分析法的基础上又发展了土的极限分析法, 尤其是上限法。不过上述方法都是在采用正交流动法则的基础上进行的。滑移线法与极限分析法只研究力的平衡, 未涉及土体的变形与位移。20世纪50年代开始, 人们致力于岩土本构模型的研究, 力求获得岩土塑性的应力-应变关系, 再结合平衡方程与连续方程, 从而求解岩土塑性问题。70年代发现, 用一个塑性势面和屈服面很难使计算结果与实际吻合;采用正交流动法则既不符合岩土实际情况, 还会产生过大的体胀。由此, 双屈服面与多重屈服面模型
真正的土力学必须建立在符合土本身特性的本构模型的基础上, 而本构模型的建立必须有符合岩土材料变形机制的建模理论。岩土塑性力学是一门新兴学科, 也是建立岩土本构模型的基础。
《2 经典塑性力学用于岩土材料问题》
2 经典塑性力学用于岩土材料问题
《2.1 岩土类材料的特点》
2.1 岩土类材料的特点
岩土类材料是由颗粒材料堆积或胶结而成, 属于摩擦型材料。摩擦材料的特点是抗剪强度中含有摩擦力项, 它的抗剪强度随压应力的增大而增大, 因而岩土材料的屈服条件与金属材料明显不同, 称此为岩土的压硬性, 即随压应力的增大岩土的抗剪强度与刚度增大。
岩土为多相材料, 岩土颗粒间有孔隙, 在各向等压作用下, 岩土颗粒间的水、气排出, 产生塑性体变, 出现屈服, 称为岩土的等压屈服特性。金属材料在各向等压作用下是不会产生塑性体变的。
岩土的体应变还与剪应力有关, 即在剪应力的作用下岩土会产生塑性体变 (剪胀或剪缩) , 一般称为岩土的剪胀性 (含剪缩) 。这在力学上表现为球张量与偏张量的交叉作用, 即球应力会产生剪变 (负值) , 这也是压硬性的一种表现;反之, 剪应力会产生体变。纯塑性金属材料是不具有这一特性的。
《2.2 经典塑性力学用于岩土类材料出现的问题》
2.2 经典塑性力学用于岩土类材料出现的问题
岩土塑性力学脱胎于经典塑性力学, 然而经典塑性力学只适应于金属材料, 当用于岩土类摩擦材料时就会出现一些不符合实际的情况, 理论计算结果与土工试验结果出现诸多矛盾
1) 按照经典塑性力学中的传统塑性势理论, 塑性应变增量的方向唯一取决于应力状态, 而与应力增量无关。文献
2) Poorooshasb
3) 基于传统塑性位势理论的单屈服面模型, 用莫尔-库仑一类剪切屈服面作屈服面时, 若用关联流动法则, 将会导致远大于实际的剪胀变形。
4) Matsouka等人试验证实
上述岩土变形机制与经典塑性力学相矛盾。
传统塑性位势理论是经典塑性力学的核心, 可表示为
式中 dε
式 (2) 是传统塑性位势理论的基本特征, 即各塑性应变增量分量存在比例关系。由此还可推证塑性主应变增量与主应力增量的关系为
塑性矩阵Ap元素a1i, a2i, a3i (i=1, 2, 3) 有
式 (2) 和式 (4) 表明, 各塑性主应变增量或Ap中的各行元素成比例, 只有一个势函数就可求出3个塑性主应变或Ap, 这是传统塑性位势理论的特点。由于塑性应变增量分量互成比例, 因而塑性应变增量方向不随塑性应变增量分量大小而变, 导致传统塑性位势理论与岩土材料变形机制的矛盾。
在德鲁克公设提出前后, 经典塑性力学一直沿用关联流动法则, 即塑性势函数与屈服函数相同。实际土工试验表明, 岩土材料不服从关联流动法则, 它只是一种假设。图2的简单摩擦系统, 在一定程度上也说明岩土类材料不符合正交流动法则
德鲁克公设一直是关联流动法则的理论支柱, 本来是作为弹塑性稳定材料的定义提出来的, 因此不是所有客观材料的力学行为都必须满足这个公设所导出的结论
在经典塑性力学中, 将屈服函数写成3个主应力或3个应力张量不变量的函数, 这就忽略了应力增量中3个剪应力增量dτij所引起的塑性变形。即经典塑性力学中不考虑应力主轴的旋转, 假设应力主轴与应力增量主轴始终共轴, 只有主应力增量dσ1, dσ2, dσ3, 而dτ12=dτ23=dτ13=0。实际岩土工程中, 应力主轴会发生旋转, 即存在主轴旋转的应力增量分量dτij, 并由此产生相应的塑性变形。不考虑应力主轴旋转也是经典塑性力学的一个假设, 无法算出应力主轴旋转产生的塑性变形。
《2.3 经典塑性力学特有的3条假设》
2.3 经典塑性力学特有的3条假设
经典塑性力学属于连续介质力学范畴, 一般具有各向同性、均质、连续、小变形等基本假设。此外, 经典塑性力学还具有下面3条特有的假设:
1) 假设应力空间中只存在一个满足式 (1) 的塑性势函数, 导致塑性应变增量分量成比例;塑性应变增量方向只与应力有关, 而与应力增量无关。
2) 假设应力与应力增量主轴共轴, 不考虑应力主轴旋转。
3) 材料服从关联流动法则。
上述假设, 不适应岩土材料的变形机制。消除上述假设, 将经典塑性力学改造成更一般的塑性力学, 称它为广义塑性力学
《3 不计应力主轴旋转的广义塑性位势理论》
3 不计应力主轴旋转的广义塑性位势理论
广义塑性力学从固体力学原理直接导出广义塑性位势理论。在不计应力主轴旋转情况下 (应力主轴、应力增量主轴及塑性应变增量主轴共轴) , 引入张量定律, 从理论上导出了广义塑性位势理论
式中 σk与εpk分别为3个主应力和3个塑性主应变。当张量Aij (如dε
则式 (5) 得证。
根据梯度定义, 有
式中 Qk为3个线性无关的任意势函数;dλk为3个塑性因子。将式 (7) 代入式 (5) , 则有
式 (8) 称为不计应力主轴旋转的广义塑性位势理论, 它与传统塑性位势理论有如下3点区别:
1) 广义塑性位势理论有3个塑性势面, 且3个塑性势面必须线性无关;而传统塑性位势理论只有一个塑性势面。
2) 广义塑性位势理论中, 塑性应变增量方向由3个塑性应变增量分量的方向和大小来定, 而3个分量既与塑性势面有关, 又与屈服面及应力增量有关。传统塑性位势理论是其特例, 此时塑性应变增量分量成比例, 塑性应变增量方向由此势函数唯一地确定, 而与应力增量无关。由此表明, 传统塑性力学中可事先确定塑性应变增量总量的势面。而广义塑性力学中, 因塑性应变增量总量方向与应力增量有关, 无法事先确定塑性应变增量总量方向 (即势面) 。但可事先确定塑性应变增量的3个分量方向, 即知3个分量的势面。
3) 3个塑性因子dλk (k=1, 2, 3) 不要求都大于零或等于零。dλk与屈服面有关, 当屈服面与塑性势面同向, dλk>0;屈服面与塑性势面反向, 则dλk<0。岩土材料的体积屈服面既可与塑性势面同向 (体缩) , 也可与塑性势面反向 (体胀) 。而传统塑性力学中只有一个塑性势面, 因而dλ一定大于零或等于零。
式 (8) 中3个塑性势函数是可任选的, 但必须保持线性无关, 最符合这一条件并应用最方便的是选用主应力空间中的3个座标轴作塑性势函数, 如选σ1, σ2, σ3或p, q, θσ不变量为势函数。这种情况下构造屈服函数也最为方便。这说明势函数可采用任何一种形式的3个张量不变量。
取σ1, σ2, σ3的等值面为3个塑性势函数即Q1=σ1, Q2=σ2, Q3=σ3时, 式 (8) 变为
式中 dλ1, dλ2, dλ3分别为相应上述3个势面的塑性因子, 将σ1=Q1, σ2=Q2, σ3=Q3代入式 (9) 或按其物理意义均能得到
可见dλk具有明确的物理意义。
如果取p, q, θσ为塑性势函数, 有
同理有
式中 dε
塑性应变增量可分解为塑性体应变增量与塑性剪应变增量
塑性剪应变可分为q方向上的塑性剪应变增量d
从实际情况来看, 无论是岩土或金属材料, d
这就是国内常用的“南水”双屈服面模型
对于金属材料, dε
《4 塑性势面与屈服面的关系》
4 塑性势面与屈服面的关系
经典塑性力学中采用总量势面;而广义塑性力学必须采用分量势面。广义塑性力学采用了塑性力学中的分量理论。屈服面主要用来确定塑性应变增量的大小。经典塑性力学中确定塑性因子dλ;而广义塑性力学中确定3个塑性因子dλk。
塑性应变增量矢量的方向由塑性势面确定, 而大小按其相应的屈服面确定。这就表明塑性势面与屈服面必须相关。例如求塑性应变增量分量dε
屈服面一般应由试验确定。在等向硬化模型的情况下, 体积屈服面、q方向剪切屈服面与θσ方向剪切屈服面可表达为如下形式
微分式 (15) , 可得
式中
当塑性应变总量与应力具有唯一性关系时, 即上述A1=A2=A3=1时, 并略去θσ对ε
微分式 (17) 得
由上看出, 塑性势面与屈服面存在如下关系:
1) 塑性势面可以任取, 但必须保证各势面间线性无关, 屈服面则不可任取, 它必须与塑性势面相对应, 并有明确的物理意义, 如取σ1为势面, 则对应的屈服面必为塑性主应变ε
2) 取σ1, σ2, σ3或p, q, θσ为塑性势面, 相应的屈服面最简单, 并具有明确的物理意义, 即为3个塑性主应变的等值面或为塑性体应变、q方向塑性剪切应变与θσ方向塑性剪应变的等值面。
3) 由于3个塑性势面线性无关, 则相应的3个屈服面也必然互相独立。例如体积屈服面与q方向上及θσ方向上的剪切屈服面都各自独立。这表明体积屈服面只能用来计算塑性体积变形, 而与塑性剪切变形无关, 反之亦然。
4) 通常的剪胀, 指剪应力dq, dθσ所引起的体胀, 即式 (16) 中的
5) 对于采用米赛斯屈服条件的金属材料, 式 (16) 中只保留 (∂fq/∂q) dq一项, 其余均为零。
《5 岩土材料的加载条件 (屈服条件)》
5 岩土材料的加载条件 (屈服条件)
广义塑性力学中, 相应3个塑性势函数有3个屈服条件:体积屈服条件、q方向及θσ方向的剪切屈服条件。
《5.1体积屈服面》
5.1体积屈服面
以ε
5.1.1 压缩型体积屈服面
松砂、正常固结土等土体, 受力后土体体积压缩。其体积屈服面常用的是椭圆型曲线, 如图4a所示。其表达式为
式中η=q/p, M为极限状态线的斜率, pc为加载面与p轴的右交点。
5.1.2 硬化压缩剪胀型体积屈服面
中密砂、弱超固结土等土体, 应力应变曲线处于应变硬化状态, 土体体变先压缩后剪胀。这类屈服面一般应用不多, 近年来有所进展。通过中密砂试验
5.1.3 软化压缩剪胀型体积屈服面
密砂、超固结土、岩石等岩土体, 应力应变曲线先处于应变硬化状态, 后处于应变软化状态, 其体变也是先压缩后剪胀。这类岩土的体积屈服面目前研究不多, Hvorslev面可认为是软化压缩剪胀型体积屈服面。图4c上人为地把极限状态线与状态变化线合为一条, 而实际试验是2条曲线, 这是其不足之处。
《5.2剪切屈服面》
5.2剪切屈服面
剪切屈服面是以
式中β, α1, k为与γp有关的系数, n=2或1,
从试验结果看, 多数岩土的剪切屈服曲线在子午平面上是双曲线或抛物线。当加载面Φγ发展到与
表1 不同屈服条件的系数值
Table 1 Value of coefficients for various yield conditions
《表1》
屈服条件 | n | β | α1 | k | g (θσ) | ||
广义米赛斯条件 | 米赛斯条件 | 1 | 0 | 0 | c | 1 | |
外角外接圆锥 | 1 | 0 | tgˉφ | ˉc | 1 | ||
内角外接圆锥 | 1 | 0 | 6sinφ√3(3+sinφ) | 6ccosφ√3(3+sinφ) | 1 | ||
内切圆锥 (德鲁克-普拉格条件) | 1 | 0 | 3sinφ√3√3+sin2φ | 3ccosφ√3√3+sin2φ | 1 | ||
Dimaggio-sandler | 1 | 0 | (a-d)pexp(-pclnaa-d) | c | 1 | ||
屈瑞斯卡条件 | 1 | 0 | 0 | 2c/31/2 | 31/2/2cos θσ | ||
广义屈瑞斯卡条件 | 1 | 0 | 按广义米赛斯条件确定 | 31/2/2cos θσ | |||
莫尔-库仑条件 | 1 | 0 | tgˉφ | ˉc | 3-sinφ2(√3cosθσ-sinθσsinφ) | ||
莫尔-库仑等面积圆锥 | 1 | 0 | 6√3sinφ√2√3(9-sin2φ) | 6√3ccosφ√2√3(9-sin2φ) | 1 | ||
平面应变下关联流动法则莫尔-库仑匹配条件 | 1 | 0 | 3sinφ√3√3+sin2φ | 3ccosφ√3√3+sin2φ | 1 | ||
平面应变下非关联流动法则莫尔-库仑匹配条件 | 1 | 0 | sin φ | c cos φ | 1 | ||
双剪应力条件 | 1 | 0 | 0 | 2c/31/2 | [ (31/2cos θσ) /2- (sin θσ) /2]-1, θσ≤0° [ (31/2cos θσ) /2+ (sin θσ) /2]-1, θσ≥0° | ||
广义双剪应力条件 | 1 | 0 | tgˉφ | ˉc | K[ (31/2cosθσ) /2- (sinθσ) /2]-1, θσ≤tg-1 (3-1/2sinφ) | ||
1 | 0 | tgˉφ | ˉc | K[ (31/2cos θσ) /2+ (sin θσ) /2]-1, θσ≥tg-1 (3-1/2sin φ) | |||
辛克维兹-潘德条件 | 双曲线型 | 2 | -tg2ˉφ | 2ˉc | ˉc2-a2tgˉφ | 2K/[ (1+K) - (1-K) sin 3θσ], K= (3-sin φ) / (3+sin φ) | |
抛物线型 | 2 | 0 | 1/a | d/a | (同上双曲线型行) | ||
椭圆型 | 2 | tg2ˉφ | -2(a-a1)tg2ˉφ | tgˉφ(2a-ˉctgˉφ)ˉc | (同上双曲线型行) | ||
注 ˉc=6ccosφ/[31/2(3-sinφ)],tgˉφ=6sinφ/[31/2(3-sinφ)] |
应当指出, 由徐干成、郑颖人提出的莫尔-库仑等面积圆条件是按偏平面上圆面积与莫尔-库仑不等边六角形面积相等得出的。它不但消除了莫尔-库仑条件的奇异性, 而且计算结果与莫尔-库仑条件计算结果十分接近, 因而完全可等代莫尔-库仑条件。平面应变下莫尔-库仑匹配条件实际上就是平面应变下的莫尔-库仑条件, 这种情况下, 六角形锥体已转化为圆锥体。系数α1, k是在关联流动法则与非关联流动法则 (设体变为零, 即膨胀角为φ/2) 下导出。关联流动法则下为内切圆锥。
文献
式中
式中 D=[4 (1-e2) cos (π/6+θσ) +5e2-4e]1/2; e=q1/qc;ql, qc为受压与受拉时的偏应力。
偏平面上的剪切屈服曲线, 一般取偏平面上破坏条件的形状函数g (θσ) 。式 (20) 中的g (θσ) 可写成
式中 rc, qm为三轴压缩π平面上的半径和q值;rσ, q为π平面上相应任一θσ的半径与q值。
形状函数的选择, 通常应根据试验结果来定, 应用真三轴试验就能确定g (θσ) 。此外, 它还应满足如下条件:
1) 剪切屈服曲线应是凸曲线, 即要求
2) g (30°) =1, rσ (30°) =rc, g (-30°) =K,
式中 rl为三轴拉伸时π平面上的半径。K=1, 最简单的凸曲线是圆;K=1/2时, 唯一不凹的曲线是直线。
3) 当θσ=30°和θσ=-30°时
满足式 (26) 是为了消除奇异性, 简化计算。
莫尔-库仑条件并不是最符合实际的, 随着真三轴土工试验仪器的出现, 可通过试验来获得偏平面上的剪切屈服曲线。Lade根据砂的真三轴试验得出Lade屈服条件
这两种屈服条件在π平面上都是不规则的形状, 近似为一曲边三角形, 如图5所示。这两种条件没有角点, 都是光滑曲线, 而且Lade屈服曲线外接莫尔-库仑破坏条件的3个内角顶点, 而Matsuoka-Nakai破坏曲线外接莫尔-库仑条件6个内角点。
图5 π平面上Lade, Matsuoka-Nakai和郑颖人-陈瑜瑶屈服曲线
Fig.5 Yield curve of Lade, Matsuoka-Nakai and Zheng-Chen in the π plane
清华大学根据砂的真三轴试验
通过对重庆红粘土的真三轴试验
式中K, n为系数, K=0.69, n=0.45。屈服曲线与Lade屈服曲线十分接近 (图5) 。当无试验数据时, K可采用三轴拉伸试验时偏平面上的半径rl与三轴压缩试验时偏平面上的半径rc之比, n=0.4~0.5内取值。
《5.3θσ方向与q方向的剪切屈服面》
5.3θσ方向与q方向的剪切屈服面
θσ方向的剪切屈服面以
文献
式中α为偏离角。
式 (32) 意味着屈服函数Φγ, Φq, Φθ都成比例, 由此可得
求得了θσ方向与q方向的剪切屈服面。
在实际计算中, 常用的近似方法是把α角视为零, 即略去θσ方向的塑性应变增量, 而增大q方向的塑性应变增量, 使
《6 广义塑性力学中的硬化定律》
6 广义塑性力学中的硬化定律
传统塑性力学的dλ与硬化参量的函数有关, 即
硬化函数或硬化模量A与硬化参量有关。建立A的表达式称为硬化定律。在等向强化情况下, 引用相容性条件或一致性条件, 可得
金属材料不产生体积应变, 因而可选用塑性剪应变γp, 塑性总应变εp或塑性功Wp作为硬化参量。广义塑性力学采用分量塑性势面与分量屈服面, 各屈服面都有各自的硬化参量, 可各自表征各分量的硬化历史。体积屈服面、剪切屈服面、q方向与θσ方向剪切屈服面都应采用各自硬化参量的硬化定律。
1) ε
其矩阵形式为
式中 δT=
2)
广义塑性力学中, 如
3)
4) εpi硬化定律 同理, 广义塑性力学中, 如Φ=-εpi, Q=σi, 则A=1;如Φ=-H (εpi) , Q=σi, 则A=∂H/∂εpi。
各种硬化参量的硬化定律见表2。
《7 广义塑性力学的应力应变关系》
7 广义塑性力学的应力应变关系
广义塑性力学中, 应力-应变关系较为复杂,
表2 采用各种硬化参量的硬化定律
Table 2 Hardening laws of different hardening parameters
《表2》
H | A一般 | A特殊 (k, i=1, 2, 3) |
塑性功Wp | -∂Φ∂Wpσij∂Q∂σij | |
塑性应变εpij | -∂Φ∂εpij∂Q∂σij | |
塑性主应变εpi | -∂Φ∂εpi∂Q∂σi | 1 (Φk=-εpi, Qk=σi) |
塑性体应变εpv | -∂Φ∂εpv∂Q∂p | 1 (Φ=-εpv, Q=p) |
q方向塑性剪应变ˉγpq | -∂Φ∂ˉγpq∂Q∂q | 1 (Φ=-ˉγpq,Q=q) |
θσ方向塑性剪应变ˉγpθ | -∂Φ∂γpθ∂Q∂θσ | 1 (Φ=-ˉγpθ,Q=θσ) |
因为3个屈服面可能全部处于加载状态, 即全部屈服, 也可能只有其中一、二个屈服面处于加载状态, 即部分屈服。当全部屈服时, 应采用全部屈服时的应力-应变关系;而当只有一个屈服面屈服时应采用单屈服面应力-应变关系;当2个屈服面屈服时采用双屈服面应力-应变关系, 单屈服面与双屈服面应力-应变关系是其特例。
《7.1 广义塑性力学中的弹塑性柔度矩阵》
7.1 广义塑性力学中的弹塑性柔度矩阵
弹塑性刚度矩阵Dep可由弹塑性柔度矩阵Cep求逆得到, 而求Cep的关键是求塑性柔度矩阵Cp。
1) 依据单屈服面模型中Cep推广求广义塑性力学中的Cep:
与单屈服面类似有
式中 Ak见表2。因此有
2) 先求主应力空间中的塑性柔度矩阵Ap, 然后通过转换矩阵求Cep
如果εpk=Fk (σi) , (k, i=1, 2, 3) , 则有
通过坐标变换, 可由Ap获得一般应力-应变关系中的塑性柔度矩阵Cp
由此有 Cep=Ce+Cp。
《7.2 广义塑性力学中Dep一般表达式》
7.2 广义塑性力学中Dep一般表达式
广义塑性力学中, 采用等值面硬化规律可直接导出Dep的一般表达式
式中
,
其中
。即有
单屈服面情况下, 式 (47) 即变化为传统塑性力学中的弹塑性矩阵表达式, 传统塑性力学是广义塑性力学的一个特例。
《8 广义塑性力学的加卸载准则》
8 广义塑性力学的加卸载准则 [14]
加卸载准则既可由加卸载定义直接作出判断, 同时还可由屈服面状态给出加卸载准则。经典塑性力学经常采用后者;广义塑性力学中采用加卸载定义来判断加卸载较为方便, 且一般分为应力型与应变型2种加卸载条件。塑性变形与应力之间没有一一对应关系, 所以应力型加卸载准则理论上存在不足, 因而采用应变型加卸载准则。
无论加载或卸载, 总应变ε是单调变化量, 硬化材料或软化材料加载时, 应变ε总是增大, 卸载时应变ε总是减少。采用应变作为硬化材料加卸载准则是一种较简单的形式 (其他形式略) , 以压缩型土的体变为例, 可写成
式中 dε
同理可用来分析剪切屈服情况。由于假设d
《9 考虑应力主轴旋转的广义塑性位势理论》
9 考虑应力主轴旋转的广义塑性位势理论 [26]
经典塑性力学中假设应力主轴不发生偏转, 即硬化过程中应力主轴方向不变。实际上, 由于应力导致的各向异性, 应力主轴方向会发生偏转, 并同时由此引起塑性变形。为了反映主应力轴旋转所产生的塑性变形, 国内外学者做了不少工作, 主要在两方面:直接建立一般应力增量分量与一般应变增量分量之间的关系模式, 不过这样做太复杂;采用运动硬化模型, 屈服面随主应力轴旋转而运动, 当应力路径复杂时, 很难给出屈服面运动规律。
通过应力增量的分解, 分解出共轴应力增量与旋转应力增量 (主轴旋转所产生的应力增量) 。建立包含旋转应力增量在内的广义塑性位势理论, 就能求出应力主轴旋转产生的塑性变形。
《9.1 应力增量的分解》
9.1 应力增量的分解
9.1.1 二维应力增量的分解
应力增量可分解为共轴部分分量dσc与旋转部分分量dσr, 即
式中 K1=dσ1, K3=dσ2, K2=dθ (σ1-σ2) , 变换矩阵
9.1.2 三维应力增量分解 若将二维应力增量分解推广到三维情况, 可得
式中 dθ1, dθ2, dθ3分别表示旋转应力增量dσr1, dσr2, dσr3引起的绕第三、第一、第二主应力轴旋转的旋转角增量;T为变换矩阵, 同式 (45) 。由式 (51) 可见, 应力主轴旋转角增量引起剪应力增量, 即
《9.2 考虑应力主轴旋转的广义塑性位势理论》
9.2 考虑应力主轴旋转的广义塑性位势理论
由土工试验可知, 在主应力和主应变空间内, 旋转应力增量dσr引起6个应变方向的塑性应变, 需引用6个塑性势函数。可以任意选择势函数, 但必须保持势函数的线性无关。一般可把6个应力分量写成6个势函数, 6个应力分量的方向就是6个势面的方向。结合式 (8) 可得考虑应力主轴旋转的广义塑性位势理论
式中 dε
《10 岩土屈服条件的确定》
10 岩土屈服条件的确定 [23]
岩土材料极其复杂, 不同的土性常有不同的屈服条件, 导致计算结果的不唯一性。笔者提出用试验、按科学方法来建立屈服条件, 尽量减少确定屈服条件的随意性。
《10.1 屈服条件 (加载条件) 的物理意义》
10.1 屈服条件 (加载条件) 的物理意义
弹性力学基本方程与塑性力学基本方程的差别在于应力-应变关系。表3中列出了不同力学状态下的应力-应变关系及其参数的影响因素, 可见, 弹性材料或是塑性材料, 弹性参数与屈服条件都是材料的状态参数, 这就是屈服条件的物理含义。
图8示出了线弹性、非线性弹性、经典塑性与广义塑性的应力-应变曲线。线弹性情况下是条直线, 非线性弹性情况下是条曲线, 由此即可确定E与Et;经典塑性情况下是一组曲线, 广义塑性情况下是几组曲线, 需要按屈服面定义将试验曲线转换为屈服曲线。
表3 不同力学状态下的应力-应变关系及其参数的影响因素
Table 3 Stress-strain relation and influencing factors of parameter under various mechanical states
《表3》
力学状态 | 应力-应变关系 | 力学参数 | 参数的影响因素 |
线弹性 | 单轴情况下 εi=σi/E | E (弹性模量) | 材性 |
非线性弹性 | 单轴情况下 εi=σi/Et | Et (切线弹性模量) | 材性与应力状态 |
经典塑性 | dεpij=1A∂Φ∂σijdσij∂Φ∂σij A=-∂Φ∂εpij∂Φ∂σij | Φ (加载面) | 材性、应力状态与应力历史 |
广义塑性 | 分量应力-应变关系 dεpij=3∑k=11Ak∂Φk∂σijdσij Ak=∂Φk∂εpij(k=1‚2‚3) | Φk (分量加载面) | 材性、应力状态与应力历史 |
图8 线弹性、非线弹性、传统弹塑性、 广义弹塑性应力应变关系
Fig.8 Stress-strain curves of linear elasticity, non-linear elasticity, classical plasticity and generalized plasticity
由图8可见, 弹性参数与屈服条件都应由试验确定。岩土的性质差异很大, 只有用当地土获得的屈服条件, 才能提高计算的准确度, 且较为经济。
《10.2 屈服条件及其参数的建立》
10.2 屈服条件及其参数的建立
10.2.1 由试验数据作屈服面
由室内试验可测出总应变, 将塑性应变从总应变中分离出来, 得塑性应变与应力的关系, 如图9所示。对塑性应变分量ε
通常, 图10所示的屈服曲线可写成:
式中 Hk (εpk) 是硬化函数。由于试验获得的屈服曲线并不完全符合等向强化, 因而宜采用式 (54) 。
10.2.2 剪切屈服面的拟合
1) 子午面上γ方向的剪切屈服曲线
常用的剪切屈服曲线是双曲线型和抛物线型。
a. 双曲线型 双曲线的方程为
根据试验, 针对不同的
则γ方向双曲线型剪切屈服面为
b. 抛物线型 其方程为
表4可见, a是
抛物线型剪切屈服面为
表4 a, b与
Table 4 Value of a, b and
《表4》
ˉγp | a | b | ˉγp | a | b |
1 | 1.899 | -0.000 506 | 7 | 1.115 | 0.000 089 3 |
2 | 1.689 | -0.000 43 | 8 | 1.058 | 0.000 135 |
3 | 1.439 | -0.000 104 | 9 | 1.021 | 0.000 181 9 |
4 | 1.293 | -0.000 105 | 10 | 0.969 | 0.000 357 |
5 | 1.17 | 0.000 169 9 | 11 | 0.944 | 0.000 347 8 |
6 | 1.123 | 0.000 141 | 12 | 0.914 | 0.000 574 |
c. 子午面上γ方向剪切屈服面的验证 将试验数据与模型曲线做了一系列对比, 由图11可见, 双曲线型屈服面是合理的。若这二种类型的曲线都拟合得不理想, 应另设曲线重新拟合。获得Φγ后, 由式 (34) 求得q方向剪切屈服曲线Φq。
2) 偏平面上的q方向和θσ方向剪切屈服面
要确定偏平面上q方向和θσ方向剪切屈服曲线, 需要进行真三轴试验。在一般情况下, θσ影响产生的塑性变形约10%, 可借鉴前人的试验结果。
10.2.3 体积屈服曲线
用等值面理论, 在p-q平面内作出等ε
1) 压缩型
这类体积屈服面是封闭形的, 一端与p轴相接, 另一端与极限线相接。其形状有子弹形、蛋形、椭圆形、直线形等。从重庆红粘土的常规试验的数据可以得出椭圆的拟合关系式为
对不同的ε
表5 a2, b2与ε
Table 5 Value of a2, b2 and ε
《表5》
εpv/% | a2 | b2 | εpv/% | a2 | b2 |
0.5 | 23 000 | 6 500 | 2 | 86 000 | 32 000 |
0.8 | 35 500 | 10 000 | 2.5 | 107 000 | 41 000 |
1.0 | 44 000 | 15 000 | 3 | 128 000 | 52 000 |
1.5 | 65 000 | 21 300 | 3.5 | 149 000 | 62 500 |
1.8 | 77 500 | 26 000 | 4 | 169 900 | 75 000 |
从a2, b2值中可拟合出a2, b2与ε
拟合的屈服条件与试验数据对比, 从图12可见, 拟合的椭圆形屈服条件与试验数据比较吻合。
图12 压缩型土体 (重庆红粘土) 的椭圆型体 积屈服条件与试验数据的验证
Fig.12 Proof of compressive type volumetric yield condition based on test data
2) 压缩剪胀型
段建立、郑颖人、陈瑜瑶提出的压缩剪胀型土体的屈服面是一个S型曲线, 在状态变化线的上部分近似为直线, 下部分近似为椭圆。因此用分段曲线来拟合这类体积屈服条件。
由福建标准砂的试验数据, 对不同的ε
从表6可以看出, a1随ε
表6 直线段a1, b1与ε
Table 6 Value of a1, b1 and ε
《表6》
εpv /% | a1 | b1 | εpv /% | a1 | b1 |
0.03 | 1.38 | -3.0818 | 0.20 | 1.38 | -26.560 |
0.05 | 1.38 | -5.843 | 0.25 | 1.41 | -33.480 |
0.10 | 1.40 | -12.746 | 0.3 | 1.40 | -40.350 |
0.12 | 1.42 | -15.507 | 0.35 | 1.37 | -47.260 |
0.15 | 1.39 | -19.650 | 0.40 | 1.40 | -54.180 |
表7 曲线段a2, b2与ε
Table 7 Value of a2, b2 and ε
《表7》
εpv /% | a22 | b22 | εpv /% | a22 | b22 |
0.03 | 2 000 | -1 700 | 0.20 | 4 000 | 2 500 |
0.05 | 1 700 | -1 400 | 0.25 | 6 500 | 5 000 |
0.10 | 1 500 | -600 | 0.30 | 10 200 | 7 800 |
0.12 | 1 700 | -200 | 0.35 | 15 000 | 11 200 |
0.15 | 2 200 | 700 | 0.40 | 20 100 | 15 000 |
将拟合的屈服条件与试验数据对比, 见图13, 拟合的屈服条件与试验曲线比较吻合。
《11 广义塑性力学的实际应用》
11 广义塑性力学的实际应用
《11.1 岩土建模方面的应用》
11.1 岩土建模方面的应用
以广义塑性力学为理论基础, 建立符合土体实
图13 压缩剪胀型土体 (福建标准砂) 的体积 屈服条件与试验数据的验证
Fig.13 Proof of hardening compressive dilative type volumetric yield condition based on test data
际变形特性的本构模型是广义塑性力学最主要的应用之一。文献
按广义塑性力学, 分别建立3个分量塑性势面与相应的3个屈服面, 构成了土体本构模型。如果略去洛德角方向的塑性应变增量分量, 就将上述三屈服面模型简化为双屈服面模型。若进一步略去塑性体应变, 就可简化为单屈服面模型。为了便于分析不同屈服面对土体变形的影响程度, 表8给出了几种基于广义塑性力学的土体模型的计算结果。
从表8可见, 在均布荷载作用下地基表面处的沉降变形:q方向上剪切变形引起的沉降占总沉降中的64.5%;θσ方向上剪切变形引起的沉降占总沉降的16.1%;总剪切变形引起的降沉占总沉降的80.6%;体积变形引起的沉降占总沉降的19.4%。采用剪切单屈服面模型时, 变形计算误差可达35.5%, 当不计θσ影响时计算误差可增至38.7%, 采用双屈服面模型时, 变形计算误差可达16%, 不计θσ影响时计算误差可增至19.4%。
表8 几种基于广义塑性力学的土体模型及其计算所得的地基表面沉降量
Table 8 Several soil models based on generalized plastic mechanics and the calculated settlements of foundation
《表8》
方案 | 屈服函数与塑性势函数 | 均载作用下/m | 所有荷载作用下/m | ||||
x=0.5 | x=2.5 | x=6.0 | x=11 | x=17 | |||
1 | f=fq (p, q) ; Q=q | 0.19 | 0.94 | 0.73 | 0.35 | 0.22 | 0.17 |
2 | f=fq (p, q, θ) ; Q=q | 0.20 | 1.11 | 0.86 | 0.37 | 0.22 | 0.17 |
3 | f1=fv (p, q) , f2=fq (p, q) ; Q1=p, Q2=q | 0.25 | 1.13 | 0.90 | 0.45 | 0.29 | 0.23 |
4 | f1=fv (p, q) , f2=fq (p, q, θσ) ; Q1=p, Q2=q | 0.26 | 1.23 | 0.93 | 0.46 | 0.29 | 0.23 |
5 | f1=fv (p, q) , f2=fq (p, q, θσ) , f3=fθ (q, θσ) ; Q1=p, Q2=q, Q3=θσ | 0.31 | 1.30 | 1.03 | 0.55 | 0.37 | 0.31 |
当偏载与均载共同作用时, 在最大变形处的沉降变形:q方向剪切变形引起的沉降占总沉降中的80.5%;θσ方向上剪切变形引起的沉降占总沉降的5.4%;总剪切变形引起的沉降占总沉降的85.9%;体积变形引起的沉降占总沉降的14.1%。采用剪切单屈服面模型时, 变形计算误差可达19.5%, 当不计θσ影响时误差可增至27.7%。采用双屈服面模型时, 变形计算误差可达5.4%;当不计θσ影响时, 变形计算误差可增至13.1%。
《11.2 极限分析方面的应用》
11.2 极限分析方面的应用
当前的滑移线理论中都采用经典塑性理论中的关联流动法则, 对于金属材料, 其塑性势面 (速度滑移线) 与屈服面一致, 都是圆形曲线。金属材料的破坏面 (应力特征线) 也是基于米赛斯屈服条件, 而获得的圆形曲线, 因而应力特征线与速度滑移线一致, 如图14a所示, 它们与x轴的夹角为π/4。对于岩土材料, 一般也采用关联流动法则, 假设塑性势面与莫尔-库仑屈服面一致。按前人研究, 此时应力特征线与速度滑移线夹角成φ角。而试验得知, 岩土材料并不服从关联流动法则。广义塑性力学的出现, 从理论上证明了塑性势面与莫尔-库仑屈服面之间成一定的角度, 因而应按非关联流动法则来研究速度滑移线。文献
基于传统塑性理论的极限分析方法, 广泛应用于金属材料并获得了成功, W. F. Chen
广义塑性力学还会对岩土动力模型、运动强化模型及各种新的岩土本构模型产生影响, 有些已经被认可的方法也会受到影响而需进行某些改造。
《12 结论》
12 结论
1) 广义塑性力学消除了经典塑性力学中的传统塑性势假设、正交流动法则假设与不考虑应力主轴旋转的假设, 从固体力学原理直接导出了广义塑性位势理论。
2) 广义塑性力学是基于分量塑性势面与分量屈服面的理论, 能反映应力路径转折的影响, 即应力增量对塑性应变增量的影响。
3) 广义塑性力学中的塑性势面是已知的, 因而它不会产生当前非关联流动法则中任意假定塑性势面引起的误差。
4) 广义塑性力学中要求屈服面与塑性势面对应, 而不要求相等, 避免了采用正交流动法则引起过大剪胀等不合理现象。由于它对屈服面硬化参量的选定有严格的规定, 保证了岩土材料在一定应力路径下求解的唯一性。
5) 广义塑性力学中, 按土性及其状态不同, 体积屈服面可分为压缩型、硬化压缩剪胀型与软化压缩剪胀型三类, 并依据试验首次提出了压缩剪胀型土体的体积屈服面, 可以科学地考虑土体的压缩与剪胀。剪切屈服面分为q方向及θσ方向的剪切屈服面, 一般情况下可略去θσ方向的剪切屈服面而只考虑q方向的剪切屈服面。
6) 广义塑性力学采用分量塑性势面与分量屈服面, 各屈服面都有各自与塑性势面相应的硬化参量。分量塑性势面被确定, 分量屈服面就唯一确定。给出了广义塑性力学的硬化定律和应力-应变关系。
7) 在应力增量分解的基础上, 建立了考虑应力主轴旋转的广义塑性位势理论, 从而可求出应力主轴旋转产生的塑性变形。
8) 分析屈服面的物理意义表明:屈服条件是状态参数, 它与应力状态、应力历史及材性等状态量有关;它也是试验参数, 只能由试验给出。
9) 广义塑性力学不仅可以作为岩土材料的建模理论, 而且还可以应用于诸如极限分析等土力学的诸多领域, 具有广阔的应用前景。
广义塑性力学为岩土塑性力学打下良好的理论基础, 但仍处于起步阶段, 研究范围只限于静力模型与应变硬化阶段, 尚有大量的后续工作要做。