【讨论】请教关于冲击波的计算过程

zxhst

问题如下：假定问题是一维的，简图如下四节点三单元，材料是普通弹塑性材料，弹性波在一个时间步长dt内恰好走过一个单元的长度dl，初始静止，在0时刻左端突加一撞击速度v（足以使其产生塑性），
            v  =〉。——。——。——。
那么在第一个dt中第一个单元应变e=0+de
            de=v*dt/dl
  下一步就是如何由应变增量de 求出应力增量 dp
  
在用户子程序中的计算方法都是，先假设是完全弹性变形，求出一个增量试探应力dp-trial,——〉p-trial再求出其等效应力，然后与相应的屈服应力p-0相比较，超过则表明发生了塑性变形，需要修正回去再得到真实应力。
  
我觉得这种方法好像只能适用于小的应变增量，如果在一个时间步长内，应变增量特别大时比如达到1%，应该如何处理？
  
通常abaqus材料定义中只是输入几个表征塑性的应力应变点，不知道具体过程是怎样的？
  
恳请高手指点，不胜感激！！！