岩石应变软化弹塑性有限元分析

D0I:10.13374/i.issn1001053x.1989.01.001 第11卷第4期 北京科技大学学报 Vol.11 No.4 1989年7月 Journal of University of Science and Technology Beijing July 1989 岩石应变软化弹塑性有限元分析 吴玉忠 (矿业研究所) 摘要：本文从岩石的变形特性和界面约束理论出发，建立了应变软化规律及其有限 元本构方程。这种方程的求解存在着解法难与收敛慢和漂移界面的两个问恶。解快方法： (1)分析各种解法，采用初应力法，并进行改进称之改进的初应力法，以加速其收敛，(2) 考查各种源移界面的修正方法，选取其中较好的一种。编制了一个相应的EPSF限序，并对 金川某巷道进行了软化有限元分析。 关键词：应变软化，初应力法，界面源移 Finite Element Analysis of Elastoplastic Rock with Strain-softening Wu Yuzhong ABSTRACT:Based upon rock deformation characteristics and bounding surface model,a strain-softening law and its constitutive equation are establi- shed.In the equation,two difficult problems on solution and surface drifting are existed.First,after numerical approaches to softening problem were evaluated,an initial stress method is employed and then improved.It is called improved initial stress method which accelerate convergence.Secondly, various surface drifting updatings are compared and a good correction is sele- cted.Correspondingly,the program called EPSF is developed.Finally a strain-softening analysis to a tunnel in Jinchuan is given. KEY WORDS:strain-softening,initial stress method,surface drifting 在地下巷道、矿岩崩落、屈服旷柱、边坡等工程分析中，应变软化引起的渐进破环是一 个极为重要的问题。早在1936年人们就已经提出了边坡渐进破坏问题，当时Terzagh假定为 1988一08一09收稿 293

第 卷第 期 北 京 科 。 气 兮 年 月 技 大 学 学 报 岩石应变软化弹塑性有限元分析 吴 玉 忠 矿业 研究 所 、 摘 要 本文从 岩石 的 变 形特性和界面 约 束理论 出 发 ， 建 立 了应 变软 化规 律及 其有限 元本构 方程 。 这 种 方程的求 解存在着 解法难与收 敛慢和漂移界面 的 两 个问题 。 解决方法 分析各种 解法 ， 采 用 初应力法 ， 并 进行改进称之改 进的初应力法 ， 以加速其收 敛 考查各种漂 移界面 的修正 方法 ， 选取其 中较好的 一种 。 编 制 了一 个相 应 的 程 序 ， 并 对 金川某巷道 进 行 了软化有限元分 析 。 关键词 应 变 软化 ， 初应 力法 ， 界面 漂移 一 、 叹 平“ 之 戈 ， 一 。 ， 。 ， ， · 手， 。 ， · 一 · 一 ， ， 一 在地下巷道 、 矿岩 崩落 、 屈 服矿柱 、 边 坡等工程分 析 中 ， 应 变软化引 起 的渐进破坏是一 个极 为重要 的 问题 。 一 早在 年人们就 已经提 出 了边坡渐进破坏 问题 ， 当时 假 定为 一 一 收稿 DOI ：10．13374／j ．issn1001－053x．1989．04．001

一种软化机理来解释导致这种含行裂隙的坚硬岩峰值强度随时间降低的现象。随后就现 了各种力图解释这种渐进破坏的机即【1」。后来，Ladanyi引入应变软化概念分析地基承载 能力时，指出它的重要性，同时不少学者也强调这种渐进破坏过程需要进行有限元分析〔1。 1972年Hog首次把增量弹塑性有限元应用于模拟地基岩层在不排水下的应变软化特性： Nayak等【2】把软化概念引入本构关系进行平面应变冲压问题的有限元分析等等。在理论 上，Prevost和Hoeg证明.了由弹性能限制软化率的相对大小以维持相关流动的应变软化的唯 一性；Maier等进一步讨论了软化材料的稳定性和解的存在与唯一性，给软化有限元分析提 供了理论依据。此后软化问题就得到了普遍重视〔3,4)。 木文首先从岩石的变形特性建立应变软化规律及其本构方程；鉴于这种方程求解上的闲 难，分析各种软化有限元解法：选取一种初应力法，并进一步改进：考察应力收敛于界面的 状况及其校正方法，进而编入程序并对金川某巷道进行软化有限元计算分析。 1应变软化有限元本构方程 岩石的破坏是一个渐进的发展过程。岩石受压到峰值时开始屈服，产生塑性变形，随着变 形的发展，强度不断降低，出现应变软化，对应有塑性体积膨胀，最终达到我余强度。这一 试验的变形特性同临界状态理论是一致的，即上约束界面或称上屈服面反映塑性应变软化， 并以临界状态面为限【1)。 应变软化的上约束界面的形式及其收敛过程且前众说纷纭，但可以采用能反映塑性软 化的塑性准则进行拓展。它的上界面形状有： 椭圆界面，部分椭圆面和部分双曲面：直线族 面等。其界面收敛过程的方法也有多种，如以 常偏心率的椭圆面收缩等〔1)。 在相关流动法则下，这些收缩过程均不能 很好地反映岩石的软化规律，不能满足扩容率 0 随变形而减小的变形过程。因此笔者建议采用 图1建议的应变软化示意用 图1的软化界面收缩过程，并采取如下界面函 Fig .1 Schematic of deformation process of strain-softening model 数1： proposed F=/J2-a(h)(In(h2)-I.)=0 (1) 在变形中，a(h,)是改变体积的参量，I,(h2)是j偏应力有的参量。基于Drucker-Prager 准则和膨胀率减小过程，取a(h:)和I。(h2)为： a(h:)=a。-（a,-a)h,1(h1+B,) (2a) In(h2)=3.ctgoC+(C.-C)h2/(h2+B2) (2b) 式中，g。=0.5(g,+g,),a:=(12+9cte'g,)2(任=p,r),g为内摩擦角，C为内聚力， P、脚标分别表示峰值状态和残余状态，B,、B:均为实验参数。 294

一种软 化机 理 来解释 导致 这 种 含 了裂隙 的坚硬 兴石峰值 强 度随 时 间降 低 的现 象 。 随 后就 出现 各种 力 图解释 这种 渐进 破坏 的机 理 〔 ’ 」 。 后来 ， 引 入 应 变软 化概念分 析 地 基承载 能 力时 ， 指 出它 的 重 要性 同时 不少 学者也 强 调这 种 渐 进 破 坏过程 需 要进行有限 元 分 析 〔 ‘ ’ 。 年 成 首 次 把 增 量 弹 塑 性 有限 元应 用于模 拟 地 基岩层 在 不排水下 的应 变软 化特 性 等 〔 ” 把 软 化 概 念弓 入本构关系进 行 平面 应 变冲压 问题 的有限 元分析等等 。 在理 论 上 ， 和 。 电证明 了 由弹性能限制软化 率前相对大小 以维持 相 关流动 的应变 软化 的唯 一性 等进 一步 讨 论 一 软 化材料 的稳定性 和解 的 存在 与唯一性 ， 给软 化 有限 元 分 析提 供 了理 论依据 。 此 后软 化问题就得到 了普遍 重视 〔 “ ， ’ 。 本文首 先从岩 石 的变形特性建 立应变软 化规 律及 其 本构 方程 鉴 于 这种 方程求解上 的困 难 ， 分析各种软化 有限 元解 法 选 取 一种初应 力法 ， 并进 一步 改进 考察应 力收敛 干界 面 的 状况 及其校 正方法 ， 进 而 编 入程 序并 对金 川某巷道 进 行软 化 有限元 计算 分析 。 应变软化有限元本构方程 岩石 的破坏是 一 个渐进 的发 展过程 。 岩 石受压 到峰 值 时开 始屈服 ， 产 生塑性 变形 随 着变 形 的发展 ， 强度不断 降低 ， 出现 应变软化 ， 对 应有塑 性 体积 膨胀 ， 最终达到残余强 度 。 这 一 试验 的变形 特性 同临界 · 状 态 理 论是一致 的 ， 即上 约束界面或 称上 屈服面 反映塑性应 变软化 ， 并 以 临界状 态面 为 限 〔 ‘ ’ 。 应变软 化 的上 约束界面 的形 式及 其收敛 过 程 目前 众 说纷 纭 ， 但可以 采用 能反 映 塑 性 软 化 的塑性准则进 行拓展 。 它 的上界 面形状有 椭 圆界面 ， 部分椭圆面 和部分双 曲面 直线族 面 等 。 其界 面收敛过程 的方法也 有多种 ， 如以 常偏 心 率的椭圆面收缩等 〔 ‘ ’ 。 在相 关 流动法则下 ， 这 些收缩过程 均不能 很 好地反映岩石 的软化规律 ， 不 能满足 扩容 率 随 变形而减 小 的变形过程 。 因此 笔 者建 议采 用 图 的软 化界面 收 缩 过程 ， 并采取如下界 面 函 数 「‘ ’ 沐吕 兀试 犷 呻 尸产 图 建 议的 应 变软化 示 意图 一 厂 ， 一 了 一 。 一 在变形 中 ， ， 是改 变体积 的 参量 ， 。 是 ，歹偏 应 力 有 又 参量 。 基于 一 叱 准 则 和膨胀 率减 小过程 ， 取 ， 和 。 为 ， 一 ， 一 ， · 。 ， 二 · 只切 ， 〔 ， ， 一 ， 〕 式 中 ， 、 。 、 ， 、 ， 。 ‘ 。 饭 、 ‘ 一 考 ‘ ， ， 为 内摩擦 角 ， 为内聚 力 ， 、 ，脚标 分别表 示峰值状 态和残 余状 态 ， 、 均 为实验 参数

利用小变形分解条件、线性相关流动和-致条件，与(1)式联立，并注意广义Hook©定 律，解得正常应变软化弹塑性有限元本构方程为： d{n}=〔Dep〕d{e} (8) 其中： D=0-D=0)-AA.07 ·{6,}'D (4) A=月}'o· (5) A=-3·a(h:).〔(1，(h)-1)(a,-a,)(B1+h) +(6h)+）·89a.〕 J2 (6) 6-{a0是，a+是. I 2√WJ3 、y-a(h:)+ 2,v万，1v,t.万}（) 3 dh,三de=3a(h)d h=fde=e+C:(C=0) dh:=de=（3a(h)+t+s）d h2=∫ae=e+C2(当d{e)=0附，h2=0即Ca=0) 以上各式中，〔D门、〔D门、〔D门分别为弹塑性、弹性、塑性矩阵，A、A。分别为软化参量 和临界软化参量，且A。<A。d>0,且与A有关的常量。e为体积应变，e”为等效塑性应 变。取B1=B2=0.0351)。h1、h2皆为软化材料内变量。 2应变软化有限元解法 有限元的非线性问题，包括几何的和物理的两种，本构关系是属于物理非线性的，其解 法很多。按对节点不平衡力的处理方式可分为：精确解法、增量法（或称初值法）和自校正 法3大类。它们相互结合，就产生了各式各样的混合法。决定采用的解法都与各种因素有 关，如本构关系、解的精度、收敛快慢等。本文则着重考虑与本构关系有关的因素。 295

利 用 小 变 形 分解 条件 、 线 性相关 流动 和 一致 条件 ， 与 式 联立 ， 并 注意广 义 定 律 ， 解得 正常 应变软 化弹塑性 有限 元本构 方程 为 汀 〔 ’ 〕 一 。 其 中 〔 ‘ ” 〕 〔 〕 一 〔 夕 〕 〔 〕 通 一 月 〔 〕 刁尸 日 》 ， 一 下 ’ 不 一丁 全 仁， 至 全 一乃 一 ‘ 少 麦口 圣 、 〔 〕 刁 日 仃 刁 一 · · 〔 。 石 一 · ， 一 户 · 人 十 二 十 ， 二 一 二 ， 一 ， 侧 。 甲狱。 〕 乙 口 一 月 口 ， 人 仃 二 一 了 侧 八 了 几 ， － ， 二二二二二二二 了 几 ， 厂丫 ， 了 几 ，几尸 “ “ 丫 二 。 厂 · 只 碑 了 。 里 ￡ 七 “ ￡ 。 一 产、 户 丽 二 · 。 ‘ ， 二 十 二 十 了二 二 丫 二 蕊 ， 子 当 ‘ “ ” 时 ， ‘ 即 一 。 ’ 以上各式 中 ， 〔 ’ 勺 、 〔 〕 、 〔刀勺分别为弹塑性 、 弹性 、 塑性矩 阵 ， 、 。 分别 为软 化参量 和 临界软化参量 ， 且 。 。 。 ， 且与 有关 的常量 。 此为 体积应 变 ， 户 为 等 效 塑 性 应 变 。 取刀 二 〔 ’ 〕 。 、 几 皆为软化材料内变量 。 应变软化有限元解法 有限元 的非线性问题 ， 包括 几何 的和物理 的两种 ， 本构关 系是属于物理非线性 的 ， 其解 法很多 。 按对节点 不 平衡力 的处理方式可分为 精确解 法 、 增量法 或称初值法 和 自校正 法 大类 。 它们相互结 合 ， 就产生 了各 式各样 的混合法 。 决定采 用 的解 法都 与 各 种 因 素有 关 ， 如本构关 系 、 解 的精度 、 收敛快慢等 。 本 文则 着重 考虑 与本构关 系 有关 的 因素

对软化有限元解法，Iog最早采用·种'与月前不同的切线刚度法：八ayak等()则采 州初应力法；Lo和LCc【5:也采川初应力法，但前者有别。在后所来采川的解法也多了，迄 今大体采用以下几种解法： （1)初应力法【1~4，]： (2)初应变法【1】 (3)切线刚度法〔3，)： (4)-阶白校正法1： (5)割线刚度迭代法1`。 此外，还有半解析解的迭代法和解析法。 这5种方法，无可厚此薄彼。理论上已证明(2)法在理想塑性条件下是发散的，且解软 化问题永远是发散的「？：甚至还有漂移结果的趋势(3。因此，Desi〔8指“采用与割线 刚度法联合的迭选代法，可成功地分析软化问题”并得到了应用【]。但这种方法的弊病是不 好处理初应力不为零的问题s,如地下工程中的原岩应力，更重要的是其解不收敛〔1)。对 (3)法(8)不但在塑性时1(1-2v)项趋向于无穷大，特别是存在着“因负的刚度在作用荷载 方向导致不真实的位移或应变”〔5)，但Hog的观点对此看法不同。无论如何，其刚度矩阵 在软化时是负定的。此外，人们还进行了增量法、迭代法和自校正法的比较后采用自校正 法，并进一步与初应变法和切线模量法比较，在计附上选用了一阶自校正法(3)。对(1)法， 尽管在理想塑性条件下收敛极慢，但最近的分析和实践表明：它的使用最广、解决软化问题 不会出现上述困难〔，4。作某些修正的初应力法收敛比(4)法更好1门。因此，本文软化问 题采用初应力法求解。 3改进的初应力法 增量初应力法不但可以克服上述问题，而且极易扩成：相关流动：：)和弹塑性辄合问题。 但它的收敛速度很慢。在理想塑性条件下的屈服单元较多附，迭代4次即收敛〔8】是不可 能的；计时上也不只是弹性：计算的10倍。因此，必须进行加速。月前加速的途径有两种： (1)借用变刚度修正：(2)对迭代过程进行加速，如单松地因广如速，引人·个超松弛对角 阵进行加速，还有Thomas()提出自校正单松弛因下加速。木文先采用变刚度法进行修正。 3,1增量初应力法[9 基予能量泛函、最小势能原理材料的儿何关系'与本构关系，可其弹塑性：行限心的平 衡方程为：〔Kg门{A8}”{AR}+{1R"}。在小变形条件下，第n级荷线的第次迭代 过程的方程为： 〔K。〕{A8}g={△R}。+{ARn}:-1 (8) 其中，〔K。)为常刚度矩阵，{△}为位移增t列阵，{△R}为荷载增量列阵，{AR}为 初应力{△σ°}引起的荷载增量列阵。 3,2改进的增量初应力法 改进的思路是：每次增量荷载之初采川]常刚度法，在每个增量小的送代过程采州变刚度 法进行荷载修正。弹塑性增量变刚度行限心平衡方程为（〔K,〕为变刚度矩阵）： 296

对于软 化 有限 元解 法 ， 。 己 最 旱 不川 · 种 ’ ， 目前 下同的切 线 刚度法 ， 等 〔 ’ 〕 则 采 用初 应力法 。 和 。 。 〔 石 也 采川初应 力法 ， 但 ‘ 前 择了湍 。 在 后来 采川 的解祛 也多 厂 ， 迄 今大 体采 用以 下几种解 法 初应 力祛 〔 ’ 一 ‘ ， “ ’ 初 应 变 法 〔 ” 切线 刚度法 〔 “ ， “ ’ 一阶 自校 正 法 广 “ 〕 割 线 刚度迭 代 法 〔 ” 。 此 外 ， 还 有半 解 析解 的迭 代法 和解 析 法 。 这 种方法 ， 无 可厚此 薄彼 。 理 论上 已证 明 法在 理 想 塑 性 条件下 是 发 散 的 ， 且解软 化 问题永远是 发散 的 〔 〕 甚 至还 有 漂移结 果的趋势 〔 ’ 。 因此 ， 。 ， 〔 ’ 指 出 “ 采用 一 与割 线 刚度法联 合的迭 代 法 ， 可成功 地 分析软 化 问题 ” 并 得到 了应 用 〔 ” ’ 。 但 这 种 方法 的弊病是 不 好 处理初应 力不为 零的问题 〔 了， 如地 下工 程 中的原岩应 力 ， 更 币 要 的是 其解 不收敛 〔 ’ 〕 。 对 法 〔 日 〕 不但在 塑性 时到 一 项趋 向 干无 穷大 ， 特别 是 存 在着 “ 因 负的刚度在作用 荷 载 方 向导致不真实的 位移或 应变 ” 〔 〕 ， 但 。 电 的观 点 对此 看法 不 同 。 无 论如何 ， 其刚度矩 阵 在软化时是 负定 的 。 此 外 ， 人们还 进 行 了增量 法 、 迭 代法 和 自校 正 法 的 比较后 采用 自校 正 法 ， 并进 一步 与初应变法 和切 线模 量 法 比 较 ， 在计时 卜选 用 了 一 阶 自校 正法 〔 “ ’ 。 对 法 ， 尽 管在理 想 塑性 条 件下 收敛 极慢 ， 但 最近 的分析 和 实践表 明 它 的使 用 最广 ， 解 决软化 问题 不会 出现上述 困难 〔 ” 碑 ’ 。 作某咋修正 的初 应 力法 收敛 比 法 更 好 〔 ” 。 因此 ， 本 文软 化 题 采 用初 应 力法求解 。 改 进 的初应 力法 增量 初应 力法 不但可以 克服上 述 问题 ， 而 且极 易扩成 】卜相 关 流动 「’ 〕 和 弹 塑性 祸 合问题 。 但它 的收敛 速 度很慢 。 在 理想塑 性 条件下 的 屈 服单 元较多 时 ， 迭 代 次 即 可收 敛 〔 〕 是 不可 能 的 计时上 也 不只是弹 性 计算 的 倍 。 因此 ， 必 须进 行 加速 。 目前 加速 的途 径 有 两 种 借 用变刚 度修 正 对迭 代过程 进 行 加速 ， 如单 松弛 因 一 加速 ， 引 入 · 个超 松 弛对 角 阵进 行加速 ， 还 有 、 〔 ‘ 〕 提 出 自校 正单 松弛 因 加速 。 木 文先 采 用 变 刚 度法 进 行修正 。 增 初 应 力法 〔 “ 墓 二能量 泛 函 、 最 小势 能原 理 和材 料 的儿何 关系 ’歹本构 关系 ， 衡 方程 为 〔 〕 入八 二 八 一 、 ’〕 。 在 小变形 条什 卜 ， 过程 的方程 为 才， 七弹 塑性 有限 元 的 、 几 第 刀 级 荷 载 的 第。 次迭 代 〔 。 〕 一 入 二 八 其 中 ， 〔 。 〕为 常 刚 度矩 阵 ， 八 玛 为 位移 增 欣列阵 ， 初 应 力 王△ 。 引 起 的荷载 增 笼 列阵 。 “ 八尸 。 丁 一 ’ 态 为 荷载 增 量 列阵 ， 入刀 们 为 改进 的增 初 应 力 法 改进 的思路是 睡次增 量荷 载之 初采 川常 刚度法 ， 在 侮个增 量 卜的迭 代过程 采 川 变刚 度 法进 行荷载修 正 。 弹 塑 性增 最变 刚 度有限 ，口卜衡 方程 为 〔 〕 为变刚度矩 阵

〔Kr]m-1·{△δ}0={AR}x (9) 分解K,:〔K=∫∫∫cBr(D-D)c8〔K-K,,代人(9)式立得 (CB)为几何矩阵)： 〔Ka)·{A8}={AR}.+〔K,-1·{△δ}？ 若迭代过程解唯一，则{△8}⊙{△6}：1，即上式为： 〔K{△}m={JR}n+K,-1.{Aδ}- (10) 联立(8)(10)两式，并令{△R′}-1=〔K-1.{18}”-1,得： 〔KJ〔△]={△R}n+{△Ra}m-1+{AR'}-1 (11) 这种考虑刚度修正项的方法称为“改进的初应力法”。改进后既保持了可以进行软化问 题求解的常刚度的优点：又兼有变刚度法收敛快的优点。 4应力漂离界面的校正方法 如上所述，方程(3)的非线性软化问题的有限元解法很多。然而，由于〔D〕需要多 次小步长迭代来求之，在每步求解中任何解法都假定了该步与{σ}有关的矩阵〔D·〕为一 常数。显然这是近似的，即使迭代步长再小，每个荷载增量后所求得的应力实际上不在该状 态的界面F。上，如图2。为此，作者进行过界 面追踪，结果发现10F>√J4)的单元并 不少见。对软化问题，这种偏差尤为显著。随 着荷载增量的增加，这种偏差就愈来愈大，使 结果漂移正确解。有限元法中，处于这类的单元 称之为过渡单元。在理想塑性条件下，常用的线 -1 图2漂离界面示意图 性插值方法（如方程(13）)使应力状态问到该 Fig.2 Sketch of surface drifting 状态的界面上。对于软化问题，{σ}不但与 〔D门有关，而且还与h有关，因此问题就更为复杂，必须认真对待。 对此，校正的方法有几种：Sloan和Randolph建议I,不变J:减小，使{o满足 F({σ}，h)=0:Owen等门提出改变累积的有效应力方法；Potts等〔10]用沿总的应变 增量方向的修正法；Nayak等[2)采用沿塑性应变增量方向的修正法；Vermeer也提出了“修 改”的校正方法；Thomas'4)则用试位迭代法进行线性插值修正。Potts等对前5种与用微 小的荷载增量多次迭代的标准结果比较表明：前3种误差甚大，不宜采用1：。作者针对后 3种作了进一步的对比后，选用了Nayak等人采用的方法（设h为常数）： (1)令：F({0}、h)=F.<0 (12) F({G}+{△o},h)=F>0J 297

〔 丁 〕 刃 一 ‘ · 众 丁 八 气 分解 〔 〕 〔 〕 了 、 丁 〔 〕 丁 ‘〔 〕 一 〔 ’ 〕 ，〔 〕 “ 〔 。 〕 一 〔 ， 〕 ， 代 人 式 立 得 〔 〕 为几何矩 阵 〔 。 〕 · 人占 丁 八刀 二 〔 〕 叹 一 · 入 叹 若迭 代过程 解 唯一 ， 则 八 贾公 人占 毋 一 ， 即上 式为 区 。 〕 · 王△占 下 一立 。 ‘ 区 ， 〕 万 一 一 八占 丁 一 了 联立 两 式 ， 并 令 刀 ‘ 份 一 ’ 二 〔 ， 〕 丁 一 ‘ ， 入 公 一 ， 得 〔 。 〕 · 〔 舀〕 分 八刀 穿 一 这种考虑 刚度修正项 的方 法称 为 “ 改 进 的初应 力法 ” 。 改 进 后既保 持 了可以 进 行软化 间 题求 解 的常 刚度 的优 狱 又兼 有变 刚 度法 收敛 快 的优 点 。 应 力漂离界面 的校正 方法 如上 所述 ， 方程 的非线性软 化 问题 的有限元 解 法很 多 。 然而 ， 由于 〔 〕需 要 多 次小步长迭 代来求之 ， 在每步 求解 中任何解 法都假 定 了该步 与 有 关 的矩 阵 〔 ’ 夕 〕 为 一 常数 。 显 然 这是 近似 的 ， 即使迭 代步长再 小 ， 每 个 荷 载增 量 后所 求得 的应 力 实 际 不在 该状 态 的界 面 。 上 ， 如 图 。 为此 ， 作 者进 行过界 面追 踪 ， 结 果发现 了 〔 〕 的 单 元 并 不少见 。 对软 化 问题 ， 这 种偏差 尤为显 著 。 随 着 荷载增量 的增加 ， 这种偏 差就愈来愈大 ， 使 结果 漂移正 确解 。 有 限 元法 中 ， 处于 这 类 的单 元 称之 为过 渡单 元 。 在理 想 塑性 条件下 ， 常 用 的线 性插值 方法 如 方程 〔 使 应 力状 态 回到 该 状 态的界面上 。 对 于软化 问题 ， 笼 不但 与 之 一次 匕 亏 。 ’ 一 八 图 漂离 界面 示 意 图 ， 么 可 又 〔 ” 〕 有 关 ， 而 且还 与 有关 ， 因此 问题就更 为 复杂 ， 必须 认真 对待 。 对此 ， 校正 的 方法有 几 种 引 和 建 议 ， 不变 减 小 ， 使 满 足 ， 等 〕 提 出改 变累积 的有效 应 力方 法 等 〔 ’ ” 〕 用沿 总的应 变 增量方 向 的修 正法 ， 等 〔 “ 〕 采 用沿塑 性 应 变增量 方 向 的修正法 也提 出 了 “ 修 改 ” 的校正 方 法 〔 ‘ 〕 则 用试 位迭 代法进 行线 性插值修 正 。 等 对前 种与 用微 小 的荷 载增量多 次迭 代 的标准结 果 比 较表明 前 种 误差甚大 ， 不宜 采 用 〔 ’ 。 〕 。 作 者针 对后 种 作 了进一步 的对 比后 ， 选 用 了 手 · 等 人采 用 的方法 设 为常数 令 。 八 ，