| 5 | 1/1 | 返回列表 |
| 查看: 2460 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
求助ABAQUS軟件中環(huán)形焊接移動(dòng)熱源FORTRAN程序問題 已有2人參與
|
||
|
大家好: 最近寫了個(gè)程序想算一下環(huán)焊的溫度場(chǎng),前期是分段加載溫度曲線可以進(jìn)行計(jì)算。 后來想用FORTRAN寫移動(dòng)熱源代入,繞圓周方向一圈。是拿直線焊縫的程序修改的,加了個(gè)平移和旋轉(zhuǎn)矩陣而已,平移是換了旋轉(zhuǎn)的圓心,旋轉(zhuǎn)是繞Y軸的。然后之后調(diào)理了一下焊槍角度又乘了個(gè)旋轉(zhuǎn)矩陣33度。這里繞Z軸的了。 結(jié)果發(fā)現(xiàn)提交的時(shí)候顯示編譯有問題,我在VS里面建了IVF工程但是按了BUILD之后毫無反應(yīng),0 SUCCESS 0 FAIL 1 UPDATE 回頭又看了下程序可惜實(shí)在看不出問題所在,請(qǐng)問有沒有同學(xué)能指導(dǎo)一下我,告訴我可能哪里會(huì)存在問題呢?謝謝啦! 附上環(huán)焊的INP和FOR文件。FOR文件沒傳上……就附在下面啦 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 |
新蟲 (初入文壇)
新蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +7 | 十閑wx 2026-03-24 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑,總分281 +3 | 橙子(勝意) 2026-03-30 | 3/150 |
|
|
[考研] 327求調(diào)劑 +5 | 小卡不卡. 2026-03-29 | 5/250 |
|
|
[考研] 085600,專業(yè)課化工原理,320分求調(diào)劑 +5 | 大饞小子 2026-03-29 | 5/250 |
|
|
[考研] 11408總分309,一志愿東南大學(xué)求調(diào)劑,不挑專業(yè) +4 | 天賦帶到THU 2026-03-29 | 5/250 |
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +8 | 26要上岸 2026-03-28 | 8/400 |
|
|
[考研] 329求調(diào)劑,一志愿西北工業(yè)大學(xué),材料工程(085601) +5 | 小小機(jī)靈蟲 2026-03-29 | 11/550 |
|
|
[考研] 材料與化工272求調(diào)劑 +21 | 阿斯蒂芬2004 2026-03-28 | 21/1050 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +6 | @taotao 2026-03-26 | 7/350 |
|
|
[考研] 290求調(diào)劑 +3 | dfffsar 2026-03-29 | 3/150 |
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動(dòng)285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
[考研] 壓國(guó)家一區(qū)線,求導(dǎo)師收留,有恩必謝! +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩,求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺_16 2026-03-24 | 4/200 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研] 303求調(diào)劑 +6 | 藍(lán)山月 2026-03-25 | 6/300 |
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
[考研] 各位老師您好:本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|