JXNC 发表于 2005-12-21 15:06:04

自定义材料重启动分析出错

我基于usrdll里的ovl30u_ve1.f编译了一个2D材料动力本构模型。 原程序里有很多关于温度应力应变的代码,且基本上是基于IF (MG.GT.0.AND.MG.LE.15) THEN ...和 IF (MK.GT.0.AND.MK.LE.15) THEN...来建立本构关系的,在CTI(1)-CTI(8)里要输入与剪切模量和体积模量以及与温度相关的参数。而我的本构不需要考虑温度的影响,希望温度应力应变都为零,因此我把参数CTI(1)-CTI(8)都定义为零,CTI(8)以后为用户材料参数。另外,原程序用到ARRAY()历史变量,当MG和MK都小于等于5时,数组ARRAY()的元素最大数目为70,若MG和MK大于5,数目还要增加。
   
    我的程序是假定MG和MK数目都等于零,这样利用空出的一部分ARRAY()的变量来存放用户自定义材料的变量(从ARRAY(9)开始,ARRAY(1)—ARRAY(8)不改变)。定义温度荷载为零,在Temperature-Dependent Properties里定义两个温度,即0和1000,将该自定义材料用于ADINA重启动瞬态动力计算时(静力计算用ADINA软件自带的各向同性弹性模型),出现如下错误:

*** E R R O R ***
IF TEMPERATURE, TEMPERATURE GRADIEND OR
PIPE INTERNAL PRESSURE LOADINGS WERE PRESENT
IN THE PREVIOUS ANALYSIS, THEY MUST BE PRESENT
IN THE RESTART ANALYSIS AND VICE VERSA
*** ERROR NO. A1014 ***                                                      
Temperature loading used in a previous run but not in restart

shaoyu1979 发表于 2005-12-21 23:42:14

Re:自定义材料重启动分析出错

看错误的提示:温度载荷存在于先前的载荷步中,但是重启动没有.
也许你应该加温度载荷试试.

JXNC 发表于 2006-1-2 21:45:48

Re:自定义材料重启动分析出错

To shaoyu,在Adina8.1/usrdll目录里的很多程序里,都有这样一段话——“NOTE THAT THE VARIABLES PASSED TO THE SUBROUTINE WHEN KEY=3,4 ARE THESE CALCULATED LAST: I.E., CALCULATED IN THE LAST SUBDIVISION WHEN KEY=2. HENCE THESE VARIABLES ARE NOT CALCULATED WHEN KEY=3,4.”
> 从这句话的意思来推测,KEY=3,4时的变量值是KEY=2时最后一个子步(即KTR=INTER时)的计算结果,因此在KEY=3计算STRESS/STRAIN MATRIX时根据最后一个子步的结果计算一次,也就是说KEY=2的程序段由于有INTER个子步,在每一时间步(DT)/荷载步里要被上部程序调用INTER次,而KEY=3时只计算一次STRESS/STRAIN MATRIX。但在KEY=3的程序段里计算STRESS/STRAIN MATRIX时,程序里出现了DDT-----一个子步的时间增量,如OVL30U_VP1里,从这里看似乎KEY=3的程序段要被调用/更新INTER次,所以我搞不清到底在一个时间步(DT)里,在KEY=3的程序段里应力应变矩阵要计算一次还是INTER次,是根据每一个子步结束时的变量计算一次还是根据最后一个子步只计算一次?

shaoyu1979 发表于 2006-1-2 23:41:46

Re:自定义材料重启动分析出错

你仔细看一下帮助文挡就可以发现并没有在每个DDT里计算stress/strain matrix,而是计算一个stress/strain matrix.DDT的作用是在计算CP = CP + 1.5*DDT*GAMA*(EST/YLD - 1.)/EST,对应帮助文挡389页,3.14-22.
粘塑性这个子程序是帮助文挡里写的最详细的一个了.子程序里的好多公式都可以在里面找到对应.

JXNC 发表于 2006-1-12 23:20:50

Re:自定义材料重启动分析出错

谢版主了。我还想问下adina在进行地震动力分析时一般是采用什么方法,是不是用增量迭代法?
就等肯-张模型二次开发来说,采用切线模量和割线模量有什么区别?是不是采用割线模量时用总应变来计算模量,然后按这个模量来计算各子步的应力赠量?
页: [1]
查看完整版本: 自定义材料重启动分析出错