简单案例分析import引起重力震荡对后续分析的影响程度

wei1012

前一文（http://forum.simwe.com/viewthread.php?tid=980406&extra=page%3D1%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D471）对abaqus软件import from Standard to explicit时发现重力作用出现震荡的问题。这些继续寻找这个问题带来的后续分析影响程度。
（1）隐式计算文件
core.inp--》表示的是隐式下完成施工模拟和地震波分析；
（2）显式计算文件
core_std.inp---》表示的是隐式下完成施工模拟；
core_xpl.inp---》表示的是显式下地震波分析，初始状态来自core_std.inp结果；
core_xpl.inp中为两个step，第一个仅为保持重力作用的step，第二个为时程分析；
这三个文件详见附件；

（3）两个结果对比
对比三种状态，可详见附件中的三个图：
图1左边为隐式状态的重力最终结果，图1右边为显式import过来的最初状态；
对比结果一致

图2左边为隐式下12s的混凝土受压损伤，图2右边为显式下12s的混凝土受压损伤；
对比结果，两者相差很大；




图3左边为隐式下12s的混凝土受拉损伤，图3右边为显式下12s的混凝土受拉损伤；
对比结果，两者相差很大，表现在数值大小和损伤分布上。



结论：根据这个结果，貌似import对后续的时程分析影响很大。
注：用隐式和显式 直接计算重力+动力分析，两者结果差异不大。


是这样吗？？？

wei1012

论坛里面有比较多关于import的好内容，Robert_Su大牛也做了非常全的汇总。汇总内容也提到对于重力作用下隐式输入到显式时仍然需要施加重力作用。但是否会出现震荡，有些案例是没有震荡，有些却有，对震荡是否可以忽略也没有一个说明。

这里用两个案例，一个是实体，另一个是框架-核心筒（梁柱和壳单元），仅受重力荷载作用下，计算结果显示实体单元是没有震荡的；但第二个案例却出现震动，具体该怎么避免，个人还没有什么好的方法。也请大家一起来讨论。
(一) 实体案例
下面两个图片，第一个左边是隐式计算最后结果，右边是显式的最初状态，表示import正确；
第二个左边是是显式的最初状态，右边是显式的1s后结果，结果不一样；


[（2） 一层框架核心筒
附上的图片含义同（1）：


附上这两个案例的inp文件：
（1）实体

（2）一层框架核心筒

最后附上odb：
（1）实体


（2）框架核心筒

Robert_Su

1# wei1012

沒有細看模型,
beam element的section是否太小,
導致explicit code計算時數值擾動對其影響極大?

zsq-w

lz的这个模型的section可能是工程实况。
我通常的做法是：如果施加重力引起的振幅远小于载荷作用所引起的挠度，我就忽略了。

wei1012

我先说说为什么要采用隐式import到显式：
在隐式里面考虑model change的自重作用，在此基础上输入到显式计算地震波。

在重力作用下如果都震荡的话，那么担心显式的结果可能会受影响，毕竟显式的计算过程没有隐式那么严密。

刚才也发现一个问题：我把显式文件中的重力作用删掉，结果和有重力作用下的几乎一样，那是不是说明在显式里面作用的重力是不其作用的？

zsq-w

我把显式文件中的重力作用删掉，结果和有重力作用下的几乎一样，wei1012 发表于 2011-4-26 13:46

你是说显式分析的时候不施加任何载荷，结果也存在震荡？

wei1012

补充说明一下：
以noestory为例，计算noestory.inp，然后计算abaqus j=onestory_xpl，oldjob=noestory  int
注意：onestory_xpl文件中将*dload关于grav定义删掉，计算。

zsq-w

6# wei1012

                                      


我刚刚做了一个一榀框架的小模型，做2种工况的显示运算（时间1s，见上图高20m，宽10m，截面为圆截面，半径为0.05）：
1）只是施加重力
2）施加横向力+重力

后处理：
模型1）输出A点竖向位移，见图Y(gravity-only).JPG ，位移的确一直在震荡，振幅0.08
模型1）输出A点竖向+横向位移，见图X-vs-Y-displacement.JPG
横向位移随着载荷增加，竖向位移和模型1）一样的震荡，竖向位移远小于横向位移，即使应力不太高的时候（应力水平：最大Mises应力78MPa）

附上2个inp

Robert_Su

1# wei1012

可否上傳CAE文件?

wei1012

2# Robert_Su

附件为隐式计算的CAE模型。计算直接采用inp，没有用CAE建模计算。所以从隐式到显式的CAE模型，就需要手工操作下。

Robert_Su

工作實在太忙了,抽不出時間來探討這個問題,提供幾點供樓主試試:
1.explicit是否使用雙精度計算.
2.explicit比較適用於1s以內的計算,樓主的explicit計算時間為12 s,能否比較0.1s, 0.5s,1s時的計算結果.

wei1012

11# Robert_Su

1. explicit采用单精度；
之前工作比较过（a）在隐式下计算重力和地震波 和（b）在显式下计算重力和地震波，显式仍然采用单精度；这两种情况不考虑model change下对比，结果相近。
在此基础上对比隐式计算重力，import到显式计算地震波，发现结果存在较大差异。
2. 比较过结果，重力一直在震荡。

sword1963

震荡的原因应该主要是弹性材料到弹塑性材料的变化。