Emp6:
- ! ANSYS 5.5 Heat Transfer Seminar Notes ! Workshop #6
- ! Low Pressure Gas Turbine Case Thermal Analysis
- ! MUST HAVE case.db AND work6excel.txt TO RUN
- RESUME,case,db ! Supplied with seminar files
- /TITLE, Low Pressure Gas Turbine Case Thermal Analysis
- !* SET PREFERENCES TO THERMAL AND FLUID/PMETH,OFF
- KEYW,PR_SET,1 $KEYW,PR_THERM,1
- KEYW,PR_FLUID,1 $KEYW,PR_MULTI,1
- /COM,
- /COM,Preferences for GUI filtering have been set to display:
- /COM, Thermal $/COM, ANSYS Fluid
- /PREP7 !DEFINE THE SURF151 AND FLUID116 TYPES
- ET,2,SURF151 $ET,3,FLUID116
- KEYOPT,2,3,1 $KEYOPT,2,4,1
- KEYOPT,2,5,1 ! KEYOPTIONS FOR SURF151s
- KEYOPT,2,6,0 $KEYOPT,2,8,2
- KEYOPT,2,9,0
- KEYOPT,3,1,1 ! KEYOPTIONS FOR FLUID116s
- KEYOPT,3,2,1 $ETLIST, ALL
- MPLIST ! LIST PREDEFINED TEMP-DEPENDENT MPs
- ! REAL SET FOR FLUID116
- PI=ACOS(-1) $R,2
- R,3,,1E-6, ! Diameter, area, # channels
- !CREATE KEYPOINTS 112,113,114 ON WP
- K,,16.7628783879,145.253628512,0
- K,,16.4392239406,146.100950574,0
- K,,16.6407907696,147.310654059,0
- !CREATE LINES FOR FLUID ELEMENTS
- LSTR, 112, 113 $LSTR, 113, 114
- !SET ELEMENT ATTRIBUTES FOR THE LINES
- LSEL, , , ,91,92 $LATT,20,3,3,0, , ,
- !SET ELEMENT DIVISIONS TO ONE FOR
- !FLUID LINES
- LESIZE,ALL, , ,1,1
- !MESH LINES WITH FLUID116
- LMESH,ALL $CMLIST
- ! LIST PREDEFINED LINE COMPONENTS
- !CREATE SUR151S FOR ZONE 101
- CMSEL,S,LZONE101
- LPLOT $NSLL,S,1 $NPLOT
- TYPE,2 $MAT,1 $REAL,2
- ESURF,380 ! NODE 380 IS POINT A AND EXTRA NODE
- ESEL,R,ENAME,,151
- !CREATE A COMPONENT FOR SURF151 ON ZONE101
- CM,SURF101,ELEM
- ESEL,ALL $ESEL,U,TYPE,,2
- ! REPEAT FOR ZONE 102
- CMSEL,S,LZONE102
- LPLOT $NSLL,S,1 $NPLOT
- ESURF,380 $ESEL,R,TYPE,,2
- CM,SURF102,ELEM
- ESEL,ALL $ESEL,U,TYPE,,2
- !REPEAT FOR ZONE 103
- CMSEL,S,LZONE103
- LPLOT $NSLL,S,1 $NPLOT
- ESURF,381 ! NODE 381 IS AT POINT B
- ESEL,R,TYPE,,2 $CM,SURF103,ELEM
- ESEL,ALL $ESEL,U,TYPE,,2
- CMSEL,S,LZONE104
- NSLL,S,1 $ NPLOT
- !NSLL,Type,NKEY , Selects those nodes associated with the selected lines
- ESURF,382 ! NODE 382 IS AT POINT C
- ESEL,R,TYPE,,2 $CM,SURF104,ELEM
- ESEL,ALL $NSEL,ALL
- !DIMENSION THE TABLE FOR HF FOR ZONE104
- *DIM,HF104,TABLE,4,5,1,TIME,Y,
- !READ IN THE EXCEL SPREADSHEET
- *TREAD,HF104,work6excel,txt,,4,
- ALLSEL,ALL $LPLOT
- !SELECT THE LINES WHERE COUPLING IS NEEDED
- FLST,5,8,4,ORDE,8
- !FLST,NFIELD,NARG,TYPE,Otype,LENG
- !Specifies data required for a picking operation .
- FITEM,5,16 !FITEM,NFIELD,ITEM,ITEMY,ITEMZ—Identifies items chosen by a picking operation
- FITEM,5,29 $FITEM,5,41 $FITEM,5,57 $FITEM,5,60
- FITEM,5,80 $FITEM,5,88 $FITEM,5,-89
- LSEL,S, , ,P51X !CREATE COMPONENT OF LINES AT INTERFACE
- CM,LINTERFA,LINES $/REPLOT
- NSLL,S,1 $NPLOT
- !CHECK THE DISTANCE BETWEEN NODES TO
- !SET TOLERANCE FOR COUPLING
- NDISTANCE,136,15 !CREATE THE COUPLING EQUATIONS
- CPINTF,TEMP,0.005, $ALLSEL,ALL
- FINI $/SOLUTION $SOLCONTROL,ON,0
- !SELECT COMPONENT SURF101
- CMSEL,S,SURF101
- SFE,ALL,1,CONV, , %HF101% ! APPLY TABULAR Hf
- SFE,ALL,1,CONV,2, !SELECT COMPONENT SURF102
- CMSEL,S,SURF102
- SFE,ALL,1,CONV, , %HF102% ! APPLY TABULAR Hf
- SFE,ALL,1,CONV,2, !SELECT COMPONENT SURF103
- CMSEL,S,SURF103
- SFE,ALL,1,CONV, , %HF103% ! APPLY TABULAR Hf
- SFE,ALL,1,CONV,2, !SELECT COMPONENT SURF104
- CMSEL,S,SURF104
- SFE,ALL,1,CONV, , %HF104% ! APPLY TABULAR Hf
- SFE,ALL,1,CONV,2,
- !SELECT COMPONENT LZONE105 (LINES IN ZONE105)
- CMSEL,S,LZONE105 !APPLY TABULAR Hf AND Tb
- SFL,ALL,CONV, %HF105% , , %TB05T09%
- !SELECT COMPONENT LZONE106
- CMSEL,S,LZONE106 !APPLY TABULAR Hf AND Tb
- SFL,ALL,CONV, %HF106% , , %TB05T09%
- !SELECT COMPONENT LZONE107
- CMSEL,S,LZONE107 !APPLY TABULAR Hf AND Tb
- SFL,ALL,CONV, %HF107% , , %TB05T09%
- !SELECT COMPONENT LZONE108
- CMSEL,S,LZONE108 !APPLY TABULAR Hf AND Tb
- SFL,ALL,CONV, %HF108% , , %TB05T09%
- !SELECT COMPONENT LZONE109
- CMSEL,S,LZONE109 !APPLY TABULAR Hf AND Tb
- SFL,ALL,CONV, %HF109% , , %TB05T09%
- ALLSEL,ALL !DEFINE TABLE ARRAYS FOR INLET FLUID
- !TEMPERATURE AND MASS FLOW RATES
- !IN THE TWO FLUID116S
- *DIM,TINLET,TABLE,4,1,1,TIME, ,
- *DIM,FLOWAB,TABLE,4,1,1,TIME, ,
- *DIM,FLOWBC,TABLE,4,1,1,TIME, ,
- TINLET(1,0,1) = 0 $TINLET(1,1,1) = 304.33
- TINLET(2,0,1) = 3 $TINLET(2,1,1) = 870.33
- TINLET(3,0,1) = 160 $TINLET(3,1,1) = 870.33
- TINLET(4,0,1) = 163 $TINLET(4,1,1) = 304.33
- FLOWAB(1,0,1) = 0 $FLOWAB(1,1,1) = 0.001
- FLOWAB(2,0,1) = 3 $FLOWAB(2,1,1) = 0.0042
- FLOWAB(3,0,1) = 160 $FLOWAB(3,1,1) = 0.0042
- FLOWAB(4,0,1) = 163 $FLOWAB(4,1,1) = 0.001
- FLOWBC(1,0,1) = 0 $FLOWBC(1,1,1) = 0.000692
- FLOWBC(2,0,1) = 3 $FLOWBC(2,1,1) = 0.002908
- FLOWBC(3,0,1) = 160 $FLOWBC(3,1,1) = 0.002908
- FLOWBC(4,0,1) = 163 $FLOWBC(4,1,1) = 0.000692
- LPLOT $/AUTO, 1 $/REP
- !SPECIFY INLET FLUID TEMPERATURE
- !AT NODE 380 AS TABULAR INPUT
- D,380,TEMP,%TINLET%
- !SPECIFY MASS FLOW RATE ON FLUID116
- !FROM A TO B, USING HFLUX
- SFE,252,,HFLUX, ,%FLOWAB%
- SFELIS,252,HFLUX !SPECIFY MASS FLOW RATE FROM B TO C
- SFE,253,,HFLUX, ,%FLOWBC% $SFELIS,253,HFLUX
- ANTYPE,0 ! ANALYSIS TYPE STEADY-STATE
- TIME,1.E-6 ! TIME AT END OF STEADY-STATE
- NSUB,1 ! ONE SUBSTEP
- ALLSEL
- SOLVE ! RUN STEADY-STATE
- SAVE ! SAVE THE DATABASE FOR RESTART
- FINISH $/POST1
- PLNSOL,TEMP ! REVIEW THERMAL CONTOURS
- FINISH $RESUME $/SOLU
- ANTYPE,,REST $TIMINT,ON
- !DT_START = DX**2/(10 x ALPHA) = 0.05 seconds
- DELTIM,0.05,0.005,10,0
- KBC,0 ! RAMP LOADS
- TIME,163
- !DEFINE THE ARRAY TO USE FOR TIME STEP SIZE RESET
- *DIM,TSTEP,ARRAY,2,1,1, , ,
- TSTEP(1,1,1)=3 $TSTEP(2,1,1)=160
- TSRES,%TSTEP% ! USE TSTEP ARRAY
- OUTRES,NSOL,ALL $SAVE
- SOLVE $FINI
复制代码 |