python参数化建模 书_基于python三维复合材料层压板参数化建模6-28.docx

基于python三维复合材料层压板参数化建模6-28.docx

基于ABAQUS/PYTHON的三维复合材料层压板有限元建模参数化开发伊鹏跃，于哲峰，汪海(上海交通大学航空航天学院，上海200240)摘要讨论了ABAQUS对象模型及网格对象的调用方法。利用PYTHON脚本，开发了可应用连续壳或实体单元的三维复合材料层压板参数化建模程序。通过了人机交互，可选择层间是否使用界面单元，材料是否预制损伤。实例表明，该程序既可完成单个层压板及损伤的参数化三维建模，也可应用于装配形成的多板连接结构的建模，适用范围广泛。关键词ABAQUS/PYTHON；三维；复合材料；层压板；参数化建模中图分类号TP319TB332PARAMETRICTHREEDIMENSIONALMODELINGOFCOMPOSITELAMINATESBASEDONABAQUS/PYTHONYIPENGYUE,YUZHEFENG,WANGHAISCHOOLOFAERONAUTICSANDASTRONAUTICS,SHANGHAIJIAOTONGUNIVERSITY,SHANGHAI200240ABSTRACTABAQUSOBJECTMODELANDINVOKINGOFMESHOBJECTWEREDISCUSSEDAPARAMETRICMODELINGPROGRAMFORTHREEDIMENSIONALCOMPOSITELAMINATEWASDEVELOPEDINWHICHCONTINUUMSHELLORSOLIDELEMENTISUSEDHUMANCOMPUTERINTERACTIONWASREALIZEDBYWHICHLAMINATEWITHINTERFACEANDDAMAGEORNOTCANBESELECTEDEXAMPLESSHOWTHATTHEPROGRAMISABLETOBEUSEDNOTONLYINTHEMODELINGOFONELAMINATEANDDAMAGEBUTALSOINTHESTRUCTURESCONSISTINGOFMULTIPLEPLATESBYASSEMBLINGTHEAPPLICATIONREGIONOFTHEPROGRAMISWIDEKEYWORDSABAQUS/PYTHONTHREEDIMENSIONCOMPOSITELAMINATEPARAMETRICMODELING引言由于复合材料高比强度、比刚度及较强可设计性等特点在航空航天以及民用结构中的用量越来越大1。层压板式纤维增强复合材料在结构中应用广泛，其力学性能成为研究重点2。由于材料各向异性，层间性能差，损伤机理复杂，凡考虑层压板厚度、厚度方向应力分布或层间作用时，基本都需要使用三维模型。ABAQUS具有丰富的材料本构和单元类型，并具备强大的二次开发能力，是复合材料结构分析的有力工具。基于ABAQUS有限元分析软件，贾建东等3通过不含界面层的三维复合材料层压板模型利用UMAT编写材料本构，研究了层压板的冲击后剩余压缩强度；BGFALZON等4等基于三维实体模型研究了层压板的层内渐进损伤并进行了算例验证；顾亦磊等5利用三维层压板模型使用USDFLD程序研究了层压板螺栓连接失效；VOLNEITITA等6利用三维模型对层压板静压痕试验进行了有限元仿真研究；张丽等7建立了含界面层的层压板三维有限元模型，利用VUMAT子程序进行了复合材料层压板低速冲击作用下损伤分析；滕锦等8通过含粘性层的三维层压板模型研究了增韧复合材料的冲击损伤过程，朱炜垚等9仅在铺层角度不同的层间应用界面层简化三维模型，研究了层压板低速冲击下分层损伤；CRAVEN等1011建立了含多处分层损伤与纤维损伤的层压板三维有限元模型，研究了带损伤层压板的拉伸、屈曲及后屈曲力学性能；王跃全等12利用层压板修补的三维有限元模型，进行了胶接贴补修理渐进损伤分析；张彦等13基于三维复合材料层压板模型对双悬臂梁试验进行仿真，验证了分层的临界能量释放率。通过ABAQUS的前处理模块进行三维层压板有限元建模是一个相当重复繁琐且耗时耗力的工作。如果开发参数化建模程序则可大大提高建模效率，促进层压板材料的性能分析和研究。本文采用PYTHON脚本语言，基于孤立网格及单元节点信息，通过查找相应单元进行网格偏置，建立三维有限元网格，并赋予单元类型、材料属性与铺层角度等。开发完成了可包含界面层及可预制损伤的复合材料层压板三维参数化建模程序，应用了人机交互界面，对ABAQUS对象模型、主要模块程序流程及实现方法进行介绍，并通过层压板三维建模例子演示程序的功能。1层压板建模与PYTHON脚本三维复合材料层压板的建模繁琐，以基于网格偏置的建模方法为例，首先建立单层的几何模型，划分网格后将部件转化为孤立网格部件，进入MESH模块，通过MESHEDIT功能对单层网格进行偏置，生成界面层或第二层复合材料单层。若为界面层，厚度值赋零，根据建模要求逐层偏置生成复合材料层压板三维模型，由于网格部件不具备几何信息，需为每一界面层与复合材料单层建立集合，利用集合为相应单元赋单元类型，材料属性，铺层角度等。如果铺层层数较多，通过用户界面进行操作，是一个重复工作量大且耗时耗力的过程。PYTHON脚本接口是ABAQUS的二次开发接口之一，它扩充了PYTHON的对象模型和数据类型。一般情况下，脚本接口主要用于前处理、后处理、自定义模块等。若前处理通过用户界面建模，需要大量的手动操作,而PYTHON脚本使用区区几十行代码则可实现上百次操作的效果，若以此开发复合材料层压板参数化建模程序，可大大减少工作量，提高效率。ABAQUS/PYTHON主要有SESSION对象、MDB对象和ODB对象，本文主要对MDB模型对象进行搜索调用。MDB对象保存于模型数据库中，调用前需使用IMPORT语句导入，图1为MDB对象结构。通过对PART对象的查询，给出了程序需用到的相关子对象，其中ELEMENTS下的LABEL、CONNECTIVITY分别为单元编号、对应的节点编号，NODES下的COORDINATES、LABEL分别为节点坐标、节点编号。图1MDB对象结构FIG1MDBOBJECTMODELTREE2含损伤的三维层压板参数化建模本程序包括层压板的参数化建模部分和损伤建模部分，流程分别如图2A和图2B所示。对于层压板的参数化建模，首先根据面内对网格的要求，例如是否需要局部加密等，可通过程序选择提前划分网格的平面部件，或者通过程序基于草图模块建立网格均匀的平面部件。网格偏置建模，需将平面部件转化为网格部件，然后根据层间是否引入界面层建立模型，包括不含界面层的层压板模型、各层间都有界面层的模型以及指定层间加入界面层的模型。由于模型中只含网格信息，为方便相应对象的选择，需要通过搜索单元编号为各层包括界面层建立单元集合。建立材料后，基于单元集合对相应的区域赋单元类型，对于铺层单元，可选择使用连续壳单元或实体单元，界面层单元使用COHESIVE单元，对于不同的分析步类型，可分别选择STANDARD和EXPLICIT单元，然后建立对应的截面属性并赋予相应的单元区域，最后根据提示对铺层角度进行赋值，至此三维层压板的建模完成。对于损伤预制部分，建立损伤材料及截面属性后，会提示预制损伤的位置，可选择某一铺层或者界面层，然后需选择