 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 2633 | 回復(fù): 5 | ||
xiaomiao126金蟲 (小有名氣)
|
[求助]
關(guān)于ansys二次開發(fā)UPFs，修改usercreep子程序的問題，求好人幫忙 已有1人參與
|
|
大神求教�。。。。�！本人最近在做關(guān)于ansys UPFs的二次開發(fā)！想把自己的蠕變方程寫入ansys中！網(wǎng)上查資料知道修改 ansys 的usercreep程序可以實(shí)現(xiàn)此功能！ ansys 提供的程序可以實(shí)現(xiàn) 這個(gè)蠕變方程的編譯 delcr := c1 * seqv ^ n * creqv ^ m * exp (-b/T) * dtime c1，n，m，b都是常熟！delcr蠕變?cè)隽�，creqv蠕變應(yīng)變 他的程序如下: #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 跪求大神看一看這個(gè)方程要寫入ansys 該怎么修改上述程序�。�！ 我試了修改了很多次，編譯連接成功了，但是ansys分析時(shí)候，輸入命令求解，但是ansys沒有反應(yīng)，我想肯定是程序輸入的有錯(cuò)誤現(xiàn)附我修改編寫的程序 希望大神可以指教�。。�！前面聲明變量都是一樣的直接如下： *** 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) !蠕變應(yīng)變?cè)隽?br /> 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 *** 蠕變應(yīng)變?cè)隽繉?duì)等效應(yīng)力的導(dǎo)數(shù) dcrda(1) = c2 * delcr / seqv c *** 蠕變應(yīng)變?cè)隽繉?duì)等效蠕變應(yīng)變的導(dǎo)數(shù) & dcrda(2) = 0 c *** 保存蠕變應(yīng)變?cè)隽?br /> Ustatev(1) = delcr Ustatev(2) = time 990 continue return end 請(qǐng)大神指教啊 好人一生平安 |
回復(fù)此樓» 猜你喜歡
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有6人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有244人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有30人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有63人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
c2009ss
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 92.5
- 帖子: 7
- 在線: 2.3小時(shí)
- 蟲號(hào): 1615270
- 注冊(cè): 2012-02-13
- 專業(yè): 塑料
2樓2014-11-25 13:12:59
胡海浪
金蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1936.7
- 散金: 200
- 紅花: 1
- 帖子: 295
- 在線: 93.9小時(shí)
- 蟲號(hào): 1842588
- 注冊(cè): 2012-05-31
- 性別: GG
- 專業(yè): 應(yīng)用數(shù)學(xué)方法
3樓2016-07-08 17:11:00
|
本帖內(nèi)容被屏蔽 |
4樓2018-03-14 16:37:29
|
本帖內(nèi)容被屏蔽 |
5樓2018-03-14 16:37:43
|
本帖內(nèi)容被屏蔽 |
6樓2018-03-14 16:37:52
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xiaomiao126 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +5 | Doleres 2026-03-19 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研]
 271材料工程求調(diào)劑
+6
|
.6lL 2026-03-18 | 6/300 |
|
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑 +3 | 一瓶苯甲酸 2026-03-14 | 3/150 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +5 | z1z2z3879 2026-03-18 | 5/250 |
|
|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
|
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +4 | 材化逐夢(mèng)人 2026-03-18 | 4/200 |
|
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|