NASA/CR--2003-212127

User's Guide for the

Software System

Ashwin Shah

Sest, Inc., Middleburg Heights, Ohio

Updated EST/BEST

February 2003
The NASA STI Program Office... in Profile

Since its founding, NASA has been dedicated to

the advancement of aeronautics and space

science. Tile NASA Scientific and Technical

Information (STI) Program Office plays a key part

in helping NASA maintain this important role.

Tile NASA STI Program Office is operated by

Langley Research Center, the Lead Center for

NASA's scientific and technical information. The

NASA STI Program Office provides access to the

NASA STI Database, the largest collection of

aeronautical and space science STI in file world.

The Program Office is also NASA's institutional

medlanism for disseminating the results of its

researd3 and development acti vities. These results

are published by NASA in the NASA STI Report

Series, which includes the following report types:

TECHNICAL PUBHCATION. Reports of

completed :research or a major significant

phase of research that present the results of

NASA programs and include extensive data

or theoretical analysis. Includes compilations

of significant scientific and technical data and

information deemed to be of continuing

reference value. NASA's counterpart of peerreviewed formal professional papers but

has less stringent limitations on manuscript

lengfl3 and extent of graphic presentations.

TECHNICAL MEMORANDUM. Scientific

and tedmical findings that are pre, liminary or

of specialized interest, e.g., quick release

reports, working papers, and bibliographies

that contain minimal annotation. Does not

contain extensive analysis.

CONTRACTOR REPORT. Scientific and

technical findings by NASA-sponsored

contractors and grantees.

CONFERENCE PUBLICATION. Collected

papers from scientific and technical

conferences, symposia, seminars, or other

meetings sponsored or cosponsored by

NASA.

SPECIAL PUBLICATION. Scientific,

technical, or historical information from

NASA programs, projects, and missions,

often concerned with subjects having

substantial public interest.

TECHNICAL TRANSLATION. Englishlanguage translations of foreign scientific

and technical material pertinent to NASA's

mission.

Specialized services that complement the STI

Program Office's diverse offerings incl ude

creating custom thesauri, building customized

databases, organizing and publishing research

results.., even providing vid eos.

For more information about the NASA STI

Program Office, see the following:

® Access the NASASTI Program Home Page

at http:lhuww.sti.nasa.gov

® E-mail your question via the Intemet to

help@sti.nasa.gov

* Fax your question to the NASA Access

Help Desk at 301-621-0134

* Telephone the NASA Access Help Desk at

301-621-0390

Write to:

NASA Access Help Desk

NASA Center for AeroSpace Information

7121 Standard Drive

Hanover, MD 21076
NASA/CR--2003-212127

User's Guide for the

Software System

Ashwin Shah

Sest, Inc., Middleburg Heights, Ohio

Updated EST/BEST

Prepared under Contract C.----76751----.J

National Aeronautics and

Spa ce Ad.minis tration

Glelm Research Center

February 2003
The Propulsion and Power Program at

NASA Glenn Research Center sponsored this work°

NASA Center for Aerospace Information

71121Standard Drive

Hanover, MD 211076

Available frorn

National Technical Information Service

5285 Port Royal Road

Springfield, VA 22100

Available electronically at http: //gltrs.grc.nasa.gov
USER'S GUIDE TO THE UPDATED EST/BEST SOFTWARE SYSTEM

A.1 Interactive Input General Information

The structure of the IPACS input file reflects the modularity of each module. The

structured format helps the user locate specific input data and manually enter or edit it. The

IPACS input file can have any user-specified filename, but must have a .DAT extension. The

input file may consist of up to six input data blocks; the data blocks must be separated by

delimiters beginning with the $ character. If multiple sections are desired, they must be arranged

in the following order:

$COB

- Automatic Mesh Generation Input -

$PIC

- Material Property Uncertainty Input -

$NES

- Finite Element Input -

$RAN

- Structure Level Random Variable Input -

$FPI

- Probabilistic Output Request Input -

$EXE2

-Execution Control Input

The delimiters may be abbreviated to the first four characters (i.e. $COB, $PIC, $NES, $FPI).

NASA/CR--2003-212127 1
Thestructureof theNESSUSinputfile is sameasthatisdescribedin theNESSUS

manual.

During interactive processing of NESSUS input blocks, scratch files may also be created

with the following extensions and formatting:

IPACS/NESSUS Body Force Block Card 1:

inputfilename.BODY Card 2 (IBODY=I):

Card 2 (IBODY=2):

Card 3 (IBODY=2):

*BODY, IBODY(I)

X-COMP(F), Y-COMP(F), Z-COMP(F)

XI(F), YI(F), ZI(F)

X2(F), Y2(F), Z2(F)

IPACS/NESSUS Boundary Conditions Block

inputfilename.BOUND Cards 1 thru N: NODE #(I), DOF #(I), DISP(F)

NESSUS Material Properties Block

inputfilename.prop Cards 1 thru N: NODE #(I), NODE(J) ,(PROP(I),I=I,5)

NESSUS Laver Temperature Block

inputfilename.temp Cards 1 thru N: NODE #(I), NODE(J) ,(TEMP(I),I=I,5)

IPACS/NESSUS Element Connectivity Block

input filename. CONN Cards 1 thru N: ELEM #(I), INODE(I), JNODE(I), KNODE(I),

NODE(I)

IPACS/NESSUS Nodal Coordinates Block

input filename. COORD Cards 1 thru N: NODE #(I), X(F), Y(F), Z(F), TH(F)

IPACS/NESSUS Random Variable Coordinates Block

inputfilename.CRDxxx Card 1: MEAN(F), STD DEV(F), DIST TYPE(I)

Cards 2 thru N: NODE #(I), X-SDF(F), Y-SDF(F), Z-SDF(F),

TH-SDF(F)

IPACS/NESSUS Distributed Loads Block

input filename.DIST Cards 1 thru N: BEG ELEM(I), END ELEM(I), ICODE(I),

TRAC(F), PRESS(F), BODY FORCE(F)

IPACS/NESSUS Random Variable Distributed Loads Block

inputfilename.DSTxxx Card 1: MEAN(F), STD DEV(F), DIST TYPE(I)

Cards 2 thru N: BEG ELEM(I), END ELEM(I), ICODE(I),

TRAC SDF(F), PRESS SDF(F), BODY

FORCE SDF(F)

IPACS/NESSUS Duplicate Nodes Block

inputfilename.DUPL Cards 1 thru N: MASTER NODE(I), SLAVE NODE(I)

NASA/CR--2003-212127 2
IPACS/NESSUS Nodal Forces Block

inputfilename.FORCE Cards 1 thru N: NODE #(I), DOF #(I), FORCE(F)

IPACS/NESSUS Random Variable Forces Block

input filename.FRCxxx Card 1:

Cards 2 thru N:

MEAN(F), STD DEV(F), DIST TYPE(I)

NODE #(I), DOF #(I), FORCE SDF(F)

IPACS Material Orientation Block

inputfilename.orie Cards 1 thru N: BEG NODE(I), END NODE(I), X-COMP(F),

Y-COMP(F), Z-COMP(F)

where xxx is a three-digit number identifying the random variable and where (I) indicates integer

and (F) indicates floating point numbers.

The IPACS/NESSUS scratch files are useful in recovering data following a program abort

or a system failure. To recover data, use the "Input Existing File" options found on many of the

menus, after reentering the input module using a filename different than the one used in the

previous run. The quantities found in each of these files are written in free format consistent

with the IPACS/NESSUS convention. This feature also allows users to import their own mesh,

boundary conditions, and loads from outside the input module. To import this data, use the

scratch file naming convention described above and simply have the data read into the program.

NASA/CR--2003-212127 3
A.2 Detailed Description of the Interactive Input Module

As described in section 3.1, the IPACS/NESSUS input file consists of five separate

sections including automatic mesh generation module, a module to compute probabilistic

composite material properties, probabilistic finite element analysis module, and probability

algorithm module.

The IPACS code contains an interactive input module that allows the user to interactively

create or edit an IPACS input file and submit it for processing. Interactive module enables an

inexperienced user to prepare the input in a user-friendly manner in a short time without making

common mistakes. This module allows creating or editing the input file with little or no

knowledge of the file's structure and format, and thus allows an inexperienced user to run the

code and come up to speed more quickly.

The interactive input module is designed for ease of use. An input session need not be

completed in one sitting: once an IPACS input file has been created it can be saved and returned

to later for further editing, or in order to enter the execution mode and run the code. (For

example, the user can enter or change the finite element mesh inputs, save the entries or edits,

exit the IPACS program, and return later to add or edit the material specifications.)

An input session is begun at the main menu, which offers six options that allow the user

quick access to input functions (refer to Figure 3.2. lofthe IPACS manual). The six input

options may be selected and run in any order. The options access a series of menus in which the

user is prompted for information to enter or change.

A.3 Batch Input

IPACS input can be prepared manually or by using any other preprocessor and processing

it in a batch mode. The procedure to prepare the input is described in the following section.

As described in section A. 1, the input file is divided into several different blocks. The

block pertaining to the automatic mesh generation data should not be used in the batch mode.

Each of the material property uncertainty input and finite element input data blocks are divided

into parameter and model data groups. The parameter group of a data block defines the size of

the problem, whereas the model group of the data block contains the actual model data. In both

the data blocks, the model group always follows the parameter group. Also, a 'END' card marks

end of the parameter group in a data block. The $PIC input data block contains the data related to

the material properties and the lay-up configuration of composites in the structure. The $NES

input data block contains the actual finite element data. Detailed input information on each block

is discussed below:

A.3.1 Material Property Uncertainty Input Data Block:

The material property uncertainty input data block begins with a '$PIC' card. The

composite lay-up configuration and material property uncertainties are defined in this block.

NASA/CR--2003-212127 4
Normally,the$PICdatablockendswiththebeginningof the$NESdatablock(i.e.,with a

'$NESSUS'card).Differentkeywords(akeywordalwaysbeginswith an'*'; asterisk)usedin

thematerialpropertyuncertaintydatablockarediscussedbelow:

A.3.1.1 Material Property Uncertainty Parameter Input Data Group:

CFEM

The above card is used to flag ifemx (MV) analysis

$PICn

where n = 0 is used if material properties PDF computation is desired

= 1 is used material properties PDF computation is to be bypassed

*PROSOL

The keyword PROSOL is used to define the probabilistic structural analysis solution

method. The general format is

*PROSOL

NSOL

NMCS

where NSOL is the probabilistic solution method. Different available probabilistic

solution methods and their NSOL values are given below:

NSOL = 21 PV Based method using fast probability integrator

= 22 Monte Carlo simulation technique for material property related

uncertainties

NMCS = No. of samples for Monte-Carlo simulation (required for NSOL = 22)

*INDZON

The keyword INDZON is used to define the total number of statistically independent

zones in the structure. The general input format is

*INDZON

NDZON

where NDZON is the total number of independent zones in a composite structure.

*PLYGRP

The keyword PLYGRP is used to specify the total number of ply groups in each

independent zone. The general input format is

*PLYGRP

INDI(1) NLGRP(1)

INDI(2) NLGRP(2)

INDI(3) NLGRP(3)

INDI(n) NLGRP(n)

NASA/CR--2003-212127 5
Where INDI is the independent zone no. NLGRP is the total number of ply groups in an

independent zone.

There will be NDZON number of lines under this keyboard.

Note: *INDZON key word card must precede *PLYGRP key word card.

*NODGRP

The above keyword is used to specify the total number of node groups. A node group

used to identify the common independent zone and the common layer group for all the nodes in a

group. The general format is:

*NODGRP

NGRP

where NGRP is the total number of node groups.

*PLIES

The keyword PLIES is used to specify the maximum number of plies that exists in an

independent zone. The general input format is:

*PLIES

1NDI(1) NMAXPL(1)

INDI(2) NMAXPL(2)

INDI(3) NMAXPL(3)

INDI(n) NMAXPL(n)

where 1NDI is the independent zone No., and NMAXPL is the maximum number of plies

in an independent zone. The maximum number of plies for all the independent zones must be

defined.

Note: *INDZON key word card must precede *PLIES key word card.

*MATSYS

The keyword MATSYS is used to specify the total number of material systems in a

structure. The general input format is:

*MATSYS

NMAT

where NMAT is the total number of material systems in a structure.

*ORISYS

The keyword ORISYS is used to specify the total number of orientation systems in a

structure. The general input format is:

*0RISYS

NORI

NASA/CR--2003-212127 6
whereNORIisthemaximumnumberof orientationsystemsin astructure.

*THKSYS

ThekeywordTHKSYSisusedto specifythetotalnumberof thicknesssystemsin a

structure.Thegeneralinputformatis:

*THKSYS

NTHK

whereNTHK isthetotalnumberof thicknesssystemsin astructure.

*RNDGRP

Thiskeywordis usedto specifythemaximumnumberofrandomfield groupsrelatedto

materialpropertiesfor eachindependentzonein astructure.

Thegeneralinputformatis:

*RNDGRP

INDI(1)

INDI(2)

INDI(3)

INDI(n)

NaFT(1)

NRFT(2)

NRFT(3)

NRFT(n)

where INDI is the independent zone number NRFT is the total number of random field

groups (*CORTAB) related to the constituent material properties in an independent zone.

The number of random field groups related to the constituent material properties for all

the independent zones must be defined.

Note: *INDZON key word card must precede *RNDGRP card.

*RANTEM

This key word is used to specify the maximum number of random variables related

thermal loads and corresponding number of perturbation sets.

The general input format is:

*RANTEM

NRTEM NTPER

where NRTEM is the number of temperature related random variables and NRPER is the

total number of perturbations related to the temperature related random variables.

The number of random variables related to the temperatures for all the independent zones

must be defined.

*END

This key word card is always required to terminate the parameter data block input. It

must always be put at the very end of parameter data input block. The general input format is

*END

A.3.1.2 Material Property Uncertainty Model Input Data Group:

*NODGRP

NASA/CR--2003-212127 7
ThekeywordNODGRPisusedto specifythenodenumbersthatbelongto a specificply

groupof aspecificindependentzone.Thegeneralinputformatis shownbelow:

*NODGRP

IG

NODB(IG) NODE(IG) INDZ(IG) ILAY(IG)

where IG isthenodegroupnumber

NODB(IG) isthebeginningnodenumberof aseriesof nodesingroupIG,

NODE(IG) istheendnodenumberof aseriesof nodesin groupIG,

ILAY(IG) istheply groupnumberto whichNODB(IG)throughNODE(IG)

belongs

INDZ(IG) istheindependentzonenumberto whichNODB(IG)through

NODE(IG)belongs

Rememberthatthenodegroupsfor all theply groupsandall theindependentzonesmust

bespecified.

*RESTAB

ThekeywordRESTABisusedto specifythepliesandtheirresponsesfor whichthe

probabilisticsimulationisdesired.Theresponsetablefor alltheplieswithin anindependent

zonemustbespecified.Thecodeusedto specifythecomputationalrequestis0(zero)or l(one).

0 meanstheprobabilisticply responsecomputationisnotdesiredand1meanstheprobabilistic

ply responsecomputationisdesired.Thegeneralinputformatis:

*RESTAB

IZ

NPB,NPE

IEX(1)IEXE(2)... (EX(23)

where

IZ is theindependentzonenumber

NPB isthebeginningplynumberof independentzoneIZ

NPE istheendplynumberof independentzoneIZ

IXE(I) isthecodeto requestithprobabilisticply response

Followingis thelist of ply responsesfor whichtheprobabilisticanalysiscanbe

performed.

1Longitudinalstrain

2Transversestrain

3Shearstrain

4Longitudinalstress

5Transversestress

6Shearstress

7Longitudinaltensilestrength

8Longitudinalcompressivestrength

9 Transversetensilestrength

10Transversecompressivestrength

11Shearstrength

NASA/CR--2003-212127 8
12 Modified distortion energy failure criterion

13 Hoffman's failure criterion

14 Interply delamination failure criterion

15 Fiber crushing criterion (compressive strength)

16 Delamination criterion (compressive strength)

17 Fiber micro buckling criterion (compressive strength)

18 Failure in longitudinal direction

19 Failure in transverse direction

20 Failure in Shear strength

21 Out of plane Shear 13

22 Out of plane Shear 23

23 Sigma zz

*PLYTAB

The keyword PLYTAB is used to specify the ply configuration of a laminate in a ply

group of any independent zone. The ply configuration of a laminate is specified by the existence

of a ply. The code used to specify the existence is 0(zero) or l(one). 0 means the ply does not

exists and 1 means the ply exists. The general input format is:

*PLYTAB

IZ, LG,

NPL(1,IZ,LG), NPL(2,IZ,LG), ... NPL(NMAXPL(IZ),IZ,LG)

where

IZ is the independent zone number

LG is the ply group number of independent zone IZ

NMAXPL(IZ) is the maximum number of plies in an independent zone IZ

NPL is the ply existence code, either 0 or 1. Code zero indicates that

the ply does not exist and 1 indicates that the ply exists

*MATSYS

The keyword MATSYS is used to specify the uncertainties of a specific material. The

general input format is:

*MATSYS

I

'KEY(I)KEYS(I)' FVRM(I) VVRM(I) FRS(I) FVRMS(I) VVRMS(I)

C(1,I) C(2,I) C(3,I)... C(29,I)

CF(I), CVV(I), CFRS(I), CFS(I), CVVS(I)

CS(1,I)CS(2,I) CS(3,I)... CS(29,I)

IM(1,I) IM(2,I)IM(3,I)... IM(29,I)

IF(I), IV(I), IFR(I), IFS(I), IVS(I)

IMS(1,I) IMS(2,I) IMS(3,I) ... IMS(29,I)

NASA/CR--2003-212127 9
where I is the material system number KEY(I) and KEYS(I) are the acronyms defining

the I th primary and secondary composite material system. The acronym must be specified in

quotes ('). The acronym is used to pull the mean values of composite material properties from

the data bank. Refer to Table 3.4.1.1 and Table 3.4.1.2 of the IPACS manual for the description

of composite material acronyms. FVRM(I), VVRM(I), FRS(I), FVRMS(I), and VVRMS(I) are

mean values of the primary system fiber volume ratio, the primary system void volume ratio,

fraction of secondary system, the secondary system fiber volume ratio and the secondary system

void volume ratio respectively of the I th material system. C(n,I) is the coefficient of variation

for the n th primary material property random field of the I th material system. The list of the

random field numbers is given in section 3.5 oflPACS manual. CF(I), CVV(I), CFRS(I),

CFS(I), CVVS(I) are the coefficient of variation for the primary system fiber volume ratio, the

primary system void volume ratio, fraction of secondary system, the secondary system fiber

volume ratio and the secondary system void volume ratio respectively of the I th material system.

IM(n,I) is the distribution type for the n-th material property random field number of I th primary

material system. IF(I), 1V(I), IFR(I), IFS(I), 1VS(I) are the distribution types for the primary

system fiber volume ratio, the primary system void volume ratio, fraction of secondary system,

the secondary system fiber volume ratio and the secondary system void volume ratio respectively

of the I th material system.

IMS(n,I) is the distribution types for the n th secondary system material property number.

No more than ten integer numbers and eight real numbers in a given line can be specified. If

there are more numbers for a particular data type, then their input should continue in the

subsequent lines.

*ORISYS

The keyword ORISYS is used to specify uncertainties of the ply orientations. The

general input format is:

*ORISYS

OM(1) OM(2) ... OM(NORI)

OC(1) OC(2)... OC(NORI)

IO(1) IO(2)... IO(NORI)

where OM(i), OC(i) and IO(i) are the mean, a constant to specify standard deviation

(standard deviation -- 90 x constant), and the distribution type respectively of the ith orientation

angle system. NORI is the total number of orientation systems.

*THKSYS

The keyword THKSYS is used to specify uncertainties of the ply thickness. The general

input format is:

*THKSYS

TM(1) TM(2) ... TM(NTHK)

TC(1) TC(2)... TC(NTHK)

IT(l) IT(2)... IT(NTHK)

where TM(i), TC(i) and IT(i) are the mean, the coefficient of variation and distribution

type respectively of the ith thickness system. NORI is the total number of orientation systems.

NTHK is the total number of thickness systems.

NASA/CR--2003-212127 10
*MATTAB

ThekeywordMATTAB is usedto specifythematerialof eachplyof alaminatein an

independentzone.Thematerialpropertyuncertaintyspecificationof aply in anindependent

zonecorrespondsto amaterialsystemid. All thepliesof anindependentzonemusthavea

materialsystemid defined,andthematerialpropertyuncertaintyforthatid mustalsobedefined

in

*MATSYScard.Thegeneralinputformatis:

*MATTAB

IZ

NMT(IZ,1) NMT(IZ,2) NMT(IZ,3) ...NMT(IZ,NMAXPL(IZ))

whereIZ istheindependentzonenumber.NMT(IZ,n)isthematerialsystemnumberof

thenthply of independentzoneIZ, andNMAXPL(IZ)isthemaximumnumberof pliesin

independentzoneIZ.

*ORITAB

ThekeywordORITABis usedto specifytheuncertaintiesof eachplyorientationof a

laminatein anindependentzone.Theorientationuncertaintyspecificationof aply in an

independentzonecorrespondsto anorientationsystemid. All thepliesof anindependentzone

musthaveanorientationsystemid defined,andtheorientationuncertaintyforthatid mustalso

bedefinedin

*ORISYScard.Thegeneralinputformatis:

*ORITAB

IZ

NOT(IZ,1) NOT(IZ,2) NOT(IZ,3) ... NOT(IZ,NMAXPL(IZ))

whereIZ istheindependentzonenumber,NOT(IZ,n)istheorientationsystemnumberof

thenthply of independentzoneIZ, andNMAXPL(IZ)isthemaximumnumberof pliesin

independentzoneIZ.

*THKTAB

ThekeywordTHKTABis usedto specifytheuncertaintiesof eachply thicknessof a

laminatein anindependentzone.Thethicknessuncertaintyspecificationof aply in an

independentzonecorrespondsto athicknesssystemid. All thepliesof anindependentzone

musthaveathicknesssystemid defined,andthethicknessuncertaintyforthatid mustalsobe

definedin

*THKSYScard.Thegeneralinputformatis

*THKTAB

IZ

NTT(IZ,1) NTT(IZ,2) NTT(IZ,3) ... NTT(IZ,NMAXPL(IZ))

whereIZ istheindependentzonenumberandNTT(IZ,n)isthethicknesssystemnumber

of nthplyof independentzoneIZ, andNMAXPL(IZ)is themaximumnumberof pliesin

independentzoneIZ.

NASA/CR--2003-212127 11
*CORTAB

Thiskeywordis usedonlyin caseof aprimitivevariablebasedmethodto specifythe

correlationof aparticularrandomfieldbetweendifferentplies(correlationtable)of acomposite

material.RefertoAppendixA forthedefinitionof acorrelationtable.Thegeneralinputformat

fortheprimitivevariablebasedmethodis:

*CORTAB

IZ

NaFi NaFj

ICOR(1) ICOR(2) ... ICOR (NMAXPL(IZ))

Where IZ is the independent zone number

NRFi is the beginning material related random field number. NRFj is the

end material related random field number (Legal list of material related

random field numbers are given in Appendix A, Section A.3).

ICOR(i) is the correlation id number for the ith ply 0<= ICOR(i) <=NMAXPL(IZ)

NMAXPL(IZ) is the maximum number of plies in independent zone, IZ.

The cards (NRFi NRFj) and (ICOR(1), ... ,ICOR(NMAXPL(IZ)) can be repeated for

additional random fields in independent zone, IZ. For additional independent zones the data

should be specified in the same order discussed above.

Remember that the correlation table for all the material related random fields of all the

independent zones should be specified. The random field automatically becomes deterministic if

its correlation table is not defined.

*MEANTM

This key word is used to specify mean ply temperatures at different nodes. The

temperature groups and the number of node groups must be exactly same. The correlated

temperatures at different nodes can be specified using the *RELPER card. The standard

deviations for the temperatures are specified in the *STDVTM card. The first ply is the bottom

ply based on the local z-definition. The general input format for the definition of the mean ply

temperatures is:

*MEANTM

NODi NODj

PT1 PT2 ... PTn

NODi is the beginning ply number.

NRFj is the last node number in a group.

PTi is the mean ith ply temperature

The cards (NODi NODj) and (PT1, ... ,PTn) can be repeated for additional nodes.

*STDVTM

This key word is used to specify standard deviation of the ply temperatures at different

nodes. The temperature groups and the number of node groups must be exactly same. The first

ply is the bottom ply based on the local z-definition. The general input format for the definition

of the ply temperature standard deviation is:

NASA/CR--2003-212127 12
*STDVTM

NODi NODj

PSTD1 PSTD2 ... PSTDn

NODi is the beginning ply number.

NRFj is the last node number in a group.

PSTDi is the ith ply temperature standard deviation

The cards (NODi NODj) and (PSTD1, ... ,PSTDn) can be repeated for additional nodes.

*RELPER

This key word is used to specify the autocorrelation matrix for random temperature

variables. Basically, it represents the eigen-vectors of the correlation matrix. The general input

format for the definition of the ply temperature standard deviation is:

*RELPER

NR

NODi NODj

EG1 EG2 EGn

NODi is the beginning ply number.

NRFj is the end 1st ply number.

EGi is the ith ply autocorrelation coefficient

The cards (NODi NODj) and (EG1, ... ,EGn) can be repeated for additional nodes.

*IRVPER

This key word is used to specify the perturbation of particular temperature random

variables. Any random variable can be perturbed for any number of times. It also signifies in

what order the variables shall be perturbed. Input entry in this card must be used in conjunction

with the *ABSPER card which specifies the corresponding perturbation magnitude. Basically, it

represents the eigen-vectors of the correlation matrix. The general input format for the definition

of the temperature random variables perturbation definition is:

*IRVPER

NP(1) NP(2)... NP(n)

where NP(i) represents the i th perturbation of random variable NP(i).

The cards (NPi NPn) can be repeated for additional variable perturbations.

NASA/CR--2003-212127 13
*ABSPER

Thiskeywordis usedin conjunctionwiththe*IRVPERinputcard.It isusedto specify

themagnitudeof theperturbationforthecorrespondingrandomvariableinthe*IRVPERcard.

Thegeneralinputformatforthedefinitionof theperturbationmagnitudeis:

*ABSPER

PV(1) PV(2)... PV(n)

wherePV(i) representsthemagnitudeof theithperturbationofrandomvariableNP(i)

specifiedin *IRVPERcard.

Thecards(PVi PVn)canberepeatedfor additionalvariableperturbations.

*REFTEM

Thiskeywordis usedto specifythereferencetemperatureof thematerial.Thegeneral

inputformatforthedefinitionof thereferencetemperatureis:

*REFTEM

IZ

REFTM

whereIZ representstheindependentzoneandtheREFTEMthereferencetemperaturein

°F units. ThecardsIZ andREFTMcanberepeatedfor additionalindependentzones.

A.3.2 Input for Probabilistic Output Request:

This section of input always begins with $FPI card. The response locations on the

structure and the type of probabilistic analysis are defined in this section of input. The input

details are given in the following card:

*STRU

ISTAT NTYPE NODN IDUMP IDUMP NCOMP

where

ISTAT

NTYPE

NODN

IDUMP

NCOMP

= 1 for the Static Analysis

= 2 for the Buckling Analysis

= 3 for the Frequency Analysis

= 1 for Displacement Response

= 2 for Ply Strain/Stress Response

Node Number where probabilistic output is requested

Reserved for future use

displacement or stress/strain component in static analysis

mode number in case of buckling or frequency analysis

NASA/CR--2003-212127 14
A.3.3 IPACS Execution Control Input:

The input to this section always begins with $EXE2 card.

is defined in this section. The input to this section is as follows:

IX1 IX2 IX3 IX4 IX5

The IPACS execution control

where the variables IX1, IX2, IX3, IX4 and IX5 could be either 0 (zero) or 1. 0 means

suppress the execution of the analysis type discussed below. 1 means perform the analysis.

IX1 = 1 Compute the perturbed material properties at every scale.

IX2 = 1 Prepare the input file for the probabilistic finite element analysis.

IX3 = 1 Perform the probabilistic finite element analysis.

IX4 = 1 Extract desired structural responses as indicated in $FPI section.

IX5 = 1 Compute the CDF of desired structural responses

A.3.4 Output

The user generally controls IPACS output. However, by default IPACS provides minimum

necessary output required for the probabilistic assessment of composite structure. The minimum

output consists of the following: (i) an echo of the input, (ii) CDF of the material properties and

the desired structural response, (iii) the sensitivity factors of the primitive variables for the

material properties and the desired structural response. The important output files are

xxx21 .fpimov, xxx21 .nesout and xxx21 .rantab. The names of these files and their respective

contents are listed in Table A.3.4.1, on the next page.

Several other output files are created at intermediate computation stage or at the end of the

computations. Some of these files provide an interface between different modules or serve as a

database system. These files contain information for the user to study the overall problem

behavior in greater detail.

NASA/CR--2003-212127 15
Table A.3.4.1 List of Output Files

jn- jobname

IZ Independent Zone No.

LG - Ply Group No.

IPM - Method No.

$$ A unique job number assigned by the computer

File Name File Contents

jn.msp-IZ-LG Mean value and standard deviations for ply material

properties.

jn.cpp-IZ-LG CDF of ply material properties

jn.cpl-IZ-LG CDF of laminate material properties

jn-IZ-LG-IPM- If IC = 0 Mean Values of material properties. Perturbations

IC.PIC of A, C, and D matrices terms for primitive variable based

method

jn-IZ-IPM.dat Data required for finite element analysis. The file contains

the data supplied by the user and that generated by the

program

jn-IZ-IPM.pdb Perturbation database.

jn-IZ-IPM.datu Perturbation data for each desired response.

jn-IZ-IPM.dist Statistics of Primitive Variables for each independent zone.

jn-IPM.datu FPI input file.

jn-IPM, dist Statistics of primitive variables required for the probabilistic

analysis of structural response.

jn-IPM.fpinp Input file prepared by program to perform fast probability

integration.

j n-IPM.fpimov Sensitivity factors of primitive variables at cumulative

probability levels of 0.001 and 0.999.

j n-IPM, fpibin FH logfile

jn-IPM.nesout Discretized CDF of a structural response.

jn-IZ-IPM.out Echo of FEM input and Results of Unperturbed Solution.

jn-IBM.FPIMOV Sensitivity factors for independent random variables at

probability levels 0.001 and 0.999.

NASA/CR--2003-212127 16
SAMPLE INPUT FILE:

CFEM

SPIC1

* INDZON

1

*NODGRP

1

*PROSOL

21

*PLYGRP

1 1

*PLIES

1 8

*RNDGRP

1 2

*MATSYS

1

*ORISYS

4

*THKSYS

1

*RANTEM

1 4

* END

*NODGRP

1

1 9

*RESTAB

1

1 8

c o-skip fpi

0 0

0 0

0 0

* PLY TAB

1 1

1 1

*MATSYS

1

'AS--EPOKAS--EPOK '

0 05 0.05

0 05 0.05

0 05 0.05

0 05 0.05

0 05 0.05

0 05 0.05

0 05 0.05

2 2

2 2

2 2

2 2

2 2

2 2

2 2

* MAT TAB

1

1 1

*ORISYS

-45.0 0.0

0.02 0.02

2 2

*ORITAB

1

3 1

1 1

1-with fpi

0 1

0 0

0

1 0 0 0 0 0

0 0 0 0 0 0

1 1 1 1 1 1

0.60 0.02

0 05 0 05

0 05 0 05

0 05 0 05

0 05 0 05

0 05 0 05

0 05 0 05

0 05 0 05

2 2

2 2

2 2

2 2

2 2

2 2

2 2

0.000001

0.05

0.05

0.05

0.05

0.05

0.05

0.05

2

2

1

2

2

2

1

0.60 0.02

0.05 0.05 0.05 0.05

0.05 0.05 0.05 0.05

0.05 0.05 0.05 0.05

0.05

0.05

0.05 0.05 0.05 0.05 0.05

0.05 0.05 0.05 0.05 0.05

0.05 0.05 0.05 0.05

2 2 2 2 2

1 1 2 2 2

1 1 2 2

2 2 2 2

1 1 2 2

1 1 2 2

1 1 1 1 1 1

45.0 90.0

0.02 0.02

2 2

2 4 4 2 1 3

NASA/CR--2003-212127 17
*THKSYS

0.020

0.050

2

*THKTAB

1

1

*CORTAB

1

1

1

33

1

*MEANTM

1

200.0

*STDVTM

1

10.0

*RELPER

1

1

1.0

*IRVPER

1

*ABSPER

-2.0

*END

$NES

*FEM

*DISP

*FORC 1

*ELEM 4

75

*COMPOSITE

*TEMP

*PERT

*NODES 9

*BOUN 18

*CONS 0

*PRIN

*END

*INC

0

*ITER

100

*COOR

1

2

3

4

5

6

7

8

9

*ELEM

1

2

3

4

*BOUN

1 1

4 1

7 1

1 2

4 2

1 1 1 1 1 1 1

1

1

33

1

9

200.0

9

10.0

9

1.0

-1.0

1

2

0

1

2

0

1

2

75

1

2

4

5

200.0

10.0

1.0

1

1.0

0.0001

0 000

000

000

000

000

000

000

000

000

0.0001

200.0

10.0

1.0

1

2.0

0 000

0 000

0 000

1 000

1 000

1 000

2 000

2 000

2 000

2 5 4

3 6 5

5 8 7

6 9 8

200.0

10.0

1.0

0.0001

0.000

0 000

0 000

0 000

0 000

0 000

0 000

0 000

0 000

200.0

10.0

1.0

0.0001

1 0

1 0

1 0

1 0

1 0

1 0

1 0

1 0

1 0

200.0

10.0

1.0

200.0

10.0

1.0

NASA/CR--2003-212127 18
7 2 0

1 3 0

4 3 0

7 3 0

1 4 0

4 4 0

7 4 0

1 5 0

4 5 0

7 5 0

1 6 0

4 6 0

7 6 0

*FORC

5 3 1000.0

*ORIEN

1 9 1.00

1 9 0.00

C*LAMI

INCLUDE

*PRINTOPTION

REACTION

STRESS

STRAIN

TOTALDISPLACEMENT

*END

SRAN *DEFI 1

1000.0 100.0

FORCE

5 3 1.0

$FPI *STRU

1 1 5

SEKE2

1 1 1

0.00

0.00

0.00

1.00

0 0 3

1 1

NASA/CR--2003-212127 19
Form Approved REPORT DOCUMENTATION PAGE

OMB No. 0704-0188

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources,

gathering and maintaining the data needed, and completing and review#lg the collection of information. Send corrlments regarding this burden estimate or any other aspect of this

collection of information, including suggestions for reducing this burden, to Washington Headquarters Services. Dhectorate for Information Operations and Reports, 1215 Jefferson

Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget Paperwork Reduction Project (0704-.0188), Washington, DC 20503.

1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3, REPORT TYPE AND DATES COVERED

February2003 Final Contractor Report

4, TITLE AND SUBTITLE 5. FUNDING NUMBERS

User's Guide for the Updated EST/BEST SoLlware System

6. AUTHOR(S)

Ashwin Shah

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES)

Sest, Inc.

18000 Jefferson Park Road

Middleburg tteights, Ohio 44130

WBS-22-708-90-53

C-7675 l-J

8. PERFORMING ORGANIZATION

REPORT NUMBER

E----13779

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSORING/MONITORING

AGENCY REPORT NUMBER

National Aeronautics and Space Administration

Washington, DC 20546-0(101 NASA CR-------2003-212127

11. SUPPLEMENTARY NOTES

Project Manager, Christos C. Chal-nis, Research and 'I)chnology Directorate, NASA Glenn Research Center, organization

code 5000, 216---433.---3252.

12a, DISTRiBUTION/AVAILABILITY STATEMENT

Unclassified- Unlimited

Subject Categories: 05, 07, and 35) Distribution: Nonstandard

Available electronically at http://a]trs.grc.n_sa.gov

"]['his publication is available from the NASA Center for AeroSpace lnfonnatiom 301-621-0390.

12b. DISTRNBUTION CODE

13. ABSTRACT (Maximum 200 words)

This User's Guide describes the structure of the IPACS input file that reflects the modularity of each module. The

structured format helps the user locate specific input data and manually enter or edit it. The IPACS input file can have

a W user-specified filename, but must have a DAT extension. The input file may consist of up to six input data blocks;

the data blocks must be separated by delimiters beginning with the $ character. If multiple sections are desired, they

must be arranged in the order listed.

14. SUBJECT TERMS

Input; Output; Sample cases

17. SECURITY CLASSIFICATION

OF REPORT

Unclassified

NSN 7540-01-280-5500

15. NUMBER OF PAGES

25

16. PRICE CODE

18, SECURITY CLASSiFiCATiON 19. SECURITY CLASSiFiCATION 20. LiMiTATiON OF ABSTRACT

OF THIS PAGE OF ABSTRACT

Unclassified Uncl assified

Standard Form 298 (Rev. 2-89)

Prescribed by ANSI Std. Z39-18

298-102