基于有效应力法的不排水条件下大直径管桩水平循环动力响应数值分析 Numerical analysis of horizontal cyclic responses of large-diameter pipe piles under undrained conditions based on effective stress method 1 2 first-author 1 2 corresp 欧强，博士(后)，从事加筋土和复合地基的研究；E-mail：ouq126@cqu.edu.cn 欧强，博士(后)，从事加筋土和复合地基的研究；E-mail：ouq126@cqu.edu.cn ouq126@cqu.edu.cn 1 2 1 2 重庆大学 土木工程学院，重庆，400045 重庆大学 土木工程学院，重庆，400045 College of Civil Engineering, Chongqing University, Chongqing 400045, China College of Civil Engineering, Chongqing University, Chongqing 400045, China 重庆大学 山地城镇建设与新技术教育部重点实验室，重庆，400045 重庆大学 山地城镇建设与新技术教育部重点实验室，重庆，400045 Key Laboratory of New Technology for Construction of Cities in Mountain Area, Chongqing University, College of Civil Engineering, Chongqing 400045, China Key Laboratory of New Technology for Construction of Cities in Mountain Area, Chongqing University, College of Civil Engineering, Chongqing 400045, China 刘锦伟 基于有效应力原理和黏土边界面塑性本构模型，开发适用于ABAQUS/Explicit的VUMAT子程序。在此基础上，研究不同循环幅值与桩径对软黏土大直径单桩水平响应的影响。研究结果表明：本文所开发程序可有效模拟循环荷载下土体的孔压和位移累积特性与刚度退化行为；循环次数与桩泥线处的桩累积变形大致呈线性关系，受加载幅值比影响；卸载刚度退化与桩累积变形成幂函数关系，且桩径增加，退化趋于严重；对于大直径海洋单桩，桩底易出现超孔压累积，其累积量受循环幅值控制；桩径对超孔压累积量不敏感，但桩径增加后其影响区域显著增加。 Based on the effective stress principle and clayed bounding surface plastic constitutive model, a material subroutine VUMAT suitable for the undrained analysis in ABAQUS/Explicit was developed, and the effects of the cyclic amplitude and pile diameter on the horizontal responses of large-diameter monopiles in soft clay were studied. The results show that the subroutine can effectively simulate the accumulation characteristics of the pore pressure and displacement of the soil and the stiffness degradation behavior of the clay under cyclic loads. The number of cycles has an approximately linear relationship with the pile's cumulative deformation at the mudline, and the deformation is affected by the load amplitude ratio, the degeneration of unloading stiffness has a power function relationship with the cumulative deformation, and the degeneration increases with the increase of pile diameter. For large-diameter monopiles, the excess pore pressure tends to accumulate at the pile bottom, whose cumulative amplitude is controlled by the cyclic amplitude ratio. However, the pile diameter is not sensitive to the accumulation of excess pore pressure, and the influence area increases significantly when the pile diameter increases. 数值模拟 二次开发 循环荷载 大直径管桩 numerical simulation secondary development cyclic loading large-diameter pipe monopole http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703350&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703337&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703361&type= 14.66796875 14.66793823 大直径单桩基础是海上风电基础的首选 [ 1 1 - 2 2 1-2 ] 。随着风机发电功率的不断增大，风电桩基础逐渐大直径化，其直径甚至可达10 m，然而，目前针对大直径海洋管桩的设计方法尚不成熟。我国东部沿海广泛分布深厚软黏土 [ 3 3 ] ，欧洲现有设计经验无法直接用于我国风电建设，水平动载下黏土大直径管桩的基础变形及承载力特性有待进一步研究。此外，与陆地不同，海床地基始终处于饱和状态，需考虑孔压对土体力学性能的影响：1) 波浪直接作用于海床，在海床表面形成动孔压边界，造成土体强度弱化 [ 4 4 ] ；2) 波浪作用在桩身，桩周土受到扰动，引起土体孔压累积；此外，随桩径增加，排水路径变长，这种累积更为明显 [ 5 5 ] 。鉴于此，基于孔压-有效应力分析的风电基础变形及承载力特性研究就显得尤为重要 [ 6 6 - 7 7 6-7 ] 。 数值模拟作为重要的研究手段在海上风电基础领域得到了广泛应用 [ 2 2 - 4 4 2-4 , 8 8 - 11 11 8-11 ] ，但往往受限于土体本构模型的选择。以ABAQUS为例，其内置的 M-C模型、D-P模型和MCC模型均无法有效地反映土体动力循环特性(如刚度退化和变形累积)，不适用于水平循环荷载下的桩-土非线性问题。此外，该软件的孔压-有效应力分析功能仅限于Standard静力分析，但Standard隐式模块对复杂接触和大变形问题(如桩-土接触和桩的贯入)收敛性极差。ABAQUS/Explicit子块非常适用于强烈非线性和大变形问题的研究。目前，ABAQUS/Explicit仅具有总应力分析功能，无法直接进行基于孔压-有效应力的不排水分析。 在以往的不排水分析中，通常采用M-C模型或D-P模型的总应力进行分析，并取泊松比 υ= 0.5来考虑不排水过程体变为0的条件。但该方法无法反映动孔压累积及土体有效应力变化情况，不满足有效应力分析的要求。为此，YI等 [ 12 12 - 13 13 12-13 ] 基于MCC模型开发了适用于ABAQUS/Explicit的不排水分析VUMAT子程序。MCC模型不能有效地反映土体动力循环特性(刚度退化和变形累积)，因此，上述VUMAT子程序得到土体孔压是不真实的，也不适用于研究海洋大直径管桩基础的变形累积问题。有研究表明，海上风电基础在其生命周期内需要承受约10 7 次负载循环 [ 14 14 ] 。风电结构对基础刚度和变形有着极其严格的要求，因此，在孔压-有效应力分析条件下，土体动本构模型二次开发对近海风电工程具有非常显著的现实意义。本文基于有效应力原理和黏土有效应力本构模型 [ 15 15 ] ，开发适用于不排水下的有效应力分析的VUMAT子程序。通过单元试验和离心机试验的验证，发现所开发程序可有效模拟循环动载下土体孔压、累积变形及刚度退化特性，适用于复杂海洋环境荷载下风电基础的受力分析。在此基础上，对水平循环动载下软黏土大直径单桩响应进行数值研究，以进一步揭示水平循环动载下大直径单桩基础变形及刚度退化规律，并给出相应的经验公式。 1 Explicit 模块下不排水条件有效应力分析的数值方法 1 1.1 　黏土本构模型 2 目前考虑黏土循环特性的本构模型中，以DAFALIAS [ 16 16 ] 为代表的边界面模型因参数意义明确和土动力学行为表征合理等优点而备受关注 [ 17 17 - 18 18 17-18 ] 。本文采用的模型为ZHOU等 [ 15 15 ] 提出的饱和黏土边界面本构模型。该模型可正确预测正常固结与超固结饱和黏土在小应变与大应变下的弹塑性响应及循环荷载下刚度退化和变形累积特征。在本文研究中，弹性剪切行为取决于当前的平均应力 p' 和当前的空隙比 e ，弹性应变增量可表示为： <math display="block" id="M1"><mi mathvariant="normal">d</mi><msubsup><mrow><mi>ε</mi></mrow><mrow><mi mathvariant="normal">v</mi></mrow><mrow><mi mathvariant="normal">e</mi></mrow></msubsup><mo>=</mo><mfrac><mrow><mi>p</mi><mi>'</mi></mrow><mrow><mi>K</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703059&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703012&type= 12.02266598 7.28133297 <math display="block" id="M2"><mi mathvariant="normal">d</mi><msubsup><mrow><mi>ε</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow><mrow><mi mathvariant="normal">e</mi></mrow></msubsup><mo>=</mo><mfrac><mrow><mi>q</mi></mrow><mrow><mn mathvariant="normal">3</mn><mi>G</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703050&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703061&type= 13.80066681 7.28133297 <math display="block" id="M3"><mi>K</mi><mo>=</mo><mfrac><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>+</mo><mi>e</mi><mo stretchy="false">)</mo><mi>p</mi><mi>'</mi></mrow><mrow><mi>κ</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703035&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703033&type= 18.96533203 7.28133297 <math display="block" id="M4"><mi>G</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">3</mn><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>-</mo><mn mathvariant="normal">2</mn><mi>υ</mi><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>+</mo><mi>υ</mi><mo stretchy="false">)</mo></mrow></mfrac><mi>K</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703056&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703066&type= 22.94466591 8.04333401 式中： <math id="M5"><mi mathvariant="normal">d</mi><msubsup><mrow><mi>ε</mi></mrow><mrow><mi mathvariant="normal">v</mi></mrow><mrow><mi mathvariant="normal">e</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703071&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703038&type= 4.65666676 4.06400013 和 <math id="M6"><mi mathvariant="normal">d</mi><msubsup><mrow><mi>ε</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow><mrow><mi mathvariant="normal">e</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703103&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703075&type= 4.48733330 4.06400013 分别为弹性体变增量和弹性剪变增量； q 和 p' 分别为剪切应力及有效围压； <math id="M7"><mi>K</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703093&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703079&type= 2.79399991 2.62466669 为土体体积模量； κ 为 e - p 空间中回弹线再压缩曲线斜率； <math id="M8"><mi>υ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703082&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703096&type= 1.69333339 2.62466669 为泊松比； G 为土体剪切模量。 边界面方程是控制土体屈服和塑性模量的关键，在 p' - q 平面内的边界面方程表示为 <math display="block" id="M9"><mi>F</mi><mo>=</mo><msup><mrow><mo stretchy="false">(</mo><mfrac><mrow><mi>q</mi></mrow><mrow><mi>M</mi><mi>p</mi><mi>'</mi></mrow></mfrac><mo stretchy="false">)</mo></mrow><mrow><mi>n</mi></mrow></msup><mo>+</mo><mfrac><mrow><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>p</mi><mi>'</mi><mo>/</mo><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">)</mo></mrow><mrow><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703085&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703099&type= 34.29000092 8.12800026 式中： <math id="M10"><mi>F</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703087&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703101&type= 2.62466669 2.62466669 为屈服面函数； n 和 r 为边界面形状参数，当分别取为1.6和2时与修正剑桥模型一致； p 0 为先期固结压力； M 为临界状态应力比。 边界面初始大小由 p 0 控制。 p 0 的增量取决于塑性体积应变，并应与修正剑桥模型保持一致。基于非关联流动法则和一致性条件，塑性应变增量可表示为： <math display="block" id="M11"><mi mathvariant="normal">d</mi><msubsup><mrow><mi>ε</mi></mrow><mrow><mi mathvariant="normal">v</mi></mrow><mrow><mi mathvariant="normal">p</mi></mrow></msubsup><mo>=</mo><msub><mrow><mi>D</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow></msub><msub><mrow><mi>Λ</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow></msub><mo>=</mo><mfrac><mrow><msup><mrow><mo stretchy="false">(</mo><mi>M</mi><mi>ρ</mi><mo>/</mo><mover accent="true"><mi>ρ</mi><mo>¯</mo></mover><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><msup><mrow><mi>η</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mn mathvariant="normal">2</mn><mi>η</mi></mrow></mfrac><msub><mrow><mi>Λ</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703138&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703134&type= 41.65599823 8.55133343 <math display="block" id="M12"><mi mathvariant="normal">d</mi><msubsup><mrow><mi>ε</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow><mrow><mi mathvariant="normal">p</mi></mrow></msubsup><mo>=</mo><msub><mrow><mi>Λ</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">p</mi></mrow></msub></mrow></mfrac><mo stretchy="false">(</mo><mfrac><mrow><mo>∂</mo><mi>F</mi></mrow><mrow><mo>∂</mo><mi>p</mi><mi>'</mi></mrow></mfrac><mi mathvariant="normal">d</mi><mi>p</mi><mo>+</mo><mfrac><mrow><mo>∂</mo><mi>F</mi></mrow><mrow><mo>∂</mo><mi>q</mi></mrow></mfrac><mi mathvariant="normal">d</mi><mi>q</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703122&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703117&type= 45.63533020 8.55133343 <math display="block" id="M13"><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">p</mi></mrow></msub><mo>=</mo><mfrac><mrow><mi>υ</mi></mrow><mrow><mo stretchy="false">(</mo><mi>λ</mi><mo>-</mo><mi>κ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mi>r</mi></mrow></mfrac><mo stretchy="false">[</mo><mfrac><mrow><msup><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">(</mo><mi>ρ</mi><mo>/</mo><mover accent="true"><mi>ρ</mi><mo>¯</mo></mover><mo stretchy="false">)</mo><mo>-</mo><msup><mrow><mi>η</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mn mathvariant="normal">2</mn><mi>η</mi></mrow></mfrac><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703152&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703124&type= 45.88933182 8.55133343 式中： <math id="M14"><mi mathvariant="normal">d</mi><msubsup><mrow><mi>ε</mi></mrow><mrow><mi mathvariant="normal">v</mi></mrow><mrow><mi mathvariant="normal">p</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703147&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703146&type= 4.65666676 4.06400013 和 <math id="M15"><mi mathvariant="normal">d</mi><msubsup><mrow><mi>ε</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow><mrow><mi mathvariant="normal">p</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703131&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703156&type= 4.65666676 4.06400013 分别为塑性体变增量和塑性剪变增量； Λ s 和 D s 分别为非负加载因子和剪胀因子； K p 为塑性模量； <math id="M16"><mi>λ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703179&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703132&type= 1.77800000 2.62466669 为 e - p 空间中初始固结线斜率； η 为应力比，定义为 q / p' ； ρ 和 <math id="M17"><mover accent="true"><mi>ρ</mi><mo>¯</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703183&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703180&type= 2.03200006 3.38666677 分别为当前应力点到原点和映射点到原点的距离。 1.2 　不排水有效应力分析 2 虽然ABAQUS/Explicit求解得到的均是土体总应力，但总应力 <math id="M18"><mfenced open="{" close="}" separators="|"><mrow><mi>σ</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703188&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703186&type= 4.91066647 3.72533321 为有效应力 <math id="M19"><mfenced open="{" close="}" separators="|"><mrow><mi>σ</mi><mi>'</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703197&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703177&type= 5.75733376 3.72533321 与土体孔压 <math id="M20"><mfenced open="{" close="}" separators="|"><mrow><mi>u</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703217&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703204&type= 4.74133301 3.72533321 之和。因此，有效应力 <math id="M21"><mfenced open="{" close="}" separators="|"><mrow><mi>σ</mi><mi>'</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703203&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703202&type= 5.75733376 3.72533321 与土体孔压 <math id="M22"><mfenced open="{" close="}" separators="|"><mrow><mi>u</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703217&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703204&type= 4.74133301 3.72533321 可分别通过定义土骨架的弹塑性刚度矩阵 <math id="M23"><mo stretchy="false">[</mo><msubsup><mrow><mi mathvariant="bold-italic">D</mi></mrow><mrow><mi mathvariant="normal">e</mi><mi mathvariant="normal">p</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msubsup><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703228&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703207&type= 7.53533268 4.74133301 和水体积模量矩阵[ D w ]进行求解，进而求得土体总应力 <math id="M24"><mfenced open="{" close="}" separators="|"><mrow><mi>σ</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703230&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703239&type= 4.91066647 3.72533321 ，见式( 9 9 )~( 11 11 )。将上述过程编译到用户自定义的材料子程序VUMAT中，便可实现不排水条件下的孔压-有效应力分析 [ 12 12 ] 。 <math display="block" id="M25"><mfenced open="{" close="}" separators="|"><mrow><mi mathvariant="normal">d</mi><mi>σ</mi></mrow></mfenced><mo>=</mo><mfenced open="{" close="}" separators="|"><mrow><mi mathvariant="normal">d</mi><mi>σ</mi><mi>'</mi></mrow></mfenced><mo>+</mo><mfenced open="{" close="}" separators="|"><mrow><mi>u</mi></mrow></mfenced><mo>=</mo><mo stretchy="false">(</mo><msubsup><mrow><mi mathvariant="bold-italic">D</mi></mrow><mrow><mi mathvariant="normal">e</mi><mi mathvariant="normal">p</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msubsup><mo>+</mo><msub><mrow><mi mathvariant="bold-italic">D</mi></mrow><mrow><mi mathvariant="normal">w</mi></mrow></msub><mo stretchy="false">)</mo><mfenced open="{" close="}" separators="|"><mrow><mi mathvariant="normal">d</mi><mi>ε</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703222&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703231&type= 55.20266724 4.31799984 <math display="block" id="M26"><mo stretchy="false">[</mo><msubsup><mrow><mi mathvariant="bold-italic">D</mi></mrow><mrow><mi mathvariant="normal">e</mi><mi mathvariant="normal">p</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msubsup><mo stretchy="false">]</mo><mo>=</mo><mo stretchy="false">[</mo><msub><mrow><mi mathvariant="bold-italic">D</mi></mrow><mrow><mi mathvariant="normal">e</mi></mrow></msub><mo stretchy="false">]</mo><mo>-</mo><mfrac><mrow><msub><mrow><mi mathvariant="bold-italic">D</mi></mrow><mrow><mi mathvariant="normal">e</mi></mrow></msub><mfrac><mrow><mo>∂</mo><mi>g</mi></mrow><mrow><mo>∂</mo><mi>σ</mi><mi>'</mi></mrow></mfrac><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><mo>∂</mo><mi>F</mi></mrow><mrow><mo>∂</mo><mi>σ</mi><mi>'</mi></mrow></mfrac></mrow></mfenced></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><msub><mrow><mi mathvariant="bold-italic">D</mi></mrow><mrow><mi mathvariant="normal">e</mi></mrow></msub></mrow><mrow><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><mo>∂</mo><mi>F</mi></mrow><mrow><mo>∂</mo><mi>σ</mi><mi>'</mi></mrow></mfrac></mrow></mfenced></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><msub><mrow><mi mathvariant="bold-italic">D</mi></mrow><mrow><mi mathvariant="normal">e</mi></mrow></msub><mfrac><mrow><mo>∂</mo><mi>g</mi></mrow><mrow><mo>∂</mo><mi>σ</mi><mi>'</mi></mrow></mfrac><mo>+</mo><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">p</mi></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703243&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703233&type= 49.95333481 19.47333336 <math display="block" id="M27"><mo stretchy="false">[</mo><msub><mrow><mi mathvariant="bold-italic">D</mi></mrow><mrow><mi mathvariant="normal">w</mi></mrow></msub><mo stretchy="false">]</mo><mo>=</mo><mfrac><mrow><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">w</mi></mrow></msub></mrow><mrow><mi>n</mi></mrow></mfrac><msub><mrow><mfenced open="[" close="]" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><mi mathvariant="bold-italic">E</mi><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mn mathvariant="normal">0</mn></mtd></mtr><mtr><mtd><mn mathvariant="normal">0</mn><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mn mathvariant="normal">0</mn></mtd></mtr></mtable></mrow></mfenced></mrow><mrow><mn mathvariant="normal">6</mn><mo>×</mo><mn mathvariant="normal">6</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703236&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703235&type= 30.39533424 8.46666718 式中：{d ε }为应变增量；[ D e ]为弹性刚度矩阵； g 为塑性势函数； K w 为水的体积模量，本文取 2.18 GPa； n 为土体孔隙率； E 和 0 分别表示为各元素为1和0的3阶矩阵。 1.3 　数值实现 2 VUMAT子程序的主要任务是基于上一步应力状态及当前给定的应变增量，直接对材料点的应力进行更新并反馈给主程序，无需求解材料的雅可比矩阵 [ 19 19 ] 。基于此，在VUMAT子程序内部，对所有与土体骨架相关参数均按有效应力状态进行更新求得有效应力，再对孔压进行单独求解，并根据叠加原理求解总应力。对于土体单元的孔压与有效应力输出可通过输出相关状态变量来实现。 数值积分算法是求解土体有效应力的关键。目前，数值积分算法主要有2大类：显式积分算法和隐式积分算法。考虑到对高度非线性问题的收敛性及计算效率，本文采用带有误差控制的改进欧拉应力积分算法 [ 20 20 ] 来求解积分点处的应力。该算法引入子步概念，根据主程序输入参数的初值，将问题划分成多个子步进行逐步求解。当应力预测值误差满足给定的控制误差时(为尽可能控制求解误差，本文取控制误差为10 -6 )，对子步中的材料点应力和状态变量进行更新。当所有子步的总步长等于1时，计算土体孔压，更新材料积分点处的总应力和状态变量，并反馈至主程序，具体流程见 图1 图1 ，其中：steptime为ABAQUS软件中分析时间；Δ T 为子步步长； T 为累积步长； T n 为第 n 次修正后子步步长； ξ 为步长修正系数。本程序采用Fortran语言进行双精度编程，共有9个材料常数和19个状态变量，可实现显式模块下饱和黏土的完全排水及不排水的动力分析。 2 数值方法验证 1 2.1 　单元试验 2 2.1.1　Itsukaichi海洋黏土不排水动三轴 3 为了评估本程序对饱和黏土静力与动力特性的描述性能，对HYODO等 [ 21 21 ] 的Itsukaichi高塑性海洋黏土动三轴试验结果进行数值重现。土样直径为50 mm，高度为100 mm，在平均围压200 kPa下固结，其土体参数见 表1 表1 。三轴不排水单调加载试验为应变控制，加载应变速率为0.001 min -1 ；单向循环试验为应力控制，加载频率为0.02 Hz。 图2 图2 所示为土体单调与循环响应的试验结果与模拟结果对比，其中： p c 为初始固结压力； q s 为初始偏应力； q cyc 和 q d 均为循环偏应力。由 图2 图2 可以看出：偏应力预测曲线与试验曲线基本吻合，VUMAT子程序可较好地模拟Itsukaichi黏土在单调荷载下不排水过程中的有效应力。对于循环荷载，偏应力计算曲线可以很好地反映黏土在加卸载过程中的位移累积特性，与试验曲线吻合较好。 2 土类试样 2 初始孔隙比 e 0 2 泊松比 υ 3 临界状态参数 M κ λ Itsukaichi黏土 2.238~2.632 0.3 1.56 0.057 0.349 高岭黏土 1.021~1.075 0.2 0.78 0.034 0.173 (a) 单调荷载下土体应力路径； (b) 循环荷载下土体应力-应变曲线 2.1.2　饱和高岭黏土不排水动三轴 3 为了进一步验证程序的有效性，分别对LI等 [ 22 22 ] 得出的高岭黏土单向与双向循环三轴试验结果进行了模拟。黏土材料基本参数见 表1 表1 。试样直径为98 mm，高为110 mm，循环加载频率为 0.1 Hz。单向加载试样平均围压为450 kPa，循环动应力幅值为116 kPa；双向加载试样平均围压为350 kPa，循环动应力幅值为130 kPa。 图3 图3 所示为孔压-循环次数关系曲线。由 图3 图3 可以看出：无论是双向还是单向循环加载，孔压均随循环次数的增加逐渐累积并趋于稳定，数值曲线与试验曲线差异较小，总体基本一致。 (a) 单向循环加载；(b) 双向循环加载 2.2 　单桩离心机试验 2 为验证本程序对侧向荷载下桩-土不排水行为预测的性能，选取文献[ 23 23 ]中的水平加载单桩离心机试验作为研究对象。试验缩尺比为1꞉48，原型桩外径为0.91 m，壁厚为50.8 mm，单桩预埋长度为20.2 m，位于泥线以上的偏心距为4.3 m。试验所使用土为Alwhite高岭土，其材料特性见 表2 表2 。原型钢桩密度为7 800 kg/m 3 ，弹性模量为206 GPa，屈服强度为414 MPa。由于模型试验的对称性，取桩土模型的一半进行建模。土体采用增强沙漏控制的C3D8R实体单元建模，单桩采用C3D8I单元。为消除边界的影响，土体半径取20倍桩径。土体底部采用固定边界，采用外侧弧面约束侧向位移；此外，限制桩和土对称平面的 Y 向位移。对浅层和靠近桩身(5 D 内， D 为单桩直径)的土体网格进行加密处理，保证桩-土分析的计算精度。采用库仑摩擦法则和考虑接触后分离的硬接触的面-面接触属性表征桩-土界面行为。根据文献[ 24 24 ]推荐公式计算桩-土界面的摩擦因数，本文取0.3。试验循环加载为单向循环位移控制，荷载频率为1 Hz。 图4 图4 所示为单调与循环荷载下数值与试验结果对比。从 图4 图4 (a)可以看出：数值曲线可以很好地模拟单桩荷载-位移曲线的加工硬化特性；此外，在深度为6 D 处的归一化循环 P - y 曲线的模拟曲线中，随着循环次数的增长，滞回曲线逐渐下移，即存在刚度与承载力退化现象，这与试验结果是一致的。由于文献[ 23 23 ]并未给出高龄黏土的边界面形状参数 n 与 r ，因此，本文采用了与修正剑桥一致的形状参数(即 n =1.6和 r =2)，这可能也是引起计算与试验误差的主要原因。虽然模拟荷载-位移曲线略低估了单桩承载性能且单圈滞回圈面积比试验结果略大，但数值与试验结果的整体趋势基本一致。 2 相对密度 2 孔隙比 e ( p' =1 kPa) 2 静止土压力 K 0 (正常固结) 2 泊松比 υ 3 临界状态参数 M κ λ 2.64 3.58 0.64 0.3 0.8 0.05 0.25 (a) 单调荷载-位移曲线；(b) 归一化的循环 P - y 曲线 3 计算结果及分析 1 在涉及桩-土相互作用的风机动力分析中，常采用子结构法，因此，水平循环荷载下大直径管桩在泥线处的刚度和位移累积规律研究是十分重要的。本节针对不同循环幅值和桩径的大直径管桩，讨论桩的位移累积、刚度退化及孔压累积规律。桩身锚固长度 L p 为30 m，壁厚取桩径的1.1%。桩土模型的建模方法与2.2节的相同，土体与管桩的基本材料参数也与文献[ 23 23 ]中的一致。桩顶荷载采用力控制的单向循环水平荷载，加载频率为 0.5 Hz，最大循环次数为100次。桩土相对刚度系数 K r 的计算及类型采用Poulos准则 [ 25 25 ] (即刚性和柔性桩的 K r 上限和下限分别为0.208 0和0.002 5)进行判别。桩研究参数见 表3 表3 。循环幅值比定义 ξ 为加载幅值 F cyc 比单桩极限荷载 F u ，其中极限荷载 F u 取泥线0.2 D 侧向位移所对应荷载 [ 2 2 ] 。 桩径 D /m 循环幅值比 ξ /% 桩土相对刚度 K r 桩型 4 20 0.048 半刚性 6 6 0.270 刚性 6 20 0.270 刚性 6 40 0.270 刚性 8 20 0.769 刚性 3.1 　荷载-位移曲线 2 图5 图5 所示为循环荷载下桩泥面处荷载-位移曲线(考虑到原荷载-位移过于密集，对其进行了过滤处理)。由 图5 图5 可以看出：在整个加载过程中，桩泥线处的侧向位移出现不可恢复的塑性变形，有明显的棘轮效应。循环幅值比对单桩的荷载-位移响应影响非常显著，即循环荷载比越大，由加载引起的土体扰动越大，桩泥线处的位移累积越大。类似地，对于同一循环荷载比下不同直径的桩，桩径越大，累积位移越大。 (a) 不同加载幅值比；(b) 不同桩径 3.2 　累积位移 2 图6 图6 所示为循环荷载下桩泥面处位移随循环次数的演化过程。由 图6 图6 可以看出：桩泥线处的侧向位移 y m 随循环次数 N 的增加而增加；桩的累积位移与循环次数基本呈线性趋势，且在100个循环内无衰减迹象，这与文献[ 26 26 ]中软黏土的离心机试验结果较为类似；另外，随循环幅值的增加，桩泥线处的位移累积速率显著增加；单桩直径越大，累积速率也越大，这与3.1节中的结果一致。使用首次循环卸载后的残余变形(荷载卸至0时的位移) y res,1 对残余侧向位移进行归一化处理，结果见 图7 图7 。从 图7 图7 可知：归一化后的累积位移受幅值比控制，但不受直径影响。为预测水平荷载下大直径桩的累积变形情况，采用对数坐标对数据进行拟合。式( 12 12 )~( 14 14 )为循环次数与累积位移的经验公式，预测公式所得结果与数值结果吻合较好。 (a) 加载幅值比；(b) 桩径 <math display="block" id="M28"><mi mathvariant="normal">l</mi><mi mathvariant="normal">g</mi><mfrac><mrow><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">s</mi><mo>,</mo><mi>N</mi></mrow></msub></mrow><mrow><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">s</mi><mo>,</mo><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfrac><mo>=</mo><mi>α</mi><msup><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">g</mi><mi>N</mi><mo stretchy="false">)</mo></mrow><mrow><mi>β</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703296&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703295&type= 25.99266624 8.80533314 <math display="block" id="M29"><mi>α</mi><mo>=</mo><mo>-</mo><mn mathvariant="normal">0.14</mn><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><msub><mrow><mi>F</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">y</mi><mi mathvariant="normal">c</mi></mrow></msub></mrow><mrow><msub><mrow><mi>F</mi></mrow><mrow><mi mathvariant="normal">u</mi></mrow></msub></mrow></mfrac></mrow></mfenced></mrow><mrow><mo>-</mo><mn mathvariant="normal">0.5</mn></mrow></msup><mo>+</mo><mn mathvariant="normal">1.01</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703297&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703310&type= 38.01533508 10.24466610 <math display="block" id="M30"><mi>β</mi><mo>=</mo><mn mathvariant="normal">0.97</mn><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><msub><mrow><mi>F</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">y</mi><mi mathvariant="normal">c</mi></mrow></msub></mrow><mrow><msub><mrow><mi>F</mi></mrow><mrow><mi mathvariant="normal">u</mi></mrow></msub></mrow></mfrac></mrow></mfenced></mrow><mrow><mn mathvariant="normal">0.18</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703321&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703311&type= 26.16200066 10.24466610 式中： y res,1 和 y res, N 分别为第1次和第 N 次循环的桩泥线处的水平累积位移； F cyc 和 F u 分别为循环荷载和静力极限荷载； N 为循环次数。 (a) 加载幅值；(b) 桩径 3.3 　卸载刚度 2 对于风机的设计，基础刚度是一个非常重要的参数。定义循环荷载-位移曲线的卸载段最高点与最低点(0荷载点)之间连线的斜率为卸载刚度。考虑到刚度的退化主要是土体塑性变形引起的，因此，相比于循环次数，桩基础的残余塑性位移更适合用于表征单桩刚度的退化规律。 图8 图8 所示为归一化后的卸载刚度 k un 与 y res 的演化关系曲线。由 图8 图8 可以看出：在循环初期，桩的卸载刚度随累积变形急剧衰减，然而，随累积变形的持续增长，刚度衰减逐渐趋于平缓，这点与文献[ 23 23 ]的报道一致。此外，不同幅值下的归一化曲线基本一致，归一化的卸载刚度退化受循环幅值比的影响很小；对于不同直径的大直径单桩，直径增加会使卸载刚度衰减更迅速，存在“直径效应”。造成这一现象的原因是：桩径越大，加载引起的土体扰动越大，刚度衰减越迅速。 式(15) 式(15) 反映了卸载刚度与累积残余位移之间的经验关系，可有效地反映出单桩直径对于刚度退化的影响，与数值结果整体吻合度较高。 <math display="block" id="M31"><mfrac><mrow><msub><mrow><mi>k</mi></mrow><mrow><mi mathvariant="normal">u</mi><mi mathvariant="normal">n</mi><mo>,</mo><mi>N</mi></mrow></msub></mrow><mrow><msub><mrow><mi>k</mi></mrow><mrow><mi mathvariant="normal">u</mi><mi mathvariant="normal">n</mi><mo>,</mo><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfrac><mo>=</mo><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">s</mi><mo>,</mo><mi>N</mi></mrow></msub></mrow><mrow><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">s</mi><mo>,</mo><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfrac></mrow></mfenced></mrow><mrow><mfenced separators="|"><mrow><mo>-</mo><mn mathvariant="normal">0.039</mn><mfrac><mrow><mi>D</mi></mrow><mrow><msub><mrow><mi>D</mi></mrow><mrow><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">f</mi></mrow></msub></mrow></mfrac><mo>-</mo><mn mathvariant="normal">0.068</mn></mrow></mfenced></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703331&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44703330&type= 42.33333206 12.86933327 式中： k un,1 和 k un, N 分别为第1次和第 N 次循环的桩泥线处的卸载刚度； D ref 为参考直径，本文取4 m。 (a) 加载幅值比；(b) 桩径 3.4 　土体孔压 2 图9 图9 给出了直径为4 m和8 m桩在100个荷载循环后的土体超孔压(状态变量SDV5)分布。由 图9 图9 可以看到：单桩底部出现明显的超孔压累积，且桩径变大后孔压区域显著增加，尤其是桩底。造成这一现象的原因是：在水平荷载下，刚性桩(或半刚性桩)桩体发生刚体转动，桩底出现“踢脚”位移，对桩底土造成扰动使土体超孔压增大。因此，在工程实际中，应考虑桩底超孔压对土体的弱化。 图10 图10 和 图11 图11 所示分别为深度3 m和桩底处(见 图9 图9 红色方框)的超孔压与循环次数关系曲线。由 图10 图10 和 图11 图11 可以看出：加载初期孔压增长迅速，但后期逐渐趋于稳定；另外，荷载增加，桩周土扰动程度加剧，超孔压幅值及其累积显著增加，且初期增长速率增加，但土体超孔压对桩径不敏感。虽然土体超孔压累积量不受桩径控制，但桩径增加后孔压影响区域显著增加(见 图9 图9 )，这意味着桩径增加会加剧整个土体性能下降，进而影响单桩的受力性能。 (a) D =4 m；(b) D =8 m (a) 加载幅值比；(b) 桩径 (a) 加载幅值比；(b) 桩径 4 结论 1 1) 不排水单元试验及水平单桩离心机试验结果表明所开发VUMAT子程序可有效模拟循环荷载下土体在不排水过程中的孔压和位移累积特性及刚度退化行为且具有较高精度。 2) 在水平循环荷载下，桩侧向位移累积随加载次数的增加而增加，有明显的棘轮效应，并与循环次数大致呈线性关系。归一化的位移累积速率会随着循环幅值的增加而显著增加，但与直径无关。 3) 随累积变形增加，单桩卸载刚度衰减速率先大后小，大致呈幂函数关系；归一化的卸载刚度和累积变形演化关系与循环幅值无关，但与直径有关，即卸载刚度的衰减存在“直径效应”，桩径越大，衰减速率越大。 4) 对于大直径海洋单桩，单桩底部会出现显著的孔压累积，工程实际应考虑桩底超孔压对土体弱化的影响。此外，超孔压累积量受循环幅值控制；超孔压对桩径不敏感，但桩径增加后其影响区域显著增加。 参考文献 CAO Guangwei , CHEN Zhixiong , WANG Chenglong , et al . Dynamic responses of offshore wind turbine considering soil nonlinearity and wind-wave load combinations [J]. Ocean Engineering , 2020 , 217 : 108155 . Dynamic responses of offshore wind turbine considering soil nonlinearity and wind-wave load combinations CAO Guangwei , DING Xuanming , YIN Zhenyu , et al . A new soil reaction model for large-diameter monopiles in clay [J]. Computers and Geotechnics , 2021 , 137 : 104311 . A new soil reaction model for large-diameter monopiles in clay 尚景宏 . 海上风力机基础结构设计选型研究 [D]. 哈尔滨 : 哈尔滨工程大学 , 2010 : 26 - 64 . 海上风力机基础结构设计选型研究 SHANG Jinghong . Study on the design and selection of offshore wind turbine foundation structure [D]. Harbin : Harbin Engineering University , 2010 : 26 - 64 . Study on the design and selection of offshore wind turbine foundation structure ZHU Bin , REN Jie , YE Guanlin . Wave-induced liquefaction of the seabed around a single pile considering pile-soil interaction [J]. Marine Georesources & Geotechnology , 2018 , 36 ( 1 ): 150 - 162 . Wave-induced liquefaction of the seabed around a single pile considering pile-soil interaction 任杰 . 波浪-大直径单桩-地基相互作用超重力模型试验与数值分析 [D]. 杭州 : 浙江大学 , 2020 : 98 - 120 . 波浪-大直径单桩-地基相互作用超重力模型试验与数值分析 REN Jie . Centrifuge and numerical studies on ocean wave-monopile-seabed interactions [D]. Hangzhou : Zhejiang University , 2020 : 98 - 120 . Centrifuge and numerical studies on ocean wave-monopile-seabed interactions LUAN Lubao , DING Xuanming , ZHENG Changjie , et al . Dynamic response of pile groups subjected to horizontal loads [J]. Canadian Geotechnical Journal , 2020 , 57 ( 4 ): 469 - 481 . Dynamic response of pile groups subjected to horizontal loads LUAN Lubao , ZHENG Changjie , KOURETZIS G , et al . Dynamic analysis of pile groups subjected to horizontal loads considering coupled pile-to-pile interaction [J]. Computers and Geotechnics , 2020 , 117 : 103276 . Dynamic analysis of pile groups subjected to horizontal loads considering coupled pile-to-pile interaction FAN Shengsheng , BIENEN B , RANDOLPH M F . Effects of monopile installation on subsequent lateral response in sand. i: pile installation [J]. Journal of Geotechnical and Geoenvironmental Engineering , 2021 , 147 ( 5 ): 04021021 . Effects of monopile installation on subsequent lateral response in sand. i: pile installation ZHANG Youhu , ANDERSEN K H . Soil reaction curves for monopiles in clay [J]. Marine Structures , 2019 , 65 : 94 - 113 . Soil reaction curves for monopiles in clay 张海洋 , 刘润 , 袁宇 , 等 . 海上大直径单桩基础 p - y 曲线修正 [J]. 水利学报 , 2020 , 51 ( 2 ): 201 - 211 . 海上大直径单桩基础p-y曲线修正 ZHANG Haiyang , LIU Run , YUAN Yu , et al . A modified p - y curve method for offshore large-diameter monopile foundations [J]. Journal of Hydraulic Engineering , 2020 , 51 ( 2 ): 201 - 211 . A modified p-y curve method for offshore large-diameter monopile foundations 费康 , 张建伟 . ABAQUS在岩土工程中的应用 [M]. 北京 : 中国水利水电出版社 , 2010 : 350 - 387 . ABAQUS在岩土工程中的应用 FEI Kang , ZHANG Jianwei . Application of ABAQUS in geotechnical engineering [M]. Beijing : China Water & Power Press , 2010 : 350 - 387 . Application of ABAQUS in geotechnical engineering YI Jiangtao , LEE F H , GOH S H , et al . Eulerian finite element analysis of excess pore pressure generated by spudcan installation into soft clay [J]. Computers and Geotechnics , 2012 , 42 : 157 - 170 . Eulerian finite element analysis of excess pore pressure generated by spudcan installation into soft clay 孟振 , 陈锦剑 , 王建华 . 基于修正剑桥模型的软黏土中沉桩过程欧拉-拉格朗日耦合模拟分析 [J]. 上海交通大学学报 , 2017 , 51 ( 3 ): 263 - 268 . 基于修正剑桥模型的软黏土中沉桩过程欧拉-拉格朗日耦合模拟分析 MENG Zhen , CHEN Jinjian , WANG Jianhua . The coupled eulerian-Lagrangian analysis of pile jacking process in saturated soft clay by using modified cam-clay model [J]. Journal of Shanghai Jiao Tong University , 2017 , 51 ( 3 ): 263 - 268 . The coupled eulerian-Lagrangian analysis of pile jacking process in saturated soft clay by using modified cam-clay model SCHAUMANN P , LOCHTE-HOLTGREVEN S , STEPPELER S . Special fatigue aspects in support structures of offshore wind turbines [J]. Materialwissenschaft und Werkstofftechnik , 2011 , 42 ( 12 ): 1075 - 1081 . Special fatigue aspects in support structures of offshore wind turbines ZHOU C , NG C W W . A thermomechanical model for saturated soil at small and large strains [J]. Canadian Geotechnical Journal , 2015 , 52 ( 8 ): 1101 - 1110 . A thermomechanical model for saturated soil at small and large strains DAFALIAS Y F . Bounding surface plasticity. I: mathematical foundation and hypoplasticity [J]. Journal of Engineering Mechanics , 1986 , 112 ( 9 ): 966 - 987 . Bounding surface plasticity. I: mathematical foundation and hypoplasticity 杨召焕 , 王建华 . 循环荷载作用下饱和软土的各向异性边界面模型 [J]. 岩土力学 , 2016 , 37 ( S1 ): 63 - 71 . 循环荷载作用下饱和软土的各向异性边界面模型 YANG Zhaohuan , WANG Jianhua . Anisotropic bounding surface model for saturated soft clay under cyclic loading [J]. Rock and Soil Mechanics , 2016 , 37 ( S1 ): 63 - 71 . Anisotropic bounding surface model for saturated soft clay under cyclic loading 胡存 , 刘海笑 . 考虑饱和黏土初始各向异性的循环边界面塑性模型 [J]. 岩土力学 , 2014 , 35 ( 10 ): 2807 - 2814, 2830 . 考虑饱和黏土初始各向异性的循环边界面塑性模型 HU Cun , LIU Haixiao . A cyclic bounding surface plasticity model for saturated clay with initial anisotropy [J]. Rock and Soil Mechanics , 2014 , 35 ( 10 ): 2807 - 2814, 2830 . A cyclic bounding surface plasticity model for saturated clay with initial anisotropy 曹光伟 . 考虑地震动力效应的混凝土结构连续倒塌破坏过程模拟 [D]. 沈阳 : 沈阳建筑大学 , 2018 : 9 - 13 . 考虑地震动力效应的混凝土结构连续倒塌破坏过程模拟 CAO Guangwei . Simulation on the progressive collapse process of concrete structure considering strain rate effect [D]. Shenyang : Shenyang Jianzhu University , 2018 : 9 - 13 . Simulation on the progressive collapse process of concrete structure considering strain rate effect SLOAN S W . Substepping schemes for the numerical integration of elastoplastic stress-strain relations [J]. International Journal for Numerical Methods in Engineering , 1987 , 24 ( 5 ): 893 - 911 . Substepping schemes for the numerical integration of elastoplastic stress-strain relations HYODO M , YAMAMOTO Y , SUGIYAMA M . Undrained cyclic shear behaviour of normally consolidated clay subjected to initial static shear stress [J]. Soils and Foundations , 1994 , 34 ( 4 ): 1 - 11 . Undrained cyclic shear behaviour of normally consolidated clay subjected to initial static shear stress LI Tao , MEISSNER H . Two-surface plasticity model for cyclic undrained behavior of clays [J]. Journal of Geotechnical and Geoenvironmental Engineering , 2002 , 128 ( 7 ): 613 - 626 . Two-surface plasticity model for cyclic undrained behavior of clays JEANJEAN P . Re-assessment of P - y curves for soft clays from centrifuge testing and finite element modeling [C]// Offshore Technology Conference . Houston, USA , 2009 : 1 - 23 . Re-assessment of P-y curves for soft clays from centrifuge testing and finite element modeling RANDOLPH M F . The response of flexible piles to lateral loading [J]. Géotechnique , 1981 , 31 ( 2 ): 247 - 259 . The response of flexible piles to lateral loading POULOS H G , HULL T S . The role of analytical geomechanics in foundation engineering [C]// Foundation Engineering: Current Principles and Practices . New York, USA : ASCE , 1989 : 1578 - 1606 . The role of analytical geomechanics in foundation engineering 何奔 . 软黏土地基单桩和复合桩基水平受荷性状 [D]. 杭州 : 浙江大学 , 2016 : 99 - 180 . 软黏土地基单桩和复合桩基水平受荷性状 HE Ben . Lateral behaviour of single pile and composite pile in soft clay [D]. Hangzhou : Zhejiang University , 2016 : 99 - 180 . Lateral behaviour of single pile and composite pile in soft clay 10.11817/j.issn.1672-7207.2023.05.024 TU473 A 1672-7207（2023）05-1905-11 曹光伟, 欧强, 丁选明, 等. 基于有效应力法的不排水条件下大直径管桩水平循环动力响应数值分析[J]. 中南大学学报(自然科学版), 2023, 54(5): 1905-1915. 曹光伟, 欧强, 丁选明, 等. 基于有效应力法的不排水条件下大直径管桩水平循环动力响应数值分析[J]. 中南大学学报(自然科学版), 2023, 54(5): 1905-1915. CAO Guangwei, OU Qiang, DING Xuanming, et al. Numerical analysis of horizontal cyclic responses of large-diameter pipe piles under undrained conditions based on effective stress method[J]. Journal of Central South University(Science and Technology), 2023, 54(5): 1905-1915. CAO Guangwei, OU Qiang, DING Xuanming, et al. Numerical analysis of horizontal cyclic responses of large-diameter pipe piles under undrained conditions based on effective stress method[J]. Journal of Central South University(Science and Technology), 2023, 54(5): 1905-1915. 国家自然科学基金资助项目(51878103)；中国博士后科学基金资助项目(2021M693740)；重庆市自然科学基金资助项目(cstc2020jcyj-cxttX0003) Project(51878103) supported by the National Natural Science Foundation of China; Project(2021M693740) supported by the China Postdoctoral Science Foundation; Project(cstc2020jcyj-cxttX0003) supported by the Innovation Group Science Foundation of the Natural Science Foundation of Chongqing, China 2022-06-07 2022-09-28 2023-05-26 2023-06-20 38909468 中南大学学报(自然科学版) 5 54