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

thdong

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

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

thdong

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

sysh320

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

sysh320

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

thdong

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

thdong

以下是模型建立和划分网格的命令流，加载后就可以算算看（外压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

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

/PREP7  
ET,1,SOLID45
MPTEMP,,,,,,,,  
MPTEMP,1, ...

你认为常理的变形是什么样的？

venture

看起来也还正常。

thdong

谢谢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

thdong

thdong

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

venture

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

thdong

谢谢venture的解答，用八分之一模型计算的结果是对的，太谢谢了。
另外您说的虚拟约束是什么啊，如何用？