abaqus_后处理

abaqus_后处理
abaqus_后处理

22 后处理

22.1 显示局部坐标系上的结果

问:我前处理用的是直角坐标系,但是我想在后处理中输出关于柱坐标的位移分量

是不是要设计局部的坐标系 ?怎样设计?

答:后处理时点菜单tools / coordinates system / create, 创建柱坐标系(例如使用默认的名称csys-1). 菜单result / options, 点tranformation, 点user-specified, 选中csys-1, 点OK. 窗口左上角显示的变量如果原来是U, U1,现在就变为U,U1(CSYS-1).

22.1 显示局部坐标系上的结果

问:我前处理用的是直角坐标系,但是我想在后处理中输出关于柱坐标的位移分量

是不是要设计局部的坐标系 ?怎样设计?

答:后处理时点菜单tools / coordinates system / create, 创建柱坐标系(例如使用默认的名称csys-1). 菜单result / options, 点tranformation, 点user-specified, 选中csys-1, 点OK. 窗口左上角显示的变量如果原来是U, U1,现在就变为U,U1(CSYS-1).

22.2 绘制曲线(X–Y data)

问:例如我想用odb文件建立这样一个曲线:

x y(自行指定)

currentmax("my-xy01") 1.0

currentmax("my-xy02") 3.3

………………

搜索了半天也找不到,在此向用过的前辈请教,或者有第三方软件也请指点。(虽然我手工

excel也能做,但是那个时间就……-_-b)

答:

后处理菜单 tools / XY Data / Create, 先创建你的"my-xy01",再选择Operate on XY Data,

来建立公式。

详见 Getting Started with ABAQUS 附录D.11 Operating on X–Y data

问:我在学习forming a channel例子的时候,发现帮助是将冲头力和位移画在一张图里了(图13-17,见getting started with abaqus),并且将explicit 和standard的结果也都画在一张图中了,可是我只能将field output的结果,也就是随时间变化的结果画出来,

怎么修改x轴的输出量呢?

答:后处理菜单tools / xy-data / manager,点create, 选odb field output,为每条曲线定义这样一个xy-data (两个不同的odb文件也没关系),选中多个xy-data, 点plot,

就可以把它们画在一起。

问:非线性分析,包括n-r,rik法,各增量步的载荷显示。我想画载荷位移曲线

还有,时间步长,弧长和载荷的关系是什么呢。

答:如果你没有为载荷设专门的Amplitude, 载荷大小就是你说得那样。

Field output中的CF就是集中载荷。如果你的载荷是pressure,恐怕不能直接输出,可以用distributing coupling 把受力面和一个参考点耦合起来,然后在参考点上施加等效的

集中载荷。

载荷位移曲线: Getting Started with ABAQUS 附录D.11 Operating on X–Y data 有一

个应力应变曲线的例子。

22.3 field output / history output

问:field variable 是什么概念,什么时候应用?请赐教

答:field output 的输出结果来自于整个模型或模型的大部分区域,被写入ODB文件的频率相对较低,用来在Visualization功能模块中生成云纹图、变形位移图、矢量图和XY图。

例如,使用field output 来在一个分析步结束时输出整个模型的位移场。

history output 的输出结果来自于模型的一小部分区域,被写入ODB文件的频率相对较高,用来在Visualization功能模块中生成XY图. 例如,使用history output 来输出某个节

点在所有时间增量步上的位移。

简言之,field output是某个量随空间位置的变化,history output 是某个量随时间的变

化。

问:桩土1/4模型,施加的symm边界约束,桩顶表面施加pressure,桩体内 embed 了 rebar,结果运行之后桩顶的四个顶点被疯狂地往下压,其他地方都没反应!

这是什么错误啊!

答:是显示的问题。你的field output只输出桩顶部的位移,不输出其他节点的位移,造成后处理显示异常。改为默认的field output就正常了。

如果桩底面和侧面都用接触,等于在底面的圆弧上定义了两次接触,同一个节点有两个法线方向,会出现收敛问题。应该在此处倒一个小的圆角(土也倒圆角),圆角处不定义接触,

就没有这种问题。

不倒角的话,土在拐角处的应力在理论上是无穷大的,可能也会造成收敛问题。

22.4 显示结果

问:不同的set能否用不同的颜色表示出来?同样被赋予不同材料属性的单元能否也用不同

的颜色或其他方式相互区别开来?

答:后处理: 菜单tools / color code

ben_shi2答:云图叠加详细步骤

云图叠加详细步骤:(分别用于应力、变形位移、等)

前提:同一几何模型建在同一cae文件, 在step1, 定义load1; step2,定义load2 (同时inactive load1), run 得到同一模型在两重载荷状态下的云图(分别存在step1和step2).

仿真分析,有限元,模拟,计算,力学,航空,航

天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent;m j$v i G+e

叠加步骤:(应力叠加)

“visualization”à “tools”à “create field outputs”à “from field” à in the opened window: give name, for example “AA” àchoose “step1” and last frameà click “s1f0-s”àclick “+”àchoose “step1” and last frameà click “s1f0-s”àclick

“apply”https://www.360docs.net/doc/b115623834.html, L9D Z1K q:w

To see the combined results, move on to “Result” à”field output”à “AA”|https://www.360docs.net/doc/b115623834.html,|有限元|模拟|计算|力学|航空|航天

|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluent j [^9f;|#W

j n E$N

应力叠加不叠加变形,若看变形叠加结果,重复以上叠加步骤,选择位移即可。

欢迎批评、补充。

问:abaqus得到的是积分点处的应力,但是由于不知道从属于那个面,这时候的应力方向

如何确定呢?现在我非常着急

答:菜单plot / symbols, 可以画矢量图。

问:谢谢各位。现在明白了几点:https://www.360docs.net/doc/b115623834.html, R;i"[#r O2\

I

仿真分析,有限元,模拟,计算,力学,航空,航

天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM c8g

Z

a0XX,l

1. 比较rpt文件输出的Element Nodal 和 Unique Nodal结果,Unique Nodal的结果是同一节点相邻单元在此节点处应力结果的平均值,即后者通过前者在同一节点处平均所得,包括mises,s11,s22,s33,t12,t13,t23都是平均所得的。是否可以说Element Nodal 的结果就是通过高斯点外推出的单元节点处的应力,Unique Nodal就是那个加权平均值?

https://www.360docs.net/doc/b115623834.html,8Z

X){J:KY:M&[ L

r!q V4Z x A

i J s [

2. 比较rpt文件输出的Unique Nodal结果和plot->contours view出的mises应力结果,在我的模型中,两者最大mises应力值差别很大(20%),最大mises应力值出现的位置(节

点号)也不一样。前者是加权平均的话,后者是如何得出的呢?

仿真分析,有限元,模拟,计算,力学,航空,航

天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent:Q'V0J#v0t2y7O

\ B

3. 我猜想plot->contours view出的mises应力结果好像与积分点的应力值关系很大,但不知道如何把显示的mises应力结果输出到结果文件中保存,在rpt文件输出选项中没有这一项,rpt输出的只是积分点上应力。由于在强度校核时,经常要关心某些区域甚至某些点处的最大应力值,所以想把contours view出的mises应力结果输出到文件中保存下来,不知能否办到?我觉得这个应力值和rpt输出的Unique Nodal应力值有差别,但不知道差别

在哪里?

答:看帮助吧。ABAQUS/CAE User's Manual

24.4.1 Understanding how results are computed

24.4.2 Understanding result value averaging

Element-based field output variables are written to the output database at the integration points, the element centroid, or the element nodes, depending on the variable. For the display of element-based field output variables, ABAQUS/CAE reads values from the output database for all elements connected to all nodes included in the plot. Computations are then applied to these values to produce contours, nodal probe results, display groups or color coding based on results, or X–Y data along

a path.

问:要在后处理里面查看两个点的相对位移,应该如何操作?谢谢!

答:

在后处理中可以直接看,菜单Tools / Query , 选 Distance。 其含义见:

问:question1:

用abaqus/cae, visualization, Query看形变后的长度,scaled 和 unscaled有什么不同

答:

实际形变后的长度是Deformed distance (unscaled): 2.02557e+002.

变形前的长度 Base distance: 200

图形显示时的缩放系数 Scale: 8.8011

变形后的长度(不考虑缩放系数)Deformed distance (unscaled): 202.557

(所以变形量为 202.557 - 200 = 2.557 )

变形后的长度 * 缩放系数 Deformed distance (scaled): 222.5

(即 200 + 2.557 * 8.8011 = 222.5)

注意:3个方向都有变形时,变形量是Deformed distance (unscaled)的magnitude 减去 Base distance 的magnitude ,而不是Relative displacement (unscaled)的magnitude 。

问:在选项deformation scale factor里有三个选项

auto-compute

uniform

nonuniform

他们的用法是是什么

是不是这只是一个显示的选项 我选不同的选项 变形的结果不一样

是不是uniform=1才是真正的变形情况

答:auto-compute: ABAQUS自动确定变形的缩放系数。

uniform: 用户自己输入变形的缩放系数,x,y,z方向上的缩放系数相同。

nonuniform: 用户自己输入变形的缩放系数,x,y,z方向上的缩放系数可以不相同。

问:我把一個複合材料薄版

用拉的

但是結果沒出現應力圖示

很奇怪

可否幫我看看哪裡出了錯

謝謝

答:这种模型好像是不能在后处理中显示应力,例题也是这样,见ABAQUS Benchmarks Manual 1.1.3 Composite shells in cylindrical bending中的

compositeshells_c3d20r.inp。

例题是把应力输出到dat文件:

*EL PRINT

1,2,3,4,5,6,7,8,9,10,11,12,13,14,15

S,MISES,PRESS

C3D8单元精度很差,例题中都是用C3D20R或C3D20, 有接触或塑性的话就要用C3D8I,

你的单元太薄,计算结果可能会有问题。

问:变形体用云图,而没有结果的刚体用透明,如何设置呢?

答:菜单options / contour, 把visible edges设为feature edges, 点other标签页,

选translucency.

22.5 输出shell积分点上的应力和应变

问:我用EXPLICIT的SHELL作模拟

怎样输出每个积分点上的应力和应变呀?

學生答:*node print

1,2,3,4,5,6

s

1~~6表示積分點由最下面一層算起

一層有三個(看你用什麼積分法就有幾個積分點我選用三個的)

所以在層與層的交界處會有兩點在同一平面重複

像是我舉兩層例子

1,2,3 4,5,6

3,4要分開但他們在同一平面上

henryzhang答:在Field output 中可以选择所要输出的积分点

(完整)总结Abaqus操作技巧总结(个人),推荐文档

Abaqus操作技巧总结 打开abaqus,然后点击file——set work directory,然后选择指定文件夹,开始建模,建模完成后及时保存,在进行运算以前对已经完成的工作保存,然后点击job,修改inp文件的名称进行运算。切记切 记!!!!!! 1、如何显示梁截面(如何显示三维梁模型) 显示梁截面:view->assembly display option->render beam profiles,自己调节系数。 2、建立几何模型草绘sketch的时候,发现画布尺寸太小了 1)这个在create part的时候就有approximate size,你可以定义合适的(比你的定性尺寸大一倍); 2)如果你已经在sketch了,可以在edit菜单--sketch option ——general--grid更改 3、如何更改草图精度 可以在edit菜单--sketch option ——dimensions--display——decimal更改 如果想调整草图网格的疏密,可以在edit菜单--sketch option ——general——grid spacing中可以修改。 4、想输出几何模型 part步,file,outport--part 5、想导入几何模型? part步,file,import--part 6、如何定义局部坐标系 Tool-Create Datum-CSYS--建立坐标系方式--选择直角坐标系or柱坐标系or球坐标 7、如何在局部坐标系定义载荷

laod--Edit load--CSYS-Edit(在BC中同理)选用你定义的局部坐标系 8、怎么知道模型单元数目(一共有多少个单元) 在mesh步,mesh verify可以查到单元类型,数目以及单元质量一目了然,可以在下面的命令行中查看单元数。 Query---element 也可以查询的。 9、想隐藏一些part以便更清楚的看见其他part,edge等 view-Assembly Display Options——instance,打勾 10、想打印或者保存图片 File——print——file——TIFF——OK 11、如何更改CAE界面默认颜色 view->Grahphic options->viewport Background->Solid->choose the wite colour! 然后在file->save options. 12、如何施加静水压力hydrostatic load --> Pressure, 把默认的uniform 改为hydrostatic。这个仅用于standard,显式分析不支持。 13、如何检查壳单元法向 Property module/Assign/normal 14、如何输出单元体积 set步---whole model ----volume/Tickness/Corrdinate-----EVOL 15、如何显示最大、最小应力 在Visualization>Options>contour >Limits中选中Min/Max:Show Location,同样的方法可以知道具体指定值的位置。 16、如何在Visualization中显示边界条件 View——ODB display option——entity display——show boundary conditions 17、后处理有些字符(图例啊,版本号啊,坐标系啊)不想显示, viewport-viewport annotation option ,选择打勾。同样可以修改这些字体大小、位置等等。

abaqus 后处理

22 后处理 22.1 显示局部坐标系上的结果 问:我前处理用的是直角坐标系,但是我想在后处理中输出关于柱坐标的位移分量 是不是要设计局部的坐标系?怎样设计? 答:后处理时点菜单tools / coordinates system / create, 创建柱坐标系(例如使用默认的名称csys-1). 菜单result / options, 点tranformation, 点user-specified, 选中csys-1, 点OK. 窗口左上角显示的变量如果原来是U, U1,现在就变为U,U1(CSYS-1). 22.1 显示局部坐标系上的结果 问:我前处理用的是直角坐标系,但是我想在后处理中输出关于柱坐标的位移分量 是不是要设计局部的坐标系?怎样设计? 答:后处理时点菜单tools / coordinates system / create, 创建柱坐标系(例如使用默认的名称csys-1). 菜单result / options, 点tranformation, 点user-specified, 选中csys-1, 点OK. 窗口左上角显示的变量如果原来是U, U1,现在就变为U,U1(CSYS-1). 22.2 绘制曲线(X–Y data) 问:例如我想用odb文件建立这样一个曲线: x y(自行指定) currentmax("my-xy01") 1.0 currentmax("my-xy02") 3.3 ……………… 搜索了半天也找不到,在此向用过的前辈请教,或者有第三方软件也请指点。(虽然我手工 excel也能做,但是那个时间就……-_-b) 答: 后处理菜单 tools / XY Data / Create, 先创建你的"my-xy01",再选择Operate on XY Data,

ABAQUS后处理中各应力解释个人收集

ABAQUS中的壳单元S33代表的是壳单元法线方向应力,S11 S22 代表壳单元面内的应力。因为壳单元的使用范围是“沿厚度方向应力为0”,也即沿着法相方向应力为0,且满足几何条件才能使用壳单元,所以所有壳单元的仿真结果应力查看到的S33应力均为0。 S11 S22 S33 实体单元是代表X Y Z三个方向应力,但壳单元不是,另外壳单元只有S12,没有S13,S23。 LE----真应变(或对数应变) LEij---真应变 ... 应变分量; PE---塑性应变分量; PEEQ---等效塑性应变 ABAQUS Field Output Stresses S stress components and invariants 应力分量和变量 SVAVG volume-averaged stress components and invariants (Eulerian only) MISESMAX 最大 Mises 应力 TSHR transverse shear stress(for thick shells)横向剪切应力 CTSHR transverse shear stress in stacked continuum shells 连续堆垛壳横向剪切应力 TRIAX stress triaxiality 应力三轴度 VS stress in the elastic-viscous network 弹粘性网格应力PS stress in the plastic-viscous network 塑粘性网格应力 SFABRIC stress components in fabric constitutive measure 纤维本构应力分量 SSAVG average shell section stress 平均壳节点应力 Strain E total strain components 总应变分量 VE viscous strain in the elastic-viscous network 黏性应变弹黏性网格PE plastic strain components 塑性应变分量 PEVAVG volume-averaged plastic strain components (Eulerian only) VEEQ equivalent viscous

abaqus后处理中各应力解释个人收集修订版

a b a q u s后处理中各应力 解释个人收集修订版 IBMT standardization office【IBMT5AB-IBMT08-IBMT2C-ZZT18】

ABAQUS中的壳单元S33代表的是壳单元法线方向应力,S11 S22 代表壳单元面内的应力。因为壳单元的使用范围是“沿厚度方向应力为0”,也即沿着法相方向应力为0,且满足几何条件才能使用壳单元,所以所有壳单元的仿真结果应力查看到的S33应力均为0。 S11 S22 S33 实体单元是代表X Y Z三个方向应力,但壳单元不是,另外壳单元只有S12,没有S13,S23。 LE----真应变(或对数应变) LEij---真应变 ... 应变分量; PE---塑性应变分量; PEEQ---等效塑性应变 ABAQUS Field Output Stresses S stress components and invariants 应力分量和变量 SVAVG volume-averaged stress components and invariants (Eulerian only) MISESMAX 最大 Mises 应力 TSHR transverse shear stress(for thick shells)横向剪切应力 CTSHR transverse shear stress in stacked continuum shells 连续堆垛壳横向剪切应力 TRIAX stress triaxiality 应力三轴度 VS stress in the elastic-viscous network 弹粘性网格应力 PS stress in the plastic-viscous

总结Abaqus操作技巧总结

打开abaqus,然后点击file——set work directory,然后选择指定文件夹,开始建模,建模完成后及时保存,在进行运算以前对已经完成的工作保存,然后点击job,修改inp文件的名称进行运算。切记切 记!!!!!! 1、如何显示梁截面(如何显示三维梁模型) 显示梁截面:view->assembly display option->render beam profiles,自己调节系数。 2、建立几何模型草绘sketch的时候,发现画布尺寸太小了 1)这个在create part的时候就有approximate size,你可以定义合适的(比你的定性尺寸大一倍); 2)如果你已经在sketch了,可以在edit菜单--sketch option ——general--grid更改 3、如何更改草图精度 可以在edit菜单--sketch option ——dimensions--display——decimal更改 如果想调整草图网格的疏密,可以在edit菜单--sketch option ——general——grid spacing中可以修改。 4、想输出几何模型 part步,file,outport--part 5、想导入几何模型? part步,file,import--part 6、如何定义局部坐标系 Tool-Create Datum-CSYS--建立坐标系方式--选择直角坐标系or柱坐标系or球坐标 7、如何在局部坐标系定义载荷

laod--Edit load--CSYS-Edit(在BC中同理)选用你定义的局部坐标系 8、怎么知道模型单元数目(一共有多少个单元) 在mesh步,mesh verify可以查到单元类型,数目以及单元质量一目了然,可以在下面的命令行中查看单元数。 Query---element 也可以查询的。 9、想隐藏一些part以便更清楚的看见其他part,edge等 view-Assembly Display Options——instance,打勾 10、想打印或者保存图片 File——print——file——TIFF——OK 11、如何更改CAE界面默认颜色 view->Grahphic options->viewport Background->Solid->choose the wite colour! 然后在file->save options. 12、如何施加静水压力hydrostatic load --> Pressure, 把默认的uniform 改为hydrostatic。这个仅用于standard,显式分析不支持。 13、如何检查壳单元法向 Property module/Assign/normal 14、如何输出单元体积 set步---whole model ----volume/Tickness/Corrdinate-----EVOL 15、如何显示最大、最小应力 在Visualization>Options>contour >Limits中选中Min/Max:Show Location,同样的方法可以知道具体指定值的位置。 16、如何在Visualization中显示边界条件 View——ODB display option——entity display——show boundary conditions

利用Python对Abaqus进行后处理结果输出

利用Python对Abaqus进行后处理结果输出 -----中大_戚超_2016.10.3 1 概述 在Abaqus的二次开发过程中,通常需要采用Python脚本语言将Abaqus的计算结果进行输出,然后再进行处理。Python使Abaqus的内核语言,使用较为方便,Abaqus运行Python语言的方式有多种,可以直接命令窗口,也可以读入脚本,还可以采用类似批处理的方式。 本次以一个例子细说Python语言在Abaqus后处理中的应用,模型的计算结果云图如图1所示。 图1 计算结果 2 输出所有积分点上的Mises应力 直接上Python代码: import os from odbAccess import* from textRepr import* myodb=openOdb(path='Job-1.odb') cpFile=open('artlcF1.txt','w') RF=myodb.steps['Step-1'].frames[1].fieldOutputs['S'].values for i in range(len(RF)) :

cpFile.write('%.3F\n' %(RF[i].mises)) else: cpFile.close() #引入模块,因为需要打开结果文件 #打开结果文件,并复制给变量myodb #打开一个txt文件 #将输出场赋值给RF #循环语句,向txt文件逐行写入mises应力 Abaqus的结构层次分的很细,比如结果文件下分如下: 使用过Abaqus的都知道step表示载荷步,frame表示载荷子步,因而在读取Mises应力时需要详细地指定输出哪一步的应力,而应力结果是输出场数据(fieldOutput)的中一种,需要指定是何种应力,程序才知道怎么读取并写入。 由于Abaqus里面涉及的变量特别多,通常很难记清楚那一项下面都有哪些量可以调用,此时比较好的方式是采用print 函数查看,例如查看myodb.steps['Step-1'].frames[1].fieldOutputs 下面有哪些变量可以调用,在窗口输入: print myodb.steps['Step-1'].frames[1].fieldOutputs 显示:

【2018年整理】abaqus的后处理

1 如何显示最大、最小应力 在Visualization>Options>contour >Limits中选中Min/Max:Show Location 2 后处理有些字符(图例啊,版本号啊,坐标系啊)不想显示 view-view annotation option ,选择打勾 3 后处理中显示边界条件 Viewport--ODB Display Options 边界条件处打勾 4 在模型上只显示云图,不显示网格 option菜单-common-visible edges--free edges 顺便窜一下,在前处理mesn之后不想显示网格,只显示模型的话,更简单:工具栏有显示mesh、显示线框、显示实体连续的三个按钮。 5 你想调大变形放大系数(Deformation Scal Factor)让变形显示更明显一点? 注意:非线性问题,这个默认为是1(也就是不放大),只有线性的才能改。 option菜单-common-visible edges--Deformation Scale factor 6 如何在后处理中移动图例 选Gneral--legend 可以隐藏 选Legent--Upper left corner 可以移动 7 对ODB结果处理以后,如何保存结果 從visualization模塊開啟odb檔(不要直接從job manager開啟),並取消勾 選 read only,對已有的XY-data作copy to odb。 8 job步提交运算后警告信息出现setwarning 或者nodewarning,(类似这样 的语句: ....have been identified in element set ErrElemVolSmallNegZero .... identified in element set ErrElemDistorted WarnNodeUnconstrained WarnNodeMissMasterIntersect

ABAQUS前_后处理模块二次开发的应用_朱兆华

文章编号:1001-2265(2009)01-0030-04 收稿日期:2008-08-11  *基金项目:江西省科技厅科技支撑项目[2007];江西省教育厅2009年度一般科技项目(G J J 09025);江西省教育厅科技研究项目(G J J 08429) 作者简介:朱兆华(1979—)男,江苏涟水人,南昌大学机电学院硕士研究生,助理工程师,研究方向为材料加工新技术,(E-m a i l )b y j 245@163. c o m 。 A B A Q U S 前、后处理模块二次开发的应用 * 朱兆华1 ,黄菊花1 ,张庭芳1 ,谢世坤2 ,白引娟 3 (1.南昌大学机电工程学院,南昌 330031;2.井冈山大学工学院,江西吉安 343009;3.江西协中汽车 内饰有限公司,南昌 330032) 摘要:文章以实例说明了P y t h o n 脚本语言和A B A Q U S G U I T o o l k i t 在A B A Q U S 的前、后处理模块二次开发中的应用,并阐述了不同模块之间的调用流程。通过二次开发程序控制A B A Q U S 的建模和装配的过程,有效地解决了模型装配时的繁琐、易错等问题,提高了前处理的效率;因A B A Q U S 对板料拉深进行数值 模拟的后处理功能不够全面,为扩展后处理的功能,更好地查看和分析模拟的结果,文章对A B A Q U S 后处理进行二次开发来达到这一目的。 关键词:A B A Q U S ;P y t h o n ;A B A Q U S G U I T o o l k i t ;二次开发中图分类号:T H 16;T G 65 文献标识码:A A p p l i c a t i o no f S e c o n d -d e v e l o p e d o n A B A Q U S P r e -p r e c e s s a n d P o s t -p r o c e s s Z H UZ h a o -h u a 1,H U A N GJ u -h u a 1,Z H A N GT i n g -f a n g 1,X I ES h i -k u n 2,B A I Y i n -j u a n 3 (1.S c h o o l o f M e c h a n i c a l a n dE l e c t r i c a l E n g i n e e r i n g ,N a n c h a n g U n i v e r s i t y ,N a n c h a n g 330031;2.C o l l e g e o f E n g i n e e r i n g ,J i n g g a n g s h a n U n i v e r s i t y ,J i a n g x i J i 'a n 343009,C h i n a ) A b s t r a c t :T h i s a r t i c l e u s i n g a n e x a m p l e t o e x p l a i n P y t h o n a n d A B A Q U S G U I T o o l k i t p r o g r a m m e f o r A B A Q U S p r e -p r o c e s s a n d p o s t -p r o c e s s s e c o n d -d e v e l o p m e n t a n de x p o u n dt h e p r o c e s s a m o n g d i f f e r e n t m o d u l e s 。T h o u g h s e c o n d -d e v e l o p m e n t t h e p r o c e s s o f m o d e l a n da s s e m b l e c a n b e c o n t r o l l e d ,t h e t e d i o u s 、e r r o r -p r o n e a n do t h e r p r o b l e m s c a n b e e f f e c t i v e l y s o l v e d ,s o t h ew o r k i n g e f f i c i e n c y c a nb e h i g h l y i m p r o v e di nt h e p r e -p r o c e s s ;F o r d r a w i n g s h e e t ,p o s t -p r o c e s s f u n c t i o n s p r o v i d e db y A B A Q U S a r ei n s u f f i c i e n t .I no r d e r t o e x t e n dp o s t -p r o c e s s f u n c t i o n s a n d e x a m i n e s i m u l a t i o n r e s u l t e f f e c t i v e l y s e c o n d d e v e l o p o n A B A Q U S p o s t -p r o c e s s i s u s e d i n t h i s p a -p e r . K e y w o r d s :A B A Q U S ;P y t h o n ;A B A Q U S G U I T o o l k i t ;s e c o n d -d e v e l o p m e n t 0 引言 A B A Q U S 是国际上最先进的大型通用有限元计算分析软件之一,可以模拟绝大部分工程材料的线性和非线性行为。A B A Q U S 自带的C A E 是进行有限元分析的前后处理模块,也是建模、分析和后处理的人机交互平台,它具有良好的人机对话界面,因此A B A Q U S 软件在工程中得到了广泛的应用。 P y t h o n 是一种面向对象的脚本语言,它功能强大,既可以独立运行,也可以用作脚本语言。特别适用于 快速的应用程序开发。 1 A B A Q U S /C A E 处理过程和二次开发接口 介绍 A B A Q U S /C A E 处理有两个程序:①内核程序;②G U I 程序。内核程序实际上就是它的脚本语言,它采用的是P y t h o n 语言,同时扩展了P y t h o n 语言,额外提供了大约500个对象模型,对象模型之间的关系复杂,它们间部分关系如图1所示。 图1中,C o n t a i n e r 表示容器,里面包括其他对象, · 30·

abaqus--后处理

1 如何显示最大、最小应力 在V isualization>Options>contour >Limits中选中Min/Max:Show Location 2 后处理有些字符(图例啊,版本号啊,坐标系啊)不想显示 view-view annotation option ,选择打勾 3 后处理中显示边界条件 V iewport--ODB Display Options 边界条件处打勾 4 在模型上只显示云图,不显示网格 option菜单-common-visible edges--free edges 顺便窜一下,在前处理mesn之后不想显示网格,只显示模型的话,更简单:工具栏有显示mesh、显示线框、显示实体连续的三个按钮。 5 你想调大变形放大系数(Deformation Scal Factor)让变形显示更明显一点? 注意:非线性问题,这个默认为是1(也就是不放大),只有线性的才能改。 option菜单-common-visible edges--Deformation Scale factor 6 如何在后处理中移动图例 选Gneral--legend 可以隐藏 选Legent--Upper left corner 可以移动 7 对ODB结果处理以后,如何保存结果 從visualization模塊開啟odb檔(不要直接從job manager開啟),並取消勾 選 read only,對已有的XY-data作copy to odb。 8 job步提交运算后警告信息出现setwarning 或者nodewarning,(类似这样 的语句: ....have been identified in element set ErrElemV olSmallNegZero .... identified in element set ErrElemDistorted WarnNodeUnconstrained WarnNodeMissMasterIntersect

完整word版ABAQUS后处理一些知识

提取剪力:即可RF 立SET 提取是否应该在计算之前,把某一层所有柱子底部的单元定义为一个集合,然后在历程输出中,定义这个集合的剪力输出?这样是不是就是某层的剪力时程曲线?关于剪力问题,大家是怎么提取的,是提取RF?Sxx?还是别的。。。 直接把某层的柱子定义成一个集合,然后提取这个集合的水平力。不知行么? 在abaqus里(猜测LZ应该问的是这个软件) 建模型时,采用不同的单元,输出层间剪力的方式是不同的; 对于beam和shell单元,可以直接在History output 或者field output里面输出SF,即可构件截面剪力,同理,SM为截面弯矩。 对于实体solid单元,则就没有那么简单了,需要在inp文件中定义相应的截面cutting surface,从而在后处理dat文件(stangard隐式计算)或整体截面输出(Explicit显式计算)输出SOF,即可直接输出截面剪力值,同理,相应的弯矩值可采用SOM。 补充一下: 采用beam单元模拟柱时,要注意柱的n1方向与之对应的剪力,否则可能会出现混乱, 即,最好是柱的n1和t方向均相同,以避免出现错误, 举的简单的例子, 某两根柱子的x向剪力时程是v1(t)和v2(t), 但如果v1(t)和v2(t)的正负方向定义反了,出来的合力v1(t)+v2(t)肯定不是所需的。 再有,beam单元输出的内力是积分点的值而非端节点的值,应当引起注意, 假如层间柱就划为2个单元,此时输出来的V、M内力和柱端差距是比较大的, 讨论。 abaqus 做高层结构分析后,如何提取某一层的剪力,基底剪力可以通过RF输出然后求和,其它层如何处理?SF(SF1-SF6)对柱单元和墙单元来说分别表示什么含义?对于柱,请参看所用单元的element output,直接将sfi求和即可,help中对sfi的说明是很清楚的; 对于壳(即你所用墙)的单元输出说明没有图,理解起来有点费劲。

基于PYTHON的ABAQUS后处理开发

基于PYTHON的ABAQUS后处理开发 ABAQUS的后处理功能不能完全提供我们在分析过程中所需的数据,为更好的扩展后处理功能,查看和分析结果数据,本文提出了使用Python语言对ABAQUS进行二次开发来达到这一目的的方法。文中讨论了ABAQUS的脚本接口和对象模型在二次开发中的作用和调用流程,以及文件的读写与复制、数据读取与处理、结果输出与查看等关键技术。以共轨管锥面密封性的分析为例,使用Python语言提取了分析结果数据并将结果作为初始条件加载于新的分析中,最终得到所需的分析数据。 标签:ABAQUS Python 后处理有限元分析 0 引言 ABAQUS是目前国际上最为先进的通用非线性有限元分析软件之一,软件包括种类丰富的材料库和单元库,可以模拟绝大部分工程材料的线性和非线性行为,而且材料库和单元库分开,材料和单元之间的组合能力很强,可以胜任复杂结构的静态与动态分析[1]。ABAQUS自带的CAE模块是一个完整的ABAQUS环境,提供—个简单一致的接口,可以用于创建、提交、监视和评价模拟所得到的结果[2]。 Python是一种简单易学、功能强大的编程语言,它有高效率的高级数据结构,可以简单而有效地实现面向对象编程[2]。ABAQUS有限元程序就通过集成脚本语言Python向二次开发者提供了很多库函数,通过Python语言调用这些库函数来增强ABAQUS的交互式操作能力。本文通过Python脚本语言来提取ABAQUS 的后处理结果,并对结果数据进行相应的计算和转换,并将处理结果作为初始条件应用于其后的开发过程中。 1 ABAQUS、PYTHON接口程序 ABAQUS的二次开发可以通过Python脚本语言控制ABAQUS内核实现前处理建模和后处理的计算分析。ABAQUS脚本接口是Python语言的一个扩展,可以使用Python语言编制脚本接口可以执行的程序,从而实现自动化重复性的工作、创建和修改模型数据库、访问数据库的功能。ABAQUS脚本接口中的每个对象都拥有相应的数据成员(data)和函数,对象中的函数专门用来处理对象中的数据成员,被称为相应对象的方法,用来生成对象的方法被称为构造函数(constructor)[3]。在对象创建以后,可以使用该对象提供的方法来处理对象中的数据成员。ABAQUS的对象模型是指ABAQUS脚本接口中的各种对象类型之间的层次关系。ABAQUS脚本接口对Python语言进行了扩展,额外提供了大约500个对象模型,这些对象模型之间的关系比较复杂,大致将这写对象分为3类,如图1所示。其中,session对象用来定义视图、远程队列、用户定义的视图等;mdb对象包括计算模型对象和作业对象,其中odb对象包括计算模型和计算结果数据,每一类对象下面又包括各类子对象,比如mdb对象下面的计算模型models对象又包括很

ABAQUS后处理一些知识

提取剪力: 立SET 提取RF 即可 是否应该在计算之前,把某一层所有柱子底部的单元定义为一个集合,然后在历程输出中,定义这个集合的剪力输出?这样是不是就是某层的剪力时程曲线? 关于剪力问题,大家是怎么提取的,是提取RF?Sxx?还是别的。。。 直接把某层的柱子定义成一个集合,然后提取这个集合的水平力。不知行么? 在abaqus里(猜测LZ应该问的是这个软件) 建模型时,采用不同的单元,输出层间剪力的方式是不同的; 对于beam和shell单元,可以直接在History output 或者field output里面输出SF,即可构件截面剪力,同理,SM为截面弯矩。 对于实体solid单元,则就没有那么简单了,需要在inp文件中定义相应的截面cutting surface,从而在后处理dat文件(stangard隐式计算)或整体截面输出(Explicit显式计算)输出SOF,即可直接输出截面剪力值,同理,相应的弯矩值可采用SOM。 补充一下: 采用beam单元模拟柱时,要注意柱的n1方向与之对应的剪力,否则可能会出现混乱,即,最好是柱的n1和t方向均相同,以避免出现错误, 举的简单的例子, 某两根柱子的x向剪力时程是v1(t)和v2(t), 但如果v1(t)和v2(t)的正负方向定义反了,出来的合力v1(t)+v2(t)肯定不是所需的。 再有,beam单元输出的内力是积分点的值而非端节点的值,应当引起注意, 假如层间柱就划为2个单元,此时输出来的V、M内力和柱端差距是比较大的, 讨论。 abaqus 做高层结构分析后,如何提取某一层的剪力,基底剪力可以通过RF输出然后求和,其它层如何处理?SF(SF1-SF6)对柱单元和墙单元来说分别表示什么含义?对于柱,请参看所用单元的element output,直接将sfi求和即可,help中对sfi的说明是很清楚的; 对于壳(即你所用墙)的单元输出说明没有图,理解起来有点费劲。

总结Abaqus操作技巧总结(个人)

Abaqus操作技巧总结 打开abaqus,然后点击file——setwork directory,然后选择指定文件夹,开始建模,建模完成后及时保存,在进行运算以前对已经完成的工作保存,然后点击job,修改inp文件的名称进行运 算。切记切记!!!!!! 1、如何显示梁截面(如何显示三维梁模型) 显示梁截面:view->assembly display option->render beam prof iles,自己调节系数。 2、建立几何模型草绘sketch的时候,发现画布尺寸太小了 1)这个在create part的时候就有approximatesize,你可以定义合适的(比你的定性尺寸大一倍); 2)如果你已经在sketch了,可以在edit菜单--sketch option ——general--grid更改 3、如何更改草图精度 可以在edit菜单--sketch option——dimensions--display——decimal更改 如果想调整草图网格的疏密,可以在edit菜单--sketch option ——genera l——grid spacing中可以修改。 4、想输出几何模型 part步,file,outport--part 5、想导入几何模型? part步,file,import--part 6、如何定义局部坐标系 Tool-Create Datum-CSYS--建立坐标系方式--选择直角坐标系or柱坐标系or球坐标 7、如何在局部坐标系定义载荷

laod--Edit load--CSYS-Edit(在BC中同理)选用你定义的局部坐标系 8、怎么知道模型单元数目(一共有多少个单元) 在mesh步,meshverify可以查到单元类型,数目以及单元质量一目了然,可以在下面的命令行中查看单元数。 Query---element 也可以查询的。 9、想隐藏一些part以便更清楚的看见其他part,edge等 view-Assembly Display Options——instance,打勾 10、想打印或者保存图片 File——print——file——TIFF——OK 11、如何更改CAE界面默认颜色 view->Grahphicoptions->viewport Background->Solid->choose the wite colour! 然后在file->save options. 12、如何施加静水压力hydrostatic load--> Pressure, 把默认的uniform 改为hydrostatic。这个仅用于standard,显式分析不支持。 13、如何检查壳单元法向 Propertymodule/Assign/normal 14、如何输出单元体积 set步---whole model----volume/Tickness/Corrdinate-----EVOL 15、如何显示最大、最小应力 在Visualization>Options>contour>Limits中选中Min/Max:ShowLocation,同样的方法可以知道具体指定值的位置。 16、如何在Visualization中显示边界条件 View——ODB display option——entitydisplay——show boundary conditions 17、后处理有些字符(图例啊,版本号啊,坐标系啊)不想显示, viewport-viewport annotation option ,选择打勾。同样可以修改这

相关主题
相关文档
最新文档