| 查看: 557 | 回复: 0 | ||
[求助]
生死单元
|
|
我的本意是让他从下到上可是运行的时候他从上到下下面是命令流/Clear /FILENAME,WELD /TITLE,Ttemperaturre and stress analysis during welding /UNIT,SI /PREP7 ET,1,13,4 ET,2,13,4 MPTEMP,1,30,500,1000,1500,2000 MPDATA,EX,1,1,1.93E11,1.5E11,0.7E11,0.1E11,0.01E11 MPDATA,EX,2,1,1.02E11,0.5E11,0.08E11,0.001E11,0.0001E11 MPDATA,EX,3,1,1.17E11,0.9E11,0.3E11,0.05E11,0.005E11 TB,BKIN,1,5 TBTEMP,20,1 TBDATA,1,1.2E9,1.93E10 TBTEMP,500,2 TBDATA,1,0.933E9,1.5E10 TBTEMP,1000,3 TBDATA,1,0.435E9,0.7E10 TBTEMP,1500,4 TBDATA,1,0.07E9,0.1E10 TBTEMP,2000,5 TBDATA,1,0.007E9,0.01E10 TB,BKIN,2,5 TBTEMP,20,1 TBDATA,1,0.8E9,1.02E10 TBTEMP,500,2 TBDATA,1,0.4E9,0.5E10 TBTEMP,1000,3 TBDATA,1,0.07E9,0.08E10 TBTEMP,1500,4 TBDATA,1,0.0001E9,0.001E10 TBTEMP,2000,5 TBDATA,1,0.0001E9,0.00001E10 TB,BKIN,3,5 TBTEMP,20,1 TBDATA,1,0.9E9,1.17E10 TBTEMP,500,2 TBDATA,1,0.7E9,0.9E10 TBTEMP,1000,3 TBDATA,1,0.23E9,0.3E10 TBTEMP,1500,4 TBDATA,1,0.04E9,0.05E10 TBTEMP,2000,5 TBDATA,1,0.004E9,0.005E10 MP,DENS,1,7800 MP,DENS,2,4850 MP,DENS,3,8900 MP,ALPX,1,1.78E-5 MP,ALPX,2,9.36E-6 MP,ALPX,3,1.66E-5 MP,NUXY,1,0.29 MP,NUXY,2,0.30 MP,NUXY,3,0.29 MP,KXX,1,16.3 MP,KXX,2,7.44 MP,KXX,3,393 MP,C,1,502 MP,C,2,544 MP,C,3,385 K,1,0,0,0 K,2,0.5,0,0 K,3,1,0,0 K,4,0,0.26,0 K,5,0.35,0.26,0 K,6,0.65,0.26,0 K,7,1,0.26,0 LENGTH=DISTKP(2,5) K,8,0.5,LENGTH,0 L,1,2 L,2,5 L,5,4 L,4,1 L,2,3 L,3,7 L,7,6 L,6,2 LARC,6,5,8 AL,1,2,3,4 AL,2,8,9 AL,5,6,7,8 ESIZE,0.025 TYPE,2 MAT,2 MSHAPE,0,2D MSHKEY,1 AMESH,2 ESIZE,0.05 TYPE,1 MAT,1 AMESH,1 MAT,3 AMESH,3 AMESH,3 FINISH /SOLU ANTYPE,4 TRNOPT,FULL NSEL,S,LOC,X,0 D,ALL,UX D,ALL,UY NSEL,ALL NSEL,S,LOC,X,0 NSEL,A,LOC,X,1 D,ALL,TEMP,20 NNA=2 ESEL,ALL *GET,EMAX,ELEM,,NUM,MAX ASEL,S,,,NNA ESLA,S *GET,NSE,ELEM,,COUNT *DIM,NE,,NSE *DIM,NEX,,NSE *DIM,NEY,,NSE *DIM,NEORDER,,NSE MINE=0 *DO,II,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(II)=MINE *ENDDO MAX_TEMP=1500 DT1=1E-3 DT=5 T=0 ESEL,ALL /AUTO,1 *DO,I,1,NSE EKILL,NEORDER(I) ESEL,S,LIVE EPLOT *ENDDO ALLSEL,ALL OUTRES,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+DT1 TIME,T NSUBST,1 *DO,J,1,4 D,NELEM(NEORDER(I),J),TEMP,MAX_TEMP *ENDDO SOLVE T=T+DT1 TIME,T SOLVE *DO,J,1,4 DDELE,NELEM(NEORDER(I),J),TEMP *ENDDO T=T+DT-2*DT1 TIME,T NSUBST,NSUB1 SOLVE *ENDDO T=T+50000 TIME,T NSUBST,NSUB2 SOLVE FINISH /POST1 @wuming524 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
filecode
已经有11人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有33人回复
-
FileCode能看出啥？
已经有19人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有5人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有13人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
