用ansys分析轨道结构

LIEF

各位高手：
  我的毕业论文马上就要交了  可是现在结果还是不对，希望高手指点，我的分析是用ANSYS分析轨道稳定，步骤就是建立模型，对钢轨施加温度荷载，然后看钢轨的横向位移何时达到2mm，这时的温度就是临界温度，我的命令如下，请各位高手检查一下，小弟是新手上路。

/PREP7  
CSYS,1
DSYS,1 !(柱面坐标)
!*  
ET,1,BEAM3  !   
!*  
ET,2,SOLID45  
!*  
ET,3,COMBIN39   !   
KEYOPT,3,1,0
KEYOPT,3,2,0
KEYOPT,3,3,2
KEYOPT,3,4,0
KEYOPT,3,6,0 !  

R,1,0.007745,0.524e-5,0.176,0,0,0,
r,2,0,0,2e-005, 1278.4100644,0.0004,229.8621007,
rmore,0.0005,4782.7153825,0.0006,5297.8746595,0.001,7105.2951183,
rmore,0.0012,7903.4881363,0.0014,8651.3549236,0.0016,9357.7885437,
rmore, 0.002,10670.402263  

  !(单元关键项实常数)

toffst,273!             （温度单位）

MP,EX,1,2.059e11
MP,PRXY,1,0.3   
MP,DENS,1,7.83e3
MP,ALPX,1,1.18e-5  !
MP,EX,2,3.0e+10
MP,PRXY,2,0.2   
MP,DENS,2,2500   !  （材料属性）

K,1,0,0,0,  

K,2,599.25,-2.8662420,0,
K,3,599.25,2.8662420,0,  
LARC,2,3,1,599.25,   !(建立弧线)

k,4,600,2.8662420,0,
k,5,600,-2.8662420,0,
LARC,4,5,1,600,

k,6,600.75,-2.8662420,0,
k,7,600.75,2.8662420,0,
LARC,6,7,1,600.75

k,8,598.75,2.88057324,0
k,9,598.75,2.85191082,0
k,10,601.25,2.85191082,0
k,11,601.25,2.88057324,0
k,12,598.75,2.88057324,-0.3
k,13,598.75,2.85191082,-0.3
k,14,601.25,2.85191082,-0.3
k,15,601.25,2.88057324,-0.3

n,100,598.25,2.8662420,-0.15
NGEN,101,1,100, , ,0,-0.05732484,0,1,

v,8,9,10,11,12,13,14,15
lsel,,,,all,,,  
mat,2
real,
type,2
lesize,4,,,2,,,,,   
lesize,5,,,10,,,,,  
lesize,8,,,2,,,,,  
vmesh,1

VGEN,101,1, , ,0,-0.05732484,0, ,0   !(建立轨枕复制五十个)

nsel,,,,all,,,

*do,i,0,99
mat,3
real,2
type,3
e,100+i,99*i+239
*enddo            

nsel,,,,all,,,
*do,i,0,99
mat,1
real,1
type,1
e,226+i*99,226+(i+1)*99
e,232+i*99,232+(i+1)*99
*enddo                          !(建立钢轨)
NSEL,ALL
NROTAT,ALL
finish

/solu   
antype,transient
tref,0    !
nsel,,,,all,,,  
*do,i,0,100
d,100+i,,0,,,,all,,,,,         !(弹簧远离轨枕段全约束)
*enddo

d,226,,0,,,,all,,,,,
d,273,,0,,,,all,,,,,
d,10126,,0,,,,all,,,,,
d,10132,,0,,,,all,,,,,

D,ALL,UZ,
D,ALL,UY
*DO,I,4141,4340

bfE,i,temp,,0
*ENDDO
time,5  
kbc,1   
nsubst,1
allsel,all  
lswrite,1,  

  
*DO,I,4141,4340
bfE,I,temp,,50
*ENDDO
time,30
kbc,0   
nsubst,25
allsel,all  
lswrite,2      !(设置温度荷载)                        !(创建载荷部)


OUTPR,ALL,ALL
LSSOLVE,1,2,1

[ 本帖最后由 rodge 于 2007-6-10 19:31 编辑 ]