不同应力路径下粉煤灰本构模型研究.docVIP

不同应力路径下粉煤灰本构模型研究 胡德金　高正中　张青云(四川联合大学水利工程系) 摘　要　本文研究粉煤灰在不同等应力增量比路径下的应力-应变-体变特性，提出考虑土的应变软化、剪胀(剪缩)性和应力路径影响的非线性弹性模型. 通过与试验及有限元计算成果对比，验证了模型的适用性.关键词　应变软化，剪胀(剪缩)，应力路径，应力比，本构模型，有限元. 　　土的力学特性是复杂的，具非线性、弹塑性、剪胀性、各向异性和流变性，同时应力状态、应力路径及土的物理状态组成、结构等对其均有影响. 土的本构模型就是用数学手段来体现试验中所揭示的土的这些复杂力学特性.现常用的本构模型有弹性模型和弹塑性模型，后者较为复杂. 本文结合粉煤灰提出考虑土的应变软化、剪胀(剪缩)性和应力路径影响的非线性弹性模型. 不少国内外学者指出，土工问题不但多属于三维和平面应变问题，且土石坝工程中的应力路径更接近于比例加载的路径. 本文结合某发电厂粉煤灰(分类属中液限粉土)，进行等应力增量比(η0=C) (η0为偏应力q与平均有效应力p增量之比)的三轴剪切试验，以揭示土的应力-应变特性. 试验中分别进行以下几种应力路径： p=C(η0→∞), σ3=C(η0=3.0), η0=2.591,2.250,1.714. 试样按等压固结，以等应力增量比控制，直至剪切破坏. 试验结果表明，应力路径对土的变形特性有显著影响. 1　模型的提出　　根据大量试验成果，通过分析计算，提出考虑土体应变软化、剪胀(剪缩)性及应力路径影响的非线性弹性模型为， 　　　　　　　　　　　　　(1) 　　　　　　　　　　　　(2) 式中：Kt为平均有效应力p与偏应力q藕合体变模量，Ktp为平均有效应力体变模量，Ktp为偏应力体变模量，Gt为剪切模量. 2　模型参数的确定 2.1　剪切模量Gt的确定　试验土样共进行密实(筑坝灰，γd=12.2kN/m3)和中密(库灰，γd=11.0kN/m3)两种干容重的等应力增量比(η0=C)加载试验，其应力-应变关系存在不同程度的软化特征(如图1)，采用的应力-应变关系表达式为： 图1　土样应变软化曲线 　　　　　　　　　　　　　　(3) 当b=c时，转化为 　　　　　　　　　　　　　　　　(4) 即为硬化型应力-应变关系. 　　　　　　　　　(5) a,b,c——模型参数，与固结围压、应力路径有关. 　　由应力-应变关系曲线特征点值可求解模型参数a,b,c. 　　　　　　　　　　　　(6) c=b2qr????????????? (7) α=（b-2c)εsf (8) qf——峰值强度， qr——残余强度， εsf——峰值应变. 　　qf,qr,εsf分别与初始固结围压σ30、应力路径(η0)及材料特性有关. 在一定应力路径下(η0=C), qf,qr与σ30在双对数座标图中具良好线性关系. εsf与σ30亦呈良好线性关系，相关系数均在0.99以上. lgqf=c1+c2lgσ30 } (9) lgqr=c3+c4lgσ30 εsf=c5+c6lgσ30 　　c1、c2、c3、c4、c6随η0的变化规律可由1/η0的线性关系表示： c1=c11+c12/η0 } (10) c2=c21-c22/η0 c3=c31+c32/η0 c4=c41-c42/η0 c6=c61-c62/η0 c5主要与材料特性有关. 2.2　平均有效应力体变模量Ktp的确定　Ktp参照成都科大简化K-G模型［1］ ktp=ki+αk p.　　　　　　　　　　　　(11) 表1　试验土样参数ki,αk 参数 密实饱和 密实非饱和 中密饱和 K/100kPaαk 31.3977.0 40.5923.23 121.7282.11 　　根据一个试样多级加荷等压固结试验的p-εvp曲线，采用数解法计算曲线各切点斜率与p的线性关系表达式，或采用程序利用牛顿法均可求解参数Ki、αk. 试验土样参数ki,αk见表1.2.3　偏应力体变模量Ktq的求取　从p=C三轴试验成果可以看出，纯粹的剪切也会引起体积变化，这说明体应变与平均有效应力p和偏应力q有关. 总的体应变应由p和q分别引起的两部分体应变组成.即 　　　　　　　　　　　　　(12) 　　根据有关研究文献［2-4］，塑性体应变与塑性剪应变之比是应力比的函数. 　　　　　　　　　　　　　　　(13) 　　假定一切剪应变都不可恢复(弹性剪应变很小，可忽略不计)，即剪应变中的弹性部分等于零，从而推导得， 　　　　　　　　　　(14) 　　　　　　　　　　　　(15) 　　弹性体应变为 　　　　　　　　　　　(16) 　　　　　　　　　　　　　　(17) Cs——土体膨胀模量. 　　　　　　　　　(18) 　　整理应力比η-