本帖最后由 zhuizhu18 于 2012-10-10 09:13 编辑
花了几个小时的时间给你重新划分了一下网格。
1)、你这个问题应该是个对称问题,所以我只取了1/8球体划分的网格。
2)、如果计算的话,你可以在四周加上对称约束就行了,能减少不少工作量。
3)、如果就想用1/2球体,那你就再镜像吧。
4)、因为你底部的梁网格那块划分的比较细,虽然对其周边的网格进行了过渡处理,可以网格密度还是比较大,你可以根据自己的需要重新调一下。
下面的是改过的命令流,你自己看看吧,我用的13.0的版本,不知道能不能在你那运行,哈哈。
网格密度比较大,网格如果有划分的不合适的地方,或者可以改进的地方,各位大神还望不吝赐教,小弟在此先谢过了。(真想有个积分啊,盼了好长时间了.....)
fini
/clear
/title,11
/prep7
MP,EX,1,2.1E11
MP,DENS,1,7800
MP,NUXY,1,0.3
ET,1,SOLID45
blc4,,,0.01,1,0.01
alls
vgen,20,all,,,0.01
FLST,3,1,6,ORDE,1
FITEM,3,20
VGEN,8,P51X, , , , ,-0.010, ,0
FLST,3,1,6,ORDE,1
FITEM,3,27
VGEN,4,P51X, , ,-0.010, , , ,0
vsel,s,loc,z,0,0.010
vgen,8,all,,,20*0.01,,
vsel,s,loc,z,0,-0.100
vgen,7,all,,,20*0.01,,
wpoffs,0.01*20*8/2,1/2,0.01
wprot,0,90,0
SPHERE,0,4,0,180 !SPHERE, RAD1, RAD2, THETA1, THETA2
alls
vsbw,all
wprota,,,90
vsbw,all
csys,4
vsel,s,loc,z,0,4
vsel,r,loc,x,0,4
vsel,inve
vdele,all,,,1
csys,1
wpcsys,-1
allsel
vglue,all
numcmp,all
kwpave,4
wprota,,90,0
vsbw,all
kwpave,162
wprota,,,90
vsbw,all
wprota,,-90,0
wpoff,,,0.3
vsbw,all
kwpave,226
csys,4
wprota,-90
vsel,s,loc,z,-0.100,0
VATT,1,,1
esize,0.01
VMESH,all !从这里往上的部分基本与原命令流相同 所以底部中间的网格密度没有更改
lsel,s,,,910
lsel,a,,,924
lsel,a,,,928
lsel,a,,,930
lesiez,all,0.03,,,10
allsel
mshape,0
vsweep,122
lsel,s,,,441
lsel,a,,,556
lsel,a,,,908
lsel,a,,,915
lesize,all,,,50,5
allsel
lsel,s,,,442
lsel,a,,,920
lesize,all,,,50,0.2
lsel,s,,,894
lsel,a,,,895
lesize,all,,,50,5
allsel
vsel,s,loc,z,0,0.3
vsel,u,,,113
vsel,u,,,122
mshape,0
vsweep,all
allsel
vsel,s,,,113
mshape,0
mshkey,1
vmesh,all
allsel
vsel,s,loc,z,-0.1,0.3
vsel,inve
lsel,s,,,558
lsel,a,,,919
lsel,a,,,913
lsel,a,,,925
lesize,all,,,50,0.5
mshape,0
mshkey,1
vmesh,120
vmesh,117
vmesh,118
vmesh,112
allsel
save
|