ABAQUS若干技巧

后处理中create X-Ydata

1. 选择ODB field output，对话框里选择加载点的RF，建立reference force和加载步之间的关系，保存；

, J3 K: z, o0 {( p' O1 `, R. V2. 在ODB field output里，对话框里选择需要位移的节点和U，建立位移U和加载步间的坐标，保存。

3. 选择operate on XY data，combined上两步里的结果。

Initial Conditions和Excel的使用

: l( N  K, P) H6 R- o

By 小梦

关键字格式：

$ n# A- p+ K7 E+ I5 _! q

*initial conditions, type=stress,input=bb.dat

上面的关键字插于*STEP语句之前，两语句之间不能有空格。

施加预应力场只是initial conditions关键字的一个应用，详见abaqus6.8帮助文档，《ABAQUS Analysis User’s Manual》的第28.2节“initial conditions”。

实例：平衡初始地应力

" l  {% x! z# \\  |# [平衡条件：由应力场形成的等效节点载荷要和外载荷相平衡，如果平衡条件得不到满足，将不能得到一个位移为0的初始状态，此时所对应的应力场也不再是所施加的初始应力场。

解决方法：首先将重力载荷施加于土体上，施加符合工程实际条件的边界条件，计算得到在重力载荷下的应力场，再将得到的应力场定义为初始应力场，和重力载荷一起作用于原始的有限元模型，就可以得到既满足平衡条件又不违背屈服准则的初始应力场，可以保证各节点的初始位移近似为0。

步骤：

1 P$ R5 f4 ?, d" e7 O) O2 E) z5 u1、  1、建立有限元模型，部件类型为轴对称模型，网格类型为CAX4R。

下载 (12.33 KB)

2009-7-12 14:48

) c. C  A! m& ^2、建立分析步：Geostatic

下载 (57.63 KB)

2009-7-12 14:48

; L) A7 V# p( Y# r5 M# f1、  3、建立载荷，在Geostatic分析步中，只需要施加重力载荷。

  E6 d! v- o0 g 

下载 (8.12 KB)

2009-7-12 14:48

4

4、创建工作，进行分析

! p1 Y! Y2 V% R' ?- |, q" u5 b5、将分析得到的应力场保存为一个文本文件，inp格式。

: ~, n& d8 b$ M6 V' _  L 

下载 (22.33 KB)

2009-7-12 14:56

下载 (15.91 KB)

2009-7-12 14:56

在工作目录下会生成一个bb.inp文件

1、  6、用excel打开bb.inp文件。注意：先打开一个空白的excel文件，再导入bb.inp文件。

删去多余的行号和列号，只各个应力分量列。

3 ]- l3 T4 t3 f! [* |7 G将各个应力分量的科学计数法格式改为显示小数点后5位。如下图所示，将第一列添加：Soil-1.1、Soil-1.2、…，先在第一列第一行填入Soil-1.1，然后使用excel的拖拽功能。

下载 (5.99 KB)

2009-7-12 14:58

下载 (15.37 KB)

2009-7-12 14:58

" d8 }/ _$ d8 w7 F* \ 

这是应力场

下载 (10.37 KB)

2009-7-12 15:01

* v$ i; T/ K& i) T! e: h 

1、  7、将文件另存为“CSV(逗号分隔)”，出现的提示信息，单击“是”即可。

2、  8、将CSV文件重命名为dat文件。

9、返回原来的Abaqus模型，单击edit keywords。

下载 (10.35 KB)

2009-7-12 15:05

$ L4 q. W0 Q& t3 q 

下载 (18.03 KB)

2009-7-12 15:05

10、重命名工作job，进行分析。可以发现，应力结果和原先一样，但是位移结果接近于零了。

7 ^: i8 J& f) B# d3 S4 K% n& o 

我也没做过土木工程案例，initial conditions, type=stress主要是用来施加初始应力场，平衡地应力只是它应用的一个方面，不知道下面的说明能不能解您的疑惑~侧向压力可以通过在下面的关键字中添加系数来定义。

定义海拔相关的初始应力场：*INITIAL CONDITIONS,TYPE=STRESS,GEOSTATIC

您可以定义海拔相关的初始应力场。当一个地应力初始条件指定于一个单元集，垂直方向的应力在垂向坐标上是分段线性变化的（三维模型中垂直方向是Z方向，轴对称模型和二维模型中垂直方向是Y方向）

  ]5 n- |3 R, j; Q对于垂向应力分量，你应该给定两对应力-海拔（高度）关系来定义单元集中的应力分布。对于这两个给定海拔之间的材料点，Abaqus将使用线性插值方法来确定其初始应力场；对于位于这两个给定海拔之外的点，Abaqus将使用线性外插外推法来确定。另外，水平线方向（侧压力）分量由两个“侧压力系数”来给定，垂向压力乘以给定的系数。在轴对称模型中，只需要给定一个侧压力系数。

单元：

耐地压初始应力只能应用于实体单元（包括三维实体单元和二维实体单元）。在Abaqus/Standard中随海拔变化的的初始应力场应该使用用户子程序SIGINI来为梁单元和壳单元定义。在Abaqus/Explicit中不支持对于梁单元和壳单元的定义。

& p& P. h: S1 o4 x$ ^耐地压应力提案件应该和给定的载荷（如重力）和边界条件保持平衡。初始分析步（initial）将会检查这种平衡。

; E% R( |6 W# ^, h! a2 n& K; i+ }文件输入格式：*INITIAL CONDITIONS,TYPE=STRESS,GEOSTATIC

9 p' T- o" [" b' ^* k% A不支持CAE中定义。

关键字的数据行：

3 P( W0 T7 y! b2 A/ c. t1、单元号或者单元几何标签

, P4 H% \\, I& l7 H2、垂向应力分量的第一个值

3、对应于上面值的点的垂向坐标

3 W) M4 n! a7 H, K5 i4 Y1 G. Q4、垂向应力分量的第二个值

5、对应于上面值的点的垂向坐标

6、侧压力的第一个系数。这个系数定义了X方向的应力分量。

3 V5 U9 E0 `3 p( O) N, T0 j& A9 T7、侧压力的第二个系数。这个系数在三维模型中定义了Y方向的应力分量，在平面（平面应力和平面应变）和轴对称模型中定义了厚度方向的应力分量。如果这个值没定义，假定和第一个给定的系数一样。

如果有必要的话，可以继续定义新的数据行来定义不同的单元或者单元几何的初始耐地压应力。& l4 b# j. d" q1 c" ~6 _6 d  N* |/ s  m: {+ W