ansys时间载荷步

ANSYS各种时间步求解方法比较



ANSYS 2009-11-06 07:49:36 阅读257 评论0 字号：大中小

下面以一个弹塑性立方块求解为例.

基本建模

/prep7

Et,1,solid45

Mp,ex,1,3e4

Mp,nuxy,1,0.3

TB,BKIN,1,1,2,1

TBDATA,,20,10,,,,

Block,,10,,10,,10

Esize,1

Vmesh,1

Al,s,loc,y,0

Da,all,all





!求解部分共三种方法比较



方法1：采用do循环批处理

*DO,I,1,2

/SOLU

TIME,I

nl,s,loc,y,10

D,all,uy,-1*I

alll,all

Solve

*ENDDO



I=1时： 时间TIME从0 算到1

I=2时： 时间TIME从0算到2 ，都是基于初始构型计算的,不具有对前一步计算结果的继承性,相当于时间函数只有一个点(1,2)

最终循环计算的结果是number step= 1 ； time=2





方法2：step by step求解步

ALLSEL,ALL

/SOLU

TIME,1 !一个时间对应一个荷载步，后一步是在前一步继承计算的

NSEL,S,LOC,Y,10

D,all,uy,-1

alll,all

Solve



TIME,2

NSEL,S,LOC,Y,10

D,all,uy,-2

alll,all

Solve

该方法是在前一步计算结果的基础上继承计算的,时间函数上有两个点(1,1),(2,2)

最终计算的结果是NUMBER STEP=2, TIME=2

!如果TIME=1计算顺利完毕，采用重启动是多余的，如果中途结束，需要重启动把TIME=1算完



方法3：

用lswire文件的方法，设置荷载步(在线性计算时相当于荷载工况)和时间，此方法对前一步结果有

继承性,计算结果同方法2



方法4：

对第一种方法，要具有继承性，采取重启动的方法:

/solu

RESCONTRL,DEFINE,ALL,1,1 !保存每一子步计算结果，用于重启动

*DO,I,1,2

*IF,I,GT,1,THEN

ANTYPE,,REST !重启动以后，变量I恢复到上一步了,新的荷载步需要I=I+1,重启动文件将变量也保存了

!I=I+1

*ENDIF

TIME,I

nl,s,loc,y,10

D,all,uy,-1*I

alll,all

Solve

SAVE

*ENDDO

计算结果同方法2



方法5：是对方法4的一点补充，如果不明示TIME

/solu

RESCONTRL,DEFINE,ALL,1,1 !保存每一子步计算结果，用于重启动

*DO,I,1,2

*IF,I,GT,1,THEN

ANTYPE,,REST,I-1,LAST,0 !重启动以后，变量I恢复到上一步了,新的荷载步需要I=I+1,重启动文件将变量也保存了

I=I+1

*ENDIF

nl,s,loc,y,10

D,all,uy,-1*I

alll,all

Solve

*ENDDO

计算结果，time的值自动+1,计算结果同上



注意：1. 每个时间输入的荷载，是荷载-时间曲线上的点，也就是荷载水平，不是荷载增量，这一点与ADINA相同

2. 后继时间点一定1,3要大于前面的时间点，如果小于，则系统在前面时间自动+1

3. NUMBER STEP(在线性分析中相当于荷载工况)的数目，是通过solve来确定的，每执行solve一次，就会有一个记录，与时间没什么关系

4. 关于*do循环与Fortran语法的区别

*do,i,

1,3

i=i-1 !这里的变量不参与循环,开始就把循环次数算完了，没有递归

*enddo



!即当：i=1, i=i-1=0,此时i=0没有返回循环结构*DO中

! i=2, i=i-1=1

! i=3, i=i-1=2

但在FORTRAN中，这种写法会出现错误Error: An assignment to a DO variable within a DO body is invalid. [I],循环体中不能出现