小弟最近在做管道凹痕分析,当我在管道端部施加弯矩(用MPC184单元),在管道内侧施加内压(均布压力),在管道另一端施加凹痕(接触)。程序仅不计算我施加的弯矩荷载对模型的影响。望各位大小不吝赐教,此问题困扰小弟太久,若能解决,不胜感激!!!
以下是我做的部分命令流,我在接触部分暂不会命令,我用GUI路径做的。
前处理:(除接触)
/prep7
et,1,solid95
mp,ex,1,2.1e5
mp,nuxy,1,0.3
TB,BKIN,1
TBDATA,1,485,769
cylind,480,500,0,1000,-90,90
sph4,0,540,0,40
vsbw,2
kwpave,9
wprota,0,90
vsbw,3
vsbw,4
wprota,0,0,90
vsbw,2
vsbw,3
vsbw,4
vsbw,5
vsbw,6
wpcsys,-1,0
wpoffs,0,0,300
vsbw,1
wprota,0,90
vsbw,6
vsbw,11
wprota,0,0,60
vsbw,12
vsbw,13
lesize,6,,,2
lesize,49,,,2
lesize,35,,,2
lesize,8,,,2
lesize,46,,,2
lesize,33,,,2
lesize,1,,,2
lesize,11,,,2
lesize,37,,,2
lesize,59,,,2
lesize,54,,,2
lesize,3,,,2
lesize,36,,,30
lesize,34,,,30
lesize,56,,,30
lesize,58,,,30
lesize,45,,,30
lesize,43,,,30
lesize,31,,,30
lesize,32,,,30
lesize,7,,,30
lesize,5,,,30
lesize,57,,,30
lesize,55,,,30
lesize,65,,,30
lesize,66,,,30
lesize,40,,,15
lesize,41,,,15
lesize,50,,,15
lesize,64,,,15
lesize,38,,,15
lesize,51,,,15
lesize,42,,,15
lesize,39,,,15
lesize,60,,,10
lesize,61,,,10
lesize,62,,,10
lesize,63,,,10
lesize,47,,,10
lesize,48,,,10
lesize,4,,,10
lesize,2,,,10
lesize,12,,,10
lesize,44,,,10
lesize,9,,,10
lesize,10,,,10
lesize,26,,,6
lesize,23,,,6
lesize,22,,,6
lesize,19,,,6
lesize,16,,,6
lesize,24,,,6
lesize,15,,,6
lesize,17,,,6
lesize,20,,,6
lesize,25,,,6
lesize,21,,,6
lesize,13,,,6
lesize,14,,,6
lesize,29,,,6
lesize,18,,,6
lesize,30,,,6
lesize,27,,,6
lesize,28,,,6
vmesh,all
et,2,184
keyopt,2,1,1
n,10000000,0,0,2000
nsel,s,loc,z,1000
type,2
*get,nnum,node,0,count
*get,ND,node,0,num,min
*do,i,2,nnum
E,10000000,ND
ND=NDNEXT(ND)
*enddo
allse
/solu
antype,0
nlgeom,on
autots,on
deltim,1e-4,1e-8,0.25,on
time,1
eqslv,sparse,,,
da,2,symm,,
da,46,symm,,
da,3,symm,,
da,29,symm,,
da,33,symm,,
da,13,symm,,
da,35,symm,,
da,50,ux,0,,
da,50,uy,0,,
da,50,uz,0,,
da,37,ux,0,,
da,37,uy,0,,
da,37,uz,0,,
da,1,ux,0,,
da,1,uy,0,,
da,1,uz,0,,
dk,9,ux,0
dk,9,uz,0
dk,10,ux,0
dk,10,uz,0
dk,14,uy,-10
dk,14,ux,0
dk,14,uz,0
sfa,32,,pres,4
sfa,52,,pres,4
sfa,48,,pres,4
sfa,39,,pres,4
sfa,5,,pres,4
sfa,30,,pres,4
f,10000000,mx,1e9
|