 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 5999 | 回復(fù): 2 | ||
ee522522新蟲 (初入文壇)
|
[求助]
幫忙啊specified range of 0 to 0 is not permitted. the esel command is ignored！
|
|
/clear,start /filename,weld /title,weld Analysis by "Element Birth and Death" /prep7 /unit,si!采用國(guó)際單位制 ！ et,1,13,4!13號(hào)二維耦合單元，同時(shí)具有溫度和位移自由度 et,2,13,4 ! mp,ex,1,2.3e11 mp,ex,2,2.05e11 ! tb,bkin,1,5 tbtemp,20,1 tbdata,1,2800e6,9.055e10 tbtemp,300,1 tbdata,1,2800e6,9.055e10 tbtemp,600,1 tbdata,1,2800e6,9.055e10 tbtemp,900,1 tbdata,1,2800e6,9.055e10 tbtemp,1200,1 tbdata,1,2800e6,9.055e10 ! tb,bkin,2,5 tbtemp,20,1 tbdata,1,335e6,8.2e10 tbtemp,300,1 tbdata,1,335e6,8.2e10 tbtemp,600,1 tbdata,1,335e6,8.2e10 tbtemp,900,1 tbdata,1,335e6,8.2e10 tbtemp,1200,1 tbdata,1,335e6,8.2e10 ! !材料密度 mptemp,1,20,300,600,900,1200 mpdata,dens,1,1,7700,7700,7700,7700,7700 mpdata,dens,2,1,7800,7800,7800,7800,7800 !熱膨脹系數(shù) mptemp,1,20,300,600,900,1200 mpdata,alpx,1,1,1.2e-5,1.2e-5,1.2e-5,1.2e-5,1.2e-5 mpdata,alpx,2,1,1.159e-5,1.31e-5,1.467e-5,1.65e-5,1.70e-5 !泊松比 mptemp,1,20,300,600,900,1200 mpdata,nuxy,1,1,0.27,0.27,0.27,0.27,0.27 mpdata,nuxy,2,1,0.25,0.25,0.25,0.25,0.25 !熱傳導(dǎo)系數(shù) mptemp,1,20,300,600,900,1200 mpdata,kxx,1,1,7.23,7.23,10.5,13.8,17 mpdata,kxx,2,1,50.2,50.2,50.2,50.2,50.2 !比熱 mptemp,1,20,300,600,755,1000 mpdata,c,1,1,3.73,3.73,3.73,3.73,3.73 mpdata,c,2,1,480,524,700,1064,602 !由于該13號(hào)單元有磁自由度，此處假設(shè)一個(gè)磁特性，但并不影響我們關(guān)心的結(jié)果 mptemp,1,20,300,600,900,1200 mpdata,murx,1,1,1,1,1,1,1 mpdata,murx,2,1,1,1,1,1,1 !假設(shè)涂層的初始溫度為1000 mp,reft,1,20 mp,reft,2,1000 mp,reft,3,20 ! !下面建立幾何模型 csys,0 k,1,0,0,0 k,2,0.06,0,0 k,3,0.06,0.005,0 k,4,0,0.005,0 k,5,0,0.010,0 k,6,0.06,0.010,0 a,1,2,3,4 a,3,4,5,6 !劃分單元 esize,0.0005 type,2 mat,2 MSHAPE,0,2D mshkey,1 amesh,1 esize,0.00005 type,1 mat,1 amesh,2 eplot ! /solu antype,4 trnopt,full ! !在模型下面施加位移約束 nesl,all *get,miny,node,,mnloc,y nsel,s,loc,y,miny d,all,uy,0 ! !假設(shè)模型的下邊界溫度始終為20 nsel,all *get,miny,node,,mnloc,y nsel,s,loc,y,miny d,all,temp,20 ! ！殺死涂層區(qū)域單元 nna=2 esel,all *get,emax,elem,,num,max asel,s,area,,nna esla *get,nse,elem,,count *dim,ne,,nse *dim,nex,,nse *dim,ney,,nse *dim,neorder,,nse mine=0 !下面的do循環(huán)用于將涂層區(qū)域的單元按其形心Y坐標(biāo)排序 !模擬涂層由下向上“生長(zhǎng)”的過程 *do,i1,1,nse esel,u,elem,,mine !此句有問題？ *get,nse1,elem,,count ii=0 *do,i,1,emax *if,esel(i),eq,1,then ii=ii+1 ne(ii)=i *endif *enddo *do,i,1,nse1 *get,ney(i),elem,ne(i),cent,y *get,nex(i),elem,ne(i),cent,x *enddo miny=1e20 minx=1e20 *do,i,1,nse1 *if,ney(i),lt,miny,then miny=ney(i) minx=nex(i) mine=ne(i) *else *if,ney(i),eq,miny,then *if,nex(i),lt,minx,then miny=ney(i) minx=nex(i) mine=ne(i) *endif *endif *endif *enddo neorder(i1)=mine *enddo ! max_tem=1000 dtl=1e-3 dt=8.33*20e-8 t=0 esel,all eplot /auto,1 /replot *do,i,1,nse ekill,neorder(i) esel,s,live eplot *enddo allsel,all outress,all,all ic,all,temp,20 kbc,1 timint,0,struct timint,1,therm timint,0,mag tintp,0.005,,,1,0.5,0.2 ! ! nsub1=2 nsub2=40 *do,i,1,nse ealive,neorder(i) esel,s,live eplot esel,all !下面的求解用于建立溫度的初始條件 t=t+dtl time,t nsubst,1 *do,j,1,4 d,nelem(neorder(i),j),temp,max_tem *enddo solve !下面的求解用于保證初始的升溫速度為零 t=t+dtl time,t solve !下面的步驟用于求解溫度分布 *do,j,1,4 ddele,nelem(neorder(i),j),temp *enddo t=t+dt-2*dtl time,t nsubst,nsub1 solve *enddo t=t+50000 time,t nsubst,nsub2 solve finish [ Last edited by ee522522 on 2013-11-8 at 10:27 ] |
回復(fù)此樓» 猜你喜歡
0854電子信息求調(diào)劑 324
已經(jīng)有3人回復(fù)
-
303求調(diào)劑
已經(jīng)有3人回復(fù)
-
336求調(diào)劑
已經(jīng)有3人回復(fù)
-
306求調(diào)劑
已經(jīng)有9人回復(fù)
-
工科0856求調(diào)劑
已經(jīng)有4人回復(fù)
-
工科材料085601 279求調(diào)劑
已經(jīng)有10人回復(fù)
-
081700 調(diào)劑 267分
已經(jīng)有4人回復(fù)
-
276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷
已經(jīng)有9人回復(fù)
-
263求調(diào)劑
已經(jīng)有9人回復(fù)
-
求老師收我
已經(jīng)有3人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 大神們，求幫忙啊，謝謝
已經(jīng)有3人回復(fù)
- 新手求助高斯輸入文件寫法，求大神幫忙
已經(jīng)有8人回復(fù)
- 求兩個(gè)值中最大的值但編寫的代碼不知道哪里錯(cuò)了，求大神幫忙。
已經(jīng)有10人回復(fù)
- 求助fortran大神，我是初學(xué)者，幫我分析一下下面的代碼
已經(jīng)有5人回復(fù)
- 求大神幫忙設(shè)計(jì)一個(gè)工裝夾具
已經(jīng)有6人回復(fù)
- 求各位大神幫忙下一篇文章
已經(jīng)有3人回復(fù)
- 跪求各位大神解答~本人高三學(xué)生 高考想考建筑學(xué) 求各位大神幫忙解答一下
已經(jīng)有9人回復(fù)
- 磷酸鐵鋰電池40度的容量居然比23度低，求大神幫忙解釋！有金幣感激！
已經(jīng)有12人回復(fù)
- 跪求大神幫忙解決A類錯(cuò)誤
已經(jīng)有3人回復(fù)
- 求大神幫忙研究一下關(guān)于鋁合金的抗彎強(qiáng)度的問題�。。。。。�！
已經(jīng)有13人回復(fù)
- 求各位大神幫忙 檢測(cè)豬肝中銅鐵鋅鉛的前處理方法有哪些�。。�
已經(jīng)有8人回復(fù)
- 哪位大神能幫忙總結(jié)下常見天然產(chǎn)物的預(yù)檢識(shí)方法啊，膜拜~~~~
已經(jīng)有3人回復(fù)
- 初學(xué)者請(qǐng)各位高手幫忙看一下GSAS精修的結(jié)果，請(qǐng)指點(diǎn)
已經(jīng)有14人回復(fù)
張大俠要學(xué)習(xí)
新蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 114.5
- 紅花: 4
- 帖子: 20
- 在線: 9.6小時(shí)
- 蟲號(hào): 3860437
- 注冊(cè): 2015-05-10
2樓2015-05-10 16:11:24
3樓2018-04-09 09:29:43
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 ee522522 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0854電子信息求調(diào)劑 324 +3 | Promise-jyl 2026-03-23 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 284求調(diào)劑 +6 | Zhao anqi 2026-03-22 | 6/300 |
|
|
[考研] 0703 調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-16 | 6/300 |
|
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
|
[考研] 307求調(diào)劑 +11 | 冷笙123 2026-03-17 | 11/550 |
|
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +5 | ￣^￣゜汗 2026-03-19 | 7/350 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
|
[考研] 0856材料專碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|
|
[考研] 280求調(diào)劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|
|
[考研] 求調(diào)劑 +4 | 要好好無(wú)聊 2026-03-21 | 4/200 |
|
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-16 | 3/150 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
|
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
|
[考研]
 320求調(diào)劑0856
+3
|
不想起名字112 2026-03-19 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|