基于Abaqus的有限元渐进损伤强度分析_孔祥宏

2910Vol.29No．10201210ComputerApplicationsandSoftwareOct．2012Abaqus2100162011－12－22。。。。Hill-TsaiTsai-Wu。Abaqus/CAEPython。。、。TP301ADOI10．3969/j．issn．1000-386x．2012．10．062FINITEELEMENTINCREMENTALDAMAGESTRENGTHANALYSISBASEDONABAQUSKongXianghongWangZhijinCollegeofAerospaceEngineeringNanjingUniversityofAeronauticsandAstronauticsNanjing210016JiangsuChinaAbstractTakingthedamagedcompositesandwichpanelwithfoamcoreasanexampleweintroducetheapplicationoffiniteelementin-crementaldamagestrengthanalysisinstructuralstrengthcalculation．ForcompositeskinsHill-TsaistrengththeoryandTsai-Wutensortheoryareadoptedasthestrengthjudgementbasisandthemethodofmaterialperformancedegradationisusedinstiffnessreductionofthedamagedelementsinfiniteelementmodel．ParametricmodellingisimplementedusingAbaqus/CAE．Pythonprogramisutilisedtocontrolwholeprocessofincrementaldamageanalysis．Theresultsoffiniteelementincrementaldamagestrengthanalysisagreeverywellwiththeresultsoftestthroughtestingandnumericalexamplesverification．Consequentlytheincrementaldamageanalysismethodinconjunctionwithsuitablestrengththeoryandmaterialperformancedegradationmethodcanbeusedtothestrengthanalysesofcomplexstructure．KeywordsCompositesandwichpanelwithfoamcoreIncrementaldamageStrengthanalysisMaterialperformancedegradationStiff-nessreductionParametricmodelling0。。AbaqusPython。1、。1、2、3Abaqus。1．11。110Abaqus2371．2、、Hill-Tsai、Tsai-Wu。。1．2．1Hill-TsaiX、Y、SHill-Tsai1σ11X2－σ1σ2X2+σ22Y2+τ212S2=111X1Y2S1－2。Hill-TsaiX、Y、S。Hill-Tsai、。、Hill-Tsai1X、Y2X=XTσ1＞0XCσ1＜0{Y=YTσ2＞0YCσ2＜0{2XT1XC1YT2YC2。1．2．2Tsai-WuTsai-WuFiσi+Fijσiσj+Fijkσiσjσk+…=1ijk…=12…633σiσjσkFiFijFijkσ4=τ23σ5=τ31σ6=τ12。σ=σ1σ2σ3τ23τ31τ12T43F1σ1+F2σ2+F3σ3+F11σ21+F22σ22+F33σ23+2F12σ1σ2+2F23σ2σ3+2F31σ3σ1+F44σ24+F55σ25+F66σ26=15512F1=1XC－1XTF2=1YC－1YTF3=1ZC－1ZTF11=1XTXCF22=1YTYCF33=1ZTZCF44=1S223F55=1S231F66=1S212F12=－12XTXCYTY槡CF23=－12YTYCZTZ槡CF31=－12ZTZCXTX槡C6σ3=τ23=τ31=05F1σ1+F2σ2+F11σ21+F22σ22++2F12σ1σ2+F66σ26=1776。1．3。23。2E=0．01EE=0．01E31σ1＞XTE11=E22=E33=0G12=G13=G23=0v12=v13=v23=01－σ1＞XCE11=E22=E33=0G12=G13=G23=0v12=v13=v23=02σ2＞YTE22=0．01E22G12=0．2G12G13=0．2G132－σ2＞YCE22=0．01E22G12=0．2G12G13=0．2G133σ3＞ZTE33=0．01E33G12=0．2G12G13=0．2G133－σ3＞ZCE33=0．01E33G12=0．2G12G13=0．2G131－2|τ12|＞S12E22=0．01E22G12=0．01G121－3|τ13|＞S13E33=0．01E33G13=0．01G132－3|τ23|＞S23E22=0．01E22E33=0．01E33G12=0．01G12G13=0．01G13G23=0．01G232。、。2382012、、。、、。2．1。。。、、。AbaqusPython。1。1PythonAbaqus、、、、、、。、。PythonAbaqus．odb。。。0．7。。2．22．2．1Python、。、、、、、、、、、、、。2。2、、CompMdCompMdCompMdobjectParentclassCompMd。。2CompMd、。myCmpMddmdp。crtElemMarkEmMkEmMk、。initProperty。myPropertycgEmSetsEmSets、、、。saveMdbnmnm。clsMdbAbaquscae。newMdbModel-1。10Abaqus239myJobnmnm。reLoadldld。getVal．odbEmValEmVal。myCrACScrEmValcrEmVal。2．2．2Python。3、、CompMd。3CrACS。CrA、CrC、CrS、、。MJname．cae、．odb、．txt。wrtFilerpt1rpt2nufmxfrePt。AnsCompSWHdmpsplprpd、、、。34180mm120mm9．6mm。、0．8mm45°/－45°s0．2mm8mm。。43．1。60mm、8．8mm0．4mm7．5mm75mm60mm60mm。5、6、745°－45°。。。5617224020123．24Ply1、Ply275mmPly3、Ply460mm。X0°Z90°。4±45°。43．35。。5/kN17．8518．55Hill-Tsai18．33Tsai-Wu18．304。Hill-TsaiTsai-Wu。。、。1．Tsai-WuJ．2006341274-76．2．J．20072851371-1375．3．J．2008404517-520．4．J．1998154．5．J．2007253342-346．6．J．2010423335-339．2295Android。2%。8。。。1CarlijnVCBoutenKarelTMKoekkoekMaartenVerduinetal．ATriaxialAccelerometerandPortableDataProcessingUnitfortheAs-sessmentofDailyPhysicalActivityJ．IEEEtransactionsonbiomed-icalengineering19973．2ZhuYiRuanXingyunShengYietal．ResearchfortheRelationshipsbetweenBodyAccelerationandEnergyExpenditureJ．IEEEtransac-tionsonbiomedicalengineering20043．3LiMeinaJangMyoungKimYounTaeKim．ACombinedHeartRateandMovementIndexSensorforEstimatingtheEnergyExpenditureC//IEEESENSORS2010Conference2010．4YuichiFujiki．Makingsenseofaccelerometermeasurementsinperva-sivephysicalactivityapplicationsC//CHI2009StudentResearchCompetition20094．5．MEMSJ．2010．6KangJ．BioenergeticsprimerforexercisescienceM．HumanKinet-ics2008．7．J．2002252．8．J．1994174．9．J．2009291．10．J．2004．11DavidMoore．Areal-worldsystemforhumanmotiondetectionandtrackingD．CaliforniaInstituteofTechnology2003．