 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 2449 | 回復(fù): 5 | ||
[求助]
求助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ù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : Job-6.inp
2015-11-25 14:07:09, 21.93 M
» 猜你喜歡
281求調(diào)劑（0805）
已經(jīng)有8人回復(fù)
-
環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名
已經(jīng)有3人回復(fù)
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復(fù)
-
301求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有7人回復(fù)
-
302求調(diào)劑
已經(jīng)有6人回復(fù)
-
26博士申請(qǐng)
已經(jīng)有3人回復(fù)
-
268求調(diào)劑
已經(jīng)有3人回復(fù)
-
311求調(diào)劑
已經(jīng)有10人回復(fù)
-
被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了
已經(jīng)有4人回復(fù)
2樓2015-11-26 08:38:16
qinglove220
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 55.8
- 帖子: 35
- 在線: 13.6小時(shí)
- 蟲號(hào): 4241207
- 注冊(cè): 2015-11-23
- 專業(yè): 金屬材料的力學(xué)行為
3樓2016-01-05 08:45:22
4樓2017-03-31 07:01:04
wanghuanjie
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 80.3
- 散金: 5
- 帖子: 28
- 在線: 7.4小時(shí)
- 蟲號(hào): 4154533
- 注冊(cè): 2015-10-18
- 性別: GG
- 專業(yè): 制造系統(tǒng)與自動(dòng)化
5樓2018-04-26 21:40:53
tang2233
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 18
- 帖子: 3
- 在線: 33分鐘
- 蟲號(hào): 23811777
- 注冊(cè): 2020-10-15
- 專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域
6樓2022-03-05 10:50:37
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 owen2011 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 293求調(diào)劑 +7 | zjl的號(hào) 2026-03-16 | 12/600 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +3 | 困于星晨 2026-03-17 | 3/150 |
|
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[基金申請(qǐng)] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6 | zlingli 2026-03-13 | 6/300 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧， 2026-03-12 | 4/200 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
|
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +5 | 瑪卡巴卡啊哈 2026-03-11 | 5/250 |
|
|
|
[考研] 301求調(diào)劑 +6 | Liyouyumairs 2026-03-11 | 6/300 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
|
[考研] 314求調(diào)劑 +7 | 無懈可擊的巨人 2026-03-12 | 7/350 |
|
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|