瞬态分析：请高手指点该怎么做？
各位大侠，我是做一个瞬态分析的工程，需要在一个时间步长内进行迭代，也就是说在一个节点加一次荷载（荷载是与结构的位移、速度、加速度有关），算出位移、速度、加速度等以后，在时间不变的情况下，利用上一次计算出来的位移等算出荷载，重新在此节点加算出来的荷载，如此循环，不知道该怎么做，请高手指点！回复 楼主 ydcow 的帖子
/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
这是我做的一个温度载荷瞬态分析的例子，你看看对你是否有帮助。 !先进行后处理，然后重启动
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 编辑 ]
回复 板凳 jxxansys 的帖子
谢谢，瞬态分析里面还有一个时间问题，不知道该这样考虑？ /soluantype,TRANS
TIME,1
solve
fini
/solu
ANTYPE,,REST,1
time,2
solve
fini
/solu
ANTYPE,,REST,2
time,2.1
solve
fini
页:
[1]