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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

owen2011

新蟲 (初入文壇)

應助: 7 (幼兒園)
金幣: 102.7
紅花: 1
帖子: 40
在線: 13.5小時
蟲號: 1451222
注冊: 2011-10-19
專業(yè): 機械結構強度學

[求助] 求助ABAQUS軟件中環(huán)形焊接移動熱源FORTRAN程序問題已有2人參與

大家好：

最近寫了個程序想算一下環(huán)焊的溫度場，前期是分段加載溫度曲線可以進行計算。

后來想用FORTRAN寫移動熱源代入，繞圓周方向一圈。是拿直線焊縫的程序修改的，加了個平移和旋轉矩陣而已，平移是換了旋轉的圓心，旋轉是繞Y軸的。然后之后調(diào)理了一下焊槍角度又乘了個旋轉矩陣33度。這里繞Z軸的了。

結果發(fā)現(xiàn)提交的時候顯示編譯有問題，我在VS里面建了IVF工程但是按了BUILD之后毫無反應，0 SUCCESS 0 FAIL 1 UPDATE

回頭又看了下程序可惜實在看不出問題所在，請問有沒有同學能指導一下我，告訴我可能哪里會存在問題呢？謝謝啦！

附上環(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

回復此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權問題，其責任在于用戶本人，如對版權有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : Job-6.inp

2015-11-25 14:07:09, 21.93 M

» 猜你喜歡

277調(diào)劑已經(jīng)有6人回復
341求調(diào)劑已經(jīng)有7人回復
328求調(diào)劑，英語六級551，有科研經(jīng)歷已經(jīng)有8人回復
293求調(diào)劑已經(jīng)有12人回復
070300化學319求調(diào)劑已經(jīng)有4人回復
268求調(diào)劑已經(jīng)有9人回復
301求調(diào)劑已經(jīng)有4人回復
0703化學336分求調(diào)劑已經(jīng)有5人回復
301求調(diào)劑已經(jīng)有4人回復
材料專碩306英一數(shù)二已經(jīng)有9人回復

1樓 2015-11-25 14:08:08

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

owen2011

新蟲 (初入文壇)

應助: 7 (幼兒園)
金幣: 102.7
紅花: 1
帖子: 40
在線: 13.5小時
蟲號: 1451222
注冊: 2011-10-19
專業(yè): 機械結構強度學

好像發(fā)現(xiàn)了一點問題，

判定條件是 arc1=ATAN(COORDS(2)/COORDS(3)) 這里錯了，XYZ都搞錯了。

這個還是初始沒變化過相對著（0,0）原點的絕對角度，是不是也需要進行改寫 arc1=ATAN(COORDS(3)/COORDS(1)-R)

俺再試試看。。

贊一下(3人)

回復此樓

2樓2015-11-26 08:38:16

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

qinglove220

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 55.8
帖子: 35
在線: 13.6小時
蟲號: 4241207
注冊: 2015-11-23
專業(yè): 金屬材料的力學行為

【答案】應助回帖

請教一下要是建立軸對稱模型焊接的熱源應該如何施加

贊一下

回復此樓

3樓2016-01-05 08:45:22

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

luyizong

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 5
帖子: 1
在線: 1.3小時
蟲號: 6187900
注冊: 2017-03-29

【答案】應助回帖

引用回帖:

3樓: Originally posted by qinglove220 at 2016-01-05 08:45:22
請教一下要是建立軸對稱模型焊接的熱源應該如何施加

和整體一樣施加。只有半個熱源有效作用在半個體上。所以和原來一樣

發(fā)自小木蟲IOS客戶端

贊一下

回復此樓

4樓2017-03-31 07:01:04

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

wanghuanjie

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 80.3
散金: 5
帖子: 28
在線: 7.4小時
蟲號: 4154533
注冊: 2015-10-18
性別: GG
專業(yè): 制造系統(tǒng)與自動化

樓主，請教下，如何讓熱源沿圓周移動呢？能指出具體程序嗎？謝謝。

贊一下(1人)

回復此樓

5樓2018-04-26 21:40:53

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

tang2233

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 18
帖子: 3
在線: 33分鐘
蟲號: 23811777
注冊: 2020-10-15
專業(yè): 工程熱物理相關交叉領域

引用回帖:

5樓: Originally posted by wanghuanjie at 2018-04-26 21:40:53
樓主，請教下，如何讓熱源沿圓周移動呢？能指出具體程序嗎？謝謝。

老哥環(huán)形移動熱源咋整呀

贊一下

回復此樓

6樓2022-03-05 10:50:37

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

相關版塊跳轉我要訂閱樓主 owen2011 的主題更新

返回列表

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

[考研] 268求調(diào)劑 +8	一定有學上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 梁成偉老師課題組歡迎你的加入 +8	一鴨鴨喲 2026-03-14	10/500	2026-03-17 15:07 by 一鴨鴨喲
[考研] 材料與化工專碩調(diào)劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考研] 08工科 320總分求調(diào)劑 +4	梨花珞晚風 2026-03-17	4/200	2026-03-17 13:38 by houyaoxu
[考研] 環(huán)境工程調(diào)劑 +6	大可digkids 2026-03-16	6/300	2026-03-16 17:16 by barlinike
[考研] 304求調(diào)劑 +4	ahbd 2026-03-14	4/200	2026-03-16 16:48 by 我的船我的海
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 一志愿華中師范071000，325求調(diào)劑 +6	RuitingC 2026-03-12	6/300	2026-03-16 14:50 by 可淡不可忘
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[考研] 341求調(diào)劑 +4	番茄頭--- 2026-03-10	4/200	2026-03-13 23:12 by JourneyLucky
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[考研] 26調(diào)劑/材料科學與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 310求調(diào)劑 +3	【上上簽】 2026-03-11	3/150	2026-03-13 16:16 by JourneyLucky
[考博] 2026年博士申請 +3	QwQwQW10 2026-03-11	3/150	2026-03-12 17:58 by gxch43
[考研] 大連大學化學專業(yè)研究生調(diào)劑 +3	琪久. 2026-03-10	8/400	2026-03-11 10:02 by 琪久.
[考研] 279求調(diào)劑 +3	莫xiao 2026-03-10	4/200	2026-03-11 08:06 by 斬魂滴兔子！

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

24小時熱門版塊排行榜

owen2011

[求助] 求助ABAQUS軟件中環(huán)形焊接移動熱源FORTRAN程序問題 已有2人參與

» 本帖附件資源列表

» 猜你喜歡

owen2011

qinglove220

【答案】應助回帖

luyizong

【答案】應助回帖

wanghuanjie

tang2233

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

[求助] 求助ABAQUS軟件中環(huán)形焊接移動熱源FORTRAN程序問題已有2人參與