壳体所围体积怎么计算？

w1h2y3

用shell41单元（设置了壳体厚度）模拟一空心球壳在外力作用下的变形，怎么计算变形后球壳所围的体积？多谢指教。
另外，壳体设置厚度以后为什么ansys不能识别它是volume？

venture

壳体设置厚度以后为什么ansys不能识别它是volume，壳体设置的厚度它还是壳体，还是壳单元属性。

w1h2y3

也就是说壳体设置厚度以后ansys还是识别其为面，而不是体？

另外，空心球壳所围的体积应该怎么计算呢？

93758902

我想法是这样子的供你参考：
计算完变形后，通过upgeom命令修改Factor=1，原有模型更新的为变形后的模型位子，然后再做一个大的球体或任意实体包围着这个壳体（比这个壳体大就行），通过partition命令把大实体分为壳内部的部分和壳外部的部分，再通过Vsum命令求出壳内实体的体积，也就求出了壳所包围的体积

[ 本帖最后由 93758902 于 2008-10-20 12:46 编辑 ]

16443

这个方法可以试试。

w1h2y3

非常谢谢你的建议。我尝试了一下，用upgeom处理过以后只是改变了node和element的值，所以用plot->nodes或者elements能得到变形后的模型，但是用plot->areas或者lines都还是变形前的模型，所以不能用partition命令处理，而且partition只能处理两个volume之间的重叠，但是对于壳体，ansys只是把它识别为area而不是volume，不知道该如何处理？

93758902

你建的这个壳应该是个封闭的壳，下面我的两个想法都是基于upgeom命令更新后的node coordinate，
1.你可以先把Node和element信息存起来，再创建一个新的项目加载上面的信息，如果原来Node编号不合理，可以先合并并且压缩所有编号形成合理的编号，然后使用KNODE命令循环地在原有Node上形成Kp，再通过A命令多选几个形成光滑的面，用循环*do形成一个多个包围的面，将所有的面相加就形成一个面，然后你在建个体包围着这个壳体，在通过divide命令（divide不行再找一个别的命令，总有一个行）就把体分成壳内和壳外部分，再用vsum一算内部的体积。
2.在壳中间创建一个node，用E命令，与壳体上的node形成四面体单元，然后求所有单元的体积和。如果单元Jacobi太大，可能算起来的体积不准。不过你可以建一个比较大圆壳，划分圆壳上的单元，圆壳上的节点连接外壳的节点形成四面体单元，然后求所有单元的体积。
希望对你有用

w1h2y3

非常谢谢你的建议，我试试。但是这样由node集合组成的面所围成的体积可能很大程度上取决于面的光滑度，可能会对最后的体积计算带来误差。

[ 本帖最后由 w1h2y3 于 2008-10-20 16:34 编辑 ]

93758902

多少会有误差的把，不知道是否有好点办法

[ 本帖最后由 93758902 于 2008-10-20 18:04 编辑 ]

cxw_luck

要是不复杂的话，你重建个模型，用实体单元做一下行不

w1h2y3

因为问题本身是空心壳体在内压和外力作用下变形的问题，如果用实体单元做会不会凸现不出薄壳效应呢？

93758902

有薄壳和中厚壳单元，也有相对自由度壳单元，你要你选对了单元，薄壳的薄膜效应本身就考虑了，也不会发生薄膜锁死