版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

lifenhong

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 2.2
紅花: 1
帖子: 31
在線: 62.8小時
蟲號: 2979222
注冊: 2014-02-19
專業(yè): 金屬材料的合金相、相變及

[求助] marc模擬晶粒尺寸變化，能夠提交運算，其他的程序可以，這個不行，結(jié)果不顯示晶粒云圖

subroutine ugrain(m,n,nn,kc,mats,eplas,erate,dt,dtdl,
   * ignmod,grndat,grnsiz,time,deltime)
c
   implicit real*8 (a-h,o-z)
c
   dimension dt(1),dtdl(1),grndat(1)
c
c 計算晶粒大小的用戶子程序
c
c m 單元號
c n elsto number
c nn 積分點號
c kc 層號
c mats 材料號
c eplas 等效塑性應(yīng)變
c erate 等效塑性應(yīng)變率
c dt 增量步開始時的狀態(tài)變量(state variables)
c dtdl 狀態(tài)變量增量
c grndat 材料數(shù)據(jù)
c grndat(1) 指初始晶粒尺寸大小
c grnsiz 當(dāng)前晶粒尺寸大小(grain size)
c time 增量步開始時的時間
c deltime 時間增量
c ignmod 不同的模型的輸入模式 (如果是用戶定義的，則值為-1)
c
   real Nbeff,z,dsrx,dmrx,tsrx,tmrx,tip,a1,a2,a3,a4,a5,a6,epsc,Ep
   real Nb,Mn,Si,b1,bb
   include 'C:\MSC.Software\Marc\2010\marc2010\common/concom'
   include 'C:\MSC.Software\Marc\2010\marc2010\common/creeps'
c
cc
c 初始晶粒尺寸大小 grndat(1)
cc
   grndat(1)=205.9
cc
c 判斷等效塑性應(yīng)變率 erate 是否小于零
cc
   if (erate<0) then
   erate=0
   end if
cc
c 求解參數(shù) Z 和臨界應(yīng)變 epsc
cc
   z=erate*exp(306085/(8.314*(dt(1)+273)))
   epsc=5.6E-4*grndat(1)**0.3*z**0.17
c
   if (eplas>=0.04) then
c 當(dāng)?shù)刃苄詰?yīng)變 eplas>=0.04 時，
c 才會有再結(jié)晶過程發(fā)生
   if (eplas<epsc) then
c 當(dāng)?shù)刃苄詰?yīng)變小于臨界應(yīng)變時，即 eplas<epsc ，
c 發(fā)生靜態(tài)再結(jié)晶
   tip=abs(time+deltime-2.5)
c
   a1=exp(-45000.0/(8.31*(dt(1)+273)))
   a2=exp(230000.0/(8.31*(dt(1)+273)))
cc
c 靜態(tài)再結(jié)晶晶粒尺寸(微米) dsrx
c 發(fā)生 50%再結(jié)晶所需要的時間 tsrx
   dsrx=343.0*eplas**(-0.5)*grndat(1)**0.4*a1
   tsrx=2.3E-15*eplas**(-2.5)*grndat(1)**2.0*a2
cc
c 計算再結(jié)晶百分?jǐn)?shù) xv ，
c 并且判斷 xv 是否合理(應(yīng)該介于 0～ 1 之間)
   tt=tsrx ! t0.5
c xv=1-exp(-log(2.0)*(tip/tt)**1.5)
   xv=1-exp(-0.693*(tip/tt)**1.5)
c
   if (xv<0) then
   xv=0
   end if
   if (xv>1) then
   xv=1
   end if
cc
c 當(dāng)再結(jié)晶百分?jǐn)?shù) xv>=0.95 ，
c 晶粒尺寸為再結(jié)晶完成后的再結(jié)晶晶粒大�。�
c 當(dāng)再結(jié)晶百分?jǐn)?shù) xv<0.95 ，
c 晶粒尺寸為平均晶粒尺寸
   if (xv>=0.95) then
   grnsiz=dsrx
   else
   grnsiz=grndat(1)*(1-xv)+dsrx*xv
   end if
   else
c 當(dāng)?shù)刃苄詰?yīng)變大于或等于臨界應(yīng)變時，即 eplas>=epsc ，
c 發(fā)生動態(tài)再結(jié)晶
   tip=abs(time+deltime-10)
cc
c 動態(tài)再結(jié)晶晶粒尺寸(微米) Dmrx
c 發(fā)生 50%再結(jié)晶所需要的時間 tmrx
   dmrx=2.6E4*z**(-0.23)
   tmrx=0.4*z**(-0.8)*exp(240000.0/(8.31*(dt(1)+273)))
cc
c 計算再結(jié)晶百分?jǐn)?shù) xv ，
c 并且判斷 xv 是否合理(應(yīng)該介于 0～ 1 之間)
   tt=tmrx ! t0.5
   xv=1-exp(-0.693*(tip/tt))
c
   if (xv<0.0) then
      xv=0.0
   end if
   if (xv>1) then
   xv=1
   end if
cc
c 當(dāng)再結(jié)晶百分?jǐn)?shù) xv>=0.95 ，
c 晶粒尺寸為再結(jié)晶完成后的再結(jié)晶晶粒大�。�
c 當(dāng)再結(jié)晶百分?jǐn)?shù) xv<0.95 ，
c 晶粒尺寸為平均晶粒尺寸
   if (xv>=0.95) then
   grnsiz=dmrx
   else
   grnsiz=grndat(1)*(1-xv)+dmrx*xv
   end if
c
   end if
c
   else
c 當(dāng)?shù)刃苄詰?yīng)變 eplas<0.04 時，
c 沒有再結(jié)晶過程發(fā)生，
c 晶粒大小等于初始晶粒尺寸大小
   grnsiz=grndat(1)
   end if
c
   return
c
   end

回復(fù)此樓

» 猜你喜歡

英國巴斯大學(xué) University of Bath 博士后及博士招聘已經(jīng)有5人回復(fù)
求調(diào)劑收留已經(jīng)有1人回復(fù)
電氣科學(xué)與工程論文潤色/翻譯怎么收費? 已經(jīng)有195人回復(fù)
機械352分求調(diào)劑，本科太原理工大學(xué) 已經(jīng)有1人回復(fù)
085500機械專碩 331分求調(diào)劑已經(jīng)有1人回復(fù)
機械專碩311求調(diào)劑已經(jīng)有0人回復(fù)
(Nuclear) Physics of Nuclear Reactors-Academic(2021) 已經(jīng)有0人回復(fù)
(Nuclear) The fourth generation of nuclear reactors-Springer(2022) 已經(jīng)有0人回復(fù)
Heat Transfer Principles and Applications-Elsevier (2020) 已經(jīng)有0人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

1樓 2015-07-23 17:17:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hh12345bb

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 72
帖子: 3
在線: 12.1小時
蟲號: 3432576
注冊: 2014-09-22
專業(yè): 聚合物共混與復(fù)合材料

請問樓主解決了么，如果解決了還望不吝賜教啊

贊一下

回復(fù)此樓

2樓2015-10-08 15:37:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 lifenhong 的主題更新