24小时热门版块排行榜    

查看: 520  |  回复: 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客户端
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 刘春杨 的主题更新
信息提示
请填处理意见