thdong 发表于 2007-8-20 10:26

曲面施加压强需注意哪些坐标系啊?

空心球体内外均加压强,请问各位大哥,施加压强时应该是加在单元的法向对吧,这里面需要考虑的问题很多,建模时总体坐标系、划分网格时单元坐标系的选取,还有单元的节点坐标系怎么设定呢?
请各位指教!

[ 本帖最后由 sysh320 于 2007-8-20 10:30 编辑 ]

thdong 发表于 2007-8-20 10:29

模型很简单,空心球外径2cm,内径1.8cm,模量2e11,泊松比0.3,外压5个气压,内压一个气压。我用的是solid45单元。
以上涉及很多ansys的基本概念,模型变成椭球体呢,单元坐标系等又如何考虑?

sysh320 发表于 2007-8-20 10:35

我感觉你是自己把问题搞复杂了,首先应该搞清楚各个坐标系的用途。
总体坐标系:用来确定几何图形参数的空间位置;
局部坐标系:用户自定义的坐标系,用来建模;
节点坐标系:定义每个节点的自由度方向和节点内力与位移的方向;
单元坐标系:确定材料特性主轴和单元内力与位移的方向;
显示坐标系:用于几何形状参数的列表和显示;
结果坐标系:可以在post1中将节点或单元结果转换到另一个坐标系中

sysh320 发表于 2007-8-20 10:40

根据你给定的材料参数,你用的是钢材,一般的材料属性就可以了,没有必要改变单元坐标系,除非你这种材料很特别。
还有就是节点坐标系跟加载也没有什么必然的关系,及时你这个模型很复杂。转换节点坐标系只有在加约束或载荷时需要,比如加扭矩时。
所加压强,无非就是面荷载(pressure),这个单元说明里都有的,就是那个几何图(gemotry),不同面有不同的加载序列号,是根据右手螺旋定义的

thdong 发表于 2007-8-21 00:22

谢谢sysh320的解答,我确实被这个空心球的静力分析搞糊涂了,模型简单,载荷也简单,又是线性情况下的静力分析,惟独分析出来的变形结果违反常理。
仁兄可以试着做做看,变形的结果不是想象中的球对称了。

thdong 发表于 2007-8-21 01:00

以下是模型建立和划分网格的命令流,加载后就可以算算看(外压1e6帕,内压1e5帕),看看各位仁兄分析出来的变形结果会不会违反常理。
欢迎一起探讨!

/PREP7
ET,1,SOLID45
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2e11   
MPDATA,PRXY,1,,0.3
CSWPLA,11,2,1,1,
CSWPLA,11,2,1,1,
SPH4, , ,0.02,0.019
CSYS,0
SMRT,6
SMRT,3
MSHAPE,1,3D
MSHKEY,0
CM,_Y,VOLU
VSEL, , , ,       1
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y
VMESH,_Y1   
CMDELE,_Y   
CMDELE,_Y1
CMDELE,_Y2
!*

venture 发表于 2007-8-22 11:27

本帖最后由 wdhd 于 2016-3-9 13:33 编辑

原帖由 thdong 于 2007-8-21 01:00 发表
以下是模型建立和划分网格的命令流,加载后就可以算算看(外压1e6帕,内压1e5帕),看看各位仁兄分析出来的变形结果会不会违反常理。
欢迎一起探讨!

/PREP7
ET,1,SOLID45
MPTEMP,,,,,,,,
MPTEMP,1, ...
你认为常理的变形是什么样的?

venture 发表于 2007-8-22 11:43

看起来也还正常。

thdong 发表于 2007-8-22 17:33

谢谢venture百忙中的计算,我把节点坐标系改成沿着曲面的法向后,也得出合理的结果了,以前球体是向一侧有刚体位移的。
但是现在我把模型改成椭球体用同样的方法计算,又出现向一侧的刚体位移了。模型中椭球长轴2厘米,短轴1厘米,其它条件不变。命令流如下。
/PREP7
ET,1,SOLID95
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2e11   
MPDATA,PRXY,1,,0.3
CSWPLA,11,1,0.5,1,
K,1,-0.02,0,0,
K,2,0.02,0,0,   
K,3,-0.019,0,0,
K,4,0.019,0,0,
L,       1,       2
L,       3,       4
LSTR,       1,       3
LSTR,       4,       2
FLST,2,4,4
FITEM,2,1   
FITEM,2,2   
FITEM,2,3   
FITEM,2,4   
AL,P51X
FLST,2,1,5,ORDE,1   
FITEM,2,1   
FLST,8,2,3
FITEM,8,1   
FITEM,8,2   
VROTAT,P51X, , , , , ,P51X, ,90, ,

SMRT,6
SMRT,3
MSHAPE,1,3D
MSHKEY,0
CM,_Y,VOLU
VSEL, , , ,       1
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y
VMESH,_Y1   
CMDELE,_Y   
CMDELE,_Y1
CMDELE,_Y2

/UI,MESH,OFF
FLST,2,2,5,ORDE,2   
FITEM,2,2   
FITEM,2,-3
NORA,P51X,1

CSYS,0,
FINISH
/SOL
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,4   
DA,P51X,SYMM
FLST,2,1,5,ORDE,1   
FITEM,2,2   
/GO
SFA,P51X,1,PRES,1e6
FLST,2,1,5,ORDE,1   
FITEM,2,3   
/GO
SFA,P51X,1,PRES,100000

CSYS,0
ALLSEL,ALL
SOLVE

thdong 发表于 2007-8-22 17:53

thdong 发表于 2007-8-22 17:54

thdong 发表于 2007-8-22 18:02

各位高手帮我分析一下哈,四分之一椭球变形后,上面倒数第2张是应力图,倒数第一张是X方向的位移图,位移图全是红色,表明整个模型沿X方向位移一致,即发生刚体位移,有违常理呢。

venture 发表于 2007-8-23 08:38

约束不对,加虚拟约束,或用1/8模型3个对称面都加对称约束。

thdong 发表于 2007-8-23 13:56

谢谢venture的解答,用八分之一模型计算的结果是对的,太谢谢了。
另外您说的虚拟约束是什么啊,如何用?
页: [1]
查看完整版本: 曲面施加压强需注意哪些坐标系啊?