| 查看: 2657 | 回复: 5 | ||
xiaomiao126金虫 (小有名气)
|
[求助]
关于ansys二次开发UPFs，修改usercreep子程序的问题，求好人帮忙 已有1人参与
|
|
大神求教啊！！！！！本人最近在做关于ansys UPFs的二次开发！想把自己的蠕变方程写入ansys中！网上查资料知道修改 ansys 的usercreep程序可以实现此功能！ ansys 提供的程序可以实现 这个蠕变方程的编译 delcr := c1 * seqv ^ n * creqv ^ m * exp (-b/T) * dtime c1，n，m，b都是常熟！delcr蠕变增量，creqv蠕变应变 他的程序如下: #include "impcom.inc" DOUBLE PRECISION ZERO PARAMETER (ZERO = 0.0d0) c c INTEGER ldstep, isubst, matId , elemId, & kDInPt, kLayer, kSecPt, nstatv, & impflg, nprop DOUBLE PRECISION dtime , time , temp , dtemp , toffst, & creqv , seqv , pres DOUBLE PRECISION prop(*), dcrda(*), Ustatev(nstatv) c c DOUBLE PRECISION c1 , c2 , c3 , c4 , & con1 , delcr , t if (seqv.LE.ZERO.AND.creqv.LE.ZERO) GO TO 990 t = temp + toffst c1 = prop(1) c2 = prop(2) c3 = prop(3) c4 = prop(4) con1 = ZERO if(c4.ne.ZERO .and. t.gt.ZERO) con1 = c4/t if (creqv .le. TINY) creqv = sqrt(TINY) delcr = ZERO IF(c1.gt.ZERO) delcr = (exp( log(c1) + c2 * log(seqv) + & c3 * log(creqv) - con1 )) * dtime dcrda(1)= c2 * delcr / seqv dcrda(2)= c3 * delcr / creqv Ustatev(nstatv) = creqv 990 continue return end 而我的方程是 delcr = (c1*(seqv**c2)*(((c3+1)**2)* (time**c3)+c4*(time**c5)+c6*(time**c7))*exp(-con1))*dtime 跪求大神看一看这个方程要写入ansys 该怎么修改上述程序！！！ 我试了修改了很多次，编译连接成功了，但是ansys分析时候，输入命令求解，但是ansys没有反应，我想肯定是程序输入的有错误现附我修改编写的程序 希望大神可以指教啊！！！前面声明变量都是一样的直接如下： *** skip when stress and creep strain are all zero if (seqv.LE.ZERO.AND.creqv.LE.ZERO) GO TO 990 c *** add temperature off set t = temp + toffst c *** Primary creep function c delcr := (c1*(seqv**c2)*(((c3+1)**2)*(time**c3)+c4*(time**c5)+c6*(time**c7))*exp(-c8/T)) * dtime c1 = prop(1) c2 = prop(2) c3 = prop(3) c4 = prop(4) c5 = prop(5) c6 = prop(6) c7 = prop(7) c8 = prop(8) c *** user need to make sure if c4 has nonzero value, temperature should be also nonzero. con1 = ZERO if(c8.ne.ZERO .and. t.gt.ZERO) con1 = c8/t delcr = Ustatev(1) !蠕变应变增量 c *** calculate incremental creep strain if (creqv .le. TINY) creqv = sqrt(TINY) delcr = ZERO IF(c1.gt.ZERO) delcr = (c1*(seqv**c2)*(((c3+1)**2)* & ((time+dtime)**c3)+c4*((time+dtime)**c5)+ & c6*((time+dtime)**c7))*exp(-con1))*dtime c *** 蠕变应变增量对等效应力的导数 dcrda(1) = c2 * delcr / seqv c *** 蠕变应变增量对等效蠕变应变的导数 & dcrda(2) = 0 c *** 保存蠕变应变增量 Ustatev(1) = delcr Ustatev(2) = time 990 continue return end 请大神指教啊 好人一生平安 |
回复此楼» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
-
物理学I论文润色/翻译怎么收费?
已经有229人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有0人回复
2楼2014-11-25 13:12:59
胡海浪
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1936.7
- 散金: 200
- 红花: 1
- 帖子: 295
- 在线: 93.9小时
- 虫号: 1842588
- 注册: 2012-05-31
- 性别: GG
- 专业: 应用数学方法
3楼2016-07-08 17:11:00
|
本帖内容被屏蔽 |
4楼2018-03-14 16:37:29
|
本帖内容被屏蔽 |
5楼2018-03-14 16:37:43
|
本帖内容被屏蔽 |
6楼2018-03-14 16:37:52
