Python-Matlab联合编程Abaqus高级后处理技术研究

367201572014－11－25513280501011978—1957—。【】doi10．11809/scbgxb2015．07．034Python-MatlabAbaqus264001、AbaqusPython-MatlabAbaqusAbaqusAbaqus。Python-MatlabAbaqus．Python-MatlabAbaqusJ．20157133－138．CitationformatRENHai-fengGAOMing．ResearchonAbaqusAdvancedPostProcessingTechnologyBasedonPython-MatlabProgrammingJ．JournalofSichuanOrdnance20157133－138．V435．1A1006－0707201507－0133－06ResearchonAbaqusAdvancedPostProcessingTechnologyBasedonPython-MatlabProgrammingRENHai-fengGAOMingNavalAeronauticalandAstronauticalUniversityYantai264001ChinaAbstractThefurtherapplicationsoffiniteelementanalysisputforwardhighrequesttoflexibilityandcus-tomizationmodeofresultoutput．ForitisdifficulttooutputnodestressinbatcheswithAbaqus/CAEthecustomizationoutputmethodswereproposedutilizingprogramsofPython-Matlabbyreconfigurationcalcula-tionandcallingkernelofAbaqus．Themasscustomizationoutputofnodestressforthesimulationresultsofbondstresssensorusinginhealthmonitoringofsolidrocketmotorwererealizedbyproposedmethods．Theoutputresultsindicatethatthesuggestedmethodsarecorrectquiteeffectivehigherprecisiongoodflexi-bilityandunderstandability．KeywordssolidrocketmotorbondstresssensorPython-MatlabProgrammingoutputinbatchesad-vancedpostprocessingofAbaqusreconfigurationcalculationkernelcallAbaqus、Ansys。。。Abaqus1－3。。、。ChaoXingAbaqus。Abaqus。AbaqusPython“”Python、、、PythonAbaqusAbaqus/CAE、、24。、、PythonAbaqus。1PythonAbaqusAbaqus。Python、、ExcelMatlab。PythonAbaqusExcelMatlab。。1EPart－1Set-1Abaqus\CAEPartTools→Set→CreateSetSet-1type=Geome-try。odbJob-31．OdbStep-1C3D8RSet-1。1ESet-12Abaqus5。。16Abaqus1。Abaqus、、Mises。。2．1AbaqusAbaqus。AbaqusSESSIONODBMDB17－10ODBMDB2．1．1AbaqusODB2ODBModelDataResultDatarootAssenbly、parts、sec-tionCategories、materialsstepsstep、frame、historyoutputsfieldoutputs。odb．setps．frames．fieldOutputs—'S'—'E'—'U'odb．setps．historyRegions．historyOutputs2．1．2AbaqusMDB3MDBModelDataJobDatamodelparts、materials。partselements、setsnodeslabel、coordinate、connectivity。mdb．models．parts．setsodb．rootAssembly．instances．elementSets、mdb．models．parts．sets．elements．elementLabelodb．rootAssembly．instances．elementSets．elements．elementLabelmdb．models．parts．sets．elements．connectivityodb．rootAssembly．instances．elementSets．elements．connectivitymdb．models．parts．sets．nodes．coordinates431http//scbg．qks．cqut．edu．cn/ChaoXing2ODB3MDB2．2AbaqusPythonAbaqusJob－31．odbSet－1Scomponent．txt。Set－1elmementofset．txt。Set－1node-ofset．txt。1、#odbAccessfromodbAccessimport*#odbmyodb=openOdb'Job－31．odb'#Step－1frames－1stressfield=myodb．steps'Step－1'．frames－1．fieldOut-puts'S'#Set－1skinset=myodb．rootAssembly．instances'PART－1－1'．ele-mentSets'SET－1'#Set－1field1=stressfield．getSubsetregion=skinsetposition=INTE-GRATION_POINTelementType='C3D8R'#Set－1val=field1．values#Scomponent．txtSet－1okScomponent．txtwithopen'D/Scomponent．txt''w'asf1forvinvalf1．writelinesstrv．elementLabel+''f1．writelinesstrv．data+'\n'elseprint'ok'f1．close##elmementofset．txtSet－1#ok2elme-mentofset．txtwithopen'D/elmementofset．txt''w'asf2forvinvalelementnode=myodb．rootAssembly．instances'PART－1－1'．elementsv．elementLabel－1．connectivityf2．writelinesstrv．elementLabel+''f2．writelinesstrelementnode+'\n'elsef2．closeprint'ok2'#Set－1#nodeofset．txtok3nodeofset．txtnd=mdb．models'Model－1'．parts'Part－1'．sets'Set－1'．nodeswithopen'D/nodeofset．txt''w'asf3v=1whilev＜=lenndnodescoors=mdb．models'Model－1'．parts'Part－1'．sets'Set－1'．nodesv－1．coordinatesnodeslabel=mdb．models'Model－1'．parts'Part－1'．sets'Set－1'．nodesv－1．labelf3．writelinesstrnodeslabel+''f3．writelinesstrnodescoors+'\n'v=v+1else531Python-MatlabAbaqusChaoXingf3．closeprint'ok3'#odbmyodb．close2．3PythonAbaqusMatlabExcel。ExcelPython*．txt。Excel。Excel*．xls*．txtMatlab。2．4AbaqusS11S22S33S12S13S23S11S22S33xxyyzzS12xyS13xzS23yzxyxzyzyxzxzy16xyzxyzS11S12S13S12S22S23S13S23S[]331σ1σ2σ3VonMises6σ=12σ1－σ22+σ2－σ32+σ3－σ12槡22clearclc%loadD\report．txt%fori=11lengthScomponent%s11=Scomponenti2s22=Scomponenti3s33=Scomponenti4s12=Scomponenti5s13=Scomponenti6s23=Scomponenti7A=s11s12s13s12s22s23s13s23s33%ps=eigA%mpofseti1=Scomponenti1%mpofsetmpofseti2=ps1%mpofsetmpofseti3=ps2mpofseti4=ps3mpofseti5=maxps%mpof-setmpofseti6=sqrt0．5*ps1－ps2^2+ps2－ps3^2+ps3－ps1^2%misesmpofsetend2．5Abaqus。。C3D8R13。3load'D\elmementofset．txt'%load'D\nodesofset．txt'%elmementofset1=%nodestress=%、fori=11lengthnodesofset%kn=findelmementofset==nodesofseti1%sum=0sum2=0forij=11lengthk%sum=sum+mpofsetkij15%sum2=sum2+mpofsetkij16%Misesendnodestress=nodestressnodesofseti14sum/lengthksum2/lengthk%End3AbaqusAbaqus。Abaqus\CAEvisualizationnodelistPathPathXYDatalabel。XYDataX。PythonAbaqusAbaqus。3．1AbaqusSessionAbaqusSession4session、。Abaqus\CAEVisualizationSession631http//scbg．qks．cqut．edu．cn/ChaoXingAbaqus。4SESSION3．2AbaqusAbaqusPythonsessionnodelistpathsession．PathpathxyDataObjectsession．XYDataFromPath。4Abaqus#Job－31．odbvisualization。fromodbAccessimport*fromabaqusConstantsimport*myodb=openOdb'Job－31．