一种考虑非理想界面的颗粒增强复合材料微观结构的几何建模方法技术

本发明专利技术公开了一种考虑非理想界面的颗粒增强复合材料微观结构的几何建模方法，属于颗粒增强复合材料结构分析领域。基于UG/OPEN API二次开发，在三维代表性体积元（Representative Volume Element，RVE）中生成多面体颗粒增强相和非理想界面相，包括：通过在随机点之间构造剪切平面，对立方体进行切割，生成随机多面体；然后判断该多面体是否发生几何干涉，若干涉则选取干涉平面切割多面体，保留非干涉部分，否则保留该多面体；接着，对多面体进行几何抽壳，生成界面相，将界面相与RVE进行布尔求差并保留界面相；重复此过程，直至增强相和界面相的总体积比达到指定值。本发明专利技术设计思路简单，易于编程实现，几何模型考虑了微观组织界面相，可直接导入有限元软件，便于后续有限元分析。

【技术实现步骤摘要】
一种考虑非理想界面的颗粒增强复合材料微观结构的几何建模方法
本专利技术涉及颗粒增强复合材料结构分析领域，具体地说是一种在RVE基体中建立互不干涉的随机增强体颗粒和非理想界面相的一种几何建模方法，利用UG/OPENAPI二次开发构建颗粒增强复合材料微观结构几何模型。
技术介绍
颗粒增强复合材料是一种以金属或非金属为基体，陶瓷、金属或碳素材料等作为增强相复合制备而成的复合材料。颗粒增强复合材料综合了基体和增强体的性能，以金属为基的复合材料具有高强度、高弹性模量、耐磨损和导电导热性能好等特点，广泛应用于航空航天、电子、汽车及建筑等行业，而且具有无法替代的优势。由于金属基复合材料的广泛应用，以及计算机技术的发展，以有限元为代表的数值模拟方法逐步成为研究复合材料的有力工具。有限元模拟的方法是对材料宏观力学性能和材料微观结构、组织成份之间定量关系进行分析，模拟出建立复合材料模型的宏观力学性能。该分析方法通常先建立代表复合材料细观结构的体积单元，通过模拟该单元在一定载荷和边界条件下的力学行为，分析材料组织结构对宏观性能的作用。由于平面结构的有限元模型近似处理过于粗糙，模拟结果偏差较大，本专利仅对颗粒增强金属基复合材料的三维有限元模型进行分析讨论。在构建代表体元模型时一般有两类方法：一是基于复合材料的局部真实结构来构建模型，称为真实结构模型法；二是通过提取复合材料结构的典型特征模拟出复合材料的结构，称为虚拟结构模型法。3D真实结构模型是基于颗粒增强复合材料的真实结构建立的模型。一般使用连续切片法建立，其原理是将样品分成很多层，取得每一层的平面结构图像输入计算机，然后利用三维重构技术，还原出样品三维真实结构。样品分的层数越多，三维结构也越接近真实。但3D重构技术工作量大、成本高，且为降低局部性影响，必须使模型包含足够多的颗粒，颗粒越多，样品分层越多，工作量越大。虚拟结构模型是基于复合材料的结构特征构建的，虽然不反映真实结构，但模型结构具有一定的代表性，且可以通过周期性边界条件和模型的对称性，大量减少模拟计算量。虚拟结构模型主要分为单颗粒单胞模型和多颗粒立方单胞模型，单颗粒单胞模型反映的是均匀规则分布的颗粒增强复合材料，单颗粒单胞模型不能体现颗粒间的相互作用，而且由于所做近似较多，模拟结果的准确性较差。多颗粒立方单胞模型在一个立方体单胞内可以包含多个增强颗粒，颗粒分布可以是规则的也可以随机分布。对于颗粒分布有一定规律的模型，可以在确定颗粒的位置坐标后进行生成；随机分布的颗粒，一般可以采用随机序列吸附方法来生成。SuY等人的随机顺序吸附算法是通过向RVE（代表性体积元）所占据的几何空间中逐一投放颗粒,并判断当前颗粒是否与已存在颗粒之间发生干涉,如不干涉则保存当前颗粒的相关数据,投放终止条件是达到体积分数要求或循环次数达到最大设定值。随机顺序吸附法的缺点是难以生成增强体体积含量较高的单胞模型，很难实现30%体积分数以上的复合材料模型，因为增强体达到一定数量后，剩余的空隙难以容纳一个完整的新增强体颗粒，而已存在的增强体又不能移动位置，所以容易达到一个饱和状态；而且新生成的颗粒很容易与已生成的颗粒发生干涉，导致计算量较大，工作效率低。根据随机顺序吸附法的缺点，ZhangJ等提出下落堆积法，将模型高度提升到原来的n倍，则模型体积也提高到原来的n倍，颗粒随机分布在n倍高的模型中，因为体积的增大，颗粒间相交的可能性大大降低。通过相交判断算法保证任何颗粒都不相交后，对颗粒施加自由落体运动，模型的外表面构成一个容器，颗粒落下后在模型底部开始堆积。通过对颗粒施加一个虚拟厚度，可以控制颗粒间的最小间距，实现颗粒在空间的均匀分布。这种方法虽然可以实现较大体积分数模型建立，但是该方法建立颗粒增强体的形状为圆球、圆柱等规则形状，未考虑是否与实际材料的相符。虽然与随机顺序吸附法相比大大减少了判断颗粒相交的计算量，但与本专利相比计算时间依旧较长。Sheng等采用了多面体颗粒生成法。先生成一个球形，然后在球面上随机取点作为多面体的顶点，取到一定的点后将点连起来构成了一个凸多面体，Su等实际观察到的颗粒尺寸和形状特点，对颗粒形状进行了虚拟重构，他们所生成的颗粒也非常接近真实形状。但是真实颗粒增强金属基复合材料一般由颗粒相、基体相和界面相三部分组成，他们与本专利相比未考虑界面相的建立。而且本专利中界面相可以直接由颗粒增强体抽壳建立，无需借助其他软件，操作简单，大大减少了界面相模型建立的工作难度。本专利技术运用多颗粒随机分布代表性体积元模型，对传统随机顺序吸附算法中发生干涉的颗粒多面体不是直接舍去，而是对其进行几何切割，保留其非干涉部分，大大减少了判断颗粒相交的计算量，可以获得较高的体积比模型；进一步通过对颗粒多面体进行几何抽壳，可以建立有厚度的非理想界面相，无需借助其他软件，大大减小了整体建模难度。本专利技术采用C语言编程，利用UG/OPENAPI二次开发，通过简单的几何切割，建立颗粒增强复合材料微观结构几何模型，方法简单直观，易于编程实现。生成的几何模型考虑了有厚度的非理想界面相，可以在后续的复合材料力学分析中，研究界面相对复合材料力学特性的影响。相比于传统随机顺序吸附算法，该方法可以获得较大体积比的RVE几何模型，几何模型可以直接导入到ABAQUS、ANASYS等有限元软件中进行仿真分析。
技术实现思路
为了提供一种有效降低建模难度，减少计算量，提升工作效率的考虑非理想界面相的颗粒增强复合材料建模方法，本专利技术利用UG/OPENAPI二次开发，在UG三维建模软件中实现颗粒增强复合材料微观几何模型的建立。本专利技术采用技术方案如下：1.一种考虑非理想界面的颗粒增强复合材料微观结构的几何建模方法，其特征在于：该方法包括如下步骤：1）根据指定的复合材料细观结构代表性体积元（RVE）的边长LR，生成复合材料基体的RVE立方体Q，并计算其体积VR=LR*LR*LR；2）在基体的RVE立方体Q内随机生成种子点i，坐标值为（xi，yi，zi），判断种子点i是否位于已生成的随机多面体（P1，P2，…，Pi-1）内部，若在内部则重新生成种子点i，若种子点i位于所有已生成随机多面体外部，执行下述步骤，随机生成边长Li,以种子点i为中心，边长为Li，生成种子立方体Ci；3）在种子立方体Ci中，随机生成n个散乱点，然后在该种子立方体Ci内利用种子点i，分别与n个散乱点之间分别构造两点间的中心对称面，作为该种子点的切割平面R，利用切割平面R对该种子立方体Ci进行几何切割，得到随机多面体Pi；4）依次判断随机多面体Pi与已生成随机多面体（P1，P2，…，Pi-1）是否发生干涉，若发生干涉，则在已生成多面体Pi-1上选取与随机多面体Pi发生干涉的，并与种子点i距离最近的面作为切割平面M，对随机多面体Pi进行几何切割，并保留非干涉部分；若不发生干涉，则直接保留随机多面体Pi；5）通过UG/OPENAPI编程得到每个随机多面体Pi体积Vi，计算体积比w=（V1+V2+……+Vi）/VR，根据本文档来自技高网...

【技术保护点】
1.一种考虑非理想界面的颗粒增强复合材料微观结构的几何建模方法，其特征在于：该方法包括如下步骤：/n1）根据指定的复合材料细观结构代表性体积元（RVE）的边长L

【技术特征摘要】
1.一种考虑非理想界面的颗粒增强复合材料微观结构的几何建模方法，其特征在于：该方法包括如下步骤：
1）根据指定的复合材料细观结构代表性体积元（RVE）的边长LR，生成复合材料基体的RVE立方体Q，并计算其体积VR=LR*LR*LR；
2）在基体的RVE立方体Q内随机生成种子点i，坐标值为（xi，yi，zi），判断种子点i是否位于已生成的随机多面体（P1，P2，…，Pi-1）内部，若在内部则重新生成种子点i，若种子点i位于所有已生成随机多面体外部，执行下述步骤，随机生成边长Li，以种子点i为中心，边长为Li，生成种子立方体Ci；
3）在种子立方体Ci中，随机生成n个散乱点，然后在该种子立方体Ci内利用种子点i，分别与n个散乱点之间分别构造两点间的中心对称面，作为该种子点的切割平面R，利用切割平面R对该种子立方体Ci进行几何切割，得到随机多面体Pi；
4）依次判断随机多面体Pi与已生成随机多面体（P1，P2，…，Pi-1）是否发生干涉，若发生干涉，则在该已生成多面体上选取与随机多面体Pi发生干涉的，并与种子点i距离最近的面作为切割平面M，对随机多面体Pi进行几何切割，并保留非干涉部分；若不发生干涉，则直接保留随机多面体Pi；
5）通过UG/OPENAPI编程得到每个随机多面体Pi体积Vi，计算体积比w=（V1+V2+……+Vi）/VR，根据事先指定的体积比标准w0，判断体积比w是否超过w0，若w<wo,则重新执行步骤2）3）4）生成一个新的随机多面体，并重新...

【专利技术属性】
技术研发人员：郑兴帅，孙腾飞，周吉星，尚晨阳，明平美，秦歌，闫亮，张新民，牛屾，李欣潮，张振，
申请(专利权)人：河南理工大学，
类型：发明
国别省市：河南;41