Transcription
NASA-CR-198795//tCeramicANSYSMaterialLife PredictionResults- A Programto CARES/LIFEto TranslateReliabilityAnalysisFinal ReportDr. PieterPeriodVon Hermann/AdamCovered:Cleveland3/1/88Pintzto 4/26/94State University1983 E. 24th -CR-]9B?95)CERAMICMAIERIALLIFEPRFOICTION:A ITYANALYSISFinal1994Report,(ClevelandiMar.I988- 26StateUniv.)N95-29367UnclasApr.29 pG31270055114
ANSCARESA ProgramTranslatingANSYSfor CARES/LIFELynnFiniteReliabilityM. Powers"Clevelandand JohnStateCleveland,ElementResultsAnalysisW. Spring s)used to readfile using Resident ResearchVersion5. 0resultssubroutinesto CARES)Associate5.0, 5 in the formatfor use by experiencedFORTRANLewisprogramsResearchof existingto prepareneededusersa neutralandThefileofby the CARES/LIFE'of ANSYSis also required.Centerroutinesthe f'mite elementFORTRANManufacturingprogramand o extractFORTRANANSCARESReleaseIt is intendedis madethe ANSYS(MSC/NASTRANt AdvancedStateanalysis.ANSYSby NASA'suseof theand linkingdevelopedat ClevelandANSYSthe usefromof nmanual- NASAandfile.resultsmachineare nearlyother CARESand esthe neutralas describedare computedidenticalinterfaceto thoseandusedprogramsandfile for inputto thein Reference6 arewrittento a neutralin the NASCARESinterface.Lewis Research CenterCenter- 1 -May,1995
CeramicsReliabilityIncreasingadvantagenumbersof highermaterialsAnalysis:temperaturedo not followCARES/LIFEProgram)of componentsprogramand reliabilitysimpleAnalysishas beenof ceramicProbabilisticand corrosionthe same(Ceramicare beingresistancestrengthandmadeby NASALewisUnfortunately,that characterizeEvaluationofResearchin order to evelopedof ceramicductileStructuresCenterceramicmetals.Lifeto predictThePredictionthe dedare performedmany simpleor cyclicload is lityof failureconcurrentTheseissuesTheof a finitecommercialversionFEMdistributionelement5. 0versionswithin(FEM)ABAQUS,3, andpreparetheANSYS,a neutralis ue)and the ndanacceptablethe statisticalto allownaturefor multiaxialsubcriticalby a 4 (ANSYS/CARES)file, in standardized- 2 -and lations,the CARES/LIFEmethodprograms67 ic,for a complexareexperimentsA static,strengthparameterschangesof crackthesetest tyThis designoccurring)withinor tensileUsingthe mechanics(simultaneouslyare addresseddata.Appropriatewithrupture data.until fracture.thesehas been reached.flawsflexuralthe time-dependentpredicted.of strength-controllingstates,geometryfromanalysis,is thenfrom specimento each specimenare then determineda finiteresultscomponentfromdevelopedtheto the2, MSC/NASTRAN4.The ABAQUSis acceptedandby the mostMay,1995
recentreformulatedCARES/LIFE program,while the ANSYS version usesthe older conceptofreading the FEM results through subroutinescalled from inside the CARES program. Thiscurrenteffortdevelopsand CARES/LIFE.NeutralDataTheto countingat eachGaussfilesversionsto CARES/LIFEof ANSYSis alsofor the interfacingthe r theadvantagefile)of thecomputergroupsto be consideredof resystems.aref'mitesoftwareused.into theof transferreliabilityinputviaanin eachof thisInterpreterin theusuallybyassembledis the easeto differentin theinterpretedelementevaluationof this approach(availableto beto the finitereliabilityof CARES/LIFEanalysishavefor outputhas to be adaptedrequiredfor sA bution5. centcontrolanalysis,an interpreterfile.A subelementfinitebase.as theTheresultsdatabasesuchwith the meters,to theelementneutralneutralor binaryDuefiniteformattedfor ANSYSpackages.outputprogram,is compatibleConceptprogram.Thethatof a formattedof standardelementflagsA translatorBasefiniteinterfacean f usingare determinedevaluated.correspondsthehas beenTheto a subelement.pointto theelementreliabilityaveragingat the GaussiansubelementTheelement-3 -analysisis now iessubelementvolumevolumein theIn the contextthator areacourseofwherethe localeachGaussianis definedof theas thenumericalMay,1995
integrationprocedure.Theas welllocationof the Gaussianas correspondingis seforsuppliesin the ement,subelementelementin stressreliabilityof the f'mitethein eachelementof thethespaceto be consideredof subelementsin the accuracydatapointas a on,chosen,integrationon fortheireliabilityanalysis,structureof the neutralis arrangedwithinand s),two-dimensionalconditionsor exploitingusingsubstantiallyTheon theVersion5. spectdata containfor volumeelementsareof the structure,mesh.Thus,thetheelementdataByof degreesof CARESfor surfaceassumingelementelementsof elementsversionelementssize of the finitenumberThethe numbermembraneimplemented.points.groups,in the leflawThehierarchy:informationtypesat the integrationto memory.the followingto the dentificationis optimizeda two-dimensionalthe areaare not takentensorsflawplanestressmodelof freedomcan becouldbereduced.methodmultiplyingstressdataIn addition(three-dimensionalreducedbasethe neutralThethe group.themodelingto determineof a evolume.Forby the thicknessthe stressto the elementsFinally,availabletheANSYSsubelement- 4 -loading5.0,theof the element.is assumedinclude(in planegroupvolumeTherefore,to be constantthe numberdatavs. bending)throughof subelementscontaindependis calculatedbendingbyeffectsthe thickness.and the materialinformationregardingMay,1995
subvolume, stress tensor (i.e., resataretheassumedintegrationpoint),to be reliabilitygroups(bothanalysisshellon the globalevaluation.volumetheTheas welland the a flagoutlinedf'u'stnextspecifiesthethetitlenumberto be consideredas well as surfaceinto recordsfortheof elementin theflaw ementtypegroups,of thetypeindicatorelementcodes,havea loopthe numberto be input.pertainingof elementsthe elementgroupThe elementnumberanalysis,elementlevelof subelementsthematerialthicknessin the finiteThe numberVersionof the5. 0numberelementelementof shellof theelementof subelementsin the f'mitecenterin elementon thecaseis readelementflawbe thein NASTRAN- 5 -type.group,in a loopandtoa shellthe indicatorover the numberof the elementanalysis,and the averageddependscorrespondingin the analysis,identificationif further postprocessingthis elementtypethis elementnumberfor surfaceshould(format415,3F15.0).the materialelements,statementelementare usedlevelthe identificationnumbernextwithinto the elementthis inthetheof elementsconsidereddata containwithconsidered,If shellInformationwithinis startedgrouphas to be set to 1.TheA80)215) containselements)volumefile data is arranged(format(format(shell)combinedthe neutralrecordentryas surfaceindicatingabove,samenumbertheelementon the Gaussianintegrationtheflawtheof this element.resultsdata is onlypreparedvolume,correspondingof the reliabilityanalysesfor volumetemperatureas theIf stressconsidered,elementare planned.order chosenavailablefor previousat theCARESMay,1995
versions),one subelementthe materialcontrolnumbersonlyevaluationand the fullcoordinateis startedfor readinghaveand cr y are necessaryvolume,thethe volumeinterpretersubelementvolumesusingof finiteelementpyramidor eme,5. 0theoperatorintegralFor thearereliabilitytype elementsthe subelementtemperatureo ,err, and oin theInsteadonlyo ,, Cryyof the subelement(I5,5F15.0).with standard finite elementto a Gaussinherentintegrationto the elementpoint)type.element,into the naturaloutput.Thus,is calculatedIn the usual(i.e., brick,coordinatein thecontextwedge,space,Idet-ItemperatureFor that reasontensor.of a three-dimensionalv fffd is the Jacobianvolume,of the subelementIon theFor volumeby cr,, , o ,, o ,after transformation-1data.are assumed.functionsis F15.0).are not includedIwhere(format(correspondingthe shapefile.correspondnumbersand elementin the printout(definedarea has to be specifiedof each subelementprogrampointstressfor the definitionsubelementTypically,planesummariesmustas the materialvolumethe subelementto be specifiedFor shell type elementsElementthe subelementat the integrationsystem)as wellnumbersused.of the .data are directlynumberstresselementgroup informationthe subelementthe identificationlocalthe Weibullfor elementA third loopThe materialused in the finitefile (containingnecessaryis specified.J (r,s,t)(1)dr ds dt-Iand r,s,t are the naturalcan be expressed- 6 -coordinates.By applyinga GaussianasMay,1995
nV wheren is thenaturalspaceor triangle)integrationandnEEi lj lk lorder,the W'snEin a similarfunctions.Thevolumeof an axisymmetricyi lj lANSYS/CARESAnalysisFractureprogramusingof volumeprobabilitya f'miteelementssurfacedatashouldhaveof thecomponent.descriptionskeywordin the5. ,J(r/,sj)of an elementpointsin(i. e., quadrilateralW Wj(3)be determinedsimilarlydet J(ri,sj)W i Wj(4)j lpointon a ceramicmodelof the partelementsin the sin question.attachedSupportedIt is vitalagreecoordinates,elementmembranedetof the em.Procedure:to i lR is the radialThecan alsonwhereof the(2)ItelementV W i Wj W tmannerI1A J(ri,sj tk)the coordinatesri,sj,t k areare the weightis determineddetat theelementtypesand thicknessesmaterialThissurfacesarenearin themodelzeroin Tableto avoidusedCARES/LIFEof the ANSYSwill typicallyto providelistedidentificationscontainedwith the executionboth1.in thevolumeShellaffectingconsistelementsthe romand stressthe ANSYSlinearcomponentsat all nodes.- 7 -elasticanalysisresultsThe ANSCARESfile includenodalprogramis thenMay,1995
usedto computesubelementstresses,volumes,andprogram.All neededis typeconsistsanalysisof a mainrenamedcalledandthef'mitetypedefaults,doesnot handleusingfracturethe neutralproperties.The1 reandand no userif desired.computed.elementby the CARES/LIFEfile dThismustfileANSCARESThesecodecommand5. 0andis providednineTheseloadcases.file preparedtemplatea flowbyfile alsodiagramsubroutinesneededall of whichbe copiedbe copiedfromtheareof aANSYSmachineBINCLO,languageto read dataBININI,in the ectorycodeanscares2.f:includedwiththe ANSYSBINRD,fileinand compiled.subroutinesBINIQR,to theANSCARES77 source(anscaresl.finto a suitableare designedare includedproblem.in FORTRANANSCARES2.FORshouldroutinesan exampleall supplieduse of system-dependent,tape.SYSIQR,andThis modulemakeson disk nanalysisthe CARES/LIFEprogramfor unix systems).the ANSYStheresultsANSCARESfile containingusingprogramANSCARESobjectandthereforeis to executea templatethe ANSYSnamingper run,ANSCARESmodulesfile.fromas neededANSCARES:Thetwofileof analysisreliabilityInstallingcasestepandtheare readto changeloadThe to runtheprogram.machineto createlanguagean executableto do thisANSYSversionon a VAX/VMSroutinesmustbe linkedof ANSCARES.computer- 8 -withIf theythe compiledFORTRANare all in the samelibrary,isMay,1995
LINK ANSCARES,BINLIB/LIBRARYfora silicongraphicsf77 ossibly-g -o anscaresIRIS.a onsare cares2.o/ansys50/objects inlib.acommandon otherthewhichwereusedspecific,systems.to linkprogramand can be expectedHowever,theseunderStand-Alonefor your systemtheare describedon a SGIto makeuseofin the e .theToexampleconfirmbe executedneutralfileproblem,thein turnproperusingthese(ANSCARES.NEU)andfor comparison.considerations:At sprogram.analyses.Tableprogram.To provideelementdescriptionVersion5. 01 describesemployingTheseANSYSelementthe ANSYSthe dedelementallowtypestypesby ANSCARES,- 9 -beenthree-dimensional5.0 elementset to a valuehasof twosupportedall element(2).Elementimplementedandintoaxisymmetricby the ANSCAREStypesmusttypeshavethethat usetheMay,1995
rogrampromptENTERkey willacceptmessageson nvalueAll needed5. 0thanothers.and elementof nodesITOTin thetheof "FILE.RST"of the ANSYSand"CARES.NEU".Again(ICARES 2)willallowwhento bepressingandthe ssingand warning/progressis highlyof nodeselementselementsizesize of the A arrayparameter,namesanddata is storedandThe declaredbe re-compiledparameter(0).thedisplayto specifyofthatis set to 9 insteadto changethe defaultofvaluefrom 6.nodeand is relatedfile to be generated.it may be necessaryon the numbersomeneutralfile instead,system,differentOf course,routine)mustof zerofor the namesof "n" to the promptlimitnumbermaximumMTOTvaluepromptof subelement"ANSCARES.OUT"to somein this array.spacedefaultanalysisanalysisvalueson the computermaximummoreSubelementA responseis no specificelementthe defaultor subelementdefaultin thedisplayThereavailabletheDependingfor monitorat theand the CARESwill acceptor to a file.savedremainwill interactivelyare for elementsent to the monitor(default)its valuefile to be readkey for eachAdditionalhaveexecutionANSCARESand geometryENTERshouldtypesof theelementspermittedin a generalpurposeis limitedrequireonlymoreamountANSCARESmainand all programsANSCARESof computerprogram,re-linkedprogram- 10 -MTOTamount(the Aof spacebutprogram)andor decreasedasin the ANSCARESresourcesto establish(mainarrayfinitedata for each elementarray can be increasedis set by the parameterto theby thein thenotavailable.Afternecessarilythe newthelimits.is used to defineMay,the1995
size of thevariabletemporaryRECLNGneededcontainsby the VAXprinta size of bufferif thesystemparameterThethe minimumbufferwill resultthenthe ANSCARESmainprogramwillspace.Tomoreneeded.in a messagebe giventhatprogramresultsset to the sizeThe ANSCARESrequirescanANSYSis currentlyin useto 1 whichITOTfile.ITOTwords).to determinevariableresultsneeded.(7680computeris executed.the ITOTthe ANSYSof ANSYSthe ITOTANSCARESchangingreadingit is necessaryby settingwhenwhenthe minimummessagebe cangivingthevalue.mustAswithbe re-compiledtand all programsANSCARESthe ANSYSroutinesresultsthis filefile.ANSCARESmainresultsfile.the major2 showsANSYSresultsroutinesare basedof Ref.calledVersionbuffersize.Thisfile is describedare readon andin Chapterprogramfile,with the ANSYSservesby theremain5, alsoin Refto establishof the varioussubroutinesveryEANSYS,similarare calledin Chaptersizesand call the subroutinesfunctionsprogram6 of RefRDGEOM,to the example[6].Theto readroutines[6].for arrays,that preparemodulesas neededcall the subroutinesthe CARES/LIFEneutralin the ANSCARESand RDSOLU.programsbeginningprogram.Thesethreeon pageA-21[6].Thenumbersuppliedare documentedthat read the ANSYSFigurethe newDescription:languageanalysisto accessTheto establishProgramMachineusedre-linkedEANSYSof nodesby EANSYS5. 0subroutineand elementsto extractopensthe ANSYSin the model,coordinateand otherlocations,- 11 -resultsfile ity,determinesRDGEOMandproperties.May,theis thenThis1995
datais storedin the A array for later use. A countis alsomadeof the variouskinds of elementsin the model. EANSYS thencalls subroutineRDSOLU to extractstressandtemperatureresults,which arealso storedin the A array.This completesreadingof the datafrom the ANSYSresultsfile.The HEXA,PENTA,TETRA,by ANSCARESas needed,are veryto the same-namedsimilarQUAD8,dependingTRIA6,on the kindssubroutinesand TRIAXof elementssubroutinesare then calledin the model.in the previousTheseNASCARESsubroutinesprogram. 3 ExampleProblem:Sampleto verifyproblemsModelsthe accuracyincludingaxisymmetricReferenceof theelement4, is describeda sectorsuitablymakingeach sectorelementselementandsix sidesandTheVersioncircularmid-sideprogramloadedto preserveelements,operation.for earliernodesandversionsusingOne of the models,useduplicatedThe remainingfour edges.platehas an outsidesector,plate.In ANSYSis dividedof ANSYS.solid,shell,EXAMPLE1andfromnodes.Thenode numbersof 2.00the element- 12-inchesthis is treated asinto 48 equal sectors,into 40 three-dimensionalfour solidto make36 solid and 9 shell elementsdiametershowingcircularThe plate is subdividedof the sectorall with mid-sidewithof theceramicsymmetry.The modelof the modelshell elements.shellswithoutto interfacessequencebelow.7.5 wide.3 is a side view5.0confu'medrestrained10 shellas comparedandis a transverselyat the centerand triangularwithtypesusing the igurewere executedelementsfivenumbers.Thesidedremainand a thicknessandsolidthe shellprismsolidas solidswithof 0.071four layersinches.shownMay,are1995
madeup of the solidlayer.Figures4 and 5 showis simplysupportedinches).A uniforma uniformTheelements,70 Fbetweenwith the ten shelladditionalelementspressuredetails9 andin all elements.MaterialNumber 300 (Must- Young'sModulus 5.87- Density- templatefor thisCircularDisk.Version5. 0agreefileexamplein this problemwith CARESat a radius of 0.9239.Temperatureare givenisbelow:template)inchin Appendixto demonstratefile is supplied31 throughThe sectorx 107 psiis includedis providedusedof the lowestnumbering.to a ring supportto elementspropertiesto the bottomand element10 (equivalent 0.5 pounds/cubic 0.25inputattachedof the nodeof 220 psi is applied- MaterialentireelementsA.modelgenerationto run CARES/LIFEis describedin reference- 13 -,using[7]:andneutralthe generatedExamplefile design.TheneutralThe1 - Transverselyfile.LoadedMay,1995
ANSYSElementNumberTypeof Nodes4Axisymmetric- 4 EdgesNot3Axisymmetric- 3 EdgesKEYOPT(3) 186Brick - 6 SidesPrism5 Sides4Tetrahedron4Shell- 4 Edges3Shell- 3 EdgesKEYOPT(1) 072-SOLID4Tetrahedron82-PLANE8Axisymmetric- 4 Edges6Axisymmetric- 3 hell- 4 Edges3Shell- 3 Edges20Brick- 6 Sides15Prism- 5 Sides1310Pyramid- 5 SidesTetrahedron- 4 SideselementtypesKEYOPT(5)KEYOPT(6)TableVersion- 4 Sides10For all of the ANSYS5. 01 - ANSYSelementsupported- 4 Sides92-SOLID95-SOLIDNote:CommentDescriptionNot SupportedKEYOPT(3) 1- 4 SidesNotSupportedabove: 2-- 0types- 14 -supportedby ANSCARESMay,1995
eractiveMODEL[ PATRANGENERATION] [preandpostprocessingiHEATTRANSFERANALYSISat ementFNA output :databaseI ES CERAMICSANALYSISEVALUATIONOFI[C3LIFE]{CARES.I COMPONENTVOLUME,Figure1 -Version5. ELEMENTRELIABILITYSURFACEANDINTENSITIESRISKOFBlock diagram for the analysis and time-dependentreliabilityevaluation ofceramic components.This analysis includes the CARES/LIFEintegrated designprogram, finite element analysis and, PATRAN PLUS pre- and post-processing.- 15 -May,1995
Begin ANSCARES programCheckfor ANSYS results fileCall EANSYSSubroutineExtractdisk pointersCall RDGEOMandmodelparametersSubroutineExtractandFor eachsavenodecoordinatesand saveelementelement:ExtractExtractelementSaveAddFor eachconnectivityin temporaryto elementelementtypeRDSOLUeachandSaveFor eachelementstressestemperaturestemperaturesat nodesin nodeorderfile nodesto save arrayelement:ExtractWritecounterSubroutineFor eachForspacetype:EliminateduplicateMove connectivityCallpropertiestypeappropriateForin the model:elementtype 8,TRIA6,or TRIAX)of this omputeaveragetemperatureWrite results to neutral fileCloseEndfilesANSCARESprogramFigureVersion5. 02 - ANSCARES16 -ProgramLogicMay,1995
UniformPressureon Elementsof 220 psi31 through39--Element1.071"Number cp3233435363722232425" 262712131415161734563839402829i30.18192089i0EuLLJ S."OU').f.e-.- . .J0 S.0Q.C).021 .aJuE-t,,o. istance3 - SideviewfromCenterof a sector- 17 -ofDiscof the ceramicSupport{in.)plateApril,1994
FigurePreliminaryVersion4 - Solidelements-18-in the exampleApril,1994
FigurePreliminaryVersion5 - Elementplot of bottom- 19-shelllayerApril,1994
CaresPaperNo.and l",Ohio,J.P., "CeramicAnalysis(CARES/LIFE)andUsersandto be published.44601,A., and nger,L.A.,TM-106316,'ABACARESthe RestructuredTurbineNASADivision,L.M.,Janosik,of y&Program",Researchand7, 1993.J.P., "CeramicSponsoredWilcox,ComponentPresentedby ASME,ReliabilityWithat the 37th Nemeth,Reliability5. s a Pre-ProcessorSwansonAnalysisSystemsManualVol I, II, and ,Rcv.ANSYS"User'sManualJuly,5.0 EngineeringProgrammer's- 'sBeta1992.N.N.,Manderscheid,J.M.,and Gyekenyesi,Evaluationof Structures (CARES)",NASA- 20 -J.P.,"CeramicAnalysisTP-2916,Aug. 1990.andMay,1995
7/TITLE,STRESSET, 1,95KEYOPT, 1,5,2ET,2,93KEYOPT,2,5,2ANALYSISOF A 87E07DENS,300,0.5NUXY,300,0.25TREF,70csys, ,0.22,0.,0.N,9,0.30,0.,0.N, 10,0.42,0.,0.N,11,0.54,0.,0.N,I 2,0.66,0.,0.N,13,0.78,0.,0.N, 14,0.82,0,0.N,15,0.86,0.,0.N,16,0.88,0.,0.N, 17,0.90,0.,0.N, 18,0.91,0 .,0.N,I sion5.0- A1-May,1995
NGEN,2,72,196,206, 0,0.92,3.75,0.N,3 I,1.00,3.75,0.NGEN,2,72,22,3 I,I,,,0.002NGEN,2,72,94,103, I,,,0.008NGEN,2,72,166,175, I,,,0.0255NGEN,2,72,238,247, I,,,0.0355NGEN,2,30,2,2 I,1,7.5,NGEN,2,72,32,5 I,I,,,0.002NGEN,2,72,104,123, ,267, I,,,0.0355NGEN,2,10,53,62, I,,7.5,NGEN,2,72,63,72, I,,,0.005NGEN,2,72,135,144, I,,,0.01675NGEN,2,72,207,216, I,,,0.03025mot, 1,339,1TYPE,1MAT,300E, 3,113,111EMORE, 3,41,83,85,115,113EMORE, 12,27,42,26,84,99,114,98EMORE,57,58,68,67E, 13,15,45,43,85,87,117,115EMORE, 14,28,44,27,86,100,116,99EMORE,58,59,69,68E,I 5,17,47,45,87,89,119,117EMORE, 16,29,46,28,88,101,118,100EMORE,59,60,70,69E,l 7,19,49,47,89,91,121,119EMORE, ,51,49,91,93,123,121EMORE,20,31,503Version5. 00,92,103,122,102- A2-May,1995
E,197,198,208,207E,149,151,1 201,211,210E,155,157,187,185,227,229,259,257EMORE, 11E,157,159,189,187,229,231,261,259EMORE, 12E,159,161,19I,I89,231,233,263,261EMORE, 13E,16 I,163,193,191,233,235,265,263Version5. 0- A3-May,1995
EMORE, 14E,163,165,195,193,235,237,267,265EMORE, ,31 RE,276,277,2 47,266,246,308,319,338,31 E,2,2,20WMORE,32,32WMORE,32,32,1WSTART,3,3 ,2,2,30WMORE,4,4,30WMORE,22,22, IWSTART,5,5,2,2,30Version5.0- A4-May,1995
ORE,25,25,IWSTART, 11,11,2,2,30WMORE, 12,12,30WMORE,26,26,1WSTART, I 3,13,2,2,30WMORE,14,14,30WMORE,27,27,1WSTART, 15,15,2,2,30WMORE,16,16,30WMORE,28,28,1WSTART, 17,17,2,2,30WMORE, 18,18,30WMORE,
element type indicator have to be input. If shell elements are used in the analysis, the indicator has to be set to 1. Information pertaining to the element level is read in a loop over the number of elements within the element group considered (format 415,3F15.0). The element level
