 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (1134)
- >導(dǎo)師招生 (116)
- >考研 (105)
- >蟲(chóng)友互識(shí) (97)
- >考博 (31)
- >休閑灌水 (20)
- >碩博家園 (14)
- >博后之家 (13)
- >招聘信息布告欄 (10)
- >基金申請(qǐng) (10)
- >教師之家 (9)
- >論文投稿 (9)
- >找工作 (5)
- >醫(yī)學(xué) (5)
- >文獻(xiàn)求助 (4)
- >公派出國(guó) (4)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 764 | 回復(fù): 1 | ||
| 【懸賞金幣】回答本帖問(wèn)題，作者ganglu131將贈(zèng)送您 50 個(gè)金幣 | ||
ganglu131新蟲(chóng) (小有名氣)
|
[求助]
vumat求助
|
|
|
關(guān)于vumat求助 現(xiàn)在是報(bào)這個(gè)錯(cuò)誤，我網(wǎng)上查的都說(shuō)是程序里面有錯(cuò)誤，但是我找不出文件里面哪里有問(wèn)題，哪位大神可以幫忙找一下，非常感謝 Bad Material definition in element number 11 instance DERMIS-1: zero or negative initial dilatational modulus caused by bad material data. Please check your material input and any initial conditions if necessary. Abaqus/Explicit Packager exited with an error - Please see the status file for possible error messages if the file exists. parameter (zero = 0.d0, one = 1.d0, two = 2.d0, * three = 3.d0, third = 1.d0 / 3.d0, half = 0.5d0,four=4.d0) C10=0.2d0 D1=1e-5 ak=two/D1 amu=two*C10 twomu=four*C10 alamda=(three*ak-twomu)/three if (steptime.EQ.zero) then do i=1,nblock IF(cmname.eq.'B') THEN C10=0.2d0 D1=1e-5 ak=two/D1 amu=two*C10 twomu=four*C10 alamda=(three*ak-twomu)/three trace1=strainInc(i,1)+strainInc(i,2)+strainInc(i,3) stressNew(i,1)=stressOld(i,1)+ * twomu*strainInc(k,1)+alamda*trace1 stressNew(i,2)=stressOld(i,2)+ * twomu*strainInc(i,2)+alamda*trace1 stressNew(i,3)=stressOld(i,3)+ * twomu*strainInc(i,3)+alamda*trace1 stressNew(i,4)=stressOld(i,4)+ * twomu*strainInc(i,4) ELSEIF(cmname.eq.'A') THEN C10=7.2d0 D1=1e-5 ak=two/D1 amu=two*C10 twomu=four*C10 alamda=(three*ak-twomu)/three trace1=strainInc(i,1)+strainInc(i,2)+strainInc(i,3) stressNew(i,1)=stressOld(i,1)+ * twomu*strainInc(k,1)+alamda*trace1 stressNew(i,2)=stressOld(i,2)+ * twomu*strainInc(i,2)+alamda*trace1 stressNew(i,3)=stressOld(i,3)+ * twomu*strainInc(i,3)+alamda*trace1 stressNew(i,4)=stressOld(i,4)+ * twomu*strainInc(i,4) end if enddo else C do k=1,nblock IF(cmname.eq.'B') THEN C10=props(1) D1=props(2) ak=two/D1 amu=two*C10 twomu=four*C10 alamda=(three*ak-twomu)/three det=stretchNew(k,3)*(stretchNew(k,1)*stretchNew(k,2)- * stretchNew(k,4)**two) det2=(V(2,2)*V(1,1)-V(1,2)**two)*V(3,3) scale=det**( -ONE/THREE) stretchNewBar(k,1)=stretchNew(k,1)*scale stretchNewBar(k,2)=stretchNew(k,2)*scale stretchNewBar(k,3)=stretchNew(k,3)*scale stretchNewBar(k,4)=stretchNew(k,4)*scale BBar(k,1)=stretchNewBar(k,1)**two+stretchNewBar(k,4)**two BBar(k,2)=stretchNewBar(k,2)**two+stretchNewBar(k,4)**two BBar(k,3)=stretchNewBar(k,3)**two BBar(k,4)=stretchNewBar(k,1)*stretchNewBar(k,4) * +stretchNewBar(k,4)*stretchNewBar(k,2) TRBBar=BBar(k,1)+BBar(k,2)+BBar(k,3) EG=two*C10/det C PR=two/D1*(det-one) stressNew(K,1)=EG*(BBar(k,1)-TRBBar/Three)+PR stressNew(K,2)=EG*(BBar(k,2)-TRBBar/Three)+PR stressNew(K,3)=EG*(BBar(k,3)-TRBBar/Three)+PR stressNew(K,4)=EG*BBar(k,4) C Statenew(k,2)=sqrt(half*((StressNew(k,1)-StressNew(k,2))**two+ C * (StressNew(k,2)-StressNew(k,3))**two+(StressNew(k,3)- C * StressNew(k,1))**two)+three*StressNew(k,4)**two) s1=(stressNew(K,1)+stressNew(K,2))/two+ * sqrt(((stressNew(K,1)-stressNew(K,2))/two)**two+ *stressNew(K,4)**two) s2=(stressNew(K,1)+stressNew(K,2))/two- * sqrt(((stressNew(K,1)-stressNew(K,2))/two)**two+ *stressNew(K,4)**two) Statenew(k,2)=sqrt(((s1-s2)**two+(s1-stressNew(K,3))**two+ * (s2-stressNew(K,3))**two)/two) if(Statenew(k,2).gt.(7d0)) then statenew(k,1)=zero else statenew(k,1)=one end if if (stateold(k,1).eq.zero) then statenew(k,1)=zero end if C Update the enerInelasNew C C enerInelasNew(i) = enerInelasOld(i) + C * stressPower/density(i)-DEE enerInelasNew(k) = 0 ELSEIF(cmname.eq.'A') THEN C10=props(1) D1=props(2) ak=two/D1 amu=two*C10 twomu=four*C10 alamda=(three*ak-twomu)/three det=stretchNew(k,3)*(stretchNew(k,1)*stretchNew(k,2)- * stretchNew(k,4)**two) scale=det**(-ONE/THREE) stretchNewBar(k,1)=stretchNew(k,1)*scale stretchNewBar(k,2)=stretchNew(k,2)*scale stretchNewBar(k,3)=stretchNew(k,3)*scale stretchNewBar(k,4)=stretchNew(k,4)*scale BBar(k,1)=stretchNewBar(k,1)**two+stretchNewBar(k,4)**two BBar(k,2)=stretchNewBar(k,2)**two+stretchNewBar(k,4)**two BBar(k,3)=stretchNewBar(k,3)**two BBar(k,4)=stretchNewBar(k,1)*stretchNewBar(k,4) * +stretchNewBar(k,4)*stretchNewBar(k,2) C TRBBar=BBar(k,1)+BBar(k,2)+BBar(k,3) C BBar(k,1)=stretchNew(k,1)**two+stretchNew(k,4)**two-one C BBar(k,2)=stretchNew(k,2)**two+stretchNew(k,4)**two-one C BBar(k,3)=stretchNew(k,3)**two-one C BBar(k,4)=stretchNew(k,1)*stretchNew(k,4) C * +stretchNew(k,4)*stretchNew(k,2) TRBBar=BBar(k,1)+BBar(k,2)+BBar(k,3) EG=two*C10/det C EG=two*C10 PR=two/D1*(det-one) stressNew(K,1)=EG*(BBar(k,1)-TRBBar/Three)+PR stressNew(K,2)=EG*(BBar(k,2)-TRBBar/Three)+PR stressNew(K,3)=EG*(BBar(k,3)-TRBBar/Three)+PR stressNew(K,4)=EG*BBar(k,4) C Statenew(k,2)=sqrt(half*((StressNew(k,1)-StressNew(k,2))**two+ C * (StressNew(k,2)-StressNew(k,3))**two+(StressNew(k,3)- C * StressNew(k,1))**two)+three*StressNew(k,4)**two) s1=(stressNew(K,1)+stressNew(K,2))/two+ * sqrt(((stressNew(K,1)-stressNew(K,2))/two)**two+ *stressNew(K,4)**two) s2=(stressNew(K,1)+stressNew(K,2))/two- * sqrt(((stressNew(K,1)-stressNew(K,2))/two)**two+ *stressNew(K,4)**two) Statenew(k,2)=sqrt(((s1-s2)**two+(s1-stressNew(K,3))**two+ * (s2-stressNew(K,3))**two)/two) if(Statenew(k,2).gt.(26.d0)) then statenew(k,1)=zero else statenew(k,1)=one end if if (stateold(k,1).eq.zero) then statenew(k,1)=zero end if |
回復(fù)此樓» 猜你喜歡
2026級(jí)碩士研究生招生/調(diào)劑
已經(jīng)有5人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有7人回復(fù)
-
新疆大學(xué)地質(zhì)與礦業(yè)工程學(xué)院招生
已經(jīng)有32人回復(fù)
-
0856材料與化工353分求調(diào)劑
已經(jīng)有4人回復(fù)
-
282分材料專業(yè)求調(diào)劑院校
已經(jīng)有13人回復(fù)
-
315食品工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
環(huán)境調(diào)劑
已經(jīng)有3人回復(fù)
-
接受26屆調(diào)劑生
已經(jīng)有20人回復(fù)
-
考研材料與化工，求調(diào)劑
已經(jīng)有3人回復(fù)
-
歡迎211本科同學(xué)，過(guò)A區(qū)國(guó)家線，A區(qū)非偏遠(yuǎn)一本，交叉學(xué)科課題組
已經(jīng)有35人回復(fù)
ganglu131
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 517
- 帖子: 118
- 在線: 28分鐘
- 蟲(chóng)號(hào): 21357433
- 注冊(cè): 2020-03-05
2樓2022-06-16 19:34:54
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 ganglu131 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 264求調(diào)劑 +4 | thext 2026-03-03 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料與化工 一志愿山大 321分 求調(diào)劑 +4 | 每天散步 2026-03-09 | 5/250 |
|
|
|
[考研] 0817一志愿蘇州大學(xué)280復(fù)試調(diào)劑 +4 | kk扛 2026-03-07 | 4/200 |
|
|
|
[考研] 0832食品科學(xué)與工程293調(diào)劑 +3 | 東東不東 2026-03-07 | 3/150 |
|
|
|
[考研] 334求調(diào)劑 +5 | Trying] 2026-03-07 | 5/250 |
|
|
|
[考研] 求調(diào)劑，數(shù)一英一274分 +4 | 小菲會(huì)努力 2026-03-08 | 4/200 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-08 | 6/300 |
|
|
|
[考研] 材料調(diào)劑 +4 | xxxcm 2026-03-08 | 7/350 |
|
|
|
[考研] 一志愿武理314求調(diào)劑 +4 | （￣～￣;） 2026-03-08 | 5/250 |
|
|
|
[考研] 物理學(xué)求調(diào)劑 +3 | sx宋 2026-03-05 | 3/150 |
|
|
|
[考研] 求調(diào)劑 一志愿蘇州大學(xué)，0856化工323分 | 本科應(yīng)化 | 有專利/競(jìng)賽/科研助手經(jīng)歷 | +5 | 橙子cyx 2026-03-06 | 6/300 |
|
|
|
[考研] 316求調(diào)劑 +3 | 林小星發(fā)大財(cái) 2026-03-05 | 5/250 |
|
|
|
[考研] 328求調(diào)劑 +6 | vuzhdkfjkx 2026-03-04 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +4 | 曼殊2266 2026-03-05 | 4/200 |
|
|
|
[考研] 301求調(diào)劑 +4 | 李LJR 2026-03-04 | 4/200 |
|
|
|
[考研] 085701環(huán)境工程 求調(diào)劑 +7 | xiiiia 2026-03-04 | 7/350 |
|
|
|
[考研] 0856材料專碩274能調(diào)劑去哪里？ +3 | 22735 2026-03-04 | 4/200 |
|
|
|
[考研] 一志愿985材料與化工 326分求調(diào)劑 +3 | Hz795795 2026-03-04 | 3/150 |
|
|
|
[考研] 085601 材料305分求助 +4 | 泡泡郵件 2026-03-03 | 6/300 |
|
|
|
[考研] 331求調(diào)劑 +3 | zzZ&zZ 2026-03-03 | 3/150 |
|