| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 2462 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | ||
[求助]
求助ABAQUS軟件中環(huán)形焊接移動(dòng)熱源FORTRAN程序問(wèn)題 已有2人參與
|
||
|
大家好: 最近寫(xiě)了個(gè)程序想算一下環(huán)焊的溫度場(chǎng),前期是分段加載溫度曲線(xiàn)可以進(jìn)行計(jì)算。 后來(lái)想用FORTRAN寫(xiě)移動(dòng)熱源代入,繞圓周方向一圈。是拿直線(xiàn)焊縫的程序修改的,加了個(gè)平移和旋轉(zhuǎn)矩陣而已,平移是換了旋轉(zhuǎn)的圓心,旋轉(zhuǎn)是繞Y軸的。然后之后調(diào)理了一下焊槍角度又乘了個(gè)旋轉(zhuǎn)矩陣33度。這里繞Z軸的了。 結(jié)果發(fā)現(xiàn)提交的時(shí)候顯示編譯有問(wèn)題,我在VS里面建了IVF工程但是按了BUILD之后毫無(wú)反應(yīng),0 SUCCESS 0 FAIL 1 UPDATE 回頭又看了下程序可惜實(shí)在看不出問(wèn)題所在,請(qǐng)問(wèn)有沒(méi)有同學(xué)能指導(dǎo)一下我,告訴我可能哪里會(huì)存在問(wèn)題呢?謝謝啦! 附上環(huán)焊的INP和FOR文件。FOR文件沒(méi)傳上……就附在下面啦 subroutine usdfld(field,statev,pnewdt,direct,t,celent,time,dtime, 1 cmname,orname,nfield,nstatv,noel,npt,layer,kspt,kstep,kinc, 1 ndi,nshr,coord,jmac,jmtyp,matlayo,laccfla) c include 'aba_param.inc' c character*80 cmname,orname character*3 flgray(15) dimension field(nfield),statev(nstatv),direct(3,3),t(3,3),time(2) 1 ,coord(*),jmac(*),jmtyp(*) dimension array(15),jarray(15) c call getvrm('temp',array,jarray,flgray,jrcd, 1 jmac, jmtyp, matlayo, laccfla) c statev(1)=max(statev(1),array(1)) c return end c SUBROUTINE DFLUX(FLUX,SOL,KSTEP,KINC,TIME,NOEL,NPT,COORDS,JLTYP, 1 TEMP,PRESS,SNAME) c INCLUDE 'ABA_PARAM.INC' c DIMENSION FLUX(2),TIME(2),COORDS(3) CHARACTER*80 SNAME c v=10 t=TIME(1) FLUX(1)=0 PI=3.14159 R=50 h=0.5 waittime=0 angle=(v*TIME(1))/R ang0=0 theta=-angle+ang0 x=(COORDS(1)-R)*COS(theta)-COORDS(3)*SIN(theta) y=COORD(2)+h z=(COORDS(1)-R)*SIN(theta)+COORDS(3)*COS(theta) X1=(COORDS(1)-R)*COS(ang0)-COORDS(3)*SIN(ang0) X2=X2 X3=(COORDS(1)-R)*SIN(ang0)+COORDS(3)*COS(ang0) x=x*COS(-33./180.*PI)+y*SINC(-33./180.*PI) y=-x*COS(-33./180.*PI)*+y*COS(-33./180.*PI) z=z c ff=0.6 fr=2-ff a=6 b=2*a bf=b*ff*0.5 br=b*fr*0.5 c=3 s=v*t depth=1.5 heat=240*24*0.8*1000 arc1=ATAN(COORDS(2)/COORDS(3)) d01=exp(-3*(X1**2/(a)**2+X2**2/(bf)**2+X3**2/(c)**2)) d02=exp(-3*(X1**2/(a)**2+X2**2/(br)**2+X3**2/(c)**2)) d1=exp(-3*(x**2/(a)**2+y**2/(bf)**2+z**2/(c)**2)) d2=exp(-3*(x**2/(a)**2+y**2/(br)**2+z**2/(c)**2)) if((y.ge.0) $ .and.(TIME(2).gt.0).and.(TIME(2).le.waittime)) then FLUX(1)=heat*d01*0.6*6.0*sqrt(3.)/(pi*(a*bf*c)*sqrt(pi)) endif if((y.lt.0) $ .and.(TIME(2).gt.0).and.(TIME(2).le.waittime)) then FLUX(1)=heat*d02*1.4*6.0*sqrt(3.)/(pi*(a*br*c)*sqrt(pi)) endif if((arc1.gt.angle).and.(TIME(2).gt.waittime) $ .and.(TIME(2).le.22.5)) then FLUX(1)=heat*d1*0.6*6.0*sqrt(3.)/(pi*(a*bf*c)*sqrt(pi)) endif if((arc1.le.angle).and.(TIME(2).gt.waittime) $ .and.(TIME(2).le.22.5)) then FLUX(1)=heat*d2*1.4*6.0*sqrt(3.)/(pi*(a*br*c)*sqrt(pi)) endif c RETURN END |
新蟲(chóng) (初入文壇)
新蟲(chóng) (初入文壇)
新蟲(chóng) (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 285求調(diào)劑 +5 | AZMK 2026-03-30 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 370求調(diào)劑 +3 | 080700調(diào)劑 2026-03-30 | 3/150 |
|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑,總分281 +4 | 橙子(勝意) 2026-03-30 | 4/200 |
|
|
[有機(jī)交流] 考研調(diào)劑 +8 | watb 2026-03-26 | 8/400 |
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +12 | 芋泥小鈴鐺 2026-03-28 | 12/600 |
|
|
[考研] 材料專(zhuān)碩調(diào)劑 +11 | 椰椰。 2026-03-29 | 11/550 |
|
|
[考研] 070300化學(xué)354求調(diào)劑 +7 | 101次希望 2026-03-28 | 7/350 |
|
|
[考研] 0856求調(diào)劑 +8 | 楒桉 2026-03-28 | 8/400 |
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +8 | 柳溪邊 2026-03-26 | 8/400 |
|
|
[考研] 340求調(diào)劑 +6 | Amber00 2026-03-26 | 6/300 |
|
|
[考研] 本科雙非材料,跨考一志愿華電085801電氣,283求調(diào)劑,任何專(zhuān)業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
[考研] 312,生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺(jué)_16 2026-03-24 | 4/200 |
|
|
[考研] 復(fù)試調(diào)劑,一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
[考研] 各位老師您好:本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
[考研]
|
13659058978 2026-03-24 | 4/200 |
|