sysweld命令結(jié)構(gòu)（課堂PPT）

1、1SysweldStructure of command files for the solver engine2F Loading a run file that contains the commandsF Interactive menuTwo possibilities to run commands: 3SEARch DATA 100DEFInition thermo-metallurgical analysisOPTIon THERmal METAllurgy PLANeRESTart GEOMetryMATErial properties .CONStrains.LOADs1 t

2、itle of loadingTABLesRETUrn4MATERIAL PROPERTIES COMMANDS Elements . / kx ky kz c ro material integ 2 Example of list elements : Elem 1 to 49 inte / Elem Groups $B1$ / Elem Groups $D1$ / Elem 1 to 100 step 5 / .5MATERIAL PROPERTIES COMMANDS Elements . / kx ky kz c ro material integ 2 Thermal conducti

3、vity KX KY KZ :1 ) E / KX KY KZ 0.04 constant E / KX KY KZ -1 =f(T)2) E / KX KY KZ 1 FT 2 =f(time)6MATERIAL PROPERTIES COMMANDS Elements . / kx ky kz c ro material integ 2 Specific heat capacity c :1 ) E / . c 500 . constant E / . c -2 . =f(T)2) E / . c 1 FT 2 =f(time)7MATERIAL PROPERTIES COMMANDS E

4、lements . / kx ky kz c ro material integ 2 Density ro :1 ) E / . ro 7800*-9 . constant E / . ro -3 . =f(T)2) E / . ro 1 FT 3 =f(time)Property depend on the phase and temperature : Elements . / material 1 integ 2Example of two materials : elem group $B1$ / mate 1 integ 2 elem group $B2$ / mate 2 inte

5、g 2 8MATERIAL PROPERTIES COMMANDS Elements . / kx ky kz c ro material integ 2 integ 2 for quadratic elements only+2D 2 x 2Example :E 1to49 inte / kx ky 1 ro 1 c 2 mate 1 inte 2 h19SEARch DATA 100DEFInition thermo-metallOPTIon THERmal METAllurgy PLANeRESTart GEOMetryMATErial properties .CONStrains.LO

6、ADs1 title of loadingTABLesRETUrn10CONSTRAINTS COMMAND Elements . / KT OR Nodes / TT Example of list elements : Elem 1 to 49 inte / Elem Groups $B1$ / Elem Groups $D1$ / Elem 1 to 100 step 5 / . Nodes Groups $B1$ / 11CONSTRAINTS COMMAND Elements . / KT Heat transfer coefficient KT : Elem / KT 20*-6

7、constant Elem / KT 1 variable 1 KT = f(T) Elem / KT 1 FT 1 KT = f(time) Elem / KT 1 variable 1 FT 2 KT = f(T,time)Example : KT = const.* f(T)* f(time)12CONSTRAINTS COMMAND Nodes . / TT Temperature at nodes TT : Nodes 1 to 100 / TT Nodes group $B1$ / TT13SEARch DATA 100DEFInition thermo-metallOPTIon

8、THERmal METAllurgy PLANeRESTart GEOMetryMATErial properties .CONStrains.LOADs1 title of loadingTABLesRETUrn14LOADS COMMAND Elem / TT QR ORNodes / TT QC . Example of elements list: Elem 1 to 49 inte / Elem Groups $B1$ / Nodes Groups $B1$ / Elem 1 to 100 step 5 / .15LOADS COMMAND Elem / TT Constrain E

9、 1 to 20 / KT 1000Loads1 example E 1 to 20 / TT 20Temperature of media : 1 ) Elem . / TT 20 constant 2) Elem . / TT 1 FT 2 TT = f(time)Example :16LOADS COMMAND Elem / QR Heat flux QR at elements W/m2 or W/m3: 1 ) Elem . / QR 20 constant 2) Elem . / QR 1 variable 2 QR = f(T) 3) Elem . / QR 1 FT 2 QR

10、= f(time) 4) Elem / QR 1 variable 1 FT 2 QR = f(T,time)Example : QR = const.* f(T)* f(time)17LOADS COMMAND Nodes / QC Heat flux QC at points W: 1 ) Nodes . / QC 20 constant 2) Nodes . / QC 1 variable 2 QC = f(T) 3) Nodes . / QC 1 FT 2 QC = f(time) 4) Nodes . / QC 1 variable 1 FT 2 QC = f(T,time)Exam

11、ple : QC = const.* f(T)* f(time)18LOADS COMMAND Nodes / TT Constraint n 1 to 20 / TTLoads1 example n 1 to 20 / TT 20Temperature of points TT : 1 ) Nodes . / TT 20 constant 2) Nodes . / TT 1 FT 2 TT = f(time)Examples :Constraint n 1 to 20 / TTLoads1 example n 1 to 20 / TT 1 ft 219SEARch DATA 100DEFIn

12、ition thermo-metallurgical analysisOPTIon THERmal METAllurgy PLANeRESTart GEOMetryMATErial properties .CONStrains.LOADs1 title of loadingTABLesRETUrn20TABLE COMMAND Tables Number of table / type of table ( pairs or specific values )Example of table No.1 : 1 / 1 ( 0.,0.03 ) (400,0.022) (900,0.027) Tv

13、alues21TABLE COMMAND Tables Number of table / type of table ( pairs or specific values )Example of FORTRAN table : 2 / FORTRAN FUNCTION F(T) E=0.8 T1=T+273.15 . F=D+T1 RETURN END22SEARch DATA 100DEFInition thermo-metallurgical analysisOPTIon THERmal METAllurgy PLANeRESTart GEOMetryMATErial propertie

14、s CONStrainsLOADs1 title of loadingTABLesRETUrnSAVE DATA 10123DATA 100 coordinates of nodes elements descriptionDATA 101 coordinates of nodes elements description material properties boundary conditions ( const.) loadings DATA-FILES ( *.TIT )24SEARch DATA 101TRANsient NON LINEar BEHAvior METAllurgy

15、4 file METALLURY.DATMETHod DIREct or ITERativeALGOrithm OPTI 1 ITER 100PRECision ABSOlute DISPlacement 0.5 FORCe 1INITial CONDition N / TT = 900 E / p = 0 0 0 1TIME INITial 0.RETUrnSAVE DATA TRANsient 102 25START OF CALCULATIONAPPLICATION MENU :FILE LOAD *.DATOR*ASSIgn 60 THERM.DAT*READ26BATCH MODE

16、:*ASSIgn 46 THERM.OUT*ASSIgn 60 THERM.DAT*READ 6027UNIX WINDOW :sysweld batch THERM.OUT THERM.DAT &28CCT DIAGRAM29LEBLONDS MODEL TPTPTfdtdPeq.P - proportion of phase t time T heating or cooling rateT temperaturePeq- proportion at phase equlibrium - delay time.30KOISTINEN- MARBURGER MODEL P - pro

17、portion of martensite T temperatureMs martensite at start temperatureb law parameter (0.011).(exp(1)(TMsbTP31EXTRACTION OF DATA FROM CCT DIAGRAMVELOCITY-21.4-9-3.8-3.4-1.7%1530354050FERRITEBAINITEVELOCITY-13.3-5.9-5.7-4.04-0.67%5378.684.610090TEMPERATURE 540 to 665 CTEMPERATURE 330 to 540 C32PHASeCA

18、LCulateVELOcity 0.5 1.7 3.4 -3.8 -9 -21.4 .FVELOcity .12 .4 .55 .PROPortion FT 1TAU FT 2TABLE1 / 1 540 1 665 02 / 1 540 18 665 18INITial TEMPerature 800 540 step 10DRAWFORMat33DRAWFORMat Ferritic transformationLIMIts xmin 540 xmax 800 ymin 0 ymax 1XHEAder temperatureYHEAder ferriteCURVeRETUrn 34META

19、LLURGY.DATmaterial 1 phase 4reaction 4 1 cooling k table 10 kp table 11 f table 12 fp table 12 4 2 cooling k table 20 kp table 21 f table 22 fp table 22 4 3 cooling k table 30 kp table 31 f table 32 fp table 32table 10 / 1 539 0 540 0.05555 665 0 11 / 1 539 0 540 0 665 0.05555 666 0 12 / 1 -1000 1 -

20、9 1 -3.8 0.5 -3.4 0.55 -1.7 0.4 -0.5 0.12 20 / 1 329 0 330 1 540 0 21 / 1 329 0 330 0 540 1 541 0 22 / 1 -1000 0.05 -50 0.05 -13.3 0.11 -5.9 0.13 -5.7 0.18 -4 1. -0.5 1. 30 / 1 0 0.011 430 0.011 431 0 31 / 1 0 0 430 0 32 / 1 -1000 1000 -0.5 0.5end k=Peq/ kp=(1- Peq)/ f=fp 35material 1 phase 4reactio

21、n 4 1 cooling Peq table 10 Tau table 11 f table 124 2 cooling Peq table 20 Tau table 21 f table 224 3 cooling MS 430 KM 0.011table 10 / 1 539 0 540 1 665 0 11 / 1 540 18 665 18 12 / 1 -1000 1 -9 1 -3.8 0.5 -3.4 0.55 -1.7 0.4 -0.5 0.12 20 / 1 329 0 330 1 540 0 21 / 1 330 1 540 1 22 / 1 -1000 0.05 -50

22、 0.05 -13.3 0.11 -5.9 0.13 -5.7 0.18 -4 1 -0.5 1. endMETALLURGY.DAT36MATERIAL 1 PHASE 4KX(1) = KX(2) = KX(3) = TABLE 1KX(4) = TABLE 2RO = TABLE 3ENTH(1) = ENTH(2) = ENTH(3) = TABLE 4ENTH(4) = TABLE 5REACTION 4 1 cooling k table 10 kp table 11 f table 12 fp table 12 4 2 cooling k table 20 kp table 21

23、 f table 22 fp table 22 4 3 cooling MS 400 KM 0.011TABLE 1 / 1 20.0 0.05 1000.0 0.03 2 / 1 20.0 0.015 1000.0 0.03 3 / 1 20.0 7.8*-6 1000.0 7.55*-6 4 / 2 0.0 0.0 600.0 1000.0 600000.0 600.0 5 / 2 0.0 70000.0 700.0 1000.0 770000.0 700.0 10 / 1 629.0 0.0 630.0 0.066667 780.0 0.0 11 / 1 629.0 0.0 630.0

24、0.0 780.0 0.066667 781.0 0.0 12 / 1 -10.00 0.01 -6.500 0.22 -3.000 0.23 -1.435 0.27 -0.77 0.27 -0.165 0.13 -0.116 1.0 20 / 1 399.0 0.0 400.0 1.0 630.0 0.0 21 / 1 399.0 0.0 400.0 0.0 630.0 1.0 631.0 0.0 22 / 1 -100.0 0.01 -40.0 0.08 -22.8 0.7 -15.4 0.8 -11.6 1.0 -5.54 1.5 -1.0 1.037CCT DIAGRAM PLOTTI

25、NGMETALLURGYCALCULATE CCT diagramFILE “METALLURGY.DAT”MATERIAL 1AUSTENITE 4CCT VMIN 1000 VMAX 0.1 CUTTING 20 LOGARITHMINITIAL TEMPERATURE 800 20 step 20DRAWABSCISA TIME SECONDSPRECISION 0.01DIAGRAM CYCLERETURN38MECHANICAL ANALYSIS39SEARch DATA 101DEFInition mechanical analysisOPTIon TWO-DIME TRANsla

26、tion THERMoelas.RESTart GEOMetryMATErial PROPerties E 1 to 49 inte / e 1 nu 2 lx ly lz 10 phase 4* yield 20 slope 30 model 3 integ 2CONStrain n 1000 to 1012 2 / uy n 4000 to 4012 4 / ux n 1 / ux uyLOADs1 quenching40 TABLE 1 / 1 0 200*3 900 120*3 2 / 1 0 0.29 900 0.30 10 / -12 12 12 11 4 phases 11 /

27、1 20 0.01024 900 0.01 12 / 1 20 0. 900 0.0132E 1 to 49 inte / e 1 nu 2 lx ly lz 10 phase 4 -41 TABLE 20 / -22 -23 -24 -21 phase 4 21 / 1 20 220 200 200 300 180 ; austenite 22 / 1 20 320 200 290 300 280 ; ferrite 23 / 1 20 460 200 420 300 400 ; bainite 24 / 1 20 660 200 620 300 550 ; martensite* yield 20