您的位置：360文档中心› Abaqus关于用户子程序(simwe)

Abaqus关于用户子程序(simwe)

合集下载

ABAQUS子程序

ABAQUS子程序

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用户子程序VUMAT

Abaqus用户子程序VUMAT
Pipe elements For pipe elements the stretch tensor and the deformation gradient tensor are not available. The axial strain, strainInc(*,1), and the shear strain, strainInc(*,4), associated with twist are provided along with the hoop stress, stressNew(*,2). The hoop stress is predefined based on your pipe internal and external pressure load definitions (PE, PI, HPE, HPI, PENU, and PINU), and it should not be modified here. The thickness stress, stressNew(*,3), is assumed to be zero and any value you assign is ignored. You must define the axial stress, stressNew(*,1), and the shear stress, stressNew(*,4). You must also define hoop strain, strainInc(*,2), and the pipe thickness strain, strainInc(*,3).
Component 1 2 3 4 5 6 7 8 9
2-D Case 3-D Case
Initial calculations and checks

abaqus UAMP用户子程序

abaqus UAMP用户子程序

OverviewUser subroutine UAMP:∙allows you to define the current value of an amplitude definition as a function of time;∙can be used to model control engineering aspects of your system when sensors are used (sensor values are from the beginning of theincrement);∙can use a predefined number of state variables in their definition; and ∙can optionally compute the derivatives and integrals of the amplitude function.Explicit solution dependenceThe solution dependence introduced in this user subroutine is explicit: all data passed in the subroutine for information or to be updated are values at the beginning of that increment.User subroutine interfaceSUBROUTINE UAMP(* ampName, time, ampValueOld, dt, nProps, props, nSvars,* svars, lFlagsInfo,* nSensor, sensorValues, sensorNames, jSensorLookUpTable,* AmpValueNew,* lFlagsDefine,* AmpDerivative, AmpSecDerivative, AmpIncIntegral, * AmpDoubleIntegral)CINCLUDE 'ABA_PARAM.INC'C time indicesparameter (iStepTime = 1,* iTotalTime = 2,* nTime = 2)C flags passed in for informationparameter (iInitialization = 1,* iRegularInc = 2,* iCuts = 3* ikStep = 4* nFlagsInfo = 4)C optional flags to be definedparameter (iComputeDeriv = 1,* iComputeSecDeriv = 2,* iComputeInteg = 3,* iComputeDoubleInteg = 4,* iStopAnalysis = 5,* iConcludeStep = 6,* nFlagsDefine = 6)dimension time(nTime), lFlagsInfo(nFlagsInfo),* lFlagsDefine(nFlagsDefine)dimension jSensorLookUpTable(*)dimension sensorValues(nSensor), svars(nSvars),props(nProps)character*80 sensorNames(nSensor)character*80 ampNameuser coding to define AmpValueNew, andoptionally lFlagsDefine, AmpDerivative, AmpSecDerivative,AmpIncIntegral, AmpDoubleIntegralRETURNENDVariable to be definedAmpValueNewCurrent value of the amplitude.Variables that can be updatedlFlagsDefineInteger flag array to determine whether the computation of additional quantities is necessary or to set step continuation requirements.lFlagsDefine(iComputeDeriv)If set to 1, you must provide thecomputation of the amplitude derivative.The default is 0, which means thatAbaqus computes the derivativeautomatically.lFlagsDefine(iComputeSecDeriv)If set to 1, you must provide thecomputation of the amplitude secondderivative. The default is 0, whichmeans that Abaqus computes thesecond derivative automatically. lFlagsDefine(iComputeInteg)If set to 1, you must provide thecomputation of the amplitudeincremental integral. The default is 0,which means that Abaqus computes theincremental integral automatically.lFlagsDefine(iComputeDoubleInteg)If set to 1, you must provide the computation of the amplitude incremental double integral. The default is 0, which means that Abaqus computes the incremental integral automatically.lFlagsDefine(iStopAnalysis)If set to 1, the analysis will be stoppedand an error message will be issued.The default is 0, which means thatAbaqus will not stop the analysis. lFlagsDefine(iConcludeStep)If set to 1, Abaqus will conclude the stepexecution and advance to the next step(if a next step exists). The default is 0. svarsAn array containing the values of the solution-dependent state variables associated with this amplitude definition. The number of such variablesis nsvars (see above). You define the meaning of these variables.This array is passed into UAMP containing the values of these variables at the start of the current increment. In most cases they should be updated to be the values at the end of the increment.AmpDerivativeCurrent value of the amplitude derivative.AmpSecDerivativeCurrent value of the amplitude second derivative.AmpIncIntegralCurrent value of the amplitude incremental integral.AmpDoubleIntegralCurrent value of the amplitude incremental double integral.Variables passed in for informationampNameUser-specified amplitude name, left justified.time(iStepTime)Current value of step time or frequency.time(iTotalTime)Current value of total time.ampValueOldOld value of the amplitude from the previous increment.dtTime increment.propsUser-specified array of material constants associated with this amplitude definition.nPropsUser-defined number of material constants associated with this amplitude definition.nSvarsUser-defined number of solution-dependent state variables associated with this amplitude definition.lFlagsInfoInteger flag array with information regrading the current call to UAMP.lFlagsInfo(iInitialization) This flag is equal to 1 if UAMP is called from theinitialization phase of the first analysis step andis set to 0 otherwise.lFlagsInfo(iRegularInc)This flag is equal to 1 if UAMP is called from aregular increment and is set to 0 if called fromthe initialization phase of the first analysis step. lFlagsInfo(iCuts)Number of cutbacks in this increment. lFlagsInfo(ikStep)Step number.nSensorTotal number of sensors in the model.sensorValuesArray with sensor values at the end of the previous increment. Each sensor value corresponds to a history output variable associated with the output database request defining the sensor.sensorNamesArray with user-defined sensor names in the entire model, left justified. Each sensor name corresponds to a sensor value provided with the output database request. All names will be converted to uppercase characters if lowercase or mixed-case characters were used in their definition.jSensorLookUpTableVariable that must be passed into the utilityfunctions IGETSENSORID and GETSENSORVALUE.Example: Amplitude definition using sensor and state variablesc user amplitude subroutineSubroutine UAMP(C passed in for information and state variables* ampName, time, ampValueOld, dt, nProps, props, nSvars,* svars, lFlagsInfo,* nSensor, sensorValues, sensorNames,* jSensorLookUpTable,C to be defined* ampValueNew,* lFlagsDefine,* AmpDerivative, AmpSecDerivative, AmpIncIntegral, * AmpDoubleIntegral)include 'aba_param.inc'C svars - additional state variables, similar to (V)UEL dimension sensorValues(nSensor), svars(nSvars),* props(nProps)character*80 sensorNames(nSensor)character*80 ampNameC time indicesparameter( iStepTime = 1,* iTotalTime = 2,* nTime = 2)C flags passed in for informationparameter( iInitialization = 1,* iRegularInc = 2,* iCuts = 3,* ikStep = 4,* nFlagsInfo = 4)C optional flags to be definedparameter( iComputeDeriv = 1,* iComputeSecDeriv = 2,* iComputeInteg = 3,* iComputeDoubleInteg = 4,* iStopAnalysis = 5,* iConcludeStep = 6,* nFlagsDefine = 6)parameter( tStep=0.18d0, tAccelerateMotor = .00375d0, * omegaFinal=23.26d0,* zero=0.0d0, one=1.0d0, two=2.0d0, four=4.0d0)dimension time(nTime), lFlagsInfo(nFlagsInfo),* lFlagsDefine(nFlagsDefine)dimension jSensorLookUpTable(*)lFlagsDefine(iComputeDeriv) = 1lFlagsDefine(iComputeSecDeriv) = 1lFlagsDefine(iComputeInteg) = 1lFlagsDefine(iComputeDoubleInteg) = 1c get sensor valuevTrans_CU1 = GetSensorValue('HORIZ_TRANSL_MOTION',* jSensorLookUpTable,* sensorValues)if (ampName(1:22) .eq. 'MOTOR_WITH_STOP_SENSOR' ) then if (lFlagsInfo(iInitialization).eq.1) thenAmpSecDerivative = zeroAmpDerivative = omegaFinal/tAccelerateMotor ampValueNew = zeroAmpIncIntegral = zeroAmpDoubleIntegral = zerosvars(1) = zerosvars(2) = zeroelsetim = time(iStepTime)c ramp up the angular rot velocity of thec electric motorc after which hold constantif (tim .le. tAccelerateMotor) thenAmpSecDerivative = zeroAmpDerivative =omegaFinal/tAccelerateMotorampValueNew =omegaFinal*tim/tAccelerateMotorAmpIncIntegral =dt*(ampValueOld+ampValueNew)/twoAmpDoubleIntegral =dt**2*(ampValueOld+ampValueNew)/fourelseAmpSecDerivative = zeroAmpDerivative = zeroampValueNew = omegaFinalAmpIncIntegral =dt*(ampValueOld+ampValueNew)/twoAmpDoubleIntegral =dt**2*(ampValueOld+ampValueNew)/fourend ifc retrieve old sensor valuevTrans_CU1_old = svars(1)c detect a zero crossing and count the number of crossings if (vTrans_CU1_old*vTrans_CU1 .le. zero .and.* tim .gt. tAccelerateMotor ) thensvars(2) = svars(2) + oneend ifnrCrossings = int(svars(2))c stop the motor if sensor crosses zero the second time if (nrCrossings.eq.2) thenampValueNew = zerolFlagsDefine(iConcludeStep)=1end ifc store sensor valuesvars(1) = vTrans_CU1end ifend ifreturnend。

ABAQUS用户子程序学习小结

ABAQUS用户子程序学习小结

1。

1 FORTRAN语言不区分大小写，文件扩展名为for，正文从第7列开始，第1～5列是标号区,第六列是续行标志区.1.2 FORTRAN语言中的“I-N规则”：I、J、K、L、M、N开头的为整型变量,其他开头为实型变量;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=RECT 0.2, 5.0。

ABAQUS学习零碎笔记(转自Simwe)

ABAQUS学习零碎笔记(转自Simwe)

1.接触中设置Adjust的理解：这个命令主要还是用来初始化接触的。

在分析开始之前，调整接触面中节点的初始位置，且不产生任何应力和应变。

在分析过程中，由于残留的初始过盈引起的应变将被施加在接触面上。

模型的尺寸往往会存在数值误差，所以设置一个位置误差限度，用来调整从面节点的初始坐标，ADJUST=位置误差限度，其含义是：如果从面节点与主面的距离小于此限度，ABAQUS将调整这些节点的初始坐标，使其与主面的距离为0.这种调整不产生任何压力、应力、应变。

Explicit不允许接触表面的初始过盈，分析开始前，接触面上的节点将被自动调整，以删除任何初始过盈，在随后的分析中，这样的调整将引起应变。

2.使用INTERFERENCE（干涉）来定义过盈接触：Edit Interaction》底部Interference fit负值表示过盈量，正值表示间隙量。

类似于载荷，只能在后续分析步中定义，不能在初始分析步中定义。

3.CLERENCE（间隙）可以定义两个接触面之间的初始过盈量和间隙量，它只适用于小滑移，并且不需要使用ADJUST来调整从面节点的初始位置。

4.特征：记录了设计目的，并包含几何信息，同时也是管理几何体的行为的规则。

ABAQUS中导入的几何体是没用特征的，要删除不重要的细节。

5.View=》ODB Display Options =》Sweep and Extrude6.CPRESS接触压强7．COPEN从面上节点与主面的距离8.ALE自适应网格：Step=>Other=>Adaptive Mesh Domain9.计算代价估算：隐式：自由度数目的平方显式：自由度正比10.软接触：用指数或者表格形式表达的应力-距离关系11.传说ABAQUS默认的幅值曲线是从1下降到0的。

是吗？我觉得6.9的版本好像不是这样。

给一个棱柱施加扭转力矩，载荷采用Ramp，变形是逐步增大的。

12.Visualization》Tools》Job Diagnostics13.修改严重不连续迭代尝试次数：Step>Other>General Solution Controls>Edit>Specify>Time Incrementation>More>Is>1214.接触问题中，90度圆角至少划分是个单元。

ABAQUS用户子程序小例子

ABAQUS用户子程序小例子

实例：最简单的杆件单轴拉伸，材料本构模型采用线弹性模型，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 的数值即可。

建模大家都会，故省略2 ABAQUS 调用用户子程序大家在论坛上也可以看到很多方法，我这里指介绍最简单的一种。

在Job Manager 中点击Edit 选项，在General 选项的最后一项中选择自己建立好的用户子程序文件。

（注：用户子程序文件可以使用文本编辑器进行编辑，当然也可以用Fortran编译器，如果对自己的用户子程序文件的语法不放心，可以先用Fortan 编译器进行编译，不过编译前要先建立Project，关于Fortran 编译，这里不再介绍）。

完成后，点击submit 即可进行分析。

3 结果以下两张图分布是用户子程序和ABAQUS 自带的材料模型得到的应力云图，可以看出两种图形完全相同。

4 UMAT 子程序对于初学者来说，需要注意的是，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 ----------------------------------------------------------------CIF (NDI.NE.3) THENWRITE (6,*) 'THIS UMAT MAY ONLY BE USED FOR ELEMENTS1 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 DOCRETURNEND5 INP 输入文件*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****** A SSEMBLY ***Assembly, name=Assembly ***Instance, name=Part-1-1, part=Part-1 *Node1, 0.0399999991, 0.0399999991, 0.4000000062, 0.0399999991, 0.029*******, 0.4000000063, 0.0399999991, 0.0199999996, 0.4000000064, 0.0399999991, 0.00999999978, 0.4000000065, 0.0399999991, 0., 0.4000000066, 0.0399999991, 0.0399999991, 0.3899999867, 0.0399999991, 0.029*******, 0.3899999868, 0.0399999991, 0.0199999996, 0.3899999869, 0.0399999991, 0.00999999978, 0.38999998610, 0.0399999991, 0., 0.38999998611, 0.0399999991, 0.0399999991, 0.37999999512, 0.0399999991, 0.029*******, 0.37999999513, 0.0399999991, 0.0199999996, 0.37999999514, 0.0399999991, 0.00999999978, 0.37999999515, 0.0399999991, 0., 0.37999999516, 0.0399999991, 0.0399999991, 0.37000000517, 0.0399999991, 0.029*******, 0.37000000518, 0.0399999991, 0.0199999996, 0.37000000519, 0.0399999991, 0.00999999978, 0.37000000520, 0.0399999991, 0., 0.37000000521, 0.0399999991, 0.0399999991, 0.36000001422, 0.0399999991, 0.029*******, 0.36000001423, 0.0399999991, 0.0199999996, 0.36000001424, 0.0399999991, 0.00999999978, 0.36000001425, 0.0399999991, 0., 0.36000001426, 0.0399999991, 0.0399999991, 0.34999999427, 0.0399999991, 0.029*******, 0.34999999428, 0.0399999991, 0.0199999996, 0.34999999429, 0.0399999991, 0.00999999978, 0.34999999430, 0.0399999991, 0., 0.34999999431, 0.0399999991, 0.0399999991, 0.34000000432, 0.0399999991, 0.029*******, 0.34000000433, 0.0399999991, 0.0199999996, 0.34000000434, 0.0399999991, 0.00999999978, 0.34000000435, 0.0399999991, 0., 0.34000000436, 0.0399999991, 0.0399999991, 0.33000001337, 0.0399999991, 0.029*******, 0.33000001338, 0.0399999991, 0.0199999996, 0.33000001339, 0.0399999991, 0.00999999978, 0.33000001340, 0.0399999991, 0., 0.33000001341, 0.0399999991, 0.0399999991, 0.31999999342, 0.0399999991, 0.029*******, 0.31999999343, 0.0399999991, 0.0199999996, 0.31999999344, 0.0399999991, 0.00999999978, 0.31999999345, 0.0399999991, 0., 0.31999999346, 0.0399999991, 0.0399999991, 0.31000000247, 0.0399999991, 0.029*******, 0.31000000248, 0.0399999991, 0.0199999996, 0.31000000249, 0.0399999991, 0.00999999978, 0.31000000250, 0.0399999991, 0., 0.31000000251, 0.0399999991, 0.0399999991, 0.30000001252, 0.0399999991, 0.029*******, 0.30000001253, 0.0399999991, 0.0199999996, 0.30000001254, 0.0399999991, 0.00999999978, 0.30000001255, 0.0399999991, 0., 0.30000001256, 0.0399999991, 0.0399999991, 0.28999999257, 0.0399999991, 0.029*******, 0.28999999258, 0.0399999991, 0.0199999996, 0.28999999259, 0.0399999991, 0.00999999978, 0.28999999260, 0.0399999991, 0., 0.28999999261, 0.0399999991, 0.0399999991, 0.28000000162, 0.0399999991, 0.029*******, 0.28000000163, 0.0399999991, 0.0199999996, 0.28000000164, 0.0399999991, 0.00999999978, 0.28000000165, 0.0399999991, 0., 0.28000000166, 0.0399999991, 0.0399999991, 0.27000001167, 0.0399999991, 0.029*******, 0.27000001168, 0.0399999991, 0.0199999996, 0.27000001169, 0.0399999991, 0.00999999978, 0.27000001170, 0.0399999991, 0., 0.27000001171, 0.0399999991, 0.0399999991, 0.2599999972, 0.0399999991, 0.029*******, 0.2599999973, 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.23999999582, 0.0399999991, 0.029*******, 0.23999999583, 0.0399999991, 0.0199999996, 0.23999999584, 0.0399999991, 0.00999999978, 0.23999999585, 0.0399999991, 0., 0.239999995 86, 0.0399999991, 0.0399999991, 0.23000000487, 0.0399999991, 0.029*******, 0.23000000488, 0.0399999991, 0.0199999996, 0.23000000489, 0.0399999991, 0.00999999978, 0.23000000490, 0.0399999991, 0., 0.230000004 91, 0.0399999991, 0.0399999991, 0.21999999992, 0.0399999991, 0.029*******, 0.21999999993, 0.0399999991, 0.0199999996, 0.21999999994, 0.0399999991, 0.00999999978, 0.21999999995, 0.0399999991, 0., 0.219999999 96, 0.0399999991, 0.0399999991, 0.20999999397, 0.0399999991, 0.029*******, 0.20999999398, 0.0399999991, 0.0199999996, 0.20999999399, 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.0900000036157, 0.0399999991, 0.029*******, 0.0900000036158, 0.0399999991, 0.0199999996, 0.0900000036159, 0.0399999991, 0.00999999978, 0.0900000036160, 0.0399999991, 0., 0.0900000036 161, 0.0399999991, 0.0399999991, 0.0799999982162, 0.0399999991, 0.029*******, 0.0799999982163, 0.0399999991, 0.0199999996, 0.0799999982164, 0.0399999991, 0.00999999978, 0.0799999982165, 0.0399999991, 0., 0.0799999982 166, 0.0399999991, 0.0399999991, 0.0700000003167, 0.0399999991, 0.029*******, 0.0700000003168, 0.0399999991, 0.0199999996, 0.0700000003169, 0.0399999991, 0.00999999978, 0.0700000003170, 0.0399999991, 0., 0.0700000003 171, 0.0399999991, 0.0399999991, 0.0599999987172, 0.0399999991, 0.029*******, 0.0599999987173, 0.0399999991, 0.0199999996, 0.0599999987174, 0.0399999991, 0.00999999978, 0.0599999987175, 0.0399999991, 0., 0.0599999987 176, 0.0399999991, 0.0399999991, 0.0500000007177, 0.0399999991, 0.029*******, 0.0500000007178, 0.0399999991, 0.0199999996, 0.0500000007179, 0.0399999991, 0.00999999978, 0.0500000007180, 0.0399999991, 0., 0.0500000007 181, 0.0399999991, 0.0399999991, 0.0399999991182, 0.0399999991, 0.029*******, 0.0399999991183, 0.0399999991, 0.0199999996, 0.0399999991184, 0.0399999991, 0.00999999978, 0.0399999991185, 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.0199999996192, 0.0399999991, 0.029*******, 0.0199999996193, 0.0399999991, 0.0199999996, 0.0199999996194, 0.0399999991, 0.00999999978, 0.0199999996195, 0.0399999991, 0., 0.0199999996 196, 0.0399999991, 0.0399999991, 0.00999999978197, 0.0399999991, 0.029*******, 0.00999999978198, 0.0399999991, 0.0199999996, 0.00999999978199, 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.400000006207, 0.029*******, 0.029*******, 0.400000006208, 0.029*******, 0.0199999996, 0.400000006209, 0.029*******, 0.00999999978, 0.400000006 210, 0.029*******, 0., 0.400000006 211, 0.029*******, 0.0399999991, 0.389999986212, 0.029*******, 0.029*******, 0.389999986213, 0.029*******, 0.0199999996, 0.389999986214, 0.029*******, 0.00999999978, 0.389999986 215, 0.029*******, 0., 0.389999986 216, 0.029*******, 0.0399999991, 0.379999995217, 0.029*******, 0.029*******, 0.379999995218, 0.029*******, 0.0199999996, 0.379999995219, 0.029*******, 0.00999999978, 0.379999995 220, 0.029*******, 0., 0.379999995 221, 0.029*******, 0.0399999991, 0.370000005222, 0.029*******, 0.029*******, 0.370000005223, 0.029*******, 0.0199999996, 0.370000005224, 0.029*******, 0.00999999978, 0.370000005 225, 0.029*******, 0., 0.370000005 226, 0.029*******, 0.0399999991, 0.360000014227, 0.029*******, 0.029*******, 0.360000014228, 0.029*******, 0.0199999996, 0.360000014229, 0.029*******, 0.00999999978, 0.360000014 230, 0.029*******, 0., 0.360000014 231, 0.029*******, 0.0399999991, 0.349999994232, 0.029*******, 0.029*******, 0.349999994233, 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*******394, 0.029*******, 0.00999999978, 0.029*******395, 0.029*******, 0., 0.029******* 396, 0.029*******, 0.0399999991, 0.0199999996397, 0.029*******, 0.029*******, 0.0199999996398, 0.029*******, 0.0199999996, 0.0199999996399, 0.029*******, 0.00999999978, 0.0199999996400, 0.029*******, 0., 0.0199999996 401, 0.029*******, 0.0399999991, 0.00999999978402, 0.029*******, 0.029*******, 0.00999999978403, 0.029*******, 0.0199999996, 0.00999999978404, 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.400000006412, 0.0199999996, 0.029*******, 0.400000006413, 0.0199999996, 0.0199999996, 0.400000006414, 0.0199999996, 0.00999999978, 0.400000006 415, 0.0199999996, 0., 0.400000006 416, 0.0199999996, 0.0399999991, 0.389999986417, 0.0199999996, 0.029*******, 0.389999986418, 0.0199999996, 0.0199999996, 0.389999986419, 0.0199999996, 0.00999999978, 0.389999986 420, 0.0199999996, 0., 0.389999986 421, 0.0199999996, 0.0399999991, 0.379999995422, 0.0199999996, 0.029*******, 0.379999995423, 0.0199999996, 0.0199999996, 0.379999995424, 0.0199999996, 0.00999999978, 0.379999995 425, 0.0199999996, 0., 0.379999995 426, 0.0199999996, 0.0399999991, 0.370000005427, 0.0199999996, 0.029*******, 0.370000005428, 0.0199999996, 0.0199999996, 0.370000005429, 0.0199999996, 0.00999999978, 0.370000005 430, 0.0199999996, 0., 0.370000005 431, 0.0199999996, 0.0399999991, 0.360000014432, 0.0199999996, 0.029*******, 0.360000014433, 0.0199999996, 0.0199999996, 0.360000014434, 0.0199999996, 0.00999999978, 0.360000014 435, 0.0199999996, 0., 0.360000014 436, 0.0199999996, 0.0399999991, 0.349999994437, 0.0199999996, 0.029*******, 0.349999994438, 0.0199999996, 0.0199999996, 0.349999994439, 0.0199999996, 0.00999999978, 0.349999994 440, 0.0199999996, 0., 0.349999994 441, 0.0199999996, 0.0399999991, 0.340000004442, 0.0199999996, 0.029*******, 0.340000004443, 0.0199999996, 0.0199999996, 0.340000004444, 0.0199999996, 0.00999999978, 0.340000004 445, 0.0199999996, 0., 0.340000004 446, 0.0199999996, 0.0399999991, 0.330000013447, 0.0199999996, 0.029*******, 0.330000013448, 0.0199999996, 0.0199999996, 0.330000013449, 0.0199999996, 0.00999999978, 0.330000013 450, 0.0199999996, 0., 0.330000013 451, 0.0199999996, 0.0399999991, 0.319999993452, 0.0199999996, 0.029*******, 0.319999993453, 0.0199999996, 0.0199999996, 0.319999993454, 0.0199999996, 0.00999999978, 0.319999993 455, 0.0199999996, 0., 0.319999993 456, 0.0199999996, 0.0399999991, 0.310000002457, 0.0199999996, 0.029*******, 0.310000002458, 0.0199999996, 0.0199999996, 0.310000002459, 0.0199999996, 0.00999999978, 0.310000002 460, 0.0199999996, 0., 0.310000002 461, 0.0199999996, 0.0399999991, 0.300000012462, 0.0199999996, 0.029*******, 0.300000012463, 0.0199999996, 0.0199999996, 0.300000012464, 0.0199999996, 0.00999999978, 0.300000012 465, 0.0199999996, 0., 0.300000012 466, 0.0199999996, 0.0399999991, 0.289999992467, 0.0199999996, 0.029*******, 0.289999992468, 0.0199999996, 0.0199999996, 0.289999992469, 0.0199999996, 0.00999999978, 0.289999992 470, 0.0199999996, 0., 0.289999992 471, 0.0199999996, 0.0399999991, 0.280000001472, 0.0199999996, 0.029*******, 0.280000001473, 0.0199999996, 0.0199999996, 0.280000001474, 0.0199999996, 0.00999999978, 0.280000001 475, 0.0199999996, 0., 0.280000001 476, 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.109999999557, 0.0199999996, 0.029*******, 0.109999999558, 0.0199999996, 0.0199999996, 0.109999999559, 0.0199999996, 0.00999999978, 0.109999999 560, 0.0199999996, 0., 0.109999999 561, 0.0199999996, 0.0399999991, 0.100000001562, 0.0199999996, 0.029*******, 0.100000001563, 0.0199999996, 0.0199999996, 0.100000001564, 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.0900000036570, 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.0799999982575, 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.0700000003580, 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.0599999987585, 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.0500000007590, 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.0399999991595, 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.0199999996605, 0.0199999996, 0., 0.0199999996 606, 0.0199999996, 0.0399999991, 0.00999999978607, 0.0199999996, 0.029*******, 0.00999999978608, 0.0199999996, 0.0199999996, 0.00999999978609, 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关于用户子程序（simwe）

Abaqus关于用户子程序（simwe）

Abaqus关于用户子程序（simwe）序，但可以调用用户自己编写的FORTRAN子程序和ABAQUS应用程序。

当用户编写FORTRAN子程序时，建议子程序名以K开头，以免和ABAQUS内部程序冲突。

2．当用户在用户子程序中利用OPEN打开外部文件时，要注意以下两点：一是设备号的选择是有限制的，只能取15－18和大于100的设备号，其余的都已被ABAQUS占用。

二是用户需提供外部文件的绝对路径而不是相对路径。

3．ABAQUS 应用程序必须由用户子程序调用。

当用到某个用户子程序时，用户所关心的主要有两方面：一是ABAQUS提供的用户子程序的接口参数。

有些参数是ABAQUS传到用户子程序中的，例如SUBROUTINE DLOAD中的KSTEP，KINC，COORDS；有些是需要用户自己定义的，例如F。

二是ABAQUS何时调用该用户子程序，对于不同的用户子程序ABAQUS调用的时间是不同的。

有些是在每个STEP的开始，有的是STEP结尾，有的是在每个INCREMENT的开始等等。

当ABAQUS调用用户子程序是，都会把当前的STEP和INCREMENT利用用户子程序的两个实参KSTEP和KINC传给用户子程序，用户可编个小程序把它们输出到外部文件中，这样对ABAQUS何时调用该用户子程序就会有更深的了解。

下面就选出几个常用的用户子程序和应用程序进行详细解释：一．SUBROUTINEDLOAD(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传到用户子程序中的。

ABAQUS的inp文件精讲

ABAQUS的inp文件精讲
列的时步层序.每一步就是一个响应(静态加载,动态响应等)，时步的定义包括过程类型(比 如静态应力分析,瞬时传热分析等)对于时间积分的控制参数或者非线性解过程,加载和输出 要求.At a minimum the model consists of the following information: geometry, element section properties, material data, loads and boundary conditions, analysis type, and output requests.
*ELEMENT, TYPE=SPRINGA, ELSET="One element" 1,1,2 *SPRING, ELSET="One eห้องสมุดไป่ตู้ement" 1.0E-5, *NSET, ELSET="One element", NSET=NODESET *BOUNDARY nodeset,1,2 4．数据行重复: 数据行可以重复,也就是说每行数据可以有一行响应的变量也可以有几行。同样也可以有多 行数据行,对应各自的变量行.如： *ELASTIC, TYPE=ISOTROPIC 200.E3, 0.3, 20. 150.E3, 0.35, 400. 80.E3, 0.42, 700. 定义了一个零件的材料性质,均质,线弹性,在不同应力下的扬式模量和泊松比. 三、从外存储器中引入模型或者历史数据 The *INCLUDE option can be used to specify an external file that contains a portion of the ABAQUS input file. This file can include model and history definition data, comment lines, and other *INCLUDE references. When the option is encountered, ABAQUS will immediately process the input data within the file specified by the INPUT parameter. When the end-of-file is reached, ABAQUS will return to processing the original file. A maximum of five levels of nested *INCLUDE references can be used. The option

abaqus用户子程序vusdfld单元 积分点出的极限应变

abaqus用户子程序vusdfld单元 积分点出的极限应变

在ABAQUS中，用户子程序（User Subroutine）是用于自定义模型行为的强大工具。

如果你想在VUSDFLD单元的积分点上计算极限应变，你需要编写一个用户子程序来实现这一功能。

以下是一个简单的示例，展示了如何在ABAQUS的Python脚本中定义一个用户子程序来计算极限应变。

请注意，这个示例假设你已经熟悉ABAQUS的Python API和有限元分析的基本概念。

首先，你需要导入ABAQUS的Python模块：
如果你需要更具体的帮助，例如如何编写计算极限应变的代码，或者如何将子程序集成到你的ABAQUS模型中，请提供更多详细信息。

最新ABAQUS用户子程序

最新ABAQUS用户子程序

ABAQUS用户子程序123ABAQUS/Standard subroutines：41.CREEP: Define time-dependent, viscoplastic behavior (creep and 5swelling).6定义和时间相关的、粘塑性的运动（蠕变和膨胀）2. DFLOW: Define nonuniform pore fluid velocity in a 78consolidation analysis.9在压实分析中，定义非均匀孔隙流速度3. DFLUX: Define nonuniform distributed flux in a heat transfer1011or mass diffusion analysis.在热传递和质量扩散分析中，定义非均匀的分布流量12134. DISP: Specify prescribed boundary conditions.14指定规定的边界条件5. DLOAD: Specify nonuniform distributed loads.1516指定非均匀的分布荷载176. FILM: Define nonuniform film coefficient and associated sink 18temperatures for heat19transfer analysis.20对热传递分析指定非均匀的膜层散热系数和联合的散热器温度7. FLOW: Define nonuniform seepage coefficient and associated2122sink pore pressure for23consolidation analysis.对压实分析定义非均匀的渗流系数和渗入孔隙压力24258. FRIC: Define frictional behavior for contact surfaces.26对接触面定义摩擦279. GAPCON: Define conductance between contact surfaces or nodes 28in a fully coupled29temperature-displacement analysis or pure heat transfer analysis.30在一个完全耦合的温度—置换分析或者是纯热传递分析中，定义31接触面或节点间的导热系数。

abaqus调用子程序的方法

abaqus调用子程序的方法

Abaqus运行子程序方法试了好几种调用子程序的方法最后总算找到了最简单的了，非常感谢仿真论坛上的网友的分享：首先需要装Microsoft Visual Studio（简写mvs），再装Intel Visual Fortran（IVF）。

通常的组合是Abaqus v6.10.1可使用：MVS2008+IVF10.1 MVS2008+IVF11.1 这二种组合最好都使用英文版的。

照着附录1 方法装一般没什么问题。

但是一般都是先装abaqus后才想到去装MVS以及IVF的。

很多时候都把MVS 以及IVF装好以后才想到去验证的，而且MVS很多都是中文的，MVS不仅装的慢卸载残留也很多好麻烦。

搞了好久都验证不了，看了附录二的方法下面2个网友的回复后豁然开朗，一切都是那么简单。

附录二的方法的主要作用可能就是将VC++和Fortran的环境变量一并启动了，参考过以前版上的安装教程，通过在ABAQUS 的快捷方式中添加命令，亦能够达到同样效果，同样也是不能通过verification，但是可以运行子程序，通过在ABAQUS CAE的快捷方式的目标中，添加"D:\Program Files (x86)\Intel\Compiler\Fortran\10.1.034\em64t\bin\ifortvars.bat" && （跟附录1中的修改快捷方式的目标一样）达到了同样的效果（摘自ppengine，tcboywjr网友的回复）。

只需在搜索栏中搜索ifortvars.bat，再把路径添加上去照上面的形式修改便可。

本人用的是win7系统装的是abaqus6.10.mvs2008（中文版）IVF11，abaqus好早就装了，装MVS2008之前装了mvs2010，发现2010不好装ivf，然后又卸了装mvs2008中文，验证时出现附录Ⅰ的情况，就C++项通不过还以为装C++2010就好了，又重装了一下mvs2010的C++（囧），发现还是用不了，校园网速太慢又不想下mvs2008（英文版），最后在附录2的方法下面发现了两位网友的分享非常感谢。

浅谈ABAQUS用户子程序

浅谈ABAQUS用户子程序

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子程序

1.1.37 UMATHTUser subroutine to define a material's thermalbehavior.Product: Abaqus/StandardWarning: The use of this subroutine generally requires considerable expertise. You are cautioned that the implementation of any realistic thermal model requires significant development and testing. Initial testing on models with few elements under a variety of boundary conditions is strongly recommended.References●“User-defined thermal material behavior,”Section 23.8.2 of the Abaqus Analysis User's Manual● *USER MATERIAL●“Freezing of a square solid: the two-dimensional Stefan problem,”Section 1.6.2 of the Abaqus Benchmarks Manual●“UMATHT,”Section 4.1.22 of the Abaqus Verification ManualOverviewUser subroutine UMATHT:● can be used to define the thermal constitutive behavior of the material as well as internal heat generation during heat transfer processes;● will be called at all material calculation points of elements for which the material definition includes a user-defined thermal material behavior;● can be used with the procedures discussed in “Heat transfer analysis procedures: overview,”Section 6.5.1 of the Abaqus Analysis User's Manual;● can use solution-dependent state Variables;● must define the internal energy per unit mass and its variation with respect to temperature and to spatial gradients of temperature;● must define the heat flux vector and its variation with respect to temperature and to gradients of temperature;● must update the solution-dependent state Variables to their values at the end of the increment;● can be used in conjunction with user subroutine USDFLD to redefine any field Variables before they are passed in; and● is described further in “User-defined thermal material behavior,”Section 23.8.2 of the Abaqus Analysis User's Manual.Use of subroutine UMATHT with coupled temperature-displacement elementsUser subroutine UMATHT should be used only with reduced-integration or modified coupled temperature-displacement elements if the mechanical and thermal fields are not coupled through plastic dissipation. No such restriction exists with fully integrated coupled temperature-displacement elements.User subroutine interfaceSUBROUTINE UMATHT(U,DUDT,DUDG,FLUX,DFDT,DFDG,1 STATEV,TEMP,DTEMP,DTEMDX,time,Dtime,PREDEF,DPRED,2 CMNAME,NTGRD,NSTATV,PROPS,NPROPS,COORDS,PNEWDT,3 NOEL,NPT,LAYER,KSPT,KSTEP,KINC)CINCLUDE 'ABA_PARAM.INC'CCHARACTER*80 CMNAMEDIMENSION DUDG(NTGRD),FLUX(NTGRD),DFDT(NTGRD),1 DFDG(NTGRD,NTGRD),STATEV(NSTATV),DTEMDX(NTGRD),2 time(2),PREDEF(1),DPRED(1),PROPS(NPROPS),COORDS(3)user coding to define U,DUDT,DUDG,FLUX,DFDT,DFDG,and possibly update STATEV, PNEWDTRETURNENDVariables to be definedUInternal thermal energy per unit mass, U, at the end of increment. This variable is passed in as the value at the start of the increment and must be updated to its value at the end of the increment.DUDTVariation of internal thermal energy per unit mass with respect to temperature,, evaluated at the end of the increment.DUDG(NTGRD)Variation of internal thermal energy per unit mass with respect to the spatial gradients of temperature,, at the end of the increment. The size of this array depends on the value of NTGRD as defined below. This term is typically zero in classical heat transfer analysis.FLUX(NTGRD)Heat flux vector,, at the end of the increment. This variable is passed in with the values at the beginning of the increment and must be updated to the values at the end of the increment.DFDT(NTGRD)Variation of the heat flux vector with respect to temperature,, evaluated at the end of the increment.DFDG(NTGRD,NTGRD)Variation of the heat flux vector with respect to the spatial gradients of temperature,, at the end of the increment. The size of this array depends on the value of NTGRD as defined below.Variables that can be updatedSTATEV(NSTATV)An array containing the solution-dependent state Variables.In an uncoupled heat transfer analysis STATEV is passedinto UMATHT with the values of these Variables at the beginning of the increment. However, any changes in STATEV made in usersubroutine USDFLD will be included in the values passed into UMATHT, since USDFLD is called before UMATHT. In addition, if UMATHT is being used in a fully coupled temperature-displacement analysis and user subroutine CREEP, user subroutine UEXPAN, user subroutine UMAT, or user subroutine UTRS is used to define the mechanical behavior of the material, those routines are called before this routine; therefore, any updating of STATEV done in CREEP, UEXPAN, UMAT, or UTRS will be included in the values passed into UMATHT.In all cases STATEV should be passed back from UMATHT as the values of the state Variables at the end of the current increment.PNEWDTratio of suggested new time increment to the time increment being used (DTIME, see below). This variable allows you to provide input to the automatic time incrementation algorithms in Abaqus/Standard (if automatic time incrementation is chosen).PNEWDT is set to a large value before each call to UMATHT.If PNEWDT is redefined to be less than 1.0, Abaqus/Standard must abandon the time increment and attempt it again with a smaller time increment. The suggested new time increment provided to the automatic time integration algorithms is PNEWDT ×DTIME, wherethe PNEWDT used is the minimum value for all calls to user subroutines that allow redefinition of PNEWDT for this iteration.If PNEWDT is given a value that is greater than 1.0 for all calls to user subroutines for this iteration and the increment converges in this iteration, Abaqus/Standard may increase the time increment. The suggested new time increment provided to the automatic time integration algorithms is PNEWDT ×DTIME, where the PNEWDT used is the minimum value for all calls to user subroutines for this iteration.If automatic time incrementation is not selected in the analysis procedure, values of PNEWDT that are greater than 1.0 will be ignored and values of PNEWDT that are less than 1.0 will cause the job to terminate.Variables passed in for informationTEMPTemperature at the start of the increment.DTEMPIncrement of temperature.DTEMDX(NTGRD)Current values of the spatial gradients of temperature,TIME(1)Value of step time at the beginning of the current increment.time(2)Value of total time at the beginning of the current increment.DTIMETime increment.PREDEFarray of interpolated values of predefined field Variables at this point at the start of the increment, based on the values read in at the nodes.DPREDarray of increments of predefined field Variables.CMNAMEUser-defined material name, left justified.NTGRDnumber of spatial gradients of temperature.NSTATVnumber of solution-dependent state variables associated with this material type (defined as described in “Allocating space” in “Usersubroutines: overview,”Section 15.1.1 of the Abaqus Analysis User's Manual).PROPS(NPROPS)User-specified array of material constants associated with this user material.NPROPSUser-defined number of material constants associated with this user material.COORDSAn array containing the coordinates of this point. These are the current coordinates in a fully coupled temperature-displacement analysis if geometric nonlinearity is accounted for during the step(see “Procedures: overview,”Section 6.1.1 of the Abaqus Analysis User's Manual); otherwise, the array contains the original coordinates of the point.NOEL :Element number.NPT :Integration point number.LAYER :Layer number (for composite shells and layered solids).KSPT: Section point number within the current layer.KSTEP: Step number.KINC: Increment number.Example: Using more than one user-defined thermal material modelTo use more than one user-defined thermal material model, the variable CMNAME can be tested for different material names inside user subroutine UMATHT, as illustrated below:IF (CMNAME(1:4) .EQ. 'MAT1') THENCALL UMATHT_MAT1(argument_list)ELSE IF(CMNAME(1:4) .EQ. 'MAT2') THENCALL UMATHT_MAT2(argument_list)END IFUMATHT_MAT1 and UMATHT_MAT2 are the actual user material subroutines containing the constitutive material models for each material MAT1 and MAT2, respectively. Subroutine UMATHT merely acts as a directory here. The argument list can be the same as that used in subroutine UMATHT.Example: Uncoupled heat transferAs a simple example of the coding of user subroutine UMATHT, consider uncoupled heat transfer analysis in a material. The equationsfor this case are developed here, and the corresponding UMATHT is given. This problem can also be solved by specifying thermal conductivity, specific heat, density, and internal heat generation directly.First, the equations for an uncoupled heat transfer analysis are outlined.The basic energy balance iswhere V is the volume of solid material with surface area S,is the density of the material,is the material time rate of the internal thermal energy, q is the heat flux per unit area of the body flowing into the body, and r is the heat supplied externally into the body per unit volume.A heat flux vectoris defined such thatwhereis the unit outward normal to the surface S. Introducing the above relation into the energy balance equation and using the divergence theorem, the following relation is obtained:The corresponding weak form is given bywhereis the temperature gradient andis an arbitrary variational field satisfying the essential boundary conditions.Introducing the backward difference integration algorithm:the weak form of the energy balance equation becomesThis nonlinear system is solved using Newton's method.In the above equations the thermal constitutive behavior of the material is given byAndwhereare state Variables.The Jacobian for Newton's method is given by (after dropping the subscriptson U)The thermal constitutive behavior for this example is now defined. We assume a constant specific heat for the material. The heat conduction in the material is assumed to be governed by Fourier's law.The internal thermal energy per unit mass is defined aswithwhere c is the specific heat of the material andFourier's law for heat conduction is given aswhereis the thermal conductivity matrix andis position, so thatandThe assumption of conductivity without any temperature dependence implies thatNo state variables are needed for this material, so the allocation of space for them is not necessary.A thermal user material definition can be used to read in the two constants for our simple case, namely the specific heat, c, and the coefficient of thermal conductivity, k, so thatSUBROUTINE UMATHT(U,DUDT,DUDG,FLUX,DFDT,DFDG,1 STATEV,TEMP,DTEMP,DTEMDX,time,DTIME,PREDEF,DPRED,2 CMNAME,NTGRD,NSTATV,PROPS,NPROPS,COORDS,PNEWDT,3 NOEL,NPT,LAYER,KSPT,KSTEP,KINC)CINCLUDE 'ABA_PARAM.INC'CCHARACTER*80 CMNAMEDIMENSION DUDG(NTGRD),FLUX(NTGRD),DFDT(NTGRD),1 DFDG(NTGRD,NTGRD),STATEV(NSTATV),DTEMDX(NTGRD),2 TIME(2),PREDEF(1),DPRED(1),PROPS(NPROPS),COORDS(3) CCOND = PROPS(1)SPECHT = PROPS(2)CDUDT = SPECHTDU = DUDT*DTEMPU = U+DUCDO I=1, NTGRDFLUX(I) = –COND*DTEMDX(I)DFDG(I,I) = –CONDEND DOCRETURNEND1.1.13 HETVALUser subroutine to provide internal heat generation in heat transfer analysis.Product: Abaqus/StandardReferences●“Uncoupled heat transfer analysis,”Section 6.5.2 of the Abaqus Analysis User's Manual●“Fully coupled thermal-stress analysis,”Section 6.5.4 of the Abaqus Analysis User's Manual● *HEAT GENERATION●“HETVAL,”Section 4.1.9 of the Abaqus Verification Manual OverviewUser subroutine HETVAL:● can be used to define a heat flux due to internal heatgeneration in a material, for example, as might be associated with phase changes occurring during the solution;● allows for the dependence of internal heat generation on state Variables (such as the fraction of material transformed) that themselves evolve with the solution and are stored as solution-dependent state variables;● will be called at all material calculation points for which the material definition contains volumetric heat generation during heat transfer, coupled temperature-displacement, or coupled thermal-electrical analysis procedures;● can be useful if it is necessary to include a kinetic theory fora phase change associated with latent heat release (for example, in the prediction of crystallization in a polymer casting process);● can be used in conjunction with user subroutine USDFLD if itis desired to redefine any field Variables before they are passed in; and● cannot be used with user subroutine UMATHT.User subroutine interfaceSUBROUTINE HETVAL(CMNAME,TEMP,time,DTIME,STATEV,FLUX,1 PREDEF,DPRED)CINCLUDE 'ABA_PARAM.INC'CCHARACTER*80 CMNAMECDIMENSION TEMP(2),STATEV(*),PREDEF(*),TIME(2),FLUX(2),1 DPRED(*)user coding to define FLUX and update STATEVRETURNENDVariables to be definedFLUX(1)Heat flux, r (thermal energy per time per volume: JT–1L–3), at this material calculation point.。

ABAQUS子程序

ABAQUS子程序

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 材料库中没有包含的材料进行计算，扩充程序功能。

ABAQUS初学者用户子程序小例子

ABAQUS初学者用户子程序小例子

实例：最简单的杆件单轴拉伸，材料本构模型采用线弹性模型，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用户子程序

ABAQUS用户子程序ABAQUS/Standard subroutines：1。

CREEP: Define time-dependent，viscoplastic behavior （creep and swelling）.定义和时间相关的、粘塑性的运动（蠕变和膨胀)2。

DFLOW：Define nonuniform pore fluid velocity in a consolidation analysis.在压实分析中，定义非均匀孔隙流速度3。

DFLUX：Define nonuniform distributed flux in a heat transfer or mass diffusion analysis。

在热传递和质量扩散分析中，定义非均匀的分布流量4. DISP: Specify prescribed boundary conditions。

指定规定的边界条件5。

DLOAD: Specify nonuniform distributed loads。

指定非均匀的分布荷载6. FILM：Define nonuniform film coefficient and associated sink temperatures for heat transfer analysis。

对热传递分析指定非均匀的膜层散热系数和联合的散热器温度7。

FLOW：Define nonuniform seepage coefficient and associated sink pore pressure forconsolidation analysis。

对压实分析定义非均匀的渗流系数和渗入孔隙压力8. FRIC: Define frictional behavior for contact surfaces.对接触面定义摩擦9. GAPCON: Define conductance between contact surfaces or nodes in a fully coupled temperature-displacement analysis or pure heat transfer analysis。

abaqus调用子程序的方法

abaqus调用子程序的方法

Abaqus运行子程序方法试了好几种调用子程序的方法最后总算找到了最简单的了，非常感谢仿真论坛上的网友的分享：首先需要装Microsoft Visual Studio（简写mvs），再装Intel Visual Fortran（IVF）。

通常的组合是Abaqus 可使用：MVS2008+ MVS2008+ 这二种组合最好都使用英文版的。

照着附录1 方法装一般没什么问题。

但是一般都是先装abaqus后才想到去装MVS以及IVF的。

很多时候都把MVS 以及 IVF装好以后才想到去验证的，而且MVS很多都是中文的，MVS不仅装的慢卸载残留也很多好麻烦。

搞了好久都验证不了，看了附录二的方法下面2个网友的回复后豁然开朗，一切都是那么简单。

附录二的方法的主要作用可能就是将VC++和Fortran的环境变量一并启动了，参考过以前版上的安装教程，通过在ABAQUS的快捷方式中添加命令，亦能够达到同样效果，同样也是不能通过verification，但是可以运行子程序，通过在ABAQUS CAE的快捷方式的目标中，添加"D:\Program Files (x86)\Intel\Compiler\Fortran\" && （跟附录1中的修改快捷方式的目标一样）达到了同样的效果（摘自ppengine，tcboywjr网友的回复）。

只需在搜索栏中搜索，再把路径添加上去照上面的形式修改便可。

本人用的是win7系统装的是（中文版）IVF11，abaqus好早就装了，装MVS2008之前装了mvs2010，发现2010不好装ivf，然后又卸了装mvs2008中文，验证时出现附录Ⅰ的情况，就C++项通不过还以为装C++2010就好了，又重装了一下mvs2010的C++（囧），发现还是用不了，校园网速太慢又不想下mvs2008（英文版），最后在附录2的方法下面发现了两位网友的分享非常感谢。

附录二的方法也试过好像运算完了会出错误，不知道什么原因。

ABAQUS用户子程序小例子

ABAQUS用户子程序小例子

实例：最简单的杆件单轴拉伸，材料本构模型采用线弹性模型，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用户材料子程序开发及应用

ABAQUS用户材料子程序开发及应用ABAQUS用户材料子程序开发及应用摘要：本文介绍了ABAQUS用户材料子程序的开发与应用。

首先，简要介绍了ABAQUS软件及其在工程领域的广泛应用。

然后，详细阐述了用户材料子程序的概念及作用，并介绍了子程序的开发流程和必要步骤。

接着，以一个具体的材料模型开发为例，详细介绍了子程序的实现方法和注意事项。

最后，以轴对称挤压模拟为例，展示了用户材料子程序在实际工程分析中的应用，并讨论了其优点和局限性。

关键词：ABAQUS；用户材料子程序；开发；应用一、引言ABAQUS是一款广泛应用于工程领域的有限元分析软件。

其强大的建模和分析能力使得工程师可以准确地模拟和分析各种结构和材料的行为。

然而，对于一些非标准材料或特殊材料，ABAQUS自带的材料模型可能无法满足工程师的需求。

此时，用户材料子程序的开发就显得尤为重要。

二、用户材料子程序的概念及作用用户材料子程序是指由ABAQUS用户自行编写的用于描述非标准材料行为的子程序。

它可以根据特定的材料性质和应变-应力关系，定义材料模型的行为，并将其与ABAQUS的有限元分析过程相结合。

通过用户材料子程序，工程师可以更加准确地模拟和分析特殊材料的行为，提高分析结果的可靠性和准确性。

三、用户材料子程序的开发流程和步骤用户材料子程序的开发包括以下几个基本步骤：1. 确定材料模型：根据实际需要和具体材料的性质，选择合适的材料模型。

常见的模型包括线弹性模型、塑性模型、粘弹性模型等。

2. 编写用户材料子程序：使用合适的编程语言（如Fortran）编写用户材料子程序，实现材料模型的行为。

子程序应包括材料刚度矩阵计算、应力和塑性应变更新等关键计算部分。

3. 软件接口设置：将编写好的用户材料子程序与ABAQUS 软件进行接口设置，以实现子程序与有限元分析的集成。

4. 验证与调试：使用合适的测试用例对子程序进行验证和调试，确保其计算结果与实际情况吻合。

Abaqus用户子程序VUMAT

Abaqus用户子程序VUMAT
The component ordering of the symmetric and nonsymmetric tensors for the three-dimensional case using C3D8R elements is different from the ordering specified in “Three-dimensional solid element library,” Section 25.1.4 of the Abaqus Analysis User’s Manual, and the ordering used in Abaqus/Standard.
WARNING: The use of this user subroutine generally requires considerable expertise. You are cautioned that the implementation of any realistic constitutive model requires extensive development and testing. Initial testing on a single-element model with prescribed traction loading is strongly recommended.
Component 1 2 3 4 5 6 7 8 9
2-D Case 3-D Case
Initial calculations and checks
In the data check phase of the analysis Abaqus/Explicit calls user subroutine VUMAT with a set of fictitious strains and a totalTime and stepTime both equal to 0.0. This is done as a check on your constitutive relation and to calculate the equivalent initial material properties, based upon which the initial elastic wave speeds are computed.
相关主题
  1. 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
其值为1900，前面的4为1900的位数。第三项5表示是此单元编号，第四项S4R
表示单元类型，剩余四项表示此单元四个节点的编号。
CHARACTER*(*) FUNCTION DMKNAME(FNAME,DNAME,EXTEN)
C
CHARACTER*(*) FNAME,DNAME,EXTEN
C FNAME I JOBNAME
C DNAME I DIRECTORY
C EXTEN I EXTENSION
KSTEP,KINC,DTIME,TIME是ABAQUS传给用户子程序的参数。
功能：
1． 读结果文件（.FIL）中的数据。
2． 利用LSTOP开关，用户可中止ABAQUS计算
URDFIL要用到以下两个ABAQUS应用程序：POSFIL和DBFILE
文本框: Home
END DO
LTOT = LEN(DNAME)
LD = 0
DO K1 = LTOT,2,-1
IF (LD.EQ.0.AND.DNAME(K1:K1).NE.' ') LD = K1
END DO
LTOT = LEN(EXTEN)
LE = 0
为2表示文件已读到结尾后，又从头开始读起。
2．ARRAY是用来存放从结果文件中读到的那条记录的数组。要用户注意的
是ARRAY数组的大小要能存放一条记录。
3. JRCD是返回值，如果读到文件结尾，返回值为非0的整数。
功能：
读取结果文件中的一条记录。
如果想要熟练运用SUBROUTINE URDFIL,那么用户必须对结果文件（.FIL）
！打开文件
WRITE(17,*)'Opening new user external file...'
WRITE(17,*)'Writing dummy data to this file...'
END IF
RETURN
END
c
C COMPOSE A FILENAME DIRECTORY/JOBNAME.EXTEN
CALL POSFIL (NSTEP,NINC,ARRAY,JRCD)
参数：
1. NSTEP和NINC的值都是由调用它的URDFIL用户子程序接口中的参数KSTEP
和KINC传递下来。
2. ARRAY是用来存放RECORD 2000的值
3. JRCD为返回值。如果在结果文件中找到相应的STEP和INCREMENT,返回值为
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时的荷载大小
参数：
1． LSTOP是决定ABAQUS分析是否继续的开关。如果LSTOP=1，分析中止；
否则，分析继续。
2． LOVRWRT是决定能否把上个INCREMENT的结果文件覆盖的开关。
LOVRWRT＝1，覆盖，这样可大大减少结果文件的大小；否则，不覆盖。
3． LSTEP和LOVRWRT参数是留给用户自己定义的，
录是以*I开头，数据中的实数是以D开头，输出格式为E22.15 or D22.15；整
数是以I开头，接着的数字为这个整数的位数，然后才是整数的值。字符数据以
A开头。
比如记录 *I 18I 41900I 15AS4R I 3195I 3198I 3205I 3204 这条记
录第一项18表示这条记录一共有8个数据，1为8的位数；第二项为RECORD KEY
编个小程序把它们输出到外部文件中，这样对ABAQUS何时调用该用户子程序就
会有更深的了解。
下面就选出几个常用的用户子程序和应用程序进行详细解释：
一．SUBROUTINE
DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS, JLTYP,SNAME)
例1：这个例子在每个积分点施加的荷载不仅是坐标的函数，而且是随STEP
变化而变化的。
SUBROUTINE DLOAD(P,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,
1 JLTYP,SNAME)
C
INCLUDE 'ABA_PARAM.INC'
到了GETENVVAR应用程序来获得ABAQUS的环境变量，用DMKNAME子程序来合成
所需的文件名。
SUBROUTINE UEXTERNALDB(LOP,LRSTART,TIME,DTIME,KSTEP,KINC)
INCLUDE 'ABA_PARAM.INC'
CHARACTER XINDIR*255,XFNAME*80
序，但可以调用用户自己编写的FORTRAN子程序和ABAQUS应用程序。当用户编
写FORTRAN子程序时，建议子程序名以K开头，以免和ABAQUS内部程序冲突。
2．当用户在用户子程序中利用OPEN打开外部文件时，要注意以下两点：一
是设备号的选择是有限制的，只能取15－18和大于100的设备号，其余的都已
5． COORDS为当前积分点的坐标；
6． 除F外，所有参数的值都是ABAQUS传到用户子程序中的。
功能：
1． 荷载可以被定义为积分点坐标、时间、单元编号和单元节点编号的函数。
文本框: Home
2． 用户可以从其他程序的结果文件中进行相关操作来定义积分点F的大小。
结果文件存放的是一条条的记录。为了区别不同的记录集，记录中有一项为
RECORD KEY。不同的记录集，RECORD KEY是不同的。比如节点坐标的记录集
文本框: Home
（COORD）的RECORD KEY 为8；主应力（EP）的RECORD KEY为403；用户可根
据不同的RECORD KEY来分辨不同的记录集。每条记录包括好多个数据。每个记
C
DIMENSION TIME(2),COORDS(3)
CHARACTER*80 SNAME
PARAMETER (PLOAD=100.E4)
C
IF (KSTEP.EQ.1) THEN ！当STEP=1时的荷载大小
P=PLOAD
ELSE IF (KSTEP.EQ.2) THEN ！当STEP＝2时的荷载大小
活性。
功能：
1． 可以用来和其它用户子程序及其它软件进行数据通讯。
2． 可以用来在适当的时间打开，关闭外部文件。
3． 用户可以把自己编写的ABAQUS扩充功能的程序代码通过此用户子程序嵌
文本框: Home
入到ABAQUS中。
例2．新建一个与JOB名相同但扩展名(.ALE)不同的文件，此用户子程序用
的格式有所了解。
要让ABAQUS输出结果文件，用户必须在INPUT文件或KEYWORD中加入*EL
FILE或着*EL NODE。有些数据是ABAQUS自动写道结果文件中的，如ABAQUS的
版本号，节点编号和节点坐标；有些是用户用*EL FILE或着*EL NODE定义的，
比如应力、应变、节点位移等等。
被ABAQUS占用。二是用户需提供外部文件的绝对路径而不是相对路径。
3．ABAQUS 应用程序必须由用户子程序调用。
当用到某个用户子程序时，用户所关心的主要有两方面：一是ABAQUS提供
的用户子程序的接口参数。有些参数是ABAQUS传到用户子程序中的，例如
SUBROUTINE DLOAD中的KSTEP，KINC，COORDS；有些是需要用户自己定义的，例
P=COORDS(1)**3*PLOAD
ELSE IF (KSTEP.EQ.5) THEN ！当STEP=5时的荷载大小
P=COORDS(1)**4*PLOAD
END IF
RETURN
END
二． UBROUTINE UEXTERNALDB(LOP,LRESTART,TIME,DTIME,KSTEP,KINC)
IF(LOP.EQ.0) THEN
FNAMEX=DMKNAME(XFNAME(1:LXFNAME),XINDIR(1:LXINDIR), '.ALE')
！生成所要新建文件的文件名
OPEN(UNIT=17,FILE=FNAMEX,STATUS='UNKNOWN',FORM='FORMATTED')
CHARACTER DMKNAME*255,FNAMEX*80
C
LXFNAME=0
LXINDIR=0
XFNAME =' '
XINDIR =' '
CALL GETENVVAR('FNAME',XFNAME,LXFNAME) !读取input文件名
CALL GETENVVAR ('OUTDIR',XINDIR,LXINDIR) !读取input文件所在的路径
参数：
1． F为用户定义的是每个积分点所作用的荷载的大小；
2． KSTEP,KINC为ABAQUS传到用户子程序当前的STEP和INCREMENT值；
3． TIME(1),TIME(2)为当前STEP TIME和INCREMENT TIME的值；
4． NOEL,NPT为积分点所在单元的编号和积分点的编号；
如F。二是ABAQUS何时调用该用户子程序，对于不同的用户子程序ABAQUS调用
相关文档
最新文档