版塊導(dǎo)航: 正在加載中...

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

24小時(shí)熱門版塊排行榜

返回列表

下一頁(yè)

【懸賞金幣】回答本帖問題，作者杜杜杜茄子將贈(zèng)送您 50 個(gè)金幣

杜杜杜茄子

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1764
帖子: 290
在線: 7.8小時(shí)
蟲號(hào): 28153140
注冊(cè): 2022-01-13
專業(yè): 結(jié)構(gòu)工程

[求助] ansys體熱源溫度場(chǎng)應(yīng)力場(chǎng) 已有1人參與

求助

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

回復(fù)此樓

» 猜你喜歡

0856材料與化工，270求調(diào)劑已經(jīng)有15人回復(fù)
085600 材料與化工 298 已經(jīng)有12人回復(fù)
312求調(diào)劑已經(jīng)有7人回復(fù)
304分材料專碩求調(diào)劑已經(jīng)有9人回復(fù)
320材料與化工，求調(diào)劑已經(jīng)有5人回復(fù)
成果系統(tǒng)訪問量大，請(qǐng)15分鐘后再嘗試。由此給您造成的不便，敬請(qǐng)諒解。已經(jīng)有12人回復(fù)
281求調(diào)劑已經(jīng)有5人回復(fù)
347求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)
一志愿314求調(diào)劑已經(jīng)有6人回復(fù)

1樓 2023-02-18 13:16:29

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

回帖支持 ( 顯示支持度最高的前 50 名 )

htbbzzg

鐵桿木蟲 (著名寫手)

應(yīng)助: 668 (博士)
金幣: 9199
紅花: 146
帖子: 2205
在線: 346小時(shí)
蟲號(hào): 728202
注冊(cè): 2009-03-22
性別: GG
專業(yè): 動(dòng)力學(xué)與控制

【答案】應(yīng)助回帖

你好！
  簡(jiǎn)單看了一下你的命令流，先說一個(gè)問題吧：
  你所定義的與材料有關(guān)的命令有一個(gè)共同的問題，即：每個(gè)定義材料性能隨溫度變化的命令中的溫度參數(shù)最多只能是6 個(gè)，而你在每個(gè)材料性能隨溫度變化的命令中定義了8個(gè)參數(shù)，后兩個(gè)參數(shù)將會(huì)被忽略，導(dǎo)致計(jì)算結(jié)果不正確需要修改。比如以下兩個(gè)命令：
mptemp,1,0,250,500,750,1000,1500,1700,2500                      !定義材料溫度區(qū)間(C)
mpdata,dens,1,1,7850,7700,7600,7550,7500,7350,7300,7100       !定義材料密度(kg/m3)
  需要改為：
mptemp,1,0,250,500,750,1000,1500,                   !定義材料溫度區(qū)間(C) 前6個(gè)
mptemp,7,1700,2500                                        !定義材料溫度區(qū)間(C)    后2個(gè)，第2個(gè)參數(shù)表示從第7個(gè)溫度開始
mpdata,dens,1,1,7850,7700,7600,7550,7500,7350       !定義材料密度(kg/m3)    前6個(gè)
mpdata,dens,1,7,7300,7100                                     !定義材料密度(kg/m3)    后2個(gè)，第3個(gè)參數(shù)表示從第7個(gè)溫度開始

即：原來的每個(gè)命令拆分為兩個(gè)命令，保證每個(gè)命令中的數(shù)據(jù)參數(shù)不多于6個(gè)。
因?yàn)槲业氖虑楹芏啵日f這一點(diǎn)吧。你改一下，再試試，有問題再討論。

贊一下(1人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

杜杜杜茄子

7樓2023-02-22 08:41:29

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

普通回帖

杜杜杜茄子

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1764
帖子: 290
在線: 7.8小時(shí)
蟲號(hào): 28153140
注冊(cè): 2022-01-13
專業(yè): 結(jié)構(gòu)工程

各位老師好，我在學(xué)習(xí)ansys模擬焊縫殘余應(yīng)力時(shí)，仿照一篇論文的數(shù)據(jù)做的，出現(xiàn)了一些問題：
1.在溫度場(chǎng)中，第一層熱源溫度為2000多度，第二層為1600多度，論文里1900都為左右；且第一層焊縫焊到末尾時(shí)溫度突然上升到3000多度；后處理時(shí)，取焊縫背面中心節(jié)點(diǎn)熱循環(huán)曲線時(shí)溫度比論文里的數(shù)據(jù)小很多；
2.應(yīng)力場(chǎng)中，仿照論文里的約束方式，在立板頂端施加三向約束，得到的殘余應(yīng)力云圖形狀是錯(cuò)的，請(qǐng)大佬們幫忙指點(diǎn)！

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

贊一下

回復(fù)此樓

2樓2023-02-18 13:18:26

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

杜杜杜茄子

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1764
帖子: 290
在線: 7.8小時(shí)
蟲號(hào): 28153140
注冊(cè): 2022-01-13
專業(yè): 結(jié)構(gòu)工程

引用回帖:

2樓: Originally posted by 杜杜杜茄子 at 2023-02-18 13:18:26
各位老師好，我在學(xué)習(xí)ansys模擬焊縫殘余應(yīng)力時(shí)，仿照一篇論文的數(shù)據(jù)做的，出現(xiàn)了一些問題：
1.在溫度場(chǎng)中，第一層熱源溫度為2000多度，第二層為1600多度，論文里1900都為左右；且第一層焊縫焊到末尾時(shí)溫度突然上升到 ...

這是應(yīng)力場(chǎng)圖片

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

回復(fù)此樓

3樓2023-02-18 13:19:59

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

杜杜杜茄子

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1764
帖子: 290
在線: 7.8小時(shí)
蟲號(hào): 28153140
注冊(cè): 2022-01-13
專業(yè): 結(jié)構(gòu)工程

引用回帖:

3樓: Originally posted by 杜杜杜茄子 at 2023-02-18 13:19:59
這是應(yīng)力場(chǎng)圖片
...

發(fā)不了圖片，請(qǐng)私聊

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

贊一下

回復(fù)此樓

4樓2023-02-18 13:24:04

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

htbbzzg

鐵桿木蟲 (著名寫手)

應(yīng)助: 668 (博士)
金幣: 9199
紅花: 146
帖子: 2205
在線: 346小時(shí)
蟲號(hào): 728202
注冊(cè): 2009-03-22
性別: GG
專業(yè): 動(dòng)力學(xué)與控制

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

請(qǐng)你對(duì)自己的問題做一個(gè)詳細(xì)的介紹，并提供完整的計(jì)算模型 (最好是命令流)，以便進(jìn)行分析、討論。

贊一下

回復(fù)此樓

5樓2023-02-21 08:46:40

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

杜杜杜茄子

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1764
帖子: 290
在線: 7.8小時(shí)
蟲號(hào): 28153140
注冊(cè): 2022-01-13
專業(yè): 結(jié)構(gòu)工程

引用回帖:

5樓: Originally posted by htbbzzg at 2023-02-21 08:46:40
請(qǐng)你對(duì)自己的問題做一個(gè)詳細(xì)的介紹，并提供完整的計(jì)算模型 (最好是命令流)，以便進(jìn)行分析、討論。

老師，您好，我在學(xué)習(xí)用ansys體生熱率熱源做t型鋼板的溫度場(chǎng)和應(yīng)力場(chǎng)的分析，用的是一篇博士論文里的數(shù)據(jù)，工況是兩層焊縫，第一層焊接速度為0.01m/s，第二層為0.005m/s，中間間隔10s，焊后冷卻1800s，算出來的溫度場(chǎng)結(jié)果有問題。
主要的問題就是兩層焊縫的第一個(gè)荷載步和最后一個(gè)荷載步對(duì)應(yīng)的溫度會(huì)升高，而且兩層焊縫的溫度不一樣。第一層兩端為3000度左右，中間的荷載步為2000度，第二層兩端為1800度，中間為1600度左右，而論文里都是1880左右。
下面是命令流

!******************************第一部分進(jìn)行溫度場(chǎng)分析******************************
/com, units: SI, m - degree C - W - kg -sec - J
!/title,the drive force of solidification crack

/prep7
et,1,solid70
!******************************定義熱分析物理參數(shù)************************************
mptemp,1,0,250,500,750,1000,1500,1700,2500                      !定義材料溫度區(qū)間(C)
mpdata,dens,1,1,7850,7700,7600,7550,7500,7350,7300,7100       !定義材料密度(kg/m3)
mpdata,kxx,1,1,50,45,40,30,26,35,140,145             !熱傳導(dǎo)系數(shù)(W/(m*K))(考慮熔池對(duì)流）
mpdata,c,1,1,458,482,536,677,671,661,781,825                !比熱容(J/kg/K)
*dim,conve,table,8,1,1,temp       !考慮輻射和對(duì)流的總換熱系數(shù)(W/m^2*K)，加載時(shí)按對(duì)流加載  ，定義8行、1列、1面的數(shù)組
conve(0,1,1)=1
conve(1,0,1)=0
conve(1,1,1)=5
conve(2,0,1)=250
conve(2,1,1)=30
conve(3,0,1)=500
conve(3,1,1)=50
conve(4,0,1)=750
conve(4,1,1)=97
conve(5,0,1)=1000
conve(5,1,1)=150
conve(6,0,1)=1500
conve(6,1,1)=300
conve(7,0,1)=1700
conve(7,1,1)=307
conve(8,0,1)=2500
conve(8,1,1)=351
!******************************定義應(yīng)力分析參數(shù)***************************************
mpdata,ex,1,1,20E10,  17.8E10 ,13.5E10 ,11E10 ,0.7E10 ,0.5E10 ,0.5E10,0.1E10 !定義彈性模量

tb,bkin,1,6                                                       !定義屈服強(qiáng)度和切變模量
tbtemp,20,1
tbdata,1,235e6,2.06e10
tbtemp,100,2
tbdata,1,220e6,2.03e10
tbtemp,500,3
tbdata,1,160e6,1.47e10
tbtemp,800,4
tbdata,1,75e6,0.198e10
tbtemp,1000,5
tbdata,1,12e6,0.07e10
tbtemp,1500,6
tbdata,1,1e6,0.07e10

mpdata,alpx,1,1,1.1e-5,1.25e-5,1.4e-5,1.3e-5,1.3e-5,1.3e-5,1.3e-5,1.3e-5    !定義線膨脹系數(shù)(考慮凝固收縮中的固液區(qū)）
mpdata,prxy,1,1,0.29,0.29,0.39,0.39,0.39,0.39          !定義泊松比

!******************************建模型***********************************************
K, 1,0.0875,0,0
K, 2,0.0875,0.016,0
K, 3,0.0045,0.096,0
K, 4,-0.0045,0.096,0
K, 5,-0.0875,0.016,0
K, 6,-0.0875,0,0
K, 7, 0.0045, 0.024,0
K, 8,-0.0045, 0.016,0
K, 9,-0.0025, 0.016,0
K,10,0.0115, 0.016,0
k,11,0.0075, 0.016,0
k,12,0.0025, 0.152/7,0
K,13,-0.0045, 0,0
K,14,0.0115, 0,0
K,15,-0.0045, 0.024,0
K,16,-0.006, 0.016,0
K,17,-0.006, 0,0
K,18,0.013, 0.016,0
K,19,0.013, 0,0

A,1,2,18,19
A,3,4,15,7
A,5,6,17,16
A,8,13,14,10
A,9,11,12
A,11,10,7,12
A,10,14,19,18
A,8,9,7,15
A,16,17,13,8

VEXT, 1, , ,,,0.2
VEXT, 2, , ,,,0.2
VEXT, 3, , ,,,0.2
VEXT, 4, , ,,,0.2
VEXT, 5, , ,,,0.2
VEXT, 6, , ,,,0.2
VEXT, 7, , ,,,0.2
VEXT, 8, , ,,,0.2
VEXT, 9, , ,,,0.2

/facet,norml
vglue,all                                                    !粘連所有體

!************************************體10************************************
LESIZE,20,,,2
LESIZE,22,,,2
LESIZE,100,,,2
LESIZE,102,,,2

LESIZE,18,,,5
LESIZE,21,,,5
LESIZE,63,,,5
LESIZE,101,,,5

LESIZE,61,,,80
LESIZE,66,,,80
LESIZE,45,,,80
LESIZE,67,,,80

MSHAPE,0,3D$MSHKEY,1
mat,1
type,1
VMESH,10
!************************************體5************************************
LESIZE,65,,,80
LESIZE,18,,,5
LESIZE,19,,,5
LESIZE,17,,,10
LESIZE,63,,,5
LESIZE,64,,,5
LESIZE,62,,,10

mshkey,0
mshape,1

mat,1
type,1
VMESH,5

!************************************體3************************************
LESIZE,49,,,10
LESIZE,47,,,10
LESIZE,10,,,10
LESIZE,12,,,10

LESIZE,50,,,25
LESIZE,51,,,25
LESIZE,52,,,25
LESIZE,53,,,25

LESIZE,46,,,4
LESIZE,48,,,4
LESIZE,9,,,4
LESIZE,11,,,4
MSHAPE,0,3D$MSHKEY,1
mat,1
type,1
VMESH,3
!************************************體1************************************
LESIZE,2,,,10
LESIZE,4,,,10
LESIZE,31,,,10
LESIZE,33,,,10

LESIZE,34,,,25
LESIZE,35,,,25
LESIZE,36,,,25
LESIZE,37,,,25

LESIZE,1,,,4
LESIZE,3,,,4
LESIZE,30,,,4
LESIZE,32,,,4

MSHAPE,0,3D$MSHKEY,1
mat,1
type,1
VMESH,1
!************************************體2************************************
LESIZE,7,,,3
LESIZE,5,,,3
LESIZE,40,,,3
LESIZE,38,,,3

LESIZE,6,,,9
LESIZE,8,,,9
LESIZE,39,,,9
LESIZE,41,,,9

LESIZE,42,,,80
LESIZE,43,,,80
LESIZE,44,,,80

MSHAPE,0,3D$MSHKEY,1
mat,1
type,1
VMESH,2
!************************************體13************************************
LESIZE,107,,,1
LESIZE,25,,,1

LESIZE,108,,,7
LESIZE,27,,,7

LESIZE,58,,,80
LESIZE,65,,,80

mshkey,0
mshape,1
mat,1
type,1
vmesh,13
!************************************體14************************************
LESIZE,55,,,8
LESIZE,54,,,8
LESIZE,56,,,8
LESIZE,15,,,8
LESIZE,14,,,8
LESIZE,13,,,8

LESIZE,59,,,80
LESIZE,60,,,80

mshkey,0
mshape,1
mat,1
type,1
vmesh,14
!************************************體11************************************
mshkey,0
mshape,1
mat,1
type,1
vmesh,11

!************************************體12************************************
mshkey,0
mshape,1
mat,1
type,1
vmesh,12
!**************************************合并重合點(diǎn)*************************************
nummrg,all
alls
!***************************************死亡************************************************
/solu
/auto,1
/replot
alls
eplot
vsel,s,,,10
vsel,a,,,5
eslv,s,1
eplot
ekill,all
alls
esel,s,live
eplot
!************************施加對(duì)流邊界條件******************************
esel,s,live
ASEL,S,ext
NSLA,S,1
SF,ALL,CONV,%CONVE%,30
ALLS
!************************************非線性求解選項(xiàng)****************************************
/solu
antype,4
trnopt,full
pred,on
solcontrol,on
nropt,full,,on
cutcontrol,plslimit,10
lnsrch,on
autots,on
tintp,,,,1,,
timint,on
tref,30
kbc,0
neqit,50
outres,all
!****************************************定義變量*******************************************
!HGEN=Q/A*V*DT!!
!A為焊縫的橫截面積
!V為焊接速度
!DT為每個(gè)載荷步的時(shí)間
!Q=K*U*I*DT!!
!K為焊接熱源的熱效率
!U為焊接電壓
!I為焊接電流
!*******************************************穩(wěn)態(tài)分析*********************************************
timint,off
alls
ic,all,temp,30
alls
nsla,u,1
sf,all,conv,%conve%,30
time,1e-5
deltim,1e-5,1e-5,1e-5
kbc,1
allsel,all
solve
eplot
!********************************************第一道焊縫瞬態(tài)分析*******************************************
j=0.5             !
zmin=-1E-5
zmax=0.005+1E-5             !
zd=0.005
!***************************************計(jì)算生熱率******************************************
K=0.75
U=31
I=270
A=28e-6
V=0.01
DT=0.5             !
HGENUM=(K*U*I)/(A*V*DT)
timint,on
*do,i,1,40             !
/solu
vsel,s,,,5
eslv,s,1
nsle,s,all
nsel,r,loc,z,zmin,zmax
esln,s,all
ealive,all
eplot
bfe,all,hgen,,HGENUM
time,j
deltim,0.5,0.5,0.5
alls
solve

vsel,s,,,5
eslv,s,1

nsle,s,all
nsel,r,loc,z,zmin,zmax
esln,s,all
eplot

bfedele,all,hgen

j=j+0.5
zmin=zmin+zd
zmax=zmax+zd
*enddo
alls
eplot
!*************************************第一道焊后冷卻過程*************************************************

*do,i,41,50             !
/solu
time,j
deltim,1,1,1
alls
solve
j=j+1
*enddo
!********************************************第二道焊縫瞬態(tài)分析*******************************************
j=j+0.25
zmin1=-1E-5
zmax1=0.0025+1E-5
zd1=0.0025
!***************************************計(jì)算生熱率******************************************
A1=28e-6
V1=0.005
HGENUM1=(K*U*I)/(A1*V1*DT)
timint,on
*do,i,51,130             !
/solu
vsel,s,,,10
eslv,s,1
nsle,s,all
nsel,r,loc,z,zmin1,zmax1
esln,s,all
ealive,all
eplot
bfe,all,hgen,,HGENUM1
time,j
deltim,0.5,0.5,0.5
alls
solve
vsel,s,,,10
eslv,s,1
nsle,s,all
nsel,r,loc,z,zmin1,zmax1
esln,s,all
eplot
bfedele,all,hgen
j=j+0.5
zmin1=zmin1+zd1
zmax1=zmax1+zd1
*enddo
alls
eplot
!*************************************第二道焊后冷卻過程*************************************************

*do,i,131,310
/solu
time,j
deltim,10,10,10
alls
solve
save
j=j+10
*enddo

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

贊一下

回復(fù)此樓

6樓2023-02-21 12:28:10

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

杜杜杜茄子

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1764
帖子: 290
在線: 7.8小時(shí)
蟲號(hào): 28153140
注冊(cè): 2022-01-13
專業(yè): 結(jié)構(gòu)工程

送紅花一朵

引用回帖:

7樓: Originally posted by htbbzzg at 2023-02-22 08:41:29
你好！
簡(jiǎn)單看了一下你的命令流，先說一個(gè)問題吧：
你所定義的與材料有關(guān)的命令有一個(gè)共同的問題，即：每個(gè)定義材料性能隨溫度變化的命令中的溫度參數(shù)最多只能是6 個(gè)，而你在每個(gè)材料性能隨溫度變化的命令中定 ...

謝謝老師，感謝老師的指導(dǎo)��！

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

贊一下

回復(fù)此樓

8樓2023-02-22 19:10:58

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

杜杜杜茄子

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1764
帖子: 290
在線: 7.8小時(shí)
蟲號(hào): 28153140
注冊(cè): 2022-01-13
專業(yè): 結(jié)構(gòu)工程

引用回帖:

老師，您好，按您的指導(dǎo)，溫度云圖的數(shù)據(jù)好看了一些，我后來也修改了一下網(wǎng)格劃分，但是算下來的兩層焊縫的溫度還是有差距，第一層是2200度，第二層是1450度，因?yàn)槲覄澐值膬蓪雍缚p的截面面積是一樣的，不應(yīng)該會(huì)有這么大的差距，而且論文里的兩層焊縫得溫度都是1880度，應(yīng)該還是我哪里搞錯(cuò)了，請(qǐng)老師幫我看看。
下面是修改后的命令流
!******************************第一部分進(jìn)行溫度場(chǎng)分析******************************
/com, units: SI, m - degree C - W - kg -sec - J
!/title,the drive force of solidification crack

/prep7
et,1,solid70
!******************************定義熱分析物理參數(shù)************************************
mptemp,1,0,250,500,750,1000,1500
mptemp,7,1700,2500                      !定義材料溫度區(qū)間(C)
mpdata,dens,1,1,7850,7700,7600,7550,7500,7350
mpdata,dens,1,7,7300,7100       !定義材料密度(kg/m3)
mpdata,kxx,1,1,50,45,40,30,26,35
mpdata,kxx,1,7,140,145             !熱傳導(dǎo)系數(shù)(W/(m*K))(考慮熔池對(duì)流）
mpdata,c,1,1,458,482,536,677,671,661
mpdata,c,1,7,781,825                !比熱容(Jp/K)
mpdata,enth,1,1,0,250*7700*482,500*7600*536,750*7550*677,1000*7500*671,1500*7350*661
mpdata,enth,1,7,1700*7300*781,2500*7100*825
*dim,conve,table,8,1,1,temp       !考慮輻射和對(duì)流的總換熱系數(shù)(W/m^2*K)，加載時(shí)按對(duì)流加載  ，定義8行、1列、1面的數(shù)組
conve(0,1,1)=1
conve(1,0,1)=0
conve(1,1,1)=5
conve(2,0,1)=250
conve(2,1,1)=30
conve(3,0,1)=500
conve(3,1,1)=50
conve(4,0,1)=750
conve(4,1,1)=97
conve(5,0,1)=1000
conve(5,1,1)=150
conve(6,0,1)=1500
conve(6,1,1)=300
conve(7,0,1)=1700
conve(7,1,1)=307
conve(8,0,1)=2500
conve(8,1,1)=351
!******************************定義應(yīng)力分析參數(shù)***************************************
mpdata,ex,1,1,20E10,  17.8E10 ,13.5E10 ,11E10 ,0.7E10 ,0.5E10
mpdata,ex,1,7,0.5E10,0.1E10 !定義彈性模量
tb,bkin,1,6                                                       !定義屈服強(qiáng)度和切變模量
tbtemp,20,1
tbdata,1,235e6,2.06e10
tbtemp,100,2
tbdata,1,220e6,2.03e10
tbtemp,500,3
tbdata,1,160e6,1.47e10
tbtemp,800,4
tbdata,1,75e6,0.198e10
tbtemp,1000,5
tbdata,1,12e6,0.07e10
tbtemp,1500,6
tbdata,1,1e6,0.07e10

mpdata,alpx,1,1,1.1e-5,1.25e-5,1.4e-5,1.3e-5,1.3e-5,1.3e-5
mpdata,alpx,1,7,1.3e-5,1.3e-5    !定義線膨脹系數(shù)(考慮凝固收縮中的固液區(qū)）
mpdata,prxy,1,1,0.29,0.29,0.39,0.39
mpdata,prxy,1,7,0.39,0.39          !定義泊松比

!******************************建模型***********************************************
K, 1,0.0875,0,0
K, 2,0.0875,0.016,0
K, 3,0.0045,0.096,0
K, 4,-0.0045,0.096,0
K, 5,-0.0875,0.016,0
K, 6,-0.0875,0,0
K, 7, 0.0045, 0.024,0
K, 8,-0.0045, 0.016,0
K, 9,-0.0025, 0.016,0
K,10,0.0115, 0.016,0
k,11,0.0075, 0.016,0
k,12,0.0025, 0.152/7,0
K,13,-0.0045, 0,0
K,14,0.0115, 0,0
K,15,-0.0045, 0.024,0
K,16,-0.006, 0.016,0
K,17,-0.006, 0,0
K,18,0.013, 0.016,0
K,19,0.013, 0,0

A,1,2,18,19
A,3,4,15,7
A,5,6,17,16
A,8,13,14,10
A,9,11,12
A,11,10,7,12
A,10,14,19,18
A,8,9,7,15
A,16,17,13,8

VEXT, 1, , ,,,0.2
VEXT, 2, , ,,,0.2
VEXT, 3, , ,,,0.2
VEXT, 4, , ,,,0.2
VEXT, 5, , ,,,0.2
VEXT, 6, , ,,,0.2
VEXT, 7, , ,,,0.2
VEXT, 8, , ,,,0.2
VEXT, 9, , ,,,0.2

/facet,norml
vglue,all                                                    !粘連所有體

!*******************************網(wǎng)格劃分********
!******************體10和體5
vsel,s,,,10
vsel,a,,,5
esize,0.0025
mshkey,1
mshape,0
mat,1
type,1
vmesh,all
!******************體1********體2********體3
alls
vsel,s,,,1
vsel,a,,,2
vsel,a,,,3
esize,0.01
mshkey,1
mshape,0
mat,1
type,1
vmesh,all
!******************體14
alls
vsel,s,,,14
esize,0.005
mshkey,0
mshape,1
mat,1
type,1
vmesh,all
!******************體11****體12****體13
alls
vsel,s,,,11
vsel,a,,,12
vsel,a,,,13
esize,0.0025
mshkey,0
mshape,1
mat,1
type,1
vmesh,all
alls
!**************************************合并重合點(diǎn)*************************************
nummrg,all
alls
!************************************非線性求解選項(xiàng)****************************************
/solu
antype,4
trnopt,full
pred,on
solcontrol,on
nropt,full,,on
cutcontrol,plslimit,10
lnsrch,on
autots,on
tintp,,,,1,,
timint,on
tref,30
kbc,0
neqit,50
outres,all
!***************************************死亡************************************************
!/solu
!/auto,1
!/replot
alls
eplot
vsel,s,,,10
vsel,a,,,5
eslv,s,1
eplot
ekill,all
alls
esel,s,live
eplot
!************************施加對(duì)流邊界條件******************************
esel,s,live
ASEL,S,ext
NSLA,S,1
SF,ALL,CONV,%CONVE%,30
ALLS

!****************************************定義變量*******************************************
!HGEN=Q/A*V*DT!!
!A為焊縫的橫截面積
!V為焊接速度
!DT為每個(gè)載荷步的時(shí)間
!Q=K*U*I*DT!!
!K為焊接熱源的熱效率
!U為焊接電壓
!I為焊接電流
!*******************************************穩(wěn)態(tài)分析*********************************************
timint,off
alls
ic,all,temp,30
alls
nsla,u,1
sf,all,conv,%conve%,30
time,1e-5
deltim,1e-5,1e-5,1e-5
kbc,1
allsel,all
solve
eplot
!********************************************第一道焊縫瞬態(tài)分析*******************************************
j=0.5             !
zmin=-1E-5
zmax=0.005+1E-5             !
zd=0.005
!***************************************計(jì)算生熱率******************************************
K=0.75
U=31
I=270
A=28e-6
V=0.01
DT=0.5             !
HGENUM=(K*U*I)/(A*V*DT)
timint,on
*do,i,1,40             !
/solu
vsel,s,,,5
eslv,s,1
nsle,s,all
nsel,r,loc,z,zmin,zmax
esln,s,all
ealive,all
eplot
bfe,all,hgen,,HGENUM
time,j
deltim,0.5,0.5,0.5
alls
solve

vsel,s,,,5
eslv,s,1

nsle,s,all
nsel,r,loc,z,zmin,zmax
esln,s,all
eplot

bfedele,all,hgen

j=j+0.5
zmin=zmin+zd
zmax=zmax+zd
*enddo
alls
eplot
!*************************************第一道焊后冷卻過程*************************************************

*do,i,41,50             !
/solu
time,j
deltim,1,1,1
alls
solve
j=j+1
*enddo
!********************************************第二道焊縫瞬態(tài)分析*******************************************
j=j+0.5                !
zmin1=-1E-5
zmax1=0.0025+1E-5                !
zd1=0.0025                                  !
!***************************************計(jì)算生熱率******************************************
A1=28e-6
V1=0.005
DT1=0.5
HGENUM1=(K*U*I)/(A1*V1*DT1)
timint,on
*do,i,51,130             !
/solu
vsel,s,,,10
eslv,s,1
nsle,s,all
nsel,r,loc,z,zmin1,zmax1
esln,s,all
ealive,all
eplot
bfe,all,hgen,,HGENUM1
time,j
deltim,0.5,0.5,0.5                         !
alls
solve
vsel,s,,,10
eslv,s,1
nsle,s,all
nsel,r,loc,z,zmin1,zmax1
esln,s,all
eplot
bfedele,all,hgen
j=j+0.5
zmin1=zmin1+zd1
zmax1=zmax1+zd1
*enddo
alls
eplot
!*************************************第二道焊后冷卻過程*************************************************

*do,i,131,230
/solu
time,j
deltim,10,10,10
alls
solve
save
j=j+10
*enddo

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

贊一下

回復(fù)此樓

9樓2023-02-24 09:45:08

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

htbbzzg

鐵桿木蟲 (著名寫手)

應(yīng)助: 668 (博士)
金幣: 9199
紅花: 146
帖子: 2205
在線: 346小時(shí)
蟲號(hào): 728202
注冊(cè): 2009-03-22
性別: GG
專業(yè): 動(dòng)力學(xué)與控制

【答案】應(yīng)助回帖

你好！
試算了一下你上一個(gè)命令流，有一些問題：
  1 你的模型中存在一些不合格的單元，會(huì)影響計(jì)算的精度。
本想把不合格單元的圖片傳上來，但是不成功。你可以在前處理的 meshing 下面，檢查一下單元質(zhì)量，顯示不合格的單元，相應(yīng)命令是：
      check mesh > Individual Elm > Plot Warning/Err
  2 你的對(duì)流載荷施加了兩次，后一次比前一次多施加了一些面積，不知是為什么。
   本想把二者相差的面積圖片發(fā)上來，同樣發(fā)不上來。你可以自己檢查一下。方法是：
  在每一個(gè)施加對(duì)流載荷的命令之后，緊接著分別增加一個(gè) cm, a1, area  和 cm, a2，area 命令，把兩個(gè)對(duì)流載荷所施加的面積分別存放到 a1, a2 兩個(gè) component 中，然后執(zhí)行：
   cmsel, s, a2  >  cmsel, u, a1 >  cm,a3,area  命令，則 component a3 中就是兩個(gè)對(duì)流載荷施加面積之差。
  3  你的命令流中多處使用了 /solu 命令，包括循環(huán)中也如此。而 ANSYS 中，每一個(gè) /solu 命令都意味著開始新的一次分析，即不從上一次的結(jié)果繼續(xù)計(jì)算，而是從頭開始計(jì)算，這樣一來，分析結(jié)果的準(zhǔn)確度就無(wú)法保證了。所以，除了第一個(gè) /solu 命令外，只要你不是需要離開求解模塊或開始一個(gè)新的分析，后續(xù)的 /solu 命令都應(yīng)該去掉，以保證整個(gè)分析過程是連續(xù)的。如果確有需要離開求解模塊，再回到求解模塊時(shí)，必須使用重啟動(dòng)功能以保證從上一次的分析結(jié)果繼續(xù)計(jì)算；否則就會(huì)從頭開始計(jì)算，計(jì)算結(jié)果將會(huì)是錯(cuò)誤的。
  4  對(duì)所有的載荷步，你都是使用 deltim 命令定義子步的時(shí)間長(zhǎng)度的，對(duì)不同的載荷步需要分別定義子步長(zhǎng)，而且按照你的 deltim 命令，每個(gè)載荷步都只有一個(gè)子步，這也會(huì)影響結(jié)果的精度；建議改用 nsubst 命令，可以對(duì)所有的載荷步定義相同的子步數(shù)，比較方便。同時(shí)，一個(gè)載荷步一個(gè)子步，計(jì)算精度較差，最好多一些子步。改為，比如：  nsubst,2,5,2，這樣收斂較好是，可以一直 2 個(gè)子步；收斂不好時(shí)，可以增加子步數(shù)，有利于收斂。
  先說這些吧。我因?yàn)樯眢w不好，明天將會(huì)去醫(yī)院，住院開刀，有一段時(shí)間無(wú)法回答你的問題了，請(qǐng)見諒。

贊一下

回復(fù)此樓

10樓2023-02-26 11:26:34

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

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

返回列表

下一頁(yè)

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

普通表情龍兔虎貓

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

[考研] 求調(diào)劑 +5	博斯特525 2026-03-04	5/250	2026-03-04 14:07 by kakakapanpan
[考研] 一志愿314求調(diào)劑 +6	202111120625 2026-03-03	6/300	2026-03-04 14:05 by kakakapanpan
[考研] 材料工程269求調(diào)劑 +7	白刺玫 2026-03-02	7/350	2026-03-04 08:47 by zhyzzh
[論文投稿] EST拒稿重投 5+3	15102603076 2026-03-02	3/150	2026-03-04 00:51 by bobvan
[碩博家園] 2025屆雙非化工碩士畢業(yè)，申博 +4	更多的是 2026-02-27	5/250	2026-03-03 17:44 by wenium
[考研] 085602化學(xué)工程350，調(diào)劑，有沒有211的 +5	利好利好. 2026-03-02	9/450	2026-03-03 17:06 by 利好利好.
[考研] 289求調(diào)劑 +7	BrightLL 2026-03-02	9/450	2026-03-03 16:02 by tgxtgxtgx9
[考研] 成績(jī)276，專業(yè)代碼0856求調(diào)劑 +7	小陳朵 2026-03-03	7/350	2026-03-03 15:58 by tgxtgxtgx9
[考研] 材料270求調(diào)劑 6+6	Eiiiio 2026-03-01	11/550	2026-03-03 15:04 by Phoebelala
[考研] 課題組接收材料類調(diào)劑研究生 +6	gaoxiaoniuma 2026-02-28	9/450	2026-03-03 09:27 by gaoxiaoniuma
[考研] 298求調(diào)劑 +7	axyz3 2026-02-28	8/400	2026-03-03 06:11 by tgxtgxtgx9
[考研] 298求調(diào)劑 +10	人間唯你是清歡 2026-02-28	14/700	2026-03-02 22:49 by 人間唯你是清歡
[考研] 求調(diào)劑 +11	yunziaaaaa 2026-03-01	13/650	2026-03-02 21:59 by sunny81
[考研] 材料284求調(diào)劑，一志愿鄭州大學(xué)英一數(shù)二專碩 +15	想上岸的土撥鼠 2026-02-28	15/750	2026-03-02 20:13 by hypershenger
[考研] 材料調(diào)劑 +3	恒順自然 2026-03-02	3/150	2026-03-02 18:49 by L135790
[考研] 0856化工專碩求調(diào)劑 +15	董boxing 2026-03-01	15/750	2026-03-02 15:06 by 晃晃不許晃
[考研] 一志愿鄭大材料學(xué)碩298分，求調(diào)劑 +6	wsl111 2026-03-01	6/300	2026-03-02 11:00 by ydudjddnd
[考研] 材料調(diào)劑 +6	愛擦汗的可樂冰 2026-02-28	7/350	2026-03-02 10:42 by Jy?
[考研] 307求調(diào)劑 +4	73372112 2026-02-28	6/300	2026-03-01 00:04 by ll247
[考研] 304求調(diào)劑 +3	52hz~~ 2026-02-28	5/250	2026-03-01 00:00 by 52hz~~