关于壳单元的命令流

lvgang637 · 发表于 2008-9-4 09:58

请问：1、板的厚度是通过建模见出来的还是通过赋予实常数来实现的？
2、对板进行网格划分，本文中的方式划分的网格过多，致使在进行静力分析时，导致速度太慢，有没有更好的划分网格的方法？
3、我想对这个板壳进行优化设计，以板的厚度为设计变量，以承受的最大应力和最大挠度为状态变量，以质量为目标函数，具体的命令流时什么，能不能发一些，供参考。

/FILNAM,SCRATCH
/TITLE,OPTIMIZATION OF TOWER
/PREP7
ET,1,SHELL93
MP,EX,1,2.07e11
MP,PRXY,1,0.3
MP,DENS,1,7850
R,1,0.01
*SET,D1,5.663
*SET,D2,4.593
*SET,D3,4.593
*SET,D4,4.243
*SET,D5,4.243
*SET,t1,0.018
*SET,t2,0.018
*SET,t3,0.016
*SET,t4,0.016
*SET,t5,0.014
K,1,D1/2,,,
K,2,D2/2,,20.5,
K,3,D3/2,,41,
K,4,D4/2,,61.5,
K,5,D5/2,,82,
K,6,D5/2+t5,,82,
K,7,D4/2+t4,,61.5,
K,8,D3/2+t3,,41,
K,9,D2/2+t2,,20.5,
K,10,D1/2+t1,,0,
l,1,2
l,2,3
l,3,4
l,4,5
l,5,6
l,6,7
l,7,8
l,8,9
l,9,10
l,10,1
LPLOT
al,all
K,100,0,,0,
K,101,0,,82,
vrotat,all,,,,,,100,101
AATT,1,1,1,0
LESIZE,ALL,0.25
MSHAPE,1,2D
AMESH,ALL
CSYS,0 !创建刚性区域
K,9999,0,0,82
ET,2,MASS21
R,1,0.01
KSEL,S,LOC,Z,82,82
KMESH,ALL
NSEL,S,LOC,Z,82,82
CERIG,NODE(0,0,82),ALL,ALL,,,,
/SOL
ANTYPE,0
DA,11,ALL,
DA,22,ALL,
DA,33,ALL,
DA,44,ALL,
ACEL,0,0,9.8, !定义重力加速度
F,164574,FY,385732
nsel,all
/STATUS,SOLU
SOLVE !求解
FINISH

复制代码

[ 本帖最后由 16443 于 2008-9-4 15:06 编辑 ]

16443 · 发表于 2008-9-4 15:07

1。用R参数设置就可以。
2。放大lesize就能减少网格数量。

graduates · 发表于 2008-11-12 22:35

板厚是实常数定义出来的！
虽然在模型中看不出来，实际上它是有厚度的，可以用命令/eshape,1，再/replot，你就可以看到厚度了。

93758902 · 发表于 2008-11-12 23:29

原帖由 16443 于 2008-9-4 15:07 发表

 登录/注册后可看大图

1。用R参数设置就可以。
2。放大lesize就能减少网格数量。

厚度还可以通过定义一个数组然后通过RTHINK命令赋予变厚度，这种方法是直接可以看到厚度的，不需要/eshape,1。你这个网格完全可以画成全部四面型的，计算量也少。质量可以通过单元的体积相加乘以密度得到。

[ 本帖最后由 93758902 于 2008-11-12 23:37 编辑 ]

wrmaster · 发表于 2008-11-13 18:56

目前暂时看不懂，不过很有兴趣

账号		自动登录	找回密码
密码			我要加入

[结构分析] 关于壳单元的命令流