THE UNIVERSITY OF MICHIGAN MSC/NAS TRAN INTERACTIVE TRAINING PROGRAM July 19, 1982 CON TENTS: 1. COMPUTER PROBLEM SET 2. SOLUTION TECHNIQUES 3. PRENASTRAN (INTERACTIVE PREPROCESSOR) 4. SELECTED PAGES FROM MSC/NASTRAN MANUAL (courtesy MacNeal-Schwendler Corporation) William J. Anderson

I

W. J.A. |COMPUTER PROBLEM SET | 7-6-81 Problem 1. Plane Stress Consider a thin sheet of metal under uniform tension in the x direction (Fig. 1). The sheet 4" has a circular hole which causes Y a stress concentration. - 4 — 100 mPa > (14500 psi) The sheet has a planform of -I x -- 1000 mm X 1000 mm (39. 370 --- in. X 39. 370 in.) and is 1 mm (0. 0394 in.) thick. The hole is -- 300 mm (11. 811 in.) in diameter. - The material is aluminum, with Fig. 1. Sheet metal under E = 68950. mPa (10.0 X 106 psi) plane stress. G = 26520. mPa (3.846 X 10 psi) v = 0. 3 A) Carry out a solution for stress, particularly in the vicinity of the hole. Reduce the size of the problem by using symmetry. Use isoparametric elements -- even two will suffice for a rough answer. B) Compare your answer with that given in Popov, "Mechanics of Materials, Prentice-Hall, 2nd Ed., 1976, pp. 48-50. The maximum stress is at the upper edge of the hole on the y axis. = K xx x max average K P/A effective For this geometry, K = 2. 28. Hence r x = (2. 28) (100000 N)/(700 mmZ) max = 326. mPa (47, 200 psi) The valuefor stress concentration factor is therefore 3. 26. 1

Problem 2. Beam. 50co 1 A slender cantilever beam isj 3000 mm (118. 11 in.) long. It is loaded with two concentrated / loads as shown in Fig. 1. I 200c) J The beam is made of steel, ('tq " 1.) with propertiesi bcc,-o, ICC,V^ -CCon E = 206, 840. mPa (30. X 10 psi) G = 79, 550. mPa (11. 54 X 106 psi) v = 0. 3 Fig. 1. Cantilever beam with concentrated loads. The beam is solid, with rectangular cross section (Fig. 2). The moment of inertia about the yy axis is 7 4 Z- I = 1.25 X 107 mm 14 YY (30.03 in4). 100 AM L, —~ 3.c37 3 and about the zz axis is: 7 4 JV I 2.81 X 10 mm 4 zz (67.57 in4) 15 C n - A) Carry out a finite element solution for deflection and stress in the beam. Fig. 2. Cross-section of beam. Free-end view. B) Compare the tip deflection and stress at the root with those obtained by Roark and Young in "Formulas for Stress and Strain," McGraw-Hill, 5th Ed., 1975, p. 96. W + 1.031 mm (0. 04061 in.) upward tip Co =4. 000 mPa (580. 2 psi) tension in top fiber comp. in bottom fiber ~~~root~~2

Problem 3. Plate. A uniform, elastic plate is fixed at a wall and supported by a piano hinge (drawn as a roller). Two vertical loads of 1000 N (224.8 lb.) are applied at the corners shown in Fig. 1. The plate is 10 mm (0. 394 in.) thick aluminum with properties: E = 68950. mPa (10.0 X 10 psi) G = 26520. mPa (3.846 X 106 psi) v = 0. 3 I000 N 22.48 Ib /2.6 I~ / / 100 1 00 [500 So o SO-"Fig. 1. Plate under concentrated loads. Fig. 2. Plane view of plate The piano hinge constrains all motion except rotation about its axis. Conventional plate theories do not contain a degree of freedom in rotation about the z axis. The engineer must therefore constrain this degree of freedom in a finite element program. Find the deflection under each of the loads. 3

?/ Problem 4. Solid. Consider a spur gear with "pitch circle" radius of 100 mm (3.937 in.) / \0 N/M and with 20 teeth. A tooth on the 346 gear is subjected to a line load by the mating tooth as contact is made \! (Fig. 1). This steady running load is modeled as a line load of 600 N/mm acting at the pitch circle radius. The load is assumed to act in the +y direction for simplicity. Fig. 1. Spur gear tooth under steady running load. The side view of the tooth is shown in more detail in Fig. 2. The tooth is "cylindrical, " i. e., it can be swept out by straight line generators. The gear is made of steel with -' properties:..-.. E =206, 840. mPa (30. 0 X 10 psi) 6..: —. -'. - 8.....:-. -."; t0 G =79, 550. mPa (11. 5X 10 psi) v =0. 3. A) Find the maximum stress in the -- -. - -- — o" gear tooth for the case of uni-''.-i: form line load shown in Fig. 1.- -I|Qo. B) Find the maximum stress in the - gear tooth for the case of gear u. -- -- ----- - misalignment, where the full/ - " 6000 N (1, 349 lb.) acts at one --:- --- - —' 9 side of the tooth (Fig. 3). ).. 10 C) Compare the two cases above J. -; to- see how much conservation -.. d needs to be included in the de- -. - o sign to prevent short term failure due to misalignment --- 2 and before the gears "wear - in. " ).-o..N' - 1%q Comparez e lo s e: 2 o A/ \ CoordjrKt2. mr\ Fig. 3. Eccentric load. Fig. 2. Contour of spur gear Itooth. 4

SOLUTION TECHNIQUES FOR COMPUTER PROBLEM SET Professor William J. Anderson The University of Michigan September 11, 1981

W. J.A. MSC/NASTRAN PROBLEM SET 7/1/81 Sample Problem #1 PLANE STRESS SOLUTION WITH QUAD8 I. Problem Statement Consider a thin sheet of metal 100 Pa. with a circular hole, and under uni- --- (14>l00 ) axial tension (Fig. 1). This is con- sidered a "plane stress" problem C because both the structure and the - loading on the structure lie in a 2dimensional plane. Little of interest happens in the z direction, out of the Fig. 1. Plane Stress Problem. plane. The sheet metal is aluminum with E = 68950 m Pa (10 psi) and v = 0. 3. The sheet is 1. 00 mm (0.0394 in) thick. The planform is 1000 mm by 1000 mm (39. 370 in x 39. 370 in). The hole is 300 mm (11. 811 in) in diameter. Carry out a solution for the stress in the vicinity of the hole. Use the CQUAD8 element, which is an 8-noded isoparametric element with both inplane (plane stress) and out of plane (bending) capability. We will suppress the out-of-plane behavior by constraining z displacements and all rotations. (Plane stress problems need no rotations at a node.) II. Symmetry Students have a lot of trouble with symmetry, although it is not a hard concept. Here, solve a quarter of the problem, say the first quadrant.. Note that the two cut surfaces each have a force and a displacement condition. Imposition of the displacement symmetry also removes the rigid body modes in this problem, but that it is an accidental side effect and does not always happen.

2 zero zero shear zero forces shear force force zero, -- displ. t00 m Pa (t+500S ps zero displ. zero shear force Figure 2. First Quadrant. Solved by Symmetry III. Mesh'.-. N\O OF VSGs/E t Let us use ( CQUAD elements. A crude, 2-element mesh is chosen. T, 43 A,\/- X D1rPLs CO rT?/TS7 ----— _______ 8 5 3( mm 1'K5nmfrx j,,10 \ 0 II | Figur 3. A $e No da To2 it C ht/\E o/ YDSaL. Co~c]DST[~A Z Figure 3. A Suggested Nodal Layout Chosen for Ease of Data Input

3 The x, y coordinates of the nodes are entered on "GRID" cards: mm inches 1) 0.,150. 1) 0.,59.06 2) 0., 325. 2) 0.,127.95 3) 0., 500. 3) 0.,196.85 4) 57.40, 138.60 4) 22. 60, 54.57 5) 250., 500. 5) 98.43,196.85 6) 106.10, oi6. 10 6) 41.77,41.77 7) 303.,303. 7) 119.29, 119.29 8) 500., 500. 8) 196.85, 196.85 9) 138.60, 57.40 9) 54.57, 22. 60 10) 500., 250. 10) 196.85, 98.43 11) 150.,0. 11) 59. 6,0. 12)'325., 0. 12) 127.95,0. 13) 500., 0. 13) 196.85, 0. The connectivity of the elements is entered on CQUAD8 cards: CQUAD8, 1, 1, 1, 6, 8, 3, 4, 7, 5, 2 J connectivity (norO\ oinlr CQUAD8, 2, 1, 6, 11, 13, 8, 9, 12, 10, 7 IV. Constraints We have two kinds of constraints in this problem: a) The plane stress condition that makes z displacements ("3" direction) and all rotations ("4", "5", "6" directions) zero. Use the GRDSET card for this: GRDSET,, 0,,,, 0, 3456 b) The symmetry conditions which cause some x displacements ("1" direction) and some y displacements ("2" direction) to be zero. Enter the constrained degree of freedom in the 8th field of the appropriate GRID card.

4 V. Loads Distributed loading on isoparametric elements calls for careful procedures in creating consistent equivalent nodal loads. This approach is beyond the grasp of the student at this time, therefore, a simpler "lumped" load method will be used. The lumped loads shown in Fig. 4 can be calculated based on the area exposed to pressure. This introduces a small error which will be discussed later. -- -- 1.25x \0' (zsID \t) 2.sox 10 N I~..b i.%( X I&o \13___i __ l.ZSxl01 1 (2o(0 1b) Figure 4. Lumped Loads This loading is entered on FORCE cards: FORCE, 1, 8, 0, 1., 125000., 0., 0. FORCE, 1, 10, 0, 1., 250000., 0., 0. FORCE, 1, 13, 0, 1., 125000. 0., 0. VI. Data File The above description covers the generation of the bulk data deck. When one adds the executive and the case control cards (perhaps through using the preprocessor PRENASTRAN) one should have the data file shown on pages 9 and 10.

5 VII. Execution. Plotting. Let us presume the input data for MSC/NASTRAN are in a permanent file named DATA and the output from the run is to be put in a temporary file named -OUTPUT. The run is executed by $ RUN AERO:NASTRAN SCARDS=DATA SPRINT=-OUTPUT Plot instructions (if requested) will appear in a file -PLOT, automatically created for you. A) To get a print of the tabular results, to be picked up at the Computing Center, North Campus (CCNC), $COPY -OUTPUT *PRINT* B) To view the plots on a Tektronix storage tube, and before signing off, $RUN *PLOTSEE PAR=-PLOT C) To get CALCOMP plots, to be picked up at CCNC, $ RENAME -PLOT PLOT $RUN *CCQUEUE PAR=PLOT D) The cost of this run will be $3. 00 to $ 5. 00 during the daytime (normal rates).

I_'Ente r't. erminal t. e LA36 MTS Ann Arbor (L.A36yLF26-CCOF' 00362) ~1S!G K4L7:llrE,;ter user Password.:I'T'RE:RM NORMAL y LN: V IlA*.LA ST S:(NOC)N WAS 09: 56 32:1 USEIR "IK4L7" SE(3iNED. ON AT 10:07:+03 ON SAT JUlL:8L/81:FRUN AERO: PRENAS]TRAN:tEEXECUT ION BEGINS NASTRAN FreFr rocessor by Ed Chmiel y1elen Bi..ius -anid Bill:. An ders onr:-t*.,**' VI'ERSION 4.2 *t*.**X * E e r Y for int r ucL r o duct in a n di i s rlt..::i.o s i. tr, e 1 s e t o s t a: r t' Y T'his Preprocessor proS r am allows t the iser to c:re-l;ate. (comPl.ete NASTRAN data file interactivel1 which is su:i.tabl e.For friiPut.;io AERO: NASTRAN or which can be puched orn ca rds'For e'ex!ec.ution' on a s'stem other thaln MTS i Suci7h a data file cc:ra-t bta s-nerTated for one of flour tl;es of problems: I - PFl ane St ress (us i nl CO( MEM or CQUAD8 elemenrts) 2 - Beam (usinri CBAR elements) 3 - Plate (Ius i ns CQLAD4 el:et, ents) 4 - Solid (usin rC CHEXA elemenr''t;is)'The Pro raml will Prompt the user Pfor infl''ormatl:i.on as ret.i..i red du..lrinril file xener'eat:ior)no There are ei-ht s;teps in.eere.ba.i %a f i:le: 1 Ex:ecutive Con'tr'ol Dfeck -- user su.P1plies Fr'oblem i dent i i' i cat i on 2 ) Case C ontrol DIeck.. - user su. j lies t.:L t 1. subt it le a-nd ilnuLTmber of subcases 3) G r i d F' o i D f t e:i. rn i t i o r s - r..r i,e i e s r. d P..o:i. in t i d / s I-Iand coC) o rdj. n a t'es 4 ) Grid F' o:i. n t Co n s t r a i n-t s use-. s e r u P P 1:i. e s r i d p o i n t i d / s anrd con st ra ti lt c'odes 5) El ement Def i n:t i tions - usJe. rs uFL 1. iS es eleme intr i d's p roperti s i id', a'd, rid'o i rt cor, r ec: t i. vi i' es 6) P'roperty D efi ni J. tions - user. suF:L ies materi.T'. id's anrd t, hic.'kn eesses o r sec't i ol - ro:oPe rt,i es y if' rec Q< i red 7 ) Mater i a 1 D e i:i. n s - u s e r s u P P 1 i e s Y o,..u' s Mod.. u 1 u s anld F'oisson —'s, Ratio 8 ) _Load Def i ni -Lioi nos u- u se r su.F P1.:i e s r i. d Poi n i d /. rd C a r't e s:i a n c o mr:-. o n e n t; s f,.' t':,o r c e-.. 6

R ed,Pa. r o r a i mess a e asd p' r o m iI r;, c e r e f u.1. 1. be ore e: n nr rI t er:i. n l-. ari n in format ion + Refe reri ces t;o the NASTRAN UL.JerT' s M:anu. Vol 1 are enclosed i-n Pa-renrtheses.? e-.g (24-30()7), If a messai.e is not f'ul.l understood? refer to the a PProp: C'ri at e Pa:Be0 If for some reason the Procram terminates before all ei.'ght steps of file senreration have been completed it c(anr be re-run usin r the rog ram' s Restart capab 1 i it, This allows the user to contirue erner'atinr the f:'i.e without hav:i.ng to re-enter data that has beer sluc(cessfulJ].l erntered in — th-e f:i..le, There is also a Rewrite capablili t whi. ch allows the user to editt the file It is e xtremirel:i.mporY'tan1t that the t.yp e of P r o b 1em t te o f. elemenr t y arnd t h e ni..i mb e r of.r id ro i.r ts a nr d elements be the same in Restart or Rewrite mode as itin Create. mode Choose mode of operationr 1- Create 2 - Restart 3 - Rewrite? Ex.ecutive Control Deck. Enter Problem Identification (2.2-1)t?'ANDERSONAER0510 $ JULY 18p 1981, DATA IN K4L7:ESCDATA8. Enrte r roblem type( 1=F'lane Stress 2=Beeaiyr, 3:::F'l'lat e 4::Sol: d) S? 1 Choose element tyPe for FPlane Stress: 1 - CODMEM ('for stress contour p,:lots) 2 - CQUAE). (for better acc.Juracy with f'ewer elemenrts) 2 Case Control Deck Enrter Title (2.3-137): S ESC FROBLEM *1 Enter Subtitle (2.3-117):'? PLANE STRESS. PLATE WITH HOLE. From one to rnine subcases m.a be rea..ested to allow applicatiorn of different sets of loads to the model~ FIor' each subcase, the user will be:r'rompted for' a set of loads in Step 8 of file Ysenerationi Enter des:ired number of subcases (2 3)-55 anrd 2.o 3-111).? 1 GRID Grid F'oiint D.efin itions (2.4-205) Def ine each. rid point by er',t;er:i.n.:i. t.,s id ni.umbe'r arid c o o r d:i rn a e s. R. e P e t p r o e s s f o r e a c h. ri d p o i nr'it I'F a. i e r r o r was m ade in def ini - i nd a Yr:i d Po i nt tohe same i d ar' d o r rect; coo r d i rI a ates ma ts be r e - e n t e r ed t hee' r e by re. a rc i i. i the o. d i nf'o r ma t; i or i n th e data f i 1 e When a 11. r i d. o i.r n'tv h ave beer def ined? e nt'ler $ ENDFILE or C TR L. —C t;o stop The r'.ro r.- r.am r w:i 1l ~t hen check.:if a].l rid poirlnts rnl.iumber'ed from or'le'eo t'the user' d e f i, e d m r x i um u m hve beer, de' fi:n e d E i t e r m rix:i. yi..m' r' i d i d 7 13 M. a i. ml G r i d id 13 s Y lease corisf i'ir m(Y or N) *'?7~ Y ~~? 7

Fo r: each.. ri:Ld P o:L r't e nt?. e r 1: ID X y Y? 2yO.y325. T? 300. y500? 4y57+41358.6 SyV250o5yO. 6Y106 1.?106+1? 7,303y303 T 8,500.Y500o' 9,138.6,57.4?:LO y500v 250O 11,150.,0. 1 12y325..O:. 13,500. OO0 ~? $ENDFILE Define t hle te costrainti s fos r a.1 rid F;oi nr ( GRID 2.4-205)+ Enter grid numberconstraints (anrY di:its 1-6,s no blarnk.s). Enter $ENDFILE or CTRL-C to stop.?:L 13456,Y 2,13456 ": y 13456'I' 1,23456'.I 12y23456' 13S23456 $ENDF1:LE C Q U3A D 8 El1 emert. en e f' i nit'i ons ( 2 4-1 J. 12 c ) Define each element bv enterTi.rin its id nu.mber and.Srid id's specif lin r the elemenrt's cornnleec:tivit.t If an error was made in defininr. anr element, the samle id rand correct conri'iectivi.t may be re-entered thereb. re.:l1. ac:inr- the old in'fo lrmat:i Lo in the data file. When a:.l elementsl have beenr defined, enter $ENDFILE or CTRL-C to stop, The pro<'ramt will ther chelck.:i.f a 1 elements numbered from onre to) the user defirned maxts.1 ITmum have beer def i nedr En te r max;.:i l lmu element i d'? 2 M.,iaimum eleiment id = 2y please conf:irm(Y or N)*? Y For each elee el e rt enter EID GIDi..G + + GD8 where GID1 thru GID4 define elemenlt co rners; in counrter-clock wisei order, and GID5 thru GIDI de:fiLne midrpo:i.nts co)f elemernt edi..ese in counter-clockwise order NOTE - Follow irin def i. nition of th; e f'ir>st element there will be a PromFpt for' roperty id number and an option will be iiven to either use the sam;e ProPer't;. id for all elements or to enter ea Irnew id, for eachl el ement En:'l t e e e m e n t di a t a? 1 Y y 6,8y 3,4 7 v 5 y 2 Er, ter F. rop:e'rt"L.i d:'? 1 Is.iven:. tr'er tyd id to be.used for all e leme ts (Y or N)??' y Enrter e lemenit data..? 7 2 )76,l1 3,, 89?.1.'. 2?y 7'. $1': N I F':i8 1_ E

F'`SHELL i:ro'.ertes Defiritiori ns 2. 4-326a) E r, terT mater: ri al id an-ld th ic n..r ess fos r T' P r e r t:i d I? 1?1. MAT1 Mate'rial Definitio ns (2.4-.215) Enter Younl.' s Modulus, Poisson's Ratio' for material t1? 68950..3 F:ORCE Load Definitions (2.4-179) 1 load subcase(s) have been reauested Define loads'or subcase 1 For each loads enter grid id defiinig proint at which load is to be appliedy followed b'y XYPZ coirl')rnen ts of force? ie. GIDIFXFYYFZ. To delete a load? enter -GIDY00Y0.+ I-nter $ENDFILE or CTRL-C to stop for:iven subcaset 8,12500.,0,?0.? 10,25000.,0.,0.'? 13Y12500 Y0.y0' ENDFILE Enter 1 to List datafile or 0;'P 1:$LIST DATA 10 ID ANDERSONAER0510 $ JULY 18Y 1981. DATA IN K4L7:ESCDATA8, 11 TIME 1 12 SOL 24 13 DI AG 14 14 CENI: 100 TITLE=ESC PROBLEM:1 + 101 SUBTITLE:= PLANE STRESS. PLATE WITH HOLE. 102 ECHO:=BO'H 104 DISPLACEMENT=ALL 105 STRESS=ALL 106 ELFORC E=ALL 107 S"' PCFORCES: AL L. 1. 11 SULBCASE 1 111,5 LOAD-= 140 0 UTPUT I ( FPL T ) 141 PLOTTER NASTRAN 142 SET 1 INCLUDE ALL. 153 AXES ZXY 154 VIE-W 0.00.(Y0 155 FI: ND 156 F'L OT LABEL B:)OTH 157 I'LOT STATIC D:EF'O'RMlATION 0 SE: 1 190 BEGIN BULK 191'PARAMH AUTOSFPC Y ES 192 G RDSE'T 3 456 2 01 GRID 1 0.0 150.000 0.0 1.3456 202 GRID 2 0.0 325.000 0 0 1.3.456 203 GRID 3 0.0 500. 00 0. 0 13456: 204 GRID 4 57.400 138.600 0.0: 5 GRID 5 250.000 500.00 )() 0 206 GRID 6 106100 106.100 ( 00.: 20 G7 GRID 7 303 )000 303 * 000 0 0 9

o O ~ < -..'0', ~r <,, * * o~ n'i. 4 4' - c _ - 4 000.03 3. - _'> * * o * + 0 iis 4 4,4 4 * 4 4O 4, 4 L~10 r,-r YL-a 0 0 Z C 3 3 ~.13 _ e n'0O. O O,^-4 0 ^ O O' 000.000'0 0 00 0.00 - 4,r 4 4 44. 4' 4 4 - V't' OOOOOO 0 3 0.0 0 0 " O. s0 0 4-:'!. 0. 3 o cn C ~.i - C- ~ " ~ O1 c 0,~- 4.O -'4 0 r-i i —.", *< s: -4 0 <:...- - <. a - S'*C- 0) 4 i. 4, 4 e- _- - -'_ s 2:Z; Li_. <z. a "._? -'s s".. oo-4000oo col-3 U?-.i iN C )j _S Zr Cl t - _ C Z 3 T4 S_ 0. 2 ^ 2:.- 3 r( *;- s. 2: _, i: 0',."-.""..".'"-. - -.''.:.-..'"'-..'".'-..'..' "...'". -,',;.'. 0.. * J, 10

0 0 n <Z <Z m~~~~~~~~~~~~~Cr C ~~~~~~~~~~y I rsN~~~i.,'.4 em~~~~~~~~~~~~~~~~~~~~~~~~.S~*,0 SO',0,,..0 SO SO lo r 0. ~L~ IT IT'I T'-.'.,. 1-4.. (.' I'N. <e.-4 (1.. 03 I < <; w~'.1 000 000 000 00 1 1- f........ 000 000000.4',.'' 3 ^ O O O O I 1 i' 3.. 000 00 00000. j ^ un in............. I0_. i- _ ooo 4 O 20 -4 N o AA =*-^ c- -^ r uzs -aI u^ 4 000 0 -00 02.. W1 >-. 0 F*J u f.1 C0 3 00 ) O i O (Lt3Am A t _iO O t3 0 r:3 0' _ U!I'-. Z Zj i i M _".L'.'..:.-.!o;- to X. < W 4 W i - W. -: t 3 t1.E - C — O.. J ii 1-4 4 CD <rj S: zX n S - ",-,t 13]f, 1i% T-rC i' 4n ooz" i ~U if 1 4 1-CJ H _'' _J E T t J _ C + _ L / _ Sn U' _ ~ J 0 0 J U _ 0 ^ a 4 _' _ < 4 4 \ w u r C5 -. -A s" -Z X v - O j o^ ow r r! i _S i f.d j I e sL C 3 a< l;L=Cj ~ 4'Ji i3G I_;s Qn s 5: OC C:-~!- H * ~ 3 s I i i W a: X X i3 G f n <L > u u, CL. f3: z.: s:I s U^ U t (3 fJ; U C + * + -J'.! o* t4? ^ <r ot < r O j *r tu so S <i T- o ^ sP tr^ ot3 r- c iN 3 > 3 r- r-j r^ r, T-{ c4 -j i- T ~ T T < ^ H i ~ si T _ *~- T- st q -i ~~ *-t (r 4 -. i c -4 e4 -4 s*4 ^4 Ct s C-4 nC- N 1.1

c, c C'~, ~ - 000 z: = 2: Z 2 Z Z Q2 2:2:2 X:2 m2 2 r: Q:2 * ** - CDCDCDCDCDCC CDC C.DCDCCDDCDCDCDCDC.DCD O O O.X r. (&'~.i ~i. i"i.'..D ti —ti' L i Fi.D it i{ Xi~, r.O 0i3 it c_~'30~^ i*-t3:-~'3~~1 0 4 [>i\ i',?.~ CfX XX ii > Q iN N4 _ CM C C Jl C I t -4 4 1-4 iv-4 13 0 v- $n C - in. j w 0 o. f0 O O t'3 I0 iw o,( $n X H X 3 Ur,I T! 4 r4 rt i Ci r > 1 r ~ { r t r000 if'-~(O0 CS) 0r ~ 4> C OOO 003 I:3 ~.- ^ _ s u cn ^ ^ ^(^!. Z <Z 0*;< 0 0.-:CS 0 00 Q= Zt m F m *r tI <rl-( T- C W C0 m } (. N J 3 O - \ _ i CO I- i0 0. Z uj o o io o oouo 3 i i c 0 U * -c *M} 0 o r o NSN X MX W N V M n M A N uC-> M e 0 0 3 0 H,4 n ooooO imITr00,.-'. NO N N3 0,, ^c s o VT', Q | H j j C_ <r i~ LU uj W tJ o - CJ tt N N -i - CO<oi ~ CO 0~4 r-4 W ~ } L t 4 7 ~ 8 ~ ~D f s3 O O hU 0 0 3 0 O O O 0 3 tO ri UI O o o o o O OO * U u <* U-n Z ~57 * i(j -' +r x o sj C U U iz4 co 2 Z 4 3 sM V Z j x o < w d c 3 0 m 5c % - N) U.) 54 c j+4 qo i-4 s oooI-4.T-6<^2 2 0o o o.i. * * o * * o * * * ~ * * * ~ ~* ~ ~ 4 * ~ ~ ~ = tI- i^ s_ y Ct t r =I * o o o o jO L.f S ^ - 4 4 i ^ o ro 3 rN Nt N ~r m ms m m t v3 iqr * o i " < Z L d * 0I <.: 04^ +* 2 Z Z "'4 0 "j ^ 1 0 "n 3 3 0 - 0 - r, 0,^.s 0 i3 23 O,', 30 i _- tv-a CU 4 s o o s _i. m Cl- %O M Is OC. C-4 0,;. r CO T; ^, rN co q t) 0 ilo _,s T1: O 0 Zr 1"~S * 3 -l - - - CS ii _:'' <Z Z U,)O <. cj I! t IZ — i 1:0C J. J_ &L ". —,r,~.!,', 30.. 0'..l'.0.- -1~ 30~

z z z z 2: z z z z z z ]3Bl Cfi {Si jrt C'i CIt {Q i CQ irS &i S L0 r,.DxC Cr7f..,. D r..D C. (D CD, x J E,:.I-r C -. - CD a 3 0C 0 O "3'3' fC LCL3 W i- 4 cJ.v4'.41.4 ('J (J (' (= ('4 *-, J O3 <: r a) La(.0 La CL. ll.' XLL,l~h-h M w LID N N M 0 ZZZa',0,OrQ.. W W O~nS LI W Wi LL.LL.LL.LLj Q Z ",.! 00 ci3 00 CX)-i in E (n;, W 04 w o o, 3,, ~-.. U U O _O'..O-,IN 1, 1 3 t O *-'b'', 0', b. L L. L i Z,,C,,, 2: _ 0~-0 M O M M O N M IT T- T-I C-4 C-4 w co-4 (4 IZZZZ X W ZOc C C W z WW W3 W Z i. t-G. C'- L- If4'O'' 4 04:4 000 CD 4 04 04 i04 iL O 0 Z N_ N',O'IO,0S N O r4 W z c W3 WW 0 (L Co 4 LO 4 O. Z 034 cl c'I-O' -' *- - s -, W. <I: I - 0000000000 COe En W- WI- C.'.""'.-l C-4 ******. 020200 2-0 Z0 C <G U0<E: -1O o s- WZ Z 2, = 22 s - b3i Ld L C3 0 3 r.3 ~ i-. Ld C-4 M + M ^ 0S 4: UH N O N 0 i% 0T4 M M. M M C @ U'En NNN COWCOWO0' Ul 1 S O3.O I-I-I — C-)Q:. o. -4C f *W _- _0 Z. L )O.3~ 04 O ^O O OX zL.Z Z2W-4 ZI *.C'4 q' W -- W I-3~33 COCOO'. COO".... U''. CL. ('1O('NaNaC"1C M 00 II.ii0 I- CC L O,, 0, S- -.Z *. 0Z * LZ o n Z -.I N' 2222 O 2 CZ O C L0. ~.3'WO,. }., ", U H ZI oor~.. Z W -b - - --.. - - s.- i - O - a L.J 0s. OMd O.!U",Oo.-. ~,,O l rssIo oO. * C. O O 303::: o3 C..0 CL.-'W ('.4 n. bS b' l. O, r, o ".,q.. +,3 X P + 4 _ tH * 3 * lK!' i:,,...1.'-r i — ~tft tci "o't,-J ut,!- 03 ~J ^}. ZZ *~*'J -. 2::,,.::: n o ~ t t t 4t t Sn S oo C M2 L3 2: QD Z i Z 4S2 a:Q + 0 4 *~u -I ij 2 I - O O 3 3 3 ++ OO 3 LO -: ~ o t H o a W J t LJ - i,'r'03 0cn,c 3 u X L'. s't'. Lul teX s' "'2" _ X X r='o:>:>::5. 2J - * ---.3 2';_ * i < n.' L 3 i. C = O..'t..'.,.,:'.., i,: ~. * # * o 13 8y J7 1 (2-rt

1 7/18/81 co 0IQ8 / o8 0/ 410 ~~~0-~~~~~~~~~1 C5 Co.4 -, ------------------------------- 13,9- FOM SVO PLATE WITH MOLE. < I UNOEFORHED SHfAPE 14

2 7/18/81 MX-E. m 1.20617580 2 I 3 ID STATIC oE. OA E I LOA Q- E^BtrafisSl! P~LTE WITH'ILE. 0. STI:ITIC ClftOR. USBCr1E 1 LOP,0 1 15

W.J.A. 7/25/81 MSC/NASTRAN PROBLEM SET Sample Problem 1 PLANE STRESS SOLUTION WITH QDMEM B LS I. Procedure 8 Let us solve the plane stress problem with the QDMEM ele- ment. Although giving less accuracy than the "in-plane" solution from the QUAD8 element, this element allows the generation of contour plots in MSC/NASTRAN.* Use the nodal coordinates: t 9 Fig. 1. Plane stress problem solved with 8 Q DMEM elements. GRID 1 0.0 150 000 GRID 2 0.0 325.000 GRID 3 0O0 500.000 GRID 4 57.400 138.600 GRID 5 250.000 500.000 GRID 6 106.100 106.100 GRID 7 303.000 303.000 GRID 8 500.000 500.000 GRID 9 138 600 57 400 GRID 10 500.000 250.000 GRID 11 1 50.000 0.0 GRID 12 325.000 0.0 GRID 13 500.000 0.0 J GR I: El 14 191 200 325. 000 GRID 15 325.000 191.200 In MSC/NASTRAN Version 61, contour plots are not available for many of the more advanced, isoparametric elements. 16

II. Results A) What is the stress at node 1, where the stress concentration occurs? (This may cause you some problem in interpretation. Melosh once said that for elements yielding only one stress, it is improper to interpolate that stress to the boundaries.) B) Compare your answer with that of the QUAD8 element and with the recommended solution. What do you conclude about accuracy of "straight-sided" versus isoparametric elements? 17

: i: I... I1:'ST C 1D TA.1. 3.::-:1. 0 1. A N I:it:liE R S 0 {' 2 ( C F;, ) 51 () r5 1 $ J U _L Y 2' s:19 1. 7 9: I N K 4.:C D AT'T A;;,::>...1.'T 1 I-':E 1 > 1 ": 0l1...'24:: - 1.3 II::.14:: 1 4 C IE Ni >:: 100'T I'LE.: ESC'PROBLEM i t 1..::.01. 0. SUBT ITLE E F:'L. A NE S' RES7 S SOIUT:ION W:1: TlH QDX:ME M.:: 1 02' ECHO:=BO'TH:> 104 DISPI::'_ACEMEN TAL I 1:.05 STREl S:= A Il:: 106 ELFO RIRCE:=AL.I > 107 SPCFORCES:-AI....::'> 1.1 SUBCASE 1:: 11. 5 I._OAD: 1: 1.40 C):'.J TPUT'... OT):: 1.4 1. FIL TT ER NA S TI:; AN 142 SET 1. INC:L.UX:E A.1L....: 153 AXES ZXPY:: 154 VIEW 0 0,, 0.::~ 155 FIN:N 156 PLOT ILABEL. BOTH::: 157'LOT STATI:C DX:IEF'ORtMATION 0 SET 1 158'I.P OT SIHRINK 159 CONTOC)URM MAJP,'RIN MAX: 160 PF'_ OT CONTOU1R SET 1. OJ TIN:NE::" 190 BEGIN BULK >: 191 FARAM AUTOS'PC YES: 192 GRDSET' 3456::-'201 GRID I 0 0 150. 000 0.0 1.3456 202 GRID 2 0.0 325.000 0.0 13456:: 203 GRID 3 0 + 500.000 0.0 3456::- 204 GRID 4 57.400 138o600 0.0 2: )05 GRID 5 250,000 500.000 0. 20 G6RID 6 106.100 106.100 0.0 207 GRID 7 303.000 303.000 0.0::-> 208 GR ID 8,.500.000 500.000 0.0. 209 GRID 9:138.600 57'.40 0 O0:: 210 2GRID 1. 500.000 250.000 00( 211 GRID 11 150 000 00 00 0 23456 212 GRID.12 325.000 000 0 0 2 3456 2:: 1 3 GRID 13 500.000 0.0 0 0 23456 214 GRID 14 191 200 325.000 0 0::' 2:"15 GRIX 15 325.000 191.200 0.0 301 CQDMEM 1 1 1 4 14 2::: 302 CQDMEM 2 1. 2 14 5 3: 303 CQDMEM 3 1 4 6 7 1 4 304 CQDMEM 4 1 14 7 8 5:::. 305 CQEDMI IEM 5 1 6 9 15 7 306 CQDMlEM: 6 1 7 15 1.0 8::: 307 CQDMEM 7 1 9 11 12:1.5:: 308 CQDMEM 8 1 15 12:13 10 401 F'QDMEM 1 1 1.0 5:MAT 16 i90E A 0. 4 0. )3:: 708 + 1 FORC E 1 8 0 10 1 2500. 0. 0O: 7:1. 0 ~ 1 FORCE 1 10 0 1.0 25000. o 0. 0:: 713.1 FORCE 1 13 0 1.0 12500, 0. 0.:: 1000 E' NDDATA:IFNLEND OF FILE 18

:I:RUlJN AERO'NASTRAN SCARrDS=ESCDATA'13 SP'RINT=ESCRES113 I:EXECUTION BEGINS EXECUTION BEGI::NS:$CONTRO(i -.**4NO 1 S: ZEi:. 000 ~:$CONTIROL.. -.i#t#N02 S I 1ZE:i1. OF':4$CONTRFOL —'.:lIl:NO3 S I Z'I-O. 00F:$1CONTROL. -**:IN04 Si::ZE=IO. O(P EXIECUTION B E1 GINS NAYSTRAN L4 OADE:D AT LOC:ATION 700000 CPU LF' I/O EXCF' COUNTS'TriACKS SEC SEC SEC NASTRAN TOTAL U SEDi 0 o o 0 0 0 0 S H. BEGN 0 212 8*0 05 11 0 16 XS1 0 0 452.12,0 1 0 21 0 32 IF 1 0*650 17,0 2.0 41 0 52 XGP]: 1*445 26*0 3*1 63 0 60 9 GP1 BEGN 1i575 28*0 3 +9 78 0 76 11 GP2 BEGN 1 657 29.0 4 1 83 0 76 1.5'LT HBDY BEGN 1 703 32 0 4 3 85 0 76 17'PLTSET BEGN 1 813 35 0 4 9 98 0 92 19 FRT'MSG BEGN 1i856 35.0 5, 10 100 0 92 22 PLOT BEGN 3*346 43.0 5* 7 114 0 92 23 P'RTMSG BEGN 3 422 45,0 5*8 116 0 92 26 GFP3 BEGN 3 535 47,0 6+4 128 0 92 28 TAl BEGN 3. 817 53*0 7 4 149 9' 31 EMG BEGN 4 372 56 0 7 8 157 0 92 37 EMA BEGN 4 917 66 *0 88 175 0 92 97 GP4 BEGN 5 046 67 0 9 4 189 0 96 99 GPSFP BEGN 5 *155 69.0 9+8 197 0 96 110 SCE1 BEGN 5 333 72*0 10 4 209 0 108 140 RBMG2 BEGN 5.506 75*0 10+8 215 0 108 145 SSG1 BEIGN 5,686 79*0 11 4 228 0 108 148 SSG2 BEGN 5 867 82 0 12.1 243 0 120:1.55 SSG3 BEGN 6.100 86 0 12*8 257 0 120 160 SDR1 BEGN 6.477 92,0 14 3 286 0 128 176 SDR2 BEGN 6 793 97.0 15.3 307 0:1.28 185 OI-F' BEG' N 6 882 99.0 15.5 310 0 28 12 86 SDRX BEGN 6.960 100.0 15.8 316 0 128 188 O F' B EZGN 7+.117 104,0 16.0 320 0 2 12 8 1.89 G'FDR BEGN 7 253 106 0 16 2 324 0 128 190 OFP BEGN 7.296 107*0 16.3 325 0 128 193 OFF' BEGN 7 465 111 0 16.5 330 0 1.28 205 FLT SET BEGN 7 581 113*0 17.1 343 0) 1 E8 206 PR:'TMSG BEGN 7.610 114.0 127+ 3 345 0 1.28 20 7 PLFT BEGI- i N 9+683 125.0 18.8 3.77 0 128 208 PRTMSG BE GN 9.758 126 0 18.9 379 0 128 214 EXIT BEGN 11-10O02 STOP 1. VSSCMD238W NSTN04'- O::N DSDEF BEING C..OSEDI VSSCMD238W N ST'NOt)3 -- IOPElN DSDEF BEI1:NG CLOSED YSSCMD2'38W NSTN02 -— OPEN DSDEF BEING CLOSE:'D VSSCMD238W NSTNOJ. - OP'EN DSDlEF: BEING CIL...OS:ED -P D S P F' C T ISC:''' I 0 I)N G I:E N E R A' 10 N B l G 1:: N; l: IXECt.I O':1: N'TER;M I N ATI-: D 19

41:L. ESCDlATA13 FP'RINT* >t::':RINT;' A:SS1: GNED R1CEIPT1:' I NUMTBE:R 603206 IR' PRINT* 603206 RELE:3AS El TO' CNTRI 3 PAGES. IC ESiRES13 f':PINTS::>:.F'RINT* ASSIGNED R1EC1EI:'T NUMBEER 603207 IlF'PRINT* 603207 RELEASEE TO CNTR, 19:PA'GES:tlRENAME -PLOT ESCFL.OT13 4I:ONE:RUN SICCQtUEUE PR:'ESCP'LOT 13:ITEXECUT ION IEG INS 4 I PL OTS PLO'T' TNG RIQUIRE 237 SEEC AND 47 II: N 7 PEN WAS LJUP 45% OF TIHE T IME. "' 4L7: ESCPLOT13 " IA E BE PN F'E:M I TTTIED " 1:i' 1K IE<IY:=-(EI CCQ.lF UE "'LCT ASSIGNED RECEIPT:11 72 4300.:11fEXECUUT ION TERM:INATE:D 20

I? 7/2/81 0 0 o snt r —(0 0 ~" 0 N r^\I\ ~ ~ ~ ~ ~ z.4-.RNyTBM ELUTIO WITH OWMEM 13 a>UlOLREFOED S3HIPE 0< 21

2 7/24/81 2 oo CN 0 0 0r 0 ~~~~~~~~~~2

3 7/24/81 RX-DEF. 1.06167880 3 1 00 CD 0) V 0I N 7''I.EPh ^fSS^l iOTiSN HITH OMM 5 O RSTTIC DEFAM. 916CRSE I,1 - ) r) ~) G) STATIC DEFOR. SUSCRSE I L0 1 CI) 23

/ 7/24/81 00.J I CN'- 0) \ o 7 X); ) 1 0 RTsflrC STRESS SlE 1 LSIO I -0

W. J. A. 7/16/81 MSC/NASTRAN PROBLEM SET Sample Problem 2 BAR SOLUTION I. Problem Statement +. A slender, elastic, cantilever / bar 3000 mm (118. 11 in.) long is loaded /ak X _ _- l.. as shown in Fig. 1. The beam is made20 N of steel with properties: 6 / (. E = 206, 840 mPa (30 x 10 psi) I );|-IOOOM-** —4o000+- l000 -> v = 0.3 (3 The beam has a solid rectangular cross section as shown in Fig. 2. Fig. 1. Cantilever bar with The moments of inertia about the concentrated vertical loads. axes are: 7 4 4 I = 1.25 x 10 mm (30. 03 in) yy 7 4 4 I = 2.81 x 10 mm (67.57 in ) I 1 zz eAY7 It is desired to find deflections and stresses in the beam and to plot I5o 150 ~w _M deflections. Fig. 2. Cross section of bar. Free-end view. II. Solution The BAR element will be used. Three elements are the minimum permissible and will be used. In this case, with concentrated loads at the ends of each element, the solution should be exact. (The displacement field is modelled with a cubic shape function and this exactly represents a beam with concentrated loads at the ends.) All finite element programs using beam-type elements have a problem in determining orientation LA of the principal axes in space. In general, the beam can be at some oblique angle to the coordinate.sys- - _._ (_') tem. In MSC/NASTRAN, this is solved by defining a "Plane 1" ^- — / passing through a user-specified GA GB vector V and through the beam axis (Fig. 3.). The vector v Fig. 3. Coordinate system. 25

is centered at the left end (the GA node) of the element in question. We will choose a vector v = (1., 1., 0.). For the left element, A LA 2. the global coordinate system coincides with the local system. I Bending in Plane 1 is controlled by the stiffness I. I ZZ NASTRAN designates, in this case, i = I 1 - zz Plane 2 is defined as passing through the beam axis and perpendicular to Plane 1 (Fig. 4.) Fig. 4. Definition of Plane 2. Bending in this plane is governed by I and I = I 2 yy The most general case of beam orientation is shown on page 1. 3-17 of the User' s Manual. A proposed finite element model is shown in Fig. 5. The grid cards are: GRID, l,0.,0.,0. GRID, 2, 1000.,0.,0. ( 2 3 @> 4 GRID, 3, 2000.,0.,0. GRID, 4, 3000., 0.,0. The stresses can be read at four points. We will choose the four corners of the cross section as shown in Fig. 2. The property Fig. 5. F. E. model of beam. card will be: PBAR, 1, 1, 1. 50E+04, 2. 81E+07, 1. 25E+07,,,, +A1001 +A100 1, 7550, -7, 5 5 0..,-75., -50., 75., -50. such that stresses will be printed for the four points on the cross section shown in Fig. 6. 3.... - Fig. 6. Stresses given at four points on cross section. 26

i — i i U0 *0 0 - O cx~~ ~~~~~,~~ ~C, O 0 0 0 00 *' <2: -i vs — 4 ~ s 0 0 v 0 O U~ Zji i.. 0000000 0000 -0 - h- CZ~ 4> L,._~ 0 000 -" -0 0 _. g-,,, 000 " -'3 -..; + 9 * + 3. 3 = -., 2 " 0' 0 "" 4 0 0 >C Q Vo o Z -I Z 0 0 ^ I:j3 S 5i 0 *'4 r ~ -14 T< rS T-< - 0 L - C *. * Z O ^!. _I isr- jS iS<: SU- 0 W;' ri >- * = b O L-J S <S tiL 0 U _- ) J- Zt Z u La UZ i 11 O-Z i-.,<r.sF. ej Z * 6-:- CO W... 0 W n- S t S E. Z oi 4 ii' - <rt t j t u.< tn i. u. S,*:_4 _ ii iS I v * 4 t-4m.iS- -i L -- -'.... ". 0'"-........ ~.i................. i........-.',.... 27 ~4-::Z.L.......

RUN AERO:NASTRl-AN SCARDS=ESCDATA6 SPRINT=ESCRES6:EXEC.T ION BEGINS EXECUTION BEGINS:I$CONTRI:OL -# *N01 SIZE1 0OO':i$.CONTROL -#*$N02 SIZE=1000P1:$CONTROL -:lN03 SIZE=1.000 F:$CO(NTROL -**N04 SIZE=1000P EXECUTION BEGINS NASTRAN L.OADED AT LOCATION 700000 CPU ELAP I/0 EXC P COUNTS TRAC(KS SEC SEC SEC NAST'RAN TOTAL USED,0* 000 0*0 0 0 0 SE BEGN0 163 6.0 0*5 11 0 16 XSOR 0.301 7*0 1,1 22 0 32 IFP 0.477 10.0 2,1 42 0 52 XGPI 1.248 17 0 3,2 64 0 60 9 GP1 BEGN 1.366 19.0 3.9 79 0 76 11 GF2 BEGN 1.442 20.0 4.2 84 0 76 15 PLTHBDY BEGN 1*485 22.0 4.3 86 0 76 17 PLTSET BEGN 1.585 26.0 4.9 99 0 92 19 F'RtTMSG BEGN 1,622 27*0 5,0 101 0 92 22 PLOT BEGN 2*373 31.0 5.6 112 0 92 23 F'RTMSG BEGN 2,421 32*0 5,7 114 0 92 26 GP3 BEGN 2.526 35.0 6.3 126 0 92 28 TA1 BEGN 2 777 40.0 7.3 145 0 92 31 EMG BEGN 2.958 42.0 7.6 153 0 92 37 EMA BEGN 3,431 48.0 8.5 171 0 92 97 GP4 BEGN 3.545 49,0 9.3 185 0 96 99 GPSP1'BEGN 3.642 51.0 9.6 193 0 96 110 SCE1 BEGN 3,802 54.0 10.3 205 0 108 140 RBM'2 BEGN 3.956 57.0 10.5 211 0 108 145 SSG1 BEGN 4.122 60.0 11.2 224 0 108 148 SSG2 BEGN 4.290 62.0 11.9 239 0 120 155 SSG3 BEGN 4.501 66.0 12*6 253 0 120 160 SDR1 BEGN 4 840 70.0 141 282 0 128 176 SDR2 BEGN 5.123 74.0 15.1 303 0 128 185 OFPF BEGN 5.182 75.0 15.3 306 0 128 186 SDRX BEGN 5*252 77.0 15.6 312 0 128 188 OFF' BEGN 5+396 80.0 15.8 316 0 128 189 GF'FDR BEGN 5*526 83.0 160 20 320 128 190 OFP BEGN 5,566 84.0 16.0 321 0 128 193 OFF' BEGN 5,723 88,0 16.3 326 0 128 205 F'LTSET BE-GN 5.831 90.0 16.9 339 0 128 206 PRTMSG BEGN 5.857 91.0 17+0 341 0 128 207 PLOT BEGN 6,765 97+0 17.7 354 0 128 208 PRTMSG BEGN 6.803 97.0 17.8 356 0 128 214 EXIT BEGN IH0002I STOP 1 VSSCMD238W NSTNO4 - OPEN DSDEF BEING CLOSED VSSCMD238W NSTN03 - OPEN DSDEF BEING' CLOSED VSSCMD238W NSTN02 - OPEN DSDEF BEING CLOSED VSSCMD238W NSTN01 - OPEN DSDEF BEING CLOSEI:D PFDS: P, LOT DESCRIFPTION GENERATION BEGINS. *EI..XECUTI: ON T'iRM I NA TED: $5.29 $601:LT 28

*L ESCDATA6 *PRINT*::.*PRINT* ASSIGNED RECEIPT NNUMBER 600048'*:'RINT* 600048 RELEASED TO CNTR, 3 PAGES.: $+15, $6*15T' *C ESCRES6 *PRINT* >::*(PRINT* ASSIGNED RECEIPT NUMBER 60005:1:*F'RINNT* 600051 RELEASED TO CNTR 17 PAGES+:t $+27, $6*42T *RENAME -PLOT ESCPLOT6 *DONE. * $,o01: $6*44T *RUN *CCQUEUE PAR:=ESCFLOT6 *EXECUTION BEGINS 2 PLOTSp PLOTTING REQUIRES 94 SEC. AND 24 IN. <$'36 PEN WAS UP 43% OF THE TIME. OK? "1(4L71, I:ESC'LOT6 " HAS BEEN PERMITTED "R F'IEY-== —.CCQUEUE" PLOT ASSIGNED RECEIPT * 718843. *EXECUTION TERMINATED: $ +45, $6.88T C(EkMC OF MFJ6MAL PpOcmbPc;~. 29

7/18/81 co 00 CN C( co in O a H UNDER CONCENTRATED LOADS. C"S -L Q3 o <13 UNOEFDRMEO SHAPE 30

7/18/81t IX-OEF. w.031g 0 00 I CN 0 CN ", 0 \I t0 10 T09 CN 31

W.J.A. 20 July 1981 MSC/NASTRAN PROBLEM SET Sample Problem 3 PLATE UNDER CONCENTRATED LOADS I. Modeling Use six QUAD4 elements to model the plate. This is the simplest element for data input, yet allows contour plotting (not yet available on QUAD8 elements). Fig. 1. F.E. mesh. The plate element QUAD 4 (as well as QUAD8) does not contribute stiffness to nodal rotations about an axis normal to the plate (z direction). Hence the 6th degree of freedom at every node in this problem must be T constrained. / )S The clamp at the wall constrains all degrees of freedom. The piano \ t hinge constrains all degrees of' freedom except rotation along its axis (5th d. o.f.). To visualize deformations of a 1 plate under lateral loads, it is best to take a 3/4 view, such as *- - the default view angle in MSC/ NASTRAN. Your answer for deflections at the corners of the plate should be 10.76 mm (0.424 in). Fig. 2. Degrees of freedom. 32

I:1... I: T E SCD AA rT 12:: 10 IID ANitIE R S.ONY AER0510 $ JULY 22 19.81 v DATA IN 14L7:E SCDATA:1::: 11 TIiE I 12 SOL 2.4 ]::' 1 3::]: DIA 1. 4: 14 CEND:: 00 TITLE=ESC PROBLEM.*3 FINER MESH (6 ELEMENTS).::. 101 SUBTITLE=:- LATE UNDER CONCENTRATED LOADS. 102 ECHO=BOTH 104 DIISPLACEMENT=ALL >:::1.05 STR E S = A L L 0 E::. S lF:R C EA..LL 1:. 0 7 SI"-PCFORCES ALL:::.:1.1:L SUBCASE 1::::. 5 L OAt.1::40 O UT UT( PLOT ) 1::41 PLOTTER NASTRAN: 142 SET INCLUDIE ALL:: 43 F-IND SCALE ORIGIN 1 14S4 FLOT STATIC DEFORMATION 0 SET 1::.53 AXES ZXY: 154 VIEW 0o 0.+0.:: 55 FIND SCALE ORIGIN 1: 156 PILOT LABEL BOTH: 58 FPLOT SHRINK: 159 CONTOUR ZDISP:: 160 F'LOT CONTOUR SET 1 OUTLINE 190 BEGIN BULK:: 1 PARAM AUTOSPC YES 1:: 9 2 GRDSET 6 2:: 01:L GRID 1 0.0 (0+ 0.0 123456 0()2 GRID 2 0.0 250+000 0.0 123456:: 203 GRID 3 0.0 500.000 0.0 123456 204. GRID 4 500O+00 0.0 0.0 12346 2'05 GRID 5 500.000 250.000 0.0 12346:: 0'6'GRID 6 500.000 500.000 0.0 12346 207 GRID 7 750.000 0.0 0.0 2 08 GRID 8 750.000 250.000 0.0 2. "09 GRID 9 750o000 500.000 0.0 210 GRID 10 1000000 (+O 0.0.1. 1 GRID 11 1000.000 250.000 0.0 212 GRIID 12 1000.000 500.000 0.0:: 301 CQUAD4 1 1 1 4 5 2 0:: 02 CQUAD4 2 1 2 5 6 3 3 CQUAD4 3 1 4 7 8 5 3 04 CQUAD4 4 1 5 8 9 6:::.>:305 CQULAD4 5 1 7 10 11 8 30)6 CQ JUAD4 6, 1. 8 11 12 9: 40:L F'S 1.:..I.. l -1 10. 1 501 MAT 16.90E +04 0.30:: 710.1 FORCE:1: 0 0 1.0 0. 0. -1000o:: 712.1 FOIRCE 1,12 0 1.0 0. 0. 1000.:1.000 E N D r ATA'1J:EN OI: F I L. 1: E: 33

*:U..N AER NAST''RAN SCARDS=ESCI'ATA12 SFPRINT:=ESCRES1.2.:IEXEC' UTION B[EG.INS E:XECJCUTIO(:N BEGINS:1:$tICON'TROL —:i: * N 01 SIZE=1000P:$ C 0 NTR 0 L - *:N02 SIZE= 1000P:I$COINTROL —:*N03 S:I:ZE=1000P:I$;CONTROL -*1.N0()4 SIZE=1 000' EXECUT ION DBEG INS NASTRAN L.OADED AT LO CATION 700000 CPFU ELAP I/0 EXCF' COUNTS TRACKS SEC SEC SEC NASTRAN TOTAL USED 0.0 0,0 0.0 0 0 0 SEMI BEGN 0.197 7.0 0 5 11 0 16 XSOR 0 s387 100 1.0 21 0 32 IFP 0,550 12*. 2 0 41 0 52 XGPI:1 2. 85 17+0 3+1 63 0 60 9 GP1 BEGN 1,390 19.0 3*9 78 0 76 11 GP2 BEGN 1,457 20+0 4.1 83 0 76 15 PLTHBDY BEGN 1.496 220 4*3 85 0 76 17 PLTSET BEGN 1+585 26*0 4.9 98 0 92 19 PRTMSG BEGN 1.620 27*0 5*0 100 0 92 22 PLOT BEGN 2.688 32'0 5.8 116 0 92 23 PRTMSG BEGN 2.763 33*0 5.9 118 0 92 26 GP3 BEGN 2*853 35.0 6*5 130 0 92 28 TA1 BEGN 3+092 40,0 7,5 151 0 92 31 EMG BEGN 3.277 42.0 8.0 160 0 92 37 EMA BEGN 3,760 47*0 9+0 180 0 92 97 GP4 BEGN 3 863 49 0 9.7 194 0 96 99 GPSP1 BEGN 3.956 51.0 10.1 203 0 96 110 SCE1 BEGN 4,121 54.0 10,8 216 0 108 140 RBMG2 BEGN 4,283 56.0 11*1 222 0 108 145 SSG1 BEGN 4.4432 60.0 1. *8 235 0 108 148 SSG2 BEGN 4.585 62.0 12.5 250 0 120 155 SSG3 BEGN 4,780 65,0 13,2 264 0 120 160 SDR1 BEGN 5s088 68,0 14+6 293 0 128 176 SDR2 BEGN 5+370 72.0 15.8 316 0 128 185 OFF BEGN 5 453 74 0 15*9 319 0 28 128 86 SDRX BEGN 5+513 75*0 16.3.325 0 128 188 OFF BEGN 5,675 80 0 16.4 329 0:128 189 GPFDFR BEGN 5*793 81.0 16.6 333 0 128 190 OFF' BEGN 5,828 82.0 16.7 334 0 128 193 OFF' BEGN 5 976 86 0 16+9 339 0 128 205 PLTSET BEGN 6,072 88*0 17 6 352 0 128 206 FRTMSG BEGN 6+096 89.0 17+7 354 0 128 207 PLOT BEGN 7 738 99.0 19.2 384 0 128 208 PRTMSG BEGN 7 812 100*0 19. 3 386 0 - 28 2 14 EX IT BE-GN 11-0002I STOP:1 VSSCMD238W NSTN04 - OPEN DSDEF BEING CLOSED V3.SCMD238W NSTNO3.- OPEN DSDEF BEING CLOSED VSSCMlD238W NSTNO2 - OPFEN DSDEX:1F BEING CLOSED VSYSCMD238W NSTN01 - Ol::'EN DS.DEF BEING CLOSED F1' I F:'PlO.T'DESCIRI:P:i:: TION GI:ENEFRATION BEG INS.:t1 E X E C U' I0 iN T l-: RH:1: N A T li D 34

:1::'CON TR 0: L * PR N'T I.I 0 L D:ii:.'!F'RINTT< AS.S:TIGNED RECEIPT NUMBER 601796:IL ESCD.OAl A 2 *PF'R I N T* U:IC ESCRES:L2 *1F'PRINT*::CONTROIL *F:'RINTI* RELEASE':F:'PIRINT.' 60:1796 RELEASED TO CNTRY 22 PAGES.:I:RENANE -FPLOT ESCPLOT12.: DONE:$RtUN *^CCQUEUE PAR=ESCPLOT12:EXECUTION BEGINS 4 FPL'OTSY PLOTTING REQUIRES 235 SEC. AND 47 IN; $. 79 P::'N WAS UP 41% OF THE TIME. OK? OK "K4L7:ESCPLOT12 HAS BEEN FERMITTED "R PKEY=*CCQUEUE", PLOT ASSIGNED RECEIPT. 723210,:I:X ECUT ION TERMINATED 35

1 7/aim 0 r) W Wl W (0 0 CN I1 -J 4

2 7/23/1 2 -. N ~~~~~~~~~37 N 3 c/

3 7/23/81 7 / X-F. - 1'76101W 3 0 N CD I -J z- - 0,.. D N STATIC DEFM. SmiCIZ 1 L 1 38

4 7/25/B1 ceF a -O.*7S164 00 I C) N KJ) (b CT N I&-,a ae..M -W AM-b m A39 <s ____ CN' ---- r^ qE5______^rn~s ~cao

40 MSC/NASTRAN PROBLEM SET Sample Problem 4 SPUR GEAR TOOTH I. Modeling Use two HEXA solid elements with 20 nodes each as shown in Fig. 3 1. Align the element boundary at/ the pitch circle where the running \ load is applied. This means the equivalent nodal loads for case (A) will be applied to nodes 15, 16, 17 with values of 1000, 4000 and 1000N -. - -_ _7 respectively. (These are the values |' [! consistent with a quadratic inter- - 1 polating function. See Zienkiewicz, Ed. 3, pp. 222-235). For case B, the full 6000N (1, 349 lb.) acts on 9 node 15. This will be our first example with multiple load cases. 7 6 The nodal coordinates, as \ determined from Fig. 2. of the problem statement, are: I Fig. 1. Two-element model of spur gear tooth using HEXA elements. 1 386.500 12.700 5.000 2 86.500 0.0 5.000 3 86.500 -12.700 5.000 4 86.500 -12.700 0.0 i5 86.500 -12.700 -5.000 6 86.500 0.0 -5.000 7 86.500 12.700 -5.000 8 86.500 12.700 0.0 9 93 000 8.700 5.000 10 93.000 -8 700 5.000 11 93.000 -8.700 -5.000 12 93.000 8.700 -5.000 13 99.600 7.800 5.000 14 100.000 0.0 5.000:1 5 99. 600 -7.800 5.000 16 99.600 -7.800 0. 0 17 99.600 -7.800 -5.000 18 100.000 0.0 -5.000

19 99.600 7.800 -5.000 20 99.600 7.800 0.0 21 105000 700 5700 5.000 22-. 105.000 -5.700 5.000 23 105.000 -5.700 -5000 24 105.000 5 700 -5.000 25 110.000 3.500 5.000 26 110 000 0.0 5.000 27 110.000 -3.500 5.000 28 110.000 -3.500 0.0 29 110.000 -3.500 -5.000 30 110.000 0.0 -5.000 31 11(0000 3.500 -5,000 32 110.000 3.500 0,0 The HEXA element connectivities, as seen from page 2. 4-20? in the User' s Manual, are Element 1: 3, 5, 7, 1, 15, 17, 19, 13, 4, 6,.8, 10, 11, 12,9, 16, 18, 20, 14 bottom, far side, bottom, middle, far side, corners corners midside midside midside Element 2: 15, 17, 19, 13, 27, 29. 31, 25, 16, 18, 20, 14, 22, 23, 24, 21, 28, 30, 32, 26 Plots of the gear tooth are best obtained by using: AXES Z, X, Y VIEW -20., -30.,80. This shows the tooth in the upright configuration as in the problem statement. 41

, -. pQ - -.O -C *C 3 _-..>0 CO':<" (<^?<*?^) Ce N CO C ~?.4'" —,'...r"4'. _ t'' -' C;,. J........?. i~ ~ ~~~:- - ":Z.. ~~ - Z;Q we LL C C. Z - sr~~~~~ ~~"'" 000 0 00 0..... S. -'. 0Ls - I o u - c 4 f n _ I.-. i i C^ i; *,'.. *i7'.,0n "' 2 -O."n 0.0..' 0 O " 0'-..-. 3.5- - - < o *> ^ 4 ^ 4 >. 4 4 v >:- ZO C 0 w m 0 CO 0 03 ^ > O 0 _~. f *> Ys,-..s _' -,, <, >- - _ Q # _. c:s; in *: _ o => _ f 03. ^ ]]..' *r ~ ~ r ~."- C'~. V...'., — _ i ~ S'.'-~- 0 2: _ ~ X._.-'?' l,>k. *3 c; ~J r >- z?~, * 3 *X:: - CE - i! 9- 1>~; L S _ s..2..' >-?....S'... L: L. i. _. s _s _ * s __3. XO; _ <2_ _ _ _ _ _ _3. _ _ _._ *.. —... ~ " - "^.1 Z:m -". O ~'.-4.' C' ~ s 3 -.. ^ 0- -".-. * " ^.,.; -.'.."3 42

004 —4-.00-4-. -s Wr: N1 NO. -4.-4 -4 -i i-.4.. -.00 00 W -' OO O O 3 NN 0 Sf 0-l 0 C f 4 O 4+ 4 4> 4 4 4 04 4 4 ~ ~....44 4 4+ 4 0 -000 000000000 00 _ O 5 O 000-0:L -s O 0 0 0 - - 0- II j - 40 0.4 OO-OO 0000 1- i04 0 0 0-4 0. - -. i'.-' CO en O e O v en OO Oen en r -.n o.r4 4. o 4 r\ r\ 4.) o;? in i) ^ o o^ i3 ~i. o 4f, - o ooo _t. —'o oc. oC -oo ooooo o' oo y — 3.R v0 *- *2*3 0 v 0 0v 4...> 0. ^ "' Lv u w 0 0 0 0 0 0 0. ~. J4 [! *! i- i;,I:"U *3 ^ *3 s m.v ~.w s x A t a A A 3 ^ f | ^ _ _ _ _ t._'.0'~"'" 3 ~0'-00 h 0 0.0 o O0 0 """'"'-" 0........ 4 _.'C h - 3 3h 3 3 A A O 4 * 4 4 O., -,'4 r' "',.. -C *. -,.*. C 6'.,.*.,.-.. C. 4 ~ 4 *. C". 4 C-.i?..' ~4, *0...4 0,. 0 0 0 } 0 04 0 * < 5 T 0 ~,,,,~. ~.-~ ~-',~.~.~ v.2.. ~- i v, ~ ~ ~ ~4 ~4 ~.I.-4.-i::,,~ -— " -- ~ C.. ~- ~4 "'...'~4.4~. ~'~'. ~ ~0 O 00 00. ~.'.-.4~....~~~i ".i 4~ ~,4,...'4'.24.?' ~.~.T' i'. (-'4.~..?.4 4 ~.4.~.. ~. ~ ~ = ri h: ii,q'''.f.%.'..'....'....,,...'...'....-.......,.....'...'-.,'....-...;...''....-,,'....,~..'...'...,'-..,....'....-....'....'..........'..'". 43~~~~~~~~~~~~~~~i~

2: zz z t3t3 2:z3:;:32:3z:z2zzzz:3 z:2 zzzz;:2z22zz 2Z * s i; * t * t s.s _'-,,,.. s._ s _ t- >: r..D rD r ~.t.D:,_q~. E CD. CD CZ..D D CDDrD.D..:.CO Z Z.. r LJ v' J:.I..LL L.'.'. >L 0o 0' *4._ e ZQ - "'4- sD::C 3 1...~, t.t r_ Co, I - 0 3' * t 2 t - 4i * - t _' " s s * n * 5_,. o_ _.__. _ _.' - MLJ X IC-1C_,Cs k^t.. -.^'. a -., t N- _ 0,_..- S. -. C3 _-'.' ^ _ 0 0 3 0 -O - _ - 0 0. _0. _L CO. _,D T' n. s -.t 03', >'. r *- 3. W m C 4''_ 0 4 O. X N 0 -- ^ C.- 4'. -4'< N.. 4 N.',- 4 - r 0 0 N,- 4.-. i.-.-.... C.4.-z:'',, 0 ('.4 ~.- 40,,, ~4 C'4;-' * s4 C'4 0'*4! <,s3 O ~.4'.. " ".. - -4 o~ ~4S 0 0 -<? ^C'.< C'4 ~.4.-4 ~' c4,.~i.. - -' "~0"" ^ "" "'" "-I O 0 00 r 0 O00C, 000"'"'" 0 "-" O " " - " -' ~. -. Z~ ~~-', t < S o-.' - ~4 W.. _ 03 03 0 C'.- 0. -- 0? 4 S. C C <r o a. ^ ~~- <'....-'. —4*i %-.4.-,. i0,'.'c4 4. U' 0 V I' I. r''' %4 ~:z * v __*N Ni<_ N- w~"s ^J- ~ - - 9. ~ ~ * N M.O -0 - - -.'-":...... ~ 0 0 0,-. 0 Ld S'4.C-4, - C *. -'C ~3'3 k3 b3 b3. 0 0' 0 0. —..'" C~ ~ "~.o.. —~;Z.- ~4: o o o t - 3 *- _. —' - c-4..' < O n. -,o.-3 P<3 _.,^ 4 -- * — t 0 L, O:'-'._ 2: 2 C....J s. 3 X X:',.. s_ s _ X < 44

z i z Ciz zzz. T-i:3 ~f.._ 0 Cn: _. 0:, C. 0 \!'- ~7 _j X_ J s. ~00; U U i:i " CL. _. O Z U r-sc. Z'J,,, s Z -C;-~ US;Lt i.Li *, z <. ZZ Z -*. 0 _ o en en en I " j_ c o Q ci o <C E4 - s _'''-'. *- l - ",$ * s 7 _0' -J,, - <C U,, L t D CD ID CD LZ l 4 o0!rf),4,S 0 0 * 0 22: r. i- = I=: ccX: r-' L-W- 2L. i- I 0 J0 0 1?.-<:. II l l C U J9'9! V i-Hi —H C-Ci Z.C.nia.00:d Z' Z n Z:n Z.- X' L' O r.z z: t:e. *- 2s: s { *< iL-I - LL00ZL C.LL.i ro r^ s o ao C3; D * o c en b 34 _t 0:. U( U U * M 2 Q f X i - U ooo l yrf i X-L. ~. =DX bJ " C:. X'1 Uf 2 2 2 > 2 i^^O z u:t3 OCa*S

1 7/23/81 oo =0 1 1 0 rL/ I \ 31~g17 E~LIbI4 R t'OTH UNER LINE MD CONCENTRATED LoM. UNDEFOR1ED FPE f1)

-2 7/23/81 MAX-OEF. 0.01465214 2 00 r) 0i o TLT DF.mRt TH UNCSER LINE ANM CONCENTRATED LDS. 1)- STATIC EFOR. SMCASE 1 LOAD 1 a)

3 7/23/81 MAX-DEF. - 0.02652046 0 IC) LO Oi 4-, EOB% 4taH UNDER LINE AND CONCENTRATED LOADS. STATIC DEFOR. SUBCASE 2 LOAD 2

Helen Buus William J. Anderson Edward Chmiel February 22, 1982 PRENASTRAN - A PREPROCESSOR FOR PREPARING~ MSC/NASTRAN DATA INTRODUCTION The AERO:PRENASTRAN program allows the user to create a complete NASTRAN data file interactively which is suitable for input to AERO:NASTRAN or which can be punched on cards for execution on a system other than MTS. Such a data file can be generated for one of four types of problems: 1) Plane Stress (using CQDMEM or CQUAD8 elements) 2) Beam (using CBAR elements) 3) Plate (using CQUAD4 elements) 4) Solid (using CHEXA elements) The program can be run three ways: 1) Create - to generate new NASTRAN file 2) Restart - in case of interruption during Create mode, this capability allows the user to continue generating the file at the point of interruption 3) Edit - to make minor modifications to an existing complete NASTRAN file. While running the program, read all messages and prompts carefully before entering any information. The program accepts free format input of numerical data; simply separate items with commas or semicolons. Certain steps of file generation will continue until the user gives an "end-of-file" command. This can be done by "$ENDFILE." The first time the program is run, a file named "DATA" will be created for the user and NASTRAN data will be put into it. If a new NASTRAN file is to be created later and this file still exists, the user can request to: 1) Enter a new file name to be created; 2) Empty the existing file and use it; 3) Use the Restart or Edit options on the existing file. The program will prompt the user for information as required during file generation. There are eight steps in generating a file: 1) Executive Control Deck - problem identification 2) Case Control Deck - title, subtitle, and number of subcases 3) Grid Point Definitions - grid point id's and coordinates 4) Grid Point Constraints -,grid point id's and constraint codes 5) Element Definitions - element id's, property id's and grid point connectivities 6) Property Definitions - material id's and thicknesses or section properties, if required 7) Material Definitions - Young's Modulus and Poisson's Ratio 8) Load Definitions - grid point id's and Cartesian components of forces.

2 This manual is divided into sections corresponding to each step of file generation. Each section describes what the program does and what information the user is expected to provide. EXECUTIVE CONTROL DECK The user need only provide two 8-character names, typically last name and problem name, separated by a comma. Example: Executive Control Deck Enter Problem Identification (2.2-1):..ID? ANDERSON,STRESSES (carriage return) Note the reference to the NASTRAN User's Manual, 2.2-1. This is the page number where information about the ID card can be found. The program provides thisinformation throughout to direct the user to the appropriate page in the manual,should a question ever arise on a given type of NASTRAN card. CASE CONTROL DECK Specify the problem type. Example: Enter problem type (l=Plane Stress, 2=Beam, 3=Plate, 4=Solid):? 1 For the plane stress cases only, the user must then choose between two elements, i.e., CQDMEM or CQUAD8. Example: Choose element type for Plane Stress: 1 - CQDMEM (for stress contour plots) 2 - CQUAD8 (for better accuracy with fewer elements)? 1. Once the problem type has been determined, the program sets up the appropriate output and plot requests automatically. The user must now provide a title and subtitle and specify how many different loading conditions are to be applied to the structure. Examole: Case Control Deck Enter Title (2.3-137): TITLE=? PLATE ON PIANO HINGE SUPPORT Enter Subtitle (2.3-117): SUBTITLE=? FEB 4, 1982 From one to nine subcases may be requested to allow application of different sets of loads to the model. For each subcase, the user will be prompted for a set of loads in Step 8 of file generation. Enter desired number of subcases (2.3-55 and 2.3-111):? 1

3 GRID POINT DEFINITIONS Grid points are used to define finite elements. Each grid point is identified by an integer number from 1 (one) to a user defined maximum which must be less than or equal to 99. No gaps in numbering are permitted, although grid points may be entered in any order. These points are defined by their X, Y and Z coordinates. If the problem type is plane stress, the user need only define X and Y coordinates. Example: GRID Grid Point Definitions (2.4-205) Define each grid point by entering its id number and coordinates. Repeat process for each grid point. If an error was made in defining a grid point, the same id and correct coordinates may be re-entered, thereby replacing the old information in the data file. When all grid points have been defined, enter $ENDFILE or CTRL-C to stop. The program will then check if all grid points numbered from one to the user defined maximum have been defined. Enter maximum grid id:? 6 Maximum grid id = 6, please confirm (Y or N):? Y For each grid point, enter ID,X,Y,Z:? 1,0,0,0? 2,0,10,0? 3,-8,0,0? 4,-8,10,0? 5,-16,0,0? 6,-16,0,0? $ENDFILE GRID POINT CONSTRAINTS Grid point constraints are used to specify boundary conditions for a problem. They tell NASTRAN which degrees of freedom (d.o.f.) are to be constrained Or fixed for each grid point. Any of the digits 1 through 6 may be specified corresponding to the degrees of freedom one wishes to constrain. The digits correspond to the degrees of freedom as follows: 1) Constrain translation in the X direction 2) Constrain translation in the Y direction 3) Constrain translation in the Z direction 4) Constrain rotation about the X axis 5) Constrain rotation about the Y axis 6) Constrain rotation about the Z axis. Example: Define the constraints for a grid point (GRID 2,3-205): Enter grid number, constraints (any digits 1-6, no blanks). Enter $ENDFILE or CRTL-C to stop.? 3,12346? 4,12346? 5,123456? 6,123456? $ENDFILE

4 ELEMENT DEFINITIONS Finite elements are defined by connecting them to previously defined grid points. Each element is identified by an integer number from 1 (one) to a user defined maximum (< 99) as was the case with grid points. An element is defined by specifying the grid points to which it is connected in a well defined order, i.e. element connectivity. The user should read the appropriate page in the NASTRAN User's Manual for the element being used: CQDMEM - 2.4-103 CQUAD8 - 2.4-112c CBAR - 2.4-39 CQUAD4 - 2.4-112a CHEXA - 2.4-69 The program expects the user to provide it with an element identification number or EID, for each element, followed by the grid point id's to which it is connected. If the user has specified problem type 4, he/she may choose to specify only the eight corner grid points of the solid CHEXA element or all twenty grid points for better accuracy. For all problem types, the user may choose to use the same property id (and orientation vector for CBAR elements) for all elements or he/she may choose to specify different property id's for each element. The property id's identified at this point in the program will be used later to define the appropriate property cards. Example: CQUAD4 Element Definitions (2.4-112a) Define each element by entering its id number and grid id's specifying the element's connectivity. If an error was made in defining an element, the same id and correct connectivity may be re-entered, thereby replacing the old information in the data file. When all elements have been defined, enter $ENDFILE or CTRL-C to stop. The program will then check if all elements numbered from one to the user defined maximum have been defined. Enter maximum element id:? 2 Maximum element id = 2, please confirm (Y or N):? Y For each element, enter EID,GID1,GID2,GID3,GID4 with grid id's in counter-clockwise order: NOTE - Following definition of the first element, there will be a prompt for property id number and an option will be given to either use the same property id for all elements or to enter a new id for each element Enter element data:? 1,1,3,4,2 Enter property id:? 1 Is given property id to be used for all elements (Y or N)?? Y Enter element data:? 2,3,4,6,4? $ENDFILE

5 PROPERTY DEFINITIONS Property cards are used to define thicknesses of membrane and shell elements such as CQDMEM, CQUAD8 and CQUAD4 elements, or to define section properties of CBAR elements. All property cards require material id's for use on material cards defined later. The program keeps track of property id's specified while defining elements and prompts the user for the data required to prepare each property card. This data depend on the problem type: 1 - PQDMEM (2.4-307) or PSHELL (2.4-326a) Material id and thickness 2 - PBAR (2.4-273) Material id, cross section area, principle moments of inertia, II and 12, and dimensions of section (Dimensions are used to calculate stress points in the section.) 3 - PSHELL (2.4-326a) Material id and thickness 4 - PSOLID (2.4-326c) Material id. All material id's must be in the range 1 to 99. An example for problem type 3 is given below: PSHELL Property Definitions (2.4-326a) Enter material id and thickness for property id 1:? 1,0.25 If the program is in Edit mode, the user will first be prompted for problem type. He/she may then enter the property id to be added or modified by the program. Property cards may be deleted by entering the negative of the unwanted property id. If a positive number is entered, the user will be prompted for the necessary property data described above. MATERIAL DEFINITIONS Material cards are used to specify the Young's modulus and Poisson's ratio associated with the material in the finite element, e.g., steel, aluminum. The program keeps track of material id's specified during property definitions and prompts the user for Young's modulus and Poisson's ratio. Example: MAT1 Material Definitions (2.4-215) Enter Young's Modulus, Poisson's Ratio for material 1;? 10.E6,0.3 In Edit mode, the user must first specify the material id to be added or modified. Or the nega-ive of an id may be entered to delete an unwanted card.

6 FORCE DEFINITIONS The last step in defining a NASTRAN data file is giving,the forces on the structure. For each subcase requested in the Case Control Deck, the user must define a set of loads, Each load is defined by specifying the grid point to be loaded and the X,Y,Z components of the force to be applied at the point. Example: FORCE Load Definitions (2.4-179) 1 Load subcase(s) have been requested Define loads for subcase 1 For each load, enter grid id defining point at which load is to be applied, followed by X,Y,Z components of force, i.e., GID,FX,FY,FZ. To delete a load, enter -GID,0,0,0. Enter $ENDFILE or CRTL-C to stop for given subcase;? 1,0,0,100? 2,0,0,100? $ENDFILE

s5SE N S MANJUAL.,ER. IO. QA

PREFACE MSC/NASTRAN is a large scale, general purpose digital computer program which solves a wide variety of engineering problems by the finite element method. MSC/NASTRAN, a version of the NASTRAN general purpose structural analysis program, has been developed and Is maintained by The MacNeal-Schwendler Corporation (MSC). NASTRAN is a registered trademark of the National Aeronautics and Space Administration (NASA). MSC/NASTRAN is marketed and serviced from MSC's offices in Los Angeles, Munich, and Tokyo, and is available at most major public data centers. The program is operational on most major computing systems, including IBM 360/370 series, IBM 3000 series, IBM 4300 series, Amdahl, Itel, Fujitsu M series, CDC 7600, CDC CYBER series, Univac 1100 series, Digital's VAX 11/780, and the CRAY. The User's Manual is restricted to those items relating to the use of MSC/NASTRAN that are independent of the computing system being used. Computer-dependent matters such as operating system control cards are treated in Section 7 of the MSC/NASTRAN Application Manual. The Application Manual is updated monthly and also contains Newsletters, the Current Error List, Documentation for the rigid format alter library, and other current information which supplements the MSC/NASTRAN User's Manual. In addition to the MSC/NASTRAN User's Manual, the following manuals are available for use with MSC/NASTRAN: 1. MSC/NASTRAN Primer: Static and Normal Modes Analysis 2. NASTRAN Theoretical Manual 3. MSC/NASTRAN Handbook for Linear Static Analysis 4. MSC/NASTRAN Application Manual 5. MSGMESH Analyst's Guide 6. MSC/NASTRAN Aeroelastic Supplement 7. MSC/NASTRAN Demonstration Problem Manual 8. MSC/NASTRAN Programmer's Manual The MSC/NASTRAN Primer serves as a text on the basic concepts of finite element analysis and as an introduction to the use of MSC/NASTRAN for static analysis and real eigenvalue analysis. The Theoretical Manual not only presents the analytical and numerical procedures that underlie the program, but it also describes the structure and problem-solving capabilities of NASTRAN in a narrative style. The NASTRAN Theoretical Manual presents the analytical and numerical procedures used in the program and also describes the structure of the problem solving capabilities of NASTRAN in a narrative style. i (2/2/81)

The MSC/NASTRAN Handbook for Linear Static Analysis is an introductory manual for beginning users of MSC/NASTRAN and is also a specialty manual for the most frequently used aspects of linear static analysis. Approximately to-thirds of the Handbook is devoted to detailed descriptions and classifications of input data cards with emphasis on ease of data retrieval. The remaining onethird of the Handbook is an original, easily read textual description of the capabilities of MSC/NASTRAN for linear static analysis and of the general procedures to be followed in preparing input data. The MSC/NASTRAN Application Manual contains specific information ration eltive to the following topics: 1. Estimation of resources required to execute a problem. 2. Procedures for problem execution on the analyst's computer. 3. Articles on special techniques and methods of analysis. 4. Articles on MSC/NASTRAN features. For example, there are several articles on the elements that are recommended for general use. These discussions, which include hints on modeling procedures and examples of specific applications, expand on the information available in the User's Manual. 5. The contents of the RFALTER library are described. This library, which is delivered with each new version of MSC/NASTRAN, contains DMAP Alters for the various solution sequences. These alters provide additional capability, improved efficiency and/or otherwise unavailable user convenience. 6. MSC/NASTRAN Application Notes and monthly Newsletters. The MSC/NASTRAN Application Manual is updated monthly with Newsletters, Application Notes, articles, and additions to the other MSC/NASTRAN Manuals. The MSGMESH Analyst's Guide describes the automatic data generation capabilities of MSC/NASTRAN. The MSC/NASTRAN Aeroelastic Supplement describes the theoretical aspects and the numerical techniques used to perform an aeroelastic analysis with MSC/NASTRAN. The MSC/NASTRAN Demonstration Problem Manual contains a group of small problems that illustrate the various analysis methods available in MSC/NASTRAN. The data decks for these demonstration problems are delivered as a separate file with each new release of the program. The MSC/NASTRAN Programmer's Manual contains a description of the program code, including the mathematical equations which are implemented in the functional modules. The MSC/NASTRAN Aeroelastic Supplement describes the theoretical aspects and the numerical techniques used to perform an aeroelastic analysis with MSC/NASTRAN. The Supplement is not intended to be a self-contained text since it relies extensively upon information relative to the dynamic analysis of a finite element model that is provided in the MSC/NASTRAN User's Manual and the NASTRAN Theoretical Manual. ii (2/2/81)

The MSC/NASTRAN Demonstration Problem Manual is a useful reference in the preparation of the NASTRAN Data Deck. This manual is a suite of relatively small problems that illustrate the various analysis methods available in MSC/NASTRAN. The data decks for these demonstration problems are delivered as a separate file with each new release of the program. The MSC/NASTRAN Programmer's Manual contains a description of the program code including the mathematical equations that are implemented in the functional modules. It also contains information relative to the maintenance and modification of the program. MSC/NASTRAN uses a finite element structural model, wherein the distributed physical properties of a structure are represented by a finite number of structure elements which are interconnected at a finite number of grid points to which loads are applied and for which displacements are calculated. Procedures for defining and loading a structural model are described in Section 1. The modeling procedures for heat transfer analysis, hydroelastic analysis, and acoustic cavity analysis are also described in Section 1. The Data Deck, including the details for each of the data cards is described in Section 2. This section also discusses the control cards that are associated with the use of the program. Section 3 contains a general description of solution procedures provided in the program, along with specific instructions for the use of each type of solution sequence. Each of the solution sequences is associated with the solution of problems for a particular type of static analysis, dynamic analysis, or heat transfer analysis. The procedures for using the plotting capability are described in Section 4. Both deformed and undeformed plots in the structural model, as well as contour plots of output quantities, are available. Response curves are also available for transient response and frequency response analyses. In addition to the solution sequences provided by the program, the user may choose to write his own direct matrix abstraction program (DMAP). This procedure permits the user to execute a series of matrix operations of his choice, along with any utility models or executive operations that he may need. The rules governing the creation of DMAP programs are described in Section 5. The diagnostic messages are documented and explained in Section 6. MSC/NASTRAN is intended for use by professional engineers to assist them in the solution of engineering problems. It is expected that MSC/NASTRAN will be used in combination with other analytical and experimental procedures in the practice of engineering analysis and design. A new edition of the MSC/NASTRAN User's Manual is printed for each new version of MSC/NASTRAN. New versions will include corrections to some of the errors on previous error lists. Others will be further noted on the Current Error List or moved to the General Limitations List, both of which are located in Section 6 of the MSC/NASTRAN Application Manual. iii (2/2/81)

A great many people have been associated with the development of MSC/NASTRAN and all have had some influence on the preparation. of the MSC/NASTRAN User's Manual. Much of the material for the original edition of this manual was taken from the NASTRAN User's Manual published by the National Aeronautics and Space Administration as NASA SP-222(01). Most of the members of both the programming staff and the engineering staff have made some direct contribution to the User's Manual. Particular recognition is due Mr. Steven E. Wall who rewrote Section 5 on Direct Matrix Abstraction Programs. iv (2/2/81)

TABLE OF CONTENTS Section Page No. Preface.............................................................. Table of Contents..................... v Page Status Log........................................................... xii 1. STRUCTURAL MODELING 1.1 INTRODUCTION..................................................11.2 GRID POINTS......................................... 1.2-1 1.2.1 Grid Point Definition.................................... 1.2-1 1.2.2 Grid Point Sequencing..................................... 1.2-3 1.2.3 Grid Point Properties..................................... 1.2-8 1.3 STRUCTURAL ELEMENTS..................... 1.3-1 1.3.1 Element Definition...................................... 1.3-1 1.3.2 Beam and Bar and Bend Elements.......................... 1.3-3 1.3.3 Rod Element............................................... 1.3-5 1.3.4 Shear Panel.............................................. 1.3-6 1.3.5 Shell Elements............................................ 1.3-6 1.3.6 Conical Shell Element................................. 1.3-9 1.3.7 Axisymmetric Shell Elements.............................. 1.3-12 1.3.8 Scalar Elements................................... 1.3-13 1.3.9 Solid Polyhedron Elements................................. 1.3-15 1.4 CONSTRAINTS AND PARTITIONING.................................... 1.4-1 1.4.1 Set Definition................................... 1.4-1 1.4.2 Multipoint Constraints.................................... 1.4-4 1.4.3 Single-Point Constraints................................. 1.4-6 1.4.4 Static Partitioning and Reduction......................... 1.4-7 1.4.5 Free Body Supports........................................ 1.4-8 1.4.6 Generalized Dynamic Reduction............................. 1.4-9 v (2/13/81)

TABLE OF CONTENTS (Continued) Section Page No. 1.4.7 Sets for Dynamics......................................... 1.4-10 1.4.8 Sets for Aerodynamics.................................... 1.4-10 1.4.9 Conventions for Undefined Degrees of Freedom.............. 1.4-11 1.4.10 Output Selection Via Set Specification.................... 1.4-12 1.5 APPLIED LOADS......................................... 1.5-1 1.5.1 Static Loads........................................ 1.5-1 1.5.2 Frequency-Dependent Loads................................. 1.5-5 1.5.3 Time-Dependent Loads..................................... 1.5-5 1.5.4 Superelement Load Partitioning.................. 1.5-8 1.6 DYNAMIC MATRICES........................................ 1.6-1 1.6.1 Direct Formulation....................................... 1.6-2 1.6.2 Modal Formulation....................................... 1.6-3 1.7 HYDROELASTIC MODELING..................................... 1.7-1 1.7.1 Solution of the NASTRAN Fluid Model...................... 1.7-1 1.7.2 Hydroelastic Input Data.................................. 1.7-3 1.7.3 Rigid Formats............................................. 1.7-6 1.7.4 Hydroelastic Data Processing.............................. 1.7-8 1.7.5 Sample Hydroelastic Model................................. 1.7-9 1.7.6 The Virtual Mass Method.................................. 1.7-10 1.7.7 Example of the Incompressible Virtual Mass Method......... 1.7-12 1.8 HEAT TRANSFER PROBLEMS.............................................. 1.8-1 1.8.1 Introduction to MSC/NASTRAN Heat Transfer................. 1.8-1 1.8.2 Heat Transfer Elements................................... 1.8-2 1.8.3 Constraints and Partitioning.............................. 1.8-3 1.8.4 Thermal Loads............................................. 1.8-4 1.8.5 Linear Static Analysis................................... 1.8-4 1.8.6 Nonlinear Static Analysis................................. 1.8-5 1.8.7 Transient Analysis........................................ 1.8-6 1.8.8 Compatibility with Structural Analysis..................... 1.8-7 vi (2/13/81)

TABLE OF CONTENTS (Continued) Section Page No. 1.9 ACOUSTIC CAVITY MODELING......................................... 1.9-1 1.9.1 Data Card Functions....................................... 1.9-1 1.9.2 Assumptions and Limitations.............................. 1.9-2 1.9.3 Acoustic Cavity Example Problem........................... 1.9-3 1.10 (Deleted) 1.11 CYCLIC SYMMETRY..................................................... 1.11-1 1.11.1 Cyclic Symmetry Solution Sequences....................... 1.11-1 1.11.2 Solution of Axisymmetric Problems........................ 1.11-4 1.12 SUPERELEMENT ANALYSIS............................................... 1.12-1 1.13 NONLINEAR STRUCTURAL ANALYSIS...................................... 1.13-1 1.13.1 Geometric Nonlinearity (SOL 64)......................... 1.13-2 1.13.2 Nonlinear Static Analysis (S0L 66)....................... 1.13-3 1.14 LAYERED COMPOSITE MATERIAL ANALYSIS................................. 1.14-1 2. NASTRAN DATA DECK -" 2.1 GENERAL DESCRIPTION OF DATA DECK.................................... 2.1-1 - 2.2 EXECUTIVE CONTROL DECK.............................................. 2.2-1 2.2.1 Executive Control Card Descriptions....................... 2.2-1 2.2.2 Executive Control Deck Examples........................... 2.2-7 - b 2.3 CASE CONTROL DECK................................................... 2.3-1 2.3.1 Subcase Definition........................................ 2.3-1 2.3.2 Case Control Card Summary................................. 2.3-3c 2.3.3 Case Control Card Descriptions........................... 2.3-8b 2.4 BULK DATA DECK......................................... 2.4-1 2.4.1 Format of Bulk Data Cards................................. 2.4-1 2.4.2 Bulk Data Card Summary.................................... 2.4-4 2.5 (Deleted) vii (2/13/81)

TABLE OF CONTENTS (Continued) Section Page No. 3. SOLUTION SEQUENCES 3.1 GENERAL DESCRIPTION OF SOLUTION SEQUENCES........................... 3.1-1 3.1.1 Restart Procedures........................................ 3.1-4 3.1.2 Solution Sequence Output.................................. 3.1-11 3.1.3 Parameter Descriptions................................. 3.1-16 3.1.4 Use of Parameters........................................ 3.1-50 3.2 CASE CONTROL DECK...................................... 3.2-1 3.2.1 Case Control Structure............................. 3.2-1 3.2.2 Output Options.................................... 3.2-6 3.2.3 Superelement Case Control Structure.................... 3.2-11 3.2.4 Grid Point Stress Requests.......................... 3.2-13 3.3 SOLUTION SEQUENCE OPERATIONS................................... 3.3-1 3.3.1 Geometry Processors...........................3.3-1 3.3.2 Matrix Assembly Operations.................... 3.3-2 3.3.3 Multipoint Constraint Operations................. 3.3-3 3.3.4 Single-Point Constraint Operations....................... 3.3-4 3.3.5 Partitioning Operations.................... 3.3-5 3.3.6 Rigid Body Operations................................... 3.3-6 3.3.7 Static Solution................................. 3.3-8 3.3.8 Data Recovery Operations.............................. 3.3-11 3.3.9 Real Eigenvalue Analysis................................ 3.3-16 3.3.10 (Deleted) 3.3.11 (Deleted) 3.3.12 Formulation of Dynamics Equations......................... 3.3-19 3.3.13 Complex Eigenvalue Analysis................ 3.3-22 3.3.14 Frequency Response and Random Analysis................... 3.3-22 3.3.15 Transient Response...................................... 3.3-23 3.3.16 Linear Heat Transfer Analysis............................. 3.3-24 3.3.17 Nonlinear Heat Transfer Analysis....................... 3.3-24 viii (10/30/81)

TABLE OF CONTENTS (Continued) Section Page No. 3.3.18 Transient Heat Transfer Analysis.......................... 3.3-24 3.3.19 Plotting Operations....................................... 3.3-25 3.3.20 Cyclic Symmetry....................................... 3.3-25 3.3.21 Superelement Processing....................... 3.3-27 3.3.22 Superelement Model Checkout.............................. 3.3-33 3.3.23 Superelement Alternate Statics Solution................... 3.3-35 3.3.24 Dynamic Reduction and Component Mode Synthesis............ 3.3-36 3.3.25 Use of Superelement Solution Sequences Without Superelements......................................... 3.3-41 3.4 COMPUTER RESOURCE REQUIREMENTS...................................... 3.4-1 3.4.1 Main Storage Requirements................................. 3.4-1 3.4.2 Secondary Storage Requirements........................ 3.4-2 3.4.3 CPU Time Estimation....................................... 3.4-3 4. PLOTTING -> 4.1 PLOTTING........................................ 4.1-1 4.1.1 General Capability...................................... 4.1-1 4.1.2 Superelement Plotting.................................... 4.1-2 - 4.2 STRUCTURE PLOTTING.................................................. 4.2-1 4.2.1 General Rules.................................... 4.2-2 4.2.2 Plot Request Card Descriptions........................... 4.2-4 4.2.3 Examples of Structure Plot Requests...................... 4.2-22 4.2.4 Summary of Structure Plot Request Packet Cards............ 4.2-25 4.3 X-Y OUTPUT.............................................. 4.3-1 4.3.1 X-Y Plotter Terminology.................................. 4.3-1 4.3.2 Parameter Definition Cards............................... 4.3-2 4.3.3 Command Operations Cards.................................. 4.3-8 4.3.4 Examples of X-Y Output Request Packets.................... 4.3-17a ix (2/13/81)

TABLE OF CONTENTS (Continued) Section Page No. 4.3.5 Summary of X-Y Output Request Packet Cards................ 4.3-20 4.3.6 XY-Plots for SORT1 Output................................. 4.3-21 4.3.7 Dynamic Response Output on TABLED1 Bulk Data Cards........ 4.3-23 4.3.8 Superelement Plot Packet Example.......................... 4.3-24 5. DIRECT MATRIX ABSTRACTION - 5.1 GENERAL.............................................. 5.1-1 5.2 DMAP RULES........................................................ 5.2-1 5.2.1 DMAP Module Calling Sequences.......5.................... 5.2-1 5.2.2 Sequences of Modules................................... 5.2-10 5.2.3 Executive DMAP Instructions............................... 5.2-18 5.2.4 ALTERs................................... 5.2-30 5.3 EXAMPLES......................................................... 5.3-1 5.3.1 DMAP Example to Display Outputs........................... 5.3-1 5.3.2 DMAP Example to Solve a System of Equations.............. 5.3-5 5.3.3 Example of DMAP to Run a Static Analysis................... 5.3-7 5.3.4 Example of DMAP to Solve an Iteration Problem.......... 5.3-11 5.3.5 DMAP Example of Real Eigenvalue Extraction................ 5.3-12 5.3.6 DMAP Example to Raise the Matrix [Q] to a Power R......... 5.3-15 5.3.7 DMAP Example....................................... 5.3-17 5.4 DMAP MODULE DESCRIPTIONS............................................ 5.4-1 6. DIAGNOSTIC MESSAGES 6.1 RIGID FORMAT DIAGNOSTIC MESSAGES................................... 6.1-1 6.1.1 (Deleted) 6.1.2 (Deleted) 6.1.3 (Deleted) 6.1.4 (Deleted) 6.1.5 Rigid Format Error Messages for Buckling Analysis......... 6.1-1 6.1.6 (Deleted) x (2/13/81)

TABLE OF CONTENTS (Continued) Section Page No. 6.1.7 Rigid Format Error Messages for Direct Complex Eigenvalue Analysis....................................... 6.1-2 6.1.8 (Deleted) 6.1.9 (Deleted) 6.1.10 Rigid Format Error Messages for Modal Complex Eigenvalue Analysis.................................... 6.1-2 6.1.11 (Deleted) 6.1.12 (Deleted) 6.1.13 Variable Parameters Used for Error Exits............... 6.1-3 6.2 NASTRAN SYSTEM AND USER MESSAGES................................ 6.2-1 6.2.1 Preface Messages............................ 6.2-2 6.2.2 Executive Module Messages............................. 6.2-17 6.2.3 Functional Module Messages................................ 6.2-21 6.3 UNNUMBERED NASTRAN SYSTEM AND USER MESSAGES................. 6.3-1 6.3.1 (Deleted) 6.3.2 Superelement Processing Messages.......................... 6.3-1 6.3.3 Element Generation Messages............................. 6.3-9 6.3.4 Matrix Assembly and Update Messages....................... 6.3-15 6.3.5 Matrix Operations Messages............................... 6.3-16 6.3.6 Data Recovery Messages.................................... 6.3-19 6.3.7 Data Base Messages........................................ 6.3-20 6.3.8 Plotter Messages........................................ 6.3-21 6.3.9 NASTRAN Card Messages.................................... 6.3-23 6.3.10 Input File Processor Messages...................... 6.3-24 6.3.11 Grid Point Singularity Processor Messages.......6..... 6.3-25 6.3.12 Machine-Dependent Error Messages.......................... 6.3-25 xi (2/13/81)

STRUCTURAL ELEMENTS Zelem | a - ma ofse _aele m, b,nbX Gb) BeapeenintPlane a fre or (a) Bea1a element coordinate system a Een Fore (O ~~~ ~1.-7 (2,0,0281 w offset AGA Yb o fset Beam eleen Grid pointnts I LV Figure 1. Beam, Bar and Bend Element Coordinate Systems and Element Forces. 1.3-17 (2/2/81)

STRUCTURAL MODELING x (c) Bar element coordinate system. ^^^^^~Plane 2 a y \^ ^v T Mla Mlb x x FX Plane I vi (d) Bar element forces z v M2a M^b Plane 2 V2 Figure 1. Beam, Bar and Bend Element Coordinate Systems and Element Forces. (Cont.) 1.3-18 (2/2/81)

STRUCTURAL ELEMENTS G3 G6 G5 / 05 \ Yelem n = const. G4 G1l Xelem G2 (c) TRIA6 Yelem nG81 G7 Gi G6 ti1 2Ji~~ f._ xelem id) QUAD8 Figure 4. Shell Element Coordinate Systems. (Cont.) 1.3-23 (6/15/81)

STRUCTURAL MODELING y /IyFy y F / xy/r xy (b) Moments Figure 5. Forces, Moments, and Stresses in Plate Elements. 1.3-24 (10/30/81)

STRUCTURAL ELEMENTS y c~X i, jIf a a1 x x txy xy (c) Stresses Figure 5. Forces, Moments, and Stresses in Plate Elements. (Cont.) z Uz - Displacement Coordinates W __~RINGAX B. N. ^^^^U - Element Coordinates " RINGAX A Figure 6. Geometry for Conical Shell Element. 1.3-25 (10/30/81)

STRUCTURAL MODELING z,w,Fz ~~~w299, (A 0', F - e (Material Orientation) Figure 7. Triangular Ring Element Coordinate System. zw' IZ 0 x ~ a (Material Orientation) r,u Figure 8. Trapezoidal Ring Element Coordinate System. 1.3-26 (2/2/81)

STRUCTURAL ELEMENTS z z - Zbasic GS rm G GG3 r Xbasic radial Figure 9. The TRIAX6 Solid of Revolution Element. 1.3-27 (2/2/81)

STRUCTURAL MODELING HS 3 s''':. kSaue 9SI: x x x x x x x.alsue4 21,aH |: x x x. x | x |x x x x x SS3j1S9SW x> x x _x Z.LB0S xxx x xx x x sassalS u i0d p._9_.Jno.uo0 old e x x x x x x x x x Ix sa6p3.Ua.. u 3. I.^. I. |'x: O urns do..I I'x I - | I LvqOL9, urs aco. I x x x x x x x x x x xaeLd|3 as |...,. ^ ei r.- r-. x LU t ~ Avwy.re 3 3 00. aJoj:. _ __.. _'ssasj s I, 3 CI | S V ZI - r ea m, cn en - _d. O V. cm x..-... SSajS co cm.. p,4- c trI; LoruralU i wU i i. i i ii 2- o oau xx x x x x C s Y I x >C W I' i ~ ~,...: ii. ~ ~ ~ ~ 0. L/)iii <CQlOLe1 4 I, W Laa W We d ~ 3d04^0S~l x x x x x x x xx x oldoaufOs ) X x x:D pzamoaQ..... x > x - ssau dd X x-,- x ~ l~^uau....xx x ~x x xxxl;'w ~ o i''' ssau~) |PS x x x x x x - x * N ~ * S' m m",~,......'>C "d~i ~ acau:~2 < ^ i,: i iii "'1.3-28 (10/30/,1'. t

STRUCTURAL ELEMENTS HS3W9SW x x x x x x AJSUeJ.,9aH x x x x x x SS3U1S9SW x x Z0lHS x x x x x x sassa^,S'40 d x OOLd Jnoauoo x.0...... _ ajn | outSx x x x xx X x x x c m lL,,,,,de, sa6p3.uauOWL3 E w F 1eaH x x x x x x w IP~oL9 >n In urns 3saJj x x x x x x _, C ~ o0 ^ xai~ulo, r..% o m n~s'g c 0 __ __ 0. _ -,' --------------------- tA j (O - adO' i -. —-. ---—, —-..- - ---- - - --- x u- EC u _. Pxa l9dm3 a F) - in C., (v) iU) C, n E.4. sse.-s - - o ~ = e.. eo r, a.,-, S 0g o o i 1. ssa8 deS C. Io III U Ix Y:;oCosi x > xu-, X < ~ o x x x x x x g 5Po uo onpuoo.. 1| I= ^9 x_ x x x x g& Luau3 _ I I i I III i i ssauss x x1x x.x *' ^ uo I dOXa X X X o tdoosua v L 9a t dASSe.A dI x x C X 6 M S X X X 4J ______a _______ _ - JvutuONII r,J dssoaui s L u u a (. >C' S. En- -r- a) 3~doawosI x x xX x x 3~ I~~~rta~~~~J~~s C~~xv~ x x4-) x (,C 0'C x?~v o J,-2 deaUILUON ar SOUau43 EUC5 En eaU I. LUO,. cC4 ~~d3'ma~ u~(.3-29 3 10/ 3( / 8. x )x x x ~s..S.. -........... @3 4'a ad,(. - -- 4,.-, 1.3-29 (10/30/81)

~L__Y_~~~tY)_~IIIIIIyV hr A s rpA h, DA-Tp~ DE CK i

NASTRAN DATA DECK 2.1 GENERAL DESCRIPTION OF DATA DECK The input deck begins with the required resident operating system control cards. The type and number of these cards will vary with the installation. Instructions for the preparation of these control cards are given in Section 7.6 of the Application Manual. The operating system control cards are followed by the MSC/NASTRAN Data Deck, which consists of the following three sections: - 1. Executive Control Deck' 2. Case Control Deck -6 3. Bulk Data Deck The Executive Control Deck and the Case Control Deck both have free-field formats. Only columns 1 thru 72 are used for data. Any information in columns 73 thru 80 will appear in the printed echo, but the data will not be used by the program. As explained in Section 2.4.1, limited use is made of the data in columns 73 thru 80 for the Bulk Data Deck. The NASTRAN card is used to change the default values for certain operational parameters, such as buffer size or the number of data lines printed per page. The NASTRAN card is optional, but, if present, it must be the first card of the NASTRAN Data Deck. The NASTRAN card is a freefield card (similar to cards in the Executive Control Deck). Its format is as follows: NASTRAN keywords = value, keyword2 = value,... These keywords set specific cells in the /SYSTEM/ common block. These same cells may be set dynamically (between DMAP statements) by the PARAM module (see Section 5). The number in parentheses after the keywords is the associated number of the system cell. The most frequently used keywords are as follows: 1. BUFFSIZE(1) - Defines the number of words in a GIN0 buffer. Usually this value is standardized at any particular installation. However, the desired value may be different from the default value. In any event, in a series of related runs, such as restarts, the same BUFFSIZE must be used for all of the runs. 2.1-1 (7/21/80)

NASTRAN DATA DECK 2. NLINES(12) - Defines the number of data lines per page. This value is automatically set for a particular installation. However, when nonstandard paper is used, the use of this keyword will produce proper paging. 3. HIC0RE(57) - Defines the maximum region request in words for UNIVAC and VAX. 4. DAYLIMIT(7) - Defines the number of entries allowed in dayfile (CDC). 5. PREF0PT(31) - Selects preprocessors for execution in the Preface. The following option is available: PREF0PT Preprocessor 2 MSGMESH (limited release) 6. FILES(45) - Used to declare permanent GIN0 disk files on CDC machines. 7. MPYAD(66) - Deselects multiply/add methods as follows: MPYAD Selection O Default value - all methods are available for program selection. 1 Deselect Method 1 Nontranspose 2 Deselect Method 1 Transpose 4 Deselect Method 2 Nontranspose 8 Deselect Method 2 Transpose 16 Deselect Method 3 Nontranspose 32 Deselect Method 3 Transpose 64 Deselect Method 4 Nontranspose 128 Deselect Method 4 Transpose Several methods can be deselected by using the sum of the values for deselecting the individual methods. 8. SDC0MP(69) - Controls options for the symmetric decomposition routine during the review of the matrix prior to beginning the actual decomposition operation as follows: SDC0MP Action O Default - Print up to 50 messages for null columns and zero diagonals. Do not perform decomposition if there are any null columns or zero diagonal terms. 1 Terminate execution when first null column is encountered. 2 Suppress printing of message when a null column is encountered. 4 Terminate execution when first zero diagonal term is encountered. 2.1-2 (7/21/80)

GENERAL DESCRIPTION OF DATA DECK SDCOMP Action 8 Suppress printing of message when a zero diagonal term is encountered. 16 Place 1.0 in diagonal position for all null columns and proceed with the decomposition. 32 In case of zero diagonal terms, proceed with the decomposition unless the leading minor is zero. 64 Exit after execution of preface for symmetric decomposition. Combinations of program action are indicated by using the sum of the values for the individual actions. If the ninth parameter of the DECOMP module is set to a non-default value, it overrides any value used with the SDCOMP keyword. This parameter is set to non-default values in some of the solution sequences. It can be changed by a DMAP alter. See the DECOMP module description in Section 5.4. 9. REAL(81) - Specifies directly the amount of OPEN CORE (single precision words) to be used. This information is used on virtual machines to prevent thrashing. The value for REAL should never be set to more than the amount of real address space available to the user. If REAL a 0, the OPEN CORE is calculated from the REGION request on the operating system JOB card. Default value for IBM machines is 100000. The default value on the VAX 11/780 depends on the working set (WSL). REAL - (WSL-200) 128. 10. HEAT(56) - the integer value of 1 specifies that the current problem is heat transfer analysis rather than structural analysis. 11. DBSET(None) - defines data base subsets as follows: DBSET a = (DBij, DBij,...), where A is an integer from 1 thru 15 and DBij are data base file names. Default values are i = 1 and ij = 01. 12. I0RATE(84) - I/0 operations per second. MSC/NASTRAN will, by default, consider I/O time as well as CPU time in the internal decision logic used to select the most cost-effective method of matrix multiplication (MPYAD), forward-backward substitution (FBS), etc. The decision logic utilizes the following relation to estimate the number of equivalent CPU seconds required to perform the desired matrix operation by each of the available methods: iOBLOCKS T - CPU + IRABLCKS where T = Equivalent CPU time CPU - CPU time for matrix operation IOBLOCKS = Number of transfers between memory and disk IORATE = 20 blocks/second (default) The default value of 20 blocks/second for IORATE is approximately the read/write rate for commonly used disk devices. Other values of IORATE may be more effective in minimizing computer charges because of the characteristics of local charge algorithms. If IORATE is explicitly set to zero on the NASTRAN card, method selection will be based solely on estimated CPU time. 2.1-3 (12/31/80)

NASTRAN DATA DECK 13. UDC0MP(69) - Controls options for the unsymmetric decomposition routine during the review of the matrix prior to beginning the actual decomposition operation as follows: UDC0MP Action 0 Default - Print up to 50 messages if both row and column are null. Do not perform decomposition if there are any null rows and columns. 1 Terminate execution if both row and column are null. 2 Suppress printing of message if both row and column are null. 16 Place 1.0 in diagonal position if both row and column are null and proceed with the decomposition. 64 Exit after execution of preface for unsymmetric decomposition. Combinations of program action are indicated by using the sum of the values for the individual actions. 14. FBS0PT(70) - Selects Forward/Backward Substitution methods as follows: FBS0PT Selection -1 Select Method 1-FBS. 0 Program.selects method based on minimum of I/O + CPU time. +1 Select Method 2 FBS. The effect of using keywords on the NASTRAN card is to place information in selected locations within a common block called SYSTEM. Information can also be placed in the SYSTEM common block by specifying the location directly. Section 2.4.1.8 of the MSC/NASTRAN Programmer's Manual gives a table correlating the keywords with their location number. For example, specifying SYSTEM (81) - 30000 is equivalent to REAL = 30000. The location number method can also be used to change information in the SYSTEM common block during the course of a run using DMAP statements. (See the PARAM module description in Section 5.4). The checkpoint dictionary and other punched output are transmitted to the MSC/NASTRAN PUNCH file. The PUNCH file can be sent to an alternate FORTRAN unit by use of SYSTEM(64)=(n), where (n) is the number of an alternate unit. See the 0UTPUT2 module description in Section 5.4 for a discussion of the machine-dependent FORTRAN units. The checkpoint dictionary can be separated from the other punched output by sending it to the F0RTRAN unit defined by SYSTEM(83)=(m). The remainder of the punched output, such as MSGVIEW data, STRESS(PUNCH) files, etc., will go to the MSC/NASTRAN PUNCH file or the unit defined on SYSTEM(64). 2.1-4 (12/31/80)

GENERAL DESCRIPTION OF DATA DECK The Executive Control Deck begins with the ID card and ends with the CEND card, as indicated in Figure 1. It identifies the job and the type of solution to be performed. It also declares the general conditions under which the job is to be executed, such as maximum time allowed, type of system diagnostics desired, restart conditions, and whether or not the job is to be checkpointed. If the job is to be executed with a rigid format, the actual rigid format is declared along with any alterations to the rigid format that may be desired. If Direct Matrix Abstraction is used, the complete DMAP sequence must appear in the Executive Control Deck. The Executive Control Cards and examples of their use are described in Section 2.2. The Case Control Deck begins with the first card following CEND and ends with the card, BEGIN BULK, as indicated in Figure 1. It defines the subcase structure for the problem, makes selections from the Bulk Data Deck, and makes output requests for printing, punching and plotting. A general discussion of the functions of the Case Control Deck and a detailed description of the cards used in this deck are given in Section 2.3. The special requirements of the Case Control Deck for each solution sequence are discussed in Section 3. The Bulk Data Deck begins with the card following BEGIN BULK and ends with the card, ENDDATA, as indicated in Figure 1. It contains all of the details of the structural model and the conditions for the solution. The BEGIN BULK and ENDDATA cards must be present even though no new bulk data is being Introduced into the problem or all of the bulk data is coming from an alternate source, such as user-generated input. The format of the BEGIN BULK card is free-field. The ENDDATA card must begin in column 1 or 2. Generally speaking, only one structural model can be defined in the Bulk Data Deck. However, some of the bulk data, such as cards assocated with loading conditions, constraints, direct input matrices, transfer functions and thermal fields may exist in multiple sets. Only sets selected in the Case Control Deck will be used in any particular solution. The bulk data cards are described in Section 2.4. Comment cards may be inserted in any of the parts of the Data Deck. These cards are identified by a $ in column one. Columns 2-72 may contain any desired text. All data cards must be punched using the character set shown in the table on the next page. The EBCDIC character set may also be used. Any EBCDIC characters are automatically translated into the character set shown in the table on the next page. The EBCDIC character card punch configurations are shown in parentheses for the five characters that differ from the standard character set. 2.1-5 (12/31/80)

NASTRAN DATA DECK Character Card Punch(es) Character Card Punch(es) EBCDIC Punch(es) blank blank N 11-5 0 0 0 11-6 1 1 I p 11-i 2 2 Q 11-8 3 3 R 11-9 4 4 S 0-2 5 5 T 0-3 6 6 U 0-4 7 7 V 0-5 8 8 W 0-6 9 9 X 0-7 A 12-1 Y 0-8 B 12-2 Z 0-9 C 12-3 $ 11-3-8 D 12-4 / 0-1 E 12-5 + 12 (12-6-8) F 12-6 11 G 12-7 ( 0-4-8 (12-5-8) H 12-8 ) 12-4-8 (11-5-8) I 12-9 (6-8) 4-8 (5-8) J 11-1. 3-8 (6-8) K 11-2 0-3-8 L 11-3 12-3-8 M 11-4 * 11-4-8 2.1-6 (12/31/80)

GENERAL DESCRIPTION OF DATA DECK Figure 1. MSC/NASTRAN Data Deck. 2.1-7 (12/31/80)

C oJTeOU t ( O Ct ( ^^- --— " — c

NASTRAN DATA DECK 2.2 EXECUTIVE CONTROL DECK The format of the Executive Control cards is free field. The name of the operation (e.g., CHKPNT) begins in column 1 and is separated from the operand by one or more blanks. The fields in the operand are separated by commas, and may be integers (Ki) or alphanumeric (Ai) as indicated in the following control card descriptions. The first character of an alphanumeric field must be alphabetic followed by up to 7 additional alphanumeric characters. Blank characters may be placed adjacent to separating commas if desired. The individual cards are described in Section 2.2.1 and examples follow in Section 2.2.2. 2.2.1 Executive Control Card Descriotions - ID AI, A2 Required. Al, A2 -- Any legal alphanumeric fields chosen by the user for problem identification. RESTART Al, A2, Kl/K2/K3, K4, Required for Restart. Al, A2 -- Fields taken from ID card of previously checkpointed problem. Kl/K2/K3, K4 -- Month/Day/Year, Time that Problem Tape was generated. The complete restart dictionary consists of the RESTART card followed by one physical card for each file checkpointed and each reentry point. The restart dictionary is automatically punched when operating in the checkpoint mode. The complete restart dictionary is one logical card. Each continuation card begins with a sequence number. Each type of continuation card will be documented separately. 1. Basic continuation card NO, DATABLOCK,FLAG=Y,REEL=Z,FILE=W where:;NO is the sequence number of the card. The entire dictionary must be in sequence by this number. DATA8LOCK is the name of the data block referenced by this card. FLAG=Y defines the status of the data block where Y = 0 is the normal case and Y = 4 implies this data block is equivalenced to another data block. In this case (FLAG=4) the file number points to a previous data block which is the "actual" copy of the data. 2.2-1

NASTRAN DATA DECK REEL=Z specifies the reel number. Z=1 even if the Problem Tape is multi-reel. FILE=W specifies the GINO (internal) file number of the data block on the Problem Tape. A zero value indicates the data block is purged. For example: 1,GPL,FLAGS=O,REE=l,FILE=7 says data block GPL occupies file 7 of reel 1. 2,KGG,FLAGS=4,REEL=l,FILE=20 says KGG is equivalenced to the data block which occupies file 20. (Note that FLAGS-4 cards usually occur in at least pairs at the equivalenced operation is at least binary.) 3,USETD,FLAGS=O,REEL=1,FILE=O implies USETD is purged. 2. Reentry point card: NO,REENTER AT OMAP SEQUENCE NUMBER N where NO is the sequence number of. the card. N is the sequence number associated with the OMAP instruction at which the problem will restart. This value may be changed by adding a final such card (i.e., only the last such card is operative). This may be necessary when restarting from a Rigid Format to a DMAP sequence (to print a matrix, for example). The REENTER card provides an easy way to accomplish this. There are four types of restarts: Unmodified Restart; Modified Restart; Rigid Format Switch; and Pseudo Modified Restart. The function of the reentry point is different in each case. On an unmodified restart the program continues from the reentry point. On a modified restart, modules which must be run to process the modified data but which are ahead of the reentry point are executed first. The program then continues from the reentry point. On a Rigid Format Switch (going from a Rigid Format to another), the reentry point is meaningless in that it was determined for another DMAP sequence. In this case, the data blocks available are consulted to determine the proper sequence of modules to run. A Pseudo Modified Restart (defined by the existence of only changes to output producing data such as plotter requests) is treated like a modified restart. The type of restart is implied by the changes made in the NASTRAN Data Deck. No explicit request for a Darticular kind of restart is required. See Section 3.1 for additional information. 3. End of dictionary card: $ END OF CHECKPOINT DICTIONARY T:iis card is simply a comment card but is punched to signal the end of the dictionary for user convenience. The program does not need such a card. Terminations associated with non-NASTRAN failures (operator intervention, maximum time, etc.) will not have such a card punched. 2.2-2

EXECUTIVE CONTROL DECK REENTER K1 Optional. This card directs the program to reenter at DMAP statement K1. This card is only effective on restart and overrides the directive given in the restart dictionary. Note that this card is ignored (as is any reentry point) on Rigid Format Switches. TRUNCATE K1 Optional. The restart dictionary will be truncated (ignored) after card number K1. The truncate card must precede the restart dictionary. CHKPNT Al or CHKPNT Al, A2 Optional. Al — YES if problem is to be checkpointed, NO if problem is not to be checkpointed. Defaul tNO. A2 — DISK if checkpoint file is on direct access device. If the DISK option is used, the user must instruct the resident operating system to permanently catalog the checkpoint file. APP A Optional. A -- DISPLACEMENT indicates one of the Displacement solution sequences (default). A -- HEAT indicates one of the Heat Transfer Approach solution sequences. A - DMAP indicates Direct Matrix Abstraction Approach (DMAP). Default if a DMAP sequence is submitted. -* SOL K1 [,Ki] or SOL An [,Ki] Required when using a solution sequence. K1 — Number of Solution Sequence (see table below). Ki — Subset numbers for solution K1, default = O. Multiple subsets can be used by separating the subset numbers with commas. (See Section 3.1 for allowable subsets.) An — Name of Rigid Format (see table below) Rigid Formats Ki An Parameter 3 NORMAL MODES M0DES 5 BUCKLING BUCKLING -> 24 STATICS STATICS 25 None - old Normal Modes None 26 DIRECT FREQUENCY RESPONSE DFREQ 27 DIRECT TRANSIENT RESPONSE DTRAN 28 DIRECT C0MPLEX EIGENVALUES DCEIG 29 MODAL COMPLEX EIGENVALUES MCEIG 30 MODAL FREQUENCY RESPONSE MFREQ 31 MODAL TRANSIENT RESPONSE MTRAN 45 AERODYNAMIC FLUTTER FLUTTER 46 AERODYNAMIC RESPONSE AERO 47 CYCLIC STATICS CYC STATICS 48 CYCLIC MODES CYC MODES 2.2-3 (12/31/80)

NASTRAN DATA DECK Solution Sequences K1 An Parameter 60 SUPERELEMENT M0DEL CHECKOUT CHECKOUT 61 SUPERELEMENT STATICS SE STATICS 62 ALTERNATE SUPERELEMENT STATICS ALT SE STATICS 63 SUPERELEMENT NORMAL MODES SE MODES 64 GEOMETRIC NONLINEAR GN0LIN 65 SUPERELEMENT BUCKLING SE BUCKLING 66 MATERIAL NONLINEAR MNOLIN 67 SUPERELEMENT DIRECT COMPLEX EIGENVALUES SE OCEIG 68 SUPERELEMENT DIRECT FREQUENCY RESPONSE SE DGREQ 69 SUPERELEMENT DIRECT TRANSIENT RESPONSE SE DTRAN 70 SUPERELEMENT M0DAL COMPLEX EIGENVALUES SE MCEIG 71 SUPERELEMENT MODAL FREQUENCY RESPONSE SE MFREQ 72 SUPERELEMENT MODAL TRANSIENT RESPONSE SE MTRAN 74 STEADY NONLINEAR HEAT TRANSFER None 75 SUPERELEMENT AERODYNAMIC FLUTTER SE FLUTTER 76 SUPERELEMENT AEROELASTIC RESPONSE SE AERO 77 CYCLIC BUCKLING CYC BUCKLING 78 CYCLIC FREQUENCY RESPONSE CYC FREQ 81 SUPERELEMENT CYCLIC STATICS None 83 SUPERELEMENT CYCLIC MODES None 88 SUPERELEMENT CYCLIC FREQUENCY RESPONSE None 89 SUPERELEMENT TRANSIENT HEAT TRANSFER None ALTER K1, K2 Optional. K1,K2 -- First and last DMAP instructions of series to be deleted and replaced with any following DMAP instructions. Alter numbers do not need to be in sort; overlaps are resolved on a first come basis. ALTER K Optional. K -- Input any following DMAP instructions after statement K. > TIME K, M Required if K is greater than one minute. K -- Maximum allowable execution time in minutes. K is a real or integer number; thus, 1.5 is equivalent to 90 seconds. The default for K is one minute. M - Reserved for an I/0 limit. ENDALTER Required when using ALTER unless the alter package ends with the CEND or RESTART card. Indicates end of DMAP alters. BEGIN Required when using DMAP approach. Indicates beginning of DMAP sequence and trnrlilates any DMAP alters. END Required when using DMAP approach. Indicates end of DMAP sequence. This card must begin in column 1 and contain a blank in column 4. 2.2-4 (6/1/81)

EXECUTIVE CONTROL DECK - DIAG K Optional request for diagnostic output. K = 1 Dump memory when non-preface fatal message is generated. K - 2 Print File Allocation Table (FIAT) following each call to the File Allocator. K = 3 Print status of the Data Pool Dictionary (DPD) following each call to the Data Pool Housekeeper. K = 4 Print the DMAP cross reference. K = 5 Print BEGIN time on-line for each functional module on the operator's console. K = 6 Print END time on-line for each functional module on the operator's console. K = 7 Print eigenvalue extraction diagnostics for complex determinant method. K = 8 Print matrix trailers as the matrices are generated and data blocks used from 0PTP. K = 9 Echo checkpoint dictionary as cards are punched. K = 10 Use alternate nonlinear loading in TRD. Replace Nn+l by 1 (Nn+ + Nn + Nn_). K = 11 Not used. K = 12 Print eigenvalue extraction diagnostics for complex inverse power. K = 13 Print open core length. (REAL on the VAX). -%bK = 14 Print the DMAP Sequence (NASTRAN S0URCE PR0GRAM C0MPILATI0N). K = 15 Trace GIN0 0PEN/CL0SE operations. Print table trailers. K = 16 Trace real inverse power eigenvalue extraction operations. K = 17 Punch the DMAP sequence that is compiled. K = 18 Do not trace Heat Transfer iterations in module SSGHT. In Aeroelastic Analysis, print internal grid points picked by SET2 cards. K = 19 Print data for MPYAD and FBS method selection. K = 20 Print Data Base Manager fetch/store messages. K = 21 Not used. K = 22 Not used. K = 23 Compute strains rather than stresses for QUAD4, QUAD8, TRIA3, and TRIA6. Stresses will not be calculated for any elements. If ELFORCE is used with DIAG 23 the calculated forces will be incorrect. K = 24 Print EXCP counts for files on IBM. Print secondary file allocation messages in run log for UNIVAC. K = 25 0utput internal plot diagnostics. K = 26 Print material orientation angle 8 for the QUAD4, QUAD8, TRIA3 or TRIA6 element corresponding to the user-supplied material orientation coordinate system ID. K = 27 Dump Input File Processor (IFP) tables. 2.2-5 (6/1/81)

NASTRAN DATA DECK DIAG K (Cont.) K = 28 Punch the link specification table (Deck XBSBD.) K = 29 Process link specification table update deck. A module may be forced to run in a specified link by placing the following two cards after the ENDDATA card at the end of the data deck. M0DNAME, EPNAME, I ENDDATA where M0DNAME is the name of the module, EPNAME is the entry point name and I is the link number. K = 30 Punch the XSEMii decks (i.e., set ii via DIAG 1-15). (After Link 1, this turns on BUG output.) K = 31 Print link specification table and module properties list data. K = 32 Not used. K = 33 Print Hencky-von Mises equivalent stresses instead of maximum shear stress. K = 34 Turn off plot line optimization. K = 35 Turn off automatic RFL on CDC 6000 and Univac computers. See Section 4.3.7 for effect on XYTRAN module. K = 36 Turn off copying to FIAT on RESTARTs. K = 37 Disable superelement congruence test option. K = 38 Print material angles for QUAD4, QUAD8, TRIA3, and TRIA6 elements. K = 39 Trace module FA1 operations. K = 40 Show cards modified by 1VARY cards. K 41 Trace buffer manager on CDC 7600. K = 42 Turn on page locking on the VAX. K = 43 Not used. K = 44 Print mini dump for fatal errors and suppress user message exit. K = 45 Print FIAT before and after every module execution. K = 46 Set by program if CHKPNT is YES. K = 47 Generate automatic CHKPNT for all outputs if CHKPNT is YES. (Not needed for Solution Sequences.) K = 48 Do not calculate stresses for elements using composite materials. K = 49 Print stress output for composite materials from SDR2. K = 50 Print nonlinear convergence values per iteration. K = 51 Print nonlinear g-size displacement, unbalanced, and error vectors per iteration. K = 52 thru Not used. K = 64 2.2-6 (6/1/81)

EXECUTIVE CONTROL DECK Multiple options may be selected by using multiple integers separated by commas. Other options and other rules associated with the DIAG card which primarily concern the programmer can be found in Section 6.11.3 of the Programmer's Manual. ECH00FF Optional. Suppresses echo of following cards such as a large ALTER package or a restart dictionary. ECH00N Optional. Reactivates the echo after an ECH00FF card. CEND Required. Indicates end of executive control cards. The ID card must appear first and CEND must be the last card of the Executive Control Deck. C0MPILER = Al, A2, A3, A4 Optional. This card is an alternate for DIAGs 4, 14, and 17. Al -- LIST if the DMAP Library is to be printed. Default is NOLIST. A2 -- DECK if a DMAP with any embedded ALTERs is to be punched. Default is N0DECK. A3 -- REF if a cross reference table of the DMAP sequence is to be printed. Default is N0REF. A4 -- N0G0 if the job is to be terminated after the preface even when no errors are detected. Default is GO. 2.2.2 Executive Control Deck Examples 1. Cold start, no checkpoint, rigid format, diagnostic output. ID MYNAME, BRIDGE23 S0L 24 TIME 5 DIAG 8,13 CEND 2. Cold start, checkpoint, rigid format. ID PERS0NZZ, SPACECFT CHKPNT YES S0L 24 TIME 15 CEND 3. Restart, no checkpoint, rigid format. The restart dictionary indicated by the brace is automatically punched on previous run in which the CHKPNT option was selected by the user. ID J0ESCHM0E, PR0JECTX S0L3 TIME10 RESTART PERS0NZZ, SPACECFT, 05/13/67, 15, 1, XVPS, FLAGS=O, REEL=1, FILE=6 2, REENTER AT DMAP SEQUENCE NUMBER 7 3, GPL, FLAGS=O, REEL=1, FILE=7 $ END 0F CHECKP0INT DICTI0NARY CEND 2.2-7 (12/31/80)

NASTRAN DATA DECK 4. Cold start, no checkpoint, DMAP. User-written DMAP program is indicated by braces. ID IAMO07, TRYIT BEGIN $ {DMAP statements go heret END $ TIME 8 CEND 5. Restart, checkpoint, altered rigid format, diagnostic output. ID GOODGUY, NEATDEAL CHKPNT YES C0MPILER LIST, REF SOL 3 TIME 15 ALTER 56 MATPRN KGGX// $ RESTART BADGUY, NOSHOW, 05/09/68, 25, 1, XVPS, FLAGS=O, REEL=1, FILE=6 2, REENTER AT DMAP SEQUENCE NUMBER 7 3, GPL, FLAGS=O, REEL=1, FILE=7 $ END OF CHECKPOINT DICTIONARY CEND 6. Superelement statics cold start or restart, diagnostic output.* ID SUPER, COLD $ DIAG 8 TIME 5 S0L 61 CEND *The 0PTP-NPTP files may be used in the superelement rigid formats, but they are not required. If an OPTP file is used, its checkpoint dictionary will consist of two cards, the first starting with the word "RESTART" and the second with "1,XVPS" etc. Only the first card should be used on restart. The second card causes the values of the parameters at the time of termination of the prior run to influence the present run, and may lead to restart errors. 2.2-8 (12/31/80)

CAsE OAJTTOL 1 ) C ( %VW- quo

NASTRAN DATA DECK 2.3 CASE CONTROL DECK The Case Control Deck defines the subcase structure for the problem, makes selections from the Bulk Data Deck, and makes output requests. A summary of all Case Control cards is given inSection 2.3.2. The individual cards are described in Section 2.3.3. Structure plotter output requests and curve plot requests are treated separately in Sections 4.2 and 4.3 respectively. 2.3.1 Subcase Definition In general, a separate subcase is defined for each loading condition and/or each set of constraints. Subcases may also be used in connection with output requests, such as requesting different output for each mode in a real eigenvalue problem. Only one level of subcase definition is provided. All items placed above the subcase level (ahead of the first subcase) will be used for all following subcases unless overridden within the individual subcase. In statics problems, provision has been made for the combination of the results of several subcases. This is convenient for studying various combinations of individual loading conditions and for the superposition of solutions for symmetrical and antisymmetrical boundaries. The following examples of Case Control Decks indicate typical ways of defining subcases: 1. Static analysis with multiple loads. DISPLACEMENT - ALL MPC a 3 SUBCASE 1 SPC = 2 TEMPERATURE(L0AD): 101 L0AD = 11 SUBCASE 2 SPC = 2 DEFORM = 52 L0AO = 12 SUBCASE 3 SPC = 4 LOAD = 12 SUBCASE 4 M1PC = 4 SPC = 4 Four subcases are defined in this example. The displacements at all grid points will be printed for all four subcases. MPC: 3 will be used for the first three subcases and will be overridden by MPC = 4 in the last subcase. Since the constraints are the same for subcases 1 and 2 and the subcases are contiguous, the static solutions will be performed simultaneously. In subcase 1, thermal load 101 and external load 11 are internally superimposed, as are the external and deformation 2.3-1

NASTRAN DATA DECK loads in subcase 2. In subcase 4, the static loading will result entirely from enforced displacements of grid points. 2. Linear combination of subcases. SPC = 2 SET 1 - 1 THRU 10,20,30 DISPLACEMENT = ALL STRESS = 1 SUBCASE 1 L0AD = 101 0L0AD = ALL SUBCASE 2 LOAD = 201 OL0AD = ALL SUBC0M 51 SUBSEQ = 1.0,1.0 SUBC0M 52 SUBSEQ = 2.5,1.5 Two static loading conditions are defined in subcases 1 and 2. SUBC0M 51 defines the sum of subcases 1 and 2. SUBCOM 52 defines a linear combination consisting of 2.5 times subcase 1 plus 1.5 times subcase 2. The displacements at all grid points and the stresses for the element numbers in SET will be printed for all four subcases. In addition, the nonzero components of the static load vectors will be printed for subcases 1 and 2. 3. Statics problem with one plane of symmetry. SET 1 a 1,11,21,31,51 SET 2 = 1 THRU 10, 101 THRU 110 DISPLACEMENT 1 ELF0RCE = 2 SYM 1. SPC = 11 L0AD - 21 0L0AD - ALL SYM 2 SPC 12 L0AD = 22 SYMC0M 3 SYMC0M 4 SYMSEQ 1.0,-1.0 Two SYM subcases are defined in subcases 1 and 2. SYMC0M 3 defines the sum and SYMC3M 4 the difference of the two SYM subcases. The nonzero components of the static load will be printed for subcase 1 and no output is requested for subcase 2. The displacements for the grid point numbers in set 1 and the forces for elements in set 2 will be printed for subcases 3 and 4. 2.3-2

CASE CONTROL DECK 4. Use of REPCASE in statics problems. SET 1 = 1 THRU 10, 101 THRU 110, 201 THRU 210 SET 2 = 21 THRU 30, 121 THRU 130, 221 THRU 230 SET 3 = 31 THRU 40, 131 THRU 140, 231 THRU 240 SUBCASE 1 L0AD = 10 SPC = 11 DISP = ALL SPCF0RCE = 1 ELFORCE = 1 REPCASE 2 ELFORCE = 2 REPCASE 3 ELFORCE = 3 This example defines one subcase for solution and two subcases for output control. The displacements at all grid points and the nonzero components of the single-point forces of constraint along with forces for the elements in SET 1 will be printed for SUBCASE 1. The forces for elements in SET 2 will be printed for REPCASE 2 and the forces for elements in SET 3 will be printed for REPCASE 3. 5. Use of MODES in eigenvalue problems. METHOD = 2 SPC = 10 SUBCASE 1 DISP = ALL STRESS = ALL M0DES = 2 SUBCASE 3 DISP = ALL In this example, the displacements at all grid points will be printed for all modes. The stresses in all elements will be printed for the first two.modes. 2.3.1.1 Superelement Subcase Structure Subcases must be defined for each superelement to specify load, constraint, and output selection. The syntax rules described above also apply for superelements, i.e.,requests above the subcase level are over-ridden by requests inside the subcase level. 1. Static analysis with multiple loads, two superelements. DISP = ALL ESE = ALL SPC = 1000 SUBCASE 11 SUPER = 10,1 LOAD = 1000 SUBCASE 12 SUPER = 10,2 L0AD = 1002 SUBCASE 21 SUPER = 20,1 L0AD = 1000 2.3-3 (2/15/78)

NASTRAN DATA DECK Assume that the model generated in step 4 was used as a starting point. Superelement 100 is an alternate version of superelement 10. This run generates superelement 100, then reduces it and assembles and solves the residual structure. Superelement 10 is excluded (see SET 3), so that the alternate version of the structure Is assembled, while maintaining the matrices for superelement 10 in the data base. If it is decided to return to the original structure on the next run, superelement 10 can be included and superelement 100 excluded by the sequence SET 2 = 0 SET 3 = 100 SEMA = 2 SELA = 2 SEEX = 3 2.3.2 Case Control Card Summary This section contains a summary of all Case Control cards under the following headings: 2.3.2.1 Subcase Definition 1. Output Request Delimiters 2. Subcase Delimiters 3. Subcase Control 2.3.2.2 Data Selection 1. Static Load Selection 2. Dynamic Load Selection 3; Constraint Selection 4. Thermal Field, Selection 5. Dynamic Solution Conditions 6. Direct Input Matrix Selection 7. Nonlinear Analysis 2.3.2.3 Output Selection 1. Output Control 2. Set Defintion 3. Physical Set Output Requests 4. Solution Set Output Requests 2.3.2.4 Superelement Control 2.3-4 (10/30/81)

CASE CONTROL DECK 2.3.2.1 Subcase Definition 1. Output Request Delimiters 0UTPUT Beginning of printer output request (optional) 0UTPUT(PL0T) Beginning of structure plotter output request 0UTPUT(P0ST) Beginning of grid point stress output requests OUTPUT(XY0UT).or Beginning of curve plotter output request 0UTPUT(XYPL0T) OUTPUT(CARDS) Suppresses processing of Bulk Data cards ENDCARDS Reactivates processing of Bulk Data cards 2. Subcase Delimiters SUBCASE Defines beginning of subcase SUBC0M Defines beginning of subcase which is a linear combination of preceding subcases SYM Defines beginning of symmetry subcase SYMC0M Defines beginning of subcase for making symmetry combinations REPCASE Defines beginning of subcase for additional output requests 3. Subcase Control SUBSEQ Defines coefficients for linear combination in SUBC0M SYMSEQ Defines coefficients for symmetry combination in SYMC0M M0DES Specifies repetition of subcase for eigenvalue problems MASTER Defines a new MASTER subcase 2.3.2.2 Data Selection 1. Static Load Selection L0AD Selects static loading condition DEF0RM Selects element deformation set 2. Dynamic Load Selection DL0AD Selects dynamic loading conditions L0ADSET Selects static load sets for use in dynamic loading N0NLINEAR Selects nonlinear loading condition for transient response 2.3-5 (10/30/81)

- I L' r- " " * NASTRAN DATA DECK 3. Constraint Selection - SPC Selects set of single-point constraints MPC Selects set of multipoint constraints AXISYMMETRIC Specifies boundary conditions for conical shell problems or hydroelastic problems DYSM Selects symmetry option in cyclic symmetry 4. Thermal Field Selection TEMPERATURE(LOAD) Selects temperature set for static thermal load TEMPERATURE(MATERIAL) Selects temperature set for temperature-dependent material properties TEMPERATURE Selects temperature set for both static load and material properties TEMPERATURE(ESTIMATE) Selects temperature for heat transferproblems 5. Dynamic Solution Conditions DYNRED Selects conditions for dynamic reduction METHOD Selects conditions for real eigenvalue analysis CMETHOD Selects conditions for complex eigenvalue analysis SDAMPING Selects table for determination of modal damping FREQUENCY Selects frequency set for frequency response analysis RANDOM Selects power spectral density functions for random analysis IC Selects initial conditions for direct transient response TSTEP Selects time steps to be used for integration in transient response 6. Direct Input Matrix Selection B2GG Selects damping matrices before constraints are applied B2PP Selects damping matrices K2GG Selects stiffness matrices before constraints are applied K2PP Selects stiffness matrices M2GG Selects mass matrices before constraints are applied M2PP Selects mass matrices TFL Selects transfer functions P2G Selects load matrices MFLUID Selects MFLUID Bulk Data cards 2.3-6 (10/30/81)

CASE CONTROL DECK 7. Nonlinear Analysis NLPARM Selects parameters for nonlinear analysis 2.3.2.3 Output Selection 1. Output Control - TITLE Specifies text for first line on each printed page SUBTITLE Specifies text for second line on each printed page LABEL Specifies text for third line on each printed page LINE Sets the number of data lines per printed page - default is installation dependent, usually 50 lines/page MAXLINES Sets maximum number of output lines - default is 100,000 - ECH0 Selects echo options for Bulk Data Deck - default is sorted echo SKIP0N Defines cards in the Case Control Deck that are not to be processed SKIP0FF Resumes processing of cards in the Case Control Deck ECH00FF Suppresses echo of Case Control Deck ECH00N Reactivates echo of Case Control Deck PAGE Causes a page eject in echo of Case Control Deck PL0TID Specifies BCD identification for first frame of plotter output 2. Set Definition SET Defines lists of point numbers, element numbers, frequencies, or time steps for use in output requests 0FREQUENCY Selects a set of frequencies for output requests in frequency response problems TSTEP Selects time steps for output requests in transient response problems 0TIME Selects a subset of time steps for output requests in transient response problems PARTN Specifies a list of grid points for partitioning operations 3. Physical Set Output Requests --. b DISPLACEMENT Requests displacements for a set of physical points or VECTOR SVECT0R Requests solution set modal eigenvector output VEL0CITY Requests velocities for a set of physical points ACCELERATI0N Requests accelerations for a set of physical points 2.3-7 (10/30/81)

NASTRAN DATA DECK ELF0RCE Requests the forces for a set of structural elements or F0RCE STRESS Requests the stresses for a set of structural elements - or ELSTRESS GPSTRESS Requests the stresses at grid points STRFIELD Requests the stresses at grid points for postprocessing GPF0RCE Requests grid point force balance for a set of grid points ESE Requests strain energy for a set of elements FLUX Requests the flux and gradient for a set of heat transfer elements - SPCF0RCES Requests the single-point forces of constraint for a set of points - 0L0AD Selects a set of applied loads for output in static analysis PRESSURE Requests hydroelastic pressure output for a set of points THERMAL Requests temperatures for a set of points in heat transfer analysis NOUTPUT Requests physical output in cyclic symmetry problems 4. Solution Set Output.Requests HARM0NICS Specifies harmonics for output requests for conical shell and hyrdroeleastic problems H0UTPUT Requests harmonic output in cyclic symmetry problems SDISPLACEMENT Requests the displacements of the independent components for a set of or SVECT0R points or modal coordinates SVELOCITY Requests the velocities of the independent components for a set of points or modal coordinates SACCELERATI0N Requests the accelerations of the independent components for a set of points or modal coordinates NLL0AD Selects a set of nonlinear loads for output in transient response MPRES Requests the pressure for selected surface elements 2.3.2.4 Superelement Control SUPER Specifies superelement identification number and the load sequence number SEMGENERATE Specifies the superelement identification numbers for which stiffness, mass, and damping matrices will be generated SEMASSEMBLE Specifies the superelement identification numbers for which stiffness matrices will be assembled and reduced SELGENERATE Specifies the superelement identification numbers for which load vectors will be generated 2.3-8 (10/30/81)

CASE CONTROL DECK SELASSEMBLE Specifies the superelement identification numbers for which the static load, mass or damping matrices will be assembled and reduced SEALL Combines the functions of SEMGENERATE, SEMASSEMBLE, SELGENERATE and SELASSEMBLE SEKREDUCE Specifies the superelement identification numbers for which stiffness matrices will be assembled and reduced' SEMREDUCE Specifies the superelement identification numbers for which the mass and damping matrices,will be assembled and reduced SELREDUCE Specifies the superelement identification numbers for which the static load matrices will be assembled and reduced SEFINAL Specifies the superelement identification number of the last superelement to be assembled SEEXCLUDE Specifies the superelements which are not to be assembled ("excluded") intotheir downstream superelement SEDR Specifies the superelement identification numbers for which data recovery will be performed 2.3.3 Case Control Card Descriptions The format of the Case Control cards is free-field. In presenting general formats for each card embodying all options, the following conventions are used: 1. Upper-case letters must be punched as shown. 2. Lower-case letters indicate that a substitution must be made. 3. Braces { } indicate that a choice of contents is mandatory. 4. Brackets [ ] contain an option that may be omitted or included by the user. 5. Underlined options or values are the default values. 6. Physical card consists of information punched in columns 1 through 72 of a card. Most Case Control cards are limited to a single physical card. 7. Logical card may have more than 72 columns with the use of continuation cards. If the first four characters of a mnemonic are unique relative to all other Case Control cards, the characters following can be omitted. 2.3-8a (10/30/81)

NASTRAN DATA DECK 2.3.3 Case Control Card Descriptions The format of the case control cards is free-field. In presenting general formats for each card embodying all options, the following conventions are used: 1. Upper-case letters must be punched as shown. 2. Lower-case letters indicate that a substitution must be made. 3. Braces { } indicate that a choice of contents is mandatory. 4. Brackets [ ] contain an option that may be omitted or included by the user. 5. Underlined options or values are the default values. 6. Physical card consists of information punched in columns 1 thru 72 of a card. Most case control cards are limited to a single physical card. 7. Logical card may have more than 72 columns with the use of continuation cards. If the first four characters of a mnemonic are unique relative to all other case control cards, the characters following can be omitted. 2.3-8b (2/2/81)

CASE CONTROL DECK Case Control Data Card DISPLACEMENT - Displacement Output Request. Description: Requests form and type of displacement vector output. Format and Examole(s): F/ST PRINT REAL \ ALL DISPLACEMENT RI PRINT IMAG, PLOT - LSORT2 PUNCH PHASE /J (NE) DISPLACEMENT = 5 DISPLACEMENT(REAL) = ALL DISPLACEMENT(SORT2, PUNCH, REAL) = ALL Option Meaning SORT1 Output will be presented as a tabular listing of grid points for each load, frequency, eigenvalue, or time, depending on the rigid format. SORT1 is not available on Transient problems (where the default is SORT2). S0RT2 Output will be presented as a tabular listing of frequency or time for each grid point. SORT2 is available only in Transient and Frequency Response problems. PRINT The printer will be the output media. PUNCH The card punch will be the output media. REAL or IMAG Requests real and imaginary output on Complex Eigenvalue or Frequency Response problems. PHASE Requests magnitude and phase (0.0O < phase < 360.0~) on Complex Eigenvalue or Frequency Response problems. PLOT Generate, but do not print, displacement data. ALL Displacements for all points will be output. NONE Displacement for no points will be output. n Set identification of previously appearing SET card. Only displacements of points whose identification numbers appear on this SET card will be output (integer > 0). Remarks: 1. Both PRINT and PUNCH may be requested. 2. On a Frequency Response problem any request for Physical Set SORT2 causes all Physical Set output to be SORT2. 3. VECTOR and PRESSURE are alternate forms and are entirely equivalent to DISPLACEMENT. 4. DISPLACEMENT = NONE allows overriding an overall output request. 5. The PL0T option is used when curve plots are desired in the magnitude/phase representation and no printer output request is present for magnitude/phase representation. 2.3-19

CASE CONTROL DECK Case Control Data Card ECHO Bulk Data Echo Request Description: Requests echo of Bulk Data Deck. Format and Example(s): ECHO = [S0RT,UNS0RT, BTH,NONE,PUNCH,N0S0RT,SORT(cdnl,cdn2,... ) ] ECHO = NOSORT,UNS0RT ECHO = BOTH ECHO = SORT(MAT1,PARAM) Option Meaning SORT Sorted echo will be printed. UNSORT Unsorted echo will be printed. B0TH Both sorted and unsorted echo will be printed. NONE No echo will be printed. PUNCH The sorted Bulk Data Deck will be punched. N0SORT Blocks the printing of the sorted echo. This allows the user to obtain a listing of changed cards with CHKPNT but not get a sorted echo. S0RT(cdnl,cdn2) This requests that card types cdnl and cdn2 only be echoed (along with their line numbers). Up to 26 card types may be echoed. (Continuation cards are included.) Remarks: 1. If no ECHO card appears, a sorted echo will be printed. 2. If CHKPNT YES, a sorted echo will be printed unless ECHO = NONE or ECHO = NOSORT. 3. Comment cards will be punched at the front of the sorted deck if ECHO = PUNCH. 4. Portions of the unsorted deck can be selectively echoed by including the cards ECH00N and ECH00FF at various places within the Bulk Data Deck and requesting an unsorted echo. ECHOOFF stops the unsorted echo until an ECHOON card is encountered. Many such pairs of cards may be used. 5. If the SORT(cdnl,cdn2,...) option is used, and the card images are recovered from the 0PTP, the continuation cards will not be printed. 6. If the SORT(cdnl,cdn2,...) option is used, it must be last in the list. 2.3-25 (5-15-80)

CASE CONTROL DECK Case Control Data Card. ELFORCE - Element Force Output Request. Description: Requests form and type of element force output. Format and Examole(s): ELF0CRCE L'S0RT1 PRINT REAL ALL ELF0RCE SRT' PUNCH Ig - <n O PHASE ELFORCE = ALL ELF0RCE(REAL, PUNCH, PRINT) = 17 ELF0RCE = 25 Option Meaning SORT1 Output will be presented as a tabular listing of elements for each load, frequency, eigenvalue, or time, depending on the rigid format. SORTI is not available on Transient problems (where the default is S0RT2). S0RT2 Output will be presented as a tabular listing of frequency or time for each element type. SORT2 is available only in Transient and Frequency Response problems. PRINT The printer will be the output media. PUNCH The card punch will be the output media. REAL or Requests real and imaginary output on Complex Eigenvalue or Frequency Response IMAG problems. PHASE Requests magnitude and phase (0.00 < phase < 360.0~) on Complex Eigenvalue or Frequency Response problems. ALL Forces for all elements will be output. NONE Forces for ino elements will be jutput. n Set identification of a previously appearing SET card. Only forces of elements whose identification numbers appear on this SET card will be output (Integer > 0). Remarks: 1. Both PRINT and PUNCH may be requested. 2. ALL should not be used in a Transient problem. 3. On a Frequency Response problem, any request for Physical Set S0RT2 output causes all Physical Set output to be SORT2. 4. FORCE is an alternate form and is entirely equivalent to ELFORCE. 5. ELF0RCE = NONE allows overriding an overall request. 2.3-31

CASE CONTROL DECK Case Control Data Card ELSTRESS - Element Stress Output Request. Description: Requests form and type of element stress output. Format and Example(s): [/S0RT1 PRINT REAL ] ALL ELSTRESS SRT2' PUNL PHASE T] n. PHASE / N0NES ELSTRESS = 5 ELSTRESS = ALL ELSTRESS(SORT1, PRINT, PUNCH, PHASE) 15 ELSTRESS(PLOT) = ALL Option Meaning SORT1 Output will be presented as a tabular listing of elements for each load, frequency, eigenvalue, or time, depending on the rigid format. S0RT1 is not available on Transient problems (where the default is S0RT2). 50RT2 Output will be presented as a tabular listing of frequency or time for each element type. S0RT2 is available only in Transient and Frequency Response problems. PRINT The printer will be the output media. PUNCH The card punch will be the output media. REAL or IMAG Requests real and imaginary printout on Complex Eigenvalue or Frequency Response problems. PHASE Requests magnitude and phase (0.0~ < phase < 360.0~) on Complex Eigenvalue or Frequency Response problems. PLOT Generates stresses for requested set but no printer output. ALL Stresses for all elements will be output. n Set identification of a previously appearing SET card (Integer > 0). Only stresses for elements whose identification numbers appear on this SET card will be output. NONE Stress for no elements will be output. Remarks: 1. Both PRINT and PUNCH may be requested. 2. ALL should not be used in a Transient problem. 3. On a Frequency Response problem, any request for Physical'Set S0RT2 output causes all Physical Set output to be S0RT2. 4. STRESS is an alternate form and is entirely equivalent to ELSTRESS. 5. ELSTRESS = N0NE allows overriding an overall output request. 6. The PLOT option is used when contour plots of stresses are requested, but no printer output of stresses is desired. 2.3-33

CASE CONTROL DECK Case Control Data Card LABEL - Output Label Description: Defines a BCD label which will appear on the third heading line of each page of printer output. Format and Example(s): LABEL = {Any BCD data} LABEL = DEMONSTRATION PROBLEM Remarks: 1. LABEL appearing at the subcase level will label output for that subcase only. 2. LABEL appearing before all subcases will label any outputs which are not subcase dependent. 3. If no LABEL card is supplied, the label line will be blank. 4. LABEL information is also placed on plotter output as applicable. 2.3-51 (2/2/81)

CASE CONTROL DECK Case Control Data Card LOAD - External Static Load Set Selection. Description: Selects the external static load set to be applied to the structural model. Format and Example(s): L0AO = n L0AD = 15 Option Meaning n Set identification of at least one external load card and hence must appear on at least one F0RCE, F0RCE1, F0RCE2, F0RCEAX, GRAY, M0MAX, M0MENT, M0MENTlj_M0MENT2, L0AD, PL0AD, PL0ADi, PL0AD2, PL0AD3, PL0ADX, QV0L, QVECT, QHBDY, QBDY1, QBDY2, PRESAX, RF0RCE, or SLOAD card (Integer > 0) Remarks: 1. The above static load cards will not be used by NASTRAN unless selected in Case Control. 2. A GRAV card cannot have the same set identification number as any of the other loading card types. If it is desired to apply a gravity load along with other static loads, a LOAO bulk data card must be used. 3. L0AD is only applicable in statics, inertia relief, differential stiffness, buckling, piecewise linear problems, linear and nonlinear heat transfer problems. 4. The total load applied will be the sum of external (L0AO), thermal (TEMP(L0AD)), element deformation (DEFORM) and constrained displacement (SPC) loads. 5. Static, thermal and element deformation loads should have unique set identification numbers. 2.3-55

CASE CONTROL DECK Case Control Data Card MPC - Multipoint Constraint Set Selectiun. Descriotion: Selects the multipoint constraint set to be applied to the structural model. Format and Examole(s): MPC = n MPC = 17 Ootion Meaning n "n" is the set identification of a Multipoint-Constraint Set and hence must appear on at least one MPC or MPCADD card. (Integer > 0). Remarks: MPC or MPCADD cards will not be used by NASTRAN unless selected in Case Control. 2.3-65

CASE CONTROL DECK Case Control Data Card 0L0AD Applied Load Output Request Description: Requests form and type of applied load vector output. Format and Example(s): (L0AD S0RT1I PRINT REAL nALL f O 0RT27 *' PUNCH' i AG N 0NE PHASE/_J 0LOAO = ALL 0L0AD(SRT1,PHASE) = 5 Option Meaning S0RT1 Output will be presented as a tabular listing of grid points for each load, frequency, eigenvalue, or time, depending on the rigid format. S0RT1 is not available on transient problems (where the default is S0RT2). S0RT2 Output will be presented as a tabular listing of frequency or time for each grid point. SORT2 is available only in transient and frequency response problems. PRINT The printer will be the output media. PUNCH The card punch will be the output media. REAL or IMAG Requests real and imaginary output on complex eigenvalue or frequency response problems. PHASE Requests magnitude and phase (0.0~ < phase < 360.0~) on complex eigenvalue or frequency response problems. ALL Applied loads for all points will be output. (SORT1 will only output nonzero values.) N0NE Applied load for no points will be output. n Set identification of previously appearing SET card. Only loads on points whose identification numbers appear on this SET card will be output (Integer > 0). Remarks: 1. Both PRINT and PUNCH may be requested. 2. On a frequency response problem any request for S0RT2 causes all output to be S0RT2. 3. In a statics problem a request for S0RT2 causes loads at all points (zero and nonzero) to be output. 4. 0L0AO = NONE allows overriding an overall output request. 5. In the statics superelement rigid formats, only externally applied loads are printed, and not loads transmitted from upstream superelements. Transmitted loads can be obtained with GPF0RCE requests. 2.3-73 (2/15/78)

CASE CONTROL DECK Case Control Data Card OUTPUT - Output Request Delimiter Description: Delimits the various output requests, structure plotter, curve plotter, and pri nter/punch Format and Example(s): (FPLOT \ POST \ OUTPUT iXYeUT 1 XYPLOTI L CARDS / OUTPUT OUTPU-T(PLOT) OUTPUT(XYOUT) Option Meaning No qualifier Beginning of printer output request - this is not a required card. PLOT Beginning of structure plotter request. This card must precede all structure plotter control cards. POST Beginning of grid point stress output request. This card must precede all postprocessor control cards. XTOUT or Beginning of curve plotter request. This card must precede all curve plotter XYPL0T control cards. XYPLOT and XYOUT are entirely equivalent. CARDS Beginning of the deck of cards which will be placed on Data Block FORCE in 20A4 format. These cards have no format rules applied. This package must terminate with the card ENDCARDS (starTTng in column 1). Remarks: 1. The structure plotter request and the curve plotter request must be at the end of the Case Control Deck. 2. The delimiting of a printer request its completely optional. 3. The OUTPUT(CARDS) packet is currently used by the MSGSTRES Module. See the MSGMESH Analyst's Guide for details. 2.3-77 (10/30/81)

CASE CONTROL DECK Case Control Data Card SET - Set Definition Card Description: 1. Lists identification numbers (point or element) for output requests. 2. Lists the frequencies for which output will be printed in Frequency Response Probl ems. Format and Example(s): 1) SET n = {il[,12, 13 THRU i4 EXCEPT i5, 169 17, i8 THRU ig]} SET 77 - 5 SET 88 = 5, 6, 7, 8, 9, 10 THRU 55 EXCEPT 15, 16, 77, 78, 79, 100 THRU 300 SET 99 = 1 THRU 100000 2) SET n u {r1[, r2, r3, r4]} SET 101 = 1.0, 2.0, 3.0 SET 105 = 1.009, 10.2, 13.4, 14.0, 15.0 3) SET n = ALL Option Meaning n Set identification (Integer > 0). Any set may be redefined by reassigning its identification number. Sets inside SUBCASE delimiters are local to the SUBCASE. i!, i2, etc. Element or point identification number at which output is requested. (Integer > 0).- If no such identification number exists, the request is ignored. i3 THRU i4 Output at set identification numbers 13 thru 14 (i4 > 13). EXCEPT Set identifiction numbers following EXCEPT will be deleted from output list as long as they are in the range of the set defined by the immediately preceding THRU. r1, r2, etc. Frequencies for output (Real > 0.0). The nearest solution frequency will be output. EXCEPT and THRU cannot be used. ALL All members of the set will be processed. Remarks: 1. A SET card may be more than one physical card. A comma (,) at the end of a physical card signifies a continuation card. Commas may not end a set. 2. Set identification numbers following EXCEPT within the range of the THRU must be in ascending order. 2.3-103 (7/1/80)

CASE CONTROL DECK Case Control Data Card SPC - Single-Point Constraint Set Selection. Description: Selects the single-point constraint set to be applied to the structural model. Format and Example(s): SPC = n SPC = 10 Option Meaning n Set identification of a single-point constraint set and hence must appear on a SPC, SPC1 or SPCADD card (Integer > 0). Remarks: SPC, SPC1 or SPCADD cards will not be used by NASTRAN unless selected in Case Control. 2.3-105

CASE CONTROL DECK Case Control Data Card SPCF0RCES - Single-Point Forces of Constraint Output Request. Description: Requests form and type of Single-Point Force of constraint vector output. Format and Example(s): SPCF0RCES (S0RT1 PRINT REAL ALL PHASE NONE SPCF0RCES =5 SPCF0RCES(S0RT2, PUNCH, PRINT, IMAG) - ALL SPCF0RCES(PHASE) = N0NE Option Meaning S0RT1 Output will be presented as a tabular listing of grid points for each load, frequency, eigenvalue, or time, depending on the rigid format. S0RT1 is not available on Transient problems (where the default is S0RT2). S0RT2 Output will be presented as a tabular listing of frequency or time for each grid point. S0RT2 is available only in Transient and Frequency Response problems. PRINT The printer will be the output media. PUNCH The card punch will be the output media. REAL or Requests real and imaginary output on Complex Eigenvalue or Frequency Response IMAG problems. PHASE Requests magnitude and phase (0.0~ < phase < 360.0~) on Complex Eigenvalue or Frequency Response problems. ALL Single-Point forces of constraint for all points will be output. (SORT1 will only output nonzero values.) NONE Single point forces of constraint for no points will be output. n Set identification of previously appearing SET card. Only single-point forces constraint for points whose identification numbers appear on this SET card will be output (Integer > 0). Remarks: 1. Both PRINT and PUNCH may be requested. 2. On a Frequency Response problem any request for solution set S0RT2 output causes all solution set output to be S0RT2. 3. In a Statics problem'a request for S0RT2 causes loads at all points (zero and nonzero) to be output. 4. SPCF0RCES = NONE allows overriding an overall output request. 2.3-107

CASE CONTROL DECK Case Control Data Card STRESS - Element Stress Output Request. Description: Requests form and type of element stress output. Format and Example(s): RSTRESS Rc PRINT REAL \1 IALL 1 STRESS \ |UIMAG, PL0T I PHASE l NONEI STRESS = 5 STRESS = ALL STRESS(S0RT1, PRINT, PUNCH, PHASE) = 15 STRESS(PL0T) = ALL Option Meaning S0RT1 Output will be presented as a tabular listing of elements for each load, frequency, eigenvalue, or time, depending on the rigid format. SORT1 is not available on Transient problems (where the default is SORT2). SORT2 Output will be presented as a tabular listing of frequency or time for each element type. SORT2 is available only in Transient and Frequency Response problems. PRINT The printer will be the output media. PUNCH The card punch will be the output media. REAL or IMAG Requests real and imaginary printout on Complex Eigenvalue or Frequency Response problems. PHASE Requests magnitude and phase (0.0~ < phase < 360.0~) on Complex Eigenvalue or Frequency Response problems. PLOT Generates stresses for requested set but no printer output. ALL Stresses for all elements will be output. n Set identification of a previously appearing SET card (Integer > 0). Only stresses for elements whose identification numbers appear on this SET card will be output. NONE Stress for no elements will be output. Remarks: 1. Both PRINT and PUNCH may be requested. 2. ALL should not be used in a Transient problem. 3. On a Frequency Response problem, any request for solution set 5a0Rt2 output causes alTsoluti-on set -output- to0'be -RT2-. 4: ELSTRESS is an alternate form and is entirely equivalent to STRESS. 5. STRESS = NONE allows overriding an overall output request. 6. The PL0T option is used when contour plots of stresses are requested, but no printer output of stresses is desired. 2.3-109

NASTRAN DATA DECK Case Control Data Card STRFIELD Grid Point Stress Output Request for Postprocessing Description: Requests grid point stresses for postprocessing Format and Example(s): STRFIELD = {ALL} STRFIELD = ALL STRFIELD = 21 Option Meaning ALL Grid point stress requests for all SURFACEs defined in the 0UTPUT(POST) section will be saved for postprocessing. n Set identification number of previously appearing SET card. Only SURFACEs whose identification numbers appear on this SET card and in the 0UTPUT(P0ST) section will be included in the grid point stress output request for postprocessing. 2.3-110 (10/30/81)

CASE CONTROL DECK Case Controi Data Card SUBCASE - Subcase Delimiter. Description: Delimits and identifies a subcase. Format and Examole(s): SUBCASE n SUBCASE 101 Option Meani n n Subcase identification number (Integer > 0). Remarks: 1. The subcase identification number, n, must be strictly increasing (i.e., greater than all previous subcase identification numbers). 2. Plot requests and RAND0M requests refer to n. 2.3-111

CASE CONTROL DECK Case Control Data Card SUBTITLE - Output Subtitle Description: Defines a BCD subtitle which will appear on the second heading line of each page of printer output. Format and Example(s): SUBTITLE = {Any BCD data} SUBTITLE = PROBLEM NO. 5-1A Remarks: 1. SUBTITLE appearing at the subcase level will title output for that subcase only. 2. SUBTITLE appearing before all subcases will title any outputs which are not subcase dependent. 3. If no SUBTITLE card is supplied, the subtitle line will be blank. 4. SUBTITLE information is also placed on plotter output as applicable. 2.3-117 (2/13/81)

CASE CONTROL DECK Case Control Data Card TITLE - Output Title. Description: Defines a BCD title which will appear on the first heading line of each page of NASTRAN printer output. Format and Example(s): TITLE = { Any BCD data } TITLE = **$// ABCDEFGHI.... $ Remarks: 1. TITLE appearing at the subcase level will title output for that subcase only. 2. TITLE appearing before all subcases will title any outputs which are not subcase dependent. 3. If no TITLE card is supplied, the title line will contain data and page numbers only. 4. TITLE information is also placed on NASTRAN plotter output as applicable. 2.3-137

>ECIc

NASTRAN DATA DECK 2.4 BULK DATA DECK The primary NASTRAN input medium is the Bulk Data card. These cards are used to define the structural model and various pools of data which may be selected by Case Control at execution time. For large problems the Bulk Data Deck may consist of several thousand cards. In order to minimize the handling of large numbers of cards, provision has been made in NASTRAN to store the Bulk Data on the Problem Tape, from which it may be modified on subsequent runs. See the MSGMESH Analyst's Guide for information on the use of the MSGMESH Program for the generation of Bulk Data cards. For any cold start, the entire Bulk Data Deck must be submitted. Thereafter, if the original run was checkpointed, the Bulk Data Deck exists on the Problem Tape in sorted form where it may be modified and reused on restart. On restart the Bulk Data cards contained in the Bulk Data Deck are added to the Bulk Data contained on the Old Problem Tape. Cards are removed from the Old Problem Tape by the use of a delete card. Cards to be deleted are indicated by inserting a Bulk Data card with a slash (/) in column one and the sorted bulk data sequence numbers in fields two and three. All bulk data cards in the range of the sequence numbers in fields two and three will be deleted. In the case where only a single card is deleted, field three may be left blank. The Bulk Data Deck may be submitted with the cards in any order as a sort is performed prior to the execution of the Input File Processor. It should be noted that the machine time to perform this is minimized for a deck that is already sorted. The sort time for a badly sorted deck will become significant for large decks. The user may obtain a printed copy of either the unsorted or sorted bulk data by selection in the Case Control Deck. A sorted echo is necessary in order to make modifications on a secondary execution using the Problem Tape. This echo is automatically provided unless specifically suppressed by the user. 2.4.1 Format of Bulk Data Cards Bulk Data cards may be prepared either in the form of 8 or 16 column fields, or as free-field cards. The 8-column format is Indicated in the following diagram. 2.4-1 (5-15-80)

NASTRAN DATA DECK Small Field Bulk Data Card la 2 3 4 5 6 7 8 9 10a ^-8- 8- 8- 88-8 -8- -8 -8- l-88- — W 8The mnemonic is punched in field 1 beginning in column 1. Fields 2-9 are for data items. The only limitations in data items are that they must lie completely within the designated field, have no inbedded blanks, and must be of the proper type, i.e., blank, integer, real, double precision, or BCD. All real numbers, including zero, must contain a decimal point. The SEQGP and SEQEP cards use the decimal point in a special way. A blank will be interpreted as a real zero or integer zero as required. Real numbers may be encoded in various ways. For example, the real number 7.0 may be encoded as 7.0,.7E1, 0.7+1, 70.-1,.70+1, 7+0, etc. A double precision number must contain both a decimal point and an exponent with the character D such as 7.0DO. Double precision data values are allowed only in a few situations, such as on the PARAM card. BCD data values consist of one to eight alphanumeric characters, the first of which must be alphabetic. Normally, field 10 is reserved for optional user identification. However, in the case of continuation cards, field 10 (except column 73 which is not referenced) is used in conjunction with field 1 of the continuation card as an identifier and, hence, must contain a unique entry. The continuation card contains the symbol + in column 1 followed by the same seven characters that appeared in columns 74-80 of field 10 of the card that is being continued. This allows the data to be submitted as an unsorted deck. BCD values used as continuation mnemonics cannot contain the symbols *, =, or $. Continuation mnemonics may be generated automatically for continuation cards which are in sorted order. This option is triggered by placing a + in column 73 of the parent card and leaving field 1 of the continuation cards blank. The program will then generate the unique continuation mnemonics. This option may not be used for bulk data cards submitted on restart runs. 2.4-2 (12/15/80)

BULK DATA DECK The small field data card should be adequate for most.applications. Occasionally, however, the input is generated by another computer program or is available in a form where a wider field would be desirable. For this case, the larger field format with a 16-character data field is provided. Each logical card consists of two physical cards as indicated in the following diagram: Large Field Bulk Data Card la 2 3 4 5 10a *-8 *zJ 16 * 16 16 16 8 lb 6 7 8 9 10b t - o ~16 16 16 -- -216 — ( — 2.4-2a (12/15/80)

BULK DATA DECK The large field card is denoted by placing the symbol * after the mnemonic in field la and some unique character configuration in the last 7 columns of field 10a. The second physical card contains the symbol * in column 1 followed by the same seven characters that appeared after column 73 in field 10a of the first card. The second card may in turn be used to point to a large or small field continuation card, depending on whether the continuation card contains the symbol * or the symbol + in column 1. The use of multiple and large field cards are illustrated in the following examples: Small Field Card with Small Field Continuation Card TYPE I I I I I 10.QED123 +ED123. Large Field Card TYPE* OED124 *X3124 - Large Field Card with Large Field Continuation Card TYPE* ___ ED301 *~EDl ________________________________________________ _0 * ED302 ____________________________ _________________ED305 ~*ED305 _ Large Field Card Followed by a Small Field Continuation Card and a Large Field Continuation Card TYP,*.......... __________________ __________ED4G2 *f46 Z it EDD4L 1 +ED421| 1 IEgSl *ED361 [, ED291 *ED291 __________________ Small Field Card with Larce Field Continuation Card TYPE' ______________I' ED632 *ED0632______________________ __________________ ) ED204 fED204 | | |. l l l_________________________o 2.4-3

NASTRAN DATA DECK The free-field format also includes a limited data generation capability. The following rules apply to the use of the free-field format: 1. Only-small-field cards can be created. 2. Continuation cards are not generated automatically. 3. Data items must be separated with a comma or one or more blanks. A comma or equal sign must appear in the first ten columns of the card. 4. Duplication of fields from the preceding card is accomplished by coding the symbol. 5. Duplication of all trailing fields from the preceding card is accomplished by coding the symbol == 6. Incrementing a value from the previous card is indicated by coding *(i), where i is the value of the increment. 7. Repeated replication is indicated by coding =(n), where n is the number of card images to be generated using the values of the increments on the preceding card. The following is an example of the use of free-field cards: GRID,101,17,1.0,10.5,,17,3456:,*(1 ),=,*(o.2) $ =(3) The above free-field cards will generate the following bulk data cards in the 8-column format: 1 2 3 4 5 6 7 8 9 10 GRID 101 17 1.0 10.5 __ 17 3456 __ GRID 102 17 1.2 10.5 __ 17 3456 GRID 103 17 1.4 10.5 ___ 17 3456 GRID 104 17 1.6 10.5 __ 17 3456 GRID 105 17 1.8 10.5 __ 17 3456 __ 2.4.2 Bulk Data Card Summary This section contains a summary of all Bulk Data cards under the following headings: Section 2.4.2.1 Geometry - 1. Grid Points 4. Fluid Points - 2. Coordinate Systems 5. Axisymmetry 3. Scalar Points 6. Cyclic Symmetry 7. Superelement/Substructuring 2.4-4 (1/15/77)

BULK DATA DECK Section 2.4.2.2 Elements - 1. Elastic Line Elements 7. Mass Elements " 2. Elastic Surface Elements 8. Damping Elements - 3. Elastic Solid Elements 9. Fluid Elements -- 4. Elastic Scalar Elements 10. Heat Transfer 5. Axi symmetric Elements 11. Dummy Elements 6. Rigid Elements 12. Nonlinear Elements Section 2.4.2.3 Materials - 1. Isotropic 4. Stress Dependent 2. Anisotropic 5. Fluid 3. Temperature Dependent Section 2.4.2.4 Constraints - 1. Single-Point Constraints 4. Free-Body Supports 2. Multipoint Constraints 5. Component Mode Boundary Conditions 3. Partitioning 6. User Sets Section 2.4.2.5 Loads - 1. Static Loads 2. Dynamic Loads 3. Heat Transfer Section 2.4.2.6 Problem Control 1. Buckling Analysis 5. Frequency Response 2. Eigenvalue Analysis 6. Random Response 3. Cyclic Symmetry 7. Transient Response 4. Dynamics 8. Nonlinear Analysis Section 2.4.2.7 Miscellaneous -* 1. Comments 5. Tabular Input 2. Delete 6. Output Control --- 3. Parameters 7. Matrix Assembly -> 4. Direct Matrix Input 2.4-5 (1/15/81)

NASTRAN DATA DECK Section 2.4.2.8 Bulk Data Generator (MSGMESH) 1. Grid Point Generation 2. Element Generation 3. Temperature Definition 4. Load and Constraint Definition 5. MSGMESH Control See the MSGMESH Analyst's Guide for information regarding the generation of Bulk Data cards. Section 2.4.2.9 Variance Analysis See Section 2.8 of the MSC/NASTRAN Application Manual for a discussion of variance analysis. Section 2.4.2.10 Aeroelastic Analysis 1. Aerodynamic Elements 2. Aerodynamic Data 3. Aerodynamic-Structure Connection Section 2.4.2.11 Deleted bulk data cards This section contains a summary of all bulk data cards which have been deleted. 2.4.2.1 Geometry 1. Grid Points -- GRID Grid point location, coordinate system selection GRIDB Grid point location on boundary of axisymmetric fluid problem - GRDSET Default options for GRID cards - SEQGP Grid and scalar point number resequencing 2. Coordinate Systems C0RDiC Cylindrical coordinate system definition C0RDiR Rectangular coordinate system definition C0RDiS Spherical coordinate system definition BEAM0R Orientation default for CBEAM BAROR Orientation default for CBAR 2.4-6 (5-15-80)

BULK DATA DECK 3. Scalar Points SP0INT Scalar point definition EPOINT Extra point definition for dynamics GRIDS Scalar degree of freedom for acoustic cavity analysis GRIDF Scalar degree of freedom for acoustic cavity analysis -- SEQGP Grid and scalar point number resequencing SEQEP Extra point number resequencing 4. Fluid Points GRIDB Grid point location on RINGFL RINGFL Circle (fluid point) definition GRIDF Scalar degree of freedom for acoustic cavity analysis GRIDS Scalar degree of freedom for acoustic cavity analysis FREEPT Surface point location for data recovery PRESPT Pressure point location for data recovery FSLIST List of fluid points (RINGFL) on free surface boundary SLBDY List of slot points on interface between fluid and radial slots 5. Axisymmetry AXIC Defines number of harmonics for conical shell problem AXIF Defines parameters for axisymmetric fluid analysis FLSYM Axisymmetric symmetry control AXSL0T Defines parameters for acoustic cavity analysis RINGAX Ring location for conical shell problem SECTAX Sector location on RINGAX P0INTAX Point location on RINGAX 6. Cyclic Symmetry CYAX Defines grid points on axis of symmetry CYJOIN Defines boundary points of a segment 7. Superelement Analysis GRID Defines interior points for a superelement SESET Defines interior points for a superelement CSUPEXT Defines exterior points for a superelement 2.4-7 (1/15/81)

NASTRAN DATA DECK CSUPER Defines grid point connections for identical or mirror-image superelement or superelements from an external source SEQSEP Used with CSUPER card to define internal sequence of degrees of freedom for mirror-image or identical superelements RELEASE Defines released degrees of freedoms for superelement exterior grid points SEELT Changes superelement membership of elements 2.4.2.2 Elements 1. Elastic Line Elements CBAR Connection definition for prismatic beam;- - PBAR Property definition for CBAR BAROR Default for orientation and property for CBAR BEAM0R Default for orientation and property for CBEAM - CBEAM Connection definition for general beam element > PBEAM Property definition for CBEAM CBEND Connection definition for curved beam PBEND Propery definition for CBEND C0NROD Connection and property definition for rod with axial and torsional stiffness -^ CROD Connection definition for rod with axial and torsional stiffness - PR0D Property definition for CR0D CTUBE Connection definition for tube with axial and torsional stiffness PTUBE Property definition for CTUBE 2. Elastic Surface Elements CTRIA3 Connection definition for an isoparametric triangle with bending and membrane stiffness PSHELL Property definition for CTRIA3, CTRIA6, CQUAD4, and CQUAD8 - CTRIA6 Connection definition for curved triangular shell element with six grid points — > CQUAD4 Connection definition for an isoparametric quadrilateral with bending and membrane stiffness -- CQUAD8 Connection definition for curved quadrilateral shell element with eight grid points CSHEAR Connection definition for shear panel 2.4-8 (1/15/81)

BULK DATA DECK PSHEAR Property definition for CSHEAR PC0MP Property definition for composite material laminate _- CQDMEM Connection definition for quadrilateral membrane composed of constant strain triangles _- PQDMEM Property definition for CQDMEM 3. Elastic Solid Elements CTETRA Connection definition for constant strain tetrahedron CPENTA Connection definition for five-sided solid element with from six to fifteen grid points CHEXAi Defines six-sided solid element composed of tetrahedra 4 CHEXA Connection definition for six-sided solid element with from eight to twenty grid points CHEX20 Connection definition for isoparametric hexahedron with a maximum of twenty nodes PHEX Property definition for CHEX20 PS0LID Property definition for CHEXA and CPENTA 4. Elastic Scalar Elements — * CELASI Connection definition for scalar spring, also property definition for 1=2 or 4 PELAS Property definition for CELAS1 or CELAS3 GENEL Element definition in terms of stiffness coefficients or deflection influence coefficients 5. Axisymmetric Elements CC0NEAX Connection definition for conical shell PC0NEAX Property definition for CC0NEAX CTRIARG Connection and property definition for constant strain triangular ring CTRAPRG Connection and property definition for trapezoidal ring CTRIAX6 Connection and property definition for linear strain triangular ring 6. Rigid Elements - RBAR Defines rigid bar with six degrees of freedom at each end RBE1 Defines rigid body connected to an arbitrary number of grid points RBE2 Defines rigid body with independent degrees of freedom at a single grid point RBE3 Defines the motion at a "reference" grid point as the weighted average of the motions at a set of other grid points RROD Defines pin-ended rigid rod 2.4-9 (1/15/81)

NASTRAN DATA DECK RSPLINE Defines multipoint constraints for interpolation of displacements at grid points RTRPLT -Defines rigid triangular plate 7. Mass Elements CMASSi Connection definition for scalar mass, also property definition for i = 2 or 4 PMASS Property definition for CMASS1 or CMASS3 C0NM1 Defines 6x6 mass matrix at a grid point C0NM2 Defines concentrated mass at a grid point 8. Damping Elements CDAMPi Connection definition for a scalar damper, also property definition for i = 2 or 4 PDAMP Property definition for CDAMP1 or CDAMP3 CVISC Connection definition for viscous damper element PVISC Property definition for CVISC 9. Fluid Elements CAXIFi Defines axisymmetric fluid element for acoustic cavity analysis CFLUIDi Defines axisymmetric fluid element CSLOTi Defines slot element for acoustic cavity analysis ELIST Defines wetted side of structural elements 10. Heat Transfer Elements CHBDY Connection definition for boundary element PHBDY Property definition for CHBDY CFTUBE Connection definition for FTUBE element PFTUBE Property definition for FTUBE element The following elastic elements may also be used as heat conduction elements: Linear BAR, R0D, C0NR0D, TUBE, BEAM, BEND Membrane TRIA3, TRIA6, QUAD4, QDMEM, QUAD8 Axisymmetric TRIARG, TRAPRG Solid TETRA, HEXA, PENTA, HEXA1, HEXA2 11. Dummy Elements ADUMi Attributes for CDUMi element CDUMi Connection definitions for user-defined element 2.4-10 (1/15/81)

BULK DATA DECK PDUMi Property definition for CDUMi PL0TEL Defines plot element 12. Nonlinear Elements CGAP Defines gap or frictional element PGAP Property definition for CGAP 2.4.2.3 Materials 1. Isotropic MAT1 Defines elastic material properties MAT4 Defines thermal material properties 2. Anisotropic MAT2 Defines anisotropic material properties for two-dimensional element MAT3 Defines orthotropic material properties for TRIARG, TRAPRG and TRIAX6 elements MATS Defines anisotropic thermal material properties MAT8 Defines orthotropic material properties MAT9 Defines anisotropic material properties for isoparametric solid elements 3. Temperature Dependent MATT1-5 Table references for temperature-dependent MAT1 - MATS materials MATT9 Table references for temperature-dependent MAT9 materials TEMP Defines temperature at grid points TEMPD Specifies default temperature at grid points TEMPPi Defines temperature field for surface elements TEMPRB Defines temperature field for line elements TEMPAX Defines temperature field for conical shell problem TABLEMi Tabular functions for generating temperature-dependent material properties 4. Stress Dependent MATS1 Table references for stress-dependent MAT1 materials TABLES1 Defines tabular stress-strain function 2.4-11 (1/15/81)

NASTRAN DATA DECK 5. Fluid AXIF Includes default values for mass density and bulk modulus AXSL0T Includes default values for mass density and bulk modulus BDYLIST Includes mass density at boundary CFLUIDi Includes mass density and.bulk modulus CSL0Ti Includes mass density and bulk modulus FSLIST Includes mass density at free surface SLBDY Includes mass density at interface between fluid and radial slots MFLUID Defines properties for incompressible fluid 2.4.2.4 Constraints and Partitioning 1. Single-Point Constraints SPC Defines single-point constraints and enforced displacements -- SPC1 Defines single-point constraints SPCADD Defines a union of single-point constraint sets on SPC of SPC1 cards SPCAX Defines single-point constraints for conical shell problems -- GRID Includes single-point constraint definition GRID8 Includes single-point constraint definition GRDSET Includes default for single-point constraints FLSYM Symmetry control for boundary in axisymmetric fluid problem 2. Multipoint Constraints -,>MPC Defines a linear relationship for two or more degrees of freedom MPCADD Defines a union of multipoint constraint sets on MPC cards MPCAX Defines multipoint constraints for conical shell problems P0INTAX Program generates MPC equations for point on conical shell -RBAR Program generates MPC equations for rigid bar RBE1-3 Program generates MPC equations for RBE1, RBE2, RBE3 RR0D Program generates MPC equations for rigid rod RSPLINE Program generates MPC'equations for spline element RTRPLT Program generates MPC equations for rigid triangular plate 2.4-12 (1/15/81)

BULK DATA DECK 3. Partitioning ASET Defines independent degrees of freedom in analysis set ASET1 Defines independent degrees of freedom in analysis set 0MIT Defines dependent degrees of freedom in omitted set 0MIT1 Defines dependent degrees of freedom in omitted set 0MITAX Defines omitted degrees of freedom for conical shell problems GRID Defines interior points for a superelement SESET Defines interior points for a superelement CSUPEXT Defines exterior points for a superelement 4. Free Body Supports SUP0RT Defines coordinates for determinate reactions SUPAX Defines determinate reactions for conical shell problems CYSUP Defines determinate reactions for cyclic symmetry problems 5. Component Mode Boundary Conditions BSET Defines fixed boundary points for residual structure BSET1 Defines fixed boundary points for residual structure CSET Defines free boundary points for residual structure CSET1 Defines free boundary points for residual structure QSET Defines generalized coordinates for residual structure QSET1 Defines generalized coordinates for residual structure SEBSET Defines fixed boundary points for superelement SEBSET1 Defines fixed boundary points for superelement SECSET Defines free boundary points for superelement SECSET1 Defines free boundary points for superelement SEQSET Defines generalized coordinates for superelement SEQSET1 Defines generalized coordinates for superelement SESUP Defines free boundary points used as references for rigid body modes for superel ement 6. User Sets DEFUSET Defines names for user sets SEUSET Defines superelement user sets 2.4-13 (1/15/81)

NASTRAN DATA DECK SEUSET1 Defines superelement user sets USET User set definition USET1 User set definition 2.4.2.5 Loads 1. Static Loads DEFORM Enforced axial deformation for line elements -v FORCE Defines concentrated load at grid point -4 FORCEi Defines concentrated load at grid point - GRAV Defines gravity load vector MOMENT Defines moment at grid point MOMENTi Defines moment at grid point PLEAD Defines pressure load on an area -- PLOAD1 Defines distributed and concentrated loads on BAR, BEAM, and BEND elements — 4>'PL0AD2 Defines pressure loads on surface elements "-*) PL0AD3 Defines pressure loads on surfaces of HEX20 element - PL0AD4 Defines pressure loads on surfaces of HEXA,. PENTA, TRIA3, QUAD4 and QUAD8 elements -- PLOADX Defines pressure load on TRIAX6 RFORCE Defines load due to centrifugal force field SPCD Defines value for enforced displacement SLOAD Defines load on scalar point -- LOAD Linear combination of static load sets TEMP Defines temperature at grid points TEMPD Specifies default temperature at grid points TEMPPi Defines temperature field for surface elements TEMPRB Defines temperature field for line elements F0RCEAX Defines concentrated load for conical shell problem MOMAX Defines moment for conical shell problem PRESAX Defines pressure load for conical shell problem 2.4-14 (1/15/81)

BULK DATA DECK TEMPAX Defines temperature field for conical shell problem L0ADCYH Defines harmonic coefficients of static load in cyclic symmetry analysis L0ADCYN Defines physical load input in cyclic symmetry analysis L0ADCYT Defines load input for AXI type cyclic symmetry problems 2. Dynamic Loads DAREA Dynamic load scale factor definition DELAY Dynamic load function time delay definition DPHASE Dynamic load function phase lead definition RL0ADI Frequency dependent load definition TL0ADi Time dependent load definition N0LINi Nonlinear transient load definition TABLEDi Tabular functions for generating dynamic loads DL0AD Linear combination of dynamic load sets L0ADCYH Defines harmonic coefficients of dynamic load in cyclic symmetry analysis L0ADCYN Defines physical load input in cyclic symmetry analysis LSEQ Defines static load sets for dynamic analysis 3. Heat Transfer NFTUBE Defines nonlinear transient load for heat convection N0LINS Defines nonlinear transient radiant heat transfer QHBDY Defines uniform heat flux into a set of grid points QBDY1 Defines uniform heat flux into HBDY element QBDY2 Defines grid point heat flux into HBDY element QVECT Defines thermal vector flux from distant source into HBDY element QV0L Defines internal heat generation RADLST Defines list of radiation areas RADMTX Defines matrix of radiation exchange coefficients TEMP Defines temperature at grid points TEMPO Specifies default temperature at grid points VIEW Defines shading and submesh sizes for radiation exchange calculations 2.4-14a (1/15/81)

NASTRAN DATA DECK 2.4.2.6 Problem Control I 1. Buckling Analysis EIGB Defines eigenvalue extraction data 2. Eigenvalue Analys s EIGR Defines real eigenvalue extraction data EIGC Defines complex eigenvalue extraction data EIGP Defines poles for complex eigenvalue analysis TABDMP1 Defines structural damping as a tabular function of frequency 3. Cyclic Symmetry CYSYM Defines cyclic symmetry parameters 4. Dynamics DYNRED Defines data for dynamic reduction 5. Frequency Response FREQ List of frequencies for problem solution FREQi Defines a set of frequencies for problem solution TABDMP1 Defines structural damping as a tabular function of frequency 6. Random Response RANDPS Defines load set power spectral density factors RANDT1 Defines time lag constants for use in autocorrelation function computation TABRND1 Defines power spectral density as a tabular function of frequency 7. Transient Response TIC Specifies initial values for displacement and velocity TSTEP Specifies time step intervals for solution and output 8. Nonlinear Analysis NLPARM Defines parameters for nonlinear analysis 2.4.2.7 Miscellaneous 1. Comments -—, $ For inserting comment cards in unsorted echo of Bulk Data Deck 2.4-15 (1/15/81)

NASTRAN DATA DECK 2. Delete / For removing cards from the Bulk Data Deck on restart 3. Parameters' — PARAM Specifies values for parameters used in DMAP sequences or rigid formats 4. Direct Matrix Input C0NM1 Define 6x6 mass matrix at a grid point DMI.User-defined direct matrix input DMIG Direct matrix input related to grid points DMIAX Direct matrix input for fluid analysis TF Defines dynamic transfer function 5. Tabular Input DTI User-defined direct table input TABDMP1 Defines structural damping as a tabular function of frequency TABLEDI Tabular functions for generating dynamic loads TABLEMi Tabular functions for generating temperature-dependent material properties TABLES1 Defines tabular stress-strain function TABRND1 Defines power spectral density as a tabular function of frequency 6. Output Control CBARA0 Defines auxiliary output points along axis of BAR element FREEPT Surface point location for data recovery in hydroelastic problems PL0TEL Defines dummy element for plotting P0INTAX Defines point on conical shell ring (RINGAX) for recovery of displacements PRESPT Pressure point location for data recovery in hydroelastic problems SET1 Defines a set of grid points TSTEP Specifies time step intervals for data recovery in transient response 7. Matrix Assembly CNGRNT Specifies congruent elements 2.4.2.8 Bulk Data Generator (MSGMESH) See the MSGMESH Analyst's Guide for the MSGMESH bulk data cards. 2.4-16 (1/15/81)

BULK DATA DECK 1. Grid Point Generation ACTIVE Replaces (activates) a group of grid points in a grid' point field by grid points in another grid point field DELETE Deletes a set of grid points within a grid field DIST0RT- Modifies grid point location entries on generated GRID card EGRID Defines grid points for use by MSGMESH program EQUIY Equivalences (replaces) grid point identifications or activates automatic equivalence option GRIDG Defines a field of grid points and generates GRID cards GRIDM0D Modifies entries on generated GRID cards that do not affect grid point locations GRIDU Defines a user-defined grid point field INSECT Generates EGRID cards on the intersection of surfaces LIST Describes an unequal spacing of grid points in a GRIDG field MAREA Defines an active grid point area for automatic equivalencing MLINE Defines an active grid point line for automatic equivalencing MP0INT Defines a common grid point for automatic equivalencing RENAME Replaces (renumbers) a group of grid points in a grid point field by grid points in another or the same grid point field 2. Element Generation CBARG Generates CBAR cards CGEN Generates element connection cards PGEN Specifies property identificaiton numbers for generated connection cards PL0TE Generates PL0TEL cards which form the outlines of element sets PL0TG Generates PL0TEL cards which form the outlines of grid point fields RSPLING Generates a logical RSPLINE card that interconnects grid fields of different mesh sizes TRIC0N Generates triangular elements that fill the space between two fields of different mesh sizes 3. Temperature Definition ETEMP Defines temperatures at the vertices of a grid point field TEMPG Generates TEMP cards for a grid point field 2.4-16a (2/2/81)

NASTRAN DATA DECK 4. Load and Constraint Definition PL0ADG Generates PL0AD or PL0AD3 cards for a grid point field SPCG Generates SPC1 cards for a subset of grid points within a grid point field 5. MSGMESH Control MESH0FF Suppresses the generation of bulk data cards MESH0N Overrides a previous MESH0FF card MESH0PT Controls MSGMESH output options and selects location numbers for HEX fields PL0T0PT Selects MSGMESH plot options 2.4.2.9 Variance Analysis VARIAN Deviation values for variance analysis 1PARM Parameter values for variance analysis 1VARY Functional relationships for variance analysis 2.4.2.10 Aeroelastic Analysis See the Aeroelastic Supplement for the aeroelastic bulk data cards. 1. Aerodynamic Elements AEFACT Specifies lists of real numbers CAER01 Connection definition for aerodynamic panel CAER02 Connection definition for aerodynamic body CAER04 Connection definition for aerodynamic macro-strip PAER01 Defines associated bodies for panels PAER02 Property definition for aerodynamic bodies PAER04 Property definition for strip elements 2. Aerodynamic Data AERO Specifies aerodynamic parameters FLFACT Specifies aerodynamic data for flutter analysis FLUTTER Specifies aerodynamic data for flutter analysis 2.4-16b (2/2/81 )

BULK DATA DECK Input Data Card $ Comment Description: For user convenience in inserting commentary material into the unsorted echo of his input Bulk Data Deck. The $ card is otherwise ignored by the program. These cards will not appear in a sorted echo nor will they exist on the New Problem Tape. Format and Example: 1 2 3 4 5 6 7 8 9 10 |$ followed by any'egitimats characters in ca-d column; 2-80 $ TI THIS IS A REMARK (*,'$$)-T/ 2.4-17

BULK DATA DECK Input Data Card CBAR Simple Beam Element Connection Description: Defines a simple beam element (BAR) of the structural model. Format and Example: 1 2 3 4 5 6 7 8 9 10 CBAR EID PID. GA GB X1,GO X2 X3 CBAR 2 39 7 3' 13 __23 PA PB W1A W2A W3A W1 B W2B W3B _ +23 513 Field Contents EID Unique element identification number (Integer > 0). PID Identification number of a PBAR property card (Default is EID unless BAR0R card has nonzero entry in field 3)(Integer > 0 or blank *). GA,GB Grid point identification numbers of connection points (Integer > O; GA f GB). X1,X2,X3 Components of vector v, at end A, (Figure l(a) in Section 1.3) measured at end A, parallel to the components of the displacement coordinate system for GA, to determine (with the vector from end A to end B) the orientation of the element coordinate system for the BAR element (Real, 0 or blank*). GO Grid point identification number to optionally supply Xl, X2, X3 (Integer > 0 or blank*). Direction of orientation vector is GA to GO. PA,PB Pin flags for bar ends A and B, respectively (up to 5 of the unique digits 1 - 6 anywhere in the field with no imbedded blanks; Integer > 0). Used to remove connections between the grid point and selected degrees of freedom of the bar. The degrees of freedom are defined in the element's coordinate system (see Figure l(a), Section 1.3). The bar must have stiffness associated with the pin flag. For example, if PA=4 is specified, the PBAR card must have a value for J, the torsional stiffness. W1A,W2A,W3A Components of offset vectors wa and wb, respectively (see Figure l(a), Section 1.3), W1B,W2B,W3B in displacement coordinate systems at points GA and GB, respectively (Real or blank). *See the BAR0R card for default options for fields 3 and 6 - 8. Remarks: 1. Element identification numbers must be unique with respect to all other element identification numbers. 2. For an explanation of BAR element geometry, see Section 1.3.2. 3. If there are no pin flags or offsets, the continuation card may be omitted. 2.4-39 (1/30/81)

NASTRAN DATA DECK CBAR (Cont.) 4. The old CBAR card used field 9 for a flag, F, which was used to specify the nature of fields 6 - 8 as follows: FIELD 6 7 8 F= Xl X2 X3 F=2 GO Bl ank Bl ank or 0 or 0 F=blank Provided by BAR0R card. This data item is no longer required but may continue to be used if desired (See Remark 5). If F=l in field 9, a zero (0) in field 6, 7, or 8 will override entries on the BAR0R card, but a blank will not. 5. For the case where-field 9 is blank and not provided by the BAR0R card, if Xl,GO is integer, then GO is used; if Xl,GO is blank or real, then Xl, X2, X3 is used. 2.4-40 (5-15-80)

BULK DATA DECK Input Data Card CBEAM Beam Element Connection.Description: Defines a beam element (BEAM) of the structural model. Format and Example: 1 2 3 4 5 6 7 8 9 10 CBEAM EID PID GA GB X1,GO X2 X3 -. _ _ CBEAM 2 39 7 3 13., __. 123_ _ P PB WlA BW2A W3A WlB W2B W3B +23 513 3.0 | |.. 1234 SA SB +34 8 5i Field Contents EID Unique element identification number (Integer > 0) PID Identification number of PBEAM property card (default is EID)(Integer > 0 or blank*) GA.GB Grid point identification numbers of connection points (Integer > O; GA A GB) X1,X2,X3 Components of vector v, at end A (shown in the following figure), measured at the offset point for end A, parallel to the components of the displacement coordinate system for GA, to determine (with the vector from offset end A to offset end B) the orientation of the element coordinate system for the beam element (Real or blank;* see Remark 3) GO Grid point identification number to optionally supply X1, X2, X3 (Integer > 0 or blank;* see Remark 3) PA,PB Pin flags for beam ends A and B respectively (Up to five of the unique digits 1 - 6 with no imbedded blanks; integer > 0). Used to move connections between the grid point and selected degrees of freedom of the beam. The degrees of freedom are defined in the element's coordinate system and the pin flags are applied at the offset ends of the beam (see the following figure). The beam must have stiffness associated with the pin flag. For example, if PA-4, the PBEAM card must have a nonzero value for J, the torsional stiffness. W1A,W2A,W3A Components of offset vectors, measured in the displacement coordinate systems at W1B,W2B,W3B grid points A and B, from the grid points to the end points of the axis of shear center (Real or blank) SA,SB Scalar or grid point identification numbers for the ends A and B, respectively. The degrees of freedom at these points are the warping variables, d9/dx (Integers > 0 or both blank). *See the BEAM0R card for default options for fields 3 and 6 - 8. (Continued) 2.4-42c (5-15-80)

NASTRAN DATA DECK CBEAM (Cont.) Remarks: 1. Element identification numbers must be unique with respect to all other element identification numbers. 2. For an explanation of beam element geometry, see Section 1.3.2. 3. If X1,GO is integer, then GO is used. If Xl,GO is blank or real, then X1,X2,X3 is used. 4. GO f GA or GB. 5. If there are no pin flags or offsets or warping variables, both continuation cards may be omitted. 6. The first continuation card must be included, even if all fields are blank, if the second continuation card is used. 7. If the second continuation card is omitted, torsional stiffness due to warping of the cross-section will not be considered. 8. If warping is allowed (SA and SB > 0), then SA and S3 must be defined with SP0INT or GRID cards. If GRID cards are used, the warping degree of freedom is attached to the first (T1) component. elem /wfnaft ^ el 7a mas Gd..pn.t. G A I -,' y, ^ (, op (0100,O) V( oWA ~ ~ ~ ~ ~ ~ ~. offfsete Element Coordinate System Grid point Gr 2.4-42d (5-15-80)

BULK DATA DECK Input Data Card CELAS2 Scalar Spring Property and Connection Description: Defines a scalar spring element of the structural model without reference to a property card. Format and Example: 1 2 3 4 5 6 7 8 9 10 CELAS2 EID K G1 C1 G2 C2 GE S CELAS'2 28 6.2+3 32 19 4; Field Contents EID Unique element identification number (Integer > 0) K The value of the scalar spring(Real) G1, G2 Geometric grid point identification number (Integer > 0) C1, C2 Component number (6 > Integer > 0) GE Damping coefficient (Real) S Stress coefficient (Real) Remarks: 1. Scalar points may be used for G1 and/or G2 in which case the corresponding Ct and/or C2 must be zero or blank. Zero or blank may be used to indicate a grounded terminal G1 or G2 with a corresponding blank or zero C1 or C2. A grounded terminal is a point whose displacement is constrained to zero. If only scalar points and/or ground are involved, it is more efficient to use the CELAS4 card. 2. Element identification numbers must be unique with respect to all other element identification numbers. 3. This single card completely defines the element since no material or geometric properties are required. 4. The two connection points (G1, C1) and (G2, C2), must be distinct. 5. For a discussion of the scalar elements, see Section 5.6 of the Theoretical Manual. 2.4-57

BULK DATA DECK Input Data Card CHEXA Six-sided Solid Element with from Eight to Twenty Grid Points Description: Defines the connections of the HEXA solid element. Format and Example: 1 2 3 4 5 6 7 8 9 10 CHEXA EID |PID G1 IG2 | G3 G4 1G5 G6 1 CHEXA 71 4 3 14 5 6 7 8 _ A8C +6C 9'1 0 I 0 0I 30 31 53 54 DEF G15 G16 G17 G18 G19 G20 +EF J55 56 57 58 59 160 Field Contents EID Element identification number (Integer > 0) PID Identification number of a PSOLID property card (Integer > 0) G1,...,G20 Grid point identification numbers of connection points (Integer > 0 or blank) GG8 G7/ G G1 G2~G Remarks. Element identification numbers must be unique with respect to all other elent Remarks: 1. Element identification numbers must be unique with respect to all other element identification numbers. 2. Grid points G1,.., G4 must be given in consecutive order about one quadrilateral face. G5,..., G8 must be on the opposite face with GS opposite G1, G6 opposite G2, etc. 3. The edge points, G9 to G20, are optional. Any or all of them may be deleted. If the ID of any edge connection point is left blank or set to zero (as for G9 and GO1 in the example), the equations of the element are adjusted to give correct results for the reduced.number of connections. Corner grid points ca.1 not be deleted. The element is an isoparametric element (with shear correction) in all cases. 4. Components of stress are output in the material coordinate system. 5. The second continuation card is not reauired. (Continued) 2.4-69.(5-15-80)

NASTRAN DATA DECK CHEXA (Cont.) 6. The element coordinate system for the HEXA element is defined in terms of the three face-to-face lines. The axes are chosen by the following rules: Taxis: Longest line joining centroids of opposite faces. XY plane: Plane containing longest and next longest lines joining centroids of opposite faces. XYZ axes: Permute XYZ axes to make X-axis approximately parallel to edge 1-2 and Yaxis approximately parallel to edge 1-4. In the example shown, X = Y 7 Z Z G7 348"CentCenroid id G3 4G.7. It is recommended that the edge points be located within the middle third of the edge. If the edge point is located at the quarter point, the calculated stresses will be meaningless. 2.4-70 (5-15-80)

BULK DATA DECK Input Data Card C0RD2C Cylindrical Coordinate System Definition, Form 2 Description: Defines a cylindrical coordinate system by reference to the coordinates of three points. The first point defines the origin. The second point defines the direction of the z-axis. The third lies in the plane of the azimuthal origin. The reference coordinate system must be independently defined. z U z A Z r Format and Example: 1 2 3 4 5 6 7 8 9 10 CORD2C C ID RID Al A2 A3 B1 B2 B3 ABC 0RD2C 3 17 -2.9' 1.0 0.0 3.6 0.0.0 123 BC C1 C2 C3 23 1 5.2 1.0 -2.9 f____ I I _ I___ Field Contents CID Coordinate system identification nunber (Integer > 0) RID Reference to a coordinate system which is defined independently of new coordinate system (Integer > 0 or blank) A1,A2,A3 B1,B2,B3 Coordinates of three points in coordinate system defined in field 3 (Real) C1,C2,C3 (continued) 2.4-97

NASTRA DATA DECK C3RD2C (cont.) Remarks: 1. Continuation card must be present. 2. The- three points (A1, A2, A3), (B1, B2, B3), (C1, CZ, C3) must be unique and nonco31 near. Nncollinearity is checked by the geometry processor. 3- Coordinate system identification numbers on all C0RD1R, CROlC, CWRD1S, C0RD2R, CORDZC, and CORDZS cards must all be unique. 4. An RID of zero references the basic coordinate system. 5. The location of a grid point (P in the sketch) in this coordinate system is given by (R, e, Z) where 8 is measured in degrees. 5. The displacenent coordinate directions at P are dependent- on the tocation of P as shcwn above by (Ur u, u z). 7.. Points on the z-axis may not have their displacment direction defined irr this coordinate system since an ambiguity results. 2.4-98

BULK DATA DECK Input Data Card CPENTA Five-sided Solid Element with from 6 to 15 grid points Description: Defines the connections of the CPENTA element. Format and Example: 1 2 3 4 5 6 7 8 9 10 CPENTA EID PID G1 G2 G3 |G4 I G5 G6 CPENTA 112 1 2 T 3 15 14 14 I 103 115 ABC | G7 G8 G9 G10 I Gl G12 G13 G14 +8C 5 116 8 1 120 125 8CD G15 [CD 130 G5 G Field Contents EID Element identification number (Integer > 0) PID Identification number of a PSOLID property card (Integer > 0). G1 - G15 Identification numbers of connected grid points (Integer > 0 or blank) G6 G15 G14 G131 G5 G4 Gl3M G5 ~/ I ~*Gl2 GlO 10_ 4 G11 / GG9r G3 "G1 8 \ Gl G7 G2Re Remarks: 1. Element ID numbers must be unique with respect to all other element ID numbers. 2. The topology of the diagram must be preserved, i.e., G1, G2, G3 define a triangular face, G1, GlO and G4 are on the same edge, etc. 3. The edge grid points, G7 to G15, are optional. Any or all of them may be deleted. In the example shown, GlO, G11 and G12 have been deleted. The continuation cards are not required if all edge grid points are deleted. 4. Components of stress are output in the material coordinate system. 5. The element coordinate system is defined as follows: (Ccntinucd) 2.4-102a (5-15-80)

NASTRAN DATA DECK CPENTA (Cont.) G6 nes joning G and G 6 G4 G69 x G 3 / GI where the x-axis joins the midpoints of straight lines joining points G1-G4 and G2G5, and the z-axis is normal to a plane passing through the midpoints of straight lines joining G1-G4, G2-G5 and G3-G6. 6. It is recommended that the edge points be located within the middle third of the edge. If the edge point is located at the quarter point, the calculated stresses will be meaningless. 2.4-102b (5-15-80)

BULK DATA DECK Input Data Card CQDMEM Quadilateral Element Connection Description: Defines a quadrilateral membrane element (QDMEM) of the structural model consisting of four overlapping triangular membrane elements. Format and Example: 1 2 3 4 5 6 7 8 9 10 CQDMEM EID PID G1 G2 G3 G4 TH CQDMEM 72 13 3 14 15 16 29.2 Field Contents EID Element identification number (Integer > 0) PID Identification number of a PQDMEM property card (Default is EID) (Integer > 0) G1,G2,G3,G4 Grid point identification numbers of connection points (Integer > 0; G1 # G2 f G3 ) G4) TH Material property orientation angle in degrees (Real). The sketch below gives the sign convention for TH. G3 G4 GiL- --- - 2 Remarks: 1. Element identification numbers must be unique with respect to all other element identification numbers. 2. Grid points Gi thru G4 must be ordered consecutively around the perimter of the el ement. 3. All interior angles must be less than 180~. 4. This element is for heat transfer. The QUAD4 element should be used for structural problems. 5. This element may also be used for differential stiffness applications with membrane stiffness only. 2.4-102c (5-15-80)

BULK DATA DECK Input Data Card CQUAD4 Quadrilateral Element Connection Description: Defines a quadrilateral plate element (QUAD4) of the structural model. This is an isoparametric membrane-bending element. Format and Example: 1 2 3 4 5 6 7 8 9 10 CQUAD4 EID PID Gi G2 G3 G4 e CQUAD4 111 203 31 74 75 32 2.6 __BC_ i 2 T3 T4 Lm! ~ _ - _. _ +BC 1.77 2.04 2.09 1.80 Field Contents EID Element identification number (Unique integer > 0). PID Identification number of a PQUAD4 property card (Integer > 0 or blank, default is EID). G1,G2,G3,G4 Grid point identification numbers of connection points (Integers > 0, all unique). e Material property orientation specification (Real or blank; or 0 < Integer < 1,000,000). If Real or blank, specifies the material property orientation angle in degrees. If Integer, the orientation of the material x-axis is along the projection onto the plane of the element of the x-axis of the coordinate system specified by the integer value. The sketch below gives the sign convention for e. T1,T2,T3,T4 Membrane thickness of element at grid points G1 through G4 (Real or blank, see PQUAD4 for default). Yelement G3 Remarks: 1. Element identification numbers must be unique with respect to all other element el ement. 3. All the interior angles must be less than 180~ ( Continued) 2.4-112a (8/4/80)

NASTRAN DATA DECK CQUA04 (Cont.) 4. The continuation card is optional. If it is not supplied, then T1 through T4 will be set equal to the value of T on the PSHELL data card. 5. Stresses are output in the element coordinate system. 6. The PQUAD4 card may be used in place of the PSHELL card. If a single PQUAD4 card is in the Bulk Data Deck, there must not be any references on CQUAD4 cards to PSHELL cards. The PSHELL card is recommended. 2.4-112b (8/4/80)

BULK DATA DECK Input Data Card CQUAD8 Quadrilateral Element Connection Description: Defines a curved quadrilateral shell element (QUAD8) with eight grid points. Format and Example: 1 2 3 4 5 6 7 8 9 10 CQUAD8 EID PID G1 G2 G3 G4 G5 G6 CQUAD8 207 3 31 33 73 71 32 51 ABC i Fi G7 G8 I T2 I T3 T4 +BC [ 53 1 72 0.125 0.025 0.030.025 30. Field: Contents EID Element identification number (Integer > 0) PID Identification number of a PSHELL property card (Integer > 0) G1,G2,G3,G4 Identification numbers of connected corner grid points (Unique integers > 0). Required data for all four grid points. G5,G6,G7,G8 Identification numbers of connected edge grid points (Integer > 0 or blank). Optional data for any or all four grid points. TI,T2,T3,T4 Membrane thickness of element at corner grid points. See Remark 4 for default. a Material property orientation specification (Real or blank; or 0 < Integer < 1,000,000). If Real or blank, specifies the material property orientation angle in degrees. If Integer, the orientation of the material x-axis is along the projection on to the plane of the element of the x-axis of the coordinate system specified by the integer value. The sketch below gives the sign convention for 9. Ye lem ~/-^ 1 f G7 G4 =-O=a. 6 Ym G3 G( ~~~8~6 Xelem (Continued) 2.4-112C (8/4/80)

NASTRAN DATA DECK CQUAD8 (Cont.) Remarks: 1. Element identification numbers must be unique with respect to all other element ID's of any kind. 2. Grid points G1 to G8 must be numbered as shown. 3. The material property orientation angle, e, is defined locally at each interior integration point as the angle made by x-axis of the material coordinate system with a line, nrconst. If the shape of the element is a parallelogram and if edge points are located at midpoints of the sides, the lines n=const. are parallel to side 1-2. 4. T1, T2, T3 and T4 are optional. If not supplied they will be set equal to the value of T on the PSHELL card. 5. It is recommended that the edge points be located within the middle third of the edge. If the edge point is located at the quarter point, the program may fail with a divide by zero or the calculated stresses will be meaningless. 2.4-112d (5-15-80)

BULK DATA DECK Input Data Card CR0D Rod Element Connection DescriDtion: Defines a tension-compression-torsion element (ROD) of the structural model. Format and Example: 1 2 3 * 4 5 6 7 8 9 10 CR0D EID PID G1 G2 CR0D 12 13 21 23. Field Contents EID Element identification number (Integer > 0) PID Identification number of a PROD property card (Defauit is EID) (Integer > 0) G1, G2 Grid point identification numbers of connection points (Integer > O; G1 f G2) Remarks: i. Element identification numbers must be unique with respect to all other element identification numbers. 2. See C0NR0D for alternative method of rod definition. 3. Only one RVD element may be defined on a single card. 2.4-113 (5-15-79)

BULK DATA DECK Input Data Card CTRAPRG Trapezoidal Ring Element Connection Description: Defines an axisymmetric trapezoidal cross-section ring element (TRAPRG) of the structural model without reference to a property card. Format and Example: 1 2 3 4 5 6 7 8 9 10 CTRAPRG EID G 2 G1 3 G4 TH MID_ CTRAPRG 72 13 14 15 16 29.2 13 -- Field Contents EID Element identification number (Integer > 0). G1,G2,G3,G4 Grid point identification number of connection points (Integers > 0; G1 # G2 # G3 # G4). TH Material property orientation angle in degrees (Real) - The sketch below gives the sign convention for TH. MID Material property identification number (Integer > 0). z Material Axis G3 Axis G4 of TH Symmetry Gl -G2-L G - - x Remarks: 1. Element identification numbers must be unique with respect to all other element identification numbers. 2. The four grid points must lie in the x-z plane of the basic coordinate system and to the right of the axis of symmetry (the z-axis). Either (X1 thru X4 > 0.) or (X l X4 - 0 with X2, X3 > 0.). 3. Grid points G1, G2, G3 and G4 must be ordered counterclockwise around the perimeter of the element as in the above sketch. 4. The line connecting grid points GI and G2 and the line connecting grid points G3 and G4 must both be parallel to the x-axis. 5. All Interior angles must be less than 180~. 6. For structural problems, the material property identification number must reference only a MAT1 or MAT3 card. 7. For heat transfer problems, the material property identification number must reference only a MAT4 or MATS card. 8. See Section 5.11 of the Theoretical Manual for mathematical details. 9. For a core element, X1 - X4 - 0, the x-component should be constrained at G1 and G4 to avoid matrix singularity. 2.4-123 (5-15-79)

BULK DATA DECK Input Data Card CTRIA6 Triangular Element Connection Description: Defines a curved triangular shell element (TRIA6) with six grid points. Format and Example: 1 2 3 4 5 6 7 8 9 10.._.___,.,-......... CTRIA6 EID PID G1 G2 G3 G4 G5 G6 CTRIA6 302 3 31 33 71 32 51 52 ABC +BC 45..020.025.025 Field Contents EID Element Identification number (Integer > 0). PID Identification number of PSHELL property card (Integer > 0). G1 to G3 Identification numbers of connected corner grid points (Unique integers > 0). G4- to G6 Identification number of connected edge grid points (Integer > 0 or blank). Optional data for any or all three points. a Material property orientation specification (Real or blank; or 0 < Integer < 1,000,000). If Real or blank, specifies the material property orientation angle in degrees. If Integer, the orientation of the material x-axis is along the projection on to the plane of the element of the x-axis of the coordinate system specified by the integer value. The sketch below gives the sign convention for 9. T1,T1,T3 Membrane thickness of element at corner grid points. See Remark 4 for default. G3 G G G2 g 2.4-1 (8 ) const

NASTRAN DATA DECK CTRIA6 (Cont.) Remarks: 1. Element identification numbers must be unique with respect to all other element ID's of any kind. 2. Grid points G1 to G6 must be numbered as shown. 3. The material property orientation angle, 9, is defined locally at each interior integration point as the angle made by x-axis of the material coordinate system with a line, n =const. If the sides are straight and if edge grid points are located at the midpoints of the sides, the lines n=const. are parallel to side 1-2. 4. T1, T2 and T3 are optional. If not supplied they will be set equal to the value of T on the PSHELL card. 5. It is recommended that the edge points be located within the middle third of the edge. If the edge point is located at the quarter point, the program may fail with a divide by zero and the calculated stresses will be meaningless. 2.4-128d (7/21/80)

BULK DATA DECK Input Data Card CTRIAX6 Connections for TRIAX6 Element Description: Defines a linear strain axisymmetric triangular cross section ring element (TRIAX6) with mi dside grid points. Format and Example: 1 2 3 4 5 6 7 8 9 10 CTRIAX6 EID MID G1 G2 G3 G4 G5 G6 CTRIAX6 22 999 10 11 12 21 1 22 32 +C22 J TH +C22 9.0 Field Contents EID Element identification number (Integer > 0) MID Material identification number (Integer > 0) G1 thru G6 Grid point identification numbers of connected points (unique Integers > 0) TH Material orientation angle (Real, in degrees) z mZm ~ ~ ~ ~ ~ ~ rm G5 G4 G6 G3 TH G2 r 5 xBASIC Remarks: 1. The grid points must lie in the x-z plane of the basic coordinate system, with x = r > 0. The grid points must be listed consecutively beginning at a vertex and proceeding around the perimeter in either direction. 2. For structural problems, the material may be MAT1 or MAT3. 3. The continuation card is not required. 4. Material properties (if MAT3) and stresses are given in the (rm,Zm) coordinate system shown in the sketch. 5. Concentrated loads on grid circles for this element must be computed for 360~, i.e., multiply load per unit length by 2 irr. 6. G2, G4, and G6 are assumed to lie at the midpoints of the sides. The locations of these grid points (on GRID bulk data cards) are used only for global coordinate system definition, GPWG (weight generator module), centrifugal forces, and deformed structure plotting. 2.4-131(10/1/80)

BULK DATA DECK Input Data Card DMI Direct Matrix Input Description: Used to define matrix data blocks directly. Generates a matrix of the form All A12............... n EA) A21 A22...............A2n Am.. * —-.......... -..** Amn_ where the elements Aij may be real or complex numbers. Formats and Example: (The first logical card is a header card.) 1 2 3 4 5 6 7 8 9 10 DMI [ NAME " "0, FORM TIN T0UT UM N DMI QQQ 0 2 3 3 4.2 DMI NAME J 11 A(I1,J) A (I1+1,J etc. 12 DMI QQQ 1 1 1.0 2.0 3.0 4.0 3 1 A(I12,J) -etc.-________ +1 5.0 6.0__ DMI _QQQ 2 2 6.0 7.0 4 0 9.0 (etc. for each non-null column) Field Contents NAME Any NASTRAN BCD value (1-8 alphanumeric characters, the first of which must be alphabetic) which will be used in the DMAP sequence to reference the data block. FORM 1 Square matrix (not symmetric) 2 General rectangular matrix 3 Diagonal matrix (M = number of rows, N = 1) 4 Lower triangular factor 5 Upper triangular factor 6 Symmetric matrix 7 Row matrix (M = number of columns, N: 1) 8 Identity matrix (M = number of rows, N = M) TIN Type of matrix being input, as follows: 1 Real, S.P. (one field used/element) 3 Complex, S.P. (two fields used/element) 2 Real, D.P. (one field used/element)'4 Complex, D.P. (two fields used/element) (Continued) 2.4-155 (7/1/80)

NASTRAN DATA DECK DMI (Cont.) TOUT Type of matrix being output, as follows: 0 Set by precision cell 1 Real, single-precision 3 Complex, single-precision 2 Real, double-precision 4 Complex, double-precision M Number of rows in A (Integer > 0). Except for a row matrix, FORM 7. N Number of columns in A (Integer > 0). Except for FORMS 3, 7, AND 8. J Column number of A (Integer > 0). I1, 12, etc. Row number of A (Integer > 0). A(Ix,J) Values of element of A (see TIN) (Real) Remarks: 1. The user must write a DMAP (or make alterations to a solution sequence) in order to use the DMI feature since he is defining a data block. All of the rules governing the use of data blocks in DMAP sequences apply. In the example shown above, the data block QQQ is defined to be the complex, single-precision rectangular 4 x 2 matri x (1.0, 2.0) (0.0, 0.0) 1 (3.0, 4.0) (6.0, 7.0) [QQQ] (5.0, 6.0) (0.0, 0.0) (0.0, 0.0) (8.0, 9.0) 2. A limit to the number of DMI's which may be defined is set by the size of the Data Pool Dictionary. The total number of DMIs and DTIs may not exceed 79. 3. There are a number of reserved words which may not be used for DMI names. These are P00L, NPTP, 0PTP, DBOJ, INPT, PLT2, GE0M1, GE0M2, GE0M3, GE0M4, GE0MS, EDT, MPT, EPT, DIT, DYNAMICS, IFPFILE, AXIC, FORCE, MATPO0L, PCDB, XYCDB, CASECC, any DTI names, and SCRATCH1 thru SCRATCH9. 4. Field 3 of the header card must contain an integer of zero (0). 5. For symmetric matrices, the entire matrix must be input. 6. Only nonzero terms need be entered. 7. A blank field on this card is not equivalent to a zero. If zero input is desired, the appropriate type zero must be entered (i.e., 0.0 or O.ODO). 8. Complex input must have both the real and imaginary parts entered if either part is nonzero. 9. If A(Ix,J) is followed by "THRU" in the next field and an integer row number "IX' after the THRU, then A(Ix,J) will be repeated in each row thru IX. The "THRU" must follow an element value. For example, the entries for a real matrix RRR would appear as follows: (Continued) 2.4-156 (12/15/80)

BULK DATA DECK DMI (Cont.) 1 2 3 4 5 67 89 10 DMI NAME J I1 A(I1,J) I2 A(I2,J) DMI RRR I 2 1.0 THRU 10 12 2.0 |+1 These entries will cause the first column of the matrix RRR to have a zero in row 1, the values 1.0 in rows 2 through 10, a zero in row 11, and 2.0 in row 12. 10. Each column must be a single logical card. The terms in each column must be specified in increasing row number order. 11. The "F0RM" options 3, 4, 5, 7, and 8 are nonstandard forms and may be used only in conjunction with the modules indicated in the following table. Modules FORM Matrix Description ADD FBS MATPRN MPYAD 3 Diagonal X X X 4 Lower Triangular Factor X X 5 Upper Triangular Factor X X 7 Row X X 8 Identity X X X X 12. The maximum number of DMI and DTI data blocks is 79. 2.4-156a (12/15/80)

BULK DATA DECK Input Data Card F0RCE Static Load Description: Defines a static load at a grid point by specifying a vector. Format and Example: 1.2 3 4 5 6 7 8 9 10 F0RCE SID G CID F N1 N2 N3 FORCE 2 5 6 2.9 0.0 1.0 0.0 Field Contents SID Load set identification number (Integer > 0) G Grid point identification number (Integer > 0) CID Coordinate system identification number (Integer > 0, or blank) (Default 0) F Scale factor (Real) N1,N2,N3 Components of Vector measured in coordinate system defined by CID (Real; must have at least one nonzero component) Remarks: 1. The static load applied to grid point G is given by FN f = F N where N is the vector defined in fields 6, 7 and 8. 2. Load sets must be selected in the Case Control Deck (L0AD=SID) to be used by MSC/NASTRAN. 3. A CID of zero references the basic coordinate system. 2.4-179 (10/1/80)

BULK DATA DECK Input Data Card FORCEl Static Load, Alternate Form 1 Description: Used to define a static load by specification of a value and two grid points which determine the direction. Format and Example: 1 2 3 4 5 6 7 8 9 10 F0RCEI SID G F G1 G2 lllll F0RCE1 6 13 -2.93 16 13 Field Contents SID Load set identification number (Integer > 0) G Grid point identification number (Integer > 0) F Value of load (Real) G1, G2 Grid point identification numbers (Integer > O; G1 f G2) Remarks: 1. The direction of the force is determined by the vector from G1 to G2. 2. Load sets must be selected in the Case Control Deck (L0AD=SID) to be used by NASTRAN. 2.4-181

BULK DATA DECK Input Data Card GRAV Gravity Vector Description: Used to define gravity vectors for use in determining gravity loading for the structural model. Format and Example: 1 2 3 4 5 6 8 9 10 GRAV SID CID G N1 N2 N3 GRAV 1 3 32.2 0.0 0.0 -1.0 Field Contents SID Set identification number (Integer > 0) CID Coordinate system identification number (Integer > 0) G Gravity vector scale factor (Real) N1, N2, N3 Gravity vector components (Real; at least one nonzero component) Remarks: 1. The gravity vector is defined by g = G(N1, N2, N3). The direction of g is the direction of free fall. 2. A CID of zero references the basic coordinate system. 3. Gravity loads may be combined with "simple loads" (e.g., FORCE, MOMENT) only by specification on a L0AD card. That is, the SID on a GRAV card may not be the same as that on a simple load card. 4. Load sets must be selected in the Case Control Deck (LOAD=SID) to be used by MSC/NASTRAN. 5. At most nine GRAV cards can be selected in a given run either by Case Control or the LOAD Bulk Data card. Multiples or reflections of a given gravity load can be economically accomplished by use of the L0AD Bulk Data card. 6. In Type II solution sequence cyclic symmetry analyses, the T3 axis of the coordinate system referenced in field 3 must be parallel to the axis of symmetry. In the DIH type of cyclic symmetry the T1 axis must, in addition, be parallel to Side 1 of segment 1R of the model. 2.4-201 (5-15-80)

BULK DATA DECK Input Data Card GRDSET Grid Point Default Descriotion: Defines default options for fields 3, 7 and 8 of all GRID cards. Format and ExamDle: I 1 2 3 4 5 6 7 8 9 10 GRDSET CP CD PS GRDSET \ _16 ___ ___i3 1 32 345611m Field Contents CP Identification number of coordinate system in which the location of the grid point is defined (Integer > 0) CD Identification number of coordinate system in which displacements are measured at grid point (Integer > 0) PS Permanent single-point constraints associated with grid point (any of the digits 1-6 with no imbedded blanks) (Integer > 0). Remarks: 1. The contents of fields 3, 7 or 8 of this card are assumed for the corresponding fields of any GRID card whose field 3, 7 and 8 are blank. If any of these fields on the GRID card are blank, the default option defined by this card occurs for that field. If no permanent single-point constraints are desired or one of the coordinate systems is basic, the default may be overridden on the GRIn card by making one of fields 3, 7 or 8 zero (rather than blank). Only one GRDSET card may appear in the user's Bulk Data Deck. 2. The primary purpose of this card is to minimize the burden of preparing data for problems with a large amount of repetition (e.g., two-dimensional pinned-joint problems). 3. At least one of the entries CP, CD, or PS must be nonzero. 2. 4-203

BULK DATA DECK Input Data Card GRID Grid Point Description: Defines the location of a geometric grid point of the structural model, the directions of its displacement, and its permanent single-point constraints. Format and Example: 1 2 3 4 5 6 7 8 9 10 GRID ID |CP xi X3 | CD _v S -SET | GRID 2 13 1.0 -2.0 3.0 $ 316 111 Field Contents ID Grid point identification number (1,000,000 > Integer > 0) CP Identification number of coordinate system in which the location of the grid point is defined (Integer > 0 or blank*) X1,X2,X3 Location of the grid point in coordinate system CP (Real) CD Identification number of coordinate system in which displacements, degrees of freedom, constraints, and solution vectors are defined at the grid point (Integer > 0 or blank*) PS Permanent single-point constraints associated with grid point (any of the digits 1-6 with no imbedded blanks) (Integer > 0 or blank*) SEID Superelement identification number (Integer > 0 or blank) Remarks: 1. All grid point identification numbers must be unique with respect to all other structural, scalar and fluid points. 2. The meaning of Xl, X2 and X3 depend on the type of coordinate system, CP, as follows: (see C0RDi card descriptions) Type X X2 X3 Rectangular X Y Z Cylindrical R e(degrees) Z Spherical R' 9(degrees) b(degrees) 3. The collection of all CD coordinate systems defined on all GRID cards is called the Global Coordinate System. All degrees-of-freedom, constraints, and solution vectors are expressed in the Global Coordinate System. 4. The SEID entry can be overridden by use of the SESET Bulk Data card. * See the GRDSET card for default options for fields 3, 7 and 8. 2.4-205 (2/15/70)

BULK DATA DECK I.:put Data Card L3AD Static Load Combination (Superposition) Dascriotion: Cefines a static load as a linear combination of load sets defined via F0RCE,',ME:'iT, F0RCE1, MfITENT1, F-RCE2, M0MENT2, PL0AD, PL0AD1, PL5AD2, PL0AD3, PL0AD4, PL^AOX, SL^AD, RF0RCE, and GRAV cards. Fnrmat and Examsl e: 1 2 3 4 5 6 7 8 9 10 LJAO SID S S Sl LI S2 L2 S3 L3 ['LD 101 L -O.5 1.0 3 6.2 4 ___ IABC... S4 - L4. -etc.- (etc.) Ffeld Contents SfD Load set identification number (Integer > 0) S Scale factor (Real) Sf Scale factors (Real) Lf Load set identification numbers defined via card types enumerated above (Integer > 0) Renarks: 1. The load vector defined is given by * {PI= Si {PL1} 2. The Li must be unique. The remainder of the physical card containing the last entry must be blank. 3. This card must be used if gravity loads (GRAV) are to be used with any of-the other types. 4. Load sets must be selected in the Case Control Deck (LOAD=SID) if they are to be applied to the structural model. 5. A L/AD card miay not reference a set identification numnber defined by another LVAD card. 6. There may be at most 300 (Si, Li) pairs. 2.4-213 (5-15-80)

BULK DATA DECK Input Data Card MAT1 Material Property Definition, Form 1 Description: Defines the material properties for linear, temperature-independent, isotropic materials. Format and Example: 1 2 3 4 5 6.7 8 9 10 MAT1 MID E G NU RHO A TREF GE MAT1 17 3.+7 ____0.33 4.28 6.5-6 5.37+2 0.23 ABC ST SC SS MCSID +8C 20.+4 15.+4 12. +4 1003 Field Contents MID Material identification number (Integer > 0) E Young's modulus (Real or blank) G Shear modulus (Real or blank) NU Poisson's ratio (-1.0 < Real < 0.5 or blank) RH0 Mass density (Real) A Thermal expansion coefficient (Real) TREF Thermal expansion reference temperature (Real) GE Structural element damping coefficient (Real) ST,SC,SS Stress limits for tension, compression, and shear (Real). (Used only to compute margins of safety in certain elements; they have no effect on the computational procedures.) MCSID Material Coordinate System identification number (Integer > 0 or blank) Remarks: 1. The material identification number must be unique for all MAT1, MAT2, MAT3 and MAT9 cards. 2. MAT1 materials may be made temperature dependent by use of the MATT1 card. 3. The mass density, RHO, will be used to automatically compute mass for all structural el ements. 4. Weight density may be used in field 6 if the value 1/g is entered on the PARAM card WTMASS, where g is the acceleration of gravity (see Section 3.1.5). 5. MCSID must be nonzero if the CURV module is used to calculate stresses or strains at grid points. 6. To obtain the damping coefficient, GE, multiply the critical damping ratio C/Co, by 2.0. (Continued) 2.4-215 (12/15/80)

BULK DATA DECK Input Data Card MPC Multipoint Constraint Description: Defines a multipoint constraint equation of the form ~jj uj j. Format and Example: 1 2 3 4 5 6 7 8 9 10 MPC ISID G C A __G C A ab MPC 3 28 3_____6.2 2 __________4.29 +8_________ +bc I< G C A -etc. - +B | 1 4 -2.91 _ __ Field Contents SID Set identification number (Integer > 0). G Identification number of grid or scalar point (Integer > 0). C Component number - any one of the digits 1-6 in the case of geometric grid points; blank or zero in the case of scalar points (Integer). A Coefficient (Real; the first A must be nonzero). Remarks: 1. The first coordinate in the sequence is assumed to be the dependent coordinate. A dependent degree of freedom assigned by one MPC card cannot be assigned dependent by another MPC card or by a rigid element. 2. Forces of multipoint constraint are not recovered. 3. Multipoint constraint sets must be selected in the Case Control Deck (MPC=SID) to be used by NASTRAN. 4. The m-set coordinates specified on this card may not be specified on other cards that define mutually exclusive sets. See Section 1.4.1 for a list of these cards. 2.4-245 ( 5-15-79)

BULK DA7A DECK Input Data Card PARAM Parameter Description: Specifies values for parameters used in solution sequence. Format and Examole: 1 2 3 4 5 6 7 8 9 TO PARAM N V V2 PARAM IRES 1 Field Contents N Parameter name (one to eight alphanumeric characters, the first of which is alphabetic). V1,V2 Parameter value based on parameter type as follows: Type VI V2 Integer Integer Blank Real, single-precision Real Blank BCD BCD Blank Real, double-precision Double-precision Blank Complex, single-precision Real Real Complex, double-precision Double-precision Double-precision Remarks: 1. Only parameters for which assigned values are allowed may be given values via the PARAM card. Section 5 describes parameters as used in CMAP. 2. See Section 3.1.3 for a list of parameters used in solution sequences which may be set by the user on PARAM cards. 2.4-271 (5-15-80)

BULK DATA DECK Input Data Card PBAR Simple Beam Property Description: Defines the properties of a simple beam (bar) which is used to create bar elements via the CBAR card. Format and Example: 1 2 3 4 5 6 7 8 9 10 PBAR PID MID A II 1 2 J NSM PSAR 39 6 2..9 5.97 I 1 I I 123 rl|C1 "C2 1~0 IJ _ 02 El I E2 I F1 F2 +23 Cl 2 2,0 4.0 K1 K2 12I I Field Contents PID Property identification number (Integer > 0) MID Material identification number (Integer > O) A Area of bar cross-section (Real) II, 12, 112 Area moments of inertia (Real) (11 > 0., 12., 2 > I2) J Torsional constant (Real) NSM Nonstructural mass per unit length (Real) K1, K2 Area factor for shear (Real) Ci,Di,Ei,Fi Stress recovery coefficients (Real) Remarks: 1. For structural problems, PBAR cards may only reference MAT1 material cards. 2. See Section 1.3.2 for a discussion of bar element geometry. 3. For heat transfer problems, PBAR cards may only reference MAT4 or MATS material cards. 4. The transverse shear stiffnesses in planes 1 and 2 are (K1)AG and (K2)AG, respectively. The default values for Kl and K2 are infinite; in other words, the transverse shear flexibilities are set equal to zero. K1 and K2 are ignored if 112 f O. 5. The stress recovery coefficients Cl and C2, etc., are the y and z coordinates in the BAR element coordinate system of a point at which stresses are computed. Stresses are computed at both ends of the BAR. 2.4-273 (5-15-80)

BULK DATA DECK Input Data Card PBEAM Beam Property Description: Defines the properties of a beam which is used to create beam elements via the CBEAM card. Format and Example: 1 2 3 4 5 6 7 8 9 10 PBEAM PID MID A(A) I1(A) I2(A) I12(A) J(A) NSM(A) PBEAM 39 6 _ 2.9.__.5.5,_97 - 123 C1(A) C2(A) D1(A) D2(A) E1(A) E2(A) F1(A) F2(A) +23 1 1 2.0 -4.0 _ r T 1234 ___ I S0 X/XB A Il 12 112 J NSM__ J +34 YES 1.0 5.3 56.2 78.6 __ 1 1345 [i C2 Dl D2 El E2 Fl F2 +45 2.5 -5.0 456 K1 K2 S1 S2 NS I(A) NSI(B) CW (A) CB) M1(A) M2(A) M1(B) M2(B) N1(A) N2(A) N1(B) N2(B) 67 0.5 0.05 Field Contents Default Values PID Property identification number (Integer > 0) Required MID Material identification number (Integer > 0) Required A(A) Area of beam cross section at end A (Real > 0.0) Required I1(A) Area moment of inertia at end A in plane 1 about the neutral Required axis (Real > 0.0) I2(A) Area moment of inertia at end A in plane 2 about the neutral Required axis (Real > 0.0) I12(A) Area product of inertia at end A (Real (1 I I 2 > 0) 0.0 1 2- 12 0) 0.0 J(A) Torsional stiffness parameter at end A (Real) J > 0.0 if 0.0 warping is present) NSM(A) Nonstructural mass per unit length at end A (Real) 0.0 Ci(A),Di(A) The y,z locations in element coordinates (see diagram y =zi=0.0 Ei(A),Fi(A) following Remarks) at end A for stress data recovery (Real) (Continued) 2.4.274a (12/15/80)

NASTRAN DATA DECK PBEAM (Cont.) Field Contents Default Values S0 Stress output request option (BCD) Required* YES: Stresses recovered at points C,D,E,F on next continuation card YESA: Stresses recovered at points with the same y,z location as end A. NO: No stresses or forces-are recovered (see Remark 9). X/XB Distance from end A in the element coordinate system divided by Required* the length of the element (Real > 0.0). See following figure. See Remark 5 A,I1,I2,I12, Area, moments of inertia, torsional stiffness parameter and See Remark 6 J,NSM nonstructural mass for the cross-section located at x (Real) (J > 0.0 if warping is present.) Ci,Di,Ei,Fi The y,z locations in element coordinates (see diagram following Remarks) for the cross-section located at x/xb. The values are fiber locations for stress data recovery (Real) K1,K2 Shear stiffness factor K in KAG for plane 1 and plane 2 (Real) 1.0,1.0 S1,S2 Shear relief coefficient due to taper for plane 1 and plane 2 (Real) 0.0,0.0 NSI(A),NSI(B) Nonstructural mass moment of inertia per unit length about 0.0, same as nonstructural mass center of gravity at end A and end B (Real) end A (see following figure) CW(A),CW(B) Warping coefficient for end A and end B (Real) 0.0, same as end A M1(A),M2(1), (y,z) coordinates of center of gravity of nonstructural mass 0.0 (no offset M1(B),M2(B) for end A and end B (see following figure) from shear center), same values as end A N1(A),N2(A) (y,z) coordinates of neutral axis for end A and end B (see 0.0 (no offset following figure) from shear center) same values as end A Remarks: 1. For structural problems, PBEAM cards may reference only MAT1 material cards. 2. For heat transfer problems, PBEAM cards may reference only MAT4 and MATS material cards. 3. If no stress data at end A is to be recovered, and a continuation card with the BCD stress option is present, the continuation card +23 which contains the C,D,E and F (y,z) coordinates may be omitted. *Required only if card is present. (Continued) 2.4-274b (12/15/80)

BULK DATA DECK PBEAM (Cont.) 4. If SO is VESA or NO, the following continuation card for the C,D,E and F locations must be omitted. If SO is YES, the continuation card for C,D,E and F must be the next card. 5. Continuation cards +34 and +45 shown above may be repeated nine (9) more times for intermediate x/sb values. The order of the +34, +45 pairs is independent of the x/xb value; however, one value of x/xb must be 1.0, corresponding to end B. 6. If any fields 4 through 9 are blank on the continuation card with the value of x/xb = 1.0, then the values for A, I1, 12, 112, J and NSM are set to the values given for end A. For the continuation cards which have intermediate values of x/xb between 0.0 and 1.0, and which use the default option (any of the fields 4 through 9 are blank), a linear interpolation between the values at ends A and B is performed to obtain the missing section properties. 7. If SO is YES, blank fields are defaulted to O.0. 8. Blank fields for K1, K2 are defaulted to 1.0. If a value of 0.0 is used for K1 and K2, the transverse shear flexibilities are set to infinity. 9. If end B forces are desired, put YESA in the BCD field even when no end A stress points are input. Zelem y.'l0 Gri_ point ma i S, na zma SS /wA offset ^eem \ S^-Ymb I Grid point u^~ <^ He., A\W offset X Grid point GB 2.4-274c (12/15/80 sffset Gr poinint GB 2.~-274c (12/15/80

BULK DATA DECK Input Data Card PL0ADX Pressure Load on TRIAX6 Element Description: Defines a pressure load to be used with the TRIAX6 axisymmetric element. Format and Example: 1 2 3 4 5 6 7 8 9 10 PLAODX |SID P1i' P3. G1.G2 G3... PL)ADX 200 13.5 10.5 10 20 130 Field Contents SID Load set identification number (integer > 0). P1^ P3 Pressure at first and third listed points (real). G1, G2, G3 Grid point identification numbers (integers > 0). Remarks: 1. This load card is designed for use with the TRIAX6 element. G1 and G3 should be vertices of a TRIAX6 element, and G2 the corresponding midside grid point. 2. The pressure is assumed to vary linearly between grid points 1 and 3. If G1, G2, G3 appear in the order determined by going around an element in counterclockwise direction, a positive pressure will produce a force directed toward the element interior. P3 G3 Axis of Revolution G4 G2 G5 -G6 3. The pressure is input as force per unit area. 2.4-293

BULK DATA DECK Input Data Card PQDMEM1 Quadrilateral Membrane Property Description: Used to define the properties of a quadrilateral membrane. Referenced by the CQDMEM card. No bending properties are included. Format and Example: 1 2 3 4 5 6 7 8 9 10 PQDMEM PID MID T NSM 1 PQDMEM 235 2 0.5 0.0 Field Contents PID Property identification number (Integer > 0) MID Material identification number (Integer > 0) T Thickness of membrane (Real > 0.0) NSM Nonstructural mass per unit area (Real) Remarks: 1. All PQDMEM cards must have unique property identification numbers. 2.4-306a (5-15-80)

BULK DATA DECK Input Data Card PR)D Rod Property Description: Defines the properties of a rod which is referenced by the CR0D card. Format and Examole: 1 2 3 4 5 6 7 8 9 10 PRoD PID'|. ID A J C NSM i JFR0D j 17 | 23 42.6.17.92 4.236 0.5.. Field Contents PID Property identification number (Integer > 0),MID Material identification number (Tnteger > 0) A Area of rod (Real) J Torsional constant (Real) C Coefficient to determine torsional stress (Real) NSM Nonstructual mass per unit length (Real) Remarks: 1. PR0D cards must all have unique property identification numbers. 2. For structural problems, PR0D cards may only reference MAT1 material cards. 3. For heat transfer problems, PR0D cards may only reference MAT4 or MAT5 cards. 4. The formula used to compute torsional stress is J where Mt is the torsional moment. 2.4-323

BULK DATA DECK Input Data Card PSHELL Shell Element Property Description: Defines the membrane, bending, transverse shear, and coupling properties of thin shell elements. Format and Example: 1 2 3 4 5 6 7 8 9 10 PSHELL PID MIDi T MID2 121/T3 MID3 TS/T NSM PSHELL 203 204 1.90 205 1.2 206 0.8 6.32 BCD Zi. Z2 MID4 +CD +.95 -.95 Field Contents PID Property identification number (Integer > 0) MID1 Material identification number for membrane (Integer > 0 or blank) T Default value for membrane thickness (Real) MID2 Material identification number for bending (Integer > 0 or blank) 12I/T3 Bending stiffness parameter (Real or blank, default = 1.0) MID3 Material identification number for transverse shear (Integer > 0 or blank), must be blank unless MID2 > O) TS/T Transverse shear thickness divided by membrane thickness (Real or blank, default.833333) NSM Nonstructural mass per unit area (Real) Z1,Z2 Fiber distances for stress computation. The positive direction is determined by the righthand rule and the order in which the grid points are listed on the connection card. (Real or blank, defaults are -1/2 T for Z1 and 1/2 T for Z2.) MID4 Material identification number for membrane-bending coupling (Integer > 0 or blank, must be blank unless MID1 > 0 and MID2 > 0, may not equal MID1 or MID2) Remarks: 1. All PSHELL property cards must have unique identification numbers. 2. The structural mass is computed from the density using the membrane thickness and membrane material properties. 3. The results of leaving an MID field blank are: MID1 No membrane or coupling stiffness. MID2 No bending, coupling, or transverse shear stiffness. MID3 No transverse shear flexibility. MID4 No bending-membrane coupling. (Continued) 2.4-326a (12/15/80)

NASTRAN DATA DECK PSHELL (Cont.) 4. The continuation card is not required. 5. The structural damping (for dynamics rigid formats) uses the values defined for the MID1 material. 6. The MID4 field should be left blank if the material properties are symmetric with respect to the middle surface of the shell. 7. This card is used in connection with the CTRIA3, CTRIA6, CQUAD4 and CQUAD8 cards. 8. For structural problems, PSHELL cards may reference MAT1 or MAT2 material property cards. 9. For heat transfer problems, PSHELL cards may reference MAT4 or MAT5 material property cards. 2.4-326b (12/15/80)

BULK DATA DECK Input Data Card PSOLID Properties of HEXA and PENTA Solid Elements Description: Defines the properties of solid elements. Referenced by CHEXA and CPENTA cards. Format and Example: 1 2 3 4 5 6 7 8 9 10 PS0LID PID MID CoROM IN PS0LID 2 100 6 Field Contents PID Property identification number (Integer > 0) MID Identification number of a MAT1, MAT4, MAT5, or MAT9 card (Integer > 0) C'RDM Identification number of material coordinate system (Integer, Default - -1) IN Integration network (Integer, 0, 2, 3, or blank) Remarks: 1. PS0LID cards must have unique ID numbers. 2. Either isotropic (MAT1,, OR MAT4) or anisotropic (MATS or MAT9) materials may be referenced. 3. See the CHEXA or CPENTA card for the definition of t;e elemrent coordinate system. The material coordinate system may be the basic system (0), any defined system (Integer > 0) or the element coordinate system (-1 or blank). 4. Stress components are output in the material coordinate system for the corner points and the midpoint of the element. 5. There is a choice of tao integration networks for each element. HEXA PENTA IN=2 8 point 6 point IN=3 27 point 9 point The recommendation is to use the default (zero or blank) which selects IN=2 if all edge points are deleted and IN=3 if any edge points are present. 2.4-326c (5-15-80)

BULK DATA DECK Input Data Card RBAR Rigid Bar Description: Defines a rigid bar with six degrees of freedom at each end. Format and Example: 1 2 3 4 5 6 7 8 9 10 RBAR IEID IGA GB CNA CNB CMA ICMB RBAR 5 11 2 34 123 Fie d Contents EID Identification number of rigid element. GA, GB Grid point identification number of connection points (integer > 0). CNA, CNB Independent degrees of freedom in the global coordinate system for the element at grid points GA and GB, indicated by any of the digits -6 with no imbedded blanks (integer > 0 or blank). See Remark. 1. CMA, CMB Component numbers of dependent degrees of freedom in the global coordinate system assigned by the element at grid points GA and GB, indicated by any of the digits 1 - 6 with no imbedded blanks (integer > 0 or blank). See Remarks 2 and 3. Remarks: 1. The total number of components in CNA and CNB must equal six; for example, CNA = 1236, CNB = 34. Furthermore, they must jointly be capable of representing any general rigid body motion of the element. 2. If both CMA and CMB are zero or blank, all of the degrees of freedom not in CNA and CNB will be made dependent; i.e., they will be made members of the {u;Tset. 3. The m-set coordinates specified on this card may not be specified on other cards that define mutually exclusive sets. See Section 1.4.1 for a list of these cards. 4. Element identification numbers must be unique. 5. Rigid elements, unlike MPC's, are not selected through the Case Control Deck. 6. Forces of constraint are not recovered. 7. Rigid elements are ignored in heat transfer problems. 8. See Section 2.10 of the Application Manual for a discussion of rigid elements. 2.4-359 (5-15-80)

BULK DATA DECK Input Data Card SEQGP Grid and Scalar Point Resequencing Descriotion: Used to order the grid points and user-supplied scalar points of the problem. Tne purpose of this card is to allow the user to reidentify the formation sequence of the grid and scalar points of his structural model in such a way as to optimize bandwidth which is essential for efficient solutions by the displacement method. For-at and Examole: 1. 2 3 4 5 6 7' 8 9 10 SECGP I ID SEQID 10o SEQID 10 SEQID ID SEQID SEQG? 5392 15.6 1.596 0.2 2 1.9.2.6 3 2 Field Contents ID Grid or scalar point identification number (Integer > 0). SEQID Sequenced identification number (a special number described below). Remarks: 1. ID is any grid or scalar point identification number which is to be reidentified for sequencing purposes. The grid point sequence number (SEQID) is a special number which may have any of the following forms where X is a decimal integer digit - XXXX.X.X.X, XXXX.X.X, XXXX.X or XXXX where any of the leading X's may be omitted. This number must contain no fmbedded blanks. The leading character must not be a decimal point. 2. If the user wishes to insert a grid point between two already existing grid points, such as 15 and 16, for example, he would define it as, say 5392, and then use this card to insert grid point number 5392 between them by equivalencing it to, say 15.6. All output referencing this point will refer to 5392. 3. The SEQID numbers must be unique and may not be the same as a point ID which is not being changed. No grid point ID may be referenced more than once. 4. No continuation cards (small field or large field) are allowed with either the SEQGP or SEQEP card. 5. From one to four grid or scalar points may be resequenced on a single card. 6. If a grid point ID is referenced more than once, the last reference will determine its sequence. 2.4-389 (5-15-79)

BULK DATA DECK Input Data Card SPC Single-Point Constraint Description: Defines sets of single-point constraints and enforced displacements. Format and Example: 1 2 3 4 5 6 7 8 9 10 SPC SID GI C D G C 0 SPC 2 1 32 1436 -2.6 5 + 2.9 Field Contents SID Identification number of single-point constraint set (Integer > 0). G Grid or scalar point identification number (Integer > 0). C Component number (6 > Integer > O; up to six unique digits may be placed in the field with no imbedded blanks). D Value of enforced displacement for all coordinates designated by G and C (Real). Remarks: 1. Coordinates specified on this card form members of a mutually exclusive set. They may not be specified on other cards that define mutually exclusive sets. See Section 1.4.1 for a list of these cards. 2. Single-point forces of constraint are recovered during stress data recovery. 3. Single-point constraint sets must be selected in the Case Control Deck (SPC=SID) to be used by NASTRAN. 4. From one to twelve single-point constraints may be defined on a single card. 5. SPC degrees of freedom may be redundantly specified as permanent constraints on the GRID card. 6. Continuation cards are not allowed. 2.4-399 ( 5-15-79)

BULK DATA DECK Input Data Card SPC1 Single-Point Constraint, Alternate Form Description: Defines sets of single-point constraints. Format and Example: SPC1 SID C G1 G2 G3 G4 G5 G6 abc SPC1 3 2 1 3 10 9 6 5 ABC +bc G7 G8 G9 -etc+BC 2 8 Alternate Form SPC1 SID C GID1 "THRU" GID2 SPC1 313 12456 6 THRU 32 Field Contents SID Identification number of single-point constraint set (Integer > 0). C Component number (any unique combination of the digits 1-6 (with no imbedded blanks) when point identification numbers are grid points; must be null if point identification numbers are scalar points). Gi,GIDi Grid or scalar point identification numbers (Integer > 0). Remarks: 1. Note that enforced displacements are not available via this card. As many continuation cards as desired may appear when "THRU" is not used. 2. Coordinates specified on this card form members of a mutually exclusive set. They may not be specified on other cards that define mutually exclusive sets. See Section 1.4.1 for a list of these cards. 3. Single-point constraint sets must be selected in the Case Control Deck (SPC=SID) to be used by NASTRAN. 4. SPC degrees of freedom may be redundantly specified as permanent constraints on the GRID card. 5. If the alternate form is used, points in the sequence GID1 thru GID2 are not required to exist. Points which do not exist will collectively produce a warning message but will otherwise be ignored. 2.4-401 (5-15-79)

.BULK DATA DECK Input Data Card SPCADD Single-Point Constraint Set Combination Description: Defines a single-point constraint set as a union of single-point constraint sets defined via SPC or SPC1 cards. Format and Example: 1 2 3 4 5 6 7 8 9 10 SPCADD SID Si I12 S3 S4 S5 S6 i 7aS7 abc SPCADO 101 3 92 9 11 +bc S8 S9 -etc. I l -etc.Field Contents SID Identification number for new single-point constraint set (Integer > 0). Si Identification numbers of single-point constraint sets defined via SPC or SPC1 cards (Integer > O; SID # Si). Remarks: 1. Single-point constraint sets must be selected in the Case Control Deck (SPC=SID) to be used by NASTRAN. 2. No Si may be the identification number of a single-point constraint set defined by another SPCADD card. 3. The Si values must be unique. 4. SPCADD cards take precedence over SPC or SPC1 cards. If both have the same set ID, only the SPCADD card will be used. 2.4-403 ( 5-15-79)

PLoT 1^i - e \ OLgO"D

PLOTTING 4.1 PLOTTING 4.1.1 General Capability MSC/NASTRAN provides the capability for generating the following kinds of plots: 1. Undeformed geometric projections of the structural model. - 2. Static deformations of the structural model by either displaying the deformed shape (alone or superimposed on the undeformed shape), or displaying the displacement vectors at the grid points (superimposed on either the deformed or undeformed shape). 3. Modal deformations resulting from real or complex eigenvalue analysis by the same options stated in 2 above. Complex modes for flutter analysis may be plotted for any user-chosen phase lag. 4. Deformations of the structural model for transient response or frequency response by displaying either vectors or the deformed shape for specified times or frequencies. 5. X-Y graphs of transient response or frequency response. 6. V-F and V-G graphs for flutter analysis. - 7. Contour plots of displacements, temperature and stress on the structure. Structural plots (items 1-4) are discussed in Section 4.2, while X-Y plots (item 5) are discussed in Section 4.3. Requests for structure plots or X-Y plots are accomplished in the Case Control Deck by submitting a structure plot request or an X-Y output request. The optional PL0TID card is considered to be part of the plot request (although it must precede any OUTPUT(PLPT), 0UTPUT(XY0UT) or 0UTPUT(XYPL0T) cards. See Section 2.3 for a description of the PL0TID card. Plot requests are separated from case control by the 0UTPUT(PLOT), 0UTPUT(XYPLOT) or 0UTPUT(XY0UT) cards. Data above this card (except PL0TID) will not be recognized by the plotter, even though it may have the same name (for example, the SET card). 4.1-1 (5-15-80)

PLOTTING 4.1.2 Superelement Plotting Plotting may take place at four different places in the superelement solution sequences. There are two plot commands, SEPL0T and SEUPPL0T, used with other case control and PARAM cards to control the type of plot to be prepared. A flow chart of the solution process is given in Figure 1. Undeformed structure plots are made early during the SEMG process. They can be made for either one superelement only, or for a superelement and all of its upstream members, as controlled by PARAM, PLOTSUP (see Figure 1, block 1). Plots are made for superelements selected for the SEMG process, and listed in a plot request headed by SEPL0T SEID, where SEID is the superelement identification number (see Section 4.3.8 for an example of a complete superelement plot packet). Undeformed plots can be used to check geometry and connectivity, and do not require the presence of property or material cards. A branch to the end of the loop immediately after the plot module can be made by use of the PARAM,PL0T,-1 card. In dynamic analysis, solution set XY-plots are requested by use of SEPL0T 0 (see Figure 1, block 2). In the data recovery phase, XY-plots and deformed structure plots for elements in one superelement only are requested by the SEPL0T command (see Figure 1, block 3). Deformed structure plots for a superelement and all its upstream plots are requested with the SEUPPL0T command (see Figure 1, block 4). Undeformed plots require the presence of the SEMG case control card for all superelements listed on SEPL0T requests. XY-plots and deformed structure plots are regarded as output requests, and will result in automatic execution of the data recovery loop as is required to produce the plots requested. For example, the command SEUPPL0T 0 will result in data recovery being performed on the entire model, even in the absence of any other output requests. For SEUPPL0T requests, if SUBC0M or SYMC0M subcases are used, each superelement must have identical SUBC0M and SYMC0M structure. 4.1-2 (10/1/80)

PLOTTING Case Control I. GENERATION, ASSEMBLY, REDU(TION LOOP Commands PARAM,PLOTSUP,O! PARAM,PL0TSUP,-1 Prepare plotting data da r on (default) superelement on ly "1. Undeformed Plots SEMG, PARAM,PL0)T,-1 _ SEPL0T SEID PARAM,PL0T,O (Oeimult) S SEUPPLOT SEID Generate, assemble and reduce one superelement II. RESIDUAL STRUCTURE ONLY P ASE Eigensolutions, Forced Response 2. Solution Set XY Plots SEPL0T 0 Continued Figure 1. Superelement Plot Control 4.1-3 (5-15-80)

PLOTTING III. DATA RECOVERY LOOP Recover interior di slacements, stresses 3. XYPLOTS, Deformed structure SEPL0T SEID plots (Local) IV. UPSTREAM PLOT LOOP Assemble Upstream Plot Vectors 4. Deformed Structure SEUPPL0T SEID Plots (Upstream) FINIS Figure 1. Superelement Plot Control (Cont.). 4.1-4 (5-15-79)

PLOTTING 4.2 STRUCTURE PLOTTING In order to assist users both in the preparation of the analytical model and in the interpretation of output, the structure plotter provides the following capabilities for undeformed structures: 1. Connect the grid points in a predetermined manner using the structural elements or PL0TEL elements. 2. Place a symbol at connected grid point locations. (optional) 3. Identify connected grid points by placing the grid point identification number to the right of the grid point locations. (optional) 4. Identify elements by placing the element identification number and element symbol at the center of each element. (optional) 5. Identify connected grid points that are constrained to zero with single point constraints by placing the grid point number and integers defining the constrained degrees of freedom to the right of the grid point locations. 6. Identify elements by placing the element identification number, the element symbol, and the property identification number at the center of each element. 7. Draw an outline of the structure. (optional) 8. Shrink one- and two-dimensional elements by a fraction. (optional) The following capabilities are provided for deformed structures: 1. Connect the deflected grid points in a predetermined manner using the structural elements or PL0TEL elements. 2. Place a symbol at the deflected grid point location. (optional) 3. Identify the deflected grid points by placing the grid point identification number to the right of the deflected grid point location. (optional) 4. Draw lines originating at the undeflected or deflected grid point location, drawn to user-specified scale, representing the X, Y, Z components or resultant summations of the grid point deflections. 5. Draw an outline of the structure. (optional) 6. Shrink one- and two-dimensional elements by a fraction. (optional) 7. Contours of element stresses. 8. Contours of displacements or temperatures. 4.2-1 (10/1/80)

PLOTTING The above plots are available in either orthographic, perspective, or stereoscopic projections. Stereoscopic plots are normally made only on microfilm plotters since a stereoscopic viewer or projector must be used to obtain the stereoscopic effect. A request for structure plotting is made in the Case Control Deck by means of a plot request which includes all cards from an 0UTPUT(PL0T) card to either a BEGIN BULK or 0UTPUT(XYOUT) [or 0UTPUT(XYPL0T)] card. It should be noted that only elements can be plotted. Grid points that are not associated with elements cannot be plotted. Grid points may be connected with PLOTEL elements for plotting purposes. The data card format is free-field, subject to rules in paragraphs below. The cards are basically sequence-dependent even though some interchanging in sequence of defining parameters is permissible. The elements and grid points to be plotted may be defined anywhere in the plot package, but the parameters describing the characteristics of the plot are evaluated on the current basis every time a PL0T or FIND card (see Section 4.2.2.2) is encountered. In order to minimize mistakes, it is suggested that a strict sequence dependency be assumed. 4.2.1 General Rules 4.2.1.1 Rules for Free-Field Card Specification 1. Only columns 1 thru 72 are available. Any information specified in columns 73 thru 80 will be ignored. 2. If the last character on a card is a comma (not necessarily in column 72), the next card is a continuation of this physical card. Any number of continuation cards may be specified, and together they form a logical card. 3. The mnemonics or values can be placed anywhere on the card, but must be separated by delimiters. 4. The following delimiters are used: a. blank b. comma (,) c. left parenthesis ( ] d. right parenthesis [)] e. equal sign (=) All of these delimiters can be used as needed to aid the legibility of the data. 4.2-2 (2/2/81)

STRUCTURE PLOTTING 4.2.1.2 Plot Request Packet Card Format In the plot request packet card descriptions presented in Section 4.2.2, the following notation will be used to describe the card format: 1. Upper case letters must be punched exactly as shown. 2. Lower-case letters indicate that a substitution must be made. 3. Braces { } indicate that a choice of the contents is mandatory. 4. Brackets [ ] contain an option that may be omitted or included by the user. 5. Underlined options or values are those for which a default option or an initialized (or computed) value was programmed. 6. A physical card consists of information punched in columns 1 thru 72 of a card. 7. A logical card may consist of more than one physical card through the use of conti nuati on cards. 8. Numerical values may usually be either integer or real numbers, even though a specific type is at times suggested in order to conform to the input in other sections of the program. If the first four characters of a mnemonic are unique relative to all other case control cards, the characters following can be omitted. 4.2.1.3 Plot Titles Up to four lines of title information will be printed in the lower left-hand corner of each plot. The text for the top three lines is taken from the TITLE, SUBTITLE, and LABEL cards in the Case Control Deck. (See Section 2.3 for a description of the TITLE, SUBTITLE, and LABEL cards.) The text for the bottom line may be of two forms depending on the type plot requested. One form contains the word UNDEF0RMED SHAPE. The other form contains the type of plot (statics, modal, etc.), subcase number, load set or mode number, frequency or eigenvalue or time, and (for complex quantities) the phase lag or magnitude. 4.2-3 (2/2/81)

PLOTTING The sequence number for each plot is printed in the upper corners of each frame. The sequence number is determined by the relative position of each PL0T execution card in the plot package. The information on the PTITLE card will be printed on the line below the sequence number. The date and (for deformed plots) the maximum deformation are also printed at the top of each frame. 4.2.2 Plot Request Card Descriptions The general form for each card of the plot request is shown enclosed in a box. Description of the card contents then follows for each card. 4.2.2.1 SET Definition Cards These cards specify sets of elements, corresponding to portions of the structure, which may be referenced by PL0T and FIND cards. The SET card is required. Each set of elements defines by implication a set of grid points connected by those elements. The set may be modified by deleting some of its grid points. The elements are used for creating the plot itself and element labeling while the grid points are used for labeling, symbol printing, and drawing deformation vectors. SET i CINCLUDE][ELEMENTS] jl, J2, j3 THRU J4, j5, etc. INCLUDE EXCLUDE ELEMENTS EXCEPT GRID P INTS k1, k2 3 k3 THRU k4, k5 etc. i = set identification number (positive integer, unique for each set) j = element identification numbers or element types k = element identification numbers or grid point identification numbers or element types 4.2-4 (2/2/81)

STRUCTURE PLOTTING Permissible element types and plotting labels are: Element Type Plot Label Element Type Plot Label AER01 AE PENTA PA AXIF2 A2 PL0TEL PL AXIF3 A3 QUAD4 Q4 AXIF4 A4 QUAD8 Q8 BAR BR ROD RD BEAM BM SHEAR SH BEND BD SLOT3 S3 C0NE CN SLOT4 S4 C0NR0D CR TETRA TE DUMi Di TRAPRG TA FLMASS FM TRIAX6 D1 FLUID2 F2 TRIA3 T3 FLUID3 F3 TRIA6 T6 FLUID4 F4 TRIM6 02 HBDY* HB TRIARG Ti HEXA HA TUBE TU HEX20 XQ VISC VS *The boundary condition elements for the MSC/NASTRAN heat conduction analyzer can be plotted for undeformed structure plots. There are several types of HBDY elements, as follows: Type P0INT LINE AREA3 AREA4 REV ELCYL FTUBE No. of Primary Grid Pts. 1 2 3 4 2 2 2 Normals Available yes yes yes yes no no yes The secondary grid points are used for "ambient" conditions and are ignored by the plotter. Type P0INT must have a nonzero associated area (see AF on the associated PHBDY data card) and a defined normal direction (see V1, V2, V3 on the CHBDY data card) to be plotted. It is plotted as a hexagon with approximately the correct area. Type LINE must have a nonzero width (see AF on the associated PHBDY data card) and a normal defined in order to plot. 4.2-4a (2/2/81)

PLOTTING ALL may be used to select all permissible element types. INCLUDE may be used at any time for element information. EXCLUDE can be used to delete elements or element types. All grid points that are associated with deleted elements are also deleted. INCLUDE or EXCLUDE, when used for grid points, will have the first reference to a grid point used. EXCEPT is a modifier to an INCLUDE or an EXCLUDE statement for elements. 4.2-4b (2/2/81)

STRUCTURE PLOTTING THRU is used to indicate all of the integers in a sequence of identification numbers, starting with the integer preceding THRU and ending with the integer following THRU. The integers in the range of the THRU statement need not be consecutive, e.g., the sequence 2, 4, 7, 9 may be specified as 2 THRU 9. THRU is not applicable if element types are specified. THRU should not be used for large open sets because of high computation cost. Each SET must, be a logical card. Redefinition of sets previously defined is not permitted; however, there is no restriction on the number of sets. The sets of identification numbers can be assembled by use of the word ALL, or by individually listing the integers in any order such as 1065, 32, 46, 47, 7020, or by listing sequences using THRU, EXCLUDE, and EXCEPT such as 100 THRU 1000 EXCEPT 182 EXCLUDE 877 THRU 911. Examples of SET cards: 1. SET 1 INCLUDE 1, 5, 10 THRU 15 EXCEPT 12 (Set will consist of elements 1, 5, 10, 11, 13, 14 and 15) 2. SET 2 = ROD, CONROD, EXCEPT 21 (Set will consist of all ROD and CONROD elements except element 21) 3. SET 10 SHEAR INCLUDE GRID POINTS 35, 36 EXCLUDE GRID POINTS 20, 30 THRU 60, INCLUDE ELEMENTS 70 THRU 80. (This set will include all shear elements plus elements 70 thru 80, and the associated grid point set will contain all grid points connected by these elements. Grid points 20, 30 thru 34 and 37 thru 60 will not appear on the plots with their symbols, labels, or vectors. 4. SET(15) = (15 THRU 100) EXCEPT (21 THRU 25) (This set will include all elements from 15 to 20 and from 26 to 100) 5. SET 2 = ALL EXCEPT BAR (This set will include all elements except bars) Notes: The equal signs, commas, and parentheses above are delimiters and are not required, because blanks also serve as delimiters. Plotter sets must be specified after the 0UTPUT(PLOT) card. Case Control sets are ignored. 4.2.2.2 Cards Defining Parameters These cards specify how the structure will be plotted, i.e., type of projection, view angles, scales, etc. All the multiple choice parameters are defaulted to a preselected choice if not specified. Each parameter requiring a numerical value that is not specified by the user can either be established internally in the program by means of the FIND card or can assume default values. The FIND card is used to request that the program select a SCALE, ORIGIN, and/or VANTAGE POINT to allow the construction of a plot in a user-specified region of the paper or film. The FIND card is described at the end of this section, following the discussion of the associated parameters. 4.2-5

PLOTTING The parameter cards are listed here in a logical sequence; however, they need not be so specified. Any order may be used, but if a parameter is specified more than once, the value or choice stated last will be used. Each parameter may be either an individual card, or any number of them may be combined on one logical card. All the parameters used in the generation of the various plots will be printed out as part of the output, whether they are directly specified, defaulted or established using the FIND card. Initialization of parameters to default values occurs only once. Subsequently, these values remain until altered by a direct input. The only exceptions are the view angles, scale factors, vantage point parameters, and the origins. Whenever the plotter or the method of projection is changed, the view angles are reset to the default values, unless they are respecified by the user. In addition, the scale factors, vantage point parameters, and the origin must be redefined by the user. PL0TTER { AST The default option (NAST) is the NASTRAN general purpose plotter. The NASTRAN plotter option prepares a general purpose plot file which must be interpreted by a postprocessing program. The format of the NASTRAN general purpose plot file is described in Section 6 of the Programmer's Manual. See Section 2.3 of the Application Manual for a description of NASPL0T, a utility routine for interpreting the NASTRAN plot file. SC selects the Stromberg Carlson 4020 microfilm plotter. 0RTH0GRAPHIC PERSPECTIVE STERE0SC0PIC The default option is orthographic projection. See Section 13 of the Theoretical Manual for a discussion of the various projections. This card is optional. SYMMETRIC AXES r,s,t ANTISYMMETRIC VIEW y,,a r, s, t = X or MX, Y or MY, Z or MZ (where "M" implies the negative axis) y, B, a = three angles of rotation in degrees (real numbers) 4.2-6 (7/1/80)

STRUCTURE PLOTTING These two parameter cards define the orientation of the object in relation to the observer, that is the angles of view. Both of these cards are optional. Defining the observer's coordinate system as R, S, T and the basic coordinate system of the object as X, Y, Z, the angular relationship between the two systems is defined by the three angles y, 8 and a as follows: T / r — n of 0<^ ^^\ -- Direction of View /_Ja 7 (Always in negative R-direction. 7/ / / The projection plane is always in, A/^ /- / or parallel to, the S-T plane.) R Using the above convention, y and 8 represent the angles of turn and tilt. The default values are: y = 34.27~ B = 23.17~ for orthographic and stereoscopic projections, 0.0~ for perspective projection a = 0.0~ The order in which y, 8, and a are specified is critically important as illustrated in Figure 1, at the end of this section. Also, see Section 13.1.1 of the Theoretical Manual. The AXES card can be used to preposition the object in 90~ increments in such a manner that only rotations less than 90~ are required by the VIEW card to obtain the desired orientation. This is accomplished by entering X, Y, Z, MX, MY or MZ in the fields corresponding to R, S, T axes, where MX, MY and MZ represent the negative, X, Y and Z axis directions respectively. The default values are X, Y, Z. Tle 14 values can also be used to define left-handed cooruinate systems. (Note that the default system is right-handed.) An undeformed or deformed plot of the symmetric portion of an object can be obtained by reversing the sign of the axis that is normal to the plane of symmetry. In the case of multiple planes of symmetry, the signs of all associated planes should be reversed. The ANTISYMMETRIC option should be specified when a symmetric structure is loaded in an unsymmetric manner. This will cause the deformations to be plotted antisymmetrically with respect to the specified plane or 4.2-7

PLOTTING planes. Since the AXES card applies to all parts (SETS) of a single frame, symmetric and antisymmetric combinations cannot be made with this card (see the symmetry option on the the PL0T execution card in Section 4.2.2.3). MAXIMUM DEF0RMATI0N 5% structure dimension t The value of d represents the length to which the maximum displacement component is scaled in each subcase. The maximum deformation must be specified in units of the structure (not inches of paper). / actual displacement d actual maximum displacement The default value for d is 5% of the maximum dimension of the structure. This data is necessary since the actual deformations are usually too small to be distinguishable from the undeformed structure if they were plotted to true scale. If FIND card parameters are to be based on the deformed structure, the FIND card must be preceded by the MAXIMUM DEF0RMATI0N card. SCALE aC, b] a = real number representing scale to which the model is drawn b = ratio of model size/real object size (stereoscopic projection only) For orthographic or perspective projections, the scale "a" is the ratio of the plotted object in inches to the real object in the units of the structural model, i.e., one inch of paper equals one unit of structure. For stereoscopic projection, the stereoscopic effect is enhanced by first reducing the real object to a smaller model (scale "b"), and then applying scale "a." The ratio of plotted/real object is then the product a x b. A scale must be defined in order to make a plot. However, the SCALE card is not recommended for general use. See the FIND card described at the end of this section in order to have the scale determined automatically. 4.2-8 (1/15/81)

STRUCTURE PLOTTING If the NASTRAN general purpose plotter is used in combination with the MSC/NASPL0T postprocessing routine, the user may determine a scale as follows: 20x 1 a = p x-x 1 7 SF(NASPL0T) where p = ratio of object size to plot size 20 = ratio of default PAPER SIZE to default NASPL0T frame size -7 SF(NASPLOT) - scale factor on NASPL0T data card (default = 1.0) (1.0, 1.0, 1.0 I DIST0RTI0N -,Y i This card specifies the distortion factor in the X, Y and Z directions. The X, Y and Z coordinates are multiplied by the distortion factors prior to any other scaling operations. The default values are all 1.0. If this card is used, all three values must be specified even though one or two of the values are the default. CSCALE 0.5 This card is used to control the spacing of characters when plots are made with the NASTRAN plotter and they are postprocessed with the MSC/NASPL0T routine. For example, if the SCALE FACT0R on the NASPL0T data card is 2.0, a value for cs of 0.5 will result in characters of default size (.07 inches) at the regular spacing. On the other hand, if the user wishes to double the size of both the plot and the characters, the SCALE FACTOR and the CSCALE FACT0R on the NASPL0T data card should both be set equal to 2.0. The CSCALE card must immediately precede the PLOTTER selection card. If a second CSCALE card is used, a second PL0TTER selection crd must also be used. ORIGIN i, u, v 4.2-9 (1/15/81)

PLOTTING i - origin identification number (any positive integer) u = horizontal displacement of paper origin from RST origin v = vertical displacement of paper origin from RST origin In the transformation performed for any of the three projections, the origins of both the object (XYZ system) and of the observer (RST system) are assumed to be coincident. This card refers to the paper origin. It represents the displacement of the paper origin (lower left hand corner) from the RST origin. The units are inches and are not subject to the scaling of the plotted object. The ORIGIN card is not recommended for general use. See the FIND card described at the end of this section in order to have the origin located so as to place the plotted object in the center of the image area. Ten (10) origins are permitted to be active at one time. However, any one can be redefined at any time. An eleventh origin is also provided if more than ten origins are erroneously defined (i.e., only the.last of these surplus origins will be retained). CAUTION: When a new projection or plotter is called for, all previously-defined origins are deleted. VANTAGE P0INT ro, so, to, do [, sor] (perspective and stereoscopic projections only) ro = R-coordinate of the observer so = S-coordinate of the observer in perspective projection or S-coordinate of the left eye of the observer in the stereoscopic projection to = T-coordinate of the observer do = R-direction separation of the observer and the projection plane sor = S-coordinate of the right eye of the observer in the stereoscopic projection (required in stereoscopic projection) This card defines the location of the observer with respect to the structural model. A vantage point is required for either perspective or stereoscopic projections. The VANTAGE POINT card is not recommended for general use. See the FIND card described at the end of this section. A theoretical description of the vantage point is contained in Section 13 of the Theoretical Manual. 4.2-10 (1/15/81)

STRUCTURE PLOTTING PROJECTION PLANE SEPARATION { 2.0 (stereoscopic projections only) This card specifies the R-direction of the observer and the projection plane when a VANTAGE P0INT card is not used. The default value is 2.0. The PR0JECTI0N PLANE SEPARATI0N card is not recommended for general use. See the FIND card described at the end of this section. A theoretical description of projection plane separation is contained in Section 13 of the Theoretical Manual. 0CULAR SEPARATION j 2.756 n (stereoscopic projection only) Ocular separation - S-coordinate separation of the two vantage points in the stereoscopic projection is defaulted to 2.756 inches, which is the separation used in the standard stereoscopic cameras and viewers (70mm). It is recommended that the default value be used. FILMPAPER CAMERA PAPR [, BLANK FRAMES - BOTH ) (n) (microfilm plotters only) This card offers three options of different cameras or combinations: FILM - 35mm or 16mn film (positive or negative images) PAPER - positive prints B0TH - positive prints and 35mm or 16mm film 4.2-11 (1/15/81)

PLOTTING The request for a 35mm or 16mm camera and positive or negative images must be communicated to the plotter operator through normal means of communications at the installation. Insertion of blank frames between plots is optional and is applicable only to plots generated on film. The type option must be FILM or B0TH if blank frames are desired. The plotter must be operated in the manual mode in order to have blank frames inserted between positive prints. If blank frames are desired only on film, and not on paper, the plotter must be operated in the automatic mode. The default values are PAPER and n = 0. This card is completely optional. PAPER SIZE [ TYPE VELLUM ] ( ~a ) (Dl b )name (SC plotters) 20.0 X 20.0 VELLUM PAPER SIZE {2.}{ }20 [ TYPE nELaM} (NASTRAN plotters) a = horizontal size of paper in inches b = vertical size of paper in inches name = any BCD value desired by user for identification purposes. The default parameter for SC plotters is 7.5 x 7.5, type VELLUM. The default paper size for the NASTRAN plotter is 20 x 20 inches and this is converted to a 7 x 7 inch plot frame by the NASPL0T postprocessor.. The PAPER SIZE parameter can be used along with the NASPL0T postprocessor to create rectangular plots. For example, the card will result in a 14 x 7 inch plot frame if the default value of 1.0 is used for the SCALE FACT0R on the NASPL0T data card. The SCALE FACT0R on the NASPL0T data card can be used to make larger plots having the shape defined with PAPER SIZE. The PAPER SIZE parameter is also used to control the raster count for the NASTRAN plotter. The default raster count is 1000 for a paper size of 20 x 20. Doubling the paper size to 40 x 40 will double the raster count to 2000. 4.2-12 (1/15/81)

STRUCTURE PLOTTING PEN i [COLOR name] i = pen designation number (1 to 4) name = color desired This card generates a message on the printed output which may be used to inform the plotter operator as to what size and which color pen point to mount in various pen holders. The actual number of pens available will depend on the plotter hardware configuration at each installation. This card does not control the pen used in generating the plot (see the PEN option on the PL0T execution card in Section 4.2.2.3). The PEN card is optional, and is not appropriate for microfilm plotters. PTITLE = { tblanks plot title text The plot title text on this card Is printed at the top of the plot on the line below the sequence number. The complete text must be on a single physical card. The default value for the text is all blanks. FIND [SCALE], [ORIGIN i],[VANTAGE POINT],[SET j], REGION 0.0,0.0,1.0,1.0 le, be, re, te i = origin identification number (any positive integer) j = set identification number (any positive integer) le = fractional distance of left edge of plot region from the lower left corner of the image area (default = 0.0) be = fractional distance of bottom edge of plot region from the lower left corner of the image area (default = 0.0) re = fractional distance of right edge of plot region from the lower left corner of the image area (default = 1.0) te = fractional distance of top edge of plot region from the lower left corner of the image area (default = 1.0) 4.2-13 (1/15/81)

PLOTTING The FIND card requests the structure plotter to compute any of the parameters SCALE, 0RIGIN1, and/or VANTAGE P0INT indicated by the user based on 1) the plotter requested on the PL0TTER card, 2) the projection requested on the PR0JECTI0N card, 3) SETJ and REGI0N le, be, re, te requested on the FIND card, 4) the orientation requested on the VIEW and/or AXES card(s), 5) the deformation scaling requested on the MAXIMUM DEF0RMATI0N card, and 6) the paper size as requested on the PAPER SIZE card. All dependencies on which a FIND card is based must precede the FIND card. It is recommended that the FIND card be placed just before the PL0T card or the C0NT0UR card, if present. Any one, two or all three parameters may be computed by the program by using this card, provided that the parameters not requested have already been defined. If no set is specified on this card, the first set defined is used by default. If no options are specified on the FIND card, a SCALE and VANTAGE POINT are selected and ORIGIN 1 is located, using the first defined SET, so that the plotter object is located within the image area. The plot region is defined as some fraction of the image area (image area = 0.,, 0., 1. and first quadrant -.5,.5, 1., 1.). The image area is located inside the margins on the paper. Each FIND card must be one (1) logical card. The FIND card is recommended for general use. MAJPRIN MINPRIN MAXSHEAR XN0RMAL (EVEN 10 Z1 ) C0NT0UR YN0RMAL EVEN n 77 C0MM0N ZN0RMAL LIST a,b,... MAX L0 XYSHEAR MIDJ XZSHEAR YZSHEAR (stress contour plots only) 4.2-14 (1/15/81)

STRUCTURE PLOTTING 1. Type of Stress MAJPRIN - Major principal stress (default) MINPRIN - Minor principal stress MAXSHEAR - Maximum shear stress XNORMAL) YN0RMAL - X, Y, Z components of normal stress ZNORMAL XYSHEAR) XZSHEAR - Components of shear stress YZSHEAR 2. Number of Stress Contours EVEN - Integer following is number of stress contours (default is 10, maximum is 50). LIST - Real numbers following are values of stress contours. 3. Location of Stress (Z1 and Z2 specified on property card) Z1 - Distance Z1 from neutral plane Z2 - Distance Z2 from neutral plane MAX - Maximum of stress at Z1 and Z2 MID - Average of stress (membrane stress) at Z1 and Z2 4. Coordinate System C0MM0N - Plot stress contours in basic coordinate system (default) L0CAL - Plot stress contours in element coordinate system Stress contour plots are available for the following elements: TRIA3, QUAD4, SHEAR, and TRIAX6. The connection cards for all elements must list the grid points in either clockwise or counterclockwise order. Some in each order will result in meaningless plots. When selecting contour options, note: 1. MAJPRIN, MINPRIN, MAXSHEAR are the same in COMMON and LOCAL. 2. ZNORMAL, XZSHEAR, YZSHEAR, if selected in LOCAL, will be changed to COMMON. 3. SHEAR elements only have the MAXSHEAR value. 4.2-15 (1/15/81)

PLOTTING The TRIAX6 stress contour plots are different in that they must be selected as C0MM0N, then the following equivalence table applies: XN0RMAL = RADIAL YN0RMAL = AZIMUTH ZNORMAL = AXIAL XYSHEAR = SHEAR XZSHEAR MAX PRINCIPAL YZSHEAR = MAX SHEAR MAXSHEAR = OCTAHEDRAL The CONTOUR card should be placed immediately before the associated PL0T execution card in order to be certain that all applicable parameters have been previously defined. A STRESS request must appear in the Case Control Deck for all elements included in a C0NT0UR request. If printed output is not desired, the request can be of the form STRESS(PL0T) = { ALL4 C0NT0UR ZD I VSP }V MAGNITSP LIST a,b,... (displacement or temperature contour plots only) 1. Type of Displacement XDISP YDISP X,Y,Z components of displacement in global coordinate system ZDISP ) MAGNIT Magnitude of deformation, also used for contours of temperatures 2. Number of Displacement Contours EVEN - Integer following is number of displacement contours (default is 10, maximum is 50). LIST - Real numbers following are values of displacement contours. 4.2-16 (1/15/81)

STRUCTURE PLOTTING 4.2.2.3 PLOT Execution Card STATIC MODAL I (DEFORMATION RANGE fl, f2 PLOT ICMODAL I VEL0CITY f I({CONTOUR} [il, i2 THRU 13, 14. etc.] RANGE X1, X2 I TRANSIENTI L-ACCELERATION (TIME tl, t2 / FREQUENCY l blank {PHASE LAG i MAGNITUDE [MAXIMUM DEFORMATION d] [SET jl] [ORIGIN kl] [ASYMMETRY [ENSTY] [SYMB0LS m [ n]] GRID POINTS ELEMENTS SHAPE LABEL BOTH UTLIN [VECTOR v], [PRINT] EPID - GSPC [SHRINK t,o] [NORMALS] [ SET j2] [0RIGIN k2]..., etc. This logical card will cause one picture to be generated for each subcase, mode or time step requested using the current parameter values. If only the word PLOT appears on the card, a picture of the undeformnned structure will be prepared using the first defined set- and the first defined origin. Options 1 thru 6, if used, must be in the order shown. The available plot options and their meanings are: 1. STATIC - Plot static deformations MODAL - Plot mode shapes CM0DAL - Plot complex mode shapes TRANSIENT - Plot transient solutions FREQUENCY - Plot frequency solutions blank - No entry is a request for either undeformed plots or undeformed underlay for contour plots 4.2-17 (1/15/81)

PLOTTING 2. DEF0RMATI0N - Plot displacements or temperatures (in the Z direction) (default value) VEL0CITY - Plot velocities ACCELERATI0N - Plot accelerations 3. C0NT0UR - Request for contour plots 4. 1, i2,... - Nonzero integers following refer to subcase numbers that are to be plotted. Default is all subcases. See SHAPE and VECTOR for use of "O" (underlay) command. In eigenvalue problems, the list of subcases must refer to subcase IDs whenever the number of modes recovered is equal to or less than the number of subcases defined. If the number of modes recovered is more than the subcases defined, the plot request for those modes associated with the subcases must refer to subcase IDs. After the mode associated with the last defined subcase, higher modes will be identified by incrementing the last defined subcase ID by one for each such higher mode. In eigenvalue problems using Type II cyclic symmetry solution sequences (Rigid Format 48 and buckling analysis solution sequence DMAP5CY) the plot requests for segments of the model must refer to the coded subcase IDs (see Section 1.11 for definition of coded subcase IDs). All eigenvectors recovered for the segment will be plotted. The RANGE option can be used to select a subset of all eigenvectors for plotting. 5. RANGE - Refers to range of eigenvalues (buckling) or frequencies using requested subcases, for which plots will be prepared. TIME - Refers to time interval, using requested subcases and output time steps, for which plots will be prepared. 6. PHASE LAG - Real number, +, in degrees (default is 0.0). The plotted value is UR cos + U sin *, where uR and uI are the real and imaginary parts of the response quantity. MAGNITUDE - plotted value. is u2 + U. 7. MAXIMUM DEFORMATION - Real number following is used instead of the actual maximum displacement component in scaling the displacements for all subcases. This option is not recommended for general use. Each subcase is separately scaled according to its own maximum if this item is absent. If d is omitted, the set will be scaled to the maximum within the set being plotted. 8. SET - Integer following identifies a set which defines the portion of the structure to be plotted. Default is first set defined. 9. ORIGIN - Integer following identifies the origin to be used for the plot. Default is first origin defined. 10. SYMMETRY w - Prepare an undeformed or deformed plot of the symmetric portion of the object which is defined by SET j. This symmetric portion will be located in the space adjacent to the region originally defined by ORIGIN k, and will appear as a reflection of the antisymmetrically deformed structure about the plane whose normal is oriented parallel to the coordinate direction w. 4.2-18 (1/15/81)

STRUCTURE PLOTTING ANTISYMMETRY w - Prepare a deformed plot of the symmetric portion of the antisymmetrically loaded object which is defined by SET j. This symmetric portion will be located in the space adjacent to the region originally defined by ORIGIN k, and will appear as a reflection of the antisymmetrically deformed structure about the plane whose normal is oriented parallel to the coordinated direction w. The symbol w may specify the basic coordinates X, Y or Z or any combination thereof. This option allows the plotting of symmetric and/or antisymmetric combinations, provided that an origin is selected for the portion of the structure defined by the bulk data that allows sufficient room for the complete plot. This does not permit the combination of symmetric and antisymmetric subcases, as each plot must represent a single subcase. In the case of a double reflection, the figure will appear as one reflected about the plane whose normal is parallel to the first of the coordinates w, followed by a reflection about the plane whose normal is oriented parallel to the second of the coordinates w. This capability is primarily used in the plotting of structures that are loaded in a symmetric or an antisymmetric manner. The plane of symmetry must be one of the basic coordinate planes. 11. PEN - Integer following controls the internal NASTRAN pen number that is used to generate the plot. DENSITY - Integer following specifies line density for film plotters. A line density of d is d times heavier than a line density of 1. 12. SYMBOLS m[,n] - All the grid points associated with the specified set will have syibol m overprinted with symbol n printed at its location. If n is not specified, only symbol m will be printed. Grid points excluded from the set will not have a symbol. Grid points in an undeformed underlay will be identified with symbol 2. Following is a table of symbols available on each plotter. Symbols that are not available on a given plotter are defaulted to a similar symbol indicated in parentheses. SYMBOL AVAILABILITY NO. m or n SYMBOL SC4020 NAST 0 no symbol X X 1 X X X 2 * X X 3 + X X 4 X X 5 X X 6 6 X X 7 D X X 8 (7) X 9 / (7) X 13. LABEL GRID - All the grid points included in the specified set have their i-ntftification number printed to the right of the undeflected or deflected location (undeflected location in the case of superimposed plots). 4.2-19 (1/15/81)

PLOTTING LABEL ELEMENTS - All the elements included in the specified set are identified by the element identification number and type at the center of each element (undeflected location in the case of superimposed plots). LABEL B0TH - Label both the grid points and elements. LABEL GSPC - Label those degrees of freedom that are constrained to zero through permanent single point constraints on GRID and GRDSET bulk data cards or are constrained through SPC and SPC1 bulk data cards. The label consists of the grid point ID number and the constrained degrees of freedom. For superelements, this type of plot is available only after restarting from a run which completes the SEMG operation. For non-superelements, RFxx020 is required. LABEL EPID - Label elements with their respective property card identification (PID) numbers. The label consists of the standard element labels and element PID. For non-superelements, RFxxD2O is required. 14. SHAPE - All the elements included in the specified set are shown by connecting the associated grid points in a predetermined manner. OUTLINE - Only the outline of all the elements in the specified set are shown. Elements not supported by contour plots are ignored. Outlines are always drawn with PEN 1. Both deformed and undeformed shapes or outlines may be specified. All the deformed shapes relating to the subcases listed may be underlaid on each of their plots by including "0" with the subcase string on the PL0T card. The undeformed plot will be drawn using PEN 1 or DENSITY 1 and symbol 2 (if SYMB0LS is specified). Omitting both on contour requests gives undeformed shape. 15. VECTOR v - A line will be plotted at the grid points of the set, representing in length and direction the deformation of the point. Vectors representing the total deformation or its principal components may be plotted by insertion of the proper letter(s) for variable v. Possible vector combinations are: X or Y or Z - requesting individual components XY or XZ or YZ - requesting two specified components XYZ - requesting all three components RXY or RXZ or RYZ - requesting vector sum of two components R - requesting total vector deformation N - used with any of the above combinations to request no underlay shape be drawn All plots requesting the VECTOR option will have an underlay generated of the undeformed shape using the same sets, PEN 1 or DENSITY 1, and symbol 2 (if SYMB0LS is specified). If "SHAPE" and "VECT0R" are specified, the underlay will depend on whether "0" is used with DEF0RMATI0N. It will be the deformed shape when not used and will be both deformed and undeformed shapes when it is used. The part of the vector at the grid point will be the tail when the underlay is undeformed and the head when it is deformed. If the "N" parameter is used, no shape will be drawn but other options such as SYMB0LS will still be valid. 4.2-20 (1/15/81)

STRUCTURE PLOTTING 16. PRINT - A list of the average stresses at the interior grid points in the set will be printed for contour stress plots 17. SHRINK - The real number (t) following is used to shrink two-dimensional elements with three or four grid points (0 < t < 1.0). If t is omitted, the default value is 0.1 which results in a 10% reduction. The second real number (o) is used to shrink onedimensional elements. If o is omitted, the one-dimensional elements are not reduced. Both t and o must be nonzero in order to shrink one-dimensional elements. QUAD8 and TRIA6 do not shrink, regardless of number of deleted midside nodes. 18. NORMALS - Plot vector normal to HBDY elements Examples of PLOT Cards 1. PLOT Undeformed SHAPE using first defined SET, first defined ORIGIN and PEN 1 (or DENSITY 1). 2. PLOT SET 3 ORIGIN 4 PEN 2 SHAPE SYMBOLS 3 LABEL Undeformed SHAPE using SET 3, ORIGIN 4, PEN 2 (or DENSITY 2) with each grid point of the set having a + placed at its location, and its identification number printed adjacent to it. 3. PLOT MODAL DEFORMATION 5 SHAPE Modal deformations as defined in subcase 5 using first defined SET, first defined ORIGIN, and PEN 1 (or DENSITY 1). 4. PLOT STATIC DEFORMATI0N 0, 3 THRU 5, 8 PEN 4, SHAPE STATIC deformations as defined in subcases 3, 4, 5 and 8, deformed SHAPE; drawn with PEN 4, using first defined SET and ORIGIN, underlayed with undeformed SHAPE drawn with PEN 1. This command will cause for plots to be generated. 5. PLOT STATIC DEFORMATION 0 THRU 5, SET 2 ORIGIN 3 PEN 3 SHAPE, SET 2 ORIGIN 4 PEN 4 VECTORS XYZ SYMBOLS O, SET 35 SHAPE Deformations as defined in subcases 1, 2, 3, 4, and 5, undeformed underlay with PEN 1, consisting of SET 2 at ORIGIN 3, SET 2 at ORIGIN 4 (with an * placed at each grid point location), and SET 35 at ORIGIN 4. Deflected data as follows: SHAPE using SET 2 at ORIGIN 3 (PEN 3) and SET 35 at ORIGIN 4 (PEN 4); 3 VECTORS (X, Y and Z) drawn at each grid point of SET 2 at ORIGIN 4 (PEN 4) (less any excluded grid points), with o placed at the end of each vector. 6. PLOT STATIC DEFORMATIONS 0, 3, 4, SET 1 ORIGIN 2 DENSITY 3 SHAPE, SET 1 SYMMETRY Z SHAPE, SET 2 ORIGIN 3 SHAPE, SET 2 SYMMETRY Z SHAPE 4.2-21 (1/15/81)

PLOTTING Static deformations as defined in subcases 3 and 4, both halves of a problem solved by symmetry using the X-Y principal plane as the plane of symmetry. SET 1 at ORIGIN 2 and. SET 2 at ORIGIN 3, with the deformed shape plotted using DENSITY 3 and the undeformed structure plotted using DENSITY 1. The deformations of the "opposite" half will be plotted to correspond to symmetric loading. This command will cause two plots to be generated. 7. PLOT TRANSIENT DEFORMATI0N 1, TIME 0.1, 0.2, MAXIMUM DEFORMATI0N 2.0, SET 1, ORIGIN 1, PEN 2, SYMBOLS 2, VECTOR R Transient deformations as defined in subcase 1 for time - 0.1 to time - 0.2, using SET 1 at ORIGIN 1. The undeformed SHAPE using PEN or DENSITY 1 with an * at each grid point location will be drawn as an underlay for the resultant deformation vectors using PEN or DENSITY 2 with an * typed at the end of each vector drawn. In addition, a plotted value of 2.0 will be used for the single maximum deformation occurring on any of the plots produced. All other deformations on all other plots will be scaled relative to this single maximum deformation. This command will cause a plot to be.generated for each output time step which lies between 0.1 and 0.2. 8. CONTBUR XNORMAL PLOT CONTOUR, SET 2, ORIGIN 4, OUTLINE Contour plot of x-component of normal stress for elements in SET 2 in basic coordinate system at a distance Z1 from neutral plane with 10 contour lines, an outline of elements in SET 2, and using ORIGIN 4. 9. CBNTOUR MAGNIT, LIST 2., 4., 6., 8., 10 PLOT CBNTOUR, SET 5, OUTLINE Contour plot of magnitude of displacements at grid points associated with: elements in SET 5 with 5 contours having values of 2., 4., 6., 8., 10., and an outline of the elements in SET 5 using ORIGIN 4. 10. PLOT CM0DAL DEFORMATION PHASE LAG 90. SET 1 VECTOR R The imaginary part of the complex mode shape will be plotted for SET 1. 4.2.3 Examples of Structure Plot Requests The BEGIN BULK card is shown in each of the examples to remind the user to place the structure plot request at the end of the Case Control Deck. 4.2-22 (1/15/81)

STRUCTURE PLOTTING Example 1 0UTPUT(PL0T) PL0TTER NAST SET 2 - ALL FIND SCALE, 0RIGIN 5, SET 2 PL0T SET 2, 0RIGIN 5 BEGIN BULK A single undeformed plot using all elements will be produced for the NASTRAN general purpose plotter. Example 2 0UTPUT(PL0T) PL0TTER NAST CSCALE = 1.8 SET 10 = R0D AXES Z,X,Y VIEW 0.0, 0.0, 0.0 FIND SCALE, 0RIGIN 5, SET 10 PL0T SET 10, 0RIGIN 5, LABEL BEGIN BULK A single undeformed plot using all R0D elements will be produced for the NASTRAN general purpose plotter. The grid points will be labeled with the grid point identification numbers. The view will be along the z-axis, or normal to the X-Y plane. 4.2-23 (1/15/81)

PLOTTING Example 3 OUTPUT(PLOT) PLOTTER NAST CSCALE = 1.8 SET 2 = 1 THRU 50 SET 3 = 101 THRU 199 FIND SCALE, ORIGIN 20, SET 2 PLOT STATIC DEFORMATION SET 2, ORIGIN 20 FIND SCALE, ORIGIN 25, SET 3 PLOT STATIC DEF0RMATION 0, 2, SET 3, ORIGIN 25 PTITLE - THERMAL LOAD PLOT STATIC OEFORMATION 3 THRU 10, SET 3, ORIGIN 25, LABEL BEGIN BULK Several deformed plots will be produced for the NASTRAN plotter. The maximum deformation will be scaled to the default value of 5% of the maximum dimension of the model. The first PL0T command card will produce one deformed plot for each subcase in the Case Control Deck. All elements having identification numbers in the range of 1 through 50 will be included on the plot. All plots will be located at 0RIGIN 20, as determined by the first FIND card. The second PL0T command card will produce one deformed plot for subcase 2 with an underlay of the undeformed structure. All elements having identification numbers in the range of 101 through 199 will be included on the plot. The plot will be located at ORIGIN 25, as determined by the second FIND card. The third PL0T command card will produce one deformed plot for each subcase having an identification number in the range 3 through 10. All elements having identification numbers in the range of 101 through 199 will be included on the plot. The grid points will be labeled with the grid point identification numbers. The plot will be located at ORIGIN 25, as determined by the second FIND card. 4.2-24 (1/15/81)

STRUCTURE PLOTTING 4.2.4 Summary of Structure Plot Request Packet Cards SET Definition - Required SET i [INCLUDE][ELEMENTS], J2' J3 THRU j4, jg, etc INEXCLUDE EXCiEPT GRID PINTS k1 k2, k3 THRU k, k5, etc. Parameter Definition - Optional, except as noted PLOTTER { NASTC} J (required if not NAST) (ORTHOGRAPHIC) PERSPECTIVE STEREOSCOPIC | AXES R, S, T TISYMTR YVIEW y, B. (| SCALE a[, b] (Required if not on FIND card) 0 RIGIN i u, v I (Required if not on FIND card) VANTAGE POINT ro,so,to,do[,sor] (Required for perspective and stereoscopic _1 ----— 0 —-— I ~projections if not on FIND card) PR0JECTION PLANE SEPARATI0N {2-O OCULAR SEPARATION 2.756 in } MAXIMUM DEF0RMATION d PEN i [C0L0R name] 4.2-25 (1/15/81)

PLOTTING FILM CAMERA )PAPER. iBLANK FRAMES i.. (BOTH ln) PAPER SIZE {^} ( XB!75 [, TYPE ( ^^j_ (SC plotters) ( a ) (BY) ( b ) [ ( name PAPER SIZE 20 20.0, TYP VELL (NASTRAN plotters) a BY ( name blanks PTITLE {. blanks -t } PTIT E = p'plot title text DIST0RTI0N x,y,z FIND [SCALE],[0RIGIN i],[VANTAGE POINT],[SET j],[REGI0N le, be, re, te] MAJPRIN MINPRIN MAXSHEAR XNORMAL EVEN 10 (Z C0NT0UR YN0RMAL EVEN n 7T GMMON ZN0RMAL (LIST a,b,...) MAX CAL XYSHEAR. MID XZSHEAR YZSHEAR.E.,A (stress contour plots only) (XDISP EVEN 1 (EVEN 10 ) CONTOUR iXD ISP V n ZDISP LIST a,b,... MAGNIT) (displacement contour plots only) 4.2-26 (1/15/81)

STRUCTURE PLOTTING PLOT Execution Card - Required STATIC 1 LMODAL 1 (DEF0RMATION! RANGE fl, f2) PLOT CMODAL V I VELOCITY C {CONT0UR} [il, 12 THRU 13, 14, etc.] IRANGE X1, X2I TRANSIENT (ACCELERATION (TIME tl, t2 FREQUENCY L J L blank J PHASE LAG )1 [MAXIMUM DEFORMATION d] MAGNITUDE [SET j1] [ORIGIN kl] [{ASYMM ETRY} w PENSiTY P] [SYMOOLS m[,n]] GRID POINTS ELEE M BN T H SHAPE ELABEL BT [HOU}TLINE'[ CVECTOR v] [PRINT] GSPC j [SHRINK t,o], [NORMALS] [SET j2] [ORIGIN k2]..., etc. 4.2-27 (1/15/81)

PLOTTING (a) y - rotation about T-axis. T z R (b) S rotation about S-axis T tt R (c) a - rotation about R-axis Figure 1. Plotter Coordinate System-model Orientation. 4.2-28 (1/15/81)

DIRECT MATRIX ABSTRACTION 5.1 GENERAL In addition to using the solution sequences provided by MSC/NASTRAN, the user may wish to execute a series of modules in a different manner than provided by a Solution Sequence. Or, he may wish to perform a series of matrix operations which are not contained in any existing Solution Sequence. If the modifications to an existing solution sequence are minor, the ALTER feature described in Sections 2 and 5.2.4 may be employed. Otherwise, a user-written Direct Matrix Abstraction Program (DMAP) should be used. DMAP is the Data Block-oriented language used by MSC/NASTRAN to solve problems. A solution sequence is basically a collection of statements in this language. DMAP, like English or FORTRAN, has many grammatical rules which must be followed to be interpretable by the MSC/NASTRAN DMAP compiler. Section 5.2 provides the user with the rules of DMAP which will allow him to understand the solution sequence DMAP sequences, write ALTER packages, and construct his own DMAP sequences using the many modules contained in the DMAP repertoire. Section 5.3 provides general examples of DMAP usage. Section 5.4 describes individually the many DMAP programmer-oriented modules contained in the DMAP library. Before proceeding with the specific rules of DMAP program creation, a few general comments on the nature of the language and its relation to the MSC/NASTRAN program are in order. DMAP began as a spur capability of a structurally-oriented computer program. Its purpose was to allow the user of MSC/NASTRAN access to the matrix routines which would support the structural solution. However, DMAP grew as a language until today it dominates MSC/NASTRAN. In fact, all operations undertaken by MSC/NASTRAN are specified via DMAP. A DMAP programmer has access to all parts of the program. The MSC/NASTRAN program is separated into functional blocks called modules each of which has a name (module name) and a specific function to perform (purpose). Modules are not orderconstrained in that the DMAP programmer can schedule any module order that is meaningful to him. Modules are constrained in their methods of communication with subsequent modules. In fact, one module can communicate with another module only through the MSC/NASTRAN Executive System (NES). In particular, modules must write their answers on disk space provided by the NES. These logical collections of data are given names and called data blocks. A good example of a data block is a matrix. The DMAP language then expresses a logical relation between data blocks as operated on by 5.1-1 (10/1/80)

DIRECT MATRIX ABSTRACTION modules. Data block names are arbitrary; the content of the data block is important, not its name. A typical DMAP program might appear as follows: BEGIN $ ADD A,B/C $ MATPRN C//$ END $ This program adds the two matrices (data blocks) A and B together forming C, which is then printed. Four modules are invoked, BEGIN, ADD, MATPRN and END. If C were changed to ANSWER in both its occurrences, the DMAP program would be unchanged. In any language another class of objects is needed. These control the flow and communicate options to operations taking place. In MSC/NASTRAN these scalar values are called parameters. Most modules provide for one or more parameters to express options. Thus, the add matrix routine allows a scalar multiplier for each matrix. Parameters, however, may have default values, and the example program used them. If the DMAP programmer wanted to multiply [A] by 5.0, he could code ADD A,B/C/(5.O,O.O) $ A DMAP module then reads its input data blocks (A,B), operates as specified by its parameters ((5.0,0.0)), and writes its outputs (C). Meanwhile, the NES is establishing and controlling the sequence of module executions (according to the DMAP program), establishing, protecting and communicating values of parameters, allocating disk space (files) to data blocks, and potentially providing a restart capability. To learn to write a DMAP program, the user must accomplish a number of steps. He must learn the basic syntax rules of the DMAP compiler and how to correctly specify the inputs and outputs to a module (Section 5.2.1). He must then learn to write sequences of modules which properly link together (Section 5.2.2). As his sequences become more complex, he will want to invoke special functions of the NES (Section 5.2.3). Examples of DMAP programs are provided in Section 5.3 as well as by the solution sequences (which can be printed by including DIAG 14 in the Executive Control Deck). 5.1-2 (10/1/80)

GENERAL The DMAP programmer must also expand his DMAP vocabulary. Section 5.4 provides a list of commonly used modules which may serve as a beginning. Of course, all modules may be used in a DMAP program. A list of all modules can be printed out by DIAG 31. Descriptions of all modules are found in the MSC/NASTRAN Programmer's Manual, Section 4. A particularly valuable use of DMAP is to modify one of the existing DMAP programs (Solution Sequences) via the ALTER package to provide for a specialized application. New capability, if expressed through new modules, is only available through the use of DMAP ALTERs. Many prewritten ALTERs are available through the MSC/NASTRAN RFALTER Library. Section 5.2.4 will discuss the ALTER process in some detail. 5.1-3 (7/28/80)

OMAP MODULE DESCRIPTIONS I. NAME: MATPRN (General Matrix Printer) II. PURPOSE: To print general matrix data block. III.. DMAP CALLING SEQUENCE: MATPRN M,M2,M3,M4,M5// $ IV. INPUT DATA BLOCKS: Mi - Matrix data blocks, any of which may be purged. (See Reniark 2). V. OUTPUT DATA BLOCKS: None VI. PARAMETERS: None VII I. OUTPUT: The nonzero band of each column of each input matrix data block is printed (see Remark 4). VIII. REMARKS: 1. Any or all input data blocks may be purged. 2. If any data block is not a matrix, it will be printed as if it were a table. 3. MATPRN prints the row index for the term which begins each line of printout. 4. MATPRN will not printout two or more consecutive lines of zeroes, but instead will issue a message of the form: ROW P0SITIONS xxxx THRU yyyy NOT PRINTED - ALL - 0.0. 5. If DIAG 30 is set by the PARAM module before MATPRN (see Example 3), and turned off after MATPRN, most of the digits of the internal representation will be printed. Normally, the output is truncated to five or six digits. 6. For large, sparse matrices with scattered terms, the user is advised to use either the MATPRT or MATGPR modules. IX. EXAMPLES: 1. MATPRN KGG// $ 2. MATPRN KGG,PL,PG,BGG,UPV// $ 3. PARAM //OIAG0N // 30 $ PRINT EXTENDED PRECISI0N MATPRN KGG// S PARAM // DIAG0FF // 30 $ 5.4-63 (5-15-80)

DMAP MODULE OESCRIPTIONS I. NAME: SOLVE (Linear System Solver) II. PURPOSE: To solve the matrix equation [A][X] +~[B] III. OMAP CALLING SEQUENCE: SOLVE A,B/X/V,Y,SYM/V,Y,SIGN $ IV. INPUT DATA BLOCKS: A - square real or complex matrix B - rectangular real or complex matrix (if purged, the identity matrix is assumed) V. OUTPUT OATA BLOCKS: X - a rectangular matrix Note: A standard matrix trailer will be written, identifying [X] as a rectangular matrix witn the same dimensions as [B] and the type specified. VI. PARAMETERS: 0 - use decomposition consistent with form of [A] SYM - Input-integer-default = 0 1 - use symmetric decomposition -1 - use unsymmetric decomposition SIGN - Input-integer-default = 1 1 - solve [A][X] = [] VII. METHOD: Depending on the SYM flag and the type of [A], one of subroutines SDC0MP or UDCjMP is called to form rA] [ CLICU]. One of FBS or UFBS is then called to solve [L]C[Y + [B] and [L][X] = CY], as appropriate. VIII. EXAMPLES: 1. Solve a system of equations CA][U] - [P] S0LVE A,P/U/ S 2. Invert [A] SOLVE A,/AINV/ $ 5.4-121 (5-15-80)

UNIVERSITY OF MICHIGAN 3 9015 02499 5584 THE UNIVERSITY OF MICHIGAN DATE DUE X? /5,//0 o Q 4s''/ i,4 ~o