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

采用0.125gtnt当量的微型炸药球作为爆炸源，利用在?1370mm×1200mm黄土样品180～1280m？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.13gpa、松弛因子?m=21μs，将gm、?m作为黄土黏弹性本构参数的描述对球面波在黄土中的传播进行了数值模拟，结果表明：模拟得到的粒子速度曲线与试验曲线吻合较好，粒子速度幅值vmax、位移幅值umax与试验结果的最大偏差分别为8%、6%；粒子速度、位移幅值对应时刻与试验结果的最大偏差分别为1%、5%；数值模拟得到不同比距离处的径向应力?r、切向应力??、径向应变?r、切向应变??与基于强间断及变模量模型假设得到的结果具有较好的一致性。