【原创】管道局部腐蚀模拟

pengfm

各位朋友：
　　您好！本人建立了一管道局部腐蚀的有限元模型，尽量采用映射网格划分，本人认为该模型基本合适，贴出来与各位网友共同讨论。命令流如下：
  
/prep7
  
ET,1,MESH200   
KEYOPT,1,1,7   
KEYOPT,1,2,0
  
et,2,solid95
  
mp,ex,1,2.0e6
mp,prxy,1,0.3
  
pi=3.14159265
pipe_ri=40  !管道内径
pipe_ro=50　!管道外径
pipe_th=pipe_ro-pipe_ri  !管道壁厚
pipe_l=200  !管道长度
corr_l=30　　!腐蚀区长度
corr_w=20    !腐蚀区宽度
corr_t=4　　!腐蚀区深度
frp_l=100   !粘贴的FRP的长度
frp_w=50　　!粘贴的FRP的宽度
  
!!!建立XZ平面内的纵向截面(Y=0)
k,1,pipe_ri,      
k,2,pipe_ro-corr_t
k,3,pipe_ro-corr_t,,corr_l-corr_t
k,4,pipe_ri,,corr_l-corr_t
k,5,pipe_ro,,corr_l
k,6,pipe_ro-corr_t,,corr_l
k,7,pipe_ri,,corr_l
k,8,pipe_ro,,corr_l-corr_t
k,9,pipe_ro,,frp_l
k,10,pipe_ro-corr_t,,frp_l
k,11,pipe_ri,,frp_l
k,12,pipe_ro,,pipe_l
k,13,pipe_ro-corr_t,,pipe_l
k,14,pipe_ri,,pipe_l
  
a,1,2,3,4
a,4,3,6,7
a,7,6,10,11
a,6,5,9,10
  
larc,3,5,8,corr_t*1.5   
al,5,14,11
  
a,11,10,13,14
a,10,9,12,13
  
theta1=(corr_w-corr_t)*180/(pi*pipe_ro)   
theta=corr_w*180/(pi*pipe_ro)      
theta2=frp_w*180/(pi*pipe_ro)      
  
NUMCMP,ALL   
*get,kpmax1,kp,,num,max     
k1=kpmax1+1                  
k2=kpmax1+2
  
k,k1                        
k,k2,,,corr_l
  
VROTAT,all, , , , , ,k1,k2,theta1,,   
NUMCMP,ALL  
  
VROTAT,12,16,20, , , , k1,k2,theta-theta1,,
VROTAT,24,26,30,33, , , k1,k2,theta-theta1,,
NUMCMP,ALL  
  
csys,1
cyl4,,,pipe_ri,theta,pipe_ro-corr_t,theta2   
cyl4,,,pipe_ro-corr_t,theta,pipe_ro,theta2   
cyl4,,,pipe_ri,theta2,pipe_ro-corr_t,180      
cyl4,,,pipe_ro-corr_t,theta2,pipe_ro,180
nummrg,kp
  
*get,kpmax2,kp,,num,max
k,kpmax2+1,pipe_ro,theta1
larc,51,18,kpmax2+1,corr_t*1.5
al,98,57,73
nummrg,kp
NUMCMP,ALL  
  
VDRAG,61,62,63,64,50, ,55,62　　　
nummrg,kp
NUMCMP,ALL  
  
VDRAG,85,89,93,96 , , ,67,79　　　
nummrg,kp
NUMCMP,ALL  
  
!vadd,12,24
!nummrg,kp
!NUMCMP,ALL
!aadd,95,124
!aadd,51,126
  
VPTN,12,24   
NUMCMP,ALL
  
vadd,31,32,33   
aadd,96,126     
aadd,52,129
nummrg,kp
NUMCMP,ALL
  
!!!给线分段!!!
LESIZE,16, , ,2, , , , ,0   
LESIZE,19, , ,2, , , , ,0   
  
LESIZE,10,5, , , , , , ,0   
LESIZE,17,5, , , , , , ,0
LESIZE,48,5, , , , , , ,0
LESIZE,80, , ,3, , , , ,0
LESIZE,151,5, , , , , , ,0
LESIZE,159,5, , , , , , ,0
  
TYPE,1
MSHAPE,0,2D  
MSHKEY,1
AMESH,6
AMESH,7
AMESH,3
AMESH,4
  
VSEL,S,LOC,Z,corr_l,pipe_l
type,2
mat,1
vsweep,all
allsel,all
  
LESIZE,7, , ,3, , , , ,0
LESIZE,4,4, , ,2, , , ,0   
LESIZE,14, , ,4, , , , ,0
  
TYPE,1
MSHAPE,0,2D  
MSHKEY,1
AMESH,2
AMESH,1
  
MSHKEY,0
AMESH,5
  
type,2
mat,1
vsweep,2   
vsweep,5
vsweep,1
vsweep,9
vsweep,8
vsweep,19
vsweep,14
vsweep,21
vsweep,16
  
TYPE,1
MSHAPE,0,2D  
MSHKEY,1
AMESH,93
AMESH,76
type,2
mat,1
vsweep,22
vsweep,17
vsweep,20
vsweep,15
  
LESIZE,73, , ,4, , , , ,0
TYPE,1
MSHAPE,0,2D  
MSHKEY,0
AMESH,50
type,2
mat,1
vsweep,18
  
LESIZE,43, , ,6, , , , ,0
type,2
mat,1
MSHAPE,1,3D  
vmesh,31
  
ESEL,S,TYPE,,1      
aclear,all
  
allsel,all

flic

请教一下，你建模中”frp_l=100   !粘贴的FRP的长度 frp_w=50　　!粘贴的FRP的宽度“是什么意思？