该怎么进行网格划分

alwayszi

!!!!!1, fr,2,copper,3,plcc,4,solder!!!!!!!!
finish
/clear
/file, pcb thermalmechanical analysis
/title, pcb thermalmechanical analysis
/prep7
et,1,70
lpcb=200
wpcb=180
uh=0.2
wplcc=30
thplcc=1
thsd=0.5
l1=85
l2=45
w1=45
hair=10
tair=25
hng=30
t1=230
t2=180
t3=260
t4=25
mp,dens,1,1870e-9                                               
mp,kxx,1,0.25
mptemp,1,30,70,100,130,220
mpdata,c,1,1,1200e6,1380e6,1452e6,1512.5e6,1625e6
mp,dens,2,8933e-9                                                  
mp,kxx,2,394
mp,c,2,390e6
mp,dens,3,1800e-9                                                  
mp,kxx,3,18
mptemp,1,40,80,160,170,220,225,230,250
mpdata,c,3,1,740e6,850e6,880e6,960e6,970e6,1050e6,920e6,900e6
mp,dens,4,8350e-9
mp,kxx,4,50.9
mp,c,4,0.15e6

   
block,0,lpcb,0,wpcb,0,0.5*uh                                       
block,0,lpcb,0,wpcb,0.5*uh,1.5*uh
block,0,lpcb,0,wpcb,1.5*uh,2.5*uh
block,0,lpcb,0,wpcb,2.5*uh,3.5*uh
block,0,lpcb,0,wpcb,3.5*uh,4.5*uh
block,0,lpcb,0,wpcb,4.5*uh,5.5*uh
block,0,lpcb,0,wpcb,5.5*uh,6.5*uh
block,0,lpcb,0,wpcb,6.5*uh,7*uh
/pnum,volu,1
vplot

block,l1,l1+wplcc,w1,w1+wplcc,0,7*uh+thsd
block,l2,l2+wplcc,wpcb-w1-wplcc,wpcb-w1,0,7*uh+thsd
block,lpcb-l2-wplcc,lpcb-l2,wpcb-w1-wplcc,wpcb-w1,0,7*uh+thsd

block,l1,l1+wplcc,w1,w1+wplcc,7*uh+thsd,7*uh+thsd+thplcc
block,l2,l2+wplcc,wpcb-w1-wplcc,wpcb-w1,7*uh+thsd,7*uh+thsd+thplcc
block,lpcb-l2-wplcc,lpcb-l2,wpcb-w1-wplcc,wpcb-w1,7*uh+thsd,7*uh+thsd+thplcc
vsel,all
vovlap,all
vglue,all
   
vsel,all
wprot,0,90
wpoffs,0,0,-w1
vsbw,all
wpave
wpoffs,0,0,-w1-wplcc
vsbw,all
wpave
wpoffs,0,0,-(wpcb-w1-wplcc)
vsbw,all
wpave
wpoffs,0,0,-(wpcb-w1)
vsbw,all
wpave
wpcsys,0

wprot,0,0,90
wpoffs,0,0,l2
vsbw,all
wpoffs,0,0,wplcc
vsbw,all
wpoffs,0,0,l1-l2-wplcc
vsbw,all
wpoffs,0,0,wplcc
vsbw,all
wpave
wpoffs,0,0,lpcb-l2
vsbw,all
wpoffs,0,0,-wplcc
vsbw,all
wpave
wpcsys
vsel,s,loc,z,0,0.5*uh
vsel,a,loc,z,1.5*uh,2.5*uh
vsel,a,loc,z,4.5*uh,5.5*uh
vsel,a,loc,z,6.5*uh,7*uh
vatt,2,,1
allsel
vsel,s,loc,z,0.5*uh,1.5*uh
vsel,a,loc,z,2.5*uh,4.5*uh
vsel,a,loc,z,5.5*uh,6.5*uh
vatt,1,,1
allsel
vsel,s,loc,z,7*uh,7*uh+thsd
vatt,4,,1
allsel
vsel,s,loc,z,7*uh+thsd,7*uh+thsd+thplcc
vatt,3,,1
allsel
lesize,all,10
vmesh,all
mshape,0
mshkey,1
这么分网不行。后来单独对体16、19、22和体17、20、23指定网格大小进行扫掠，这六个体网格划分没问题，对剩下的体进行网格划分时显示“单元尺寸violate单元尺寸限制”，大家帮我看看该怎么分网比较好呢？

ldhltw

单元尺寸限制可能是有尖角存在，单元分的细点试试

alwayszi

ldhltw 发表于 2012-2-21 11:19
单元尺寸限制可能是有尖角存在，单元分的细点试试

*** WARNING ***   SUPPRESSED MESSAGE    CP =     225.297   TIME= 11:29:28
Brick element 56 has an aspect ratio of 50, which exceeds the warning   
limit of 20.
单元size改成4，划分出来的单元全显示上买你的错误。用solid单元，厚度方向是不是不能值划分一层，也就只有一层单元？对于薄板，是不是应该改用壳单元呢？谢谢您的回复

senja

我看了一下，这是层合材料吗？   楼主控制的是分数，不是单元尺寸。由于层间间距很小，但面比较大，统一控制份数会导致单元形状不好，即单元边长比太大。

PS：楼主是需要什么样的网格？厚度为多个单元嘛？  solid单元的厚度方向的单元数并没有什么限制，只是根据计算要求认为给定的，也就是单元尺寸。

senja

划分了一次，不知道能不能满足楼主的要求。若是需要单层厚度方向的单元数大于1个，那么需要把单元尺寸设置（esize）更小一些，这样才能实现。要不然还是会出现单元形状超出软件要求的警告。



finish
/clear
/file, pcb thermalmechanical analysis
/title, pcb thermalmechanical analysis

/prep7


lpcb=200  
wpcb=180
uh=0.2
wplcc=30
thplcc=1
thsd=0.5
l1=85
l2=45
w1=45
hair=10
tair=25
hng=30
t1=230
t2=180
t3=260
t4=25

et,1,70

mp,dens,1,1870e-9                                               
mp,kxx,1,0.25
mptemp,1,30,70,100,130,220
mpdata,c,1,1,1200e6,1380e6,1452e6,1512.5e6,1625e6

mp,dens,2,8933e-9                                                  
mp,kxx,2,394
mp,c,2,390e6
mp,dens,3,1800e-9                                                  
mp,kxx,3,18


mptemp,1,40,80,160,170,220,225,230,250
mpdata,c,3,1,740e6,850e6,880e6,960e6,970e6,1050e6,920e6,900e6
mp,dens,4,8350e-9
mp,kxx,4,50.9
mp,c,4,0.15e6

   
block,0,lpcb,0,wpcb,0,0.5*uh                                       
block,0,lpcb,0,wpcb,0.5*uh,1.5*uh
block,0,lpcb,0,wpcb,1.5*uh,2.5*uh
block,0,lpcb,0,wpcb,2.5*uh,3.5*uh
block,0,lpcb,0,wpcb,3.5*uh,4.5*uh
block,0,lpcb,0,wpcb,4.5*uh,5.5*uh
block,0,lpcb,0,wpcb,5.5*uh,6.5*uh
block,0,lpcb,0,wpcb,6.5*uh,7*uh


block,l1,l1+wplcc,w1,w1+wplcc,0,7*uh+thsd
block,l2,l2+wplcc,wpcb-w1-wplcc,wpcb-w1,0,7*uh+thsd
block,lpcb-l2-wplcc,lpcb-l2,wpcb-w1-wplcc,wpcb-w1,0,7*uh+thsd

block,l1,l1+wplcc,w1,w1+wplcc,7*uh+thsd,7*uh+thsd+thplcc
block,l2,l2+wplcc,wpcb-w1-wplcc,wpcb-w1,7*uh+thsd,7*uh+thsd+thplcc
block,lpcb-l2-wplcc,lpcb-l2,wpcb-w1-wplcc,wpcb-w1,7*uh+thsd,7*uh+thsd+thplcc
vsel,all
vovlap,all
vglue,all

vsel,all
wprot,0,90
wpoffs,0,0,-wl
vsbw,all
wpave
wpoffs,0,0,-w1-wplcc
vsbw,all
wpave
wpoffs,0,0,-(wpcb-w1-wplcc)
vsbw,all
wpave
wpoffs,0,0,-(wpcb-w1)
vsbw,all
wpave
wpcsys,0

wprot,0,0,90
wpoffs,0,0,l2
vsbw,all
wpoffs,0,0,wplcc
vsbw,all
wpoffs,0,0,l1-l2-wplcc
vsbw,all
wpoffs,0,0,wplcc
vsbw,all
wpave
wpoffs,0,0,lpcb-l2
vsbw,all
wpoffs,0,0,-wplcc
vsbw,all
wpave
wpcsys

allsel
wprota,,-90
wpoffs,85,,45
vsbw,all
allsel
wpcsys,-1,0

vsel,s,loc,z,0,0.5*uh   
vsel,a,loc,z,1.5*uh,2.5*uh
vsel,a,loc,z,4.5*uh,5.5*uh
vsel,a,loc,z,6.5*uh,7*uh
vatt,2,,1
allsel

vsel,s,loc,z,0.5*uh,1.5*uh
vsel,a,loc,z,2.5*uh,4.5*uh
vsel,a,loc,z,5.5*uh,6.5*uh
vatt,1,,1
allsel

vsel,s,loc,z,7*uh,7*uh+thsd
vatt,4,,1
allsel

vsel,s,loc,z,7*uh+thsd,7*uh+thsd+thplcc
vatt,3,,1
allsel

mshape,0
mshkey,1

esize,2

vmesh,all

红色部分为动过手脚的。
  

厚度方向测试图：

alwayszi

senja 发表于 2012-2-21 21:18
我看了一下，这是层合材料吗？   楼主控制的是分数，不是单元尺寸。由于层间间距很小，但面比较大，统一控 ...

是符合材料，采用的统一单元solid70，先做温度场分析。
你分析的太到位了，就是厚度方向维度给长宽维度不是一个数量级导致划分出网格的宽总比超过系统默认的20.把网格划分的尺寸变小就可以了。嘿嘿
我现在遇见了，瞬态热分析不收敛的问题，附命令流，您帮忙看看可能是什么原因
/solu
antype,trans
kbc,0
tref,30
!!!!!!!!!!!!!!!!!!!!warming!!!!!!!!!!!!!!!!!!!!!!
deltim,6  
tim=0   
timint,on
*do,i,1,10
tim=tim+6
time,tim
vsel,all
asel,r,ext,,,,,1
nsel,r,loc,z,0
nsel,a,loc,z,7*uh
nsel,a,loc,z,7*uh+thsd+thplcc
nsel,r,loc,x,0,20*i
sf,all,conv,hng,t1
vsel,all
asel,r,ext,,,,,1
nsel,r,loc,z,0
nsel,a,loc,z,7*uh
nsel,a,loc,z,7*uh+thsd+thplcc
nsel,r,loc,x,20*i,200
sf,all,conv,hair,tair
allsel
solve
sfdele,all,conv
*enddo
deltim,6  
tim=60   
timint,on
*do,i,11,15
tim=tim+6
time,tim
vsel,all
asel,r,ext,,,,,1
nsel,r,loc,z,0
nsel,a,loc,z,7*uh
nsel,a,loc,z,7*uh+thsd+thplcc
nsel,r,loc,x,0,200
sf,all,conv,hng,t1
allsel
solve
sfdele,all,conv
*enddo

求解中从graphicl solution tracking图中明显看到计算结果不收敛。后面再做热应力分析就受阻了。you are so sweet and helpful，thanks

senja

这个我没运行，也真说不出什么东西。
不过，一般瞬态分析不收敛的因素楼主可以去查查，如子步过大等等。再到错误文件里面看看，是否能够得到些有用的信息。