马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
小弟建一模型急需计算,可是就是映射网格老是弄不成,跪求高人出手相助,不胜感激!!!以下是命令流:
FINISH
/CLEAR
/FILNAME,CRTSII banshi ballastless track on luji,1
/prep7
ET,1,BEAM188
ET,2,COMBIN39
ET,3,COMBIN14
ET,4,COMBIN14
ET,5,SOLID45
ET,6,SOLID45
ET,7,SOLID45
ET,8,SOLID45
ET,9,SOLID45
ET,10,SOLID45
ET,11,SOLID45
KEYOPT,2,3,3 !X轴
KEYOPT,3,2,2 !Y轴
KEYOPT,4,2,1 !Z轴
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义单元实常数
R,2,0.0001,1000,0.001,10000,0.01,10000 ! COMBIN39/X
R,3,6.0e7,7.5e4 ! COMBIN14/Y
R,4,3.0e7,6e4 ! COMBIN14/Z
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义材料属性
MP,EX, 1,2.059e11
MP,PRXY,1,0.3
MP,DENS,1,7.8e3
MP,ALPX,1,1.18e-5
MP,DENS,5,2.4e3
MP,EX, 5,3.9e10
MP,PRXY,5,0.1
MP,DENS,6,2e3
MP,EX, 6,3.0e8
MP,PRXY,6,0.449
MP,DENS,7,2.5e3
MP,EX, 7,2.4e10
MP,PRXY,7,0.2
MP,DENS,8,1.95e3
MP,EX, 8,1.8e8
MP,PRXY,8,0.3
TB,DP,2,,, !输入(弹塑性)DP模型的相应参数
TBMODIF,8,1,0.7e5 !输入粘聚力
TBMODIF,8,2,34 !输入摩擦角
MP,DENS,9,1.9e3
MP,EX, 9,1.2e8
MP,PRXY,9,0.3
TB,DP,2,,, !输入(弹塑性)DP模型的相应参数
TBMODIF,9,1,0.55e5 !输入粘聚力
TBMODIF,9,2,28 !输入摩擦角
MP,DENS,10,1.85e3
MP,EX, 10,0.8e8
MP,PRXY,10,0.35
TB,DP,2,,, !输入(弹塑性)DP模型的相应参数
TBMODIF,10,1,0.5e5 !输入粘聚力
TBMODIF,10,2,30 !输入摩擦角
MP,DENS,11,1.8e3
MP,EX, 11,0.5e8
MP,PRXY,11,0.35
TB,DP,2,,, !输入(弹塑性)DP模型的相应参数
TBMODIF,11,1,0.3e5 !输入粘聚力
TBMODIF,11,2,16 !输入摩擦角
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!实体建模!!!!!!!!!!!!!!!!!!!!!
k,1,0,0.04,1.435/2
k,2,13,0.04,1.435/2
k,3,13,0.04,-1.435/2
k,4,0,0.04,-1.435/2
l,1,2
l,4,3
!!!!!!!!!!!!!!!!!!!建立基础
k,,0,0,-1.275
k,,0,-0.2,-1.275
k,,0,-0.2,1.275
k,,0,0,1.275
*get,n,kp,,num,max
l,n-3,n-2
l,n-2,n-1
l,n-1,n
l,n,n-3
*get,nl,line,,num,max
FLST,2,4,4
FITEM,2,nl-3
FITEM,2,nl-2
FITEM,2,nl-1
FITEM,2,nl
AL,P51X
k,,0,-0.2,-1.275
k,,0,-0.2,1.275
k,,0,-0.23,1.275
k,,0,-0.23,-1.275
*get,n,kp,,num,max
l,n-3,n-2
l,n-2,n-1
l,n-1,n
l,n,n-3
*get,nl,line,,num,max
FLST,2,4,4
FITEM,2,nl-3
FITEM,2,nl-2
FITEM,2,nl-1
FITEM,2,nl
AL,P51X
k,,0,-0.23,-1.625
k,,0,-0.23,1.625
k,,0,-0.53,1.625
k,,0,-0.53,-1.625
*get,n,kp,,num,max
l,n-3,n-2
l,n-2,n-1
l,n-1,n
l,n,n-3
*get,nl,line,,num,max
FLST,2,4,4
FITEM,2,nl-3
FITEM,2,nl-2
FITEM,2,nl-1
FITEM,2,nl
AL,P51X
k,,0,-0.53,-4.4
k,,0,-0.53,4.4
k,,0,-1.23,5.45
k,,0,-1.23,-5.45
*get,n,kp,,num,max
l,n-3,n-2
l,n-2,n-1
l,n-1,n
l,n,n-3
*get,nl,line,,num,max
FLST,2,4,4
FITEM,2,nl-3
FITEM,2,nl-2
FITEM,2,nl-1
FITEM,2,nl
AL,P51X
k,,0,-1.23,-5.45
k,,0,-1.23,5.45
k,,0,-3.53,8.9
k,,0,-3.53,-8.9
*get,n,kp,,num,max
l,n-3,n-2
l,n-2,n-1
l,n-1,n
l,n,n-3
*get,nl,line,,num,max
FLST,2,4,4
FITEM,2,nl-3
FITEM,2,nl-2
FITEM,2,nl-1
FITEM,2,nl
AL,P51X
k,,0,-3.53,-8.9
k,,0,-3.53,8.9
k,,0,-4.53,10.4
k,,0,-4.53,-10.4
*get,n,kp,,num,max
l,n-3,n-2
l,n-2,n-1
l,n-1,n
l,n,n-3
*get,nl,line,,num,max
FLST,2,4,4
FITEM,2,nl-3
FITEM,2,nl-2
FITEM,2,nl-1
FITEM,2,nl
AL,P51X
k,,0,-4.53,-15
k,,0,-4.53,15
k,,0,-10.53,15
k,,0,-10.53,-15
*get,n,kp,,num,max
l,n-3,n-2
l,n-2,n-1
l,n-1,n
l,n,n-3
*get,nl,line,,num,max
FLST,2,4,4
FITEM,2,nl-3
FITEM,2,nl-2
FITEM,2,nl-1
FITEM,2,nl
AL,P51X
vext,all,,,13,,,,,
WPCSYS, -1
kwpave,1
vSBW, all
kwpave,4
vSBW, all
kwpave,5
vSBW, all
kwpave,8
vSBW, all
kwpave,13
vSBW, all
kwpave,14
vSBW, all
wprota,,,90
*do,i,1,19,1
wpoffs,0,0,0.65
ASBW,all
*enddo
!!!!!!!!!!!!!!!!!分配属性!!!!!!!!!!!!!!
vsel,s,loc,y,-0.2,0
type,5
mat,5
vsel,s,loc,y,-0.23,-0.2
type,6
mat,6
vsel,s,loc,y,-0.53,-0.23
type,7
mat,7
vsel,s,loc,y,-1.23,-0.53
type,8
mat,8
vsel,s,loc,y,-3.53,-1.23
type,9
mat,9
vsel,s,loc,y,-4.53,-3.53
type,10
mat,10
vsel,s,loc,y,-19.53,-4.53
type,11
mat,11
怎么才能映射划分?纠结了N天了,急求指点 |