【原创实例】地层-结构三维模型的拉伸建模和网格划分

解剖地球

finish
/clear
/prep7
et,1,42  !plane42
et,2,65  !solid65
!自底向上建模思路，首先建立关键点
lsel,none,all
*get,kst_,kp,0,num,max  !get max keypoint number  kst_=0.0000
kst_=kst_+1        !kst_=1.0000
k,kst_+0,-4.926,8.424   !kst_=1.0000
k,kst_+1,-5.245,8.424   !kst_=2.0000
k,kst_+2,-4.318,8.424
k,kst_+3,-5.880,4.924
k,kst_+4,-6.166,4.924
k,kst_+5,-5.378,4.924
k,kst_+6,-4.498,1.474
k,kst_+7,-4.874,1.474
k,kst_+8,-3.727,1.474
k,kst_+9,4.498,1.474
k,kst_+10,4.875,1.474
k,kst_+11,3.727,1.474
k,kst_+12,5.879,4.924
k,kst_+13,6.149,4.924
k,kst_+14,5.379,4.924
k,kst_+15,4.928,8.424
k,kst_+16,5.246,8.424
k,kst_+17,4.322,8.424
k,kst_+18,-50.000,-25
k,kst_+19,50.000,-25
k,kst_+20,-50.000,37.951
k,kst_+21,50.000,37.951
k,kst_+22,0.001,5.344
k,kst_+23,0.001,5.349
k,kst_+24,0.002,5.361
k,kst_+25,-0.138,5.168
k,kst_+26,1.308,4.752
k,kst_+27,0.000,7.426
k,kst_+28,0.000,7.603
k,kst_+29,0.000,8.103
k,kst_+30,-1.202,4.815
k,kst_+31,-1.299,4.757
k,kst_+32,-1.308,4.752
k,kst_+33,0.275,5.099
k,kst_+34,0.097,5.174
k,kst_+35,0.085,5.180
   
!=========================连接关键点成线，画出模型（隧道+围岩）横断面
lstr,kst_+0,kst_+1
lstr,kst_+2,kst_+0
lstr,kst_+3,kst_+4
lstr,kst_+5,kst_+3
lstr,kst_+6,kst_+7
lstr,kst_+8,kst_+6
lstr,kst_+9,kst_+10
lstr,kst_+11,kst_+9
lstr,kst_+12,kst_+13
lstr,kst_+14,kst_+12
lstr,kst_+15,kst_+16
lstr,kst_+17,kst_+15
lstr,kst_+8,kst_+11
lstr,kst_+14,kst_+5
lstr,kst_+2,kst_+17
lstr,kst_+18,kst_+19
lstr,kst_+20,kst_+18
lstr,kst_+19,kst_+21
lstr,kst_+21,kst_+20
/uis,msgpop,3   !/uis,label,value  [Controls the GUI behavior]  3=Only errors displayed
!=========================连接弧线
larc,kst_+16,kst_+1,kst_+22,6.083   !larc,p1,p2,pc,rad
larc,kst_+15,kst_+0,kst_+23,5.808
larc,kst_+17,kst_+2,kst_+24,5.296
larc,kst_+13,kst_+16,kst_+25,6.292
larc,kst_+12,kst_+15,kst_+25,6.023
larc,kst_+14,kst_+17,kst_+25,5.523
larc,kst_+10,kst_+13,kst_+26,4.845
larc,kst_+9,kst_+12,kst_+26,4.575
larc,kst_+11,kst_+14,kst_+26,4.074
larc,kst_+7,kst_+10,kst_+27,7.694
larc,kst_+6,kst_+9,kst_+28,7.603
larc,kst_+8,kst_+11,kst_+29,7.606
larc,kst_+4,kst_+7,kst_+30,4.965
larc,kst_+3,kst_+6,kst_+31,4.584
larc,kst_+5,kst_+8,kst_+32,4.074
larc,kst_+1,kst_+4,kst_+33,6.443
larc,kst_+0,kst_+3,kst_+34,5.982
larc,kst_+2,kst_+5,kst_+35,5.470
kdele,23,36,1
!=========================画面
al,1,21,11,20
al,2,22,12,21
al,15,22
al,35,3,36,1
al,36,4,37,2
al,37,14,25,15
al,25,10,24,12
al,24,9,23,11
al,32,5,33,3
al,33,6,34,4
al,34,13,28,14
al,28,8,27,10
al,27,7,26,9
al,31,13
al,6,30,8,31
al,5,29,7,30
al,17,16,18,19
aovlap,all   !aovlap,NA1,NA2,NA3,NA4,NA5,NA6,NA7,NA8,NA9  [Overlaps areas]
!=========================切分面以备网格划分
wprota,,,90   !wprota,THXY,THYZ,THZX  [Rotates the working plane]  绕y轴按右手法则旋转90°
asel,s,,,all
asel,u,,,18
!asbw [Subtracts the intersection of the working plane from areas (divides areas).]
asbw,all    !asbw,NA,(blank),(blank)  利用工作平面进行面的布尔运算（切割）
allsel
wpoffs,,-5   !WPOFFS,xoff,yoff,zoff  [Offsets the working plane]工作面y负向偏5
wprota,,90
asbw,all
wpoffs,,,-21.427
asbw,all
wpoffs,,15.5
wprota,,90
asbw,all
wpoffs,,,31
asbw,all
!=========================切分线
csys,0  !csys,kcn  [Activates a previously defined coordinate system]  0=笛卡尔default
wpcsys,0 !WPCSYS,WN,KCN  [Defines the working plane location based on a coordinate system]
lsel,s,,,29,30
lsel,a,,,71,74,3
ldiv,all,,,4  !LDIV,NL1,RATIO,PDIV,NDIV,KEEP  [Divides a single line into two or more lines]
lsel,s,,,75
ldiv,all,1/8
lsel,s,,,87
ldiv,all,3/7
lsel,s,,,88
ldiv,all,3/4
lsel,s,,,76
ldiv,all,1/8
lsel,s,,,90
ldiv,all,3/7
lsel,s,,,91
ldiv,all,3/4
lsel,s,,,77
ldiv,all,1/8
lsel,s,,,93
ldiv,all,3/7
lsel,s,,,94
ldiv,all,3/4
lsel,s,,,78
ldiv,all,1/8
lsel,s,,,96
ldiv,all,3/7
lsel,s,,,97
ldiv,all,3/4
*do,i,1,3,1
l,43+i,52+i
l,46+i,49+i
l,55+i,58+i
l,61+i,64+i
*enddo
l,2,55
l,2,67
l,5,54
l,8,53
l,8,61
l,30,60
l,11,59
l,11,50
l,14,51
l,17,52
l,17,65
l,28,66
lsel,s,,,102,110,4
asbl,35,all     !asbl,NA,NL,--,KEEPA,KEEPL  [Subtracts lines from areas]
lsel,s,,,99,107,4
asbl,6,all
lsel,s,,,101,109,4
asbl,34,all
lsel,s,,,100,108,4
asbl,16,all
lsel,s,,,111,122,1
asbl,36,all
allsel
lglue,all  !LGLUE,NL1,NL2,...,NL9  [Generates new lines by "gluing" lines.]
aglue,all
nummrg,all,,,,low
numcmp,all
wpoffs,,29.951
wprota,,90
asbw,all
wpoffs,,,-4
asbw,all
asel,s,,,2,3,1
asel,a,,,14,15,1
asel,a,,,17,34,17
aadd,all  !AADD,NA1,NA2,...,NA9  [Adds separate areas to create a single area]
lcomb,22,92
lcomb,22,91
lcomb,22,90
lcomb,22,74
lcomb,22,66
allsel
csys,0
wpcsys,0
!================================以下这部分是GUI操作后从log文件中删选出来的
FLST,5,7,4,ORDE,7
FITEM,5,63  
FITEM,5,77  
FITEM,5,80  
FITEM,5,83  
FITEM,5,109
FITEM,5,118
FITEM,5,-119
LSEL,S, , ,P51X
lesize,all,,,10
FLST,5,28,4,ORDE,22
FITEM,5,38  
FITEM,5,40  
FITEM,5,-41
FITEM,5,43  
FITEM,5,-46
FITEM,5,49  
FITEM,5,-50
FITEM,5,53  
FITEM,5,-57
FITEM,5,60  
FITEM,5,84  
FITEM,5,-85
FITEM,5,87  
FITEM,5,-88
FITEM,5,93  
FITEM,5,-94
FITEM,5,120
FITEM,5,-122
FITEM,5,137
FITEM,5,139
FITEM,5,145
FITEM,5,147
LSEL,S, , ,P51X
lesize,all,,,6
lsel,s,,,1,12,1
lsel,a,,,47,48,1
lsel,a,,,51,52,1
lesize,all,,,1
FLST,5,7,4,ORDE,5   
FITEM,5,129
FITEM,5,-132
FITEM,5,135
FITEM,5,141
FITEM,5,-142
LSEL,S, , ,P51X
lesize,all,,,10
FLST,5,16,4,ORDE,15
FITEM,5,15  
FITEM,5,-16
FITEM,5,61  
FITEM,5,-62
FITEM,5,65  
FITEM,5,67  
FITEM,5,69  
FITEM,5,96  
FITEM,5,-97
FITEM,5,100
FITEM,5,-101
FITEM,5,103
FITEM,5,-105
FITEM,5,125
FITEM,5,128
LSEL,S, , ,P51X
lreverse,all
lesize,all,,,16,2
FLST,5,7,4,ORDE,7   
FITEM,5,13  
FITEM,5,-14
FITEM,5,19  
FITEM,5,31  
FITEM,5,98  
FITEM,5,102
FITEM,5,106
LSEL,S, , ,P51X
lreverse,all    !
lesize,all,,,10,2
FLST,5,7,4,ORDE,2   
FITEM,5,149
FITEM,5,-155
LSEL,S, , ,P51X
lesize,all,,,4
FLST,5,47,4,ORDE,30
FITEM,5,20  
FITEM,5,-21
FITEM,5,23  
FITEM,5,-30
FITEM,5,32  
FITEM,5,-37
FITEM,5,39  
FITEM,5,42  
FITEM,5,58  
FITEM,5,-59
FITEM,5,64  
FITEM,5,68  
FITEM,5,71  
FITEM,5,-73
FITEM,5,75  
FITEM,5,-76
FITEM,5,78  
FITEM,5,-79
FITEM,5,81  
FITEM,5,-82
FITEM,5,86  
FITEM,5,89  
FITEM,5,95  
FITEM,5,108
FITEM,5,110
FITEM,5,-117
FITEM,5,134
FITEM,5,140
FITEM,5,143
FITEM,5,148
LSEL,S, , ,P51X
lesize,all,,,6
!================================以上部分是GUI操作后从log文件中删选出来的
!=========================通过拉伸同时实现体的建立和网格划分，首先划分源面网格
allsel
asel,u,,,57
amesh,all
k,500,50,37.951,-50
l,22,500
lsel,s,,,18
lesize,all,,,50
allsel
!vdrag [Generates volumes by dragging an area pattern along a path]
vdrag,all,,,,,,18  !VDRAG,NA1,NA2,...,NA6,NLP1,NLP2,...,NLP6
!Delete lines, as well as keypoints attached to lines but not attached to other lines
ldele,18,,,1  !LDELE,NL1,NL2,NINC,KSWP  [Deletes unmeshed lines]
!=========================建立隧道上覆既有铁路路基
vsel,s,,,51
wpoffs,0,37.951,-25
wprota,0,0,60
wpoff,0,0,15.5
asbw,227
wpoff,0,0,-31
asbw,266
voffst,267,-3  !VOFFST,NAREA,DIST,KINC  [Generates a volume, offset from a given area]
vsel,s,,,64
wprota,0,71,0
vsbw,64
wprota,0,-71,0
wpoff,0,0,31
wprota,0,-71,0
vsbw,66
vdele,64,65,,1
wprota,0,71,0
wpoff,0,0,-15.5
wpoff,0,0,1.1
asbw,282
wpoff,0,0,2.8
asbw,268
wpoff,0,0,-5
asbw,266
wpoff,0,0,-2.8
asbw,272
voffst,269,-0.5
voffst,273,-0.5
vsel,s,,,67
vsbw,all
wpoffs,0,0,2.8
vsbw,all
wpoffs,0,0,2.2
vsbw,all
wpoffs,0,0,2.8
vsbw,all
lsel,s,,,366,389,23
lsel,a,,,380,381,1
lsel,a,,,404,413,9
lsel,a,,,364,408,44
lesize,all,,,4
lsel,s,,,393,394,1
lsel,a,,,385,386,1
lesize,all,,,1
lsel,s,,,369,372,3
lsel,a,,,397,399,2
lsel,a,,,360,361,1
lsel,a,,,407,412,5
lsel,a,,,374,515,41
lesize,all,,,6
lsel,s,,,368
lesize,all,0.625
vsel,s,,,64,69
vsel,a,,,71
vmesh,all
vsel,s,,,51
mshape,1
vmesh,51
aclear,all
nummrg,all,,,,low
numcmp,all
CSYS,0        
WPCSYS,0
allsel

zhangchina

技术不怎么样

husyrynb

上个图看看？