ABAQUS用户子程序小例子
abaqus2用户单元子程序
![abaqus2用户单元子程序](https://img.taocdn.com/s3/m/2ccf1a0116fc700abb68fcfc.png)
20 ABAQUS用户单元子程序(UEL)在这一章中将列举两个在这些年里发展过的ABAQUS/Standard用户单元子程序(UEL)。
第一个例子是一个非线性的索单元,我们的目的是通过这个比较简单的例子让读者了解用户单元子程序的基本开发过程;第二个例子是一个用于计算应变梯度理论的单元,应变梯度是当今比较热点的一个科研前沿问题,有各种理论,我们为了验证新的理论,需要数值结果与实验对照来进行评价,整个例子的目的是通过它说明用户子单元可以求解的问题范围很广,但是由于内容比较艰深,程序也很长,所以这个例子我们并没有给出最后的全部程序。
另外,到目前为止,ABAQUS还只有隐式求解器ABAQUS/Standard支持用户自定义单元,而显式求解器ABAQUS/Explicit中还不支持这一功能。
20.1 非线性索单元20.1.1 背景钢索斜拉桥和斜拉索结构广泛应用于土木工程建筑上。
索力的计算分析是设计和施工的关键环节。
清华大学工程力学系在采用ABAQUS进行荆沙长江斜拉桥的计算机仿真分析(这个项目我们已在第15章“ABAQUS在土木工程中的应用(一)——荆州长江大桥南汊斜拉桥结构三维仿真分析”中讨论过)时,也曾进行了自行建立索单元的尝试。
本节介绍的就是这方面的工作。
香港理工大学土木与结构工程系采用ABAQUS有限元软件进行计算,完成了香港Ting Kau斜拉桥和Tsing Ma悬索桥的结构计算和分析。
对于钢索计算,他们采用梁单元进行模拟。
由于梁单元含有弯曲刚度,计算的高阶频率值偏高,周期较低。
一般假设索是单向受拉力的构件。
随着应变的非线性增加,索力呈非线性增加。
尽管ABAQUS单元库中有500个以上的单元类型,但是,还没有索单元。
本文发展了三维非线性索单元模型,形成ABAQUS的用户单元子程序,可以利用ABAQUS输入文件调入到具体的分析中。
通过静态和动态例题的计算比较,索单元工作良好。
20.1.2 基本公式在三维索单元计算中,如图20-1所示,坐标x 和位移u 的变量表达式为:ij ji i j ji u u u x x x -=-= (x,y,z) (u,v,w) (20-1)应变的公式为:()⎥⎦⎤⎢⎣⎡+++++=222211ji ji ji ji ji ji ji ji ji w v u w z v y u x L ε (20-2)公式(20-2)中,L 为索的长度,索的张力为:0N AE N +=ε(20-3)在公式(20-3)中,A 为截面面积,E 为弹性模量,N 0为初始张力。
ABAQUS用户子程序
![ABAQUS用户子程序](https://img.taocdn.com/s3/m/550ee9fd700abb68a982fb0b.png)
将没有任何影响的。也就是说,ABAQUS中调用的时候,实际并不存在这个错误。FT,忘了
说一句了,你把ABA_PARAM.INC.dp或ABA_PARAM.INC.sp拷到你的程序工作空间后 ,应该将
几个关于子程序的问题及相应解答
Q 本人在用umat作本构模型时,
static,
1,500,0.000001,0.1 此时要求的增量步很多,即每次增量要很小,
static
1,500 时,在弹性向塑性过度时,出现错误,增量过大,出现尖点.
A YOU CAN TRY AS FOLLOWS
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是坐标的函数
态变量矩阵中的任何矢量或者张量都必须通过旋转来考虑材料的刚体运动。
状态变量矩阵的维数,等于关键字“DEPVAR”定义的数值。状态变量矩阵的维数通过ABAQUS 输入文件中的关键字“DEPVAR”定义,关键字下面数据行的数值即为状态变量矩阵的维数。
材料常数的个数,等于关键字“USER MATERIAL”中“CONSTANTS”常数设定的值。
STRESS (NTENS)
应力张量矩阵,对应NDI 个直接分量和NSHR 个剪切分量。在增量步的开始,应力张量矩阵中的数值通过UMAT 和主程序之间的接口传递到UMAT 中,在增量步的结束UMAT 将对应力张量矩阵更新。对于包含刚体转动的有限应变问题,一个增量步调用UMAT 之前就已经对应力张量的进行了刚体转动,因此在UMAT 中只需处理应力张量的共旋部分。UMAT 中应力张量的度量为柯西(真实)应力。
Abaqus材料用户子程序UMAT基础知识及手册例子完整解释
![Abaqus材料用户子程序UMAT基础知识及手册例子完整解释](https://img.taocdn.com/s3/m/e2cd26cf89eb172ded63b7f5.png)
1、为何需要使用用户材料子程序(User-Defined Material, UMAT )?很简单,当ABAQUS 没有提供我们需要的材料模型时。
所以,在决定自己定义一种新的材料模型之前,最好对ABAQUS 已经提供的模型心中有数,并且尽量使用现有的模型,因为这些模型已经经过详细的验证,并被广泛接受。
UMAT 子程序具有强大的功能,使用UMAT 子程序:(1)可以定义材料的本构关系,使用ABAQUS 材料库中没有包含的材料进行计算,扩充程序功能。
(2) 几乎可以用于力学行为分析的任何分析过程,几乎可以把用户材料属性赋予ABAQU S 中的任何单元。
(3) 必须在UMAT 中提供材料本构模型的雅可比(Jacobian )矩阵,即应力增量对应变增量的变化率。
(4) 可以和用户子程序“USDFLD ”联合使用,通过“USDFLD ”重新定义单元每一物质点上传递到UMAT 中场变量的数值。
2、需要哪些基础知识?先看一下ABAQUS 手册(ABAQUS Analysis User's Manual )里的一段话:Warning: The use of this option generally requires considerable expertise(一定的专业知识). The user is cautioned that the implementation (实现) of any realistic constitutive (基本) model requires extensive (广泛的) development and testing. Initial testing on a single eleme nt model with prescribed traction loading (指定拉伸载荷) is strongly recommended. 但这并不意味着非力学专业,或者力学基础知识不很丰富者就只能望洋兴叹,因为我们的任务不是开发一套完整的有限元软件,而只是提供一个描述材料力学性能的本构方程(Constitutive equation )而已。
abaqus子程序umat编写hashin准则
![abaqus子程序umat编写hashin准则](https://img.taocdn.com/s3/m/58077c57b6360b4c2e3f5727a5e9856a5712266b.png)
Hashin准则是一种常用的材料失效准则,用于模拟复合材料的损伤和失效。
在ABAQUS中,你可以使用UMAT(用户自定义材料行为)子程序来实现Hashin准则。
以下是一个简单的UMAT子程序示例,用于实现Hashin准则:```c#include "umat.h"#define N 10void UMAT(const int *nIntPoints, const int *nExtPoints, const double *dS, const double *dSP, const double *dStrain, const double*dStress, double *dStrainInc, double *dStressInc, double *dStress历史, int *nFail){double S[N], S0[N];int i;double E1 = dStrain[0];double E2 = dStrain[1];double E3 = dStrain[2];double sigma1 = dStress[0];double sigma2 = dStress[1];double sigma3 = dStress[2];double sigma历史1 = dStress历史[0];double sigma历史2 = dStress历史[1];double sigma历史3 = dStress历史[2]; double m1 = 0.0;double m2 = 0.0;double m3 = 0.0;double n1 = 0.0;double n2 = 0.0;double n3 = 0.0;for (i=0; i<N; i++) {S[i] = sqrt(dS[i]);S0[i] = sqrt(dSP[i]);m1 += S[i] / sigma1;m2 += S[i] / sigma2;m3 += S[i] / sigma3;n1 += S0[i] / sigma历史1;n2 += S0[i] / sigma历史2;n3 += S0[i] / sigma历史3;}if (m1 < 1.0e-6) m1 = 1.0e-6;if (m2 < 1.0e-6) m2 = 1.0e-6;if (m3 < 1.0e-6) m3 = 1.0e-6;if (n1 < 1.0e-6) n1 = 1.0e-6;if (n2 < 1.0e-6) n2 = 1.0e-6;if (n3 < 1.0e-6) n3 = 1.0e-6;double sigmaT = sqrt(m1*m2*m3);double sigmaLT = sqrt(n1*n2*n3);if (sigmaT < sigmaLT) {*nFail = 1; // 损伤发生,将nFail设为1表示材料失效return; // 退出UMAT函数,不再进行应力和应变增量计算} else {*nFail = 0; // 没有损伤发生,将nFail设为0表示材料仍然有效return; // 退出UMAT函数,继续进行应力和应变增量计算}}```以上代码实现了Hashin准则的简单形式,即当总剪切应力小于长纤维方向的纤维应力时,复合材料发生损伤。
ABAQUS_Standard用户材料子程序实例
![ABAQUS_Standard用户材料子程序实例](https://img.taocdn.com/s3/m/e3e24f3610661ed9ad51f3c8.png)
ABAQUS/Standard 用户材料子程序实例-Johnson-Cook 金属本构模型卢剑锋 庄茁* 张帆清华大学工程力学系 北京 100084摘要:用户材料子程序是ABAQUS 提供给用户定义自己的材料属性的Fortran 程序接口,使用户能使用ABAQUS 材料库中没有定义的材料模型。
ABAQUS 中自有的Johnson-Cook 模型只能应用于显式ABAQUS/Explicit 程序中,而我们希望能在隐式ABAQUS/Standard 程序中更精确的实现本构积分,而且应用Johnson-Cook 模型的修正形式。
这就需要通过ABAQUS/Standard 的用户材料子程序UMAT 编程实现。
在UMAT 编程中使用了率相关塑性理论以及完全隐式的应力更新算法。
1 Johnson-Cook 强化模型简介Johnson-Cook (JC )模型用来模拟高应变率下的金属材料。
JC 强化模型表示为三项的乘积,分别反映了应变硬化,应变率硬化和温度软化。
这里使用JC 模型的修正形式:()()*01ln 11n mA B C T εσεε=+++−&& 并使参考应变率01ε=&,这样公式中的A 即为材料的静态屈服应力。
公式中包含,,,,A B n C m 五个参数,需要通过实验来确定。
2 ABAQUS 用户材料子程序用户材料子程序(User-defined Material Mechanical Behavior ,简称UMAT )通过与ABAQUS 主求解程序的接口实现与ABAQUS 的数据交流。
在输入文件中,使用关键字“*USER MATERIAL”表示定义用户材料属性。
子程序概况与接口UMAT 子程序具有强大的功能,使用UMAT 子程序:(1) 可以定义材料的本构关系,使用ABAQUS 材料库中没有包含的材料进行计算,扩充程序功能。
(2)几乎可以用于力学行为分析的任何分析过程,几乎可以把用户材料属性赋予ABAQUS中的任何单元;(3)必须在UMAT中提供材料本构模型的雅可比(Jacobian)矩阵,即应力增量对应变增量的变化率。
浅谈ABAQUS用户子程序
![浅谈ABAQUS用户子程序](https://img.taocdn.com/s3/m/93963d8fa0116c175f0e484b.png)
Home浅谈ABAQUS用户子程序李青清华大学工程力学系摘要本文首先概要介绍了ABAQUS的用户子程序和应用程序,然后从参数,功能两方面详细论述了DLOAD, UEXTERNALDB, URDFIL三个用户子程序和GETENVVAR,POSFIL,DBFILE三个应用程序,并详细介绍了ABAQUS的结果文件(.FIL)存储格式。
关键字ABAQUS,用户子程序,应用程序,结果文件一、前言:ABAQUS为用户提供了强大而又灵活的用户子程序接口(USER SUBROUTINE)和应用程序接口(UTILITY ROUTINE)。
ABAQUS 6.2.5一共有42个用户子程序接口,13个应用程序接口,用户可以定义包括边界条件、荷载条件、接触条件、材料特性以及利用用户子程序和其它应用软件进行数据交换等等。
这些用户子程序接口使用户解决一些问题时有很大的灵活性,同时大大的扩充了ABAQUS的功能。
例如:如果荷载条件是时间的函数,这在ABAQUS/CAE 和INPUT 文件中是难以实现的,但在用户子程序DLOAD中就很容易实现。
二.在ABAQUS中使用用户子程序ABAQUS的用户子程序是根据ABAQUS提供的相应接口,按照FORTRAN语法用户自己编写的代码。
在一个算例中,用户可以用到多个用户子程序,但必须把它们放在一个以.FOR为扩展名的文件中。
运行带有用户子程序的算例时有两种方法,一是在CAE中运行,在EDIT JOB菜单的GENERAL子菜单的USER SUBROUTINE FILE对话框中选择用户子程序所在的文件即可;另外是在ABABQUS COMMAND用运行,语法如下:ABAQUS JOB=[JOB] USER¡[.FOR]¡C用户在编写用户子程序时,要注意以下几点:1.用户子程序不能嵌套。
即任何用户子程序都不能调用任何其他用户子程Home序,但可以调用用户自己编写的FORTRAN子程序和ABAQUS应用程序。
ABAQUS用户子程序学习小结
![ABAQUS用户子程序学习小结](https://img.taocdn.com/s3/m/980ee8af312b3169a551a4cc.png)
ABAQUS用户子程序学习小结1 FORTRAN语言中的“I-N规则”:I、J、K、L、M、N开头的为整型变量,其他开头为实型变量;1.1 Fortran书写格式:1-5列:标号区;6列:续写标号区,一般就写"1";7-72列:语句区,本区内空格无效;注释行以C开头,本行内书写格式无要求;参考周煦《FORTRAN77结构化程序设计》,中国科学技术出版社,1995,38-40页 2 DIMENSION COORDS(3)表示声明一个含3个元素的数组,下标分别为1、2、3,访问形式为COORDS(n),n为1,3;3 子程序(*.for)文件中如何输出调试信息:WRITE(6,*)'COORDS(1)',COORDS(1),在*.dat文件中可看到输出,如果希望WRITE输出到msg文件中,则写为WRITE(7,*)'COORDS...;4 用户子程序DLOAD中COORDS数组的含义:COORDS(1)也是一个数组,存贮单元集合中所有单元积分点的X坐标,COORDS(2)存贮Y坐标,相应INP文件中的写法为:*DLOADPY,PYNU其中PY为单元集合名称,定义方法为:*Elset, elset=BEAM, generate1, 5, 1...*ELSET,ELSET=PYBEAM5 DLOAD中F的定义方法:F只有定义在单元积分点上才有效,例如:F=1.0*COORDS (1)附一个简单实例:beam.inp文件:*Heading** Job name: Job-1 Model name: beam *Preprint, echo=NO, model=NO, history=NO, contact=NO**** PARTS***Part, name=PART-1*End Part**** ASSEMBLY***Assembly, name=Assembly***Instance, name=PART-1-1, part=PART-1 *Node1, 0., 0.2, 20., 0.3, 40., 0.4, 60., 0.5, 80., 0.6, 100., 0. *Element, type=B311, 1, 22, 2, 33, 3, 44, 4, 55, 5, 6*Elset, elset=BEAM, generate1, 5, 1** Region: (Section-1-BEAM:BEAM), (Beam Orientation:BEAM)** Section: Section-1-BEAM Profile: Profile-1*Beam Section, elset=BEAM, material=STEEL, temperature=GRADIENTS, section=RECT0.2, 5.0.,0.,-1.*End Instance*Nset, nset=ENDS, instance=PART-1-11, 6*Nset, nset=_M4, internal, instance=PART-1-16,*Nset, nset=_M5, internal, instance=PART-1-11,*End Assembly**** MATERIALS***Material, name=STEEL*Elastic210000., 0.3*ELSET,ELSET=PYBEAM**** BOUNDARY CONDITIONS**** Name: Disp-BC-1 Type: Symmetry/Antisymmetry/Encastre*Boundary_M4, ENCASTRE** ---------------------------------------------------------------- **** STEP: Step-1***Step, name=Step-1*Static**** LOADS**** Name: CFORCE-1 Type: Concentrated force*DLOADPY,PYNU**** OUTPUT REQUESTS****** FIELD OUTPUT: F-Output-1 ***Output, field, variable=PRESELECT **** FIELD OUTPUT: F-Output-2 ***Output, field*Element OutputSF,**** HISTORY OUTPUT: H-Output-1 ***Output, history*Node Output, nset=ENDSCF1, CF2, CF3, CM1, CM2, CM3, RF1, RF2 RF3, RM1, RM2, RM3, U1, U2, U3, UR1 UR2, UR3*El Print, freq=999999*Node Print, freq=999999*End Stepbbb.for文件SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,1 JLTYP,SNAME) CINCLUDE 'ABA_PARAM.INC'CDIMENSION TIME(2), COORDS (3)CHARACTER*80 SNAMEWRITE(6,*)'COORDS(3)',COORDS(3)F=1.0*COORDS (1)RETURNEND运行方法:在Abaqus Command提示符后输入:abaqus job=beam user=bbb interactive 子程序如下~如何编写,调研子程序,子程序USDFLD中给固化度赋了一个初值1×10-4。
浅谈ABAQUS用户子程序
![浅谈ABAQUS用户子程序](https://img.taocdn.com/s3/m/93963d8fa0116c175f0e484b.png)
Home浅谈ABAQUS用户子程序李青清华大学工程力学系摘要本文首先概要介绍了ABAQUS的用户子程序和应用程序,然后从参数,功能两方面详细论述了DLOAD, UEXTERNALDB, URDFIL三个用户子程序和GETENVVAR,POSFIL,DBFILE三个应用程序,并详细介绍了ABAQUS的结果文件(.FIL)存储格式。
关键字ABAQUS,用户子程序,应用程序,结果文件一、前言:ABAQUS为用户提供了强大而又灵活的用户子程序接口(USER SUBROUTINE)和应用程序接口(UTILITY ROUTINE)。
ABAQUS 6.2.5一共有42个用户子程序接口,13个应用程序接口,用户可以定义包括边界条件、荷载条件、接触条件、材料特性以及利用用户子程序和其它应用软件进行数据交换等等。
这些用户子程序接口使用户解决一些问题时有很大的灵活性,同时大大的扩充了ABAQUS的功能。
例如:如果荷载条件是时间的函数,这在ABAQUS/CAE 和INPUT 文件中是难以实现的,但在用户子程序DLOAD中就很容易实现。
二.在ABAQUS中使用用户子程序ABAQUS的用户子程序是根据ABAQUS提供的相应接口,按照FORTRAN语法用户自己编写的代码。
在一个算例中,用户可以用到多个用户子程序,但必须把它们放在一个以.FOR为扩展名的文件中。
运行带有用户子程序的算例时有两种方法,一是在CAE中运行,在EDIT JOB菜单的GENERAL子菜单的USER SUBROUTINE FILE对话框中选择用户子程序所在的文件即可;另外是在ABABQUS COMMAND用运行,语法如下:ABAQUS JOB=[JOB] USER¡[.FOR]¡C用户在编写用户子程序时,要注意以下几点:1.用户子程序不能嵌套。
即任何用户子程序都不能调用任何其他用户子程Home序,但可以调用用户自己编写的FORTRAN子程序和ABAQUS应用程序。
abaqus用户子程序vusdfld单元 积分点出的极限应变
![abaqus用户子程序vusdfld单元 积分点出的极限应变](https://img.taocdn.com/s3/m/1fef255cc4da50e2524de518964bcf84b9d52daa.png)
abaqus用户子程序vusdfld单元积分点出的极限应变摘要:1.Abaqus 用户子程序简介2.VUSDFLD 单元概述3.积分点出的极限应变计算方法4.应用实例与结果分析5.总结与展望正文:【1.Abaqus 用户子程序简介】Abaqus 是一款广泛应用于结构分析和仿真的有限元软件。
用户子程序是Abaqus 中一种强大的功能,可以通过自定义编程实现对模型的特定处理和计算。
用户子程序可以在Abaqus 的各个分析步骤中调用,为分析过程提供更多的灵活性和可控性。
【2.VUSDFLD 单元概述】VUSDFLD 单元是Abaqus 中的一种特殊的单元类型,用于模拟材料在非常规应力状态下的行为,例如在高应变率和高应变下的材料行为。
VUSDFLD 单元可以定义材料的本构关系,即应力和应变之间的关系,以便在Abaqus 中进行分析。
【3.积分点出的极限应变计算方法】在Abaqus 中,可以通过用户子程序实现对模型中各积分点的极限应变计算。
具体步骤如下:(1)首先,在用户子程序中定义模型的几何参数、材料参数和边界条件。
(2)然后,通过循环遍历模型中的所有积分点,计算每个积分点的应力。
(3)接着,根据材料的本构关系,计算每个积分点的应变。
(4)最后,判断每个积分点的应变是否达到极限应变,如果达到则输出结果,并结束计算。
【4.应用实例与结果分析】假设有一个模型,由VUSDFLD 单元组成,需要计算在特定载荷条件下各积分点的极限应变。
可以通过编写用户子程序实现该功能。
在程序中,首先定义模型的几何参数、材料参数和边界条件,然后计算每个积分点的应力和应变,最后判断并输出极限应变结果。
【5.总结与展望】通过Abaqus 用户子程序和VUSDFLD 单元,可以实现对模型的极限应变计算。
该方法具有较强的灵活性和可控性,可以满足不同问题的分析需求。
abaqus用户子程序vusdfld单元 积分点出的极限应变
![abaqus用户子程序vusdfld单元 积分点出的极限应变](https://img.taocdn.com/s3/m/7113c8e4a48da0116c175f0e7cd184254b351b2b.png)
abaqus用户子程序vusdfld单元积分点出的极限应变(最新版)目录1.Abaqus 用户子程序概述2.VUSDFLD 单元介绍3.积分点出的极限应变计算方法4.应用实例与注意事项正文【1.Abaqus 用户子程序概述】Abaqus 是一款广泛应用于工程领域的有限元分析软件,其强大的功能和灵活的接口使得用户可以根据实际需求进行自定义和扩展。
用户子程序是 Abaqus 提供的一种编程接口,可以通过编写 C 或 Fortran 语言的程序,实现对 Abaqus 求解过程的控制和结果的后处理。
【2.VUSDFLD 单元介绍】VUSDFLD 单元是 Abaqus 中的一种通用粘结滑移软弹簧非线性几何单元,适用于模拟材料在复杂应力/应变路径下的非线性行为。
该单元可以描述材料在拉伸、压缩、剪切和扭转过程中的非线性特性,广泛应用于橡胶、泡沫、复合材料等非线性材料的分析中。
【3.积分点出的极限应变计算方法】在 Abaqus 中,通过用户子程序可以实现对 VUSDFLD 单元积分点出的极限应变的计算。
具体步骤如下:(1)首先,在用户子程序中定义一个函数,用于计算 VUSDFLD 单元的积分点出的极限应变。
该函数需要接收单元的各阶导数作为输入参数,并返回积分点出的极限应变。
(2)然后,在 Abaqus 的脚本中,通过调用该函数,获取 VUSDFLD 单元的积分点出的极限应变。
【4.应用实例与注意事项】在实际应用中,用户可以根据具体的材料特性和加载条件,编写相应的用户子程序,实现对 VUSDFLD 单元积分点出的极限应变的计算。
在编写过程中,需要注意以下几点:(1)正确定义函数输入输出参数,与 Abaqus 的接口保持一致。
(2)确保函数的计算过程正确无误,避免出现数值不稳定或误差较大的情况。
(3)在实际应用中,可能需要根据不同的材料特性和加载条件,调整用户子程序的计算参数,以获得更准确的结果。
Abaqus材料用户子程序UMAT基础知识及手册例子完整解释
![Abaqus材料用户子程序UMAT基础知识及手册例子完整解释](https://img.taocdn.com/s3/m/e2cd26cf89eb172ded63b7f5.png)
1、为何需要使用用户材料子程序(User-Defined Material, UMAT )?很简单,当ABAQUS 没有提供我们需要的材料模型时。
所以,在决定自己定义一种新的材料模型之前,最好对ABAQUS 已经提供的模型心中有数,并且尽量使用现有的模型,因为这些模型已经经过详细的验证,并被广泛接受。
UMAT 子程序具有强大的功能,使用UMAT 子程序:(1)可以定义材料的本构关系,使用ABAQUS 材料库中没有包含的材料进行计算,扩充程序功能。
(2) 几乎可以用于力学行为分析的任何分析过程,几乎可以把用户材料属性赋予ABAQU S 中的任何单元。
(3) 必须在UMAT 中提供材料本构模型的雅可比(Jacobian )矩阵,即应力增量对应变增量的变化率。
(4) 可以和用户子程序“USDFLD ”联合使用,通过“USDFLD ”重新定义单元每一物质点上传递到UMAT 中场变量的数值。
2、需要哪些基础知识?先看一下ABAQUS 手册(ABAQUS Analysis User's Manual )里的一段话:Warning: The use of this option generally requires considerable expertise(一定的专业知识). The user is cautioned that the implementation (实现) of any realistic constitutive (基本) model requires extensive (广泛的) development and testing. Initial testing on a single eleme nt model with prescribed traction loading (指定拉伸载荷) is strongly recommended. 但这并不意味着非力学专业,或者力学基础知识不很丰富者就只能望洋兴叹,因为我们的任务不是开发一套完整的有限元软件,而只是提供一个描述材料力学性能的本构方程(Constitutive equation )而已。
ABAQUS子程序
![ABAQUS子程序](https://img.taocdn.com/s3/m/bcdb1a3710661ed9ad51f3f2.png)
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' CDIMENSION TIME(2),COORDS(3)CHARACTER*80 SNAMEPARAMETER (PLOAD=100.E4)IF (KSTEP.EQ.1) THEN !当STEP=1时的荷载大小P=PLOADELSE IF (KSTEP.EQ.2) THEN !当STEP=2时的荷载大小P=COORDS(1)*PLOAD !施加在积分点的荷载P是坐标的函数ELSE IF (KSTEP.EQ.3) THEN !当STEP=3时的荷载大小P=COORDS(1)**2*PLOADELSE IF (KSTEP.EQ.4) THEN !当STEP=4时的荷载大小P=COORDS(1)**3*PLOADELSE IF (KSTEP.EQ.5) THEN !当STEP=5时的荷载大小P=COORDS(1)**4*PLOADEND IFRETURNENDUMAT 子程序具有强大的功能,使用UMAT 子程序:(1) 可以定义材料的本构关系,使用ABAQUS 材料库中没有包含的材料进行计算,扩充程序功能。
abaqus2用户单元子程序(1)
![abaqus2用户单元子程序(1)](https://img.taocdn.com/s3/m/d4fb965f700abb68a882fb98.png)
20 ABAQUS用户单元子程序(UEL)在这一章中将列举两个在这些年里发展过的ABAQUS/Standard用户单元子程序(UEL)。
第一个例子是一个非线性的索单元,我们的目的是通过这个比较简单的例子让读者了解用户单元子程序的基本开发过程;第二个例子是一个用于计算应变梯度理论的单元,应变梯度是当今比较热点的一个科研前沿问题,有各种理论,我们为了验证新的理论,需要数值结果与实验对照来进行评价,整个例子的目的是通过它说明用户子单元可以求解的问题范围很广,但是由于内容比较艰深,程序也很长,所以这个例子我们并没有给出最后的全部程序。
另外,到目前为止,ABAQUS还只有隐式求解器ABAQUS/Standard支持用户自定义单元,而显式求解器ABAQUS/Explicit中还不支持这一功能。
非线性索单元20.1.1 背景钢索斜拉桥和斜拉索结构广泛应用于土木工程建筑上。
索力的计算分析是设计和施工的关键环节。
清华大学工程力学系在采用ABAQUS进行荆沙长江斜拉桥的计算机仿真分析(这个项目我们已在第15章“ABAQUS在土木工程中的应用(一)——荆州长江大桥南汊斜拉桥结构三维仿真分析”中讨论过)时,也曾进行了自行建立索单元的尝试。
本节介绍的就是这方面的工作。
香港理工大学土木与结构工程系采用ABAQUS有限元软件进行计算,完成了香港Ting Kau斜拉桥和Tsing Ma悬索桥的结构计算和分析。
对于钢索计算,他们采用梁单元进行模拟。
由于梁单元含有弯曲刚度,计算的高阶频率值偏高,周期较低。
一般假设索是单向受拉力的构件。
随着应变的非线性增加,索力呈非线性增加。
尽管ABAQUS单元库中有500个以上的单元类型,但是,还没有索单元。
本文发展了三维非线性索单元模型,形成ABAQUS的用户单元子程序,可以利用ABAQUS输入文件调入到具体的分析中。
通过静态和动态例题的计算比较,索单元工作良好。
20.1.2 基本公式在三维索单元计算中,如图20-1所示,坐标x 和位移u 的变量表达式为:ij ji i j ji u u u x x x -=-= (x,y,z) (u,v,w) (20-1)应变的公式为:()⎥⎦⎤⎢⎣⎡+++++=222211ji ji ji ji ji ji ji ji ji w v u w z v y u x L ε (20-2)公式(20-2)中,L 为索的长度,索的张力为:0N AE N +=ε(20-3)在公式(20-3)中,A 为截面面积,E 为弹性模量,N 0为初始张力。
Abaqus用户子程序USDFLD实例详解-kxh
![Abaqus用户子程序USDFLD实例详解-kxh](https://img.taocdn.com/s3/m/095a5541e45c3b3567ec8b44.png)
NUAA
NUAA--Kong Xianghong
Page � 4
2. Solution
2.3. Create MeshPart
NUAA
NUAA--Kong Xianghong
Page � 5
2. Solution
2.4. Create Solid Layers
NUAA
2.15. Create Datum & Assign Material Orientation
NUAA
创建Datum csys-1,原点为(0,0,0),X轴上的点为(1,0,0),X-Y平面上的点为(0,1,0)。
NUAA--Kong Xianghong
Page � 19
2. Solution
2.8. Create Step
NUAA
NUAA--Kong Xianghong
Page � 10
2. Solution
2.9. Create & Edit Field Output & History Output
NUAA
NUAA--Kong ng
Page � 11
2. Solution
NUAA--Kong Xianghong
Page � 22
2. Solution
2.17. Results & Visualization (3)
NUAA
NUAA--Kong Xianghong
Page � 23
2. Solution
2.17. Results & Visualization (4)
NUAA--Kong Xianghong
ABAQUS_Standard用户材料子程序实例-Johnson-Cook金属本构模型卢剑锋庄茁
![ABAQUS_Standard用户材料子程序实例-Johnson-Cook金属本构模型卢剑锋庄茁](https://img.taocdn.com/s3/m/71a57d74ddccda38376baf9f.png)
。动运体刚的料材虑考来转旋过通须必都量张者或量矢何任的中阵矩量变态 状�外以新更据数的起引为行构本料材了除�题问变应限有于对�是的同不阵矩量张力应和 结的步量增在。中 TAMU 到递传据数的后新更将时始开步量增后然�据数新更先中 NAPXEU 或 。据数的中阵矩量变态状新更须必束
DLFDSU 序程子在可也 。中 TAMU 到递传值数将时始开步量增在 �阵矩的量变态状储存于用
3m/gK[
]
质 性
数参型模 kooC-nosnhoJ
度密
义定料材的件试 3 表 。序程合拟值数的面前于源来值数的数参中
型模 kooC-nosnhoJ 中其�示所 3 表如数参料材�料材义定自中 TAMU 在户用用采件试。3m/gk
301×58.7
为度密�3.0 和 aPG002 为别分比松泊和量模性弹�料材性弹线用使杆力出和杆力入
) V E T A T S N ( V ET A T S
的量张力应中 TAMU。分部旋共的量张力应理处需只中 TAMU 在此因�动转体刚了行进的量张 力应对经已就前之 TAMU 用调步量增个一�题问变应限有的动转体刚含包于对。新更阵矩量张 阵矩量张力应�始开的步量增在。量分切剪个 RHSN 和量分接直个 IDN 应对�阵矩量张力应 力应对将 TAMU 束结的步量增在�中 TAMU 到递传口接的间之序程主和 TAMU 过通值数的中
�候时的 TAMU 用使是但。义定项选”CITSALE*“过通中义定料材在常通�值的量模切剪始初料 材于基度刚些这。的义定质性性弹的中性属料材过通是度刚切剪向横的元单梁、壳、板和度刚制 控漏沙的元单分积缩减常通。度刚切剪向横和度刚制控漏沙的元单意注要需时 TAMU 用使 。 RHSN + IDN = SNETN �数个的量分力应总� SNETN �数个的量分力应切剪� RHSN �量增间时的步量增� EMITD �阵矩变应� ) SNETN ( NARTS �量变他其 。出输量能为作
ABAQUS用户子程序小例子
![ABAQUS用户子程序小例子](https://img.taocdn.com/s3/m/d66c70b710661ed9ad51f3c1.png)
实例:最简单的杆件单轴拉伸,材料本构模型采用线弹性模型,E=206000MPa, ν=0.3。
左端约束,右端施加均布力q=1。
Author: xueweek希望这个小例子能帮助初学者减少一两天的摸索时间1材料本构模型:σ=Eε在Property中定义材料时,在General下选中User Material,输入206000,0.3。
以上两个数值代表E、ν,在用户子程序中代表着PROPS(1)、PROPS(2)。
然后在General 下选中Depvar,由于该例子中没有使用状态变量,因此在第一项中输入大于0的数值即可。
建模大家都会,故省略2ABAQUS调用用户子程序大家在论坛上也可以看到很多方法,我这里指介绍最简单的一种。
在Job Manager中点击Edit选项,在General选项的最后一项中选择自己建立好的用户子程序文件。
(注:用户子程序文件可以使用文本编辑器进行编辑,当然也可以用Fortran 编译器,如果对自己的用户子程序文件的语法不放心,可以先用Fortan编译器进行编译,不过编译前要先建立Project,关于Fortran编译,这里不再介绍)。
完成后,点击submit即可进行分析。
3结果以下两张图分布是用户子程序和ABAQUS自带的材料模型得到的应力云图,可以看出两种图形完全相同。
4UMAT子程序对于初学者来说,需要注意的是,FORTRAN对于程序语言格式上的要求。
例如,对于FORTRAN语言,前六个字符必须空出来,等等。
检查语法最好的方法就是在FORTRAN编译器上进行编译。
SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,DDSDDT,1 DRPLDE,DRPLDT,STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,2 CMNAME,NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,3 PNEWDT,CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)include 'aba_param.inc'CHARACTER*8 CMNAMEDIMENSION STRESS(NTENS),STATEV(NSTATV),DDSDDE(NTENS,NTENS),1 DDSDDT(NTENS),DRPLDE(NTENS),STRAN(NTENS),DSTRAN(NTENS),2 TIME(2),PREDEF(1),DPRED(1),PROPS(NPROPS),COORDS(3),DROT(3,3),3 DFGRD0(3,3),DFGRD1(3,3)C UMAT FOR ISOTROPIC ELASTICITYC CANNOT BE USED FOR PLANE STRESSC ----------------------------------------------------------------C PROPS(1) - EC PROPS(2) - NUC ----------------------------------------------------------------IF (NDI.NE.3) THENWRITE (6,*) 'THIS UMAT MAY ONLY BE USED FOR ELEMENTS 1 WITH THREE DIRECT STRESS COMPONENTS'CALL XITENDIFCC ELASTIC PROPERTIESEMOD=PROPS(1)ENU=PROPS(2)EBULK3=EMOD/(1-2*ENU)EG2=EMOD/(1+ENU)EG=EG2/2EG3=3*EGELAM=(EBULK3-EG2)/3CC ELASTIC STIFFNESSCDO K1=1, NDIDO K2=1, NDIDDSDDE(K2, K1)=ELAMEND DODDSDDE(K1, K1)=EG2+ELAMEND DODO K1=NDI+1, NTENSDDSDDE(K1 ,K1)=EGEND DOCC CALCULATE STRESSCDO K1=1, NTENSDO K2=1, NTENSSTRESS(K2)=STRESS(K2)+DDSDDE(K2, K1)*DSTRAN(K1)END DOEND DOCRETURNEND5INP输入文件*Heading** Job name: Job-1 Model name: Model-1** Generated by: Abaqus/CAE 6.10-1*Preprint, echo=NO, model=NO, history=NO, contact=NO**** PARTS***Part, name=Part-1*End Part****** ASSEMBLY***Assembly, name=Assembly***Instance, name=Part-1-1, part=Part-1*Node1, 0.0399999991, 0.0399999991, 0.400000006 2, 0.0399999991, 0.029*******, 0.400000006 3, 0.0399999991, 0.0199999996, 0.400000006 4, 0.0399999991, 0.00999999978, 0.400000006 5, 0.0399999991, 0., 0.400000006 6, 0.0399999991, 0.0399999991, 0.389999986 7, 0.0399999991, 0.029*******, 0.389999986 8, 0.0399999991, 0.0199999996, 0.389999986 9, 0.0399999991, 0.00999999978, 0.389999986 10, 0.0399999991, 0., 0.389999986 11, 0.0399999991, 0.0399999991, 0.379999995 12, 0.0399999991, 0.029*******, 0.379999995 13, 0.0399999991, 0.0199999996, 0.379999995 14, 0.0399999991, 0.00999999978, 0.379999995 15, 0.0399999991, 0., 0.379999995 16, 0.0399999991, 0.0399999991, 0.370000005 17, 0.0399999991, 0.029*******, 0.370000005 18, 0.0399999991, 0.0199999996, 0.370000005 19, 0.0399999991, 0.00999999978, 0.370000005 20, 0.0399999991, 0., 0.370000005 21, 0.0399999991, 0.0399999991, 0.360000014 22, 0.0399999991, 0.029*******, 0.360000014 23, 0.0399999991, 0.0199999996, 0.360000014 24, 0.0399999991, 0.00999999978, 0.360000014 25, 0.0399999991, 0., 0.360000014 26, 0.0399999991, 0.0399999991, 0.349999994 27, 0.0399999991, 0.029*******, 0.349999994 28, 0.0399999991, 0.0199999996, 0.349999994 29, 0.0399999991, 0.00999999978, 0.349999994 30, 0.0399999991, 0., 0.349999994 31, 0.0399999991, 0.0399999991, 0.340000004 32, 0.0399999991, 0.029*******, 0.340000004 33, 0.0399999991, 0.0199999996, 0.340000004 34, 0.0399999991, 0.00999999978, 0.340000004 35, 0.0399999991, 0., 0.340000004 36, 0.0399999991, 0.0399999991, 0.330000013 37, 0.0399999991, 0.029*******, 0.330000013 38, 0.0399999991, 0.0199999996, 0.330000013 39, 0.0399999991, 0.00999999978, 0.330000013 40, 0.0399999991, 0., 0.330000013 41, 0.0399999991, 0.0399999991, 0.319999993 42, 0.0399999991, 0.029*******, 0.319999993 43, 0.0399999991, 0.0199999996, 0.319999993 44, 0.0399999991, 0.00999999978, 0.319999993 45, 0.0399999991, 0., 0.319999993 46, 0.0399999991, 0.0399999991, 0.310000002 47, 0.0399999991, 0.029*******, 0.310000002 48, 0.0399999991, 0.0199999996, 0.310000002 49, 0.0399999991, 0.00999999978, 0.310000002 50, 0.0399999991, 0., 0.310000002 51, 0.0399999991, 0.0399999991, 0.300000012 52, 0.0399999991, 0.029*******, 0.300000012 53, 0.0399999991, 0.0199999996, 0.300000012 54, 0.0399999991, 0.00999999978, 0.300000012 55, 0.0399999991, 0., 0.300000012 56, 0.0399999991, 0.0399999991, 0.289999992 57, 0.0399999991, 0.029*******, 0.289999992 58, 0.0399999991, 0.0199999996, 0.289999992 59, 0.0399999991, 0.00999999978, 0.289999992 60, 0.0399999991, 0., 0.289999992 61, 0.0399999991, 0.0399999991, 0.280000001 62, 0.0399999991, 0.029*******, 0.280000001 63, 0.0399999991, 0.0199999996, 0.280000001 64, 0.0399999991, 0.00999999978, 0.280000001 65, 0.0399999991, 0., 0.280000001 66, 0.0399999991, 0.0399999991, 0.270000011 67, 0.0399999991, 0.029*******, 0.270000011 68, 0.0399999991, 0.0199999996, 0.270000011 69, 0.0399999991, 0.00999999978, 0.270000011 70, 0.0399999991, 0., 0.270000011 71, 0.0399999991, 0.0399999991, 0.25999999 72, 0.0399999991, 0.029*******, 0.25999999 73, 0.0399999991, 0.0199999996, 0.2599999976, 0.0399999991, 0.0399999991, 0.25 77, 0.0399999991, 0.029*******, 0.25 78, 0.0399999991, 0.0199999996, 0.25 79, 0.0399999991, 0.00999999978, 0.25 80, 0.0399999991, 0., 0.25 81, 0.0399999991, 0.0399999991, 0.239999995 82, 0.0399999991, 0.029*******, 0.239999995 83, 0.0399999991, 0.0199999996, 0.239999995 84, 0.0399999991, 0.00999999978, 0.239999995 85, 0.0399999991, 0., 0.239999995 86, 0.0399999991, 0.0399999991, 0.230000004 87, 0.0399999991, 0.029*******, 0.230000004 88, 0.0399999991, 0.0199999996, 0.230000004 89, 0.0399999991, 0.00999999978, 0.230000004 90, 0.0399999991, 0., 0.230000004 91, 0.0399999991, 0.0399999991, 0.219999999 92, 0.0399999991, 0.029*******, 0.219999999 93, 0.0399999991, 0.0199999996, 0.219999999 94, 0.0399999991, 0.00999999978, 0.219999999 95, 0.0399999991, 0., 0.219999999 96, 0.0399999991, 0.0399999991, 0.209999993 97, 0.0399999991, 0.029*******, 0.209999993 98, 0.0399999991, 0.0199999996, 0.209999993 99, 0.0399999991, 0.00999999978, 0.209999993 100, 0.0399999991, 0., 0.209999993 101, 0.0399999991, 0.0399999991, 0.200000003 102, 0.0399999991, 0.029*******, 0.200000003 103, 0.0399999991, 0.0199999996, 0.200000003 104, 0.0399999991, 0.00999999978, 0.200000003 105, 0.0399999991, 0., 0.200000003 106, 0.0399999991, 0.0399999991, 0.189999998 107, 0.0399999991, 0.029*******, 0.189999998 108, 0.0399999991, 0.0199999996, 0.189999998 109, 0.0399999991, 0.00999999978, 0.189999998 110, 0.0399999991, 0., 0.189999998 111, 0.0399999991, 0.0399999991, 0.180000007 112, 0.0399999991, 0.029*******, 0.180000007 113, 0.0399999991, 0.0199999996, 0.180000007 114, 0.0399999991, 0.00999999978, 0.180000007 115, 0.0399999991, 0., 0.180000007 116, 0.0399999991, 0.0399999991, 0.170000002 117, 0.0399999991, 0.029*******, 0.170000002 118, 0.0399999991, 0.0199999996, 0.170000002 119, 0.0399999991, 0.00999999978, 0.170000002 120, 0.0399999991, 0., 0.170000002 121, 0.0399999991, 0.0399999991, 0.159999996 122, 0.0399999991, 0.029*******, 0.159999996 123, 0.0399999991, 0.0199999996, 0.159999996 124, 0.0399999991, 0.00999999978, 0.159999996 125, 0.0399999991, 0., 0.159999996 126, 0.0399999991, 0.0399999991, 0.150000006 127, 0.0399999991, 0.029*******, 0.150000006 128, 0.0399999991, 0.0199999996, 0.150000006 129, 0.0399999991, 0.00999999978, 0.150000006 130, 0.0399999991, 0., 0.150000006 131, 0.0399999991, 0.0399999991, 0.140000001 132, 0.0399999991, 0.029*******, 0.140000001 133, 0.0399999991, 0.0199999996, 0.140000001 134, 0.0399999991, 0.00999999978, 0.140000001 135, 0.0399999991, 0., 0.140000001 136, 0.0399999991, 0.0399999991, 0.129999995 137, 0.0399999991, 0.029*******, 0.129999995 138, 0.0399999991, 0.0199999996, 0.129999995 139, 0.0399999991, 0.00999999978, 0.129999995 140, 0.0399999991, 0., 0.129999995 141, 0.0399999991, 0.0399999991, 0.119999997 142, 0.0399999991, 0.029*******, 0.119999997 143, 0.0399999991, 0.0199999996, 0.119999997 144, 0.0399999991, 0.00999999978, 0.119999997 145, 0.0399999991, 0., 0.119999997 146, 0.0399999991, 0.0399999991, 0.109999999 147, 0.0399999991, 0.029*******, 0.109999999 148, 0.0399999991, 0.0199999996, 0.109999999 149, 0.0399999991, 0.00999999978, 0.109999999 150, 0.0399999991, 0., 0.109999999 151, 0.0399999991, 0.0399999991, 0.100000001 152, 0.0399999991, 0.029*******, 0.100000001 153, 0.0399999991, 0.0199999996, 0.100000001156, 0.0399999991, 0.0399999991, 0.0900000036 157, 0.0399999991, 0.029*******, 0.0900000036 158, 0.0399999991, 0.0199999996, 0.0900000036 159, 0.0399999991, 0.00999999978, 0.0900000036 160, 0.0399999991, 0., 0.0900000036 161, 0.0399999991, 0.0399999991, 0.0799999982 162, 0.0399999991, 0.029*******, 0.0799999982 163, 0.0399999991, 0.0199999996, 0.0799999982 164, 0.0399999991, 0.00999999978, 0.0799999982 165, 0.0399999991, 0., 0.0799999982 166, 0.0399999991, 0.0399999991, 0.0700000003 167, 0.0399999991, 0.029*******, 0.0700000003 168, 0.0399999991, 0.0199999996, 0.0700000003 169, 0.0399999991, 0.00999999978, 0.0700000003 170, 0.0399999991, 0., 0.0700000003 171, 0.0399999991, 0.0399999991, 0.0599999987 172, 0.0399999991, 0.029*******, 0.0599999987 173, 0.0399999991, 0.0199999996, 0.0599999987 174, 0.0399999991, 0.00999999978, 0.0599999987 175, 0.0399999991, 0., 0.0599999987 176, 0.0399999991, 0.0399999991, 0.0500000007 177, 0.0399999991, 0.029*******, 0.0500000007 178, 0.0399999991, 0.0199999996, 0.0500000007 179, 0.0399999991, 0.00999999978, 0.0500000007 180, 0.0399999991, 0., 0.0500000007 181, 0.0399999991, 0.0399999991, 0.0399999991 182, 0.0399999991, 0.029*******, 0.0399999991 183, 0.0399999991, 0.0199999996, 0.0399999991 184, 0.0399999991, 0.00999999978, 0.0399999991 185, 0.0399999991, 0., 0.0399999991 186, 0.0399999991, 0.0399999991, 0.029******* 187, 0.0399999991, 0.029*******, 0.029******* 188, 0.0399999991, 0.0199999996, 0.029******* 189, 0.0399999991, 0.00999999978, 0.029******* 190, 0.0399999991, 0., 0.029******* 191, 0.0399999991, 0.0399999991, 0.0199999996 192, 0.0399999991, 0.029*******, 0.0199999996 193, 0.0399999991, 0.0199999996, 0.0199999996 194, 0.0399999991, 0.00999999978, 0.0199999996 195, 0.0399999991, 0., 0.0199999996 196, 0.0399999991, 0.0399999991, 0.00999999978 197, 0.0399999991, 0.029*******, 0.00999999978 198, 0.0399999991, 0.0199999996, 0.00999999978 199, 0.0399999991, 0.00999999978, 0.00999999978 200, 0.0399999991, 0., 0.00999999978 201, 0.0399999991, 0.0399999991, 0. 202, 0.0399999991, 0.029*******, 0. 203, 0.0399999991, 0.0199999996, 0. 204, 0.0399999991, 0.00999999978, 0. 205, 0.0399999991, 0., 0. 206, 0.029*******, 0.0399999991, 0.400000006 207, 0.029*******, 0.029*******, 0.400000006 208, 0.029*******, 0.0199999996, 0.400000006 209, 0.029*******, 0.00999999978, 0.400000006 210, 0.029*******, 0., 0.400000006 211, 0.029*******, 0.0399999991, 0.389999986 212, 0.029*******, 0.029*******, 0.389999986 213, 0.029*******, 0.0199999996, 0.389999986 214, 0.029*******, 0.00999999978, 0.389999986 215, 0.029*******, 0., 0.389999986 216, 0.029*******, 0.0399999991, 0.379999995 217, 0.029*******, 0.029*******, 0.379999995 218, 0.029*******, 0.0199999996, 0.379999995 219, 0.029*******, 0.00999999978, 0.379999995 220, 0.029*******, 0., 0.379999995 221, 0.029*******, 0.0399999991, 0.370000005 222, 0.029*******, 0.029*******, 0.370000005 223, 0.029*******, 0.0199999996, 0.370000005 224, 0.029*******, 0.00999999978, 0.370000005 225, 0.029*******, 0., 0.370000005 226, 0.029*******, 0.0399999991, 0.360000014 227, 0.029*******, 0.029*******, 0.360000014 228, 0.029*******, 0.0199999996, 0.360000014 229, 0.029*******, 0.00999999978, 0.360000014 230, 0.029*******, 0., 0.360000014 231, 0.029*******, 0.0399999991, 0.349999994 232, 0.029*******, 0.029*******, 0.349999994 233, 0.029*******, 0.0199999996, 0.349999994236, 0.029*******, 0.0399999991, 0.340000004 237, 0.029*******, 0.029*******, 0.340000004 238, 0.029*******, 0.0199999996, 0.340000004 239, 0.029*******, 0.00999999978, 0.340000004 240, 0.029*******, 0., 0.340000004 241, 0.029*******, 0.0399999991, 0.330000013 242, 0.029*******, 0.029*******, 0.330000013 243, 0.029*******, 0.0199999996, 0.330000013 244, 0.029*******, 0.00999999978, 0.330000013 245, 0.029*******, 0., 0.330000013 246, 0.029*******, 0.0399999991, 0.319999993 247, 0.029*******, 0.029*******, 0.319999993 248, 0.029*******, 0.0199999996, 0.319999993 249, 0.029*******, 0.00999999978, 0.319999993 250, 0.029*******, 0., 0.319999993 251, 0.029*******, 0.0399999991, 0.310000002 252, 0.029*******, 0.029*******, 0.310000002 253, 0.029*******, 0.0199999996, 0.310000002 254, 0.029*******, 0.00999999978, 0.310000002 255, 0.029*******, 0., 0.310000002 256, 0.029*******, 0.0399999991, 0.300000012 257, 0.029*******, 0.029*******, 0.300000012 258, 0.029*******, 0.0199999996, 0.300000012 259, 0.029*******, 0.00999999978, 0.300000012 260, 0.029*******, 0., 0.300000012 261, 0.029*******, 0.0399999991, 0.289999992 262, 0.029*******, 0.029*******, 0.289999992 263, 0.029*******, 0.0199999996, 0.289999992 264, 0.029*******, 0.00999999978, 0.289999992 265, 0.029*******, 0., 0.289999992 266, 0.029*******, 0.0399999991, 0.280000001 267, 0.029*******, 0.029*******, 0.280000001 268, 0.029*******, 0.0199999996, 0.280000001 269, 0.029*******, 0.00999999978, 0.280000001 270, 0.029*******, 0., 0.280000001 271, 0.029*******, 0.0399999991, 0.270000011 272, 0.029*******, 0.029*******, 0.270000011 273, 0.029*******, 0.0199999996, 0.270000011 274, 0.029*******, 0.00999999978, 0.270000011 275, 0.029*******, 0., 0.270000011 276, 0.029*******, 0.0399999991, 0.25999999 277, 0.029*******, 0.029*******, 0.25999999 278, 0.029*******, 0.0199999996, 0.25999999 279, 0.029*******, 0.00999999978, 0.25999999 280, 0.029*******, 0., 0.25999999 281, 0.029*******, 0.0399999991, 0.25 282, 0.029*******, 0.029*******, 0.25 283, 0.029*******, 0.0199999996, 0.25 284, 0.029*******, 0.00999999978, 0.25 285, 0.029*******, 0., 0.25 286, 0.029*******, 0.0399999991, 0.239999995 287, 0.029*******, 0.029*******, 0.239999995 288, 0.029*******, 0.0199999996, 0.239999995 289, 0.029*******, 0.00999999978, 0.239999995 290, 0.029*******, 0., 0.239999995 291, 0.029*******, 0.0399999991, 0.230000004 292, 0.029*******, 0.029*******, 0.230000004 293, 0.029*******, 0.0199999996, 0.230000004 294, 0.029*******, 0.00999999978, 0.230000004 295, 0.029*******, 0., 0.230000004 296, 0.029*******, 0.0399999991, 0.219999999 297, 0.029*******, 0.029*******, 0.219999999 298, 0.029*******, 0.0199999996, 0.219999999 299, 0.029*******, 0.00999999978, 0.219999999 300, 0.029*******, 0., 0.219999999 301, 0.029*******, 0.0399999991, 0.209999993 302, 0.029*******, 0.029*******, 0.209999993 303, 0.029*******, 0.0199999996, 0.209999993 304, 0.029*******, 0.00999999978, 0.209999993 305, 0.029*******, 0., 0.209999993 306, 0.029*******, 0.0399999991, 0.200000003 307, 0.029*******, 0.029*******, 0.200000003 308, 0.029*******, 0.0199999996, 0.200000003 309, 0.029*******, 0.00999999978, 0.200000003 310, 0.029*******, 0., 0.200000003 311, 0.029*******, 0.0399999991, 0.189999998 312, 0.029*******, 0.029*******, 0.189999998 313, 0.029*******, 0.0199999996, 0.189999998316, 0.029*******, 0.0399999991, 0.180000007 317, 0.029*******, 0.029*******, 0.180000007 318, 0.029*******, 0.0199999996, 0.180000007 319, 0.029*******, 0.00999999978, 0.180000007 320, 0.029*******, 0., 0.180000007 321, 0.029*******, 0.0399999991, 0.170000002 322, 0.029*******, 0.029*******, 0.170000002 323, 0.029*******, 0.0199999996, 0.170000002 324, 0.029*******, 0.00999999978, 0.170000002 325, 0.029*******, 0., 0.170000002 326, 0.029*******, 0.0399999991, 0.159999996 327, 0.029*******, 0.029*******, 0.159999996 328, 0.029*******, 0.0199999996, 0.159999996 329, 0.029*******, 0.00999999978, 0.159999996 330, 0.029*******, 0., 0.159999996 331, 0.029*******, 0.0399999991, 0.150000006 332, 0.029*******, 0.029*******, 0.150000006 333, 0.029*******, 0.0199999996, 0.150000006 334, 0.029*******, 0.00999999978, 0.150000006 335, 0.029*******, 0., 0.150000006 336, 0.029*******, 0.0399999991, 0.140000001 337, 0.029*******, 0.029*******, 0.140000001 338, 0.029*******, 0.0199999996, 0.140000001 339, 0.029*******, 0.00999999978, 0.140000001 340, 0.029*******, 0., 0.140000001 341, 0.029*******, 0.0399999991, 0.129999995 342, 0.029*******, 0.029*******, 0.129999995 343, 0.029*******, 0.0199999996, 0.129999995 344, 0.029*******, 0.00999999978, 0.129999995 345, 0.029*******, 0., 0.129999995 346, 0.029*******, 0.0399999991, 0.119999997 347, 0.029*******, 0.029*******, 0.119999997 348, 0.029*******, 0.0199999996, 0.119999997 349, 0.029*******, 0.00999999978, 0.119999997 350, 0.029*******, 0., 0.119999997 351, 0.029*******, 0.0399999991, 0.109999999 352, 0.029*******, 0.029*******, 0.109999999 353, 0.029*******, 0.0199999996, 0.109999999 354, 0.029*******, 0.00999999978, 0.109999999 355, 0.029*******, 0., 0.109999999 356, 0.029*******, 0.0399999991, 0.100000001 357, 0.029*******, 0.029*******, 0.100000001 358, 0.029*******, 0.0199999996, 0.100000001 359, 0.029*******, 0.00999999978, 0.100000001 360, 0.029*******, 0., 0.100000001 361, 0.029*******, 0.0399999991, 0.0900000036 362, 0.029*******, 0.029*******, 0.0900000036 363, 0.029*******, 0.0199999996, 0.0900000036 364, 0.029*******, 0.00999999978, 0.0900000036 365, 0.029*******, 0., 0.0900000036 366, 0.029*******, 0.0399999991, 0.0799999982 367, 0.029*******, 0.029*******, 0.0799999982 368, 0.029*******, 0.0199999996, 0.0799999982 369, 0.029*******, 0.00999999978, 0.0799999982 370, 0.029*******, 0., 0.0799999982 371, 0.029*******, 0.0399999991, 0.0700000003 372, 0.029*******, 0.029*******, 0.0700000003 373, 0.029*******, 0.0199999996, 0.0700000003 374, 0.029*******, 0.00999999978, 0.0700000003 375, 0.029*******, 0., 0.0700000003 376, 0.029*******, 0.0399999991, 0.0599999987 377, 0.029*******, 0.029*******, 0.0599999987 378, 0.029*******, 0.0199999996, 0.0599999987 379, 0.029*******, 0.00999999978, 0.0599999987 380, 0.029*******, 0., 0.0599999987 381, 0.029*******, 0.0399999991, 0.0500000007 382, 0.029*******, 0.029*******, 0.0500000007 383, 0.029*******, 0.0199999996, 0.0500000007 384, 0.029*******, 0.00999999978, 0.0500000007 385, 0.029*******, 0., 0.0500000007 386, 0.029*******, 0.0399999991, 0.0399999991 387, 0.029*******, 0.029*******, 0.0399999991 388, 0.029*******, 0.0199999996, 0.0399999991 389, 0.029*******, 0.00999999978, 0.0399999991 390, 0.029*******, 0., 0.0399999991 391, 0.029*******, 0.0399999991, 0.029******* 392, 0.029*******, 0.029*******, 0.029******* 393, 0.029*******, 0.0199999996, 0.029*******396, 0.029*******, 0.0399999991, 0.0199999996 397, 0.029*******, 0.029*******, 0.0199999996 398, 0.029*******, 0.0199999996, 0.0199999996 399, 0.029*******, 0.00999999978, 0.0199999996 400, 0.029*******, 0., 0.0199999996 401, 0.029*******, 0.0399999991, 0.00999999978 402, 0.029*******, 0.029*******, 0.00999999978 403, 0.029*******, 0.0199999996, 0.00999999978 404, 0.029*******, 0.00999999978, 0.00999999978 405, 0.029*******, 0., 0.00999999978 406, 0.029*******, 0.0399999991, 0. 407, 0.029*******, 0.029*******, 0. 408, 0.029*******, 0.0199999996, 0. 409, 0.029*******, 0.00999999978, 0. 410, 0.029*******, 0., 0. 411, 0.0199999996, 0.0399999991, 0.400000006 412, 0.0199999996, 0.029*******, 0.400000006 413, 0.0199999996, 0.0199999996, 0.400000006 414, 0.0199999996, 0.00999999978, 0.400000006 415, 0.0199999996, 0., 0.400000006 416, 0.0199999996, 0.0399999991, 0.389999986 417, 0.0199999996, 0.029*******, 0.389999986 418, 0.0199999996, 0.0199999996, 0.389999986 419, 0.0199999996, 0.00999999978, 0.389999986 420, 0.0199999996, 0., 0.389999986 421, 0.0199999996, 0.0399999991, 0.379999995 422, 0.0199999996, 0.029*******, 0.379999995 423, 0.0199999996, 0.0199999996, 0.379999995 424, 0.0199999996, 0.00999999978, 0.379999995 425, 0.0199999996, 0., 0.379999995 426, 0.0199999996, 0.0399999991, 0.370000005 427, 0.0199999996, 0.029*******, 0.370000005 428, 0.0199999996, 0.0199999996, 0.370000005 429, 0.0199999996, 0.00999999978, 0.370000005 430, 0.0199999996, 0., 0.370000005 431, 0.0199999996, 0.0399999991, 0.360000014 432, 0.0199999996, 0.029*******, 0.360000014 433, 0.0199999996, 0.0199999996, 0.360000014 434, 0.0199999996, 0.00999999978, 0.360000014 435, 0.0199999996, 0., 0.360000014 436, 0.0199999996, 0.0399999991, 0.349999994 437, 0.0199999996, 0.029*******, 0.349999994 438, 0.0199999996, 0.0199999996, 0.349999994 439, 0.0199999996, 0.00999999978, 0.349999994 440, 0.0199999996, 0., 0.349999994 441, 0.0199999996, 0.0399999991, 0.340000004 442, 0.0199999996, 0.029*******, 0.340000004 443, 0.0199999996, 0.0199999996, 0.340000004 444, 0.0199999996, 0.00999999978, 0.340000004 445, 0.0199999996, 0., 0.340000004 446, 0.0199999996, 0.0399999991, 0.330000013 447, 0.0199999996, 0.029*******, 0.330000013 448, 0.0199999996, 0.0199999996, 0.330000013 449, 0.0199999996, 0.00999999978, 0.330000013 450, 0.0199999996, 0., 0.330000013 451, 0.0199999996, 0.0399999991, 0.319999993 452, 0.0199999996, 0.029*******, 0.319999993 453, 0.0199999996, 0.0199999996, 0.319999993 454, 0.0199999996, 0.00999999978, 0.319999993 455, 0.0199999996, 0., 0.319999993 456, 0.0199999996, 0.0399999991, 0.310000002 457, 0.0199999996, 0.029*******, 0.310000002 458, 0.0199999996, 0.0199999996, 0.310000002 459, 0.0199999996, 0.00999999978, 0.310000002 460, 0.0199999996, 0., 0.310000002 461, 0.0199999996, 0.0399999991, 0.300000012 462, 0.0199999996, 0.029*******, 0.300000012 463, 0.0199999996, 0.0199999996, 0.300000012 464, 0.0199999996, 0.00999999978, 0.300000012 465, 0.0199999996, 0., 0.300000012 466, 0.0199999996, 0.0399999991, 0.289999992 467, 0.0199999996, 0.029*******, 0.289999992 468, 0.0199999996, 0.0199999996, 0.289999992 469, 0.0199999996, 0.00999999978, 0.289999992 470, 0.0199999996, 0., 0.289999992 471, 0.0199999996, 0.0399999991, 0.280000001 472, 0.0199999996, 0.029*******, 0.280000001 473, 0.0199999996, 0.0199999996, 0.280000001476, 0.0199999996, 0.0399999991, 0.270000011 477, 0.0199999996, 0.029*******, 0.270000011 478, 0.0199999996, 0.0199999996, 0.270000011 479, 0.0199999996, 0.00999999978, 0.270000011 480, 0.0199999996, 0., 0.270000011 481, 0.0199999996, 0.0399999991, 0.25999999 482, 0.0199999996, 0.029*******, 0.25999999 483, 0.0199999996, 0.0199999996, 0.25999999 484, 0.0199999996, 0.00999999978, 0.25999999 485, 0.0199999996, 0., 0.25999999 486, 0.0199999996, 0.0399999991, 0.25 487, 0.0199999996, 0.029*******, 0.25 488, 0.0199999996, 0.0199999996, 0.25 489, 0.0199999996, 0.00999999978, 0.25 490, 0.0199999996, 0., 0.25 491, 0.0199999996, 0.0399999991, 0.239999995 492, 0.0199999996, 0.029*******, 0.239999995 493, 0.0199999996, 0.0199999996, 0.239999995 494, 0.0199999996, 0.00999999978, 0.239999995 495, 0.0199999996, 0., 0.239999995 496, 0.0199999996, 0.0399999991, 0.230000004 497, 0.0199999996, 0.029*******, 0.230000004 498, 0.0199999996, 0.0199999996, 0.230000004 499, 0.0199999996, 0.00999999978, 0.230000004 500, 0.0199999996, 0., 0.230000004 501, 0.0199999996, 0.0399999991, 0.219999999 502, 0.0199999996, 0.029*******, 0.219999999 503, 0.0199999996, 0.0199999996, 0.219999999 504, 0.0199999996, 0.00999999978, 0.219999999 505, 0.0199999996, 0., 0.219999999 506, 0.0199999996, 0.0399999991, 0.209999993 507, 0.0199999996, 0.029*******, 0.209999993 508, 0.0199999996, 0.0199999996, 0.209999993 509, 0.0199999996, 0.00999999978, 0.209999993 510, 0.0199999996, 0., 0.209999993 511, 0.0199999996, 0.0399999991, 0.200000003 512, 0.0199999996, 0.029*******, 0.200000003 513, 0.0199999996, 0.0199999996, 0.200000003 514, 0.0199999996, 0.00999999978, 0.200000003 515, 0.0199999996, 0., 0.200000003 516, 0.0199999996, 0.0399999991, 0.189999998 517, 0.0199999996, 0.029*******, 0.189999998 518, 0.0199999996, 0.0199999996, 0.189999998 519, 0.0199999996, 0.00999999978, 0.189999998 520, 0.0199999996, 0., 0.189999998 521, 0.0199999996, 0.0399999991, 0.180000007 522, 0.0199999996, 0.029*******, 0.180000007 523, 0.0199999996, 0.0199999996, 0.180000007 524, 0.0199999996, 0.00999999978, 0.180000007 525, 0.0199999996, 0., 0.180000007 526, 0.0199999996, 0.0399999991, 0.170000002 527, 0.0199999996, 0.029*******, 0.170000002 528, 0.0199999996, 0.0199999996, 0.170000002 529, 0.0199999996, 0.00999999978, 0.170000002 530, 0.0199999996, 0., 0.170000002 531, 0.0199999996, 0.0399999991, 0.159999996 532, 0.0199999996, 0.029*******, 0.159999996 533, 0.0199999996, 0.0199999996, 0.159999996 534, 0.0199999996, 0.00999999978, 0.159999996 535, 0.0199999996, 0., 0.159999996 536, 0.0199999996, 0.0399999991, 0.150000006 537, 0.0199999996, 0.029*******, 0.150000006 538, 0.0199999996, 0.0199999996, 0.150000006 539, 0.0199999996, 0.00999999978, 0.150000006 540, 0.0199999996, 0., 0.150000006 541, 0.0199999996, 0.0399999991, 0.140000001 542, 0.0199999996, 0.029*******, 0.140000001 543, 0.0199999996, 0.0199999996, 0.140000001 544, 0.0199999996, 0.00999999978, 0.140000001 545, 0.0199999996, 0., 0.140000001 546, 0.0199999996, 0.0399999991, 0.129999995 547, 0.0199999996, 0.029*******, 0.129999995 548, 0.0199999996, 0.0199999996, 0.129999995 549, 0.0199999996, 0.00999999978, 0.129999995 550, 0.0199999996, 0., 0.129999995 551, 0.0199999996, 0.0399999991, 0.119999997 552, 0.0199999996, 0.029*******, 0.119999997 553, 0.0199999996, 0.0199999996, 0.119999997556, 0.0199999996, 0.0399999991, 0.109999999 557, 0.0199999996, 0.029*******, 0.109999999 558, 0.0199999996, 0.0199999996, 0.109999999 559, 0.0199999996, 0.00999999978, 0.109999999 560, 0.0199999996, 0., 0.109999999 561, 0.0199999996, 0.0399999991, 0.100000001 562, 0.0199999996, 0.029*******, 0.100000001 563, 0.0199999996, 0.0199999996, 0.100000001 564, 0.0199999996, 0.00999999978, 0.100000001 565, 0.0199999996, 0., 0.100000001 566, 0.0199999996, 0.0399999991, 0.0900000036567, 0.0199999996, 0.029*******, 0.0900000036568, 0.0199999996, 0.0199999996, 0.0900000036569, 0.0199999996, 0.00999999978, 0.0900000036 570, 0.0199999996, 0., 0.0900000036 571, 0.0199999996, 0.0399999991, 0.0799999982572, 0.0199999996, 0.029*******, 0.0799999982573, 0.0199999996, 0.0199999996, 0.0799999982574, 0.0199999996, 0.00999999978, 0.0799999982 575, 0.0199999996, 0., 0.0799999982 576, 0.0199999996, 0.0399999991, 0.0700000003577, 0.0199999996, 0.029*******, 0.0700000003578, 0.0199999996, 0.0199999996, 0.0700000003579, 0.0199999996, 0.00999999978, 0.0700000003 580, 0.0199999996, 0., 0.0700000003 581, 0.0199999996, 0.0399999991, 0.0599999987582, 0.0199999996, 0.029*******, 0.0599999987583, 0.0199999996, 0.0199999996, 0.0599999987584, 0.0199999996, 0.00999999978, 0.0599999987 585, 0.0199999996, 0., 0.0599999987 586, 0.0199999996, 0.0399999991, 0.0500000007587, 0.0199999996, 0.029*******, 0.0500000007588, 0.0199999996, 0.0199999996, 0.0500000007589, 0.0199999996, 0.00999999978, 0.0500000007 590, 0.0199999996, 0., 0.0500000007 591, 0.0199999996, 0.0399999991, 0.0399999991592, 0.0199999996, 0.029*******, 0.0399999991593, 0.0199999996, 0.0199999996, 0.0399999991594, 0.0199999996, 0.00999999978, 0.0399999991 595, 0.0199999996, 0., 0.0399999991 596, 0.0199999996, 0.0399999991, 0.029*******597, 0.0199999996, 0.029*******, 0.029*******598, 0.0199999996, 0.0199999996, 0.029*******599, 0.0199999996, 0.00999999978, 0.029******* 600, 0.0199999996, 0., 0.029******* 601, 0.0199999996, 0.0399999991, 0.0199999996602, 0.0199999996, 0.029*******, 0.0199999996603, 0.0199999996, 0.0199999996, 0.0199999996604, 0.0199999996, 0.00999999978, 0.0199999996 605, 0.0199999996, 0., 0.0199999996 606, 0.0199999996, 0.0399999991, 0.00999999978 607, 0.0199999996, 0.029*******, 0.00999999978 608, 0.0199999996, 0.0199999996, 0.00999999978 609, 0.0199999996, 0.00999999978, 0.00999999978 610, 0.0199999996, 0., 0.00999999978 611, 0.0199999996, 0.0399999991, 0. 612, 0.0199999996, 0.029*******, 0. 613, 0.0199999996, 0.0199999996, 0. 614, 0.0199999996, 0.00999999978, 0. 615, 0.0199999996, 0., 0. 616, 0.00999999978, 0.0399999991, 0.400000006 617, 0.00999999978, 0.029*******, 0.400000006 618, 0.00999999978, 0.0199999996, 0.400000006 619, 0.00999999978, 0.00999999978, 0.400000006 620, 0.00999999978, 0., 0.400000006 621, 0.00999999978, 0.0399999991, 0.389999986 622, 0.00999999978, 0.029*******, 0.389999986 623, 0.00999999978, 0.0199999996, 0.389999986 624, 0.00999999978, 0.00999999978, 0.389999986 625, 0.00999999978, 0., 0.389999986 626, 0.00999999978, 0.0399999991, 0.379999995 627, 0.00999999978, 0.029*******, 0.379999995 628, 0.00999999978, 0.0199999996, 0.379999995 629, 0.00999999978, 0.00999999978, 0.379999995 630, 0.00999999978, 0., 0.379999995 631, 0.00999999978, 0.0399999991, 0.370000005 632, 0.00999999978, 0.029*******, 0.370000005 633, 0.00999999978, 0.0199999996, 0.370000005。
ABAQUS用户材料子程序开发及应用
![ABAQUS用户材料子程序开发及应用](https://img.taocdn.com/s3/m/6f6ca324974bcf84b9d528ea81c758f5f61f290c.png)
ABAQUS用户材料子程序开发及应用ABAQUS用户材料子程序开发及应用摘要:本文介绍了ABAQUS用户材料子程序的开发与应用。
首先,简要介绍了ABAQUS软件及其在工程领域的广泛应用。
然后,详细阐述了用户材料子程序的概念及作用,并介绍了子程序的开发流程和必要步骤。
接着,以一个具体的材料模型开发为例,详细介绍了子程序的实现方法和注意事项。
最后,以轴对称挤压模拟为例,展示了用户材料子程序在实际工程分析中的应用,并讨论了其优点和局限性。
关键词:ABAQUS;用户材料子程序;开发;应用一、引言ABAQUS是一款广泛应用于工程领域的有限元分析软件。
其强大的建模和分析能力使得工程师可以准确地模拟和分析各种结构和材料的行为。
然而,对于一些非标准材料或特殊材料,ABAQUS自带的材料模型可能无法满足工程师的需求。
此时,用户材料子程序的开发就显得尤为重要。
二、用户材料子程序的概念及作用用户材料子程序是指由ABAQUS用户自行编写的用于描述非标准材料行为的子程序。
它可以根据特定的材料性质和应变-应力关系,定义材料模型的行为,并将其与ABAQUS的有限元分析过程相结合。
通过用户材料子程序,工程师可以更加准确地模拟和分析特殊材料的行为,提高分析结果的可靠性和准确性。
三、用户材料子程序的开发流程和步骤用户材料子程序的开发包括以下几个基本步骤:1. 确定材料模型:根据实际需要和具体材料的性质,选择合适的材料模型。
常见的模型包括线弹性模型、塑性模型、粘弹性模型等。
2. 编写用户材料子程序:使用合适的编程语言(如Fortran)编写用户材料子程序,实现材料模型的行为。
子程序应包括材料刚度矩阵计算、应力和塑性应变更新等关键计算部分。
3. 软件接口设置:将编写好的用户材料子程序与ABAQUS 软件进行接口设置,以实现子程序与有限元分析的集成。
4. 验证与调试:使用合适的测试用例对子程序进行验证和调试,确保其计算结果与实际情况吻合。
关于用户子程序的东东194673-ABAQUS
![关于用户子程序的东东194673-ABAQUS](https://img.taocdn.com/s3/m/e8a3a96fb84ae45c3b358cb1.png)
以下内容为我最近中的一点总结和整理,希望对大家有所帮助!当用到某个用户子程序时,用户所关心的主要有两方面:一是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变化而变化的。
SUBROUTINEDLOAD(P,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,1 JLTYP,SNAME)INCLUDE 'ABA_PARAM.INC' CDIMENSION TIME(2),COORDS(3)CHARACTER*80 SNAMEPARAMETER (PLOAD=100.E4)IF (KSTEP.EQ.1) THEN !当STEP=1时的荷载大小P=PLOADELSE IF (KSTEP.EQ.2) THEN !当STEP=2时的荷载大小P=COORDS(1)*PLOAD !施加在积分点的荷载P是坐标的函数ELSE IF (KSTEP.EQ.3) THEN !当STEP=3时的荷载大小P=COORDS(1)**2*PLOADELSE IF (KSTEP.EQ.4) THEN !当STEP=4时的荷载大小P=COORDS(1)**3*PLOADELSE IF (KSTEP.EQ.5) THEN !当STEP=5时的荷载大小P=COORDS(1)**4*PLOADEND IFRETURNENDUMAT子程序具有强大的功能,使用UMAT 子程序:(1)可以定义材料的本构关系,使用ABAQUS 材料库中没有包含的材料进行计算,扩充程序功能。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
**
*Assembly, name=Assembly
**
*Instance, name=Part-1-1, part=Part-1
*Node
1, 0.0399999991, 0.0399999991, 0.400000006
2, 0.0399999991, 0.0299999993, 0.400000006
3, 0.0399999991, 0.0199999996, 0.400000006
4, 0.0399999991, 0.00999999978, 0.400000006
5, 0.0399999991,
0., 0.400000006
6, 0.0399999991, 0.0399999991, 0.389999986
C IF (NDI.NE.3) THEN WRITE (6,*) 'THIS UMAT MAY ONLY BE USED FOR ELEMENTS
1 WITH THREE DIRECT STRESS COMPONENTS' CALL XIT ENDIF
C C ELASTIC PROPERTIES
EMOD=PROPS(1) ENU=PROPS(2) EBULK3=EMOD/(1-2*ENU) EG2=EMOD/(1+ENU) EG=EG2/2 EG3=3*EG ELAM=(EBULK3-EG2)/3 C C ELASTIC STIFFNESS C DO K1=1, NDI DO K2=1, NDI DDSDDE(K2, K1)=ELAM END DO DDSDDE(K1, K1)=EG2+ELAM END DO DO K1=NDI+1, NTENS DDSDDE(K1 ,K1)=EG END DO C C CALCULATE STRESS C DO K1=1, NTENS DO K2=1, NTENS STRESS(K2)=STRESS(K2)+DDSDDE(K2, K1)*DSTRAN(K1) END DO END DO C RETURN END 5 INP 输入文件
31, 0.0399999991, 0.0399999991, 0.340000004
32, 0.0399999991, 0.0299999993, 0.340000004
33, 0.0399999991, 0.0199999996, 0.340000004
34, 0.0399999991, 0.00999999978, 0.340000004
15, 0.0399999991,
0., 0.379999995
16, 0.0399999991, 0.0399999991, 0.370000005
17, 0.0399999991, 0.0299999993, 0.370000005
18, 0.0399999991, 0.0199999996, 0.370000005
27, 0.0399999991, 0.0299999993, 0.349999994
28, 0.0399999991, 0.0199999996, 0.349999994
29, 0.0399999991, 0.00999999978, 0.349999994
30, 0.0399999991,
0., 0.349999994
39, 0.0399999991, 0.00999999978, 0.330000013
40, 0.0399999991,
0., 0.330000013
41, 0.0399999991, 0.0399999991, 0.319999993
42, 0.0399999991, 0.0299999993, 0.319999993
55, 0.0399999991,
0., 0.300000012
56, 0.0399999991, 0.0399999991, 0.289999992
57, 0.0399999991, 0.0299999993, 0.289999992
58, 0.0399999991, 0.0199999996, 0.289999992
SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,DDSDDT, 1 DRPLDE,DRPLDT,STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED, 2 CMNAME,NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT, 3 PNEWDT,CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)
include 'aba_param.inc' CHARACTER*8 CMNAME DIMENSION STRESS(NTENS),STATEV(NSTATV),DDSDDE(NTENS,NTENS), 1 DDSDDT(NTENS),DRPLDE(NTENS),STRAN(NTENS),DSTRAN(NTENS), 2 TIME(2),PREDEF(1),DPRED(1),PROPS(NPROPS),COORDS(3),DROT(3,3), 3 DFGRD0(3,3),DFGRD1(3,3) C UMAT FOR ISOTROPIC ELASTICITY C CANNOT BE USED FOR PLANE STRESS C ---------------------------------------------------------------C PROPS(1) - E C PROPS(2) - NU C ----------------------------------------------------------------
63, 0.0399999991, 0.0199999996, 0.280000001
64, 0.0399999991, 0.00999999978, 0.280000001
65, 0.0399999991,
0., 0.280000001
66, 0.0399999991, 0.0399999991, 0.270000011
35, 0.0399999991,
0., 0.340000004
36, 0.0399999991, 0.0399999991, 0.330000013
37, 0.0399999991, 0.0299999993, 0.330000013
38, 0.0399999991, 0.0199999996, 0.330000013
实例:最简单的杆件单轴拉伸,材料本构模型采用线弹性模型, E=206000MPa, ν=0.3。左端约束,右端施加均布力 q=1。
Author: xueweek 希望这个小例子能帮助初学者减少一两天的摸索时间
1 材料本构模型:σ=Eε
在 Property 中定义材料时,在 General 下选中 User Material,输入 206000,0.3。 以上两个数值代表 E、ν,在用户子程序中代表着 PROPS(1)、PROPS(2)。然后在 General 下选中 Depvar,由于该例子中没有使用状态变量,因此在第一项中输入大于 0 的数值即 可。
47, 0.0399999991, 0.0299999993, 0.310000002
48, 0.0399999991, 0.0199999996, 0.310000002
49, 0.0399999991, 0.00999999978, 0.310000002
50, 0.0399999991,
0., 0.310000002
7, 0.0399999991, 0.0299999993, 0.389999986
8, 0.0399999991, 0.0199999996, 0.389999986
9, 0.0399999991, 0.00999999978, 0.389999986
10, 0.0399999991,
0., 0.389999986
3 结果 以下两张图分布是用户子程序和 ABAQUS 自带的材料模型得到的应力
云图,可以看出两种图形完全相同。
4 UMAT 子程序 对于初学者来说,需要注意的是,FORTRAN 对于程序语言格式上的要求。例如,对于
FORTRAN 语言,前六个字符必须空出来,等等。检查语法最好的方法就是在 FORTRAN 编 译器上进行编译。
建模大家都会,故省略
2 ABAQUS 调用用户子程序
大家在论坛上也可以看到很多方法,我这里指介绍最简单的一种。 在 Job Manager 中点击 Edit 选项,在 General 选项的最后一项中选择自己建立好的用户 子程序文件。(注:用户子程序文件可以使用文本编辑器进行编辑,当然也可以用 Fortran 编译器,如果对自己的用户子程序文件的语法不放心,可以先用 Fortan 编译器进行编译, 不过编译前要先建立 Project,关于 Fortran 编译,这里不再介绍)。 完成后,点击 submit 即可进行分析。
*Heading ** Job name: Job-1 Model name: Model-1 ** Generated by: Abaqus/CAE 6.10-1 *Preprint, echo=NO, model=NO, history=NO, contact=NO ** ** PARTS ** *Part, name=Part-1 *End Part **
59, 0.0399999991, 0.00999999978, 0.289999992
60, 0.0399999991,
0., 0.289999992
61, 0.0399999991, 0.0399999991, 0.280000001