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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

【懸賞金幣】回答本帖問題，作者ganglu131將贈送您 50 個金幣

ganglu131

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 517
帖子: 118
在線: 28分鐘
蟲號: 21357433
注冊: 2020-03-05

[求助] vumat求助

關(guān)于vumat求助
現(xiàn)在是報這個錯誤，我網(wǎng)上查的都說是程序里面有錯誤，但是我找不出文件里面哪里有問題，哪位大神可以幫忙找一下，非常感謝
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

回復此樓

» 猜你喜歡

288求調(diào)劑已經(jīng)有12人回復
化工284求調(diào)劑已經(jīng)有10人回復
一志愿北京化工大學材料與化工296分求調(diào)劑已經(jīng)有8人回復
2026級碩士研究生招生/調(diào)劑已經(jīng)有5人回復
0703化學調(diào)劑已經(jīng)有7人回復
新疆大學地質(zhì)與礦業(yè)工程學院招生已經(jīng)有32人回復
0856材料與化工353分求調(diào)劑已經(jīng)有4人回復
282分材料專業(yè)求調(diào)劑院校已經(jīng)有13人回復
315食品工程求調(diào)劑已經(jīng)有5人回復
環(huán)境調(diào)劑已經(jīng)有3人回復

1樓 2022-06-16 14:28:42

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

ganglu131

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 517
帖子: 118
在線: 28分鐘
蟲號: 21357433
注冊: 2020-03-05

有人嗎

回復此樓

2樓2022-06-16 19:34:54

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

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

返回列表

不應(yīng)助 確定回帖應(yīng)助 (注意：應(yīng)助才可能被獎勵，但不允許灌水，必須填寫15個字符以上)

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[碩博家園] 2026級碩士研究生招生/調(diào)劑 +4	lbj6746988 2026-03-03	5/250	2026-03-10 05:25 by Ariana0818
[考研] 083000環(huán)境科學與工程調(diào)劑 +3	mingmingry 2026-03-09	3/150	2026-03-09 22:47 by bingxueer79
[考研] 085600材料與化工，一志愿廣州985，求調(diào)劑 +14	qqyyaill 2026-03-05	14/700	2026-03-09 20:16 by 霧散后相遇lc
[考研] 310 070300化學求調(diào)劑 +4	撲風鈴的貓 2026-03-08	5/250	2026-03-09 19:46 by 30660438
[考研] 材料與化工求調(diào)劑 +8	與冬清寧 2026-03-07	10/500	2026-03-09 13:45 by Wahxp
[考研] 0817找調(diào)劑 +6	kk扛 2026-03-08	6/300	2026-03-09 06:38 by houyaoxu
[考研] 0703化學求調(diào)劑，總分309分，一志愿華南師范 +3	花與葉@ 2026-03-08	3/150	2026-03-09 05:48 by houyaoxu
[教師之家] 交大前校長王樹國：現(xiàn)在最先進的科技并不在大學實驗室，而是在企業(yè)研究院 +4	zju2000 2026-03-08	6/300	2026-03-08 19:15 by zju2000
[考研] 一志愿211 化學305分求調(diào)劑 +3	0703楊悅305分 2026-03-05	3/150	2026-03-07 21:15 by yinhuanshun
[考研] 267化工調(diào)劑求助 +7	聰少OZ 2026-03-04	7/350	2026-03-07 17:33 by 無際的草原
[考研] 2026調(diào)劑】考試A區(qū)0703化學類323分誠求接收 +3	卷柏卷柏 2026-03-05	4/200	2026-03-07 17:00 by yuandd_2001
[考研] 328求調(diào)劑 +6	vuzhdkfjkx 2026-03-04	7/350	2026-03-05 23:28 by Leeding1356
[考研] 材料與化工，291，求調(diào)劑 +7	咕嚕咕嚕123123 2026-03-05	8/400	2026-03-05 20:08 by houyaoxu
[考研] 304求調(diào)劑 +4	曼殊2266 2026-03-05	4/200	2026-03-05 17:10 by zhukairuo
[考研] 334求調(diào)劑 +3	Trying] 2026-03-03	4/200	2026-03-05 16:39 by zhukairuo
[考研] 085701環(huán)境工程求調(diào)劑 +7	xiiiia 2026-03-04	7/350	2026-03-05 09:49 by 熱情沙漠
[考研] 085600 材料與化工 298 +14	小西笑嘻嘻 2026-03-03	14/700	2026-03-04 16:23 by sslc1985
[考研] 能動297求調(diào)劑，本科川大 +4	邵11 2026-03-04	4/200	2026-03-04 15:29 by lylylyheart
[考研] 306求調(diào)劑 +7	張張張張oo 2026-03-03	7/350	2026-03-03 21:35 by L135790
[考研] 化學0703求調(diào)劑學碩理/工科均可總分279 +3	1一11 2026-03-03	5/250	2026-03-03 12:37 by 1一11

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

ganglu131

[求助] vumat求助

» 猜你喜歡

ganglu131

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频