粘弹性人工边界在ANSYS中实现.docVIP

粘弹性人工边界在ANSYS中的实现 从半空间无限域取一4X2的矩形平面结构，顶部中间一定范围内受随时间变化的均布荷载，荷载如下 p(t)=t??????? 当0 DIV p(t)=2-t??????? 当1=t=2时 p(t)=0??????? 当t2时 材料弹性模量E=2.5，泊松比0.25，密度1 网格尺寸0.1X0.1，在网格边界上所有结点加法向和切向combin14号单元用以模拟粘弹性人工边界（有关理论可参考刘晶波老师的相关文章）。combine14单元的两个结点，其中一个与实体单元相连，另一个结点固定。网格图如图1所示 ? 时程分析的时间步长为0.02秒，共计算16秒。计算得到四个控制点位移时程图如图2所示，控制点坐标A(0,2)、B(0,1)、C(0,0)、D(2,2). 计算所用命令流如下： /PREP7L=4??????!水平长度H=2??????????!竖起深度E=2.5?????!弹性模量density=1???!密度nu=0.25????!泊松比dxyz=0.1???!网格尺寸 G??? = E/(2.*(1.+nu))????!剪切模量alfa = E*(1-nu)/((1.+nu)*(1.-2.*nu))??!若计算平面应力，此式需要修改Cp=sqrt(alfa/density)????!压缩波速Cs=sqrt(g/density)?