 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 5013 | 回復(fù): 1 | |||
zc_1981銅蟲 (小有名氣)
|
[交流]
abaqus中場變量子程序中的單、雙精度計算問題 已有1人參與
|
|
SUBROUTINE USDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT, 1 TIME,DTIME,CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER, 2 KSPT,KSTEP,KINC,NDI,NSHR,COORD,JMAC,JMATYP,MATLAYO, 3 LACCFLA) C INCLUDE 'ABA_PARAM.INC' DOUBLE PRECISION a,b,c,d,e,f C CHARACTER*80 CMNAME,ORNAME CHARACTER*3 FLGRAY(15) DIMENSION FIELD(NFIELD),STATEV(NSTATV),DIRECT(3,3), 1 T(3,3),TIME(2) DIMENSION ARRAY(15),JARRAY(15),JMAC(*),JMATYP(*), 1 COORD(*) C C Absolute value of current strain: CALL GETVRM('POR',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP, 1 MATLAYO,LACCFLA) P1=ARRAY(1) CALL GETVRM('TEMP',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP, 1 MATLAYO,LACCFLA) T1=ARRAY(1) C field variable FIELD(1)=P1 C state variable C IF(P1.GE.7.5*(10**6).AND.T1.GE.308)THEN DOUBLE PRECISION c,e,f ,DOUBLE PRECISION a,b,c,d,e,f ,END IF(T1.GE.284.0D0.AND.T1.LE.292.0D0) THEN a=5.0D-4*(T1-273.15D0)*(T1-273.15D0) b=10.0D0**(3.42D-2*(T1-273.15D0)+a+6.4804D0) c=DEXP(-9400.0D0/T1) d=(b-8.0D6)/1.0D6 e=(0.585D0*(10.0**1.0D13))*c*d f=0.5D0*(1.0D0-e)*0.4D0 STATEV(1)=(0.6*46*1.0D6+9313.264*1.0D6*f)/1.0D6 END IF END C If error, write comment to .DAT file: C IF(JRCD.NE.0)THEN C WRITE(6,*) 'REQUEST ERROR IN USDFLD FOR ELEMENT NUMBER ', C 1 NOEL,'INTEGRATION POINT NUMBER ',NPT C ENDIF C RETURN END 這是本人編的一個子程序，在默認(rèn)的單精度條件下去計算，可以順利過去，可是由于有的數(shù)值比較小，超過了單精度的問題，所以都改成了雙精度去計算，同樣的程序，為什么改成雙精度去計算就不好使了呢？出現(xiàn)下面這樣的錯誤提示“A return statement is invalid in the main programm”.是不是場變量的子程序里，不能用雙精度來計算呢？ |
回復(fù)此樓» 猜你喜歡
生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績338
已經(jīng)有6人回復(fù)
-
321求調(diào)劑
已經(jīng)有7人回復(fù)
-
343求調(diào)劑
已經(jīng)有4人回復(fù)
-
305求調(diào)劑
已經(jīng)有3人回復(fù)
-
本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑
已經(jīng)有6人回復(fù)
-
332求調(diào)劑
已經(jīng)有3人回復(fù)
-
學(xué)碩274求調(diào)劑
已經(jīng)有8人回復(fù)
-
壓國家一區(qū)線，求導(dǎo)師收留，有恩必謝！
已經(jīng)有7人回復(fù)
-
081200-11408-276學(xué)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
22408 359分調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- Abaqus的子程序問題
已經(jīng)有5人回復(fù)
- ABAQUS中的用戶子程序fortran編寫問題
已經(jīng)有4人回復(fù)
chen201023
木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 6870.2
- 散金: 80
- 紅花: 3
- 帖子: 330
- 在線: 109.8小時
- 蟲號: 1311273
- 注冊: 2011-05-31
- 專業(yè): 金屬結(jié)構(gòu)材料
2樓2020-05-15 14:35:53
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zc_1981 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
|
[考研] 085600，材料與化工321分求調(diào)劑 +9 | 大饞小子 2026-03-28 | 9/450 |
|
|
|
[考研] 283求調(diào)劑 +7 | A child 2026-03-28 | 7/350 |
|
|
|
[考研] 一志愿上海理工能源動力（085800）310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
|
[考研] 283求調(diào)劑（080500） +4 | A child 2026-03-27 | 4/200 |
|
|
|
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
|
[考研] 291 求調(diào)劑 +7 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 8/400 |
|
|
|
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
|
[考研] 考研調(diào)劑 +6 | 來好運來來來 2026-03-21 | 7/350 |
|
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
|
[考研] 化工專碩求調(diào)劑 +3 | question挽風(fēng) 2026-03-24 | 3/150 |
|
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|