WS19-1 WORKSHOP 19 NASTRAN BULK DATA ENTRIES NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation

WS19-2 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation

WS19-3 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation n Problem Description u A cooling fin is subject to heat flux and convection cooling. We will add the CHBDYP card along with CONV,PCONV,MAT4 to model convection.

WS19-4 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation n CHBDYP,LINE, BDYOR,PHBDY,CONV,PCONV,QBDY3,MAT4 n n The following is the NASTRAN input test deck: n SOL 153 n CEND n TITLE=2D PLATE Thermal analysis n Subtitle=problem 1 n analysis=heat n spc=10 n temp(init)=5 n load=20 n NLPARM=30 n thermal=all n flux=all n spcf=all n BEGIN BULK n NLPARM,30 n include 'plate.txt' n $Apply convection boundary condition on the left edge n $ We need to define CHBDYP,LINE elements n $CHBDYP,EID, PID, TYPE, IVIEWF, IVIEWB, G1, G2, G0, + n $+,RADMIDF, RADMIDB, GMID, CE, E1, E2, E3 n CHBDYP,101,1,LINE,,,1,6 n =,*1,=,=,=,=,*5,*5 n =7 n $BDYOR,TYPE,IVIEWF,IVIEWB,RADMIDF,RADMIDB,,PID,GO,+ n $,CE,E1,E2,E3 n BDYOR,LINE,,,,,,1,,+ n +,,-1.0,0.0,0.0 n $PHBDY PID AF D1 D2 n PHBDY,1,1.0 n $CONV,EID, PCONID, FLMND, CNTLND, TA1, TA2, TA3, TA4, +BC n $+BC, TA5, TA6, TA7, TA8 n CONV,101,2,,,99 n =,*1,== n =7 n $PCONV,PCONID, MID, FORM, EXPF n PCONV,2,20,0,0.0 n $MAT4,ID,K, CP, RHO, H n MAT4,20,,,,0.123 n $DEFINE AMBIENT POINT 99 n SPOINT,99 n SPC,10,99,1,70.0 n TEMP,5,99,70.0 n TEMPD,5,200.0 n $APPLY HEAT FLUX BOUNDARY CONDITION ALONG THE RIGHT EDGE n $CHBDYP,EID, PID, TYPE, IVIEWF, IVIEWB, G1, G2, G0, + n $+,RADMIDF, RADMIDB, GMID, CE, E1, E2, E3 n CHBDYP,201,1,LINE,,,55,60 n =,*1,=,=,=,=,*5,*5 n =1 n $QBDY3,SID, Q0, CNTRLND, EID1, EID2, EID3, EID4, EID5 n QBDY3,20,0.05,,201,202,203 n PARAM,POST,0 n $Thermal material properties for the plate n $K,CP,Rho n MAT4,1,6.-4,.146,.283 n ENDDATA Problem 1: Test Deck: input1.dat

WS19-5 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation Temperature Contour

WS19-6 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation 0 LOAD STEP = E+00 T E M P E R A T U R E V E C T O R POINT ID. TYPE ID VALUE ID+1 VALUE ID+2 VALUE ID+3 VALUE ID+4 VALUE ID+5 VALUE 1 S E E E E E E+01 7 S E E E E E E S E E E E E E S E E E E E E S E E E E E E S E E E E E E S E E E E E E S E E E E E E S E E S E E E E S E E E E S E E E E S E E E E S E+01

WS19-7 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation Problem 2: have thermal conductivity varies as a function of temp Test deck is input2. dat Make the thermal conductivity as a function of temperature Temp(F)Thermal conductivity(K) times 6.E mat4,1,6.-4,.146,.283 matt4,1,10 tablem2,10,0.0,,,,,,,+ +,-10.0,1.0,0.0,1.0,100.0,1.20,500.0,1.30, +,endt

WS19-8 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation T E M P E R A T U R E V E C T O R POINT ID. TYPE ID VALUE ID+1 VALUE ID+2 VALUE ID+3 VALUE ID+4 VALUE ID+5 VALUE 1 S E E E E E E+01 7 S E E E E E E S E E E E E E S E E E E E E S E E E E E E S E E E E E E S E E E E E E S E E E E E E S E E S E E E E S E E E E S E E E E S E E E E S E+01 Result from running input2. dat where k is a function of temperature

WS19-9 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation Problem 3: The right hand edge of the model is exposed to a hot source temperature at 1500 degree F. The view factor is assumed to be 1.0, and the emissivity of the material is 0.8. The left hand edge is convected to an ambient temperature at 70 F with H= Find the temperature distribution. Test deck name: input3a.dat NASTRAN bulk data entries: CHBDYP,RADBC,RADM,PARAM,SIGMA,PARAM,TABS,BDYOR,PHBDY,CONV,PONV,MAT4 Problem 3: Add radiation boundary conditions to the model

WS19-10 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation SOL 153 CEND TITLE=2D PLATE Thermal analysis Subtitle=problem 1 analysis=heat spc=10 temp(init)=5 $load=20 NLPARM=30 thermal=all flux=all spcf=all BEGIN BULK $SIGMA IS BTU/SEC/IN2*R**4 PARAM,SIGMA, PARAM,TABS,460.0 NLPARM,30 include 'plate.txt' CHBDYP,101,1,LINE,,,1,6 =,*1,=,=,=,=,*5,*5 =7 $BDYOR,TYPE,IVIEWF,IVIEWB,RADMIDF,RADMIDB,,PID,GO,+ $,CE,E1,E2,E3 BDYOR,LINE,,,,,,1,,+ +,,-1.0,0.0,0.0 $PHBDY PID AF D1 D2 PHBDY,1,1.0 $ $CONV EID PCONID FLMND CNTLND TA1 TA2 TA3 TA4 +BC $+BC TA5 TA6 TA7 TA8 CONV,101,2,,,99 =,*1,== =7 $PCONV PCONID MID FORM EXPF PCONV,2,20,0,0.0 $MAT4 ID K CP RHO H MAT4,20,,,,0.123 $DEFINE AMBIENT POINT 99 SPOINT,99 SPC,10,99,1,70.0 TEMP,5,99,70.0 TEMPD,5, Continue------

WS19-11 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation $APPLY RADIATION HEATING ALONG THE RIGHT EDGE $CHBDYP EID PID TYPE IVIEWF IVIEWB G1 G2 G0 +CHBDYP $+CHBDYPRADMIDF RADMIDB CE E1 E2 E3 CHBDYP,201,2,LINE,,,55,60,,+c1 =,*1,=,=,=,=,*5,*5,=,*1 =1 +c1,1,,,,1.0,0.0,0.0 *1,== =1 phbdy,2,1.0 radm,1,0.8,0.8 $RADBC NODAMB FAMB CNTRLND EID1 EID2 EID3 radbc,101,1.0,,201,thru,203 $ $DEFINE AMBIENT POINT 101 SPOINT,101 SPC,10,101,1, TEMP,5,101, MAT4,1,6.-4,.146,.283 PARAM,POST,0 ENDDATA

WS19-12 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation REPLICATION COMMANDS FOR DUPLICATING THE CONTINUATION LINES: CHBDYP,201,2,LINE,,,55,60,,+c1 =,*1,=,=,=,=,*5,*5,=,*1 =1 +c1,1,,,,1.0,0.0,0.0 *1,== =1 CHBDYP LINE C1 +C CHBDYP LINE C2 +C CHBDYP LINE C3 +C

WS19-13 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation LOAD STEP = E+00 T E M P E R A T U R E V E C T O R POINT ID. TYPE ID VALUE ID+1 VALUE ID+2 VALUE ID+3 VALUE ID+4 VALUE ID+5 VALUE 1 S E E E E E E+01 7 S E E E E E E S E E E E E E S E E E E E E S E E E E E E S E E E E E E S E E E E E E S E E E E E E S E E S E E E E S E E E E S E E E E S E E E E S E S E+03 We can see that the maximum temperature is on grid 60.

WS19-14 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation Problem 4: Transient Thermal analysis. We have a constant heat flux 0.05 imposed on the right hand edge from time zero to 500 sec. On the left hand edge we have convection boundary conditions to a time-varying ambient temperature with h= Initially at time equal to zero, the temperature of the plate is at 70 degree F. Ambient TemperatureTime(sec)Temp(F) TABLED2,100,0.0,,,,,,,+ +,0.0,70.0,20.0,78.0,50.0,83.0,100.0,83.0,+ +,ENDT NASTRAN bulk data entries:DLOAD,TLOAD1,TLOAD2,TABLED2,(TEMPBC,TRAN) Problem 4: Transient thermal analysis

WS19-15 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation SOL 159 CEND TITLE=2D PLATE Transient Thermal analysis Subtitle=problem 4 analysis=heat IC=5 DLOAD=200 TSTEPNL=500 SET 77=60,99,21,25 thermal=77 $flux=all spcf=all OUTPUT(XYPLOT) XGRID=YES YGRID=YES XTITLE=SEC YTITLE=DEGREE F XYPLOT TEMP/60(T1),99(T1),21(T1),25(T1) BEGIN BULK $TSTEPNLID NDT DT NO METHOD KSTEP MAXITER CONV $TSTEPNLEPSU EPSP EPSW MAXDIV MAXQN MAXLS TSTEPNL,500,200,1.0,1,ADAPT,,,U,+ +,.01 $DLOAD SID S S1 L1 S2 L2 S3 L3 +ABC $+DLOAD S4 L4 ETC DLOAD,200,1.0,1.0,20,1.0,40 $LOADSET 20 REPRESENT HEAT FLUX IS A FUNCTION OF TIME $LOADSET 40 REPRESENT THE AMBIENT TEMPERATURE 99 IS A FUNCTION OF TIME include 'plate.txt' CHBDYP,101,1,LINE,,,1,6 =,*1,=,=,=,=,*5,*5 =7 BDYOR,LINE,,,,,,1,,+ +,,-1.0,0.0,0.0 $PHBDY PID AF D1 D2 PHBDY,1,1.0 $ CONV,101,2,,,99 =,*1,== =7 PCONV,2,20,0,0.0 MAT4,20,,,,0.123 $DEFINE AMBIENT POINT 99 $TABLED2ID X1 +ABC $TABLED2X1 Y1 X2 Y2 X3 Y3S X4 Y4 +DEF $TABLED2X5 Y5 X6 Y6 X7 Y7 X8 Y9 +GHI TABLED2,100,0.0,,,,,,,+ +,0.0,70.0,20.0,78.0,50.0,83.0,100.0,83.0,+ +,ENDT

WS19-16 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation $TEMPBC SID TYPE TEMP1 GID1 TEMP2 GID2 TEMPP3 GID3 TEMPBC,400,TRAN,1.0,99 $TLOAD1 SID L M TYPE TF TLOAD1,40,400,,,100 SPOINT,99 TEMP,5,99,70.0 $INITIAL TEMP AT TIME EQUAL TO ZERO TEMPD,5,70.0 $APPLY HEAT FLUX BOUNDARY CONDITION ALONG THE RIGHT EDGE $CHBDYP EID PID TYPE IVIEWF IVIEWB G1 G2 G0 +CHBDYP $+CHBDYPRADMIDF RADMIDB CE E1 E2 E3 CHBDYP,201,1,LINE,,,55,60 =,*1,=,=,=,=,*5,*5 =1 $QBDY3 SID Q0 CNTRLND EID1 EID2 EID3 EID4 EID5 QBDY3,200,0.05,,201,202,203 $DEFINE UNIT STEP FUNCTION OF 1.0 FOR ALL TIME BETWEEN ZERO TO 500 SEC $AFTER 500 SEC, THE LOADING IS EQUAL TO ZERO $TLOAD2 SID DAREA DELAY TYPE T1 T2 F P $+TLOAD2C B TLOAD2,20,200,,,0.0,500.0 $ MAT4,1,6.-4,.146,.283 PARAM,POST,0 ENDDATA

WS19-17 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation Problem 5: Add a time-varying convective coefficient H(time)Time(sec)H (convection coefficient) TABLED2,200,0.0,,,,,,,+ +,0.0,0.123,20.0,0.150,50.0,0.156,100.0,0.156,+ +,ENDT

WS19-18 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation We introduce an additional SPOINT to represent the convection as a function of time. DLOAD,200,1.0,1.0,20,1.0,40,1.0,60 $CONV EID PCONID FLMND CNTLND TA1 TA2 TA3 TA4 +BC $+BC TA5 TA6 TA7 TA8 CONV,101,2,,201,99 =,*1,== =7 $PCONV PCONID MID FORM EXPF PCONV,2,20,0,0.0 MAT4,20,,,,1.0 spoint,201 TABLED2,200,0.0,,,,,,,+ +,0.0,0.123,20.0,0.150,50.0,0.156,100.0,0.156,+ +,ENDT $TEMPBC SID TYPE TEMP1 GID1 TEMP2 GID2 TEMPP3 GID3 TEMPBC,600,TRAN,1.0,201 TEMP,5,201,0.123 $TLOAD1 SID L M TYPE TF TLOAD1,60,600,,,200

WS19-19 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation 0 POINT-ID = 60 T E M P E R A T U R E V E C T O R TIME TYPE VALUE 0.0 S E E+00 S E E+00 S E E+00 S E E+00 S E E+00 S E E+01 S E E+01 S E E+01 S E E+01 S E E+01 S E E+01 S E E+01 S E E+01 S E E+01 S E E+02 S E E+02 S E E+02 S E E+02 S E E+02 S E E+02 S E E+02 S E+02 POINT-ID = 201 T E M P E R A T U R E V E C T O R TIME TYPE VALUE 0.0 S E E+00 S E E+00 S E E+00 S E E+00 S E E+00 S E E+01 S E E+01 S E E+01 S E E+01 S E E+01 S E E+01 S E E+01 S E E+01 S E E+01 S E E+02 S E E+02 S E E+02 S E E+02 S E E+02 S E E+02 S E E+02 S E-01

WS19-20 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation Problem 6: Add fluid nodes and apply forced convection along the left hand edge coupling with the structure nodes Forced convection with h equal to 0.01 Btu/sec/inch**2. F and mass flow rate is equal to lbm/sec Volumetric heating on the whole structure with heating equal to 0.5 Btu/sec/inch**3 Test deck is input6. dat NASTRAN bulk data entries: QVOL,(CHBDYP,FTUBE),CONVM,PCONVM,MAT4 Problem 6: Add advection(fluid flow) in the model

WS19-21 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation SOL 153 CEND TITLE=2D PLATE Thermal analysis,forced convection Subtitle=problem 6 analysis=heat spc=10 temp(init)=5 load=20 NLPARM=30 thermal=all flux=all spcf=all BEGIN BULK $Add one dimensional fluid nodes GRID,301,,-0.2,0.0,0.0 =,*1,=,=,* ,== =8 $CHBDYP EID PID TYPE IVIEWF IVIEWB G1 G2 G0 +CHBDYP $+CHBDYPRADMIDF RADMIDB CE E1 E2 E3 CHBDYP,401,120,FTUBE,,,301,302 =,*1,=,=,=,=,*1,*1 =7 $Diameter of the duct is 1.0inch PHBDY,120,,1.0 $CONVM EID PCONID FLMND CNTMDOT TA1 TA2 CONVM,401,130,,777,1,6 =,*1,=,=,=,*5,*5 =7 $PCONVM PCONID MID FORM FLAG COEF EXPR EXPPI EXPPO PCONVM,130,220,0,1,0.01 $MAT4 ID K CP RHO H $PROPERTY OF AIR MAT4,220, ,0.997,

WS19-22 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation SPOINT,777 SPC,10,777,1,0.005 TEMP,5,777,0.005 $INLET AIR TEMPERATURE SPC,10,301,1,70.0 TEMP,5,301,70.0 NLPARM,30 include 'plate.txt' TEMPD,5,200.0 $APPLY VOLUMETRIC HEATING ON THE SURFACES $QVOL SID QV CNTRLND EID1 EID2 EID3 EID4 EID5 $QVOL SID QV CNTRLND EID1 "THRU" EID2 QVOL,20,0.50,,1,thru,48 MAT4,1,6.-4,.146,.283 PARAM,POST,0 ENDDATA

WS19-23 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation POINT ID. TYPE ID VALUE ID+1 VALUE ID+2 VALUE ID+3 VALUE ID+4 VALUE ID+5VALUE 1 S E E E E E E+01 7 S E E E E E E S E E E E E E S E E E E E E S E E E E E E S E E E E E E S E E E E E E S E E E E E E S E E S E E E E S E E E E S E E E E S E E E E S E E E E E E S E E E E S E-03 The fluid nodes are from grid 301 to 310, the inlet air temperature is 70 degree F, and the exit air temperature is degree F.

WS19-24 NAS104, Workshop 19, March 2004 Copyright 2004 MSC.Software Corporation