标题

Size: px

Start display at page:

Download "标题"

折象
7 years ago
Views:

1 第卷第期年月哈尔滨工业大学学报 JOURNAL OF HARBIN INSTITUTE OF TECHNOLOGY Vol No Dec. 饱和粘土中球孔扩张问题弹塑性解析, 李镜培,, 唐剑华,, 张亚国, 钟光玉 (. 同济大学土木工程学院地下建筑与工程系, 9 上海 ;. 岩土及地下工程教育部重点实验室 ( 同济大学 ), 9 上海 ;. 上海南汇建工建设 ( 集团 ) 有限公司, 99 上海 ) 摘要 : 为了研究静力触探试验及沉桩扩孔等工程问题, 基于修正剑桥模型, 推导了不排水条件下球孔扩张问题的半解析解. 将扩张球孔周围土体分为临界状态区塑性区以及弹性区三个区域. 弹性区内, 利用弹性理论得到应力和孔隙水压力的解答 ; 临界状态区及塑性区内, 利用相关联的流动法则拉格朗日分析法建立了关于应力的一阶非线性常微分方程组, 以弹塑性界面处的应力分量作为初值, 求解微分方程组可得到应力和孔隙水压力的解答. 研究结果表明 : 各向同性超固结比对扩孔压力土体应力超孔隙水压力以及塑性区范围均具有显著影响, 且扩孔过程中土体剪切模量并非常量, 其随扩孔半径各向同性超固结比的变化而变化 ; 同时通过与已有解答进行比较, 对本文方法的可靠性进行了验证. 关键词 : 球孔扩张 ; 剪切模量 ; 修正剑桥模型 ; 各向同性超固结比中图分类号 : TU7 文献标志码 : A 文章编号 : 7-()-7-7 Elastic plastic solution of sphee cavity expansion in satuated clay LI Jingpei,,TANG Jianhua,,ZHANG Yaguo,,ZHONG Guangyu (. Depatment of Geotechnical Engineeing, Tongji Univesity, 9 Shanghai, China;.Key Laboatoy of Geotechnical and Undegound Engineeing( Tongji Univesity), Ministy of Education, 99 Shanghai, China;.Shanghai Nanhui Constuction Goup CO. LTD., 99 Shanghai, China) Abstact: An exact semi analytical solution in the undained cavity expansion can be obtained on the basis of the MCC model to eseach the cone penetation test and the pile diving. The field aound the cavity can be divided into thee zones: citical zone, plastic defomation zone and elastic defomation zone. In the elastic zone, an analytical solution fo the distibutions of stess and excess poe pessue is deduced accoding to the elastic theoy. In the citical and plastic zone, a set of fist ode nonlinea odinay diffeential equations concening stess can be obtained accoding to the associated flow ule and the lagangian analysis method. The stessss and poe poe pessue can be solved as an initial value poblem stating at the elastic plastic bounday. The esults show that the isotopic ove consolidation atio has a significant influence on the stesses and the excess poe pessue.the shea modulus vay significantly with the cavity adius and the isotopic ove consolidation in the couse of cavity expansion. Keywods: sphee cavity expansion; shea modulus; modified cambidge model; isotopic ove consolidation atio 球孔扩张理论在旁压试验静力触探沉桩及压力注浆桩基承载力等岩土工程问题中有着广泛应用 []. 然而, 由于采用的本构关系不同, 所得到的结果也各不相同.Vesic [] 考虑塑性区土体的收稿日期 : --. 基金项目 : 国家自然科学基金 (788); 浦东新区科技发展基金创新资金项目 (PKJ-C8). 作者简介 : 李镜培 (9 ), 男, 教授, 博士生导师. 通信作者 : 李镜培,lijp77@ tongji. edu. cn. 压缩性, 给出了土体服从 M-C 屈服准则状态下球 [] 孔扩张的近似解 ;Yu 等假定土体为线性理想弹塑性体, 采用非关联流动法则的 M-C 屈服准则, 在考虑土体剪胀性的条件下求解了大应变情况下 [] 球形孔的扩张 ;Banejee 等把理想的刚塑性模型应用于正常固结粘性土的扩孔问题 ; 然而理想的弹塑性模型与刚塑性模型均不能考虑土体应力历史的影响, 而剑桥模型则可以克服这些缺陷. [] Collins 等采用临界状态模型推导了不排水情

2 7 哈尔滨工业大学学报第卷况下球形孔扩张的大应变解析解, 分析了超固结 [] 比对扩孔压力以及对超孔压的影响 ;Cao 等采用修正剑桥模型对不排水状态下球孔扩张问题进行解析, 但其在求解过程中假设塑性区内的偏应力为极限状态的偏应力, 剪切模量为常量, 使得结 [7] 果存在一定误差. 胡伟等基于剑桥模型推导了球孔不排水扩张的解析解, 但其对偏应力做了线性插值的近似处理, 与土体塑性区内偏应力非线 [8] 性变化的特征不相吻合.Chen 等在不对偏应力和平均有效应力进行任何简化假设的情况下, 利用修正剑桥模型给出了圆柱孔扩张问题的精确半解析解 ; 然而, 研究沉桩扩孔, 静力触探等问题时, 往往认为桩端或探头周围土体呈球形扩张 [9], 故而基于平面应变假设的圆柱形扩孔理论并不能直接使用. 以往的球孔扩张理论假设在扩孔过程中, 剪切模量为常量, 基于此, 本文在以往研究基础上, 利用修正的剑桥模型理论, 对球孔扩张问题进行解析. 得到了饱和粘土中球形孔扩张后, 弹性区临界状态区塑性区内的应力及超孔隙水压力分布特征, 给出了球孔周围土体剪切模量的变化规律, 分析了各向同性超固结比对球孔周围应力以及超孔压的影响. 并将球孔扩张理论应用到静力触探实验中, 本文的研究结果可为静力触探等试验提供一定的理论基础. 力学模型及基本假设采用如图所示球孔扩张模型, 球孔初始孔径为 a, 初始孔隙水压力为 u, 初始平均总应力为 p, 初始平均有效应力为 p, 经过扩张后, 最终孔径为 a. 假设压应力和压应变为正, 任意一点平均总应力为 p, 有效平均总应力为 p, 偏应力为 q, 孔隙水压力为 u. f 为塑性区与临界状态区交界面的半径, p 为弹塑性区交界处的半径, 其初始半径为 p, 对于塑性区内任意一点 x, 其初始位置为 x. 对于 a f 区域内的土体, 此区域内的土体的应力状态均处于临界状态 CSL 线上, 土体达到此状态后, 应力状态不再变化, 该区域为临界状态区 ; 当 f p, 土体处于塑性状态, 服从修正的剑桥模型 ; 当 p, 土体处于弹性状态, 服从虎克定律. 假定土体饱和均质各向同性不可压缩, 将球孔扩张看成不排水过程. 根据弹塑性理论可得 p = (σ + σ θ ) /,q = σ - σ θ, () 球孔扩张前土体处于各向同性状态, 可得到初始径向应力 σ, 初始切向应力 σ θ : σ = σ θ = σ,p = σ, () σ 为初始应力, 球孔扩张过程中, 土体内部的任意一点都满足下列平衡方程 : σ + σ - σ θ =, () 若用有效应力表示, 则方程 () 可化为 σ + u + σ - σ θ =, () 式中 :σ σ θ 分别为径向切向总应力,σ σ θ 分别为径向切向有效应力,u 为孔隙水压力. σ 弹性区图塑性区临界状态区 σa a σ σ 弹性区应力分析 f x x 球孔扩张模型 p p 弹性区内, 采用小应变理论, 假设压应变为正, 因此径向应变增量与切向应变增量可由式 () 表示 : σ dε = - du, dε θ = dε ϕ = - du. () du 为径向位移增量, 由虎克定律可得弹性区的应力应变关系 : dε e = E dσ - ν E dσ θ, () dε e θ = - ν E dσ + - ν E dσ θ. (7) 式中 :ν 为泊松比,E 为弹性模量,σ σ θ 分别为径向和切向有效应力. 对于修正的剑桥模型, 泊松比 ν 为常数. G 为剪切模量, 可由比体积 υ 平均有效应力 p 表示 [] : G = ( - ν)υp ( + ν)κ, (8) 式中 κ 为 υ - ln p 平面上卸载 - 再加载曲线的斜率. 由于球孔扩张过程可看成不排水扩张, 故弹性区的体积变化为零, 因此有 : dυ = - κdp / p =, (9) 由式 (9) 可得, 弹性区内的有效应力 p 和比

3 第期李镜培, 等 : 饱和粘土中球孔扩张问题弹塑性解析 7 体积 υ 保持不变, 因此在弹性区内, 剪切模量 G 为常数. 由平衡方程应力应变关系可得弹性区内的应力和位移 [] : æ σ = σ + (σ p - σ ) p ö ç è ø, (a) σ θ = σ ϕ = σ - (σ - σ p ) æ p ö ç. (b) è ø 式中 : σ p 为弹塑性边界处的总径向应力,G 为初始剪切模量. 由式 () 可得弹性区的平均应力保持不变, 又因为弹性区的平均有效应力保持不变, 因此弹性区内的孔压也保持不变. 塑性区应力分析有效应力塑性区内, 土体单元服从修正剑桥模型, 其屈 [-] 服函数为 F(p,q,p c ) = q - M [p (p c - p )], () 式中 :M 为 p - q 平面中 CSL 线的斜率,p c 为各向同性状态条件下的屈服应力. 利用弹塑性边界处应力连续性条件可得 p p = p, () 由式 () () 可知, 在弹塑性边界处, 偏应力 [] 为 q p = Mp R -, () 式中 R 为各向同性超固结比,R = p c / p []. 效应力 : 变增量 : 联立式 () 式 () 可得弹塑性边界处的有 σ p = p + q p, () σ θp = p - q p. () 根据塑性理论的相关联流动法则可得塑性应 dε p = ψ F σ = ψ æ F p p σ + F q ö ç è q σ ø ψ ( p (M - η ) + q ), dε p = ψ F θ σ θ = ψ æ F p p σ θ + F q ö ç è q σ θ ø dε p ϕ = dεp θ. ψ ( p (M - η ) - q ), = = (a) (b) (c) 式中 : ψ 为塑性流动因子,η = q / p, 定义为应力比. 理论可知塑性区土体的应变服从大变形理论, 由塑性 dε θ = - d /, (7) 球孔周围土体任意一点的体应变为零, 即 dε v =, 因此有 dε = - dε θ = d /, (8) [8] 根据 chen 等的方法, 结合式 () (8) 可得到塑性区内关于应力的微分方程 : 其中 dσ - b d χ =, dσ θ - c d χ =. b = E - ν - ν + E( + ν)(λ - κ)(a θ - a a θ ) υp (M - η ) (9a) (9b) { }, c = E ν + ν - - E( + ν)(λ - κ)(a a θ + a ) υp (M - η ) (a) { }, χ = - + ν E (- + ν + ν ) + (b) E(- + ν)(λ - κ)a υp (M - η ) E(λ - κ)(a + νa θ a θ ), (c) υp (M - η ) a = p (M - η ) a θ = p (M - η ) + q, - q. - (d) (e) 式 (9) 是运用拉格朗日分析法建立的微分方程组, 该方程组适用于塑性区内任意一点 x, 若要通过式 (9) 求得 x 处的应力, 只需求得 x 处土体单元由弹性状态变为塑性状态的位置 xp, 以及在 xp 处的应力初始值. 下面论述中将给出应力初始值. xp 处的应力初始值与径向距离无关 [8], 可由式 () () 确定, 即 σ ( xp ) = σ p, σ θ ( xp ) = σ θp. (a) (b) 球孔扩张过程可看作不排水过程, 球孔周围的土体体积应变为零, 结合式 () 可得 xp 表达式 xp a = - σ - σ p G æ æ ç ç è è x a ö ø æ + a ö ç è a ø ö - ø, () 以弹塑性界面处一点为研究对象, 则 x xp 等于 p, 由式 () 可得 p a = a a æ æ ö ç ç - ç è ø ç æ - σ - σ p ö ç ç è è G ø ö - ø. ()

4 7 哈尔滨工业大学学报第卷因此, 塑性区内球孔扩张问题可归结为求解一系列具有初值条件的非线性常微分方程, 其中式 (9) 为控制方程, 式 () () 为初值条件 ; 该非线性微分方程组可通过数值方法求解. 超孔隙水压力扩孔过程中, 在弹性区内, 由于孔压保持不变, 不产生超孔压 ; 在塑性区内, 超孔压可通过对式 () 积分得到, 积分区间为 [ x p ], 则 p (σ Δu( x )= u( x ) - u = σ p - σ ( x ) + - σ θ ) d, x () 由于无法得到超孔隙水压力的解析解, 故可利用数值积分求解式 (). 应力应变关系与偏应力相对应, 偏应变可由下式确定 [] : ε q = ((ε - ε ) + (ε - ε ) + (ε - ε ) ), () 对于球孔扩张问题, 式 () 可简化为 ε q = (ε - ε θ ), () 又因为球孔扩张过程为不排水过程, 体积应变为零. ε v = ε + ε θ =, (7) 塑性区内, 采用对数应变 : 将式 (8) 代入式 (7) 可得将式 (9) (8) 代入式 () 可得 ε θ = ln, (8) ε = - ln, (9) ε q = ln. () 在弹性区内, 偏应力 q 可由式 () 确定, 偏应变可由式 () 确定, 因而联立式 () 式 () 可确定弹性区内的应力应变关系 ; 在塑性区内, 通过求出微分方程 (9) 的数值解后, 可得塑性区内土体的偏应力, 结合式 () 可得塑性区土体的应力应变关系. 算例分析选取文献 [8] 上的计算参数 : 初始有效应力 p = kpa, 临界状态线 CSL 的斜率 M =, 初始等向固结曲线 INCL 的斜率 λ =, 初始比体积 υ = 97, 泊松比 ν = 78, 初始孔压 u = kpa, 初始剪切模量 G = kpa, 回弹曲线的斜率 κ =, 不排水剪切强度 c u = kpa. 结合上述理论, 分析了各向同性超固结比对偏应力扩孔压力径向应力切向应力以及超孔压的影响, 并与文献 [] 的结果作对比. 图为扩孔半径 a / a =, 各向同性超固结比 R = 时, 偏应力 q 的变化规律, 随着各向同性超固结比的增大 R, 临界状态区的半径逐渐减小. 当 R <, 随着径向距离的减小, 偏应力 q 逐渐增大, 当土体达到临界状态时, 偏应力 q 保持不变, 土体在屈服之后表现硬化的性质 ; 当 R =, 土体一旦屈服, 偏应力 q 就保持不变, 土体在屈服之后表现理想弹塑性的性质 ; 当 R >, 随着径向距离的减小, 偏应力 q 先增大后减小, 当土体达到临界状态时, 偏应力 q 保持不变, 土体在屈服之后表现软化的性质. 在塑性区内, 各向同性超固结比对偏应力 q 的影响较大, 随着各向同性超固结比的增大, 偏应力 q 增大, 在弹性区内各向同性超固结比对偏应力 q 的影响较小, 可以忽略不计. 与文献 [] 的计算方法进行对比, 可以看出两种方法得到的偏应力吻合较好. q/cu R= R= R= R=. 本文方法文献 [] 8 图偏应力 q / c u 的径向分布规律当球孔扩展到某一孔径时, 土体中不同位置处的有效应力大小也不同, 由式 (8) 可知剪切模量 G 会随之发生变化. 图为 a / a = 时, 孔周土体剪切模量 G 的分布规律, 可以看出, 临界状态区与弹性区内土体的剪切模量几乎不变, 而塑性区内的剪切模量 G 变化较大. 各向同性超固结比 R 对剪切模量 G 有较大影响, 当 R < 时, 随着径向距离的增大, 剪切模量先保持不变, 然后逐渐增大, 最后保持不变 ; 当 R = 时, 随着径向距离的增大, 剪切模量 G 保持不变 ; 当 R > 时, 随着土体径向距离的增大, 剪切模量 G 先保持不变, 然后逐渐减小, 最后保持不变 ; 从图还可得到, 在临界状态区塑性区内, 各向同性超固结比越大, 剪切模量 G 越大. 当 R < 时, 随着扩孔半径的增大, 剪切模量 G 先保持不变, 然后逐渐减小, 直至稳定 ; 当 R = 时, 随着扩孔半径的增大, 剪切模量 G 保持不变 ; 当 R > 时, 随着扩孔半径的增大, 剪切模量 G 先保持不变, 然后逐渐增大, 最后保持不

5 第期李镜培, 等 : 饱和粘土中球孔扩张问题弹塑性解析 7 变. 由此可见, 扩孔过程中土体的剪切模量 G 并非常量, 而是随着球孔的扩孔半径土体的各向同性超固结比变化而变化. G/G R= R= R= R=. 图 8 剪切模量 G / G 与径向距离 / a 的关系图中, 当 a / a 时, 随着扩孔半径的增加, 扩孔压力急剧增大 ; 当 a / a 时, 随着扩孔半径的增加, 扩孔压力几乎不变, 表明扩孔压力达到极限值. 各向同性超固结比对扩孔压力的影响也比较大, 随着各向同性超固结比的增大, 扩孔压力也急剧增加. σa/cu 图 8 R= R= R=. 8 a/a 扩孔压力 σ a / c u 与扩孔半径 a / a 的关系图为扩孔半径 a / a =, 各向同性超固结比 R = 时, 超孔压 Δu 径向分布规律. 随着各向同性超固结比的增大, 孔壁处的超孔压逐渐增大, 而弹塑性边界处的超孔压减小, 当各向同性超固结比增大到一定程度时, 在弹塑性边界处将出现负孔压, 当 R = 时, 弹塑性边界处超孔压为负. 在弹性区内, 土体的超孔压为零. 文献 [] 的方法计算得到的超孔压与本文结果有一定差异. 当 R = 时, 两者差异较大. 这是因为扩孔过程中的剪切模量为初始剪切模量, 塑性区的偏应力为临界状态偏应力, 而本文认为扩孔过程中剪切模量与有效应力成正比, 且对塑性区内的偏应力没有简化, 故本文得到的超孔压更准确. 同时也说明文献 [] 的假设对超孔压的影响较大. 根据图, 当 a / a 时, 随着扩孔半径的增加, 孔壁处的超孔压 Δu 急剧增加 ; 当 a / a 时, 随着扩孔半径的增加, 孔壁处的超孔压保持不变. 各向同性超固结比对扩孔压力的影响也比较大, 随着各向同性超固结比的增大, 孔壁处的超孔压急剧增加. 当各向同性超固结比较大时, 孔壁出现负孔压. 这主要是因为当各向同性超固结比较大时, 土体表现出剪胀的特性本文方法文献 [] 7 9 (a) R = 本文方法文献 [] 7 9 (b) R = - - 本文方法文献 [] (c) R = 图超孔压 Δu / c u 与径向距离 / a 的关系 - 8 图 a/a R= R= R=. 超孔压 Δu / c u 与扩孔半径 a / a 的关系各向同性超固结比 R 对塑性区的半径有一定的影响. 图 7 表明 : 当 R 时, 随着各向同性超固结比的增大, 塑形区的半径急剧减小 ; 当 R > 时, 塑性区的半径基本保持不变.

6 7 哈尔滨工业大学学报第卷 p/a 8 8 R 图 7 塑性区半径 p / a 与 R 的关系图 8 为当扩孔半径 a / a = 时, 径向应力 σ 与切向应力 σ θ 随径向距离的变化规律. 孔壁附近土体的径向应力与切向应力基本保持不变, 这表明球孔周围土体已达到临界应力状态. 塑性区内的应力急剧增大或者减小. 在弹性区内, 随着径向距离的增大, 应力趋于稳定. 此外, 各向同性超固结比对径向应力切向应力也有显著的影响, 随着各向同性超固结比的增大, 临界状态区内的应力急剧增大, 然而在弹性区内, 径向应力与切向应力几乎不受各向同性超固结比的影响. σ/cu σ/cu σ/cu.... σ σθ 7 9 (a) R =... σ σθ. 7 9 (b) R = 7 σ σθ (c) R = 图 8 径向应力 σ / c u 切向应力 σ θ / c u 的径向分布规律图 9 为球孔周围土体中应力应变关系. 当偏应变 ε q 较小时, 此时偏应力 q 与偏应变 ε q 呈线性关系, 此时土体发生弹性变形 ; 随着球孔的扩张, 偏应变 ε q 逐渐变大, 偏应力 q 与偏应变 ε q 呈非线性关系, 此时土体已发生塑性变形. 此外, 各向同性超固结比对土体的应力应变关系具有显著的影响, 当各向同性超固结比 R = 时, 土体在屈服之后表现出硬化的性质 ; 当 R = 时, 土体在屈服之后表现理想弹塑性的性质 ; 当 R = 时, 土体在屈服之后表现软化的性质. q/kpa..... 图 9 应用分析 εq R= R= R=. 偏应力 q 与偏应变 ε q 的关系球孔扩张理论可用于桩基承载力旁压实验静力触探实验. 本文以静力触探实验为例, 分析静力触探实验锥头的极限阻力. 文中不考虑锥头的粗糙程度, 因此根据锥头的静力平衡方程可得锥头的极限阻力为 q t = σ u, () 式中 : q t 为锥头的极限阻力,σ u 为球孔的极限扩孔压力, 基于本文理论分析, 只需令 a / a, 即可得到极限扩孔压力 σ u. 为了验证该理论模型在实际应用中的可行性, 以下将以具体的静探实验为例. 本文选取文献 [] 的数据 :M = 7 v = R = ~ p = ~ kpa, 实测圆锥阻力 q c = ~ 7 kpa. 锥头贯入过程中, 由于超孔压的影响, 会使实测超孔压偏小, 因此应考虑孔压影响, 修正后的圆锥阻力 q t [7] 为 q t = q c + ( - α)u, () 式中 : q c 为实测圆锥阻力,q t 为修正后的锥头极限阻力,u 为锥头孔压,α 为净面积比值, 根据文献 [7] 的研究,α 可取为 8. 从图可看出, 利用本文理论方法得到的数据与文献 [] 的试验数据虽然有一定的误差, 但总体趋势一致, 在一定程度上, 可利用该方法预测试验数据, 因此本文理论模型有一定的实用价值.

7 第期李镜培, 等 : 饱和粘土中球孔扩张问题弹塑性解析 77 h/m 8 图 7 结论 qt/kpa 8 文献 [] 本文方法锥头极限阻力深度分布规律采用修正剑桥模型, 在不对偏应力做任何假设的条件下, 得到了不排水条件下球孔周围土体应力和孔隙水压力的半解析解. 并通过与已有解答的对比分析说明了本文研究方法的正确性以及结果的可靠性. 研究结果表明 : ) 各向同性超固结比 R 对球孔周围土体剪切模量的影响显著,R < 时, 随径向距离的增大, 剪切模量先保持不变, 后逐渐增大直至稳定 ; R = 时, 剪切模量为定值 ; 当 R > 时, 剪切模量先保持不变, 后逐渐减小直至稳定 ;R 对临界状态区塑性区内的应力及超孔隙水压力影响较大,R 越大, 切向和径向应力越大 ; 弹性区内,R 对孔隙水压力径向和切向应力几乎无影响. ) 各向同性超固结比对临界状态区塑性区内的应力超孔压影响较大, 各向同性超固结比越大, 切向应力径向应力越大 ; 在弹性区内, 各向同性超固结比对孔压径向应力切向应力几乎无影响. 此外, 各向同性超固结比越大, 塑性区的半径临界状态区的半径越小. ) 将球孔扩张理论应用到静力触探原位实验中, 对球孔扩张理论的应用性进行了说明, 可以在一定程度上促进静力触探原位实验在实际工程中的应用. 参考文献 [] 宋勇军, 胡伟, 王德胜, 等. 基于修正剑桥模型的挤密桩挤土效应分析 [ J]. 岩土力学,, ( ): 8-8. [] VESIC A S. Expansion of cavity in infinite soil mass[j]. Jounal of the Soil Mechanics and Foundation Division, ASCE, 97, 98( SM): -9. [] YU H S, HOULSBY G T. Finite cavity expansion in dilatant soils: loading analysis [ J ]. Geotechnique, 99, (): 7-8. [] BANERJEE P K, DAVIES T G, FATHALLAH R C. Behavio of axially loaded diven piles in satuated clay fom model studies [ M ]. England: Applied Science Publishes Ltd, 98. [] COLLINS I F, YU H S. Undained cavity expansions in citicalstate soils[ J]. Int J Nume Anal Meth Geomech, 99, (7): 89-. [] CAO L F, THE C I, CHANG M F. Undained cavity expansion in modified Cam clay I: theoetical analysis [ J]. Geotechnique,,(): -. [7] 胡伟, 黄义, 刘增荣, 等. 饱和粘土中挤土桩球形孔扩张的弹塑性分析 [ J]. 工程力学, 8, ( 8): [8] CHEN S L, ABOUSLEIMAN Y N. Exact undained elasto-plastic solution fo cylindical cavity expansion in modified Cam Clay soil [ J]. Geotechnique,, (): 7-. [9] CHANG M F, TEH C I, CAO L F. Undained cavity expansion in modified Cam clay II: application to the intepetation of the piezocone test [ J]. Geotechnique,, (): -. [] WOOD D M. Soil behaviou and citical state soil mechanics [ M ]. Cambidge: Cambidge Univesity Pess, 99. [] YU H S. Cavity expansion methods in geomechanics [ M]. Nethelands: Spinge,. [] 李广信. 高等土力学 [ M]. 北京 : 清华大学出版社,. []WROTH C P. The intepetation of in situ tests [ R]. Twenty Fouth Rankine Lectue Geotechnique,98, (): [ ] CHANG M F, TEH C I, CAO L F. Citical state stength paametes of satuated clays fom the modified Cam clay model [ R]. Canadian Geotechnical Jounal, 999,(): [] SILVESTRI V, ABOUSAMRA G. Application of the exact constitutive elationship of modified Cam clay to the undained expansion of a spheical cavity [ J ]. Intenational Jounal fo Numeical and Analytical Methods in Geomechanics,, (): -. []HIGHT D W, BOND A J, LEGGE J D. Chaacteization of the Bothkennaa clay: an oveview [ J]. Géotechnique, 99, (): -7. [7] 张诚厚, 施健. 孔压静力触探试验的应用 [J]. 岩土工程学报, 997, 9(): -7. ( 编辑赵丽莹 )

! Ν! Ν Ν & ] # Α. 7 Α ) Σ ),, Σ 87 ) Ψ ) +Ε 1)Ε Τ 7 4, <) < Ε : ), > 8 7

! Ν! Ν Ν & ] # Α. 7 Α ) Σ ),, Σ 87 ) Ψ ) +Ε 1)Ε Τ 7 4, <) < Ε : ), > 8 7 !! # & ( ) +,. )/ 0 1, 2 ) 3, 4 5. 6 7 87 + 5 1!! # : ;< = > < < ;?? Α Β Χ Β ;< Α? 6 Δ : Ε6 Χ < Χ Α < Α Α Χ? Φ > Α ;Γ ;Η Α ;?? Φ Ι 6 Ε Β ΕΒ Γ Γ > < ϑ ( = : ;Α < : Χ Κ Χ Γ? Ε Ι Χ Α Ε? Α Χ Α ; Γ ;