'" ~. f_ \ 

-' I 

NI\S/\ 

NatIOnal Aeronautics ;mrl 

Space AdministratIOn 

NASA_CR-174891 

19860015460 

BURNER LINER 

THERMAL/STRUCTURAL 

LOAD MODELING 

TRANCITS PROGRAM USER'S MANUAL 

By 

R. Maffeo 

August 1985 

General Electric Company 

Aircraft Engine Business Group 

Cincinnati, Ohio 45215 

Perpared For 

National Aeronautics and Space Administration 

Lewis Research Center 

21000 Brookpark Road 

Cleveland, Ohio 44135 

NASA-Lewis Research Center 

Contract NAS3-23272 

1111111111111 IIII 11111 IIIII 11111 III11 IIII lUI 

NF01212 

I

r {U).2 r n" I"' rq)y . I • I .. • ~ • t .. f , 

SEP 2 31985 

I 'I ~'l· l -. ... • '.~ . '. '. ,.f·t " lUJTER 

!.ILR/\RY, rlASA 

I WPfON, VIRCIr.jIA 
1 Report No I 2 

Government Accession No 3 Recipient's Catalog No 

4 Tille and Subtitle 5 Report Date 

August 1985 

Burner Llner Thermal/Structural Load Model1ng - 6 Performing Organization Code TRANCrrS Program Userls Manual 

7 Author(s) 8 Performing Organization Report No 

R. Maffeo R85AEB452 

10 Work Unit No 

9 Performing Organization Name and Address 

General Electrlc Company 11 Contract or Grant No 

Alrcraft Englne Buslness Group 

Clnclnnatl, NAS3-23272 OhlO 45215 

13 Type of Report and Period Covered 

12 Sponsoring Agency Name ."d Address Software User's Manual Natlonal Aeronautlcs and Space Admlnistratl0n 

Sponsoring AgencV Code Lewls Research Center I 14 , Cleveland, OR 44135 , 

, 

15 Supplementarv Notes 

I Project Manager, Dr. R. L. Thompson 

Structures and Mechanlcal Technologles Dlvlslon I 

NASA Lewls Research Center, Cleveland, OR 44135 I 

16 Abstract , 

, The software package developed under thlS contract was called TRansfer ANalysls Code to I 

Interface !hermal/~tructural Problems (TRANCITS). The TRANCITS code-satlsfles all the obJectlves 

for transferrlng thermal data between heat transfer and structural models of combustor llners; ln I 

addltion, lt can be used as a generlc thermal translator between heat transfer and stress models 

of any component, regardless of the geometry. TRANCITS can accurately and efflciently convert the : 

temperature dlstrlbutlons predlcted by the heat transfer programs to those requlred by the stress 

codes. It can be used for both llnear and nonllnear structural codes, and lt can produce nodal 

temperatures, elemental centrold temperatures, or elemental Gauss pOlnt temperatures. The thermal 

output of both the MARC and SINDA heat transfer codes can be lnterfaced dlrectly wlth TRANCITS, 

and lt wlll automatlcally produce stress model codes formatted for NASTRAN and MARC. In adda lon 

to these codes, any thermal program and structural program can be lnterfaced by uSlng the neutral 

l.nput and output forms supported by TRANCITS. 

In summary, the TRANCITS code can be used to lnterface temperature data between thermal and 

structural analytlcal models. The use of thlS transfer module allows the heat transfer analyst 

to select the thermal mesh denslty and thermal analysls code best sUlted to solve the thermal 

problem, and lt glves the same freedoms to the stress analyst wlthout the efflclency penaltles assoclated wlth common meshes and the accuracy penaltles assoclated wlth the manual transfer of 

thermal data. 

17 Key Words (Suggested by Author(s)) 18 D,strtbutlon Statement 

Interface Flnlte Dlfference 

Reat Transfer Analysls Transfer Module 

Structural Analysls Reat Transfer 

Flnlte Element Stress Analysls 

Isoparametrlc Computer-Alded 

Enlllneerln1z 

19 SecUrtty Oasslf (of thiS report) 20 Secufl tv Classll (of thiS page) I 21 No 01 Pages 22 ;>,.ce· 

I Unc lass 1 £led Unc lass 1 £led j 47 

• For sale by the National Technical Information Service Spllnefleld Vuglnla 22161 

NASA-C·168 (Rev 1 ()'75) 
This Page Intentionally left Blank 
Sectl.on 

1.0 

2.0 

3.0 

4.0 

5.0 

6.0 

TABLE OF CONTENTS 

INTRODUCTION 

GENERAL OVERVIEW 

CURRENT LIMITATIONS 

INPUT DESCRIPTION 

4.1 Neutral Heat Transfer Input 

4.2 SINDA Temperature Input 

4.3 SINDA Geometry Input 

4.4 MARC Posttape Input 

4.5 Stress Model Coordl.nate Input 

4.6 Stress Model Element Input 

4.7 User Control Fl.le Input 

OUTPUT DESCRIPTION 

5.1 Neutral Temperature Output Fl.le 

5.2 NASTRAN Temperature File 

5.3 MARC Temperature File 

5.4 Hard Copy Output Description 

RUN COMMANDS 

APPENDIX A - CONTENTS OF MARC POSTTAPE 

APPENDIX B - SAMPLES OF FILE FORMATS 

APPENDIX C - SINDA FORTRAN OUTPUT SUBROUTINES 

hi 

1 

2 

4 

5 

5 

10 

11 

12 

12 

13 

14 

17 

17 

18 

18 

18 

21 

22 

31 

41 
1.0 INTRODUCTION 

In aircraft turb1ne engine hot section components, cyclic thermal 

stresses are the most 1mportant damage mechan1sms. Consequently, an accurate 

and reliable pred1ction of thermal loads is essent1al to improved durab1l1ty. 

To ach1eve this goal, considerable effort has been expended over the past 20 

years to (1) acqu1re eng1ne temperature test data, and (2) develop accurate, 

relIable, efficIent computer codes to predict steady-state and transient temperatures and for calculating elast1c and 1nelast1c cyclic stresses and strains 

in hot sect10n components. There 1S a need for continued development of these 

codes because the avaIlabIlIty of more accurate analysis technIques for complex 

configurations has enabled engine desIgners to use more sophIsticated deSIgns 

to achieve higher cycle effICIency and reduce weight. 

It became apparent 1n recent years that a serious problem eXIsted when 

InterfaCIng output temperatures and temperature gradients from the heat 

transfer codes with the 1nput to the stress analYSIS codes. Due to the growth 

1n computer capacity and speed, and with the development of input preprocessors and output postprocessors, the analYSIS of components USIng hundreds or 

thousands of nodes In the heat transfer and stress models became economIcal 

and routIne. ThIS has aggravated the problem of manually transferrIng output 

temperatures from heat transfer nodes to stress analysis Input to the pOInt 

where the engineering effort required is comparable to that required for the 

remainder of the analysis. Furthermore, a considerable amount of approx1mation has been introduced in an effort to accelerate the process. ThIs tends 

to introduce errors into the temperature data, thereby negatIng the Improved 

accuracy in the temperature distribut10n achIeved through use of a fIner 

mesh. 

The objective of the Burner LIner/Thermal Structural Load Modeling program was to develop a transfer module that would accurately and effiCIently 

transfer three-dimenSIonal temperature information from heat transfer analysis 

codes to stress analysis codes. ThIS objective has been satIsfied by the 

development of the computer code TRANCITS (TRansfer ANalYSIS Code to Interface !hermal/~tructural Problems). -- -- - -

ThIS report is a prelIm1nary user's manual for the TRANCITS program. 

1 
2.0 GENERAL OVERVIEW 

TRANCITS is a computer code developed by the A1rcraft Engine Business 

Group of the General Electric Co. in conjunct10n w1th the NASA-Lew1s Research 

Center to address the problem of transferring three-dimensional temperature 

information from numer1cal heat transfer analyses to numer1cal stress 

analyses. The advanced features of TRANCITS include: 

• The ab1l1ty to transfer temperatures between dissim1lar mesh 

densit1es 

• The ability to handle fin1te d1fference and fin1te element heat 

transfer codes 

• Extensive user controls on the interfac1ng procedure 

• Automat1c interface to SINDA, MARC, and NASTRAN 

• Internal checks on the accuracy of the transfer 

• The ab1lity to handle stress points slightly outs1de the heat 

transfer model 

• Accurate and efficient transfer of temperatures 

• Modular program arch1tecture to fac1litate maintenance and 

enhancements. 

Figure 1 1S a schematic of how TRANCITS 1S used in the thermal interfac1ng 

procedure. 

It should be noted that although the interfaces between SINDA, MARC, and 

NASTRAN codes are Hard Wired into the TRANCITS system, the Neutral Input and 

Neutral Output f1les allow TRANCITS to be used with any comb1nat10n of heat 

transfer codes and stress analys1s codes. 

2 
t.J 

I SINDA 

I MARC -1 TEMP/GEOM :> 

Any 

HIT 

Code 

I T/G 

T 

R 

A 

N 

C 

I 

T 

S 

I TEMPS ~ NASTRAN 

I TEMPS ~ MARC~ 

I TEMPS 

Any 

Stress 

Code 

Figure 1. TRANCITS Schematic. 
3.0 CURRENT LIMITATIONS 

The follow1ng are the current lim1tations for the first release of 

TRANCITS (Version 1.0). 

./( 

l. Heat Transfer Model Lim1tat10ns 

1~ 

• Maximum number of heat transfer l1m1tations = 5000 

• Maximum heat transfer element name = 5000 

• Maximum number of tranS1ent solutions = 500 

• Maximum number of heat transfer nodes = 5000 

• Maximum heat transfer node name = 5000 

• Heat transfer element l1brary - 8-noded 

so11d element 

2. Stress Analys1s Model 

• Max1mum number of stress nodes = 5000 

• Maximum stress node name = 5000 

• Maximum number of stress elements = Un11m1ted 

• Maximum stress element name Un11m1ted 

• Stress element library - any three-d1mens1onal 

solid element 

All of the numer1cal l1mits can be easily increased by increas1ng the amount 

of temporary storage allocated to the program. 

4 
4.0 DESCRIPTION OF INPUT 

TRANCITS can accept many different input forms, 1ncluding the output of 

the SINDA and MARC programs d1rectly. The choice of the input form depends 

on which heat transfer code was used to produce the temperatures. For heat 

transfer codes other than MARC and SINDA, a neutral heat transfer 1nput f1le 

can be used to input the requ1red data into TRANCITS. 

Both f1nite d1fference and finite element codes can be used, but 1n 

e1ther case the geometry of the heat transfer mode and the computed 

temperatures from the heat transfer analys1s must be 1nput. In addition to 

this information, a user control file must be entered. This f1le contains 

data that controls the interface procedure. All of these input files are 

descr1bed in the following sections. 

4.1 NEUTRAL HEAT TRANSFER INPUT FILE 

Th1s f11e is either input or created internally. It contains all the 

geometry and temperature information for the heat transfer model. All 

processing of this 1nformation 1n TRANCITS 1S done using th1s f1le. If the 

heat transfer code used was SINDA or MARC, the output of these two programs 

is 1nput into TRANCITS and the Neutral Heat Transfer Input file is created 

internally. On the other hand, if the heat transfer code used was some other 

program, the output of this program could be externally converted to the 

Neutral Heat Transfer Input form and then interfaced by using TRANCITS. The 

use of the Neutral Heat Transfer Input f1le allows TRANCITS to be used w1th 

any heat transfer analyzer. 

are: 

This file 1S composed of six sets, or part1t1ons, of data. The S1X 

TITLE 

DIRECTORY 

NODAL 

ELEMENT CONNECTIVITY 

ELEMENT FACE 

TEMPERATURE 

Title Partition 

PARTITION 

PARTITION 

PARTITION 

PARTITION 

PARTITION 

PARTITION 

The title part1t10n allows the input file to have up to four 80-

character titles. These t1tles appear on the first four records of the 

input file 1n the form: 

TITLE 1 

TITLE 2 

TITLE 3 

TITLE 4 

5 

RECORD 1 

RECORD 2 

RECORD 3 

RECORD 4 
D1rectory Partit10n 

The d1rectory part1tion conta1ns count and control informat10n about the 

1nput. It 1S of the form: 

where 

NSHIFT, NUMN, NOME, NOMT, NTRECS, ECODE, ACODE 

NSHIFT = No. of records in the d1rectory part1tion - 1 

(currently NSHIFT = 0) 

NUMN = No. of nodes 1n the heat transfer model 

NOME = No. elements in the heat transfer model 

NOMT = No. of tranS1ent solut10ns on the input hIe 

NTRECS = No. of heat transfer elements w1th temperature 

(usually NTRECS = NOME for ACODE = 1 or NTRECS 

ACODE = 2) 

ECODE = Heat transfer element type code 

1-4 noded 2D element 

2-8 noded 2D element 

11-8 noded 3D element 

12-16 noded 3D element 

13-20 noded 3D element 

ACODE = Heat transfer analys1s type code 

1 - Finite d1fference 

2 - Fin1te element 

Nodal Part1t1on 

1nformat10n 

= NUMN for 

The next NUMN records compose the nodal d1rectory. It conta1ns the 

rectangular cartes1an coord1nates for every node 1n the heat transfer model. 

The form of th1S directory is: 

Node x y z 

NUMN Records 

where NODE = Node Name 

6 
x 

Y = Coordinates 

Z 

Element Connect1v1ty Part1tion 

The next NUME records in the file compose the element connect1v1ty part1-

tion. It conta1ns all the heat transfer element connect1vity 1nformat10n 1n 

the form: 

NE IE(L), L = 1, NUNPE 

NUME Records 

where 

NE = Element name 

IE = Connect1v1ty array 

NUNPE = No. of nodes each element 

Element Face Part1t1on 

This partition is omitted if the heat transfer analys1s type was fin1te 

element (ACODE = 2). If ACODE = 1 (f1nite difference), th1s part1t1on 

conta1ns NUME records of the form: 

where 

NE 

IFACP 

NE IFACP 

NUME Records 

= Element name 

= Packed element face number 

IFACP 1S a s1x-d1g1t 1nteger that relates the f1n1te d1fference face labels 

to the sides of the heat transfer element. It is requ1red when the f1n1te 

d1fference element has the ab1l1ty to produce face center temperatures as 

well as element centro1dal temperatures. Th1s 1nput allows TRANCITS to use 

all temperature informat10n produced by the heat transfer element in the 

interfacing procedure. Figure 2 shows how IFACP 1S def1ned for an e1ghtnoded sol1d heat transfer element. 

SPECIAL NOTE: When the heat transfer element be1ng used produces only 

centroidal temperatures, IFACP should be set to 123456. Therefore, for SINDA 

app11cations IFACP = 123456. 

7 
Side 4 

~ 

~ 

Side 3 

K 

I 

.... '----

,/ M .... 

Sid Side 1 

If Side 1 is Face 6 

2 5 

3 4 

4 3 

5 2 

6 1 

~----~-----Side 2 

N 

Then IFACP = 654321 

Figure 2. Element Face Definitions. 

8 
Temperature Part1tion 

This partit10n conta1ns all temperature informat10n from the heat 

transfer analysis. It can appear 1n two forms, depending on the type of heat 

transfer analysis. 

FINITE DIFFERENCE TEMPERATURE PARTITION 

(ACODE = 1) 

The form 1S: 

Solution 

1 

Solution 

2 

Solution 

NUMT 

TIMEN 

NE TMP(L) , L 

TIMEN 

NE TMP(L), L 

= 1, NTPE 

(NTRES + 1) Records 

= 1, NTPE 

eNTRECS + 1) Records 

(NTRECS + 1) Records 

FINITE ELEMENT TEMPERATURE DIRECTORY 

(ACODE = 2) 

The form is: 

Solution 

1 

Solution 

2 

TIMEN 

NODE 

TIMEN 

NODE 

TP 

(NTRECS + 1) Records 

TP 

(NTRECS + 1) Records 

9 
Solut1on I 

NUMT 

where 

T1me 

NE 

TMP 

NTPE 

NODE 

TP 

N 

eNTRECS + 1) Records 

= T1me value of tranS1ent solution 

= Element name 

= Temperature array 

= No. of temperatures each elemet 

= 5 for 2D 

= 7 for 3D 

= Node name 

= Nodal temperature 

An example of this file is given 1n Appendix B. 

4.2 SINDA TEMPERATURE INPUT FILE 

Th1s f1le 1S necessary 1f the heat transfer analysis code was SINDA. 

The users of the SINDA program have the capab1l1ty of formatt1ng the output in many ways. This formatt1ng 1S accomp11shed by includ1ng Fortran 

subrout1nes in the SINDA Job stream to control the output. An example of 

these output rout1nes for SINDA is g1ven in Append1x C. The SINDA output 

produced by these routines is defined as the SINDA Temperature Input F1le. 

The form of this file is: 

NUME 

NEI 

NEll 

TIME I 

TI 

T7 

TIME2 

TI 

T7 

etc. 

NE2 

NEI2 

T2 

T8 

T2 

T8 

NEIO 

NE20 

T6 

Tl2 

T6 

TI2 

10 

NUME 

Elements 

NUME 

Temperatures 

NUME 

Temperatures 
where 

NUME 

NXX 

TXX 

TlMEXX 

= 

= 

= 

= 

No. of heat transfer elements 

Element name 

Temperature for the element (in the same order as elements) 

Time value for this tranS1ent solution 

Th1s 1nput form 1S the required input 1nto TRANCITS. It can be obta1ned 

through SINDA by uS1ng the rout1nes given 1n Append1x C or by external 

formatt1ng of other SINDA output files. An example of the SINDA temperature 

1nput f11e is given in Appendix B. 

4.3 SINDA GEOMETRY INPUT 

This input is needed only 1f the heat transfer analysis code was SINDA. 

The SINDA output files do not conta1n geometry 1nformat1on; therefore, the 

SINDA GEOMETRY FILE must be 1nput separately 1nto TRANCITS. Th1s f11e 

contains the coordinate 1nformat1on for all the nodes 1n the heat transfer 

mode as well as the connectiv1ty of all the heat transfer elements. The form 

of the f11e 1S: 

TITLE 

NUMNOD 

NODE 

NUMELM 

X Y Z 

NUNPE 

NUMNOD Records 

NE, (IE(L), L = 1, NUNPE), IFACP 

NUMELM Records 

where: 

TITLE = 80-Character title 

NUMNOD = Number of nodes 1n the heat transfer model 

NUME = Number of elements 1n the heat transfer model 

NUNPE = Number of nodes each element 

NODE = Node name 

X,Y,Z = Rectangular cartes1an coord1nates of heat transfer 

NE = Element name 

IE = Connectivity array of heat transfer element 

IFACP = 123456 for SINDA 

nodes 

All data in this file are free format, but each node and element record must 

appear on a separate line of the f11e. An example of th1s f11e is g1ven 

in Appendix B. 

11 
4.4 MARC HEAT TRANSFER INPUT 

Heat transfer informat~on from a MARC posttape can be ~nput to the 

TRANCITS program. The MARC posttape can be ~n either formatted or sequent~al 

b~nary form; however, the following lim~tations apply: 

A. 

• Coord~nate and connectiv~ty data must be included on the posttape. 

This opt~on ~s turned on with the postmodel definit~on card. 

• Only eight-noded heat transfer sol~ds (Element Type 23) are read 

from the tape and used ~n TRANCITS. 

• Rezoning (that ~s, changes ~n coordinate and connectivity data) ~s 

not permitted. 

The contents of the posttape are assumed to be as described ~n Append~x 

MARC heat transfer ~nput ~s selected by sett~ng the var~able IHTIC ~n 

the namel~st ~nput f~le to 2 for a formatted posttape or 3 for a sequent~al 

binary posttape. 

An example of the formatted MARC postfile ~s g~ven in Append~x B. 

4.5 STRESS MODEL COORDINATE INPUT FILE 

Th~s file ~s not optional. It must be input. It contains the 

coordinate ~nformat~on for each stress node ~n the stress analysis model. 

The form of the f~le ~s: 

Node x y Z 

where: Node is the STRESS NODE LABEL 

X Rectangular 

y Cartes~an coord~nates of the stress node 

Z 

For an example see Append~x B. 

The data are in free formatted form, but each nodal record must appear 

on a separate l~ne of the f~le. 

Every stress node ~n this file will be ~nterfaced w~th the appropr~ate 

heat transfer element ~n the heat transfer model. 

12 
4.6 STRESS MODEL ELEMENT INPUT 

This f1le 1S opt1onal. It must be 1nput only 1f MARC Gauss p01nt 

temperatures or MARC centro1d temperatures are requested. It conta1ns the 

connect1v1ty informat10n for the stress model elements. The form of the f1le 

1S: 

NUMETP 

NETYP NUME 

Set 1 NE (IE(L), L = 1, NUNPE) 

No. of records = No. of elements +1 

of th1s type 

NETYP NUME 

Set 2 NE (IE (L) , L = 1, NUNPE) 

No. of records = No. of elements +1 

of th1s type 

Set 3 

Set 

NUMETP 

where 

NUMETP = No. of sets of element types 

NETYP = Element type code 

Code Element 

1 4 noded 2D element 

2 8 noded 2D element 

11 8 noded 3D element 

12 16 noded 3D element 

13 20 noded 3D element 

NE = Element name 

IE(L) = Element connect1vity 

NUNPE = No. of nodes each element 

These data are free formatted, but each element record must appear on a 

separate l1ne of the file. See Appendix B for an example of th1s f1le type. 

13 
If this f1le 1S not 1nput and 1f MARC temperatures are requested, only 

the temperatures of the stress points w111 be output. 

4.7 USER CONTROL FILE OUTPUT 

This f1le 1S not opt1onal. It must be 1nput. It conta1ns 1nput var1ables that control the 1nterfac1ng procedure. This f1le 1S 1n namel1st format 

which allows the user to 1nput only the variables related to the problem of 

1nterest. All other var1ables w111 default to internal values set by 

TRANCITS. The form of the file is: 

@ (starts in Column 2) 

$ZMAF 

Variables 

$ 

All var1ables currently ava1lable and the1r defaults are descr1bed 1n the 

following tabulat1on. 

Var1able 

Name 

ITRAN 

XOFF 

YOFF 

ZOFF 

ISL 

XMIN 

YMIN 

ZMIN 

XMAX 

YMAX 

ZMAX 

NUMZAP/ 

NZAP 

Vahd 

Values , 

0 

1 

0 

1 

Default 

Value 

0 

0 

0 

0 

0 

-105 

-105 

-105 

-lOS 

-105 

-105 

o 

Descriphon 

Do not transform heat transfer coord1nates. 

Transform heat transfer coord1nates. 

If ITRAN = 1, add the values of XOFF, YOFF, and 

ZOFF to the coord1nates of the heat transfer 

nodes. 

Do not sect10n the heat transfer model. 

Section the heat transfer model. 

If ISL = 1, these minimum and maximum 

coordinates determ1ne wh1ch heat transfer 

elements w111 be considered for the interfacing procedure. Note, if any node of a 

heat transfer element 1S outs1de the w1ndow 

defined by these m1n's and max's, the heat 

transfer element is not used. These coord1-

nates are 1n the transferred heat transfer 

coord1nate system. 

This 1nput 1S used to el1m1nate unwanted heat 

transfer elements. A shorthand 1nput 1S available (such as the 1nput NUMZAP/NZAP = 1, -7, 

-2, 12, 101, -104, 203, -218, -3) wh1ch would 

elim1nate heat transfer elements 1, 3, 5, 7, 

12, 101, 102, 103, 104, 203, 206, 209, 212, 215, 

and 218. 

14 
Van.able 

Name 

NTSTP/ 

THm 

NTMDO 

NSKIP 

TIMESP 

IHTIC 

ISPOC 

Valid 

Values 

o 

1 

2 

3 

0 

1 

2 

3 

Default 

Value 

-1 

999 

o 

1E-5 

o 

0 

Descr1pt1on 

Note, a maximum of 200 numbers may be entered 

using the NZAP array. 

This 1nput is used to select the des1red transient solut10ns from the heat transfer 1nput 

file. If 1t is om1tted or entered as NTSTP/ 

TIMD = -1, all solutions w1ll be 1nterfaced 

(such as the input NTSTP/TIMD = 10, 20, 31, 

would select three solut1ons 10 s, 20 s, 31 s). 

A maximum of 500 solut1ons may be selected. 

This input allows the first NTMDO solut1ons to 

be selected. For example, NTMDO = 5 would 

select the f1rst 5 solut10ns. 

This input allows the f1rst NSKIP solut1ons to 

be skipped or not processed. For example, NSKIP 

= 50, NTMDO = 3, would skip the f1rst 50 solutions and process the next 3. 

Th1s 1nput 

time value 

fer input. 

time values 

selected. 

is the tolerance between the des1red 

and the time value on the heat transIf TIMESP = 0.01 and TIMD = 10, any 

between 9.99 and 10.01 would be 

Heat Transfer Input Code 

Input 1S Neutral Heat Transfer Input f1le 

Input 1S SINDA temperature f1le 

Input is MARC temperature file (formatted) 

Input is MARC temperature file 

(b1nary sequent1al) 

Stress Program Output Code 

Output is Neutral Temperature f1le 

Output is formatted for NASTRAN 

Output is formatted for MARC 

(element controid) 

Output 1S formatted for MARC (element Gauss 

p01nt) 

2X2X2 

Note that any or all of these vara1ables can be 1nput using the User 

Control Input file. At least one var1able must be 1nput. It 1S recommended that 1f all defaults are des1red, the input f1le should be: 

15 
$ZMAF 

NTMDO = 999, 

$ 

An example of the file is given 1n Appendix B. 

16 
5.0 OUTPUT DESCRIPTION 

TRANCITS produces a variety of d1fferent output f11es depend1ng on wh1ch 

stress analysis output code has been selected. Two files are always produced, 

the hard-copy f1le and the neutral temperature output f1le. The two opt1onal 

output files are the NASTRAN temperature f11e and the MARC temperature file. 

The follow1ng sect10ns descr1be all of these var10US outputs. 

5.1 NEUTRAL TEMPERATURE OUTPUT FILE 

The Neutral Temperature Output file 1S the pr1mary output of TRANCITS. 

It contains the 1nterfaced temperature for all of the stress p01nts processed 

for every tranS1ent solution selected. The form of the f11e is: 

TIME = 

NODE 

TIME = 

NODE 

where 

~ 

NODE 

TEMP 

~ 

TEMP 

F1rst 

Solution 

~ 

TEMP 

Second 

Solution 

Last 

Solution 

= The TIME VALUE of the tranS1ent 

= Stress point name 

= Temperature at th1s p01nt 

Formats 

2X, TIME, IX, = IX, F12.6 

18, F12.6 

This output 1S always produced. It can be reformatted w1th s1mple 

conversion programs to produce temperature 1nput for almost any stress 

analys1s code. 

An example of this output 1S 1n Appendix B. 

17 
5.2 NASTRAN TEMPERATURE FILE 

Th1s opt10nal output f1le is produced if the stress analysis program 

used is NASTRAN. The stress point temperatures are formatted so that they 

can be d1rectly 1nput into NASTRAN. The form of the f1le 1S: 

$ TIME = XXX 

TEMP NC NODE Temp 

First Solution 

Last Solution 

where 

XXX = T1me value of the solut10n 

Formats 

$, IX, TIME, IX, = IX, F12.6 

TEMP, 4X, 2IB, FB.2 

NC = Case number (1 through number of solut10ns) 

NODE = Stress point name 

TEMP = Temperature of the stress p01nt 

An example of this file 1S 1n Append1x B. 

5.3 MARC TEMPERATURE FILE 

Th1s output f1le 1S produced 1f the stress analys1s code to be used 1S 

MARC. The temperature 1nformat10n is generated as a set of Change State 

cards with a header card to 1nd1cate the t1me step from the heat transfer 

analys1s. This data must be manually edited into the MARC 1nput. 

S1nce MARC 1S a somewhat un1que code 1n that 1t requ1res either centr01d 

or 1ntegrat10n p01nt temperatures, an element connectivity f1le must also be 

1ncluded as 1nput to TRANCITS. Currently, temperature 1nput 1S ava1lable 

only for e1ght-noded stress solid elements (Element Type 7). 

5.4 HARD-COPY OUTPUT DESCRIPTION 

The hard-copy output of TRANCITS has bas1c informat10n about the heat 

transfer model and the stress model to be interfaced. It starts with a 

TRANCITS header page, then it prints counts of the number of nodes and 

elements 1n the heat transfer model. It also pr1nts the total number of 

solution time steps and the number of solut10ns processed for the run. 

18 
Next, 1t pr1nts a header 11ne saying "The following stress p01nts were 

not contained by any heat transfer element" and that "SNOD XO YO 20 NEARE XN 

YN 2N DIST IWARN." 

If no numbers appear under th1s header, all stress points 1n the stress 

analys1s model werQ 1ns1de a heat transfer element. Any stress p01nt outs1de 

of the heat treat elements will appear under the header. The definit10ns of 

each of the columns are: 

SNOD = Stress node name outs1de all heat transfer elements 

XO 

YO = Or1ginal coord1nates of this stress point 

20 

NEARE = Heat transfer element name nearest to this stress p01nt 

XN 

YN 

2N 

Surface coord1nates on nearest heat transfer element for wh1ch 

the temperature of this stress p01nt was computed. 

DIST = Distance between (XO YO 20) and (XN YN 2N) 

IWARN = 0 - Surfacing procedure converged normally 

= 1 - Surfac1ng procedure had difficulty 1n converg1ng. 

(Check temperatures at this stress p01nt.) 

NOTE: The two ma1n points to look for 1n th1s pr1ntout are the magn1tude 

of the d1stance parameter (1f 1t 1S large, then the stress p01nt 

was far outs1de the heat transfer model) and the IWARN parameter. 

The last portion of the hard-copy output deals only with f1n1ted1fference heat transfer codes. The header reads: 

"Corner Node Compar1son For T1me = XXX" 

"TCRITICAL For Th1s Case Is Set to 99°" 

"NODE /I MAXOT TBAR ;'..'';COMMON HIT ELMS -l;l," 

Th1s printout 1S a consequence of the heat transfer element corner node 

calculation. S1nce the fin~te d~fference codes yield temperatures at the 

element controids and the element face centers, TRANCITS uses these temperatures to compute the element corner temperatures for each element. These 

corner temperatures are compared and if the d~fference of the temperature at 

the same corner 1S greater than 99°, then that corner 1S flagged 1n th1s 

pr~ntout. The column def1nit1ons are: 

19 
Node No. 

Max PT 

TBAR 

Co~on Heat Transfer 

Elements 

= 

= 

= 

= 

Heat transfer corner node name 

Maximum delta temperature at that corner 

Average temperature at that corner 

Heat transfer elements bound1ng the corner 

If this type of pr1ntout occurs it 1S an indicat10n of large thermal 

gradients in the heat transfer model. If the MAXDT value is large there 15 a 

potent1al problem with the adequacy of the mesh density of the heat transfer 

mesh. When this s1tuation occurs it 1S recommended that the heat transfer 

analyst be consulted about the valid1ty of the results in the h1gh gradient 

area. 

The last message printed 1S 

"THE INTERFACING PROCEDURE HAS BEEN COMPLETED." 

20 
6.0 RUN COMMANDS 

The current verS10n of TRANCITS (Version 1.0) is operational on the NASALewis IBM 370 system. 

Upon release of this user's TRANCITS manual, a formal procedure w111 not 

have been set up to run in product10n mode at NASA-Lewis. Th1s procedure 

will be established in the near future and will be communicated to all users. 

21 
This Page Intentionally left Blank 
APPENDIX A 

CONTENTS OF MARC POSTTAPE - VERSION K.l 

Th1s append1x contains the documentat10n of the MARC Heat Transfer 

Input. Refer to the MARC manuals for more information. 

23 
APPENDIX A 

CONTENTS OF MARC POSTTAPE - VERSION K.l 

The follow1ng describes the contents of the 10 blocks found on the 

posttape. This 1nformat1on is the same for e1ther the b1nary or the 

formatted posttape. 

BLOCK NUMBER 

BLOCK I 

NR;80Al;(TITLE(J),J=I,70) 

BLOCK II 

NR; 613; INUM 

LNUM 

MNUM 

NDEG 

NSTRES 

INOD 

IPSTCC 

IPSTYP 

NCRD 

NNODMX 

DESCRIPTION 

TITLE 

TITLE(J) 1S Jth character 

CONTROL INFORMATION 

Number of variables each element 

Number of nodes 1n mesh 

Number of elements 1n mesh 

Number of degrees of freedom each node 

Maximum number of integration p01nts 

each element. 

NSTRES = 1 1f ALL POINTS 1S not flagged. 

NOTE: For constant stra1n elements 

hav1ng only one 1ntegrat10n p01nt 

NSTRES = 2 1f ALL POINTS is flagged. 

Number of nodal var1ables - see note 

below. 

Connect1v1ty and coord1nate flag 

(O-none; I-given) 

Type of tape (0 - b1nary, 1 - formatted, 

2 - both) 

Number of coordinates each node 

Maximum number of nodes each element 

24 
IANTYP 

ICOMPL 

NBCTRA 

IPOSTR 

IDM4 

IDMS 

IDM6 

IDM7 

Analys1s Type Flag: 

1 = Displacement 

forces 

without react10n 

2 = Displacement - with react10n forces 

3 = Dynam1cs - w1thout react10n forces 

4 = Dynam1cs - w1th react10n forces 

5 = Heat transfer 

6 = Rig1d plast1c flow - w1thout react10n 

forces 

7 = Rig1d plast1c flow - with react10n 

forces 

8 = Harmon1C analys1s - w1thout react10n 

forces 

9 = Harmon1c analys1s - wLth react10n 

forces 

10 = Harmon1c analys1s with dynam1cs -

w1thout reaction forces 

11 = Harmon1c analys1s w1th dynam1cs -

with reaction forces 

12 = Heat transfer w1th fluxes 

13 = Joule heat 

14 = Bear1ng 

15 = E1genvector (modal) 

16 = E1genvector (buckle) 

Set to 0 if real analysis 

Set to 1 if complex analysis 

Number of nodes with transformations 

Posttape revision number, 1 for th1S 

release (J3) 

Not used; reserved for future 

expans10n 

25 
Note: 

INOD = {NDEG}*JNODE 

IF IANTYP = 1· , JNODE = 1 

2 on 2 

3 3 

4 4 

5 1 

6 1 

7 2 

8 1 

9 2 

10 3 

11 4 

12 2 

13 3 

14 5 

15 1 

16 1 

IANTYP = 15 only appears dur1ng sub1ncrements 

IANTYP = 16 only appears during sub1ncrements 

If IANTYP = 8 and ICOMPL = 1 JNODE = 2 

If IANTYP = 9 and ICOMPL = 1 JNODE = 4 

BLOCK III CODE NUMBERS ASSOCIATED WITH ELEMENT 

VARIABLES 

NR;6I3;JPLOT(J) 

BLOCK IV 

NR;6I3;LM(l) 

BLOCK V 

LM(2) 

LM(3) 

LM(NNODMX+2) 

INUM records; J=I,INUM 

Element var1able code (see MARC Manual 

Volume C, p C3.10-7) 

CONNECTIVITY LIST 

If IPSTCC 1S zero th1s block w1ll be 

omitted 

MNUM records; J=l,MNUM 

Element type 

Number of nodes in th1s element 

First node of Element J 

Last node of Element J 

COORDINATE LIST 

If IPSTCC 1S zero this block w1ll be 

omitted. 

26 
· NR;6E13.6; SUM(l) 

BLOCK VI 

SUM(2) 

SUM (NCRD) 

NR;6I13; LM(I) 

BLOCK VII 

NR;6E13.6; DC1,1) 

D(2,1) 

D(3,1) 

D(1,2) 

D(2,2) 

D(3,2) 

D(1,3) 

D(Z,3) 

D(3,3) 

LNUM Records; J=l,LNUM 

F1rst coord1nate of Jth node 

Last coordinate of Jth node 

TRANSFORMATION LIST 

If NBCTRA 1S zero or 1f IPSTCC 1S zero, 

th1s block will be om1tted. 

Binary 1 record 

Formatted (NBCTRA-1)/6 + 1 records 

List of nodes wh1ch have transformat1ons 

apphed. 

TRANSFORMATION - DIRECTION COSINES 

If NBCTRA 1S zero th1s block w111 be 

om1tted. 

One record for each node 11sted 1n 

Block VI 1f binary tape. 

Two records each node if formatted tape. 

Transformat10ns are for local to global 

system. 

Blocks VIII, IX, and X are repeated for each increment. 

BLOCK VIII 

NR;6E13.6 Xl(l) 

X1(2) 

X1(3) 

INCREMENT, TIME, AND FREQUENCY 

Trans1ent T1me 

Increment number 1S a real formed 

as I + J/100 

I 1S the static 1ncrement number 

J is either the harmonic subincrement 

or the E1genvector number. 

Frequency 

27 
BLOCK IX 

X1(4) 

Xl(S) 

X1(6) 

DO 411 II=l,MNUM 

DO 413 JJ=l,NSTRES 

NR;6E13.6; SUM(l) 

SUM (INUM) 

413 CONTINUE 

411 CONTINUE 

BLOCK X 

NR;6E13.6; SUM(l) 

SUM (NDEG) 

SUM(NDEG+1) 

SUM (2-,'·NDEG) 

SUMCINOD) 

Flag to read new Blocks II, III, IV, V, 

VI, VII. Set to 1 to read these blocks 

again. Th1S is only used 1n a rezon1ng 

analys1s. 

IANTYP analysis type flag 

Not used; reserved for future expans10n. 

VALUES OF ELEMENT VARIABLES 

If lNUM 1S zero, th1S block is omitted. 

If IANTYP=lS or IANTYP=16 , th1S block 

is omitted. 

MNUM*NSTRES Records 

Value of 1st element var1able for 

Element II, Integrat10n P01nt JJ. 

th Value of lNUM element var1able for 

Element II, Integrat10n Point JJ. 

VALUE OF NODAL VARIABLES 

If INOD 1S zero, th1S block is omitted. 

LNUM records 

Nodal displacements, veloc1t1es, 

accelerat1ons, and react10ns 

During stat1c or dynam1c increments: 

The f1rst NDEG quant1t1es are: 

IANTYP = 1, 2, 3, 4, 8, 9, 10, 11 

Total nodal displacements 

28 
IANTYP = 6, 7 

Total nodal velocity 

IANTYP = 5, 12, 13, 

Total nodal temperature 

IANTYP = 14 

Pressure 

The second NDEG quantit~es are: 

IANTYP = 2, 7, 9 

Nodal reactions at boundary 

condit~ons; res~dual load 

correction elsewhere 

IANTYP = 3, 4, 10, 11 

Nodal veloc~ty 

IANTYP = 12, 13 

Fluxes 

IANTYP = 14 

FORCE X 

The third NDEG quantit1es are: 

IANTYP = 3, 4, 10, 11 

Nodal accelerat~ons 

IANTYP = 13 

Nodal voltages 

IANTYP = 14 

FORCEY 

The fourth NDEG quantit~es are: 

IANTYP = 4, 11 

Nodal react~ons at boundary 

cond1tions; res1dual load 

correct10n elsewhere 

IANTYP = 14 

FORCEZ 

The fifth NDEG quantit1es are: 

IANTYP = 14 

MASS FLUX 

29 
Dur1ng subincrements: 

The first NDEG quant1t1es are: 

IANTYP = 15 

IANTYP = 

IANTYP = 

Nodal components of dynam1c 

mode 

16 

Nodal components of buckle 

Eigenvector 

8, 9, 10, 11 

Real components of harmon1c 

displacement 

The second NDEG quant1t1es are: 

IANTYP = 9, 11, and ICOMPL = ° 

Real components of harmon1c 

react10n force 

IANTYP = 8, 9, 10, 11, and ICOMPL = 1 

Imaginary components of 

harmon1c displacement 

The th1rd NDEG quantities are: 

IANTYP = 9, 11, and ICOMPL = 1 

Real components of harmon1c 

react10n force 

The fourth NDEG quant~t~es are: 

IANTYP = 9, 11, and ICOMPL = 1 

Imag1nary components of 

reaction force 

NOTE: There w1ll always be INOD quant1ties; the record may be padded 

w1th zeroes 1n sub1ncrements. 

IF IANTYP = 14, same as INCREMENTS. 

30 
NOTES: 

NR 

NR;Format 

Indicates the beginn~ng of a new record 

for the b~nary posttape. 

Ind~cates the beg~nn~ng of a new set of 

~nformat~on wh~ch ~s to be read w~th the 

follow~ng format. 

31 
This Page Intentionally left Blank 
APPENDIX B 

Th1s Appendix shows examples of the follow1ng file formats: 

• Neutral Heat Transfer Input 

• SINDA Temperature Input 

• SINDA Geometry Input 

• MARC Postinput 

• Stress Model Coordinate Input 

• Stress Model Element Input 

• User Control Input 

• Neutral Temperature Output 

• NASTRAN Temperature Output 

33 
NEUTRAL HEAT TRANSFER INPUT FILE 

This is the Test Case data. This is the new file form to be 

delivered for NASA's version -- we are using a very simple test TITLE 

case, as follows: 

09-15-83 (SR MARTIN) 

0 5.1 20 1 Zg 11 1 < DIRECTORY 

1 1.0~000ii3 1.000000 1.0!3!O~00 

2 1.0000100 2.00f}f}00 1.elO!24!3@12I 

3 1.13000013 3. 0QH;~Q'g 1.0:30000 

11 2.0000130 1.13e00Q'0 1.00101000 

• NODAL 

24~: 5.013000ii3 3.0ii3130130 3.00001313 

2'51 6.06130131.3 1.00130013 3.131.360(3121 

2'5z' 6.013010130 2.00013G~ 3.0013000 

253 6.0G00!240 3. 012113 13 Gk3 3.13f}01.300 

1 1 2 lZ 11 1131 1132 112 111 

2 2 3 1 ... · .J 12 1132 103 113 112 

11 11 12 ,.., ... , /.../... 21 111 112 122 121 ELEMENT 

12 12 13 23 22 112 113 1 ~~, 

.~.j 122 CONDUCTIVITY 

21 21 22 .... ,.., .J/... 31 121 122 13Z 131 ....... 

/.../... 

,.., ... , /.../... 23 33 32 12Z 123 133 132 

31 31 .... .-, .J/... 4~ 41 131 132 142 141 

141 141 1" ~/... ... , • 15: 151 241 242 252 251 

142 142 143 153 152 242 243 253 25Z, 

1 123456 ELEMENT 

z: 123456 FACE 

1 1 123456 

12 123456 

21 123456 

ZZ 123456 

132 • 12345e. 

141 12345': 

14: 123456 

0. 

1 201!1.e 1013.0 200.0 :?00.0 200.0 21313.13 2130.0 

2 200.0 10{(!.0 2100.0 3010.13 200.0 200.0 200.0 

11 400.13 300.0 400.0 500.0 400.10 4013.0 4013.0 TEMPERATURE 12 400.0 31313.0 4100.0 5010.0 400.0 4013.13 41313.13 

21 61013.13 5010.0 6010.0 7130.0 6130.0 600.0 6010.10 

22 600.13 5130.13 6100.10 7010.13 61313.10 6~0.0 600.1!1 

31 800.0 7010.0 8130.0 9130.0 800.0 80~.g 8013.0 

• 34 
SINDA TEMPERATURE FILE 

zss 

leI 1~~ .t. lv~ 1f4 liS U~ 197 He Ib~ 119 

111 112 113 lit 115 1lt. 117 119 119 lzg 

121 122 r''> 1...1 1~~ 125 12~ 127 12B 129 1~~ 

131 131 133 134 1~ 13b 137 138 1~ .... 

1~ 

141 142 1~3 H~ HS Hb 147 148 149 ISS 

151 1-' J ..153 1~~ ISS 15b 157 158 159 Ib3 

1M Ib2 163 1M u.s 1M. It.7 108 169 173 LIST OF 

171 172 173 17~ 175 In In 178 179 1~-.. ~ ELEMENTS 

181 1&2 183 1~ le5 2;1 zsz ,53 ZJ~ m ~b Z;i z~a Z19 2U 211 212 213 214 215 

2lb 217 21B 219 m 221 212 223 224 2ZS 

22b '1:17 US 229 231 231 Z3Z Z33 £2-4 Z35 

2U 237 238 239 2Y 2H 242 243 2" 245 

24b 247 2~ 249 29 251 ZSZ 253 ZS4 ZS5 

ZSb 2S7 ZSB 259 2bi Zbl Zb, 203 2M 2b5 

• ?~1 "I. 372 373 374 375 3n 3n 378 379 ~ 

381 m 383 3S~ 3e5 

'.3916F- ii' 

1.~59E 93 5.49561£ i3 '.sa=vG2£ 33 '.49&11E J3 '.5fiSbE 93 9.49754£ §3 

'.51943E 93 '.49niE 83 '.S11'Tt 13 j.~if1l1E i3 I.SJ9S~'£ 93 i.4%9iE 13 

'.5IZS7E 93 '.48827: .3 ;.s!191E 93 i.49241£ J3 i.4S9SZE .3 §.493bbE 83 

'.4874ZE 93 J.4ffiiE ,3 J.tSeS8E '3 '.~'SE 63 '.«mar: i3 g.~n .3 

'.489oLE 93 '.4a934E J3 '.49.3Z: 13 '.~weE ~3 j.4~'~ 83 J.4S97SE )3 

'.48925E 93 J.4S939E 13 '.48SSbE J3 a.~ .3 i.48949E i3 9.lSSS1E 13 LIST OF 

1.49311£ 93 J.491b3E '3 '.49~ '3 i.49:1SE .~ '.49SIfE '3 3.49240£ '3 TEMPERATURES 

'.49543E 83 '.4~bE t3 '.49b31E .3 •• 49391£ 13 '.47638£ '3 '.47~~ '3 

'.47225t '3 '.49114£ 13 •• ~;sc: .3 J.4..~7LE .3 '.49biSf: .3 1.49548£ J3 

1.49394£ .3 '.49SZSE .3 '.49749£ .3 1.49597£ 53 ~.4998bE .3 1.49SbSE .3 

'.49691£ i3 '.Sl1SbE .3 1.&093SE .3 '.49bS3E J3 1.51272£ .3 1.49954£ .3 

'.49599£ .3 1.~:sbE .3 '.49'i'91£ 13 '.495CfJt 13 i.5J.4..4~ 13 '.49964£ 13 

35 
SINDA GEOMETRY FILE INPUT 

SINOA GEOMETRY FILE - 30 COMc~STOR LI N;::R MODEL 

500 255 8 

1 o 28640 15 75804 O. 

2 0 20570 15 792.:m- O. 

3 0 53550 15 74303 O. 

4 0 53780 15.77885 O. 

!5 0 29100 15 82606 O. 

6 o ;,4010 I;, 61478 O. 

7 0.78450 15 72802 O. 

8 0.78690 15 76570 O. 

9 0.78930 15 80339 O. 

10 1 03350 I;, 71301 O. 

11 1.03600 15.75256 O. 

12 1.03840 15.79210 O. 

13 1.10760 15 70856 0 

14 1 11090 15 16~12 O. 

15 1.11410 1!5 81757 0 

16 1.16370 15 70524 0 

17 I 16690 15.75970 O. 

, t t t 

717 2 02::>20 14 86239 5 40947 

718 2.11210 14 93652 5 43645 

719 2.11230 14 91201 5 42753 

720 2 11250 14 86749 5.41860 

721 2.11270 14 853u7 5.40972 

722 2.19960 14 93681 5.43656 

723 2 19980 14 91239 5.42767 

724 2 20000 14 8SS07 5 41882 

725 2 20020 14 863Q5 5 40993 

101 203 204 202 201 3 4 2 1 123456 

102 204 206 205 202 4 6 5 2 123456 

103 207 208 204 203 7 8 4 3 123456 

104 208 209 206 204 8 9 (; 4 1234::>6 

105 210 211 208 207 10 11 8 7 123456 

106 211 212 209 208 11 12 9 8 12~456 

107 213 214 211 210 13 14 1 1 10 1234~6 

108 214 215 212 211 14 15 12 11 123456 

109 216 217 214 213 16 17 14 13 12:3456 

110 217 218 215 214 17 18 15 14 1234~6 

111 221 222 217 216 21 22 17 16 123456 

112 222 223 210 217 22 23 18 17 123456 

113 228 229 222 221 28 29 22 21 123456 

114 229 230 223 222 29 30 23 22 123456 

115 234 235 229 226 34 35 29 28 123456 

116 235 236 230 229 35 :36 30 29 123456 

117 224 220 218 223 24 20 18 23 123456 

118 220 219 215 218 20 19 15 18 123456 

119 227 226 220 224 27 26 20 24 123456 

120 226 225 219 2?0 26 25 19 20 123456 

121 233 232 226 227 33 32 26 27 123456 

122 232 231 225 226 32 31 -. ...... 26 123456 

123 239 238 232 233 39 38 32 33 123456 

124 238 237 231 232 38 37 31 32 123456 

125 238 239 242 241 38 39 42 41 123455 

126 237 238 241 240 37 39 41 40 12~456 

127 241 242 245 244 41 42 45 44 1~:3456 

128 240 241 244 243 40 41 44 43 123456 

• • , 

i 

• 

I 

Y • t • , t 

36 
MARC POSTFILE/TAPE INPUT 

CO~2>USTOR MARC POST FILE SIMULATIO', 

0;:;.OO2~~ 1 8 1 

1 I 3 8 5 0 

0 I 

43 8 3 4 2 -----f2aI29T~o 

43 8 4 6 :5 2 

12913113::127 

43 8 7 8 4 3 

13213312912~ 

43 8 8 9 6 4 

133134131129 

43 8 10 11 8 7 

1351361:;3132 

43 8 1 I 12 9 8 

136137134133 

43 a 13 14 11 10 

49549~490~91 

43 8369368354365 

4944934894')0 

43 837537437 0371 

50049S495.~S6 

43 8374373369370 

499496-"19-"495 

43 8373372368369 

498497493'::94 

o 286400E 00 0 157580E 02 0 

o 288700E 00 0 157920E 02 0 o 535500E 00 0 157430E 02 0 

o 537800E 00 0 157789E 02 O. 

O.291000E 00 0 158261E 02 O. o 540100E 00 0 158148E 02 O. o 7S4500E 00 0 157280E 02 0 

+ 

0 211230E 01 0 149J20E OZ 0 

0 211250E 01 0 146875E 02 0 

0 211270E OJ 0 148631E 02 0 

0 219960E 01 0 149368E 02 0 

0 219960:' 01 0 149124c. 02 0 

0 220000E 01 0 148881E 02 a 

0 220020E 01 0 148637E 02 0 

a 393700E 04 0 10000CE 01 a 

O. ::11 0434c. 03 o 499.a32E 03 

a 512448E 03 

0 50'J855E 03 

0 490034E oj 

0 491059E 03 

0 50977SE 03 

0 49651SE 03 

o 485630E 03 

0.!510112E 03 

0.494647E 03 

0.482880E 03 

37 

S427S3E 

541860E 

540972E 

543656E 

542761E 

541882E 

5a0993E 

01 

01 

01 

01 

01 

01 

01 

0 0 500000E 0' 
STRESS MODEL COORDINATE INPUT FILE 

1 3<:980 15 86550 O. 

:2 1. 349;:>'""0 1 "!:JQ:ro:m----O . 

3 1.34950 15 80435 1.03643 

4 1 34980 15 83179 1.03823 

5 1.39740 1'5 80455 1 03644 

6 1.39/.;)j 15 c",169 1 03822 

7 1.39740 15 83a.50 O. 

6 1.39753 1::: 66:70 O. 

9 1 .34920 15 81080 O. 

10 1. 34920 15 lib::') 1. 03463 

" 1 39727 15.77731 1.03465 

12 1.39727 15 81120 O. 

13 1.34920 15 67~54 2 06372 

14 1 • 34£150 15 10200 206"1"3"1 

15 1 34980 15 73C07 2 07090 

16 1 39727 15 67593 2 06378 

17 I 39740 15 70300 2 C6734 

16 1. 397;:>3 15 ns97 2 07089 

19 1 44533 15 67633 2 06383 

20 1 44533 15 77771 1 03465 

21 1 44533 15 81160 0 

22 1 4':530 b 70jl0 2 06735 

23 1.44527 15.72997 2 07089 

24 1.44530 '5 80465 1. 03645 

2'5 1.44527 15 83169 1.03822 

26 1. 44530 15 83e60 O. 

27 1.44527 15 8S570 0 

28 1.34920 15 50711 3.0i3399 

29 1 34950 15 53408 3 08935 

30 1 34980 15 ;:>6105 j 09472 

31 1.39727 15 50750 3.08407 

32 1 39740 15 53428 3.08940 

33 1 39753 15 56095 3 09470 

34 1 . 44533 15 50/89 3 08415 

35 1.44530 15 53437 3 08942 

:36 1.44527 15.56C95 :3 09470 

37 1 .49340 15 50828 3 08423 

38 1 49520 15 ;:>3457 ~ 08945 

39 1.49300 15 '56C86 :3 09468 

40 1 49340 1'5 67673 2 0638t 

41 1 49340 15 77811 1 03471 

42 1 . 49340 15 81200 O. 

43 1.49320 1'5 70330 2.06738 

44 1 49300 15 72967 2 07088 

45 1 49320 15 80455 1 03646 

45 1.49320 I::: a3e60 o. 

47 1.49300 15 83159 1.03821 

48 1.49300 15 86560 O. 

49 1 34920 15 27206 4 09214 

50 1. 34950 15 2961:>2 4 09925 

'51 1.34980 15 32519 4.10637 

'52 1.39727 15 2724'5 4.09224 

'53 1.39740 15 29682 4 09931 

38 
STRESS MODEL ELEMENT INPUT FILE 

1 

11 452 

1 790 789 iSe 788 787 785 783 784 

2 789 782 781 78G 785 780 779 783 

3 762 163 lez 761 IBO 761 760 /79 

4 763 726 725 762 761 724 723 760 

5 760 78G 777 778 764 783 775 776 

6 786 781 770 777 783 779 769 775 

7 181 762 759 770 779 760 75e ISS 

8 762 725 722 759 760 723 721 7S8 

9 778 777 753 754 776 775 7!1 7~2 

10 777 770 750 753 775 769 749 751 

11 770 109 736 750 769 /oB 735 749 

12 759 722 717 736 758 721 716 735 

13 78-' 785 783 784 774 773 771 772 

14 785 780 779 783 773 768 767 771 

15 7",0 751 700 119 7to 757 756 167 

16 761 724 723 760 7'57 720 719 756 

17 784 783 775 776 772 77' 765 766 

18 783 779 7G9 775 771 767 764 765 

19 779 760 108 769 167 756 750 764 

20 760 723 721 758 756 719 718 755 

21 776 775 751 752 766 765 747 748 

22 775 769 749 751 765 764 746 747 -----

427 497 49b 492 • 493 '157 455 ~2 453 

428 496 495 491 492 4'56 455 451 452 

429 658 638 637 657 630 629 627 628 

430 638 608 607 637 629 600 599 627 

431 soli 071 570 607 600 56e 50:::l 599 

432 657 637 636 656 628 627 625 626 

433 637 607 60S 636 627 599 596 &5,5 

434 607 570 569 606 ~99 559 558 598 

43$ 656 636 635 650 626 62~ 62~ 624 

436 636 606 605 635 625 598 597 623 

437 606 569 565 605 598 558 557 597 

438 6"'" ~~ 635 631 6'51 624 623 621 622 

439 605 Sao 601 631 623 591 !S9o 621 

440 605 '565 561 601 597 557 556 596 

441 6:30 629 627 628 595 59~ 552 593 

442 629 600 '599 627 594 '591 590 592 4~3 600 560 559 5M 5;1 55:; 504 590 

444 628 627 625 626 593 592 588 589 

445 627 599 598 625 592 590 587 588 

446 599 559 558 598 590 554 553 587 

447 626 625 623 624 5SS 586 580 586 

448 625 598 597 623 588 587 584 585 

449 598 558 557 ~97 587 553 552 584 

450 6N 623 621 622 586 585 :582 :583 

aSl 623 !g, 596 621 SS5 504 5al 582 

452 597 557 556 596 584 552 ~51 S81 

39 
USER CONTROL FILE 

SZMAF 

I HTle= j, 

ISPOC=l, 

NTMDO=3, 

NUMZAP/N'A~=l,-11,-2.51. 

NTSTP/iTRO=10,30,5I, 

$ 

40 
NEUTRAL TEMPERATURE OUTPUT FILE 

TIME = 3916 (l~OOOO 

-----1 ---q;-~.::;, I I 75 

2 475 362076 

3 473 061939 

4 472 964638 

~ 477 105648 

6 477 014626 

7 479 731720 SOLUTION 1 

e 476 41640£1 

9 476 380043 

10 476 101242 

11 4aO 371784 

12 481 297924 

13 477 669193 

14 475 659309 

15 473 984970 

16 482 275578 

17 479 570258 

18 0477.891411 

19 0487 512~89 

20 485 714573 

TIME = 3947 000000 

1 913 990150 

2 914 ?88368 

3 912 673592 

--.:r 912 7i::0978 

5 921 552902 

6 921 516510 

7 925 217186 

a 9i::3 021370 

9 914 720711 SOLUTION 2 

10 912 963562 

11 923 802025 

12 9':0 0435.5 

13 911 039658 

14 910 560577 

ns 910 516815 

10 ::lZZ 9Ia~t:>j 

17 919 726219 

18 918 977898 

19 937 802277 

20 937.246u3 

TIME = 4055 000000 

1 1027 800430 

2 1028 45C317 

3 1028 727051 

4 1029 7847104 

5 1044.299011 

6 104.4 251038 SOLUTION 3 7 1045 390854 

8 loa..! 329132 

9 1027 999115 

10 1029 331284 

11 1046 3469'5.4 

12 104S ,Q~!:Ib41 

13 1030.756936 

, ... 1030.794373 

15 1030 3196e7 

lb 104'3 565552 

17 1046.791626 

18 1045 744339 

19 1077 498398 

20 1074.42./83 

41 
NASTRAN TEMPERATURE FILE 

" TIME = 3916 OCOOOO 

TEMP I 472 25 

TEMP 2 475 36 

TEMP 3 473 06 

TEMP 4 472 96 

Te-A? 5 477 11 

TEM? 6 477.01 

TEM? 7 479 73 TE~P 6 476 42 

hMf' 9 476 30 

TEM? 10 ~76 10 

TEMP 1 1 480 37 

TEMP 12 481 30 

Te-M? 13 417 51 

TEMP 14 475 66 

TEMP 15 473 98 

TEMP 16 ~82 26 

IEMP 17 479 58 

TEMP 18 477 89 

TEMP 19 487.51 

TEMP 20 485 71 

TEMP 21 487 03 

TEMP 22 464 75 

TEMP 23 483 24 

TEMP 24 482 65 

Te-MP 25 J8Z 58 

TEMP 26 485 24 

TEMP 27 482 10 

TEMP 28 ~79 20 

TEMP 29 476. 49 

TEMP 30 476 07 

TEMP 31 483 04 

TEMP 32 480 08 

hHP jj 41'3 67 

TEMP '34 487 76 

TEMP 35 ~84.96 

TE"1P 36 484 60 

42 
APPENDIX C 

SINDA FORTRAN OUTPUT SUBROUTINES 

Th1s append1x contains listings of the SINDA Fortran output subrout1nes 

used to produce the SINDA Temperature Input File. Refer to the SINDA Manual 

for proper use of these routines. 

43 
SINDA FORTRAN OUTPUT SUBROUTINE 

~V.WW ••• ~VVWV~~r_.~~~_~~~~~~ .w~~~ __ •• ~~M.M.a," •• 

aCD !!~=,;~ 

~WW~~W •• W· •••• ·.~.~ ••• K •• W.~W •• M ••• '.~.W.MWW.~.MM 

ASSIGN DYNAMIC STO~~GE PE~U!~(~EH1S 

SIZE OF X DEPEHt·S ON Hun~tP OF HeOES. SOLUTIOH ROUTIHE. 

AND sunRCUTINES CALl~D 

CO~MOH'LOGIC'LHOOE.LCC~D.LCONST.LAR~AY 

DII':ENSION X(200) NDI~:200 

NTH =0 

THE FOLLCWING conE FIllS T~[ DYNA~lC STO.AGE ARPAY 'X' ('HX') WITH THE NODE NU~BERS IN THE SA~E CROEP AS THE 'T' t~rAY. 

TriEN nEGIN~I~G WITH THE FIRST NCDE HU~BER. THE H~DES ARE I 

IJl\ITT EN TO TilE FILE USED ilY 'TEr:rSCf.H'. _ f .. -,.l '" mv -rt't'IIO • 

$$$ REFER TO SINDA USER'S M~~CAL. P;CE 5-' $S$ ~. ~. 

re.b'S ;r,.,>' ~e neU'! '". 

HOTE: THIS HA~PENS ONLY ON T~E FIRST CALL TO EXECUTIO~ DUE 

TO THE CHECK ON T~E STATUS OF 'LHODE'. 

IF(LHODE.GE.O) GO TO 20C ~ Loe : HTH ... ~ 

LEN = HOlM • "..,. .r. I • 

,!;~tl t!~DEAp(1) 5,.-C. --.=-XeCII7,c.,., . 

leO FOR~:'T(15) 

tis -cl1lr~"CA lied. on'~ .. 

110 

:!OO COIIlItWE 

EXECUTE PP("GPAM 

CNrpl~D 

CSC;T:rW 

END 

¥.MW·~.¥~¥~¥¥¥~WW.W.¥WWW.WWW~WKW ••••• W"W."W".W.WW 

aco 3Vt.PIt.BL~S 2 

W¥ •• V~.-4¥W~~~_W~~·MW¥·WWW¥W.¥W.W" •• "~M •• W~ •• ".¥W 

EIID 

~¥ •• ¥.~._ ••• MWW¥M ••• M~ •• ¥ ••• M ••••••• " ••••••• WV¥WW 

BCD 3pnTPUT CAtt c 

¥~~w ••• ~~w~.~ ••••••••• 4.~* •• w~wW~W*.*MWN~WWN~.M~_ 

lPPHTF 

~oITE NCDE TE~PERATURES TO TAPE 

KTEST:t<TEST+l 

~OITE(6.l30)rTEST.Tr~E~ 

l!O fOFM~T(~X."~5$ NOTE· DATA SET =",15." WR1iTEH TO TE~?SCAH". 

• " FILE AT TI~E : ".FlO 3) IoIPlTE (23.140) lIMeN w~-:-rr: ·-r'".E -r; O~~f":#r. +'olE' 

l:;:ITE(23.140 >CTC I). I=l.HltT)," I 

140 fOR~I:,1[6El2.5) ....... 1N.r·.+~ r-d.+,'"v~ +#"Pf'!",tt.(Y'~S 

EIiD 11/';1(1. I 

44 

+0 

c 

C 

C 

C 

C 

C 

C 

C 

F 

F 

F 

F 

C 

C 

C 

C 

C 

C 

C 

C 

C 

F 

F 

F 

F 

C 

F 

C 

F 

C 

F 

F 

C 

C 

F 

c 

C 

r 

c 

c 

C 

c 

C 

0.: 

C 

C 

F 

F 

F 

F 

F 

F 

F 

C 

C 

C 
End of Document