用ansys分析轨道结构
各位高手：我的毕业论文马上就要交了可是现在结果还是不对，希望高手指点，我的分析是用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 编辑 ]
页:
[1]