提出一种研究复杂地形二维粘弹性介质中波动问题的数值模拟方法。利用标准线性固体(SLS)微分型本构方程给出计算当前应力的公式,基于三角形网格给出控制体的动力平衡方程,在时间域上逐步递推计算给出各节点的加速度、速度、位移和各格子的应力。与解析解的对比结果验证了方法的正确性。该方法自动满足复杂地形的自由表面条件。研究了地震作用下北川县城区域的剪应力、加速度场的动态集中现象,分析了震害原因。
A numerical method is presented for simulating wave propagation in a viscoelastic media with an irregular boundary.The standard linear solid (SLS) mechanism is used to characterize the viscoelastic media.The stress formulae for recursive calculation are derived from SLS constitutive relations.The dynamic equilibrium equations of controlled volumes are established by using the mesh of triangular grid.By applying the dynamic equilibrium equations of controlled volumes and the stress formulae alternately,the nodal accelerations,velocities,displacements,and the stresses in each grid can be obtained in time domain.The result comparing with the analytical solution shows that the numerical algorithm presented is an effective method.The free boundary condition of arbitrary irregularity is satisfied naturally.The dynamic concentration phenomena of shear stress and accelerations are investigated for the area of Beichuan county town with a topography surface.The seismic disaster reasons of the Beichuan county town are analyzed.