一种输电线路湿陷性黄土本构模拟方法技术

本发明专利技术公开了一种输电线路湿陷性黄土本构模拟方法，包括以下步骤:(1)根据黄土试样室内压缩试验和三轴不排水试验结果，测量并记录试验平均应力或轴向应力；(2)确定本构模型边界面方程中涉及到的材料常数，利用修正剑桥模型得出重塑土的边界面方程；(3)采用相关联流动法则，结合等向硬化和运动硬化法则来控制屈服面大小和中心位置的变化，利用插值计算求解塑性模量，确定本构模型；(4)选取合适的土体参数与状态变量，使模拟结果与天然土样三轴试验的结果相符合。

A constitutive simulation method for collapsible loess on transmission lines

The invention discloses a transmission line constitutive modeling methods of collapsible loess, which comprises the following steps: (1) according to the loess samples indoor compression test and the three axle undrained test results, test measure and record the average stress or axial stress; (2) to determine the constitutive constants related to the material boundary model in the equation, using the boundary equation of the modified Cambridge model that of remolded soil; (3) the associated flow rule, combined with isotropic hardening and kinematic hardening rule to control the change of yield surface size and center position, calculate the plastic modulus by interpolation, determining the model; (4) selecting soil parameters and state the appropriate variables, making the simulation results and the natural soil three triaxial test results of conformity.

【技术实现步骤摘要】
一种输电线路湿陷性黄土本构模拟方法
本专利技术涉及湿陷性黄土地基应力应变模拟方法领域，具体涉及一种输电线路湿陷性黄土本构模拟方法。
技术介绍
天然状态下的黄土强度较高，压缩性小，但在一定压力情况下遇水浸湿后，其结构会迅速破坏，这给黄土区域的输电线路杆塔基础带来严重安全隐患。以往在分析土的力学变化特性时，常以湿度、密度及粒度的变化去分析而忽略了土的结构性这一重要影响因素。在此情况下，很有必要研究结构性对输电线路工程中黄土地基的工程物理力学特性的影响规律。现有的大多数本构模型都是以重塑扰动黄土作为研究对象建立本构模型，并没有反映结构变形特性这一天然原状黄土特有的性质，这样建立的本构模型参数少，对峰值强度及结构损伤造成的土体强度衰减模拟效果不佳，因此在实际工程中应用显然是不合理的。
技术实现思路
有鉴于此，本专利技术针对饱和黄土固结阶段和剪切阶段结构受损的问题，提供了一种基于室内压缩试验和常规三轴试验结果构建结构性土本构模型的发方法，通过运动硬化结构模型引入一个可移动结构面中心位置的硬化屈服面。为了实现上述目的，本专利技术采用如下技术方案：一种输电线路湿陷性黄土本构模拟方法，包括以下步骤：(1)根据黄土试样室内压缩试验和三轴不排水试验结果，测量并记录试验平均应力或轴向应力；(2)确定本构模型边界面方程中涉及到的材料常数，利用修正剑桥模型得出重塑土的边界面方程；(3)采用相关联流动法则，结合等向硬化和运动硬化法则来控制屈服面大小和中心位置的变化，利用插值计算求解塑性模量，确定本构模型；(4)选取合适的土体参数与状态变量，使模拟结果与天然土样三轴试验的结果相符合。所述步骤(1)和步骤(2)中，对同一埋深取土的原状黄土试样展开系列常规三轴试验和室内试验，记录主要土体参数：含水率ω、液限IL、塑限IP、泊松比ν、静止土压力系数K0、前期固结应力、土体内摩擦角体积应变-平均应力对数值曲线中压缩曲线的斜率λ*、体积应变-平均应力对数值曲线中回弹曲线的斜率K*、压缩应力状态下的临界状态线斜率Mc、拉伸应力状态下的临界状态线斜率Me、土体结构各向异性参数η0。所述步骤(2)中，忽略土体初始各向异性，即η0＝0。所述步骤(2)中，采用《Aspectsoffiniteelementimplementationofcriticalstatemodels》(SHENGD，SLOANSW，YuHS.Aspectsoffiniteelementimplementationofcriticalstatemodels[J].Computationalmechanics，2000，26(2):185-196.)的形式：式中，Mc为压缩应力状态下的临界状态线斜率；m表示拉伸应力状态下临界状态线斜率Me与Mc之比；θ为Lode角所述步骤(3)中，等向硬化部分与剑桥模型一致，采用体积硬化，表示屈服面大小的pc受塑性体积应变控制：式中体积应变-平均应力对数值曲线中压缩直线对应斜率λ*和膨胀直线对应斜率κ*，pc为控制该面在偏应力空间中形状的一个标量。所述步骤(3)中，运动硬化方程通过几何条件和一致性条件求得：式中，和为屈服面和结构面中心应力张量，且n为屈服面在当前应力点处的单位梯度张量；σc为当前应力点在结构面上的像应力张量。所述步骤(3)中，塑性模量按照当前应力点与像应力点之间的距离进行插值，本专利对《Akinematichardeningconstitutivemodelfornaturalclayswithlossofstructure》(ROUAINIAM，WOODMD.Akinematichardeningconstitutivemodelfornaturalclayswithlossofstructure[J].Géotechnique，2000，50(2):153-164)采用的插值公式进行了改进，见下式：式中，Hc为像应力点对应的塑性模量，由一致性条件求得；B、ψ1、ψ2均为材料常数；b为当前应力点与像应力点之间的距离，按b＝n:(σc-σ)计算；bmax和b0分别为b的容许最大值和每次反向加载时的初始值。所述步骤(3)中，结构面与参考面大小的比值r值随着塑性应变的发展而衰减：式中，k为控制土体结构损伤衰减速率的材料常数；为描述土体结构损伤的塑性损伤应变率。所述步骤(3)中，土体结构损伤的塑性损伤应变率满足下式：式中，A为控制塑性体积应变和塑性剪切应变对结构损伤贡献率的常数，其值介于0和1之间；和分别为塑性体应变和塑性剪应变。本专利中的屈服面、结构面和参考面方程采用现有技术的三个边界面方程。所述步骤(4)是利用室内压缩试验和常规三轴试验的结果求出该本构模型中涉及到的全部材料参数和状态变量，不能由试验直接确定的参数需要调整至合理值使模拟结果与等压固结不排水试验的应力路径和应力应变关系相符合。所述步骤(4)用有限元软件的子程序对这种边界面模型进行二次开发，使用Fortran语言，采用伴随误差控制的子增量步显式积分算法完成对运动硬化结构模型的数值实现，初步建立圆柱形土体试样在等压固结不排水条件下的模型后，输入材料常数与状态变量的初始值，经过不断调试后确定材料常数和状态变量的合理取值范围，并选取最优值作为该天然黄土试样的运动硬化本构模型的一套参数。本专利技术的有益效果为：(1)本专利技术为准确构建天然结构性土的弹塑性本构模型提供了一种新方法，该模型基于土样三轴常规力学试验和室内压缩试验结果标定了一套预测其他应力路径下土样的应力应变关系，具有较高的准确性和广泛的适用性。(2)综合考虑了结构性土的屈服准则、相关联流动法则和硬化准则，良好反映结构性土的弹塑性力学行为，其力学意义明确，保证了本构模型的唯一性和准确性。(3)该模型涉及18个材料常数和12个状态变量，其中大部分常数和变量能直接确定，且编译有限元子程序进行仿真验证；由此认为该方法简单便捷，准确度较高，易于推广应用到实际工程运算中。附图说明图1为本专利技术实施的流程图；图2为运动硬化本构模型；图3为本专利技术所需的三轴常规力学试验仪；图4为卡萨格兰德经验作图法求前期固结应力；图5为K0固结试验测K0系数；图6为饱和黄土等压固结三轴不排水剪切试验结果与模拟结果。具体实施下面结合附图与实施例对本专利技术作进一步说明。如图1所示，一种输电线路湿陷性黄土本构模拟方法，依次包括以下步骤：(1)根据黄土试样室内压缩试验和三轴不排水试验结果，测量并记录试验平均应力或轴向应力；(2)确定本构模型边界面方程中涉及到的材料常数，利用修正剑桥模型得出重塑土的边界面方程；(3)采用相关联流动法则，结合等向硬化和运动硬化法则来控制屈服面大小和中心位置的变化，利用插值计算求解塑性模量，确定本构模型；(4)选取合适的土体参数与状态变量，使模拟结果与天然土样三轴试验的结果相符合。所述步骤(1)中将具有相同应力历史的土样制备成重塑土，重塑土土体内部无颗粒联结，相当于欠固结土，结构完全受损。通过重塑土试样得到的土体前期固结应力为附图2中的pc的2倍；同时该本构模型中的参考面是以修正剑桥模型为基础，讨论的土样为饱和正常固结土或弱超固结土。所述步骤(1)中，由重塑土样的室内压缩试验绘制e-lgp线，通过卡萨格兰德经验作图法得到原状土的前期固结应力Pc本文档来自技高网...

【技术保护点】
一种输电线路湿陷性黄土本构模拟方法，其特征包括以下步骤:(1)根据黄土试样室内压缩试验和三轴不排水试验结果，测量并记录试验平均应力或轴向应力；(2)确定本构模型边界面方程中涉及到的材料常数，利用修正剑桥模型得出重塑土的边界面方程；(3)采用相关联流动法则，结合等向硬化和运动硬化法则来控制屈服面大小和中心位置的变化，利用插值计算求解塑性模量，确定本构模型；(4)选取合适的土体参数与状态变量，使模拟结果与天然土样三轴试验的结果相符合。

【技术特征摘要】
1.一种输电线路湿陷性黄土本构模拟方法，其特征包括以下步骤:(1)根据黄土试样室内压缩试验和三轴不排水试验结果，测量并记录试验平均应力或轴向应力；(2)确定本构模型边界面方程中涉及到的材料常数，利用修正剑桥模型得出重塑土的边界面方程；(3)采用相关联流动法则，结合等向硬化和运动硬化法则来控制屈服面大小和中心位置的变化，利用插值计算求解塑性模量，确定本构模型；(4)选取合适的土体参数与状态变量，使模拟结果与天然土样三轴试验的结果相符合。2.根据权利要求1所述的输电线路湿陷性黄土地基本构模拟方法，其特征是：所属步骤(1)中，应对同一埋深取土的原状黄土进行重塑土制备。3.根据权利要求2所述的输电线路湿陷性黄土地基本构模拟方法，其特征是：所述步骤(1)中，通过K0固结试验测定静...

【专利技术属性】
技术研发人员：胡鑫，钱建固，刘湘莅，郭新菊，席小娟，汪秋雨，顾晓强，齐道坤，奚妍，梁晟，景川，
申请(专利权)人：国家电网公司，国网河南省电力公司经济技术研究院，同济大学，
类型：发明
国别省市：北京,11