球面波加载下黄土线黏弹性本构关系

采用0.125 g TNT当量的微型炸药球作为爆炸源,利用在Φ1 370 mm×1 200 mm黄土样品180～1 280 m.kt-1/3范围内实测球面波径向粒子速度数据为基础,结合强间断波相容条件及变模量本构模型假设,反演得到黄土的弹性模量E=(1.927±0.216)GPa、体积模量K=(1.284±0.144)GPa及剪切模量G=(0.771±0.086)。把黄土视为线黏弹性ZWT(朱-王-唐)材料进行数值模拟,以实测粒子速度幅值vmax、位移幅值umax及对应时刻作为数值模拟结果的对比参数,定义了误差函数,其极小值对应的松弛剪切模量GM=0.13 GPa、松弛因子θM=21μs,将GM、θM作为黄土黏弹性本构参数的描述对球面波在黄土中的传播进行了数值模拟,结果表明:模拟得到的粒子速度曲线与试验曲线吻合较好,粒子速度幅值vmax、位移幅值umax与试验结果的最大偏差分别为8%、6%;粒子速度、位移幅值对应时刻与试验结果的最大偏差分别为1%、5%;数值模拟得到不同比距离处的径向应力σr、切向应力σθ、径向应变εr、切向应变εθ与基于强间断及变模量模型假设得到的结果具有较好的一致性。

Linear viscoelastic constitutive relation of loess under spherical stress wave