/config,nres,50000


l=500    !钢轨长度m
a=1.25     !支承间距m
le=0.125      !单元长度  
ne=l/le   !ne of rail
nk=l/a     !跨数

        
jld=2006  ! 加力节点

/PREP7
!*******************element
ET,1,beam188            !establish rail element beam188
ET,2,COMBIN14,,,1       !spring under rail



!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!if rail beam188 read beam *******Timo
SECTYPE,6,BEAM,MESH,
SECOFFSET,CENT,,,   
SECREAD,beam,sect,,MESH
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!


R,2,1.78e8,1e2              !spring under ganggui 10e6,4.0e6

!materials   
MP,EX,1,2.059e11             !定义钢轨弹性模量
MP,PRXY,1,0.3                !定义钢轨泊松比
MP,DENS,1,7.83e3             !定义钢轨密度


!rail model and mesh
N,1,0,0                      !defination rail node(first point)
N,ne+1,0,0,l                     !last point of rail/total length of rail   
FILL
type,1
real,1
mat,1
SECNUM,6
*do,i,1,ne           !establish rail element(rail have 2430 (or 810)element)
E,i,i+1
*enddo

N,9001,0.2,0          !poistion of formation node
N,9001+nk,0.2,0,l
FILL
!establish formation element
TYPE,2
REAL,2
*do,i,1,1+nk
e,10*(i-1)+1,i+9000
*enddo
!*****************************************************yueshu
D,1,ALL              !bounding start point of rail
D,ne+1,ALL         !bounding end   point of rail

*DO,i,9001,9001+nk   
D,i,all
*enddo





!*****************************************************yueshu




/SOLU   
ANTYPE,HARMIC	! Harmonic response analysis
HROPT,FULL   !FULL - Full method (default).REDUC - Reduced method,MSUP - Mode superposition method
HROUT,off      !Print results as amplitudes and phase angles
LUMPM,0         !1?
EQSLV,FRONT,1e-008, 
PSTRES,0                !pre-stress

ALPHAD,0,   
BETAD,0,
DMPRAT,0.01, 
 
HARFRQ,0,2000,             
NSUBST,500,
outres,all,all           !
KBC,1                     
 
!adding force point jld 
F,jld,Fx,10000,             !!!%%!!!! value of force  (one axle load of metro train)



SOLVE 
finish

截面名称为beam.sect,但是截面文件传不上来