【讨论】瞬态波动分析，那种加载方式是对的？

Jeffjia

我想研究一下瞬态波传播的问题，给一个结构施加短时间的冲击荷载，看看他的行波效应，对于时间积分Timint开关（ON,OFF）的使用不能确定，帮助上说要在积分前确定位移，速度，加速度，唯一没有说集中力怎么处理！
下面是我的加载过程：
  
（把冲击力作为净载分析。然后看瞬态效应）
  
/SOLU     
ANTYPE,TRANS
NROPT,FULL
AUTOTS,ON
TIMINT,OFF                ! STATIC LOAD STEP - DEFINE INITIAL Force！     
T1=1e-5                         ! TIME INCREMENT  
TIME,T1
DELTIM,3.0e-6                     
NSEL,S,LOC,Y,3.0
F,ALL,FY,-10
NSEL,ALL  
NCNV,2
CNVTOL,U,1,0.001
OUTPR,ALL,NONE  
OUTRES,ALL,10
SOLVE                 ! LOAD STEP 1 - STATIC
TIMINT,ON             ! Transient dynamic!           
FDELE,ALL,ALL
TIME,2e-5
DELTIM,3.0E-6         
SOLVE   
    
NEQIT,40
CNVTOL,U
CNVTOL,F,0.01,0.001
TIME,(2E-5+2.0E-3)
DELTIM,3.0E-6
/OUT,SCRATCH2
SOLVE                           ! LOAD STEP 1 - DYNAMIC
/OUT
SAVE
FINISH
  
得到冲击点速度时程曲线如下：

Jeffjia

另外一种情况，把时间积分效应从一开始就打开，命令流如下:
  
/SOLU     
ANTYPE,TRANS
NROPT,FULL
AUTOTS,ON
TIMINT,ON  ! Transient dynamic! - DEFINE INITIAL Force！     
T1=1e-5                   ! TIME INCREMENT  
TIME,T1
DELTIM,3.0e-6                     
NSEL,S,LOC,Y,3.0
F,ALL,FY,-10
NSEL,ALL  
NCNV,2
CNVTOL,U,1,0.001
OUTPR,ALL,NONE  
OUTRES,ALL,10
SOLVE                 
  
FDELE,ALL,ALL
TIME,2e-5
DELTIM,3.0E-6         
SOLVE   
         
NEQIT,40
CNVTOL,U
CNVTOL,F,0.01,0.001
TIME,(2E-5+2.0E-3)
DELTIM,3.0E-6
/OUT,SCRATCH2
SOLVE                           ! LOAD STEP 2  
/OUT
SAVE
FINISH
得到的时程曲线如下。

Jeffjia

关于模型，这里顺便描述一下。本计算模拟锚杆（钢筋D25mm）在杆头敲击下的动力反映，一共用到了3种材料，为了模拟他们在损伤下的动态反应，我把钢筋模拟成小块面积，他们是glue到一起的，几何模型如下！

tzhuzi

第二种方法是可取的吧。既然是冲击载荷，又要考虑波动效应（波的传播速度本来就很快），我觉得一开始就应该打开时间积分效应。供参考，可以看看
ansys english help 的example （vm84等），在瞬态分析中，时间效应默认是打开的。

molen

波动问题，要用LSDYNA吧。ANTYPE,TRANS恐怕不能算真正的波动问题。
集中力其实就是给一个初始加速度，
Jeffjia兄觉得呢？

tzhuzi

就如molen说的那样，要算波的效应，还是用LS-DYNA模块好，他处理小的时间步比较好，因为涉及到波动问题，时间步必然会很小。提议，考虑考虑！

Jeffjia

   这个问题正如molen兄所说的一样，用ls-dyna模块可能更好一些，具体的效果我咨询过ansys公司的技术支持，他们也不能保证，主要涉及到几种不同材料中波的传播问题，我下面准备试一试用ls－dyna来模拟一下，今天看到femlab软件的功能说明上说能很好进行wave propagation 的分析，我还想用它试一下！，希望有熟悉femlab的朋友给我些建议！
  我之所以用trans 也是看了一篇论文这么做。看来没有我想象的那么简单了。
  
谢谢大家的热忱帮助，希望多多交流。

tsintsun

果然是高手，顶！

iamtcad

呵呵，看来大家都是ansys高手，多支持，我猜刚刚入门

雨化清秋

我可是试试。。。