声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1013|回复: 4

[结构分析] 瞬态分析:请高手指点该怎么做?

[复制链接]
发表于 2009-9-2 10:03 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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

使用道具 举报

发表于 2009-9-2 10:19 | 显示全部楼层

回复 楼主 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
这是我做的一个温度载荷瞬态分析的例子,你看看对你是否有帮助。
发表于 2009-9-2 11:57 | 显示全部楼层
!先进行后处理,然后重启动

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 编辑 ]
 楼主| 发表于 2009-9-2 12:07 | 显示全部楼层

回复 板凳 jxxansys 的帖子

谢谢,瞬态分析里面还有一个时间问题,不知道该这样考虑?
发表于 2009-9-2 12:20 | 显示全部楼层
/solu
antype,TRANS
TIME,1
solve
fini


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


/solu
ANTYPE,,REST,2
time,2.1
solve
fini
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-6 09:42 , Processed in 0.065948 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表