生成变厚度壳时的疑问

bishenghjy

一个10×10的矩形面划分20×20的shell单元网格模型，其厚度为thick(x,y)=0.5+0.2*x+0.02*y**2!初始化ansys环境
finsih
/clear
/filname,shell-thick
!创建壳体网格模型
/prep7  
et,1,63
rect,,10,,10
eshape,2
esize,,20   +
amesh,1
eplot   
!变厚度壳的实现
*SET,mxnode,ndinqr(0,14)
!提取节点总数mxnode
*dim,thick,,mxnode  
!定义厚度数组thick(mxnode)
*do,node,1,mxnode   
*if,ndinqr(node,1),eq,1,then
*SET,thick(node),0.5+0.2*nx(node)+0.02*ny(node)**2  
*else   
*SET,thick(node),0  
*endif  
*enddo  
!根据厚度分布函数计算每个节点的厚度，并将每个节点厚度记录在厚度数组thick(mxnode)
*SET,node,  
*SET,mxnode,
!删除变量node,mxnode
rthick,thick(1),1,2,3,4
!分配单元厚度
/eshape,1   
eplot   
save


小弟初学APDL领会不了红体字的命令到底起了什么作用 望大侠指点

pfsr

1. ndinqr是查询函数
2. 定义壳在各个节点的厚度

bishenghjy

我知道ndinqr是查询函数 但是不知道它在这段程序力起了什么作用，如果不进行IF判断直接thick(node),0.5+0.2*nx(node)+0.02*ny(node)**2  回造成设么结果？  请指教

东湖楼主

*if,ndinqr(node,1),eq,1,then
确认该节点是否在选择集内。
等同于
*if,nsel(node),eq,1,then
RTHICK,THICK(1),1,2,3,4  Defines variable thickness at nodes for shell elements. 可以参见ansys的帮助。

bishenghjy

谢谢  第一句我明白了  可是RTHICK,THICK(1),1,2,3,4 究竟是怎么把厚度分配到节点的呢？能否将一下具体的分配过程，为什么THICK(1),的小括号里是1 是2，3，4不行吗  小括号里的1是什么意思？

[ 本帖最后由 bishenghjy 于 2008-6-7 21:44 编辑 ]

bishenghjy

我的理解是RTHICK,THICK(1),1,2,3,4将THICK(MXNODE)记录的厚度值分配给每个单元,THICK(1)是从数组的第一个值开始.
我的理解还很浅,希望高手再来指教.

pfsr

你给出的例子就是ANSYS HELP自带的。建议查看一下HELP里是怎么解释的。

holli

RTHICK, Par, ILOC, JLOC, KLOC, LLOC

Defines variable thickness at nodes for shell elements.

PREP7: Meshing

MP ME ST DY <> PR <> <> <> PP ED

Par
Array parameter (indexed by node number) that expresses the function to be mapped. For example, func (17) should be the desired shell thickness at node 17.

ILOC
Position in real constant set for thickness at node I of the element (default 1).

JLOC
Position in real constant set for thickness at node J of the element (default 2).

KLOC
Position in real constant set for thickness at node K of the element (default 3).

LLOC
Position in real constant set for thickness at node L of the element (default 4).

Notes
After RTHICK, each selected element will have its REAL number matching its ELEM number, for example, R (ILOC) = func (I NODE), R (JLOC) = func (J NODE), etc. Any other real constants on a previously defined real constant set will remain unchanged. This command cannot be used for beam elements.

holli

对RTHICK不明白的是，节点变量par与ILOC, JLOC, KLOC, LLOC的关系？

bishenghjy

自己顶一下   望高手赐教