CN109190167A - 单向长纤维增强复合材料的横向微观结构生成方法 - Google Patents
单向长纤维增强复合材料的横向微观结构生成方法 Download PDFInfo
- Publication number
- CN109190167A CN109190167A CN201810862061.7A CN201810862061A CN109190167A CN 109190167 A CN109190167 A CN 109190167A CN 201810862061 A CN201810862061 A CN 201810862061A CN 109190167 A CN109190167 A CN 109190167A
- Authority
- CN
- China
- Prior art keywords
- fiber
- random
- hole
- rve
- composite material
- 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.)
- Granted
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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
本发明提供了一种单向长纤维增强复合材料的横向微观结构生成方法,在需生成RVE模型的目标区域内,确定RVE模型的初始参数,将规则分布的纤维位置作为初始纤维位置,结合相邻的行或列之间交叉算法,在保证纤维在边界处的周期性条件下,通过随机摄动方法,生成具有周期重复的纤维随机分布的RVE,基于获得的随机纤维位置坐标,确定微孔隙的初始位置,并随机确定孔隙的尺寸和形状,通过对孔隙的位置进行随机摄动,最终建立考虑纤维和微孔隙随机分布的复合材料横向微观模型。本发明考虑了单向长纤维复合材料横向微观结构的重构技术,针对纤维和孔隙的随机分布,均采用随机摄动方法,可以有效、高效地建立考虑纤维和孔隙随机分布的横向RVE模型。
Description
技术领域
本发明涉及一种复合材料,具体涉及一种复合材料的横向微观结构生成方法。
背景技术
复合材料层合板由于具有高比强、高比刚、设计性强等优点,现在被广泛应用于航天、汽车、土木等行业。由于单向复合材料板是层合板的组成部分,因此单向长纤维复合材料结构在外载作用下的变形和损伤评估成为复合材料领域的研究重点。
复合材料的横向力学性能是材料整体性能中薄弱的方向,使用数值分析方法对复合材料的横向性能进行预测,可以减少试验的次数,节约成本,用于分析不同因素,例如材料内部结构、组分材料性能、材料制作过程等,并且可以基于数值方法,对材料进行设计和优化,能够真实描述材料复合材料横向微观结构(代表性体积单元RVE)是进行数值分析预测的第一步。但是材料在制作过程中,会导致纤维的随机分布和孔隙的存在,现有的横向微观结构大多仅考虑纤维的随机分布特征,没有考虑孔隙的存在。
发明内容
发明目的:本发明的目的在于针对现有技术的不足,提供一种单向长纤维增强复合材料的横向微观结构生成方法,可建立同时考虑纤维和孔隙随机分布的真实复合材料横向微观结构代表性体积单元。
技术方案:本发明提供了一种单向长纤维增强复合材料的横向微观结构生成方法,包括以下步骤:
(1)在需生成RVE模型的目标区域(即RVE窗口)内,确定RVE模型的初始参数,将规则分布的纤维位置作为初始纤维位置;
(2)结合相邻的行或列之间交叉算法,在保证纤维在边界处的周期性条件下,通过随机摄动方法,生成具有周期重复的纤维随机分布的RVE;
(3)基于获得的随机纤维位置坐标,确定微孔隙的初始位置,并随机确定孔隙的尺寸和形状;
(4)通过对孔隙的位置进行随机摄动,最终建立考虑纤维和微孔隙随机分布的复合材料横向微观模型。
进一步,步骤(1)所述初始参数包括RVE窗口边长s、纤维体积含量Vf、纤维半径r、纤维个数n和孔隙的体积含量Vv。
进一步,步骤(2)包括以下步骤:
(21)行或列之间交叉算法是指,对相邻行或列之间的纤维在x或y方向上有相对整体的移动,移动增量为Δx或Δy,首先针对相邻的行和列进行整体的相对移动;
(22)在保证纤维不发生重叠的情况下,首先确定一根纤维在平面内x和y方向上可移动范围;
(23)在边界上对应纤维之间横向或纵向间距相同的条件下,即对RVE窗口进行复制并形成九宫格排列的窗口后,内部的所有纤维均是完整的且相邻两个小RVE窗口之间边界处的纤维是互补的,随机选择位置移动变化量,重新确定纤维的位置;
(24)对RVE窗口内所有的纤维依次进行步骤(22)和步骤(23);
(25)重复步骤(21)至(24),直至最大循环数;
(26)输出每个纤维的位置信息。
进一步,步骤(3)通过使用生成的纤维随机分布代表性体积单元的纤维位置,在保证四根纤维圆心组成的四边形内部不包含其他纤维的条件下,确定所有满足条件的四根纤维几何中心位置,将其作为微孔隙的初始位置;复合材料内部微孔隙的形状简化为圆形或椭圆形,因此在保证孔隙体积含量的情况下,随机选取孔隙的形状及其对应孔隙形状和尺寸。
进一步,步骤(4)针对尺寸非零的孔隙,在其初始位置周围的区域内进行位置的随机摄动,从而确定随机分布孔隙的位置信息,将具有随机尺寸特征的孔隙位置信息和纤维随机分布的位置信息结合,即可以建立纤维和孔隙均随机分布单胞模型。
有益效果:本发明考虑了单向长纤维复合材料横向微观结构的重构技术,针对纤维和孔隙的随机分布,均采用随机摄动方法,可以有效、高效地建立考虑纤维和孔隙随机分布的横向RVE模型。
附图说明
图1为本发明方法的流程框图;
图2为纤维位置初始规则分布横向RVE模型示意图;
图3为边界纤维距离约束和纤维随机移动取值范围示意图;
图4为纤维行和列的初始交叉运动后纤维分布示意图;
图5为孔隙几何特征示意图,其中(a)为圆形孔隙,(b)为椭圆形孔隙;
图6为孔隙初始位置和随机摄动流程示意图;
图7为单向长纤维增强复合材料横向微观结构RVE示意图。
具体实施方式
下面对本发明技术方案进行详细说明,但是本发明的保护范围不局限于所述实施例。
一种单向长纤维增强复合材料的横向微观结构生成方法,如图1所示,具体操作如下:
步骤(1)在RVE窗口内确定纤维的初始位置:
1.1)确定RVE模型的初始参数:模型初始参数包括RVE窗口(坐标轴x-y平面上的正方形或者矩形,这里选择正方形)、边长s、纤维体积含量Vf、纤维半径r、完整的纤维个数n和孔隙的体积含量Vv,其中s=37.6μm,Vf=50%,r=5μm,n=36, Vv=1%;
1.2)在RVE窗口内生成纤维周期性分布的坐标位置:选择使用纤维正方形分布的单胞模型,如图2所示,或者同样可以选择纤维菱形或者六边形分布的单胞模型。
步骤(2)结合相邻的行或列之间交叉算法,在保证纤维在边界处的周期性条件下,通过随机摄动方法,生成具有周期重复的纤维随机分布的RVE:
2.1)针对相邻的行和列进行整体的相对移动ΔX(ΔY);
2.2)确定纤维的随机摄动范围,即在不发生纤维重叠的条件下,确定纤维i圆心坐标xi和yi能改变的范围Δxi和Δyi,如图3所示;
2.3)在保证纤维之间不发生交叉和边界上对应纤维之间横向或纵向间距相同的条件下(图4),随机确定纤维新的新坐标;
2.4)对RVE窗口内所有的纤维依次进行步骤2.2)和步骤2.3);
2.5)重复步骤2.1)至2.4),直至最大循环数N=500;
2.6)输出每个纤维的位置信息。
步骤(3)基于获得的随机纤维位置坐标,确定微孔隙的初始位置,并随机确定孔隙的尺寸和形状:
3.1)在保证孔隙体积含量下,随机选取孔隙的形状及其对应孔隙尺寸,圆形和椭圆形孔隙示意图如图5所示,需要随机确定的孔隙的尺寸参数包括圆形的半径r或者椭圆形的长轴a、短轴尺寸b以及长轴与水平方向的夹角θ;
3.2)使用相邻四个纤维现在的位置坐标,确定几何中心,将其作为微孔隙j的初始位置,选择周围的正方形区域作为孔隙的摄动范围,区域边长为3d,d为6倍的纤维半径(r),如图6所示。
步骤(4)针对尺寸非零的孔隙,在其初始位置周围的区域内进行位置的随机摄动,如图6所示,进而建立纤维和孔隙均随机分布单胞模型,如图7所示。
Claims (5)
1.一种单向长纤维增强复合材料的横向微观结构生成方法,其特征在于:包括以下步骤:
(1)在需生成RVE模型的目标区域内,确定RVE模型的初始参数,将规则分布的纤维位置作为初始纤维位置;
(2)结合相邻的行或列之间交叉算法,在保证纤维在边界处的周期性条件下,通过随机摄动方法,生成具有周期重复的纤维随机分布的RVE;
(3)基于获得的随机纤维位置坐标,确定微孔隙的初始位置,并随机确定孔隙的尺寸和形状;
(4)通过对孔隙的位置进行随机摄动,最终建立考虑纤维和微孔隙随机分布的复合材料横向微观模型。
2.根据权利要求1所述的单向长纤维增强复合材料的横向微观结构生成方法,其特征在于:步骤(1)所述初始参数包括RVE窗口边长s、纤维体积含量Vf、纤维半径r、纤维个数n和孔隙的体积含量Vv。
3.根据权利要求1所述的单向长纤维增强复合材料的横向微观结构生成方法,其特征在于:步骤(2)包括以下步骤:
(21)行或列之间交叉算法是指,对相邻行或列之间的纤维在x或y方向上有相对整体的移动,移动增量为Δx或Δy,首先针对相邻的行和列进行整体的相对移动;
(22)在保证纤维不发生重叠的情况下,首先确定一根纤维在平面内x和y方向上可移动范围;
(23)在边界上对应纤维之间横向或纵向间距相同的条件下,即对RVE窗口进行复制并形成九宫格排列的窗口后,内部的所有纤维均是完整的且相邻两个小RVE窗口之间边界处的纤维是互补的,随机选择位置移动变化量,重新确定纤维的位置;
(24)对RVE窗口内所有的纤维依次进行步骤(22)和步骤(23);
(25)重复步骤(21)至(24),直至最大循环数;
(26)输出每个纤维的位置信息。
4.根据权利要求1所述的单向长纤维增强复合材料的横向微观结构生成方法,其特征在于:步骤(3)通过使用生成的纤维随机分布代表性体积单元的纤维位置,在保证四根纤维圆心组成的四边形内部不包含其他纤维的条件下,确定所有满足条件的四根纤维几何中心位置,将其作为微孔隙的初始位置;复合材料内部微孔隙的形状简化为圆形或椭圆形,因此在保证孔隙体积含量的情况下,随机选取孔隙的形状及其对应孔隙形状和尺寸。
5.根据权利要求1所述的单向长纤维增强复合材料的横向微观结构生成方法,其特征在于:步骤(4)针对尺寸非零的孔隙,在其初始位置周围的区域内进行位置的随机摄动,从而确定随机分布孔隙的位置信息,将具有随机尺寸特征的孔隙位置信息和纤维随机分布的位置信息结合,即可以建立纤维和孔隙均随机分布单胞模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810862061.7A CN109190167B (zh) | 2018-08-01 | 2018-08-01 | 单向长纤维增强复合材料的横向微观结构生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810862061.7A CN109190167B (zh) | 2018-08-01 | 2018-08-01 | 单向长纤维增强复合材料的横向微观结构生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109190167A true CN109190167A (zh) | 2019-01-11 |
CN109190167B CN109190167B (zh) | 2023-04-07 |
Family
ID=64937678
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810862061.7A Active CN109190167B (zh) | 2018-08-01 | 2018-08-01 | 单向长纤维增强复合材料的横向微观结构生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109190167B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112733408A (zh) * | 2021-02-23 | 2021-04-30 | 江西省科学院应用物理研究所 | 纤维微动和硬性填充结合生成高体积分数二维数值模型的方法 |
CN112800587A (zh) * | 2021-01-07 | 2021-05-14 | 南京航空航天大学 | 单向长纤维增强复合材料中孔隙随机分布生成方法 |
CN113221416A (zh) * | 2021-05-14 | 2021-08-06 | 上海工程技术大学 | 一种构建颗粒增强复合材料二维微观构型的方法 |
CN113792420A (zh) * | 2021-08-31 | 2021-12-14 | 西北工业大学 | 一种考虑界面的单向随机纤维复合材料微观结构生成方法 |
CN115374599A (zh) * | 2022-07-13 | 2022-11-22 | 河海大学 | 一种基于随机分布的3d纤维状孔隙冰数值建模方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106815408A (zh) * | 2016-12-23 | 2017-06-09 | 上海交通大学 | 长纤维复合材料随机结构生成及其弹性性能预测方法 |
CN107256298A (zh) * | 2017-06-01 | 2017-10-17 | 上海交通大学 | 单向纤维复合材料随机结构生成方法及其系统 |
-
2018
- 2018-08-01 CN CN201810862061.7A patent/CN109190167B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106815408A (zh) * | 2016-12-23 | 2017-06-09 | 上海交通大学 | 长纤维复合材料随机结构生成及其弹性性能预测方法 |
CN107256298A (zh) * | 2017-06-01 | 2017-10-17 | 上海交通大学 | 单向纤维复合材料随机结构生成方法及其系统 |
Non-Patent Citations (1)
Title |
---|
马健等: "基于RVE模型的空间天线结构热稳定性", 《空间科学学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112800587A (zh) * | 2021-01-07 | 2021-05-14 | 南京航空航天大学 | 单向长纤维增强复合材料中孔隙随机分布生成方法 |
CN112733408A (zh) * | 2021-02-23 | 2021-04-30 | 江西省科学院应用物理研究所 | 纤维微动和硬性填充结合生成高体积分数二维数值模型的方法 |
CN112733408B (zh) * | 2021-02-23 | 2023-03-31 | 江西省科学院应用物理研究所 | 纤维微动和硬性填充结合生成高体积分数二维数值模型的方法 |
CN113221416A (zh) * | 2021-05-14 | 2021-08-06 | 上海工程技术大学 | 一种构建颗粒增强复合材料二维微观构型的方法 |
CN113792420A (zh) * | 2021-08-31 | 2021-12-14 | 西北工业大学 | 一种考虑界面的单向随机纤维复合材料微观结构生成方法 |
CN113792420B (zh) * | 2021-08-31 | 2024-03-19 | 西北工业大学 | 一种考虑界面的单向随机纤维复合材料微观结构生成方法 |
CN115374599A (zh) * | 2022-07-13 | 2022-11-22 | 河海大学 | 一种基于随机分布的3d纤维状孔隙冰数值建模方法 |
CN115374599B (zh) * | 2022-07-13 | 2023-06-16 | 河海大学 | 一种基于随机分布的3d纤维状孔隙冰数值建模方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109190167B (zh) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109190167A (zh) | 单向长纤维增强复合材料的横向微观结构生成方法 | |
Lu et al. | Feature-based tool path generation approach for incremental sheet forming process | |
US10303157B2 (en) | Additive manufacturing method for three-dimensional object | |
CN106077638B (zh) | 一种用于增材制造的蜂窝式分区扫描方法 | |
CN103611646B (zh) | 喷涂机器人空间路径规划方法 | |
CN103722171B (zh) | 一种用于选择性激光烧结的蜂窝式激光扫描方法 | |
CN105528503B (zh) | 一种基于结构分解的大型构件动态优化设计方法 | |
CN102335687A (zh) | 渐进成形工件的方法 | |
CN107330148B (zh) | 三维编织复合材料六面体有限元模型自动生成方法 | |
CN112157911B (zh) | 一种自支撑的3d打印梯度弹性多孔材料微结构设计方法 | |
CN106925776A (zh) | 一种控制增材制造应力变形的分区域扫描路径生成方法 | |
CN104772463A (zh) | 一种用于金属增材制造的支撑结构及制作方法 | |
CN106599485B (zh) | 一种框梁整体件优化设计方法 | |
CN104573191A (zh) | 一种战时损伤武器装备零件快速再制造造型方法 | |
Kandemir et al. | Topology optimization of 2.5 D parts using the SIMP method with a variable thickness approach | |
CN105772721A (zh) | 基于增材制造的光斑大小可变的扫描路径生成方法 | |
CN112001060B (zh) | 一种三维五向编织复合材料内胞模型的建模方法 | |
CN204680754U (zh) | 一种太赫兹宽带随机表面 | |
CN104360429A (zh) | 导光板及其加工方法 | |
CN102663153B (zh) | 一种异型蜂窝结构的有限元建模方法 | |
CN106400295A (zh) | 三维编织管及制作该三维编织管的三维编织机和编织工艺 | |
CN112942837B (zh) | 一种悬挑结构混凝土3d打印方法及系统 | |
CN104716425A (zh) | 一种天线及fss系统的构建方法 | |
CN113792420A (zh) | 一种考虑界面的单向随机纤维复合材料微观结构生成方法 | |
CN108363828A (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 |