【思雨凝风】讨论displacement，fixed和remote的区别与应力奇异

kevin666210

同意jackie版的说法，但是双pres的displacement不对吧，都跑中间去了，不符合要求的说。
icem这么强大了？完了，看来我还不够努力，稍微不努力就要被时代落下啊。。。
另外37#这个说法我仍然感到困惑。就算是数值 ...
LION_ARTHAS 发表于 2011-5-31 19:55

inertia relief

likui285

网格划分越密应力奇异越大，我想应该是这样解释：实际分析中施加在一个节点的载荷约束，不可能就施加在一个点上，而是很小的一个面上，随着划分越密，自然面越小，从而与之成反比的应力值越来越大了。不知是不是这样理解，前辈纠正啊！

LION_ARTHAS

非常欢迎不同意见，有分歧才有讨论，有讨论才有灵感
我是把载荷施加到了面上，是pressure，不是force，而且pressure和force结果是一样的
从ansys classic角度解释比较容易。实际分析的时候如果把载荷施加到点线面这种model的特征上，计算的时候ansys会自动把这些载荷转移到节点上，而且classic有一步前处理就是这个
仅代表个人看法，欢迎批评指正

kevin666210

inertia relief
kevin666210 发表于 2011-6-1 15:16

inertia relief：惯性释放
1、要将symetry删除，inertia relief才可以用。
2、随便约束任意节点（slice的目的）。
3、四分之一模型，约束对称面。
4、inertia relief不是对所有单元都适用。
5、建议查看GB/T228-2002

LION_ARTHAS

和高人擦肩而过,正在阅读GB/T228-2002
惯性释放百度百科的解释

通常我们做线性静力分析需要保证结构没有刚体位移，否则求解器没有办法计算。但是很多分析，例如飞机在飞行时,轮船在航行时，要想计算结构上的应力分布，需要采用惯性释放(inertia relief)，在结构上施加一个虚假的约束反力来保证结构上合力的平衡。 　　惯性释放是MSC．NASTRAN或ANSYS中的一个高级应用，允许对完全无约束的结构进行静力分析。简单地说就是用结构的惯性(质量)力来平衡外力。尽管结构没有约束，分析时仍假设其处于一种“静态”的平衡状态。采用惯性释放功能进行静力分析时，只需要对一个节点进行6个自由度的约束(虚支座)。针对该支座，程序首先计算在外力作用下每个节点在每个方向上的加速度，然后将加速度转化为惯性力反向施加到每个节点上，由此构造一个平衡的力系(支座反力等于零)。求解得到的位移描述所有节点相对于该支座的相对运动。 当结构上承受一系列载荷作用时，ANSYS提供的惯性释放的功能可计算结构与所承受载荷平衡所需要的加速度，此时结构上的位移约束只能约束刚体运动，且计算的约束点上的反力应为零。 因为加速度是根据单元质量矩阵和所施加载荷计算的，包括平移和旋转加速度，所以输入数据中一定要有材料密度。下面为一小例子： 当结构上承受一系列载荷作用时，ANSYS提供的惯性释放的功能可计算结构与所承受载荷平衡所需要的加速度，此时结构上的位移约束只能约束刚体运动，且计算的约束点上的反力应为零。 因为加速度是根据单元质量矩阵和所施加载荷计算的，包括平移和旋转加速度，所以输入数据中一定要有材料密度。下面为一小例子：

zyr

44# kevin666210
kevin666210:
    你好！首先感谢你提出了移除中间节点的办法（dropped），但我还是有些疑惑望解答：
1.此处为啥要用惯性释放？
2.产生应力梯度与采用1/4模型或全模型没关系吧？
另外，我按你的做法（惯性释放，采用1/4模型，dropped）应力的确无梯度，但若无dropped就有应力梯度；无惯性释放，采用全模型，有dropped的六面体网格无应力梯度；无惯性释放，采用全模型，无dropped的四面体网格也无应力梯度。见如下应力云图

以下为3种情况的WB12.1的文件

LION_ARTHAS

回楼上的
kevin版主是针对我举得双pressure的例子说的，请看45#的inertia relief的百度百科，是专门针对无约束model分析的

zyr

47# LION_ARTHAS
我看了“45#的inertia relief”，也看了前面的双pressure，双pressure不就是要“两边都跑中间去吗”（两端位移相反且数值大小相等）？双pressure的受力物体是静止的，不是运动的，不用惯性释放就可计算啊？
      我是看到“kevin666210”的网格划分中采用了移除中间节点（dropped），所以我在6面体网格中采用了dropped，这样计算的结果的确无应力梯度。
      产生应力梯度的原因，到底是网格质量高低的问题，还是ANSYS计算保留数值精度的问题呢？所以，还希望你再研究一下，找到产生应力梯度的真正原因。

LION_ARTHAS

48# zyr


纯属个人看法，欢迎批评指正
应力释放是属于ansys的高级应用，用于求解无约束model的应力应变分析，如运动中的飞机。
至于舍掉中节点就没有应力梯度的现象，个人觉得jackie指出的说的有道理，有限元软件得到的都是近似解，不是解析解，就有个误差问题。

另外向jackie版主求助，您曾在本帖指出如何在wb调用icem，我已经会了
但是这里有个问题，不知道是我哪里弄得不对，我画完网格关闭之后他会重读模型，然后又调用一遍ICEM，非常无语，求真相。非常感谢，祝您端午节快乐
祝大家端午节快乐

附赠个笑话
<
问：如果有一天，一个非常漂亮的女孩，犹豫了很久之后终于向你表白，你会怎么做？
跟帖：先存一下档
>
愿大家天天开心

722tzz

接触分析时的应力奇异又应该如何处理呢

LION_ARTHAS

楼上的朋友可不可以讲详细点，什么工况下，为什么会出现。谢谢

linuaries

inertia relief的问题，建议参考Nastran相关解释，ANSYS在这方面还比较弱。你看到过飞机或者汽车的整车会用ANSYS分析的么？

xia8444252

inertia relief的问题，建议参考Nastran相关解释，ANSYS在这方面还比较弱。你看到过飞机或者汽车的整车会用ANSYS分析的么？
linuaries 发表于 2011-7-1 23:20

看了一晚上 还是没弄得很明白， 看来要理论要恶补一下了 ，好多都忘记了哦！

LION_ARTHAS

inertia relief的问题，建议参考Nastran相关解释，ANSYS在这方面还比较弱。你看到过飞机或者汽车的整车会用ANSYS分析的么？
linuaries 发表于 2011-7-1 23:20

兽兽版主所言极是，有 和 强 是两码事

路隐飘香

纠结呀，差距好大呀，很多都不懂，还得下大功夫啊，强烈支持技术探讨！

gaqrdfwxf

非常感谢，从中受到了很大的启发

甜瓜田瓜

哇塞，这是我看的最长的一个帖子了。。

云少侠

师父太强了……顶起

parrrr

学习了 同时也看到差距，得努力学习了

a1206901110

LION_ARTHAS 发表于 2011-5-18 17:47
副帖：
本副帖主要阐述主贴里出现的问题：莫名的应力集中，不是版本问题，而是网格质量引起的应力奇异。阐 ...

前辈，你好。这个我在哪里有看过。因为是手机不能上图，简要的说下。这根圆轴是轴向受压，要考虑受压后单元的变形。划分时，轴向单元边长要大于径向边长。当轴向受压时，单元变形后会是较规则四边形，避免奇异。单元边长比不是越大越好，有个平衡点，要不断比较不同边长比下，单元的变形情况(若划分时的单元是可以接受的，变形后的单元仍是良好的，这样就比较合理了吧)。我还是个新手，理解有错误的地方请见谅，希望您指正