瞬态分析：请高手指点该怎么做？

ydcow

各位大侠，我是做一个瞬态分析的工程，需要在一个时间步长内进行迭代，也就是说在一个节点加一次荷载（荷载是与结构的位移、速度、加速度有关），算出位移、速度、加速度等以后，在时间不变的情况下，利用上一次计算出来的位移等算出荷载，重新在此节点加算出来的荷载，如此循环，不知道该怎么做，请高手指点！

macrolea

/COM,THERMAL
! 进入前处理器
/prep7  

et,1,PLANE77

KEYOPT,1,3,1

!定义材料属性
MP,KXX,1,0.05   ！定义导热系数
MP,C,1,400  ！比热
MP,DENS,1,600   ！密度

!定义属材料属性
MP,KXX,2,160
MP,C,2,875
MP,DENS,2,2760

!定义材料属性
MP,KXX,3,0.1
MP,C,3,1254
MP,DENS,3,1500

TOFFST,273   ！定义温度为摄氏度

k,1000,0,0.005,0
k,1001,0,0.0055,0
k,1002,0,0.0275,0
k,1003,0,0.028,0

k,1,0.02825,0,0
k,2,0.03773,0,0
k,3,0.03325,0.005,0
k,4,0.03773,0.005,0
k,5,0.03325,0.028,0
k,6,0.03773,0.028,0
K,7,0.02825,0.033,0
K,8,0.03773,0.033,0

L,1,2
L,3,4
L,5,6
L,7,8
L,2,4
L,6,8
L,1,7
L,3,5
AL,ALL

kwpave,1000
!画下面
RECTNG,0.03325,0.03773,0,0.0005



kwpave,1001
!画

RECTNG,0.03325,0.03328,0,0.022
RECTNG,0.03353,0.03356,0,0.022
RECTNG,0.03381,0.03384,0,0.022
RECTNG,0.03409,0.03412,0,0.022
RECTNG,0.03437,0.03440,0,0.022
RECTNG,0.03465,0.03468,0,0.022
RECTNG,0.03493,0.03496,0,0.022
RECTNG,0.03521,0.03524,0,0.022
RECTNG,0.03549,0.03552,0,0.022
RECTNG,0.03577,0.03580,0,0.022
RECTNG,0.03605,0.03608,0,0.022
RECTNG,0.03633,0.03636,0,0.022
RECTNG,0.03661,0.03664,0,0.022
RECTNG,0.03689,0.03692,0,0.022
RECTNG,0.03717,0.03720,0,0.022
RECTNG,0.03745,0.03748,0,0.022



RECTNG,0.03328,0.03353,0,0.022
RECTNG,0.03356,0.03381,0,0.022
RECTNG,0.03384,0.03409,0,0.022
RECTNG,0.03412,0.03437,0,0.022
RECTNG,0.03440,0.03465,0,0.022
RECTNG,0.03468,0.03493,0,0.022
RECTNG,0.03496,0.03521,0,0.022
RECTNG,0.03524,0.03549,0,0.022
RECTNG,0.03552,0.03577,0,0.022
RECTNG,0.03580,0.03605,0,0.022
RECTNG,0.03608,0.03633,0,0.022
RECTNG,0.03636,0.03661,0,0.022
RECTNG,0.03664,0.03689,0,0.022
RECTNG,0.03692,0.03717,0,0.022
RECTNG,0.03720,0.03745,0,0.022
RECTNG,0.03748,0.03773,0,0.022



kwpave,1002
!画上面
RECTNG,0.03325,0.03773,0,0.0005


AGLUE,ALL
NUMCMP,LINE
NUMCMP,AREA
NUMCMP,NODE
kdele,1000,1003
WPAVE,0,0,0
ALLSEL,ALL
APLOT
!体33为铝
!*  

ASEL,S, , ,33

AATT, 2, ,1,0   !赋予所选体元属性
ALLSEL,ALL

!体17-32为光纤环
!*  

ASEL,S,,,17,32

AATT,1, ,1, 0   
ALLSEL,ALL
!体1-32,66,67为垫纸
!*  

ASEL,S,,,1,16
ASEL,A,,,34,35

AATT,3, , 1,0   
ALLSEL,ALL


!!选择光纤环
ASEL,S,,,17,32
Aplot  

!层数
LSEL,S,,,72,102,2
!LPLOT
LESIZE,ALL,,,1
ALLSEL,ALL
!匝数
LSEL,S,,,8
LSEL,A,,,11,67,4
!LPLOT
LESIZE,ALL,,,88
ALLSEL,ALL


ASEL,S,,,17,32

MSHAPE,0,2D
MSHKEY,1
AMESH,ALL
ALLSEL,ALL
!!选择中间垫纸层
ASEL,S,,,1,16
APLOT

!!层数
LSEL,S,,,70
LSEL,A,,,10,66,4
!LPLOT
LESIZE,ALL,,,1
ALLSEL,ALL
!!匝数
LSEL,S,,,71
LSEL,A,,,13,69,4
!LPLOT
LESIZE,ALL,,,88
ALLSEL,ALL

ASEL,S,,,1,16
MSHAPE,0,2D
MSHKEY,1
AMESH,ALL
ALLSEL,ALL


!
!底层垫纸
ASEL,S,,,34
APLOT

LSEL,S,LOC,Y,0.0055
LCCAT,ALL
ALLSEL,ALL

!层数
LESIZE,2,,,32

!匝数
LESIZE,105,,,2

MSHAPE,0,2D
MSHKEY,1
AMESH,34

!顶层垫纸
ASEL,S,,,35
APLOT

LSEL,S,LOC,Y,0.0275
LCCAT,ALL
ALLSEL,ALL

!层数
LESIZE,3,,,32

!匝数
LESIZE,106,,,2

ALLSEL,ALL

MSHAPE,0,2D
MSHKEY,1
AMESH,35

ALLSEL,ALL

!划分网格
ESIZE,0.0008,0,
MSHAPE,1,2D
MSHKEY,0
AMESH,33

!删除连接线
CM,_Y,LINE  
*SET,_Z1,LSINQR(0,13)   
*IF,_z1,ne,0,then   
LSEL,R,LCCA
*SET,_Z2,LSINQR(0,13)   
LDELE,ALL   
*SET,_Z3,_Z1-_Z2
*IF,_Z3,NE,0,THEN   
CMSEL,S,_Y  
CMDELE,_Y   
*ENDIF  
*ELSE   
CMSEL,S,_Y  
CMDELE,_Y   
*ENDIF

!施加温度激励
/SOL
antype,trans
lumpm,0

time,1  
timint,on   
deltim,1/60,,,1
outres,erase
outres,all,all  
autots,off  
kbc,0

tunif,298
LSEL,S,,,6
LSEL,A,,,107,108
LSEL,A,,,103
DL,ALL,,TEMP,353
ALLSEL,ALL
SFL,7,CONV,,,298


ALLSEL,ALL
LSWRITE,1
!SOLVE

time,7
timint,on   
deltim,1/60,,,1
outres,erase
outres,all,all  
autots,off  
kbc,1

LSEL,S,,,6
LSEL,A,,,107,108
LSEL,A,,,103
DL,ALL,,TEMP,353
ALLSEL,ALL
SFL,7,CONV,,,298
ALLSEL,ALL
LSWRITE,2

LSSOLVE,1,2
FINISH
这是我做的一个温度载荷瞬态分析的例子，你看看对你是否有帮助。

jxxansys

!先进行后处理，然后重启动

fini
/CLEAR        
/FILNAME, EXAMPLE

/PREP7
ET,1, BEAM3        
R,1,14.345e-4,245e-8,0.1         
MP,EX,1,2E11
MP,NUXY,1,0.3        
K,1,0,0,0        
K,2,1,0,0        
LSTR,1,2
LESIZE,1,,,50        
LMESH,1        
FINISH
/SOLU        
DK,1,UX        
DK,1,UY        
DK,1,ROTZ        
FK,2,FY,-10000        
SOLVE        
FINISH        

/POST1        
NSORT,U,Y   
*GET,UY_MIN,SORT,,MIN   !计算位移
fini

/solu
ANTYPE,,REST,1      !!!!!!!!!!!!!!!!!!!!!!!!!!!
FK,2,FY,UY_MIN*1e6   !施加与位移相关的载荷
SOLVE        
FINISH

[ 本帖最后由 jxxansys 于 2009-9-2 12:00 编辑 ]

ydcow

谢谢，瞬态分析里面还有一个时间问题，不知道该这样考虑？

jxxansys

/solu
antype,TRANS
TIME,1
solve
fini


/solu
ANTYPE,,REST,1
time,2
solve
fini


/solu
ANTYPE,,REST,2
time,2.1
solve
fini