ABAQUS-二次开发资料-UMAT

1、各个楼层及内容索引2-什么是umat3-umat功能简介4-umat开始的变量声明5-umat中各个变量的详细解释6-关于沙漏和横向剪切刚度7-umat流程和参数表格实例展示8-fortran语言中的接口程序interface9-关于umat是否可以用fortran90编写的问题10-17-fortran77的一些有用的知识简介20-2530-32-弹塑性力学相关知识简介34-37-用户材料子程序实例john-cook模型压缩包下载38-john-cook模型本构简介图40-用户材料子程序实例john-cook模型完整程序+david详细注解欢迎大家来看看,并提供意见,完全是自己的diy的,不

2、保证完全正确,希望共同探讨,以便更正,带?部分,还望各位大师同仁指教1 什么是umat?1.1 umat功能简介!-摘自庄茁老师的书umat子程序具有强大的功能，使用umat子程序：(1)可以定义材料的本构关系，使用abaqus材料库中没有包含的材料进行计算，扩充程序功能。abaqus软件2003年度用户年会论文集(2)几乎可以用于力学行为分析的任何分析过程，几乎可以把用户材料属性赋予abaqus中的任何单元；(3)必须在umat中提供材料本构模型的雅可比（jacobian）矩阵，即应力增量对应变增量的变化率。(4)可以和用户子程序“usdfld”联合使用，通过“usdfld”重新定义单元每一

3、物质点上传递到umat中场变量的数值。1.2 umat开始的变量声明由于主程序与umat之间存在数据传递，甚至共用一些变量，因此必须遵守有关书写格式，umat中常用的变量在文件开头予以定义，通常格式为：subroutine umat(stress,statev,ddsdde,sse,spd,scd,1 rpl,ddsddt,drplde,drpldt,2stran,dstran,time,dtime,temp,dtemp,predef,dpred,cmname3 ndi,nshr,ntens,nstatv,props,nprops,coords,drot,pnewdt,4 celent,dfg

4、rd0,dfgrd1,noel,npt,layer,kspt,kstep,kinc)cincludeaba_param.inc-此处是将abaqus本身自带的参量精度定义的文件包含进来后面详说ccharacter*80 cmnamedimension stress(ntens),statev(nstatv),1 ddsdde(ntens,ntens),ddsddt(ntens),drplde(ntens),2 stran(ntens),dstran(ntens),time(2),predef(1),dpred(1),3 props(nprops),coords(3),drot(3,3),dfg

5、rd0(3,3),dfgrd1(3,3)-以上是变量声明吧!user coding to define ddsdde,stress,statev,sse,spd,scdand,if necessary,rpl,ddsddt,drplde,drpldt,pnewdt-此处,看来是将用户定义材料属性的fortran程序编入return-这是返回值end-结束umat中各个变量的详细解释凡是-注明david的,都是我自己猜的,仅供参考ddsdde (ntens ,ntens)是一个ntensnumber of the tensions-david维的方阵，称作雅可比矩阵，应力增量/应变增量的偏导数，

6、ddsdde (i ,j)表示增量步结束时第j个应变分量的改变引起的第i个应力增量的变化!雅可比是一个对称矩阵，除非在“*user material”语句中加unsymm参数stress (ntens)应力张量矩阵，对应ndinumber of the direct components-david个直接分量和nshrnumber of the shear components-david个剪切分量.在增量步的开始,应力张量矩阵中的数值通过umat和主程序之间的接口传递到umat中,在增量步的结束,umat将对应力张量矩阵更新,即return.对于包含刚体转动的有限应变问题，一个增量步条用um

7、at之前就已经对应力张量进行了刚体转动,因此在umat中只需处理应力张量的共旋部分-这部分我没看明白,敬请高手指点.umat中应力张量的度量为柯西（真实）应力statev (nstatev)state variables (number of the state variables) 用于存储状态变量的矩阵，在增量步开始时将数值传递到umat中。也可在子程序usdfld或uexpan中先更新数据，然后增量步开始时将更新后的数据传递到umat中。在增量步的结束必须更新状态变量矩阵中的数据。 和应力张量矩阵不同的是：对于有限应变问题，除了材料本构行为引起的数据更新以外，状态变量矩阵中的任何矢量或者

8、张量都必须通过旋转来考虑材料的刚体运动。-不懂,请高手指教 状态变量矩阵的维数natatev，等于关键字“*depvar”定义的数值。状态变量矩阵的维数通过abaqus输入文件中的关键字“*depvar”定义，关键字下面数据行的数值即为状态变量矩阵的维数。 材料常数的个数，等于关键字“*user material”中“constants”常数设定的值。props (nprops)材料常数矩阵，矩阵中元素的数值对应于关键字“*user material”下面的数据行。sse，spd，scd分别定义每一增量步的弹性应变能elastic strain energy，塑性耗散plastic dissi

9、pation和蠕变耗散creep dissipation。它们对计算结果没有影响，仅仅作为能量输出。stran (ntens)：应变矩阵；dstran (ntens)：d-大抵代表deta,增量的意思-david应变增量矩阵；dtime：增量步的时间增量；ndi：直接应力分量的个数；nshr：剪切应力分量的个数；ntens：总应力分量的个数，ntens =ndi +nshr。1.3关于沙漏刚度控制和横向剪切刚度使用umat时需要注意单元的沙漏控制刚度和横向剪切刚度。通常减缩积分单元的沙漏控制刚度和板、壳、梁单元的横向剪切刚度是通过材料属性中的弹性性质定义的。这些刚度基于材料初始剪切模量的值，通

10、常在材料定义中通过“*elastic”选项定义。但是使用umat的时候，abaqus对程序输入文件进行预处理的时候得不到剪切模量的数值。所以这时候用户必须使用“*hourglass stiffness”选项来定义具有沙漏模式的单元的沙漏控制刚度，使用“*transverse shear stiffness”选项来定义板、壳、梁单元的横向剪切刚度。-一个问题:得不到剪切模量的数值和解决方案:用户必须使用“*hourglass stiffness”选项来定义具有沙漏模式的单元的沙漏控制刚度，使用“*transverse shear stiffness”选项来定义板、壳、梁单元的横向剪切刚度1.4关

11、于umat的流程图和参数表格实例跟大家说说所谓的接口程序interface-fortran的知识在fortan语言中,主调程序和被调程序是分别编译的.由于fortran90对过程的许多功能做了扩充,有些功能单靠简单的调用语句已经无法反应,因而系统也就无法进行正确的编译,这时需要在主调程序中加入interface接口块,通过它为主调程序和被调程序指明一个显示的接口.如果被调用中哑元含有假定形状assumed-shape数组,或可选变元,或含键盘输入的参数,就需要interface接口块说明.一般来讲,在fortran90程序之间需要提供interface块有三种方法:1.将interface接口

12、块直接写入调用程序,并复制被调用程序的参数列表这种方法简单易用,但也相应增加了维护代码的工作量,因为只要被调用程序的参数列表发生变化,就必须相应改变interface接口块和调用call语句.2.可以将一个软件包中所有程序的interface块写入一个模块中,该模块被软件包中的所有程序使用.这样做的优点是只需一个模块来检查interface定义,缺点是仍需对此模块和调用语句进行维护.3.fortran90编译器可在contains语句后自动提供过程之间的interface块,这种interface块可用于使用模块的任何程序.建议在同一个软件包中使用23的形式,在调用软件包的入口程序时采用12的

13、形式!是不是在umat中,我们所编译的带接口的fortran程序为调用程序,原abaqus主程序为被调用程序,调用程序中的第一部分我们先复制被调入程序的参数列表?似乎和abaqus主程序调用umat有些相反了?不过个人认为interface作为一个接口块,在fortran语法中应该放在主调程序中,且复制被调程序的参数列表.而umat的参数变量的声明,只不过是为了和abaqus使用一致的变量格式,当fortran程序处理完材料的本构定义之后再以这一致的形式将变量值返回到abaqus主程序中吧2 关于umat用fortran90可不可以的问题在论坛上搜索了关于这个问题,某位高手说是可以的,只要你自

14、己装的fortran编译器能成功编译你的fortran90程序即可,个人认为也是如此,不过我还没有证明过!2.1关于fortran77的一些有用的简介1:fortran77的基本结构:1. 一个fortran源程序由一个或者多个程序单位组成,每个独立的程序单位以end语句结束2. 每个程序单位包括若干行不能一行写多条语句,但是可以一个语句写成行;分为语句行执行语句行和非执行语句和非语句行注释行源程序的书写格式:1. 不区分大小写:每行只能80列以内,并把80列分为四个区2. 1-5列:标号区1-5位整数;第一列为*或者c时,为注释行3. 第6列,续行标志区非空格或非0;最多19个续行4. 7-

15、72列:语句区书写语句:一行只能写一个语句5. 73-80列:语句注释区一般做编号注释程序的编辑与运行:1. 创建源程序文件并编写源程序2. 编译并连接源文件3. 运行程序编译生成的可执行文件常量:六种1. 整型常量integer4位:正负整数和02. 实型常量real4位:小数和指数形式3. 双精度常量double precision8位4. 复型常量complex8位5. 逻辑性常量logical4位6. 字符型常量character1位变量 变量名1. 第一个字符是字母第一个字符后可跟1-5个字母或者数字2. 不区分大小写3. 空格不起作用4. 允许变量名与语言中具有特定含义的字同名,但

16、尽量不适用5. 尽量见名知义 变量类型 不同的变量类型用来存放不同类型的常量数据.变量相应的也有六种;变量在使用前必须首先对其进行类型说明,三种说明方法: 按优先级别排列 1. 类型说明语句:类型 变量名列表,多个变量名之间用逗号隔开,如 real a,b,c double precision x,y,z或者real*8 x,y,z character*5 缺省字符长度5 str1,str2*8,str3*19 此处,str1未指明长度,则默认使用缺省长度5;str2的长度为 8;str3的长度为19 2. implcit语句:implicit 类型字母表,类型字母表,. 如: implici

17、treal(a-d),interger(i-m),double precision(x,z) 3. i-n规则:fortran规定,不加说明的情况下,i-n为整型,其他都为实型几点说明类型语句说明位于所有可执行语句的前面;implicit语句位于最前面;implicit none取消implicit说明和i-n规则,所有的变量必须显式说明;只在本程序单位有效.2.2关于fortran77的一些有用的简介2:算术表达式:1. 运算符: +,-,*,/, *2. 优先级: ( ),* ,* / ,+-3. 书写问题1. * 不能省略2. 括号不分大小写,成对出现3. 多次乘方,按先右边后左边处理4

18、. 运算符不能连续出现,要用小括号隔开5. 运算顺序)-函数-*-*,/-+,-6. 运算中的类型问题:运算类型相同:结果仍为原类型;不同,则自动转换成同一类型7. 误差问题:1. 溢出:超出有效数字范围-解决:很大或者很小的数用实型的指数表示2. 误差:由于有效数字的位数限制,实型数运算存在误差-解决:转换成双精度型或者避免因为书写不当造成有效数字的丢失简单输出输入语句:输入输出三要素:对象哪些数据;格式;设备.输出语句输出语句的分类:格式输出表控格式输出系统标准格式;无格式输出二进制1. 表控输出语句:按计算机系统规定的格式输出:系统自动为每种类型的数据规定了列数1. 整数的表控输出格式与

19、具体的计算机系统有关:规定每个整数占13列,打印在右端,左补空格;2. 实数的表控输出格式:规定每个实数占17列,打印在右端,左补空格,小数部分占7列;当实数的绝对值=10*7或1时,按标准的指数格式输出,共占15列,指数4列,小数6列2. 表控格式输出语句:1. print*,: print*,56.8,1252. write(*,*)输入语句 输入输出三要素:对象;格式;设备. 分类:同上1. 表控输入语句1. 自由格式输入-语句:read*,;read(*,*)2. 输入数据以逗号或者空格作为间隔3. 变量名称为输入表4. 输入的数据应和输入表的变量个数类型次序严格地一一对应;少了,程序

20、停止,等待继续输入;多了,程序继续进行,多余的不起作用;较多的数据可以几个一组,回车,再输入几个一组,回车.5. 重复数据,可以7*3-7,7,76. 每一个read(*,*)和write(*,*)语句从一个新的记录以回车结束的一批输入输出数据开始读数输出1. 例如:read(*,*) a,b,c2. read(*,*) d,i,j3. 输入: 2.3,-63.5回车4. 6.4,91.0回车5. 5,8回车6. 结果: a=2.3,b=-63.5,c=6.4,7. 从新记录开始读数8. d=5.0,i=8,j未被赋值parameter语句作用:将程序中经常用到的参数或字符串定义成一个符号常量

21、,其值不可改变.语句:parameter(p1=c1,p2=c2,.,pn=cn)注意:1. 符号常量的命名规则与变量名相同,但在程序中其值不可改变,也不能赋值;2. 符号变量也有类型,可用前面的三种类型说明方法说明类型;3. 参数语句是非执行语句,位于所有可执行语句的前面,单位与类型说明语句的后面;4. 一条语句可以定义多个符号常量;5. 优点:方便修改程序end,stop,pause语句end语句:结束标志,有且仅有一条pausen语句:暂定执行;用于调试程序,n可以是一个字符串或不超过5位的数stopn语句:停止运行语句;用于调试程序,n可以是一个字符串或不超过5位的数2.3关于fort

22、ran77的一些有用的简介3:逻辑运算和选择结构 关系表达式1. 构成选择判断的基本式子2. 关系运算符:1. .gt.greater than 2. .ge.greater than or equal to =3. .lt.limiter than 4. .le.limiter than or equal to =5. .eq.equal to =6. .ne.not equal to 3. 一般形式:4. 运算结果:逻辑值:真.true.假.false.5. 运算顺序:算术运算关系运算 逻辑表达式1. 运算符:1. .and.2. .or.3. .not.4. .eqv.逻辑等5. .ne

23、qv.逻辑不等2. 一般形式:3. 结果:逻辑值:真.true.假.false.4. 运算顺序:算术运算-关系运算-逻辑运算5. 逻辑运算优先级:.not.-.and.-.or.-.eqv.-.neqv.关于fortran77的一些有用的简介4:if类选择结构 用块if实现选择结构:三种典型形式1. 基本形式1. if(条件) then (块if语句)2. 块1 (then块)3. else (else语句)4. 块2 (else块)5. endif (endif语句)6. 说明:if.then语句为块if结构的入口语句;endif语句为出口语句,必须一一对应,配对使用2. 简单结构1. if

24、(条件) then2. 块3. endif4. 说明:没有else块3. 嵌套结构1. if()then2. 块13. else if() then4. 块25. .6. else if()then7. 块n8. else块n+19. endif 逻辑if语句只用一行表示一种选择结构,当且仅当条件成立时执行,并且只执行一条语句;if(条件)语句 算术if语句if n1,n2,n3当算术表达式的值 0执行标号为n1的语句;关于fortran77的一些有用的简介5:循环结构 结构形式:循环体由一些可执行的语句组成+循环控制语句控制循环的开始和结束 分类:条件型循环和计数型循环do循环goto语句实

25、现循环 一般形式:goto其中:s1为语句标号 功能:程序执行到此语句时,无条件的转向标号为s1的语句do语句实现循环 当循环的初值终值循环次数都已知时,可用; 组成:一个do语句和循环体组成 一般形式: dos1 i=e1,e2 ,e3do i=e1,e2 ,e3 .s1 .enddo例如do10i=1,19,2do10i=e1,19 ,2 sum=sum+1s1continue10 sum=sum+1do i=1,19,2 sum=sum+1enddo 说明1. i为循环变量,s1为语句标号,是本程序单位中另一可执行语句的标号;2. 步长可以省略,缺省值=1;3. 循环初值e1,终值e2和

26、步长e3都可以是常量变量表达式;4. 由于实数在内存中存储的误差,i,e1,e2,e3尽可能用振型量5. e1,e2,e3都可正可负,e1,e2,可为0,但是e3不能为0. 具体执行过程1. 执行do语句,首先计算表达式e1,e2,e3的值,若他们的类型与循环变量i不一致,则自动转换成循环变量的类型2. 将e1的值赋予循环变量i,及执行赋值语句:i=e1;3. 计算循环次数:r=max0(e2-e1+e3)/e3,max0表示从多个整型变量中取最大的一个;4. 检查循环次数:若r=0则不执行循环体内的语句,跳出循环;r0 则执行循环体内的语句5. 执行循环终端语句:i=i+e3,即是循环变量获

27、得一个新值,而循环次数r自动减1;6. 返回步骤4,继续执行,直到r=0. continue语句 循环终端语句必须是可执行语句;那么,这种作为循环终端的语句具有双重作用:一是作为循环终端的标志;而是要完成自身的功能.因此影响了程序的可读性.fortran用一个专门的语句作为do循环的终端语句,即continue语句.它自身没有任何功能. 一些规定1. 循环变量在循环体内只能被引用,不能被赋值;2. 在执行do循环体期间,e1,e2,e3的值不能被改变,因为他们决定了循环的次数3. 离开do循环后,循环变量可以在循环体外被引用,它的值为脱离循环时最后一次被赋的值;4. 程序中用到转移语句,规定:

28、只允许从循环体内-体外;反之不行;5. 循环终端语句必须是除goto,块if,endif,end和stop语句外的任何可执行语句 do循环的嵌套 在一个do循环中还可以包含一个或者多个完整的do循环,这就是do循环的嵌套. 一般形式: do 10 i=1,10 . . . do20 j=1,10 . . . 20 continue . . . 10 continue 说明: 嵌套要完整,不能交叉 循环变量的名字,规定:并列的循环:循环变量的名字可以相同;嵌套的循环:循环变量的名字不可以相同 若多层循环的结束语句在同一个地方,可以共用一条continue语句 控制转向语句的使用体内-体外当型循环

29、的实现 在无法确定循环次数的情况下可以使用当型循环.当型循环是指执行循环体要依据实现给定的条件:当条件成立时执行循环,否则不执行. 用do-while语句实现当型循环 一般形式: do s1 , while(条件) . s1 用块if和goto语句实现循环 一般形式: s1if(条件) then 块 goto s1 endif直到型循环的实现所谓直到型循环,是指先执行循环体,再判断条件.如果条件为假,继续执行循环体,直到条件为真时终止循环. 用逻辑if语句实现: s1循环体 if(条件) gotos1几种循环形式的关系和比较 do循环适用于已知循环次数的情况 几种循环可以互换 do循环:条件型

30、循环可用次数作为条件 当型循环:直到型循环 当型:块if语句(单边)+goto语句(先判断后执行) 直到型:逻辑if语句+goto语句(先执行后判断) 各种循环可以相互嵌套2.4关于fortran77的一些有用的简介6:数据的输入和输出数据输入输出需要确定的三个基本要素: 输入输出的设备 输入输出的格式 输入输出的数据 系统中隐含的 输入输出的设备为:键盘显示器和打印机说明:#(设备,格式)数据列表,当设备显示为*,为默认设备输出,好像是显示器或默认设备输入,键盘吧;格式为*,默认格式输出输入-david有格式的输出输出语句的一般形式: write (*,s1) s1 format(格式说明)

31、或者 print s1, s1 format(格式说明)格式说明符: 主要介绍:i,f,e,d,g,l,a,(撇号),h,x,r(重复系数),/(斜杠)i 编辑符(integer)作用:用于整型数据的输出.一般形式:iw或iw.m其中:i表示整型输出,w为字段宽度,m表示输出数据的最少数字位数注意:数据输出时,在指定的区域内向右靠齐;如果数据的实际位数大于指定的字段宽度w,则不输出数据,而在该区域内充满*号;当m大于数据的实际位数时,前面添0,若小于数据实际位数,则不起作用f 编辑符(fixed point number)作用:用于实数的小数形式输出,一般形式:fw.d其中:f表示实数的小数形

32、式输出;w为字段宽度;d为输出数据的小数位数e 编辑符(iexponent)作用:用于实数的指数形式输出,一般形式:ew.d其中:e表示实数的指数形式输出;w为字段宽度;d为数字部分的小数位数注意:指数部分占4列,负号占1列,小数点前为0.如123.45-0.12345e+03d 编辑符(double precision)作用:用于双精度的指数形式输出,用法和e 编辑符相仿.一般形式w.dg 编辑符作用:由系统根据实际数据的大小来决定使用f编辑符还是e编辑符.一般形式:gw.dl 编辑符作用:用于逻辑型数据的输出,一般形式w其中表示整型输出,w为字段宽度a 编辑符作用:用于字符型数据的输出,一

33、般形式:aw或a其中:a表示整型输出,w为字段宽度;若不指定,则表示按实际长度输出 (撇号) 编辑符作用:用于输出字符常量,即把撇号内的字符串原样输出.注意:如果输出的字符中包含撇号,则用两个连续的撇号代表一个要输出的撇号h 编辑符作用:用于输出字符常量.一般形式:nh其中:h表示输出字符常量;n为输出字符个数;str为输出的字符串 (较少使用)x 编辑符作用:用于输出空格.一般形式:nx其中:x表示输出空格;n表示输出的空格数重复系数r在format语句中,如果出现几个(或者几组)相同的格式编辑符,则可以利用重复系数而只写一个(或者一组)编辑符.如format(a=,/,4(5(1x,f4,

34、0),/)反斜杠/编辑符作用:结束本记录的输出,开始下一个记录的输出,通常指换行.write语句和format语句的相互作用write语句的输出变量个数与format语句的编辑符(不含撇号,h和x)个数可以相等,也可以不等;如果编辑符个数多,则剩余的编辑符不起作用;如果变量的个数多,则当编辑符用完后,重新使用该格式说明,当如果格式说明含带重复系数的编辑符组,则格式说明用完后,只有最右面一个带重复系数的编辑符组(包含重复系数)及其右面的编辑符被重复使用.可以有空格式说明,如format(),用于输出一个空行.有格式的输入有格式的输入语句 一般形式: read(*,s1) s1 format(格式

35、说明由各种格式编辑符构成) 例如: read(*,100) a,b,c 100 format(f5.1,e12.2,f7.2) end 键盘输入:_15.7_2345.67e+04_705.83enter在printwriteread语句中包含格式说明 例如: print 100,k,y 100 format(18,f7.2) 也可以写成:print(18,f7.2),k,y 注意写法: (.)关于fortran77的一些有用的简介7:数组使用原则:先声明,后使用说明方法: 用类型说明语句(显式说明) 用dimension语句(隐式说明)数组的说明和数组元素的引用 用类型语句说明数组1. 一般

36、形式:类型说明数组说明符2. 其中:数组说明符的一般形式为:数组名(维数说明符,.)维数说明符,由下标下界:下标上界组成3. 例如:realx(1:10),w(1:2,1:3),k(10:20)或者integer b(1:100),py(0:2,0:3,0:5) 用dimension语句说明数组1. 一般形式:dimension 数组说明符,. 说明:1. 在数组说明符中,维数说明符(下标)的个数成为数组的维数2. 维数说明符只能使用整型常量或者整型符号常量表达式:如 paraneter(i=1,j=10) realkx(i:j+5)3. 维数说明符的下标下界为1时,可以省略.如:real x

37、(1:10)-real x(10)4. 数组说明语句必须写在所有可执行的语句之前.属于非执行语句 数组元素的引用1. 一般形式:数组名(下标,.)2. 即要有确定的数组名和下标值,如xn(5),w(1,3),kw(1,2,3)3. 引用数组元素时,下标可用算术表达式,如果算术表达式的值为实行,则自动取整.数组的逻辑结构和存储结构逻辑结构:数组所表示的实际数据结构 存储结构:数组在机器内存储时的排列结构 一维数组 逻辑结构:依次排列的一串数据 存储结构:一组连续存放的一列数据块 二维数组 逻辑结构:一张二维数据表 存储结构:一组按列连续存放的数据块 三维数组 逻辑结构:若干张二维数据表 存储结构

38、:一组按页连续存放的数据块数组的输入和输出三种方式:用do循环用隐含do循环用数组名 用do循环实现数组的输入输出1. 优点:数组元素的输入输出次序可由用户控制2. 缺点:做一次循环就换行输入或输出 用隐含do循环实现数组的输入输出1. 优点:既能控制数组元素的输入输出顺序,又能控制一行内输入输出数据的个数2. 例如:read(*,*) (g(i,j),j=1,3),i=1,2),由于是一个read语句,所以既可以一行输入,也可以多行输入,键盘输入如下:86,75,72enter87,70,83enter3. 注意:一个read语句可以多行输入;但是多个read语句时,每一个read语句必须从

39、心的一行读数.用数据名进行数组的输入输出使用时,其顺序要与数组元素在机器内的存储顺序一致 例如: dimension k(5) read *,k-对数组进行整体操纵 等价于:read*,k(1),k(2),k(3),k(4),k(5) 也等价于:read*,(k(i),i=1,5)使用data语句给数组赋初值 一般形式: data变量列表初值表,变量列表初值表. 功能:在程序编译期间给变量或者数组赋初值.其中,变量列表可以是变量名数组名数组元素隐do循法;初值表只能是常量,不允许出现任何形式的表达式 例如:data a, b/7.85,9.1/代表赋初值a=7.85,b=9.1-david,

40、i,j, k /5,10,15/代表赋初值i=5,j=10,k=15-david, 例如: dimensionk(2,3) data(k(i,j),j=1,3),i=1,2)/90,23,20,42,14,32/-初值列表2维3列,i=1: 90 23 20 i=2: 42 14 32 或 datak/90,42,23,14,20,32/-排列为按列排,排满一列之后,再排下一列; 90 23 20 -david 42 14 32 例如: dimensiona(10) data a/10*1.0/(表示10个1.0)注意:data语句属于说明语句,但是它可以放在end语句之前的任意行;当程序中有

41、多个data语句给同一个变量赋初值时,以最后一条为准;程序在编译期间给变量赋予初值,在程序执行期间,data语句不起任何作用!2.5关于fortran77的一些有用的简介8子程序fortran子程序:包括函数子程序,子例行程序,数据块子程序执行:从主程序开始执行,遇到调用语句再执行相应的子程序.不同类型的子程序,关键字不同,调用方法也不同 函数子程序:一种可以作为函数来调用的子程序(外部函数)1. 定义:一般形式,由function语句和子程序体组成类型说明function函数名(虚参表) .(子程序体) endfunction函数名(虚参表)类型说明 函数名 .(子程序体)end2.1. f

42、unction语句:是函数子程序的第一条语句,标志着该函数子程序的开始1. 类型说明function函数名(虚参表)2. 注意:1. 虚元也有类型,需在子程序体中说明1. 例如:realfunction intep(x1,x2,x3)2. integerx1,x2,x32. 函数名的命名规则和类型都和变量相同3. 虚参可以是变量名数组名子程序名,但不允许用常量和数组元素,它表示了函数自变量的个数顺序和类型.2. 子程序体:完成一个具体任务的程序段3. 注意:1. 若无虚参时,括号不能省2. 函数子程序中所有变量和标号(除函数名和虚参外),与其他程序单位无任何关系3. 函数体的说明部分包括对虚参

43、和本函数体内所用变量和数组的说明4. 函数体中可设置一条或者多条return语句,表示执行到此语句时返回调用程序.1. 当return语句和end语句紧挨着的时候,可省略return语句2. 也可以不设return语句,但需从中间返回时,必须设置return语句5. 函数名的作用:函数名在函数体中一定要被赋值,因为函数名把函数值带回调用程序. 函数子程序的调用1. 一般形式:调用方式和内部函数相似:函数名(实参数) 或 函数名( )2. 说明:1. 调用程序中函数名必须与函数子程序中定义的函数名相同2. 实参与虚参在个数类型位置上必须一一对应,但名字可以不同3. 当虚参是变量名的时候,实参可以

44、是常量变量数组元素或者表达式;但是当虚参要在函数体中被赋予初值的时候,则实参不可以是常量或者表达式因为两者共用一个存储单元4. 函数子程序是一个独立的程序单位,可以作为一个单独的源程序进行存贮和编译,并与调用程序连编后才能起作用.4 材料本构的相关力学知识虽然目前要做的工作就是用umat把一个新的材料的本构加进abaqus中调用,但是我并不想仅仅局限于此.在这个过程中,我会把自己遇到的所有相关问题学懂弄透再去编这个程序.所以估计我这个帖子会很长滴!我力学知识不是很好,所以这两天要好好学习一下力学知识,主要是弹塑性力学吧!说说弹塑性力学-1弹性力学塑性力学弹塑性力学 弹性力学和塑性力学时固体力学的两个重要分支.1. 固体力学:研究固体材料及其构成的物体结构在外部干扰(载荷温度变化等)下的力学响应的科学.按不同的研究对象区分为不同的学科分支.2. 弹性力学:研究固体材料及由其构成的物体结构在弹性变形阶段的力学行为,