- 积分
- 1
- 注册时间
- 2011-4-24
- 仿真币
-
- 最后登录
- 1970-1-1
|
悬赏10仿真币已解决
本帖最后由 清澈 于 2013-5-14 17:38 编辑
现在想通过数值模拟来验证一个解析解,涉及到的模型很简单,就是一个球体,在圆心处有个球形的空洞。
因为解析解要求任意方向都是均匀的,所以要求数值模拟中网格划分也必须沿任意一条直径,网格也对称。。
我开始想用一段命令流来尝试做个八分之一球体,想法是先画节点,再连单元。命令流如下:
finish
/clear
or=1 !设定球外径
ir=0.1 !设定球内径
n=20 !设定沿径向的单元分段数目
m=15 !设定沿环向的单元分段数目
z=15 !设定沿周向的单元分段数目
b=1.05 !径向放大因子
/prep7
et,1,185
csys,2
*do,p,0,z
*do,i,0,m
*do,j,0,n
n,,ir+(or-ir)/n*(1-b**j)/(1-b),360.0/m*i,180.0/z*p
*enddo
*enddo
*enddo !至此画完节点
*do,s,1,z-1 !沿周向通过循环生成单元
*do,q,0,m-1 !沿环向通过循环生成单元
inc=q*(n+1)
e,s*(m+1)*(n+1)+1+inc,s*(m+1)*(n+1)+2+inc,s*(m+1)*(n+1)+n+3+inc,s*(m+1)*(n+1)+n+2+inc,1+inc+(s+1)*(m+1)*(n+1),2+inc+(s+1)*(m+1)*(n+1),n+3+inc+(s+1)*(m+1)*(n+1),n+2+inc+(s+1)*(m+1)*(n+1)
*repeat,n,1,1,1,1,1,1,1,1
*enddo
*enddo !至此连成单元
但是这样的话,会出错,说后面的单元雅克比矩阵怎么地怎么地,如果在周向也设置与径向类似的缩小因子的话(这是我瞎琢磨的,就感觉后来越划分网格上面的径向与周向单元边长差的越多,),把n,,ir+(or-ir)/n*(1-b**j)/(1-b),360.0/m*i,180.0/z*p改成n,,ir+(or-ir)/n*(1-b**j)/(1-b),90.0/m*i,90.0/z*(1-(1/b)**p)/(1-(1/b)),这样画出来的倒是么有错误,也能连成单元。
问题是这样连不全呀,也就是周向本来设置的从0°到90°,根本就到不了90°就停了。。。差了一个角~~~如图。。。
这说明这样建立单元是有问题的吧。。。。
在缺的那个角那,即使能补上也不可能跟另外两个角处的网格对称哇。。。。。
所以我想请教,这个问题是出在哪里,最主要的,空心球体可以按照分别在径向R、环向θ和周向φ分别划分多少份,就可以连成完全对称的六面体吗?
谢谢大家啦!!
|
附件: 您需要 登录 才可以下载或查看,没有账号?注册
|