Python语言和ABAQUS后处理二次开发定稿

1、Python语言和ABAQUS后处理二次开发11001ABAQUS和二次开发后的计算结果，有效地解决了提取大量结果进行重Python语言二次开发ABAQUS后处理和开发PythonABAQUS问题，提高了后处理的效率。ABAQUS后处理二次开发PythonABAQUSABAQUSA后处理ABAQUSPython脚本语言,它应用程序开发1。开发PythonABAQUS1二次开发ABAQUS脚本接口是一个函数专门用于处理对象中的数建后，可以使用该对象提供的方法来处理对象ABAQUSAPI，Ap后处理等功能。接口编程采用Python语件之一,可以模拟绝大部分工程材料的线性和,也是建模、分析和后处理的

2、人机交互平台,可以独立运行,也可以用作脚本语言。特别适后的计算结果，有效地解决了提取大量结果进数据成员和函数。对象tonPoABAQUSAPython语言之间的层次结构和相互关系。ontn对on和o象包括计算模型对象和session；一mdbodtiviewportsjobsAssemblypartsfieldReportOptioris*=Ojntainersectionssketches=Singulagobjecit一parts+materialsModelDatasteps图1一framesi次和c和数据fieldOutputsrootAssemblyontsinsrABAQUS后处i

3、理二次开发sectionCategorie-SingularobjecthltyRegia.KtS一pointABAQUS后处理即的二次开发ABAQUSo和让ryOutputsABAQUS*o和tResultsD（10472）0、0大3连理GZ08)1、验室开放基金（30715）0、助10。ytnABAQUSPythonABAQUSPythonABAQUSABAQUSABAQUS四几种途径1：1用户子程序接口，通过用户子ABAQUS2通过环境初ABAQUS多缺省设置；3脚本ABAQUSScriptingInterface算结果；4图形用户接ABAQUSGUIToolkitGUI本文用到的是第3

4、种方法,通过PythonABAQUS12000Hz性的5个典型点的法向加1-2000HzAbaqus0-2000HzABAQUS取5个典NodePART-1-1.823、NodePART-1-1.541NodePART-1-1.1175NodePART-1-1.698NodePART-1-1.1388分别对应15号节点，如图3所示，同时确定各典：1号点的法x方向（水平方向），23、4点在侧面上，法向为与侧面垂直的方向，5号y,StepHistoryOutputABAQUSCAEinprandom_restart.odbABAQUSDOS2)提取拱板固有频率和振mdbviewpointodbno

5、k4库创建odbPath=gongban_damp.odbmyOdb=visualization.openOdb(path=odbPath)myViewport.setValues(displayedObject=myOdb)Frequency.txtE:houchuliFrequency.txt路径下：xyp=session.XYPlot(name=XYPlot-1)odb=session.odbsgongban_damp.odbxy1=session.XYDataFromHistory(name=Temp-1,odb=odb,outputVariableName=Eigenfrequenc

6、y:EIGFReQforWholeModel,steps=(fre,)xyp.Curve(name=Temp-1,xyData=xy1)xyp.setValues(curvesToPlot=(Temp-1,)myViewport.setValues(displayedObject=xyp)x0=session.xyPlotsXYPlot-1.curvesTemp-1session.writeXYReport(fileName=fre.txt,xyData=(x0,)getlnputsgetlnputs104户输入值。然后通过whileModei.pngnumMode=NonewhilenumM

7、ode=None:图4.提示输入需要提取振型图的数量getlnputsNonewhile输入数量的模态振型图，并以PNG格式输出E:houchulinumMode=getInput(EnteranumberforModes,10)numMode=int(numMode)i=1whilei=max:max=y1Fre=x1return1号点A1和A2OneData(A0)maxValue(x1,y1)defOneData(A0):forfrequency1,acceleration1inA0:printfile,%f%f%(frequency1,acceleration1)maxValue(fr

8、equency1,acceleration1)TwoData(A1,A2)2号、3号和4号点的法向加速度为A1和A2maxValue(x1,y1)defTwoData(A1,A2):x=pi/4.0i=0whileifile,%f%f%(frequency2,acceleration2)maxValue(frequency2,acceleration2)i=i+1odbStep-2Acceleration.txtodb=openOdb(path=random_restart.odb)step2=odb.stepsStep-2file=open(Acceleration.txt,w)forreg

9、ionNameinstep2.historyRegions.keys():region=step2.historyRegionsregionNameA1Data=region.historyOutputsA1.dataA2Data=region.historyOutputsA2.dataFre=0max=0ifregionName=NodePART-1-1.823:printfile,ThefrequencyandaccelerationforNode823inNODE1:OneData(A1Data)printfile,ThemaxAccelerationofNODE1is%fatthefr

10、equencyvalue%f%(max,Fre)elifregionName=NodePART-1-1.541:printfile,ThefrequencyandaccelerationforNode541inNODE2:TwoData(A1Data,A2Data)printfile,ThemaxAccelerationofNODE2is%fatthefrequencyvalue%f%(max,Fre)elifregionName=NodePART-1-1.1175:printfile,ThefrequencyandaccelerationforNode1175inNODE3:TwoData(

11、A1Data,A2Data)printfile,ThemaxAccelerationofNODE3is%fatthefrequencyvalue%f%(max,Fre)elifregionName=NodePART-1-1.698:printfile,ThefrequencyandaccelerationforNode698inNODE4:TwoData(A1Data,A2Data)printfile,ThemaxAccelerationofNODE4is%fatthefrequencyvalue%f%(max,Fre)elifregionName=NodePART-1-1.1388:prin

12、tfile,ThefrequencyandaccelerationforNode1388inNODE5:OneData(A2Data)printfile,ThemaxAccelerationofNODE5is%fatthefrequencyvalue%f%(max,Fre)file.close()odb.close()E:houchuliAcceleration.txtNode,maxAcceleration,FrequencyNODE823,0.132961168885,204.81187439,NODE541,6.54942031394,121.438873291,NODE1175,1Wey-hn,.Python().社,2008.2Abaqus6.8HTMLDocumentation/Documentation/DocumentatisonM/aAnbuaqlu.sScriptingU3宋岩新,杨庆,唐小微,万少石.ABAQ