渊源1983 发表于 2006-5-31 13:53

大家看看这个网格是如何划分出来的

第一个图是一本博士学位论文上有限元划分的图
第二个图是我划分的（计算不收敛）
想请教下大家
第一个图的网格是如何划分出来的？
我认为的难点就是圆柱不能Mapped划分，要划分出这种六面体只能用sweep
但是圆柱和长方体不能一起sweep，只能各自sweep
但是各自sweep的效果，明显不好，见图三
所以只有在接触处用free 然后进行单元退化
然而生成的四面体较多计算时不收敛
如何能划分出如图1这样规则而且耦合的又好的网格出来呢？
希望大家畅所欲言 ～～
小弟不胜感激

xiaozhu 发表于 2006-5-31 17:02

试着把上面的小圆柱分割下来
规则的体就可以用mapped了

渊源1983 发表于 2006-5-31 17:16

但是做起来真的很难啊 ～～
怎么分割呢？

渊源1983 发表于 2006-6-1 11:08

长方体里面还有钢筋
不能divide
哎 ～～
用四面体单元吧
不收敛
快疯了

渊源1983 发表于 2006-6-1 16:35

/PREP7
et,1,link8 ！钢筋单元
et,2,solid65 ！承台体单元
et,3,solid95 ！过渡单元
et,4,solid45 ！加载区分散荷载用
r,1,113 ！12号钢筋
r,2
KEYOPT,2,1,1
KEYOPT,2,3,0
KEYOPT,2,5,0
KEYOPT,2,6,0
KEYOPT,2,7,1
KEYOPT,2,8,0
mp,ex,1,2.1e5 ！定义钢筋属性
mp,prxy,1,0.3
tb,bkin,1
tbdata,,335
mp,ex,2,3e4 ！定义混凝土承台体属性
mp,prxy,2,0.2
*SET,fc,20.97
*SET,ft,2.2
TB,CONC,2,
TBTEMP,0
TBDATA,,0.5,0.9,ft,-1,,
TBDATA,,,,,,,
tb,miso,2
tbtemp,0
tbpt,,0,0
tbpt,,0.0002,0.0002*3e4
tbpt,,0.0004,0.36*fc
tbpt,,0.0006,0.51*fc
tbpt,,0.0008,0.64*fc
tbpt,,0.001,0.75*fc
tbpt,,0.0012,0.84*fc
tbpt,,0.0014,0.91*fc
tbpt,,0.0016,0.96*fc
tbpt,,0.0018,0.99*fc
tbpt,,0.002,fc
tbpt,,0.0025,0.95*fc
tbpt,,0.003,0.9*fc
tbpt,,0.0035,0.85*fc
tbpt,,0.004,0.8*fc
mp,ex,3,3.65e4 ！ 定义桩、柱混凝土属性
mp,prxy,3,0.2
mp,ex,4,3.3e4 ！定义加载区单元属性
mp,prxy,4,0.2

n,1,50,50,50 ！建钢筋单元
NGEN,15,1,1 , , ,50, , ,1,
NGEN,9,15,1, 15, 1, ,50, ,1,
type,1
mat,1
r,1
*do,ii,1,14,1
e,ii,ii+1
*enddo
FLST,4,14,2,ORDE,2
FITEM,4,1
FITEM,4,-14
EGEN,9,15,P51X, , , , , , , , ,50, ,
*do,ii,1,106,15
e,ii,ii+15
*enddo
FLST,5,8,2,ORDE,2
FITEM,5,127
FITEM,5,-134
ENGEN,8,15,1,P51X, , , , , , , ,50, , ,

blc4,,,800,500,400 ！建混凝土承台体模型
wpoff,,,50
vsbw,all
vglue,all
wpcsys,-1,0
wpoff,300,150,400
blc4,,,200,200,150 ！建柱体
wpcsys,-1,0
wpoff,100,100
CYL4, , ,50, , , ,-80 ！建桩体
VGEN,2,4, , , ,300, , ,0
VGEN,3,4, 5, ,300, , , ,0
allsel
vsel,s,loc,z,500,550 ！加载区划分单元
aslv,a
lsla,s
lesize,all,50
vatt,4,2,4
vmesh,all
allsel
vsel,s,loc,z,400,500 ！柱体划分单元
aslv,a
lsla,s
lesize,all,50
vatt,3,2,2
vmesh,all
allsel
vsel,s,loc,z,50,400 ！承台体划分单元
aslv,a
lsla,s
lesize,all,50
vatt,2,2,2
vmesh,all
allsel
vsel,s,loc,z,-80,50 ！划分桩及其桩承台接触层
SMRT,7
mshape,1,3d
MSHKEY,0
vatt,3,2,3
vmesh,all
allsel,all
tchg,95,92,0
allsel

asel,s,loc,z,-80 ！桩底施加约束
nsla,s
d,all,all
FINISH
/solu
ANTYPE,0 ！非线性设置
CNVTOL,F, ,0.05,2,0.01,
outres,all,all
autots,on
lnsrch,on
arclen,off
ncnv,2
neqit,20
pred,on

allsel ！施加荷载
nsubst,20
allsel
nsel,s,loc,z,550
f,all,fz,-16000
kbc,0
nsel,all
lswrite,1
allsel
fcum,add,1
nsubst,20
allsel
nsel,s,loc,z,550
f,all,fz,-16000
kbc,0
nsel,all
lswrite,2
allsel
fcum,add,1
nsubst,20
allsel
nsel,s,loc,z,550
f,all,fz,-8000
kbc,0
nsel,all
lswrite,3
allsel
fcum,add,1
nsubst,20
allsel
nsel,s,loc,z,550
f,all,fz,-8000
kbc,0
nsel,all
lswrite,4
allsel
fcum,add,1
nsubst,20
allsel
nsel,s,loc,z,550
f,all,fz,-4000
kbc,0
nsel,all
lswrite,5
save

linqus 发表于 2006-6-1 17:20

看看这个网格是如何划分出来

这个是用abq做的，
picture1：

linqus 发表于 2006-6-1 17:22

回复：（linqus）看看这个网格是如何划分出来

picture2：
一半模型

linqus 发表于 2006-6-1 17:38

回复：（linqus）回复：（linqus）看看这个网格是如何...

picture3
底部：

共32976个单元

渊源1983 发表于 2006-6-1 19:30

看起来比较不错哦
不过好多单元数
尺寸取的比我的要小很多 ～～

usefully 发表于 2006-6-1 20:28

这个网格划的真不错，希望作者能把划分的方法奉献一下，让大家也学习一下。谢谢！

渊源1983 发表于 2006-6-2 10:15

。。。。
他是用abq做的

linqus 发表于 2006-6-3 09:16

用adina自由划分做的。

尝试着用adina做一个，自由划分，
效果还行：

pic1：9792elements

linqus 发表于 2006-6-3 09:17

用adina自由划分做的。

pic2：
22464elements

linqus 发表于 2006-6-3 09:20

用adina自由划分做的。

由于没有加以控制，
所以划分时间较abq长一些，
还有，底部网格与abq划分不同。
手头没有ansys，没做。呵呵

底部：

渊源1983 发表于 2006-6-4 09:47

回复：（渊源1983）大家看看这个网格是如何划分出来的...

强人！！
不过这两种软件都不熟悉
也没时间去熟悉了
就想用ansys做下 ～～
页: [1] 2 3
查看完整版本: 大家看看这个网格是如何划分出来的