公开/公告号CN114861398A
专利类型发明专利
公开/公告日2022-08-05
原文格式PDF
申请/专利权人 中国人民解放军火箭军工程大学;
申请/专利号CN202210338967.5
申请日2022-04-01
分类号G06F30/20(2020.01);G16C60/00(2019.01);G06F113/26(2020.01);G06F119/14(2020.01);
代理机构西安方诺专利代理事务所(普通合伙) 61285;
代理人景丽娜
地址 710025 陕西省西安市灞桥区洪庆街道同心路2号
入库时间 2023-06-19 16:16:00
法律状态公告日
法律状态信息
法律状态
2022-08-23
实质审查的生效 IPC(主分类):G06F30/20 专利申请号:2022103389675 申请日:20220401
实质审查的生效
2022-08-05
公开
发明专利申请公布
技术领域
本发明属于复合材料多尺度力学分析领域,具体地涉及一种预测纤维增强复合材料宏观力学性能的方法。
背景技术
目前,预测复合材料宏观工程弹性常数的方法可以分为两大类:一是解析计算方法,以经典混合法(Rule of Mixture,RoM方法)、Chamis方法、夹杂法、Mori-Tanaka(M-T)方法、自洽法和桥联法等为代表;二是有限元方法,以夹杂法、渐近均匀化方法为代表。解析计算方法计算效率高,在工程上应用快捷便利,但精度较低;有限元方法虽然计算过程繁琐,但精度较高。目前,在预测未知复合材料宏观力学性能时,通常都是采用有限元方法得到宏观工程弹性常数。
在纤维增强复合材料多尺度计算与工程分析领域经常需要分析纤维体积分数对复合材料宏观力学性能的影响。若采用解析计算方法,除过RoM方法和Chamis方法外,其余算法如M-T方法都较为复杂或者需要先获取实验数据,但RoM方法和Chamis方法计算结果误差较大。而有限元方法针对特定结构较为便捷,计算纤维体积分数对宏观力学性能的影响时十分繁琐。因此,发展一种简便的解析计算方法十分必要。
发明内容
本发明针对上述现有技术存在的问题,本发明提出了一种预测纤维增强复合材料宏观力学性能的方法(updated Rule of Mixture,uRoM方法),通过解决经典混合法和Chamis方法带来的纤维和基体间变形不协调问题,显著提高混合法预测精度,实现纤维增强复合材料宏观力学性能的预测,对纤维增强复合材料尤其是高性能碳纤维增强复合材料的多尺度计算与工程分析提供指导。相较于Chamis方法和RoM方法,uRoM方法在预测宏观弹性模量和剪切模量时与作为基准的M-T方法计算结果误差最小,精度最高,可便捷化地实现纤维增强复合材料宏观力学性能的预测。
为了实现上述目的,本发明所采用的技术方案如下:
一种预测纤维增强复合材料宏观力学性能的方法,包括:
步骤1:从纤维增强复合材料中取一代表体积元,即RVE,在2方向施加应力,2方向即坐标系y方向,计算纤维和基体沿2方向的纤维应变和基体应变;
步骤2:在步骤1得到纤维和基体沿2方向的理论纤维应变和基体应变的基础上,根据纤维基体的实际变形,在变形区域引入一变形协调平面和拉伸变形协调因子α,计算纤维和基体沿2方向加载时的真实应变,进而求解2方向的真实宏观弹性模量;
步骤3:基于步骤2求解真实宏观弹性模量的公式,引入剪切变形协调因子β和γ,求解1-2方向、1-3方向和2-3方向的纤维增强复合材料真实宏观工程弹性常数。
优选的,所述步骤1包括:
步骤1.1:从纤维增强复合材料中取一代表体积元,即RVE,一个RVE由基体包裹着纤维组成,纤维横截面为圆形,为横观各向同性材料,基体为各向同性材料,定义f表示纤维,m表示基体,用1,2,3表示笛卡尔坐标系的x,y,z方向,其中1方向为纤维纵向,则纤维的工程弹性常数定义为,1方向纤维的弹性模量E
步骤1.2:为建模简便化,在RVE中沿2-3方向截取一个真实RVE平面,定义纤维半径为r,纤维体积分数为c
步骤1.3:将步骤1.2的真实RVE平面转换为一纤维体积分数相同的等效RVE平面,计算出c
步骤1.4:对步骤1.3的等效RVE平面的2方向施加应力σ
步骤1.5:将步骤1.4的纤维应变ε
沿2方向在l上的应变差Δε
优选的,所述步骤2具体为:
步骤2.1:根据纤维基体变形实际,对于整个等效RVE,根据式(3),求得纤维和基体沿2方向额外应变的和
步骤2.2:设纤维和基体的变形一致面为等效RVE平面的变形协调面N
步骤2.3:因为在σ
步骤2.4:设纤维名义变形为Δaf,基体名义变形为Δl
步骤2.5:设纤维名义应变为
步骤2.6:在步骤2.5的基础上,若存在纤维基体间相互作用,定义纤维真实应变为
进而计算出2方向的真实宏观弹性模量
优选的,所述步骤3具体包括:
步骤3.1:根据RoM方法和Chamis方法,G
步骤3.2:基于步骤2求解真实宏观弹性模量的公式,定义1方向的真实宏观弹性模量
与现有技术相比,本发明的有益效果是:
本发明提出的uRoM方法在步骤2和步骤3中考虑了纤维和基体的变形协调,求解了等效RVE上纤维增强复合材料的真实宏观应变,解决了Chamis方法和RoM方法预测宏观弹性模量和剪切模量精度不高的问题。uRoM方法与作为基准的M-T方法相比,在预测宏观弹性模量和剪切模量时误差最小,精度最高,便于工程计算,可便捷化地实现纤维增强复合材料宏观力学性能的预测。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。
在附图中:
图1为本发明的方法流程图;
图2为本发明纤维增强复合材料RVE示意图;
图3为本发明RVE等效转换图,(a)真实RVE,(b)等效RVE;
图4为E
图5为本发明IM7/8552宏观力学性能预测结果:(a)
图6为本发明T700S/Epoxy宏观力学性能预测结果:(a)
图7为本发明T800H/Epoxy宏观力学性能预测结果:(a)
图8为本发明E-Glass/PMR-15宏观力学性能预测结果:(a)
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
实施例:
纤维增强复合材料通常由纤维和基体两部分组成,然而,由于这两部分力学性能的差异性,需要对宏观尺度上的纤维增强复合材料单向板的力学性能进行预测。
目前常用的混合预测方法模型有:
RoM方法:
RoM方法是最早被提出的预测纤维增强复合材料宏观力学性能的方法,定义1方向的真实宏观弹性模量E
式(1)中,f表示纤维,m表示基体,逗号不表示求导。
Chamis方法:
Chamis方法基于RoM方法做了改进,将真实RVE转换为一纤维截面为正方形但体积不变的等效RVE,得到其计算纤维增强复合材料宏观工程弹性常数的模型为:
对比基准为M-T方法:
M-T方法计算模型为:
式中,Z
RoM方法是最早提出的预测纤维增强复合材料宏观力学性能的混合计算方法,计算效率高,但该方法在预测宏观弹性模量和剪切模量时预测值偏低且误差太大;Chaims方法通过将真实RVE转换为等效RVE,相较于RoM方法提升了预测精度,但预测值偏高且仍有较大误差。
本发明的方法为:
一种预测纤维增强复合材料宏观力学性能的方法,如图1所示,包括:
步骤1:从纤维增强复合材料中取一代表体积元(Representative VolumeElement,RVE),在2方向施加应力,2方向即坐标系y方向,计算纤维和基体沿2方向的纤维应变和基体应变;具体为:
步骤1.1:如图2所示,从纤维增强复合材料中取一RVE,一个RVE由基体包裹着纤维组成,纤维横截面为圆形,为横观各向同性材料,基体为各向同性材料,定义f表示纤维,m表示基体,用1,2,3表示笛卡尔坐标系的x,y,z方向,其中,1方向为纤维纵向,则纤维的工程弹性常数定义为,1方向纤维的弹性模量E
各种混合方法(如Chamis方法、RoM方法)对预测宏观1方向弹性模量E
在RoM方法预测纤维复合材料宏观力学性能的过程中,计算E
σ
根据上式,由于ε
E
计算E
ε
由于σ
然而,采用上式计算得到的E
E
但采用上式计算出的E
为了直观地分析纤维基体间的相互作用,可参考Chamis方法的原理,将纤维等效为一体积不变的六面体,则由真实RVE转换得到的等效RVE如图3所示。
由图3(a)(b)可以发现,在1方向施加σ
步骤1.2:为建模简便化,在RVE中沿2-3方向截取一个真实RVE平面,如图4(a)所示,定义纤维半径为r,纤维体积分数为c
步骤1.3:将步骤1.2的真实RVE平面转换为一纤维体积相同的等效RVE平面,计算出c
步骤1.4:对步骤1.3等效RVE平面的2方向施加σ
步骤1.5:将步骤1.4的纤维应变ε
沿2方向在l上的应变差Δε
步骤2:在步骤1得到纤维和基体沿2方向的理论纤维应变和基体应变的基础上,由于真实RVE受载过程中基体和纤维所受应力和应变均不同,因此,根据纤维基体的实际变形,在变形区域引入一变形协调平面和拉伸变形协调因子α,计算纤维基体沿2方向加载时的真实应变,进而求解2方向的真实宏观弹性模量;
步骤2.1:事实上,由于纤维和基体间的相互作用(如内聚力),在非破坏性外载荷下不会出现如此大的变形差和应变差,而是应该保持变形的一致性,这表明纤维和基体界面间的正应力分别影响了纤维和基体的应力。
对于整个等效RVE,根据式(3),求得纤维和基体沿2方向额外应变的和
步骤2.2:设纤维和基体的变形一致面为等效RVE平面的变形协调面N
步骤2.3:因为在σ
步骤2.4:设纤维名义变形为Δa
步骤2.5:设纤维名义应变为
步骤2.6:在步骤2.5的基础上,由于纤维基体间存在相互作用,定义纤维真实应变为
进而计算出2方向的真实宏观弹性模量
步骤3:引入剪切变形协调因子β和γ,求解纤维增强复合材料真实宏观工程弹性常数。具体包括:
步骤3.1:根据RoM方法和Chamis方法,G
步骤3.2:参考步骤2.6中求解等效RVE沿2方向上的真实宏观应变
据此,可求得1-2方向的真实宏观剪切模量
定义1方向的真实宏观弹性模量
仿真实验:
为验证本发明uRoM方法预测纤维增强复合材料宏观力学性能的精度,使用表2中的纤维和基体工程弹性常数,计算纤维体积分数c
表2纤维和基体工程弹性常数
对比图5~图8中IM7/8552、T700S/Epoxy、T800H/Epoxy和E-Glass/PMR-15等四种材料的宏观力学性能,可以发现,在计算
由图5(b)、图6(b)、图7(b)和图8(b)可以发现,与精度较高作为基准的M-T方法预测结果相比,Chamis方法计算得到的
由上述分析可知,考虑变形协调提出的预测纤维增强复合材料宏观力学性能的uRoM方法与M-T方法误差较小,能够高精度地预测纤维增强复合材料的宏观弹性模量和剪切模量,便于工程应用计算。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
机译: 一种减少纤维素纤维中杂质的方法,用于制造纤维增强水泥复合材料;纸浆的提取过程;低鳕鱼纤维素纤维;一种包含增强纤维配方的建筑用复合材料;
机译: 涂料剂,片状中间基材,光固化树脂膜,纤维增强复合材料中间体,纤维增强复合材料,生产纤维增强复合材料中间体的方法,以及生产纤维增强复合材料的方法
机译: 纤维增强复合材料坯料,纤维增强复合材料组分,转子叶片元件,转子叶片和风轮机以及生产纤维增强复合材料坯料的方法和生产纤维增强复合材料组分的方法