ABAQUS学习笔记
inp编码介绍
(一).ABAQUS头信息文件段(1-4)
1.*PREPRINT 输出求解过程所要求的信息(在dat文件中)
ie:*PREPRINT, ECHO=YES, HISTORY=YES, MODEL=YES
2.*HEADING 标题输出文件(出现在POST/VIEW窗口中,且出现在结果输出文件中)
ie:*HEADING
STRESS ANALYSIS FOR A PLATE WITH A HOLE
3.*RESTART 要求abaqus/standard输出其POST/view模块所需要的.res文件。其中的FREQ =?控制结果在每次迭代(或载荷步)输出的次数。
ie:*RESTART, WRITE, FREQ=1
4.*FILE FORMAT 要求abaqus/standard输出到.fil中的某些信息。它也用于post。对于在后处理中得到x-y形式的诸如应力-时间、应力-应变图有用!
ie: *FILE FORMAT, ZERO INCREMENT
(二).ABAQUS网格生成段
定义结点、单元,常用的命令有:结点定义(*NODE,*NGEN),单元定义(*ELEMENT,*ELGEN 等)。
1.*NODE 定义结点,其格式为:
*NODE
结点号,x轴坐标,y轴坐标,(z轴坐标)
2.*NGEN 在已有结点的基础上进行多个结点的生成,一般是在两结点间以某种方式(直线、圆)产生一定分布规律的结点。
如:*NGEN, LINE=C, NSET=HOLE,
119, 1919, 100, 101 在两结点(结点号为119,1919)间以圆弧形式生成多个结点,100为任意相邻结点的单元号增量,101为圆弧形成时圆心位置的结点(对于直线形式生成没有此结点)。所有这些生成的结点(包括119,1919)被命名成HOLE的集合(这样做的目的是以后的命令中使用到它,比如说对这些结点施加同等条件的边界条件或载荷等,HOLE就是这些结点的代称)。*NGEN使用的前提就是必须存在已有结点。
*NGEN, NSET=OUTER
131, 1031, 100 以线形式形成结点,结点号增量100,结点集合名为OUTER。
*NGEN, NSET=OUTER
1031, 1931, 100 同上生成结点,可以同上结点集合名,这样OUTER就包括这两次生成的所有结点
3.*NFILL 在如上生成的结点集(实际上,代表两条几何意义上的边界线)之间按一定规律(BIAS=?)填充结点。这样所有生成的结点构成一定形状的实体(面)。
如:*NFILL, NSET=PLATE, BIAS=0.8
HOLE, OUTER, 12, 1 以HOLE为第一条边界,OUTER为第二条边界(终止边),以从疏到密的规律(BIAS小于1)分布,其生成结点数在两内外对应结点间为12,1为每组结点号的增量。所有这些结点被置于PLATE的集合中。
下面以上面生成的结点来生成单元:
4.*ELEMENT
定义单元所使用的类型(TYPE=?),然后另行定义通过联结结点形成单元,其结点数目
依靠单元类型而变。
*ELEMENT, TYPE=CPS4 //采用四单元的平面应力单元
19, 119, 120, 220, 219 //定义顺序:单元号,以逆时针方向形成单元的各结点号
(三)ABAQUS单元
注意:分析前要选择合适的元素,这时要考虑的问题就是:使用什么样类型的单元?有限元的基本思路就是将实际中的连续体离散化,实际结果是将众多离散分析结果的集合,这似乎有点像积分的概念。选择元素种类最重要考虑的是分析必要的现象,满足必要的准确度基础上去掉不必要的细节与准确度。是选择1-D, 2-D or 3-D单元、用于何种分析的单元、是否高阶单元等。
(四)ABAQUS材料
ABAQUS本身提供了丰富的材料库供分析使用,并已能满足常用的分析。但对于新型本构关系的材料abaqus本身是无法体现的,UMAT则为这个问题提供了解决。自己编程将材料的应力应变本构表示出来,ABAQUS调用完成分析。
ABAQUS 的材料行为模式主要分为
弹性材料:
Linear elasticity (线弹性)(
No( compression or tension elasticity (无压缩或位伸弹性材料,即单力性材料)
Plane stress( orthotropic failure (平面应力单元)
Porous elasticity (多孔弹性)(
( Hypoelasticity (亚弹性)
Hyperelasticity (超弹性)(
Foam( elasticity (泡沫单元)
Viscoelasticity (粘弹性)(
非弹性材料
( Classical metal plasticity (塑性)
Metals subjected to cyclic( loading (受周期荷载金属单元)
Rate-dependent yield(率相关屈服单元)(
Creep( and Swelling (蠕变)
Anisotropic yield and creep (各向异性)(
( Porous metal plasticity (多孔塑性)
Deformation plasticity ( (塑变单元)
Granular materials or polymers (粒状材料或复合材料)(
( Clay plasticity (粘土塑性)
Crushable foam plasticity ( (可压泡沫塑性)
Jointed material (?……)(
Concrete ( (混凝土)
(五)ABAQUS求解
对于一个inp文件,不进入CAE时,需要这样做:
1. 检查inp文件的正确性(当然主要是指keyword的使用),自己能做检查最好,否则可以通过:ABAQUS datacheck job=yourjobname
2. 检查确认修正后进行计算:
通过:ABAQUS job=yourjobname
3. 检验分析结果的合理性:不只是会算,更要会对分析结果进行确认。首先要对整个分析及分析的并键之处成竹在心。然后可以通过以下途径作结果确认:
①自已能够得到的解析解
②实验数据
③其它数值解
④别人的求解结果(当然你得信任他)
⑤直觉与经验
4. 如果迭代无法收敛:需要通过.msg,.sta文件查看出错信息并做出判断(在CAE中submit分析时可以通过monitor查看),判断依据为:
①结构约束是否足够或过多
②材料数据是否正确
③单元是否适合此分析
④网格有没有过扭曲、奇异
⑤接触单元是否足够
⑥步长是否过大
来自:likuiming收集的资料
[ 本帖最后由yejet 于2006-9-2 14:46 编辑]
Abaqus 使用日记 Abaqus标准版共有“部件(part)”、“材料特性(propoterty)”、“装配(assemble)”、“计算步骤(step)”、“交互(interaction)”、“加载(load)”、“单元划分(mesh)”、“计算(job)”、“后处理(visualization)”、“草图(sketch)”十大模块组成。 建模方法: 一个模型(model)通常由一个或几个部件(part)组成,“部件”又由一个或几个特征体(feature)组成,每一个部分至少有一个基本特征体(base feature),特征体可以是所创建的实体,如挤压体、切割挤压体、数据点、参考点、数据轴,数据平面,装配体的装配约束、装配体的实例等等。 1.首先建立“部件” (1)根据实际模型的尺寸决定部件的近似尺寸,进入绘图区。绘图区根据所输入的近似尺寸决定网格的间距,间距大小可以在edit菜单sketcher options选项里调整。 (2)在绘图区分别建立部件中的各个特征体,建立特征体的方法主要有挤压、旋转、平扫三种。同一个模型中两个不同的部件可以有同名的特征体组成,也就是说不同部件中可以有同名的特征体,同名特征体可以相同也可以不同。部件的特征体包括用各种方法建立的基本特征体、数据点(datum point)、数据轴(datum axis)、数据平面(datum plane)等等。(3)编辑部件可以用部件管理器进行部件复制,重命名,删除等,部件中的特征体可以是直接建立的特征体,还可以间接手段建立,如首先建立一个数据点特征体,通过数据点建立数据轴特征体,然后建立数据平面特征体,再由此基础上建立某一特征体,最先建立的数据点特征体就是父特征体,依次往下分别为子特征体,删除或隐藏父特征体其下级所有子特征体都将被删除或隐藏。××××特征体被删除后将不能够恢复,一个部件如果只包含一个特征体,删除特征体时部件也同时被删除××××× 2.建立材料特性 (1)输入材料特性参数弹性模量、泊松比等 (2)建立截面(section)特性,如均质的、各项同性、平面应力平面应变等等,截面特性管理器依赖于材料参数管理器 (3)分配截面特性给各特征体,把截面特性分配给部件的某一区域就表示该区域已经和该截面特性相关联 3.建立刚体 (1)部件包括可变形体、不连续介质刚体和分析刚体三种类型,在创建部件时需要指定部件的类型,一旦建立后就不能更改其类型。采用旋转方式建立部件,在绘制轴对称部件的外形轮廓时不能超过其对称轴。 (2)刚体是不能够施加质量、惯性轴等特性的,建立刚体后必须给刚体指定一个参考点(reference point)。在加载模块里对参考点施加约束和定义其运动,对参考点施加的荷载或运动就相当于施加给了整个刚体。 4.模型装配 (1)在装配(assemble)模块里首先建立部件实例(part instance),一个部件实例可以看作部件的代表,但并不是原部件的拷贝。实例一直和原部件保持关联,当原部件几何形状发生变化时,实例也发生相应变化。不能对部件实例直接编辑,一个装配模型可以包含一个部件的多个实例。所有装配模型中的实例都是该装配模型的特征体,在创建第一个实例时所生成的装配模型总体坐标系也是该装配模型的一个实例。 同一个部件中所有特征体在装配模块中对该部件建立实例时会形成一个整体,也即形成了装配模型中一个特征体。选择该实例时,该实例在装配之前原部件中所有特征体都被选择了,原部件中所有特征体在装配后形成了一个整体。
1学习笔记 学习资料:《Python语言在Abaqus中的应用》 20世纪60年代,Ray W. Clough教授在发表的论文《The Finite Element in Plane Stress Analysis》中首次提出有限单元法,此后,有限单元法的理论得到迅速发展,并广泛应用于各种力学问题和非线性问题,成为分析大型复杂工程结构的强有力手段。 Abaqus提供两种接口: 1)用户子程序接口(User Subroutine)。该接口使用Fortran语言进行开发,主要用于自定义本构关系、自定义单元等。常用的用户子程序包括(V)UMAT、(V)UEL、(V)FRIC、(V)DLOD等。 2)Abaqus脚本接口(Abaqus Scripting Interface)。该接口是在Python语言的基础上进行的定制开发,它扩充了Python的对象模型和数据类型,使Abaqus脚本接口的功能更加强大。一般情况下,Abaqus的脚本接口主要用于前处理(例如,快速建模)、后处理(例如,创建和访问输出数据库)、自定义模块(例如,自动后处理模块等)。 Python语言是公认的功能强大的面向对象的编程语言,Abaqus脚本接口在它的基础上又添加了许多数据类型和核心模块,因此功能更加强大。即便如此,Abaqus脚本接口仍然允许读者编写自己的模块或函数,以扩展其功能。 保留字 Python的保留字:Python语言的代码简洁,易于阅读,保留字相对较少;Python语言中不包含分号(;)、begin、end等标记,而是通过使用空格或制表键缩进的方式进行代码分隔。编写程序时,尽量不要选择保留字作为变量名、函数名等
ABAQUS用户子程序 当用到某个用户子程序时,用户所关心的主要有两方面:一是ABAQUS提供的用户子程序的接口参数。有些参数是ABAQUS传到用户子程序中的,例如SUBROUTINE DLOAD中的KSTEP,KINC,COORDS;有些是需要用户自己定义的,例如F。二是ABAQUS何时调用该用户子程序,对于不同的用户子程序ABAQUS调用的时间是不同的。有些是在每个STEP的开始,有的是STEP结尾,有的是在每个INCREMENT的开始等等。当ABAQUS调用用户子程序是,都会把当前的STEP和INCREMENT利用用户子程序的两个实参KSTEP和KINC传给用户子程序,用户可编个小程序把它们输出到外部文件中,这样对ABAQUS何时调用该用户子程序就会有更深的了解。 (子程序中很重要的就是要知道由abaqus提供的那些参量的意义,如下) 首先介绍几个子程序: 一.SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS, JLTYP,SNAME) 参数: 1.F为用户定义的是每个积分点所作用的荷载的大小; 2.KSTEP,KINC为ABAQUS传到用户子程序当前的STEP和INCREMENT值; 3.TIME(1),TIME(2)为当前STEP TIME和INCREMENT TIME的值; 4.NOEL,NPT为积分点所在单元的编号和积分点的编号; 5.COORDS为当前积分点的坐标; 6.除F外,所有参数的值都是ABAQUS传到用户子程序中的。 功能: 1.荷载可以被定义为积分点坐标、时间、单元编号和单元节点编号的函数。 2.用户可以从其他程序的结果文件中进行相关操作来定义积分点F的大小。 例1:这个例子在每个积分点施加的荷载不仅是坐标的函数,而且是随STEP变化而变化的。 SUBROUTINE DLOAD(P,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS, 1 JLTYP,SNAME) INCLUDE 'ABA_PARAM.INC' C DIMENSION TIME(2),COORDS(3) CHARACTER*80 SNAME PARAMETER (PLOAD=100.E4) IF (KSTEP.EQ.1) THEN !当STEP=1时的荷载大小 P=PLOAD ELSE IF (KSTEP.EQ.2) THEN !当STEP=2时的荷载大小 P=COORDS(1)*PLOAD !施加在积分点的荷载P是坐标的函数 ELSE IF (KSTEP.EQ.3) THEN !当STEP=3时的荷载大小 P=COORDS(1)**2*PLOAD ELSE IF (KSTEP.EQ.4) THEN !当STEP=4时的荷载大小 P=COORDS(1)**3*PLOAD ELSE IF (KSTEP.EQ.5) THEN !当STEP=5时的荷载大小 P=COORDS(1)**4*PLOAD END IF RETURN END UMAT 子程序具有强大的功能,使用UMAT 子程序: (1) 可以定义材料的本构关系,使用ABAQUS 材料库中没有包含的材料进行计算,扩充程序功能。 (2) 几乎可以用于力学行为分析的任何分析过程,几乎可以把用户材料属性赋予ABAQUS 中的任何单元; (3) 必须在UMAT 中提供材料本构模型的雅可比(Jacobian)矩阵,即应力增量对应变增量的变化率。 (4) 可以和用户子程序“USDFLD”联合使用,通过“USDFLD”重新定义单元每一物质点上传递到UMAT 中场变量的数值。 由于主程序与UMAT 之间存在数据传递,甚至共用一些变量,因此必须遵守有关UMAT 的书写格式,UMAT 中常用的变量在文件开头予以定义,通常格式为: SUBROUTINE UMA T(STRESS,STATEV,DDSDDE,SSE,SPD,SCD, 1 RPL,DDSDDT,DRPLDE,DRPLDT,
ABAQUS学习笔记 一.AQUS-.inp编码介绍 (一).ABAQUS头信息文件段(1-4) 1.*PREPRINT 输出求解过程所要求的信息(在dat文件中) ie:*PREPRINT, ECHO=YES, HISTORY=YES, MODEL=YES 2.*HEADING 标题输出文件(出现在POST/VIEW窗口中,且出现在结果输出文件中)ie:*HEADING STRESS ANAL YSIS FOR A PLATE WITH A HOLE 3.*RESTART 要求abaqus/standard输出其POST/view模块所需要的.res文件。其中的FREQ =?控制结果在每次迭代(或载荷步)输出的次数。 ie:*RESTART, WRITE, FREQ=1 4.*FILE FORMAT 要求abaqus/standard输出到.fil中的某些信息。它也用于post。对于在后处理中得到x-y形式的诸如应力-时间、应力-应变图有用! ie: *FILE FORMAT, ZERO INCREMENT (二).ABAQUS网格生成段 定义结点、单元,常用的命令有:结点定义(*NODE,*NGEN),单元定义(*ELEMENT,*ELGEN等)。 1.*NODE 定义结点,其格式为: *NODE 结点号,x轴坐标,y轴坐标,(z轴坐标) 2.*NGEN 在已有结点的基础上进行多个结点的生成,一般是在两结点间以某种方式(直线、圆)产生一定分布规律的结点。 如:*NGEN, LINE=C, NSET=HOLE, 119, 1919, 100, 101 在两结点(结点号为119,1919)间以圆弧形式生成多个结点,100为任意相邻结点的单元号增量,101为圆弧形成时圆心位置的结点(对于直线形式生成没有此结点)。所有这些生成的结点(包括119,1919)被命名成HOLE的集合(这样做的目的是以后的命令中使用到它,比如说对这些结点施加同等条件的边界条件或载荷等,HOLE就是这些结点的代称)。*NGEN使用的前提就是必须存在已有结点。 *NGEN, NSET=OUTER 131, 1031, 100 以线形式形成结点,结点号增量100,结点集合名为OUTER。 *NGEN, NSET=OUTER 1031, 1931, 100 同上生成结点,可以同上结点集合名,这样OUTER就包括这两次生成的所有结点 3.*NFILL 在如上生成的结点集(实际上,代表两条几何意义上的边界线)之间按一定规律(BIAS=?)填充结点。这样所有生成的结点构成一定形状的实体(面)。 如:*NFILL, NSET=PLATE, BIAS=0.8 HOLE, OUTER, 12, 1 以HOLE为第一条边界,OUTER为第二条边界(终止边),以从疏到密的规律(BIAS小于1)分布,其生成结点数在两内外对应结点间为12,1为每组结点号的增量。所有这些结点被置于PLATE的集合中。 下面以上面生成的结点来生成单元: 4.*ELEMENT
Abaqus学习笔记 Abaqus 使用日记Abaqus 标准版共有“部件(part)”、“材料特性(propoterty)”、“装配(assemble)”、“计算步骤(step)”、“交互(interaction)”、“加载(load)”、“单元划分(mesh)”、“计算(job)”、“后处理(visualization)”、“草图(sketch)”十大模块组成。建模方法:一个模型(model)通常由一个或几个部件(part)组成,部件又由一个或几个特征体(feature)组成,每一个部分至少有一个基本特征体(base feature),特征体可以是所创建的实体,如挤压体、切割挤压体、数据点、参考点、数据轴、数据平面、装配体的装配约束、装配体的实例等等。1.首先建立部件 (1)根据实际模型的尺寸决定部件的近似尺寸,进入绘图区。绘图区根据所输入的近似尺寸决定网格的间距,间距大小可以在edit菜单sketcher options 选项里调整。 (2)在绘图区分别建立部件中的各个特征体,建立特征体的方法主要有挤压、旋转、平扫三种。同一个模型中两个不同的部件可以有同名的特征体组成,也就是说不同部件中可以有同名的特征体,同名特征体可以相同也可以不同。部件的特征体包括用各种方法建立的基本特征体、数据点(datum point)、数据轴(datum axis)、数据平面(datum plane)等等。 (3)编辑部件可以用部件管理器进行部件复制,重命名,删除等,部件中的特征体可以是直接建立的特征体,还可以间接手段建立,如首先建立一个数据点特征体,通过数据点建立数据轴特征体,然后建立数据平面特征体,再由此基础上建立某一特征体,最先建立的数据点特征体就是父特征体,依次往下分别为子特征体,删除或隐藏父特征体其下级所有子特征体都将被删除或隐藏。特征体被删除后将不能够恢复,一个部件如果只包含一个特征体,删除特征体时部件也同时被删除。 2.建立材料特性 (1)输入材料特性参数弹性模量、泊松比等 (2)建立截面(section)特性,如均质的、各项同性、平面应力平面应变等等,截面特性管理器依赖于材料参数管理器 (3)分配截面特性给各特征体,把截面特性分配给部件的某一区域就表示该区域已经和该截面特性相关联 3.建立刚体 (1)部件包括可变形体、离散刚体和解析刚体三种类型,在创建部件时需要指定部件的类型,一旦建立后就不能更改其类型。采用旋转方式建立部件,在绘制轴对称部件的外形轮廓时不能超过其对称轴。
Abaqus学习笔记 Abaqus使用日记Abaqus标准版共有"部件(part)"、"材料特性(propoterty)5\ “装配(assemble)n、“计算步骤(step)”、“交互(interaction),5、“加载(load)”、“单元划分(mesh)”、“计算(job)”、“后处理(visualization)”、“草图(sketch)” 十大模块组成。建模方法:一个模型(model)通常由一个或儿个部件(part)组成,部件乂由一个或儿个特征体(feature)组成,每一个部分至少有一个基本特征体(base feature),特征体可以是 所创建的实体,如挤压体、切割挤压体、数据点、参考点、数据轴、数据平面、装 配体的装配约束、装配体的实例等等。 1.首先建立部件 (1)根据实际模型的尺寸决定部件的近似尺寸,进入绘图区。绘图区根据所输入的近似尺寸决定网格的间距,间距大小可以在edit菜单sketcher options 选项里调整。 (2)在绘图区分别建立部件中的各个特征体,建立特征体的方法主要有挤压、旋转、平扫三种。同一个模型中两个不同的部件可以有同名的特征体组成,也就是说不同部件中可以有同名的特征体,同名特征体可以相同也可以不同。部件的特征体包括用各种方法建立的基本特征体、数据点(datum point)>数据轴 (datum axis)、数据平面(datum plane)等等。 (3)编辑部件可以用部件管理器进行部件复制,重命名,删除等,部件中的特征体可以是直接建立的特征体,还可以间接手段建立,如首先建立一个数据点特征体,通过数据点建立数据轴特征体,然后建立数据平面特征体,再山此基础上建立某一特征体,最先建立的数据点特征体就是父特征体,依次往下分别为子特征体,删除或隐藏父特征体其下级所有子特征体都将被删除或隐藏。特征体被删除后将不能够恢复,一个部件如果只包含一个特征体,删除特征体时部件也同时被删除。 2.建立材料特性 (1)输入材料特性参数弹性模量、泊松比等 (2)建立截面(section)特性,如均质的、各项同性、平面应力平面应变等等,截面特性管理器依赖于材料参数管理器 (3)分配截面特性给各特征体,把截面特性分配给部件的某一区域就表示该区域已经和该截面特性相关联 3.建立刚体 (1)部件包括可变形体、离散刚体和解析刚体三种类型,在创建部件时需要指定部件的类型,一旦建立后就不能更改其类型。采用旋转方式建立部件,在绘制轴对称部件的外形轮廓时不能超过其对称轴。
衬砌开挖对上方框架结构的影响存在的一些相关问题(2019.9.1) 1.在构建框架结构后,在装配模块要对其进行布尔运算。步骤:构建主体,布 尔运算,赋予截面和材料属性,设定方向。 2.衬砌建模时用壳。 3.线单元的部件要设定n1方向。(n1方向为,大拇指指向杆的方向时,四指为 n1方向,四指弯曲90度为n2方向) 4.普通混凝土衬砌,在视频中只设置了密度和弹性模量。但是ABAQUS结构工 程分析中,在研究混凝土简支梁时,还设置了塑性(混凝土损伤塑性模型)? 这是根据研究重点不同而设置的吗? 5.土体属性设置了密度、弹性模量、塑性中的摩尔库伦模型(摩擦角、剪胀角、 粘聚力、abs=0) 6.衬砌与土体接触,讲师的意思是谁不动,谁设定为主面(还有准则为刚度大 的为主面) 7.框架结构与土体的接触Creat constrain中“嵌入Embedded region”。 8.施加荷载,重力gravity(整体施加)。然后地应力平衡,土、衬砌、框架分别 进行地应力平衡。原来好像记得。壳存在时无法平衡要先将其隐藏。首先stress 导入ODB,不行的话用gravity stress。 9.底部的边界条件设置为三个方向都限制。只设置Z方向计算有偏差。 10.衬砌第一步用model change杀死,开挖第一步完成后,在第二步激活。 基坑开挖与支撑 1.在基坑土体建立时,就将桩体位置土体挖除(Part左边倒数第二个图标,二维) 2.支撑杆由于是线,需要定义为梁单元,设置截面和方向等 3.二维土体采用壳单元,墙和锚杆是线单元,锚索定义截面尺寸和梁属性 4.将不同步骤建立的锚索分成不同集合。 5.第一步地应力平衡中没有支护结构(将其杀死);第二步添加桩;第三步添加 锚索(激活) 6.锚索与另外一个视频中的地连墙支护都是采用“嵌入”,耦合约束后再绑定 7.约束时,要将桩体的UR3固定住,否则会旋转不收敛 8.剑桥模型,要对整个土体施加“Voids ratio孔隙比”视频中为1。此外,剑桥 模型不能使用减缩积分,直接将勾选去掉即可。
inp编码介绍 (一).ABAQUS头信息文件段(1-4) 1.*PREPRINT 输出求解过程所要求的信息(在dat文件中) ie:*PREPRINT, ECHO=YES, HISTORY=YES, MODEL=YES 2.*HEADING 标题输出文件(出现在POST/VIEW窗口中,且出现在结果输出文件中) ie:*HEADING STRESS ANAL YSIS FOR A PLATE WITH A HOLE 3.*RESTART 要求abaqus/standard输出其POST/view模块所需要的.res文件。其中的FREQ=?控制结果在每次荷步)输出的次数。 ie:*RESTART, WRITE, FREQ=1 4.*FILE FORMAT 要求abaqus/standard输出到.fil中的某些信息。它也用于post。对于在后处理中得到x-y形式的时间、应力-应变图有用! ie: *FILE FORMA T, ZERO INCREMENT (二).ABAQUS网格生成段 定义结点、单元,常用的命令有:结点定义(*NODE,*NGEN),单元定义(*ELEMENT,*ELGEN等)。 1.*NODE 定义结点,其格式为: *NODE 结点号,x轴坐标,y轴坐标,(z轴坐标) 2.*NGEN 在已有结点的基础上进行多个结点的生成,一般是在两结点间以某种方式(直线、圆)产生一定分布如:*NGEN, LINE=C, NSET=HOLE, 119, 1919, 100, 101 在两结点(结点号为119,1919)间以圆弧形式生成多个结点,100为任意相邻结点的单元为圆弧形成时圆心位置的结点(对于直线形式生成没有此结点)。所有这些生成的结点(包括119,1919)被命名集合(这样做的目的是以后的命令中使用到它,比如说对这些结点施加同等条件的边界条件或载荷等,HOLE就的代称)。*NGEN使用的前提就是必须存在已有结点。 *NGEN, NSET=OUTER 131, 1031, 100 以线形式形成结点,结点号增量100,结点集合名为OUTER。 *NGEN, NSET=OUTER 1031, 1931, 100 同上生成结点,可以同上结点集合名,这样OUTER就包括这两次生成的所有结点 3.*NFILL 在如上生成的结点集(实际上,代表两条几何意义上的边界线)之间按一定规律(BIAS=?)填充结有生成的结点构成一定形状的实体(面)。 如:*NFILL, NSET=PLATE, BIAS=0.8 HOLE, OUTER, 12, 1 以HOLE为第一条边界,OUTER为第二条边界(终止边),以从疏到密的规律(BIAS小其生成结点数在两内外对应结点间为12,1为每组结点号的增量。所有这些结点被置于PLATE的集合中。 下面以上面生成的结点来生成单元: 4.*ELEMENT 定义单元所使用的类型(TYPE=?),然后另行定义通过联结结点形成单元,其结点数目依靠单元类型而变。*ELEMENT, TYPE=CPS4 //采用四单元的平面应力单元 19, 119, 120, 220, 219 //定义顺序:单元号,以逆时针方向形成单元的各结点号 (三)ABAQUS单元 注意:分析前要选择合适的元素,这时要考虑的问题就是:使用什么样类型的单元?有限元的基本思路就是将实体离散化,实际结果是将众多离散分析结果的集合,这似乎有点像积分的概念。选择元素种类最重要考虑的是分象,满足必要的准确度基础上去掉不必要的细节与准确度。是选择1-D, 2-D or 3-D单元、用于何种分析的单元、
在ABAQUS 中必须用真实应力和真实应变定义塑性.ABAQUS 需要这些值并对应地在输入文件中解释这些数据。 然而,大多数实验数据常常是用名义应力和名义应变值给出的。这时,必须应用公式将塑性材料的名义应力(变)转为真实应力(变)。 考虑塑性变形的不可压缩性,真实应力与名义应力间的关系为: 00l A lA =, 当前面积与原始面积的关系为: 00l A A l = 将A 的定义代入到真实应力的定义式中,得到: 00 ()nom F F l l A A l l σσ=== 其中0 l l 也可以写为1nom ε+。 这样就给出了真实应力和名义应力、名义应变之间的关系: (1)nom nom σσε=+ 真实应变和名义应变间的关系很少用到,名义应变推导如下: 000 1nom l l l l l ε-==- 上式各加1,然后求自然对数,就得到了二者的关系: ln(1)nom εε=+ ABAQUS 中的*PLASTIC 选项定义了大部分金属的后屈服特性。ABAQUS 用连接给定数据点的一系列直线来逼近材料光滑的应力-应变曲线。可以用任意多的数据点来逼近实际的材料性质;所以,有可能非常逼真地模拟材料的真实性质。在*PLASTIC 选项中的数据将材料的真实屈服应力定义为真实塑性应变的函数。选项的第一个数据定义材料的初始屈服应力,因此,塑性应变值应该为零。 在用来定义塑性性能的材料实验数据中,提供的应变不仅包含材料的塑性应变,而是包括材料的总体应变。所以必须将总体应变分解为弹性和塑性应变分量。弹性应变等于真实应力与杨氏模量的比值,从总体应变中减去弹性应变,就得到了塑性应变,其关系为: /pl t el t E ε εεεσ=-=- 其中pl ε是真实塑性应变,t ε是总体真实应变,el ε是真实弹性应变。
1.接触中设置Adjust的理解: 这个命令主要还是用来初始化接触的。在分析开始之前,调整接触面中节点的初始位置,且不产生任何应力和应变。在分析过程中,由于残留的初始过盈引起的应变将被施加在接触面上。 模型的尺寸往往会存在数值误差,所以设置一个位置误差限度,用来调整从面节点的初始坐标,ADJUST=位置误差限度,其含义是:如果从面节点与主面的距离小于此限度,ABAQUS将调整这些节点的初始坐标,使其与主面的距离为0.这种调整不产生任何压力、应力、应变。 Explicit不允许接触表面的初始过盈,分析开始前,接触面上的节点将被自动调整,以删除任何初始过盈,在随后的分析中,这样的调整将引起应变。 2.使用INTERFERENCE(干涉)来定义过盈接触: Edit Interaction》底部Interference fit 负值表示过盈量,正值表示间隙量。类似于载荷,只能在后续分析步中定义,不能在初始分析步中定义。
3.CLERENCE(间隙)可以定义两个接触面之间的初始过盈量和间隙 量,它只适用于小滑移,并且不需要使用ADJUST来调整从面节点的初始位置。 4.特征:记录了设计目的,并包含几何信息,同时也是管理几何体的行为的规则。ABAQUS中导入的几何体是没用特征的,要删除不重要的细节。 5.View=》ODB Display Options =》Sweep and Extrude 6.CPRESS接触压强 7.COPEN从面上节点与主面的距离
8.ALE自适应网格:Step=>Other=>Adaptive Mesh Domain 9.计算代价估算:隐式:自由度数目的平方 显式:自由度正比 10.软接触:用指数或者表格形式表达的应力-距离关系 11.传说ABAQUS默认的幅值曲线是从1下降到0的。是吗?我觉得 6.9的版本好像不是这样。给一个棱柱施加扭转力矩,载荷采用 Ramp,变形是逐步增大的。 12.Visualization》Tools》Job Diagnostics 13.修改严重不连续迭代尝试次数: Step>Other>General Solution Controls> Edit>Specify>Time Incrementation>More>Is>12 14.接触问题中,90度圆角至少划分是个单元。 15.使用自动过盈接触限度来解决颤振收敛问题: Interaction> Contact Controls> Automatic Overclosure tolerance 再在Edit Interaction对话框中的最下角的Contact Controls 设置为已定义的接触控制名称。 16.PEEQ等效塑性应变,是塑性应变的积累。PEMAG塑性应变量, PE塑性应变分量。 17.子模型的定义方法: Model>Edit Attributes>Submodel>Read data from job
Abaqus 中安装intel visual fortran compiler v11.1不必通过子程序验证又能正常使用子程序的方法1.Pasted from
入门篇: 一.建模及后处理的一般步骤 (1)ABAQUS的一个推荐顺序 按照模块Module依次选择:Part(部件)、Property(特性,包括材料及截面)、Assembly(实体的组合)、Step(分析步,初始步Initial自动生成,后续步需自己定义,如加载步Apply load)、Interaction(相互作用)、Load(边界条件及载荷)、Mesh(网格划分及其设置)、Job(分析作业)、Visualization(后处理); 这种建模顺序是将数值模型定义在实体模型上,而非网格的单元和节点上,这样即便修改网格也无需重新定义材料及边界条件等模型参数。处理复杂问题时建议简单粗划分网格,再适当地局部细化网格。 一般创建Part常常还会用到绘图(Sketch)的操作。 以上各步骤间建议即时存档(Save),防止出现文件信息的丢失。 (2)一些快捷操作: 鼠标中键:确认结束某一步操作时,一般总在视图区中点击鼠标中键,或点击窗口底部的Done进行。 CTRL+ALT+鼠标右键/左键/中键:缩放/旋转/平移 (3)建模及分析的一些细节: 材料定义(Create Material):对话框中添加各种特性,如力学特性Mechanical→弹性Elasticity →弹性Elastic,在进行具体的设置,如弹性模量Young's Modulus及泊松比Poisson’s Ratio 的数值。 截面的定义及赋予:一般截面属性及赋予是分开的,原因是材料属性不直接赋予在单元或几何实体上,而是先在截面属性(Section)中定义材料特性,再为每个部件赋予相应的截面属性。 载荷模块(Load):这里包括加载(Create Load)及边界条件(Create Boundary Condition)等功能。注意荷载、边界条件所对应的分析步(一般加载放在加载步,边界条件步放在初始步),加载或赋予边界条件时应注意选中的区域(线、面)会有红色高亮提示,确认不要错误地定义荷载或边条的作用位置。 网格划分:一般为布种(全局种子或局部种子,后者可以在一定区域内按尺寸及数量更精确地布种),控制参数(Assign Mesh Controls),确定单元类型及划分网格四个分步骤。 控制参数主要包括划分方法(Free、Structured、Sweep等),单元形状(如二维结构化划分里包括Quad、Quad-dominated、Tri) 可以通过窗口上方的Object选框切换对象:组合体(Assembly)及部件(Part) 单元类型设置包括几何阶次(Geometric Order)等,如存在应力集中问题可使用二次单元(Quadratic)提高精度。设置一定参数后可直接在对话框内看到当前的单元类型(不是自选单元,而是根据设定自动生成单元)。 提交分析作业:通过Job Manager中创建(Create),设置分析参数,点击OK后即创建分析,然后选择待分析的作业,点击Submit,依次提示Submitted、Running、Completed,然后模型分析完毕,点击Job Manager 对话框的Results直接读取分析结果,进入后处理(Visualization)模块。 若Status提示为Aborted,表明模型存在问题,分析已终止,可以点击Monitor进行查看错误信息,对之前的建模步骤查错和修改,更正错误后再提交进行分析。
Abaqus标准版共有“部件(part)”、“材料特性(propoterty)”、“装配(assemble)”、“计算步骤(step)”、“交互(interaction)”、“加载(loa d)”、“单元划分(mesh)”、“计算(job)”、“后处理(visualization)”、“草图(sketch)”十大模块组成。 建模方法: 一个模型(model)通常由一个或几个部件(part)组成,“部件”又由一个或几个特征体(feature)组成,每一个部分至少有一个基本特征体(base featur e),特征体可以是所创建的实体,如挤压体、切割挤压体、数据点、参考点、数据轴,数据平面,装配体的装配约束、装配体的实例等等。 1.首先建立“部件” (1)根据实际模型的尺寸决定部件的近似尺寸,进入绘图区。绘图区根据所输入的近似尺寸决定网格的间距,间距大小可以在edit菜单sketcher options选项里调整。 (2)在绘图区分别建立部件中的各个特征体,建立特征体的方法主要有挤压、旋转、平扫三种。同一个模型中两个不同的部件可以有同名的特征体组成,也就是说不同部件中可以有同名的特征体,同名特征体可以相同也可以不同。部件的特征体包括用各种方法建立的基本特征体、数据点(datum point)、数据轴(d atum axis)、数据平面(datum plane)等等。 (3)编辑部件可以用部件管理器进行部件复制,重命名,删除等,部件中的特征体可以是直接建立的特征体,还可以间接手段建立,如首先建立一个数据点特征体,通过数据点建立数据轴特征体,然后建立数据平面特征体,再由此基础上建立某一特征体,最先建立的数据点特征体就是父特征体,依次往下分别为子特征体,删除或隐藏父特征体其下级所有子特征体都将被删除或隐藏。××××特征体被删除后将不能够恢复,一个部件如果只包含一个特征体,删除特征体时部件也同时被删除××××× 2.建立材料特性 (1)输入材料特性参数弹性模量、泊松比等 (2)建立截面(section)特性,如均质的、各项同性、平面应力平面应变等等,截面特性管理器依赖于材料参数管理器 (3)分配截面特性给各特征体,把截面特性分配给部件的某一区域就表示该区域已经和该截面特性相关联 3.建立刚体
ABAQUS 6.13版本安装完成后将安装目录下的6.13-1/code/bin文件中的mkl_avx2.dll 改名为mkl_avx2.dll.13.0.0.1 功能模块: 1.part:创建或者接口导入单独部件 2.Property:定义截面集合形状与材料性质,并赋予部件的某一部分 3.Assembly:将单独的部件定位于总体坐标系中构成装配件 4.Step:模拟过程,将其生成和构建一个个分析步,可以在分析步之间更改输出变量 5.Interaction:模型各区域之间或者模型区域与环境之间的相互作用,包括热力,力,约束必须规定相互作用是在哪个分析步起作用 6.Load:指定载荷,边界条件和场变量,也是需要指定具体的分析步 7.Mesh:网格自动划分和控制工具 8.Job:计算 9.Opti:优化模型(依据优化目标与限制条件) 10.Visul:模型与分析结果的图像显示 11.Sktech:(二维轮廓图)定义二维平面部件,梁,剖面,或者依据二维拉伸生成三维图 Cell是abaqus/CAE中能单独赋予材料和截面属性的最小单元,建模时,可以根据需要选择在基础上模型上增加一个cell来单独进行截面特性分配;若分配时发现需要单独分配截面特性的部分没有分离出来,可以选用工具区中适当的Partition 工具进行部件分割。 如何将一个组合的构建各个部分单独分配截面特性: abaqus梁建模: 梁建模只需要用到三维空间里的线,再对线赋予对应的梁截面特性 1.part: 创建三维线 2.property: 创建梁(beam)截面,在beam shape 中可以选定系统给定的截面特性:rectangular; I; L; box; pipe; circular 等 3.将截面特性指派给选中的区域梁(done 或鼠标中键确认) 4.指派梁的方向(assign beam oritation),done, 查看具体的梁截面三维图在View-Part Display Options-勾选render beam peofiles 5.n1 方向为梁截面宽度方向,n2 为梁截面高度方向,只需确定n1方向:选定需要修改的梁(按住鼠标左键拖拉矩形区域通过该梁),修改n1向量方向(0,0,-1)>(0,-1,0)
1、abaqus中的力载荷 集中力concentrated force、压强pressure(垂直于表面)、表面分布力surface traction (设定沿着某方向) pressure只能施加在面上(几何的面,单元的面),为垂直于表面的分布力; surface traction只能施加在面上(几何的面,单元的面),为沿着某一方向的分布力;concentrated force只能施加在点上(几何的点,节点),要使得集中力产生的效果等同于分布力,则需要将集中力施加在参考点上,然后将参考点与作用面上的节点进行耦合约束coupling(distributed coupling),而不要直接施加在节点上。 一般,如果不要求等效均布力,则集中力最好施加在几何的点上。确实需要施加节点力,则施加在节点上。 对于有限元软件,所有的力载荷本质上都由程序处理成节点力。 2、abaqus计算热电耦合出现Too many attempts made for this increment (1)调整一下计算载荷施加的速度或者调整载荷大小,要么把计算步长设置的小一点,尝试次数设的多一点。这个提示是说计算的过程中直到设定的尝试次数极限仍然求解失败。(2)分析步主要有初始分析步和后续分析步,每个分析步可以用来描述一个分析过程,例如在后续分析步中施加不同荷载,在初始分析步中施加边界条件等。 增量步是在分析步里面根据模型计算收敛情况设置的,简单模型可以设置较少的增量步,并可使初始增量为1;复杂模型设置多一点增量步,并减少初始增量值。超过设置的允许增量步数,则计算停止。 (3)检查模型,是否存在刚体位移,过约束,接触定义不当等问题 (4)分别建立四个边界条件,BC-1,BC-2,BC-3,BC-4,每一个边界条件定义板的一边固结的支承条件就行了。之前是建立了一个BC-1,四边的约束都定义在BC-1里面,就算不下去了,不清楚原因。仅供参考学习。 (5)1.可以把初始增量步最小增量步调小,最大增量步的数目调大。这个方法一般很难解决2.用不同的方法重新划分网格,并不是越细越好。3.接触有问题,适当配合方式和接触对参数。4.如果施加载荷比较复杂的话,可以细分分析步。你多尝试修改几次应该能把问题解决,多交流。
Abaqus/python入门体会(初稿) 长安大学 姜峰林 flin55@https://www.wendangku.net/doc/f515065088.html, 2009.2.1 #=========================================================== 自己的论文要用到有限元进行数值模拟分析,以前都用ansys计算,可ansys 中岩土的本构模型只有D-P模型,无法准确的反映土的硬化/软化性质,模拟计算出的结果因此也和实际差别很大。Abaqus有着丰富的材料模型,超强的非线性分析能力,岩土的模型也很多,因此才转学Abaqus。Abaqus的cae建模功能还是很好的,但科研课题一般都要进行参数分析,采用cae的建模方法有些不切实际,学了没几天就放弃cae开始学习inp,也是学了一阵子才知道inp不能建立实体模型,只能直接建节点和单元。复杂的模型inp也无法建立,但采用Python建模就可以解决这个问题。 由于Abaqus的学习资料不多,过了好些日子才知道Abaqus也可以采用Python语言进行建模计算,只是比Ansys的Apdl语言复杂得多,并且除了手册上的Script资料之外,没有较为系统的教程,刚一接触真是让人头痛。通过查看Simwe论坛上关于Python的帖子,和论坛朋友的帮助,自己在慢慢积累,现在对Python有了一点点了解,算是入了个门。 接触Abaqus也没多久,对python更是一知半解,绝大多数地方根本都不清楚,抽空写一点认识体会主要是给像自己一样刚学习Abqus Python的朋友,能少走一些弯路,节约一些时间。同时希望大家批评指正、共同讨论、补充。#-------------------------------------------------------------------------------------------------- 学习Abaqus/Python基础:Abaqus的cae建模有比较全面的认识;了解一些Python语法知识(大家都不会有太多时间单独学习Python语言本身,只需要有概念了解即可,不懂的地方可以随时查询Python script手册)Abaqus/Python学会使用不太难,可要精通应用还是要付出一定的劳动。大家所分析的课题专业不同,方向也千差万别,所用到的Abaqus的功能也就有很大的差别,能对自己的工作领域熟练应用就算成功。Abaqus毕竟只是软件,如何考虑专业知识成功建模才是最困难的。 #-------------------------------------------------------------------------------------------------- 1. Python与Abaqus 2. Abaqus/Python结构 3.模型参数分析技巧 4.几个命令的体会 5.一个Abaqus/Python例子 #--------------------------------------------------------------------------------------------------