 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 5012 | 回復(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é)與工程，一志愿華理能動285求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿華北電力大學(xué)能動專碩，293，求調(diào)劑
已經(jīng)有5人回復(fù)
-
調(diào)劑310
已經(jīng)有7人回復(fù)
-
食品工程專碩一志愿中海洋309求調(diào)劑
已經(jīng)有6人回復(fù)
-
275求調(diào)劑
已經(jīng)有14人回復(fù)
-
283求調(diào)劑
已經(jīng)有3人回復(fù)
-
085602 化工專碩 338分 求調(diào)劑
已經(jīng)有12人回復(fù)
-
312，生物學(xué)求調(diào)劑
已經(jīng)有3人回復(fù)
-
化學(xué)0703 調(diào)劑 306分 一志愿211
已經(jīng)有4人回復(fù)
-
311（085601）求調(diào)劑
已經(jīng)有3人回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分，前三科223分 +5 | dongfang59 2026-03-27 | 5/250 |
|
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
|
[考研] 07化學(xué)280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
|
[考研] 324求調(diào)劑 +8 | hanamiko 2026-03-26 | 10/500 |
|
|
|
[考研] 333求調(diào)劑 +6 | wfh030413@ 2026-03-23 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +8 | Auroracx 2026-03-22 | 8/400 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
|
[考研]
 085602化學(xué)工程求調(diào)劑。
+4
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
|
[考研]
 材料考研調(diào)劑生
+3
|
黃粱一夢千年 2026-03-24 | 3/150 |
|
|
|
[考研] 工科0856求調(diào)劑 +5 | 沐析汀汀 2026-03-21 | 5/250 |
|
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|