CN108984910B - 一种三维四向编织复合材料的参数化建模方法 - Google Patents

一种三维四向编织复合材料的参数化建模方法 Download PDF

Info

Publication number
CN108984910B
CN108984910B CN201810797369.8A CN201810797369A CN108984910B CN 108984910 B CN108984910 B CN 108984910B CN 201810797369 A CN201810797369 A CN 201810797369A CN 108984910 B CN108984910 B CN 108984910B
Authority
CN
China
Prior art keywords
fiber
section
geometry
fibre bundle
matrix
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810797369.8A
Other languages
English (en)
Other versions
CN108984910A (zh
Inventor
胡殿印
王荣桥
杨尧
张龙
曾雨琪
韩江博
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN201810797369.8A priority Critical patent/CN108984910B/zh
Publication of CN108984910A publication Critical patent/CN108984910A/zh
Application granted granted Critical
Publication of CN108984910B publication Critical patent/CN108984910B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Complex Calculations (AREA)
  • Braiding, Manufacturing Of Bobbin-Net Or Lace, And Manufacturing Of Nets By Knotting (AREA)

Abstract

本发明介绍了一种三维四向编织复合材料的参数化建模方法，实现步骤如下：针对三维四向编织复合材料的内部纤维束结构进行立体几何分析，提取几何单胞模型；根据立体几何映射关系，形成单胞内部纤维束与六面体边界的相交截面控制方程；基于Matlab平台调用单胞及其纤维束的空间几何控制方程，编写在空间内生成纤维网格信息和基体网格信息的自动化程序；将网格信息和材料信息输出到可导入Abaqus中进行仿真模拟的inp文件；在程序中将编织角、纤维束轴尺寸、网格尺寸等设为单胞模型参数化变量，以inp文件为媒介工具，实现三维四向编织复合材料的参数化建模。

Description

一种三维四向编织复合材料的参数化建模方法
技术领域
本发明属于复合材料技术领域，特别是小编织角度三维四向复合材料的参数化建模。
背景技术
复合材料因其具有较高的比强度、比刚度，在航空航天领域得到了广泛的应用。其中三维四向编织复合材料较传统的复合材料层合板而言，拥有更为优秀的抗层间分离的性能，在航空航天领域结构减重设计方面的应用前景十分广泛。三维四向编织复合材料在细观尺寸上具有复杂的结构，故不宜在有限元软件中对整体结构进行计算分析。三维四向编织复合材料是具有多尺度、离散分布的多相体，目前对其进行多尺度分析时广泛采用代表性体积单元(RVE)的方法。
传统的三维四向编织复合材料的RVE模型是由纤维段和环绕纤维段的基体组成的平行六面体，其存在以下几方面的缺点：
(1)传统的RVE模型将基体分为了在纤维段内部的和在纤维段之间的两种类型，在纤维段内部的作为横向各向同性材料，纤维段之间的作为各向同性材料。实际上基体在纤维段内部和在纤维段之间都是各向同性材料，没有本质的差别；
(2)一条三维四向编织型复合材料的纤维束拥有成千上万的不同长度的纤维段，导致传统的RVE模型在应用于有限元计算时，生成周期性网格困难，工作量巨大。
因此，目前三维四向编织复合材料的代表性体积单元的建模技术仍需要改进。
发明内容
本发明的技术解决问题是：克服了传统三维四向编织复合材料RVE模型建立时需要将基体进行区分和产生周期性网格的困难，提供了一种建立三维四向编织复合材料干纤维和基体的RVE模型周期性网格的方法，大大减少了人力和物力成本。
本发明的技术解决方案是：一种三维四向编织复合材料的参数化建模方法，其特征在于使用Matlab程序获得RVE模型的inp文件来替代Abaqus的前处理模块。
利用所述方法能够解决传统方法中建立RVE模型时无法产生周期性网格的问题，并且不需要对基体的位置进行区分。计算结果与实验数据的对比表明，在小编织角度条件下，该模型对结果的预测有很高的准确度。
本发明一种三维四向编织复合材料的参数化建模方法，包括以下步骤：
(1)针对三维四向编织复合材料的内部纤维束结构进行立体几何分析，提取六面体结构的几何单胞模型，即RVE模型；
(2)根据立体几何映射关系，形成单胞内部纤维束与六面体边界的相交截面控制方程，具体为根据(1)中进行的几何分析建立了映射关系，将纤维束主截面上的点全部映射到纤维束与椭圆截面上，得到每条纤维束主截面与在六面体边界上形成的两个椭圆截面之间的映射关系，所述映射关系即为相交截面控制方程；
(3)基于Matlab平台调用步骤(2)形成的单胞及其纤维束的空间几何控制方程，编写在空间内生成纤维网格信息和基体网格信息的自动化程序，具体为：
(31)输入编织复合材料所需要的几何和材料参数；
(32)利用“拉丁方设计”的方法，在纤维束的主截面上生成随机分布的代表纤维束中随机分布的纤维的点，纤维的主截面为椭圆面，椭圆面的长轴和短轴的长度分别为由输入的几何和材料参数来决定；
(33)取三维四向编织复合材料单胞的一束纤维，根据步骤(2)中建立的单胞内部纤维束与六面体边界的相交截面控制方程将纤维束主截面上点的坐标映射到RVE模型的两个椭圆截面上；
(34)关联每条纤维段的两个端点，并计算纤维单元的数据，在商业有限元软件Abaqus中采用B31完成，其中B31单元的尺寸由步骤(31)所决定；
(35)基于Matlab平台编写能重复进行(31)-(34)步骤运算的代码，完成对三维四向编织复合材料的RVE模型中所有纤维段单元的数据计算；
(36)在步骤(35)的基础上编写计算基体单元的数据的程序，在商业有限元软件Abaqus中采用C3D8R完成，基体单元的数据计算所需要的几何和材料参数均由步骤(31)中输入的几何和材料参数得到；
(37)在步骤(35)和(36)编写的程代码的基础上，编写将计算得到的纤维和基体的单元数据的程序，实现将纤维嵌入到基体中，最后在程序中编写输出三维四向编织复合材料RVE模型的inp文件的代码，从而完成了生成纤维网格信息和基体网格信息的自动化程序的编写；
(4)根据步骤(3)中所建立的生成纤维网格信息和基体网格信息的自动化程序，将网格信息和材料信息输出到可导入Abaqus中进行仿真模拟的inp文件，从而实现了三维四向编织复合材料的参数化建模过程。
所述步骤(1)在建立的RVE模型中纤维束间和纤维束内部的基体相同，与实际情况更加符合。
所述步骤(31)中仅需要在编制的生成纤维网格信息和基体网格信息的自动化程序中输入编织复合材料所需要的几何和材料参数，就生成RVE模型的inp文件，大大降低了建模成本。
本发明与现有技术的优点在于：首先，传统RVE模型(如图1所示)建模方法由于引入了“假定”的浸渍纤维段并将基体进行了区分，导致建立的模型与实际情况存在很大的差别；其次，传统RVE模型建模方法在建立数万个纤维段的工作上存在困难。为了解决这些问题，在有限元软件Abaqus的前处理模块引入Matlab程序对RVE模型进行建立，将所得到的RVE模型的inp文件导入Abaqus软件的后处理模块中进行数值计算，从而大大降低RVE模型建立的工作量。经过与实验结果的对比表明，所建立的模型更加符合实际情况；克服了传统模型建立周期性网格的困难，可以通过简单地改变输入值，实现任意编织角度、纤维体积分数的RVE模型，从而节省人力和物力；该模型提供了纤维行为的直观见解，包括纤维强度分散性，纤维脱粘，纤维滑移，纤维弯曲和纤维束扭曲等，有一定的工程意义和实用价值。
附图说明
图1为本发明建立RVE模型inp文件的流程；
图2为利用本发明建立的编织角分别为10°(a),20°(b),30°(b)时RVE模型的纤维束的编织情况；
图3为具体实施方式的步骤(2)中的单胞内部纤维束与六面体边界纤维束交截面示意图。
具体实施方式
下面结合附图对本发明进行详细说明。
如图1所示，对本发明的一种三维四向编织复合材料的参数化建模方法具体如下：
(1)针对三维四向编织复合材料的内部纤维束结构进行立体几何分析，提取六面体结构的几何单胞模型(RVE模型)；
(2)根据立体几何映射关系，形成单胞内部纤维束与六面体边界的相交截面控制方程：
根据(1)中进行的几何分析建立了映射关系，将纤维束主截面上的点全部映射到纤维束与椭圆截面上，如图3所示是每条纤维束主截面与在六面体边界上形成的两个椭圆截面之间映射关系示意图，可以得到每条纤维束主截面与在六面体边界上形成的两个椭圆截面之间的映射关系，该映射关系即为相交截面控制方程，该方程的具体形式如下：
在RVE模型中共有14段纤维，这些纤维段根据方向不同可以分为四类：
其中α为纤维编织角,
单位法向量为的纤维段接下来会被用来介绍如何获取RVE模型边界上的斜截面。如图3,纤维段的法截面是长轴为a，短轴为b的椭圆，法截面在局部坐标系x0Oy0中的表达式为：
在全局坐标系xyz中，x0、y0轴的单位向量分别为：
其中然后，图3中的长轴平面和短轴平面的单位法向量分别为：
1)z垂直边界上的斜截面：
z垂直边界的单位法向量为：
假定x1、y1轴是z垂直边界上的x0或y0轴映射的结果，那么x1、y1轴将是长、短轴平面与z垂直边界的相交线。
方程1-11为x1轴的单位向量的表达式：
因为x1轴是长轴平面与z垂直边界的相交线，所以同时垂直于和
联立方程11～13可得：
x1轴和x0轴的向量的夹角为：
其中
随后，x1轴的长度可以求得：
同样地，y1轴的单位向量、y1轴和y0轴向量的夹角以及y1轴的长度都可以分别获得：
可以计算得到因此x1轴垂直于y1轴，即z垂直边界上的斜截面1仍然是椭圆形，可知x1轴和全局坐标的x轴之间的向量夹角为：
方程1-16、1-19和1-20总结了P和P`之间的映射关系，其中P(x0,y0)代表了局部坐标系的点，点P`(x1,y1,z1)在全局坐标系中。映射P(x0,y0)到P`(x1,y1,z1)一共有以下四步：(1)使用x0表示全局x，y0表示全局y，0表示全局z。(2)将x坐标延伸至1/cosγ倍，将y坐标延伸至1倍。(3)通过旋转坐标(4)用(x1c,y1c,z1c)转换坐标。因此，式(4)中的函数g1可以表示为：
其中(x1c,y1c,z1c)是纤维段空间分布有关的斜段1的中心坐标。
结合等式1-5和1-21得到了全局坐标系xyz中的斜截面1的表达式：
同理可得：
斜截面2的表达式：
斜截面3的表达式：
(3)基于Matlab平台调用步骤(2)形成的单胞及其纤维束的空间几何控制方程，编写在空间内生成纤维网格信息和基体网格信息的自动化程序：
1)输入编织复合材料所需要的几何和材料参数，其中包括：编织角α、纤维束主截面长轴长度a、纤维束主截面短轴长度b、纤维束的包装系数k、纤维直径df、纤维在RVE模型中的体积分数vf、纤维的单元长度lf、基体的单元长度lm、纤维的轴向弹性模量纤维的横向弹性模量纤维在1-2平面上和1-3平面上的剪切模量纤维在2-3平面上的剪切模量纤维在1-2平面上的泊松比基体的弹性模量Em和基体的泊松比νm；
2)利用“拉丁方设计”在纤维束的主截面上生成随机分布的代表纤维束中随机分布的纤维的点，这些随机点代表了纤维在纤维束中的随机分布状态。纤维的主截面为椭圆面，其长轴和短轴的长度分别为由步骤1)中输入的纤维束主截面长轴长度a、纤维束主截面短轴长度b来决定；
3)取步骤2)中的一束纤维，根据步骤(2)中建立的单胞内部纤维束与六面体边界的相交截面控制方程将纤维束主截面上点的坐标映射到RVE模型的两个椭圆截面上；
4)关联每条纤维段的两个端点(P`到P``)，并计算纤维单元(在商业有限元软件Abaqus中采用B31)的数据，其中单元的尺寸由步骤1)输入的纤维的单元长度lf所决定；
5)基于Matlab平台编写能重复进行1)-4)步骤运算的代码，使之具有能完成对三维四向编织复合材料的RVE模型中所有纤维段单元的数据计算；
6)在步骤5)的基础上编写计算基体(在商业有限元软件Abaqus中采用C3D8R)单元的数据的程序，基体单元的数据计算所需要的几何和材料参数均由步骤1)中输入的几何和材料参数得到，其中基体单元C3D8R的尺寸由7)在步骤5)和6)编写的代码的基础上，编写将计算得到的纤维和基体的单元数据进行整合的代码，在Abaqus中可以通过关键词“embedded element”实现将纤维嵌入到基体中，最后在程序中编写输出三维四向编织复合材料RVE模型的inp文件的代码，从而完成了生成纤维网格信息和基体网格信息的自动化程序的编写。
(4)根据步骤(3)中所建立的生成纤维网格信息和基体网格信息的自动化程序，将网格信息和材料信息输出到可导入Abaqus中进行仿真模拟的inp文件，从而实现了三维四向编织复合材料的参数化建模过程。
本发明未详细阐述部分属本领域公知技术。
提供以上实施例仅仅是为了描述本发明的目的，而并非要限制本发明的范围。本发明的范围由所附权利要求限定。不脱离本发明的精神和原理而做出的各种等同替换和修改，均应涵盖在本发明的范围之内。

Claims (2)

1.一种三维四向编织复合材料的参数化建模方法，其特征在于，包括以下步骤：
(1)针对三维四向编织复合材料的内部纤维束结构进行立体几何分析，提取六面体结构的几何单胞模型，即RVE模型；
(2)根据立体几何映射关系，形成单胞内部纤维束与六面体边界的相交截面控制方程，具体为根据(1)中进行的几何分析建立了映射关系，将纤维束主截面上的点全部映射到纤维束与六面体相交产生的椭圆截面上，得到每条纤维束主截面与在六面体边界上形成的两个椭圆截面之间的映射关系，所述映射关系即为相交截面控制方程；
(3)基于Matlab平台调用步骤(2)形成的单胞及其纤维束的空间几何控制方程，编写在空间内生成纤维网格信息和基体网格信息的自动化程序，具体为：
(31)输入编织复合材料所需要的几何和材料参数；
(32)利用“拉丁方设计”的方法，在纤维束的主截面上生成随机分布的代表纤维束中随机分布的纤维的点，纤维的主截面为椭圆面，椭圆面的长轴和短轴的长度分别为由输入的几何和材料参数来决定；
(33)取三维四向编织复合材料单胞的一束纤维，根据步骤(2)中建立的单胞内部纤维束与六面体边界的相交截面控制方程将纤维束主截面上点的坐标映射到RVE模型的两个椭圆截面上；
(34)关联每条纤维段的两个端点，并计算纤维单元的数据，在商业有限元软件Abaqus中采用B31完成，其中B31单元的尺寸由步骤(31)所决定；
(35)基于Matlab平台编写能重复进行(31)-(34)步骤运算的代码，完成对三维四向编织复合材料的RVE模型中所有纤维段单元的数据计算；
(36)在步骤(35)的基础上编写计算基体单元的数据的程序，在商业有限元软件Abaqus中采用C3D8R完成，基体单元的数据计算所需要的几何和材料参数均由步骤(31)中输入的几何和材料参数得到；
(37)在步骤(35)和(36)编写的程代码的基础上，编写将计算得到的纤维和基体的单元数据的程序，实现将纤维嵌入到基体中，最后在程序中编写输出三维四向编织复合材料RVE模型的inp文件的代码，从而完成了生成纤维网格信息和基体网格信息的自动化程序的编写；
(4)根据步骤(3)中所建立的生成纤维网格信息和基体网格信息的自动化程序，将网格信息和材料信息输出到可导入Abaqus中进行仿真模拟的inp文件，从而实现了三维四向编织复合材料的参数化建模过程。
2.根据权利要求1所述的一种三维四向编织复合材料的参数化建模方法，其特征在于：所述步骤(1)在建立的RVE模型中纤维束间和纤维束内部的基体相同，与实际情况更加符合。
CN201810797369.8A 2018-07-19 2018-07-19 一种三维四向编织复合材料的参数化建模方法 Active CN108984910B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810797369.8A CN108984910B (zh) 2018-07-19 2018-07-19 一种三维四向编织复合材料的参数化建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810797369.8A CN108984910B (zh) 2018-07-19 2018-07-19 一种三维四向编织复合材料的参数化建模方法

Publications (2)

Publication Number Publication Date
CN108984910A CN108984910A (zh) 2018-12-11
CN108984910B true CN108984910B (zh) 2019-05-28

Family

ID=64550320

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810797369.8A Active CN108984910B (zh) 2018-07-19 2018-07-19 一种三维四向编织复合材料的参数化建模方法

Country Status (1)

Country Link
CN (1) CN108984910B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109657412B (zh) * 2019-01-22 2022-05-03 北京航空航天大学 一种陶瓷基复合材料的疲劳寿命的多尺度预测方法
CN110008531B (zh) * 2019-03-15 2020-01-31 北京应用物理与计算数学研究所 一种含间隙原子的无序固溶材料原子结构建模方法
CN110555279A (zh) * 2019-09-10 2019-12-10 北京航空航天大学 双随机条件下三维四向编织复合材料强度的多尺度分析方法
CN110765678A (zh) * 2019-09-20 2020-02-07 哈尔滨理工大学 一种三维四向编织复合材料的rve模型离散化方法
CN111950081A (zh) * 2020-07-29 2020-11-17 哈尔滨理工大学 一种金字塔型梯度点阵结构的参数化建模方法
CN113361147B (zh) * 2021-07-21 2023-04-11 湖北大学 三维复合材料的导热模型的构建方法、系统、终端、介质
CN116976127A (zh) * 2023-08-02 2023-10-31 南京航空航天大学 一种3d机织复合材料参数化建模的方法
CN116911059B (zh) * 2023-08-02 2024-03-15 南京航空航天大学 一种3d机织复合材料参数化建模的方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5112545A (en) * 1990-02-14 1992-05-12 Airfoil Textron Inc. Composite preforms and articles and methods for their manufacture
CN103970969A (zh) * 2014-05-30 2014-08-06 哈尔滨工业大学 一种使用有限元方法模拟编织复合材料双轴试验确定材料性能参数的方法
CN105354390A (zh) * 2015-12-01 2016-02-24 北京航空航天大学 一种预测平面斜交编织复合材料双向拉伸模量与强度的新方法
CN106815408A (zh) * 2016-12-23 2017-06-09 上海交通大学 长纤维复合材料随机结构生成及其弹性性能预测方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102938019B (zh) * 2012-10-22 2016-02-03 中材科技股份有限公司 三维编织复合材料性能自适应单元分析方法
CN106202728B (zh) * 2016-07-12 2019-04-09 哈尔滨工业大学 基于Micro-CT三维编织复合材料非均匀Voxel网格离散方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5112545A (en) * 1990-02-14 1992-05-12 Airfoil Textron Inc. Composite preforms and articles and methods for their manufacture
CN103970969A (zh) * 2014-05-30 2014-08-06 哈尔滨工业大学 一种使用有限元方法模拟编织复合材料双轴试验确定材料性能参数的方法
CN105354390A (zh) * 2015-12-01 2016-02-24 北京航空航天大学 一种预测平面斜交编织复合材料双向拉伸模量与强度的新方法
CN106815408A (zh) * 2016-12-23 2017-06-09 上海交通大学 长纤维复合材料随机结构生成及其弹性性能预测方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"Algorithm Design of Four-Step Three-Dimensional Braided Composite Structures in Matlab Environment";Deng Chao et al;《Applied Mechanics and Materials》;20130816;1144-1147
"Computational analysis of fatigue behavior of 3D 4-directional braided composites based on unit cell approach";Wenfeng Hao et al;《Advances in Engineering Software》;20150110;38-52
"一种改进的三维四向编织复合材料单胞模型及宏观弹性常数预测方法";王荣桥等;《复合材料学报》;20170930;第34卷(第9期);1973-1981
"基于ABAQUS的三维四向编织复合材料参数化有限元建模";姜慧等;《甘肃科技》;20120531;第28卷(第9期);56-58

Also Published As

Publication number Publication date
CN108984910A (zh) 2018-12-11

Similar Documents

Publication Publication Date Title
CN108984910B (zh) 一种三维四向编织复合材料的参数化建模方法
Thompson et al. Handbook of grid generation
CN100418108C (zh) 三维扫描系统中的图形重构方法
CN100468418C (zh) 由边界表示数据生成体数据的方法及其程序
CN103020358B (zh) 针对机械装备的可适应性动态设计平台构建方法
CN105242538B (zh) 基于图层的零件多加工工序mbd模型及实现方法
Anderson et al. Parametric deformation of discrete geometry for aerodynamic shape design
CN107562980B (zh) 一种基于bim模型的几何信息提取方法
Li et al. Geometry synthesis on surfaces using field-guided shape grammars
CN107145686A (zh) 一种利用AutoCAD创建支架ANSYS有限元分析模型的方法
CN101826016A (zh) 支持多核并行程序设计的可视化建模及代码框架生成方法
CN108090313A (zh) 一种复杂岩石裂隙模型建模识别方法
Peng et al. A node split method for crack growth problem
CN111898295A (zh) 一种变刚度复合材料层合板的有限元建模方法
CN103389704B (zh) 基于Petri网控制流程建模的创建石膏砌块成型生产线方法
CN113361147B (zh) 三维复合材料的导热模型的构建方法、系统、终端、介质
Shchurova A methodology to design a 3D graphic editor for micro-modeling of fiber-reinforced composite parts
Wissink et al. Validation of the strand grid approach
CN105243238B (zh) 一种一体化快速产品迭代成形装置及其方法
Mayer et al. Computational Geometry Reconstruction from 3D Topology Optimization Results: A New Parametric Approach by the Medial Axis
CN107729703A (zh) 一种基于“Ruby”语言的“Sketch Up”厂房工艺设计方法
Hui et al. Research on visual 3D assembly process design and simulation for marine diesel engine
Wei et al. Concept, method and application of computational BIM
Samareh Geometry modeling and grid generation for design and optimization
CN106528913A (zh) 一种页面模块封装重用方法及装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant