 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 6000 | 回復(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!采用國際單位制 ！ et,1,13,4!13號二維耦合單元，同時(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號單元有磁自由度，此處假設(shè)一個磁特性，但并不影響我們關(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)排序 !模擬涂層由下向上“生長”的過程 *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ù)此樓» 猜你喜歡
08工學(xué)調(diào)劑
已經(jīng)有11人回復(fù)
-
350求調(diào)劑
已經(jīng)有6人回復(fù)
-
求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601）
已經(jīng)有4人回復(fù)
-
一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑
已經(jīng)有6人回復(fù)
-
接收2026碩士調(diào)劑(學(xué)碩+專碩)
已經(jīng)有6人回復(fù)
-
081700 調(diào)劑 267分
已經(jīng)有5人回復(fù)
-
328求調(diào)劑
已經(jīng)有4人回復(fù)
-
0854 考研調(diào)劑 招生了！AI 方向
已經(jīng)有14人回復(fù)
-
招08考數(shù)學(xué)
已經(jīng)有14人回復(fù)
-
一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 大神們，求幫忙啊，謝謝
已經(jīng)有3人回復(fù)
- 新手求助高斯輸入文件寫法，求大神幫忙
已經(jīng)有8人回復(fù)
- 求兩個值中最大的值但編寫的代碼不知道哪里錯了，求大神幫忙。
已經(jīng)有10人回復(fù)
- 求助fortran大神，我是初學(xué)者，幫我分析一下下面的代碼
已經(jīng)有5人回復(fù)
- 求大神幫忙設(shè)計(jì)一個工裝夾具
已經(jīng)有6人回復(fù)
- 求各位大神幫忙下一篇文章
已經(jīng)有3人回復(fù)
- 跪求各位大神解答~本人高三學(xué)生 高考想考建筑學(xué) 求各位大神幫忙解答一下
已經(jīng)有9人回復(fù)
- 磷酸鐵鋰電池40度的容量居然比23度低，求大神幫忙解釋！有金幣感激！
已經(jīng)有12人回復(fù)
- 跪求大神幫忙解決A類錯誤
已經(jīng)有3人回復(fù)
- 求大神幫忙研究一下關(guān)于鋁合金的抗彎強(qiáng)度的問題！�。。。。。�
已經(jīng)有13人回復(fù)
- 求各位大神幫忙 檢測豬肝中銅鐵鋅鉛的前處理方法有哪些�。�！
已經(jīng)有8人回復(fù)
- 哪位大神能幫忙總結(jié)下常見天然產(chǎn)物的預(yù)檢識方法啊，膜拜~~~~
已經(jīng)有3人回復(fù)
- 初學(xué)者請各位高手幫忙看一下GSAS精修的結(jié)果，請指點(diǎn)
已經(jīng)有14人回復(fù)
張大俠要學(xué)習(xí)
新蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 114.5
- 紅花: 4
- 帖子: 20
- 在線: 9.6小時(shí)
- 蟲號: 3860437
- 注冊: 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ā)表 | |
|---|---|---|---|---|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap－ 2026-03-22 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 3/150 |
|
|
|
[考研] 352求調(diào)劑 +3 | 大米飯！ 2026-03-22 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +5 | ￣^￣゜汗 2026-03-19 | 7/350 |
|
|
[考研] 求調(diào)劑 +6 | 十三加油 2026-03-21 | 6/300 |
|
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
|
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
|
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|