版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

lifenhong

新蟲 (初入文壇)

應助: 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 等效塑性應變
c erate 等效塑性應變率
c dt 增量步開始時的狀態(tài)變量(state variables)
c dtdl 狀態(tài)變量增量
c grndat 材料數(shù)據(jù)
c grndat(1) 指初始晶粒尺寸大小
c grnsiz 當前晶粒尺寸大小(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 判斷等效塑性應變率 erate 是否小于零
cc
   if (erate<0) then
   erate=0
   end if
cc
c 求解參數(shù) Z 和臨界應變 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 當?shù)刃苄詰?eplas>=0.04 時，
c 才會有再結(jié)晶過程發(fā)生
   if (eplas<epsc) then
c 當?shù)刃苄詰冃∮谂R界應變時，即 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é)晶百分數(shù) xv ，
c 并且判斷 xv 是否合理(應該介于 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 當再結(jié)晶百分數(shù) xv>=0.95 ，
c 晶粒尺寸為再結(jié)晶完成后的再結(jié)晶晶粒大��；
c 當再結(jié)晶百分數(shù) xv<0.95 ，
c 晶粒尺寸為平均晶粒尺寸
   if (xv>=0.95) then
   grnsiz=dsrx
   else
   grnsiz=grndat(1)*(1-xv)+dsrx*xv
   end if
   else
c 當?shù)刃苄詰兇笥诨虻扔谂R界應變時，即 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é)晶百分數(shù) xv ，
c 并且判斷 xv 是否合理(應該介于 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 當再結(jié)晶百分數(shù) xv>=0.95 ，
c 晶粒尺寸為再結(jié)晶完成后的再結(jié)晶晶粒大�。�
c 當再結(jié)晶百分數(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 當?shù)刃苄詰?eplas<0.04 時，
c 沒有再結(jié)晶過程發(fā)生，
c 晶粒大小等于初始晶粒尺寸大小
   grnsiz=grndat(1)
   end if
c
   return
c
   end

回復此樓

» 猜你喜歡

0858動力工程290求調(diào)劑已經(jīng)有4人回復
求調(diào)劑，認真踏實，吃苦耐勞，有科研潛力已經(jīng)有0人回復
工程熱物理與能源利用論文潤色/翻譯怎么收費? 已經(jīng)有143人回復
0858能源動力求B區(qū)調(diào)劑，動力工程也可以！已經(jīng)有0人回復
關于Chinese Journal of Aeronautics期刊收不收版面費的投稿求助貼已經(jīng)有0人回復
TPI層壓板已經(jīng)有0人回復
MOF與高分子配位交換反應復現(xiàn) 已經(jīng)有0人回復
調(diào)劑動力工程322 一志愿中石油華東已經(jīng)有2人回復
東南大學364求調(diào)劑已經(jīng)有5人回復
求調(diào)劑，總分391 已經(jīng)有3人回復
pva水凝膠電池已經(jīng)有0人回復

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

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

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

hh12345bb

銅蟲 (初入文壇)

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

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

贊一下

回復此樓

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

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

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