 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 6004 | 回復(fù): 2 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
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號二維耦合單元，同時具有溫度和位移自由度 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坐標排序 !模擬涂層由下向上“生長”的過程 *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ù)此樓» 猜你喜歡
11408總分309，一志愿東南大學(xué)求調(diào)劑，不挑專業(yè)
已經(jīng)有6人回復(fù)
-
求調(diào)劑，一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分
已經(jīng)有8人回復(fù)
-
環(huán)境工程 085701，267求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料工程專碩求調(diào)劑
已經(jīng)有8人回復(fù)
-
297求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料科學(xué)與工程求調(diào)劑
已經(jīng)有7人回復(fù)
-
281求調(diào)劑
已經(jīng)有7人回復(fù)
-
288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行
已經(jīng)有8人回復(fù)
-
求收留
已經(jīng)有6人回復(fù)
-
085600 295分求調(diào)劑
已經(jīng)有9人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 大神們，求幫忙啊，謝謝
已經(jīng)有3人回復(fù)
- 新手求助高斯輸入文件寫法，求大神幫忙
已經(jīng)有8人回復(fù)
- 求兩個值中最大的值但編寫的代碼不知道哪里錯了，求大神幫忙。
已經(jīng)有10人回復(fù)
- 求助fortran大神，我是初學(xué)者，幫我分析一下下面的代碼
已經(jīng)有5人回復(fù)
- 求大神幫忙設(shè)計一個工裝夾具
已經(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)于鋁合金的抗彎強度的問題！�。。。。�！
已經(jīng)有13人回復(fù)
- 求各位大神幫忙 檢測豬肝中銅鐵鋅鉛的前處理方法有哪些�。�！
已經(jīng)有8人回復(fù)
- 哪位大神能幫忙總結(jié)下常見天然產(chǎn)物的預(yù)檢識方法啊，膜拜~~~~
已經(jīng)有3人回復(fù)
- 初學(xué)者請各位高手幫忙看一下GSAS精修的結(jié)果，請指點
已經(jīng)有14人回復(fù)
3樓2018-04-09 09:29:43
張大俠要學(xué)習(xí)
新蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 114.5
- 紅花: 4
- 帖子: 20
- 在線: 9.6小時
- 蟲號: 3860437
- 注冊: 2015-05-10
2樓2015-05-10 16:11:24
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行 +8 | lllllos 2026-03-30 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-30 | 7/350 |
|
|
[考研] 一志愿華中師范化學(xué)332分求調(diào)劑 +3 | Lyy930824@ 2026-03-29 | 3/150 |
|
|
|
[考研] 0703一志愿9，初試成績：338，四六級已過，有科研經(jīng)歷，求調(diào)劑！ +7 | Zuhui0306 2026-03-25 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +10 | 家佳佳佳佳佳 2026-03-29 | 10/500 |
|
|
|
[考研] 一志愿廈門大學(xué)材料工程專碩354找調(diào)劑！�。� +5 | 貝唄鋇鋇 2026-03-30 | 5/250 |
|
|
|
[考研] 求調(diào)劑323材料與化工 +10 | 1124361 2026-03-24 | 10/500 |
|
|
|
[考研] 考研調(diào)劑 +5 | Sanmu-124 2026-03-26 | 5/250 |
|
|
|
[考研] 356求調(diào)劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|
|
|
[考研] 數(shù)一英一271專碩（085401）求調(diào)劑，可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4 | 小張zxy張 2026-03-26 | 8/400 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
|
[考研] 320分，材料與化工專業(yè)，求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4 | 手機用戶 2026-03-26 | 4/200 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
|
[考研] 一志愿鄭大085600，310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
|
[考研] 0703化學(xué)338求調(diào)劑！ +6 | Zuhui0306 2026-03-26 | 7/350 |
|