CN114417634A - A multi-scale prediction method for damp-heat aging properties of plant fiber/polylactic acid composites based on mesoscopic modeling - Google Patents
A multi-scale prediction method for damp-heat aging properties of plant fiber/polylactic acid composites based on mesoscopic modeling Download PDFInfo
- Publication number
- CN114417634A CN114417634A CN202210102395.0A CN202210102395A CN114417634A CN 114417634 A CN114417634 A CN 114417634A CN 202210102395 A CN202210102395 A CN 202210102395A CN 114417634 A CN114417634 A CN 114417634A
- Authority
- CN
- China
- Prior art keywords
- composite material
- polylactic acid
- aging
- plant fiber
- fiber
- 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
- 239000000835 fiber Substances 0.000 title claims abstract description 100
- 239000002131 composite material Substances 0.000 title claims abstract description 99
- 229920000747 poly(lactic acid) Polymers 0.000 title claims abstract description 85
- 239000004626 polylactic acid Substances 0.000 title claims abstract description 85
- 238000000034 method Methods 0.000 title claims abstract description 23
- 230000003679 aging effect Effects 0.000 title description 3
- 230000032683 aging Effects 0.000 claims abstract description 76
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 45
- 238000010521 absorption reaction Methods 0.000 claims abstract description 42
- 239000000463 material Substances 0.000 claims abstract description 28
- 230000015556 catabolic process Effects 0.000 claims abstract description 23
- 238000006731 degradation reaction Methods 0.000 claims abstract description 23
- 238000012360 testing method Methods 0.000 claims abstract description 11
- 230000007613 environmental effect Effects 0.000 claims abstract description 9
- 230000008859 change Effects 0.000 claims abstract description 6
- 241000196324 Embryophyta Species 0.000 claims description 65
- 239000011159 matrix material Substances 0.000 claims description 53
- 240000000491 Corchorus aestuans Species 0.000 claims description 15
- 235000011777 Corchorus aestuans Nutrition 0.000 claims description 15
- 235000010862 Corchorus capsularis Nutrition 0.000 claims description 15
- 238000004088 simulation Methods 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 8
- 230000035882 stress Effects 0.000 claims description 7
- 238000009864 tensile test Methods 0.000 claims description 7
- 239000002657 fibrous material Substances 0.000 claims description 6
- 239000002245 particle Substances 0.000 claims description 6
- 238000013001 point bending Methods 0.000 claims description 5
- 238000012800 visualization Methods 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000001746 injection moulding Methods 0.000 claims description 3
- 230000000737 periodic effect Effects 0.000 claims description 3
- 238000007586 pull-out test Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000003325 tomography Methods 0.000 claims description 3
- 238000004140 cleaning Methods 0.000 claims description 2
- 230000001419 dependent effect Effects 0.000 claims description 2
- 239000011174 green composite Substances 0.000 abstract description 5
- 230000001808 coupling effect Effects 0.000 abstract 1
- 238000011160 research Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 229910052799 carbon Inorganic materials 0.000 description 4
- 239000003733 fiber-reinforced composite Substances 0.000 description 4
- 229920006395 saturated elastomer Polymers 0.000 description 4
- DGOBMKYRQHEFGQ-UHFFFAOYSA-L acid green 5 Chemical compound [Na+].[Na+].C=1C=C(C(=C2C=CC(C=C2)=[N+](CC)CC=2C=C(C=CC=2)S([O-])(=O)=O)C=2C=CC(=CC=2)S([O-])(=O)=O)C=CC=1N(CC)CC1=CC=CC(S([O-])(=O)=O)=C1 DGOBMKYRQHEFGQ-UHFFFAOYSA-L 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000009792 diffusion process Methods 0.000 description 3
- 229920003023 plastic Polymers 0.000 description 3
- 239000004033 plastic Substances 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 229920002488 Hemicellulose Polymers 0.000 description 2
- JVTAAEKCZFNVCJ-UHFFFAOYSA-N lactic acid Chemical compound CC(O)C(O)=O JVTAAEKCZFNVCJ-UHFFFAOYSA-N 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- KXGFMDJXCMQABM-UHFFFAOYSA-N 2-methoxy-6-methylphenol Chemical compound [CH]OC1=CC=CC([CH])=C1O KXGFMDJXCMQABM-UHFFFAOYSA-N 0.000 description 1
- 229920000049 Carbon (fiber) Polymers 0.000 description 1
- 239000006087 Silane Coupling Agent Substances 0.000 description 1
- 239000003513 alkali Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000004917 carbon fiber Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000003822 epoxy resin Substances 0.000 description 1
- 239000005007 epoxy-phenolic resin Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000013277 forecasting method Methods 0.000 description 1
- 239000003365 glass fiber Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 235000014655 lactic acid Nutrition 0.000 description 1
- 239000004310 lactic acid Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 1
- 229920001568 phenolic resin Polymers 0.000 description 1
- 229920000647 polyepoxide Polymers 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 239000011347 resin Substances 0.000 description 1
- 229920005989 resin Polymers 0.000 description 1
- 230000007281 self degradation Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 229920002994 synthetic fiber Polymers 0.000 description 1
- 239000012209 synthetic fiber Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000001291 vacuum drying Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Images
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N17/00—Investigating resistance of materials to the weather, to corrosion, or to light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N5/00—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid
- G01N5/02—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid by absorbing or adsorbing components of a material and determining change of weight of the adsorbent, e.g. determining moisture content
- G01N5/025—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid by absorbing or adsorbing components of a material and determining change of weight of the adsorbent, e.g. determining moisture content for determining moisture content
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0001—Type of application of the stress
- G01N2203/0003—Steady
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0014—Type of force applied
- G01N2203/0016—Tensile or compressive
- G01N2203/0017—Tensile
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/026—Specifications of the specimen
- G01N2203/0262—Shape of the specimen
- G01N2203/0278—Thin specimens
- G01N2203/028—One dimensional, e.g. filaments, wires, ropes or cables
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/26—Composites
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Pathology (AREA)
- Immunology (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Environmental Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Ecology (AREA)
- Biodiversity & Conservation Biology (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Reinforced Plastic Materials (AREA)
Abstract
Description
技术领域technical field
本发明属于复合材料的检测技术领域,具体涉及一种基于细观建模的植物纤维/聚乳酸复合材料湿热老化性能多尺度预测方法。The invention belongs to the technical field of detection of composite materials, and in particular relates to a multi-scale prediction method for the damp-heat aging performance of a plant fiber/polylactic acid composite material based on mesoscopic modeling.
背景技术Background technique
可降解是植物纤维增强复合材料最重要的优点之一。然而由于植物纤维的吸湿性和聚乳酸的可降解特性,其在湿度、湿热等服役环境下的耐久能力仍面临巨大挑战。目前面向汽车服役工况的复合材料老化研究主要集中在传统合成纤维增强树脂,例如碳纤维、玻璃纤维增强环氧树脂、酚醛树脂等。而对于植物纤维增强聚乳酸复合材料,由于植物纤维特有的亲水性和典型的多尺度多层次微观结构,以及聚乳酸自身的可降解特征,导致了其环境老化研究更加复杂和困难。Degradability is one of the most important advantages of plant fiber reinforced composites. However, due to the hygroscopicity of plant fibers and the degradable properties of polylactic acid, its durability in service environments such as humidity and damp heat still faces great challenges. At present, the aging research of composite materials for automobile service conditions mainly focuses on traditional synthetic fiber reinforced resins, such as carbon fiber, glass fiber reinforced epoxy resin, and phenolic resin. For plant fiber-reinforced polylactic acid composites, due to the unique hydrophilicity and typical multi-scale and multi-level microstructure of plant fibers, as well as the degradable characteristics of polylactic acid itself, the research on environmental aging is more complicated and difficult.
植物纤维增强复合材料的力学性能对湿热环境非常敏感。实验研究表明,吸湿和水热老化都会显著降低植物纤维增强复合材料的力学性能,影响其使用寿命。在植物纤维中,负责吸收水分的成分主要是半纤维素,半纤维素含量越高,水分吸收和降解程度越高,而不同的植物纤维结构形态也对水分扩散存在影响。通常在室温下的植物纤维增强复合材料的吸水规律遵循Fickian扩散规律,初始呈线性变化,长时间后吸水性逐渐变缓,趋于饱和。而在较高的温度下,吸湿行为显著加快,水分饱和时间大大缩短。The mechanical properties of plant fiber-reinforced composites are very sensitive to hot and humid environments. Experimental studies have shown that both moisture absorption and hydrothermal aging can significantly reduce the mechanical properties of plant fiber-reinforced composites and affect their service life. In plant fibers, the main component responsible for absorbing water is hemicellulose. The higher the hemicellulose content, the higher the degree of water absorption and degradation. Different plant fiber structures also have an impact on water diffusion. Usually, the water absorption law of plant fiber reinforced composites at room temperature follows the Fickian diffusion law, and the initial change is linear, and the water absorption gradually slows down and tends to be saturated after a long time. At higher temperatures, the hygroscopic behavior is significantly accelerated, and the water saturation time is greatly shortened.
由于植物纤维的亲水性和复杂结构特点,以及聚乳酸基体的可降解性,长期湿热环境作用往往会引起复合材料中各组分(纤维、基体和纤维-基体界面)性能均发生改变,而目前关于老化预测模型及方法的研究大都只针对纤维和基体,忽略了纤维/基体自身降解及其所带来的界面分离现象,不能从微-细-宏观层面多尺度充分反应分析短植物纤维增强复合材料中纤维、基体及界面各组分衰退对材料力学性能的影响。Due to the hydrophilicity and complex structural characteristics of plant fibers, as well as the degradability of the PLA matrix, the long-term humid and heat environment often causes changes in the properties of each component (fiber, matrix, and fiber-matrix interface) in the composite material. At present, most of the research on aging prediction models and methods only focus on fibers and matrix, ignoring the self-degradation of fibers/matrix and the interface separation phenomenon caused by them, and cannot fully respond and analyze short plant fiber reinforcement from the micro-fine-macro level. Effects of fiber, matrix and interface component recession on the mechanical properties of composite materials.
发明内容SUMMARY OF THE INVENTION
针对低碳汽车轻量化零部件设计需求,本发明提供一种基于细观建模的植物纤维/聚乳酸复合材料湿热老化性能多尺度预测方法,为植物纤维/聚乳酸绿色复合材料在实际工程中的设计应用中提供参考和指导。Aiming at the design requirements of lightweight parts for low-carbon automobiles, the present invention provides a multi-scale prediction method for the damp-heat aging performance of plant fiber/polylactic acid composite materials based on mesoscopic modeling, which is a useful tool for plant fiber/polylactic acid green composite materials in practical engineering. Design applications for reference and guidance.
为实现上述发明目的,本发明提供以下技术方案:一种基于细观建模的植物纤维/聚乳酸复合材料湿热老化性能多尺度预测方法,包括以下步骤:In order to achieve the above-mentioned purpose of the invention, the present invention provides the following technical solutions: a multi-scale prediction method for the wet-heat aging performance of plant fibers/polylactic acid composite materials based on mesoscopic modeling, comprising the following steps:
1)制备植物纤维/聚乳酸复合材料,分别对复合材料、聚乳酸基体和植物纤维进行老化试验:选取典型植物纤维作为填充相,可降解聚乳酸作为基体材料,制备植物纤维/聚乳酸复合材料;参考汽车行业零部件加速老化标准,把事先干燥好的植物纤维、聚乳酸基体、植物纤维\聚乳酸复合材料放入不同温度的恒温水槽中进行人工加速老化,得到不同老化温度和老化时间的老化材料;1) Prepare the plant fiber/PLA composite material, and conduct aging tests on the composite material, the PLA matrix and the plant fiber respectively: select typical plant fiber as the filling phase, and degradable PLA as the matrix material to prepare the plant fiber/PLA composite material ; Referring to the accelerated aging standard for parts in the automotive industry, put the pre-dried plant fiber, polylactic acid matrix, plant fiber\polylactic acid composite material into a constant temperature water tank with different temperatures for artificial accelerated aging, and obtain different aging temperatures and aging times. aged materials;
2)建立不同温度下的各材料吸水率随老化时间变化规律函数:对步骤1)中得到的老化材料进行吸水率测试,聚乳酸基体和植物纤维\聚乳酸复合材料的吸水率可以直接通过电子秤测量,植物纤维的吸水率由聚乳酸基体和复合材料的吸水率间接计算,对各温度下材料吸水率随老化时间变化规律进行拟合,得到三种材料的吸水率(M)与老化温度(T)和老化时间(t)之间的函数,分别为纤维Mf(t,T),聚乳酸基体Mm(t,T)和复合材料Mi(t,T);2) Establish the regular function of the water absorption rate of each material at different temperatures with aging time: The water absorption rate of the aging material obtained in step 1) is tested. The water absorption rate of the plant fiber is measured indirectly by the water absorption rate of the polylactic acid matrix and the composite material, and the variation law of the material water absorption rate with the aging time at each temperature is fitted to obtain the water absorption rate (M) of the three materials and the aging temperature. (T) and aging time (t) as a function of fiber M f (t, T), PLA matrix M m (t, T) and composite M i (t, T), respectively;
3)建立不同温度下的各组分强度随老化时间变化规律函数:对步骤1)中得到的老化材料进行强度测试:其中,对聚乳酸基体进行哑铃拉伸试验,对植物纤维进行单纤维拉伸试验,对复合材料进行单纤维拔出试验,分别得到不同温度下聚乳酸基体、植物纤维、复合材料三者强度随老化时间变化的函数:Sm(t,T),Sf(t,T)和Si(t,T);3) Establish the regular function of the strength of each component at different temperatures with aging time: The strength test is performed on the aging material obtained in step 1): wherein, the polylactic acid matrix is subjected to a dumbbell tensile test, and the plant fiber is subjected to single fiber pulling The tensile test was carried out, and the single fiber pull-out test was carried out on the composite material, and the function of the strength of the polylactic acid matrix, plant fiber and composite material with aging time at different temperatures was obtained: S m (t, T), S f (t, T) and Si (t, T);
4)分别建立聚乳酸基体、植物纤维、复合材料各组分强度与吸水率、温度之间关系函数:对前面获得的数据进行整合,通过拟合数据分别得到植物纤维、聚乳酸基体和纤维-基体界面三者强度随吸水率和老化温度变化的函数,分别为纤维Sf(M,T),基体Sm(M,T)和纤维-基体界面Si(M,T);4) Establish the relationship function between the strength of each component of polylactic acid matrix, plant fiber, and composite material, water absorption rate, and temperature: Integrate the previously obtained data, and obtain plant fiber, polylactic acid matrix and fiber- The three strengths of the matrix interface as a function of water absorption and aging temperature are respectively fiber S f (M, T), matrix S m (M, T) and fiber-matrix interface Si (M, T);
5)复合材料细观RVE模型建立:从未老化植物纤维/聚乳酸复合材料中采样,进行X射线断层扫描,获得包含材料微观结构信息的断层灰度图像;对断层图像进行三维可视化分析,建立三维视图并完成纤维及聚乳酸基体的几何简化清理,将重构几何模型导入有限元分析软件,建立植物纤维/聚乳酸复合材料的细观RVE模型;5) Establishment of the microscopic RVE model of the composite material: sampling from the unaged plant fiber/polylactic acid composite material, performing X-ray tomography to obtain a tomographic grayscale image containing the microstructure information of the material; 3D visualization analysis of the tomographic image to establish 3D view and geometric simplification and cleaning of fibers and PLA matrix are completed, the reconstructed geometric model is imported into finite element analysis software, and the meso RVE model of plant fiber/PLA composite material is established;
6)环境退化因子定义与引入:编写植物纤维、聚乳酸基体、复合材料界面模拟单元本构关系,对复合材料细观RVE模型进行属性定义,构建基于强度变化函数的各组分环境退化因子(D)函数,6) Definition and introduction of environmental degradation factors: Compose the constitutive relationship of plant fibers, polylactic acid matrix, and composite interface simulation elements, define the properties of the meso-RVE model of composite materials, and construct the environmental degradation factors of each component based on the strength change function ( D) function,
其中S(M′,T′)是温度为T′,吸湿度为M′时的强度,S(M0,T0)是材料初始强度。由上式可以得到纤维的退化因子函数Df(M,T),基体的退化因子函数Dm(M,T)和纤维-基体界面的退化因子函数Di(M,T);where S(M', T') is the strength when the temperature is T' and the moisture absorption is M', and S(M 0 , T 0 ) is the initial strength of the material. From the above formula, the degradation factor function D f (M, T) of the fiber, the degradation factor function D m (M, T) of the matrix and the degradation factor function D i (M, T) of the fiber-matrix interface can be obtained;
7)复合材料弹性性能计算:引入退化因子,对各组分弹性性能参数进行修正;在周期性边界条件下,对RVE模型施加线性不相关位移载荷并进行仿真,获得不同温度和吸水率的宏观复合材料弹性性能;7) Calculation of elastic properties of composite materials: Introduce a degradation factor to correct the elastic performance parameters of each component; under periodic boundary conditions, apply a linear uncorrelated displacement load to the RVE model and simulate it to obtain macroscopic images of different temperatures and water absorption rates. composite elastic properties;
8)复合材料失效强度计算:针对细观模型中的植物纤维、聚乳酸基体和复合材料界面三个组分,分别定义对应的初始失效准则和损伤扩展准则,引入退化因子,对失效参数进行修正,并进行数值模拟,得到不同吸水率和温度下的宏观复合材料失效强度;8) Calculation of failure strength of composite materials: For the three components of plant fiber, polylactic acid matrix and composite interface in the mesoscopic model, the corresponding initial failure criterion and damage expansion criterion are respectively defined, and the degradation factor is introduced to correct the failure parameters. , and carry out numerical simulation to obtain the failure strength of macroscopic composites under different water absorption rates and temperatures;
9)宏观复合材料湿热老化性能预测:建立复合材料拉伸、三点弯曲试件宏观模型,将基于细观RVE模型仿真得到的依赖吸水率和温度变化的复合材料宏观弹性性能函数和失效强度函数,代入到宏观模型中,设置湿热环境进行有限元分析,对湿热老化工况下不同老化程度的植物纤维/聚乳酸复合材料力学性能进行预测。9) Prediction of the wet-heat aging performance of macroscopic composite materials: establish a macroscopic model of composite tensile and three-point bending specimens, and simulate the macroscopic elastic performance function and failure strength function of composite materials that depend on water absorption and temperature changes based on the microscopic RVE model simulation. , substituted into the macro model, set up the damp-heat environment for finite element analysis, and predicted the mechanical properties of plant fiber/polylactic acid composites with different aging degrees under damp-heat aging conditions.
步骤1)中,所述植物纤维/聚乳酸复合材料由黄麻纤维和聚乳酸颗粒通过注塑成型制备得到。In step 1), the plant fiber/polylactic acid composite material is prepared from jute fiber and polylactic acid particles by injection molding.
步骤5)中,所述有限元软件为Abaqus软件;所述三维可视化分析采用Avizo9.0软件。In step 5), the finite element software is Abaqus software; the three-dimensional visualization analysis adopts Avizo9.0 software.
步骤8)中,纤维的失效准则采用最大应力准则;聚乳酸基体失效准则采用广义的Mises失效准则;对复合材料中纤维-基体界面采用二次应力准则。In step 8), the maximum stress criterion is adopted as the failure criterion of the fiber; the generalized Mises failure criterion is adopted as the failure criterion of the polylactic acid matrix; and the secondary stress criterion is adopted for the fiber-matrix interface in the composite material.
与现有技术相比,本发明的有益效果为:Compared with the prior art, the beneficial effects of the present invention are:
1)本发明围绕我国“碳达峰”、“碳中和”建设目标,针对当前节能减排和环境保护的时代需求,进行面向低碳汽车的植物纤维/聚乳酸绿色复合材料老化性能多尺度预测方法研究,探索农作物原料化新路径,破解石油能源危机困局,为植物纤维/聚乳酸绿色复合材料在汽车零部件领域的进一步应用提供参考和基础。1) The present invention focuses on the construction goals of “carbon peaking” and “carbon neutrality” in my country, aiming at the current needs of energy conservation, emission reduction and environmental protection, and conducts multi-scale aging performance of plant fiber/polylactic acid green composite materials for low-carbon vehicles. Research on forecasting methods, explore new paths for crop raw materials, solve the dilemma of oil and energy crisis, and provide reference and foundation for the further application of plant fiber/polylactic acid green composite materials in the field of auto parts.
2)目前关于短植物纤维/聚乳酸复合材料的湿热老化研究更多关注水分扩散理论模型的建立和完善、以及结构内应力分析等,尚未形成系统有效的老化后力学性能预测模型及方法。本发明从考虑各组分性能的细观建模出发,构建依赖吸水率变量的纤维、基体以及界面性能修正函数,实现多个因素耦合影响下的短植物纤维/聚乳酸复合材料老化后弹性性能及失效强度预测,为复合材料的性能预测提供新思路。2) At present, the research on the damp-heat aging of short plant fiber/polylactic acid composites pays more attention to the establishment and improvement of the theoretical model of water diffusion and the analysis of internal stress of the structure, etc., but a systematic and effective prediction model and method for mechanical properties after aging has not been formed. The invention starts from the mesoscopic modeling considering the performance of each component, constructs the fiber, matrix and interface performance correction function dependent on the water absorption rate variable, and realizes the elastic performance of the short plant fiber/polylactic acid composite material after aging under the coupling influence of multiple factors and failure strength prediction, which provides a new idea for the performance prediction of composite materials.
3)在长期服役过程中,受到环境、载荷、磨损等因素的影响,汽车零部件结构存在不同程度的老化现象,需要定期进行更新和维护。相较于金属和传统塑料等材质,植物纤维/聚乳酸复合材料对服役工况更加敏感,老化现象更加明显。如果超过部件服役寿命,则可能丧失使用功能,甚至造成严重的安全事故。而未达到使用寿命极限就提前更新维护,同样会造成人力和材料的极大浪费,显著增加使用成本。本发明所提供的老化性能模拟及寿命预测方法能够有效支撑绿色复合材料零部件结构设计理论,为制定合理的维护周期及方式提供参考,从而减少经济成本。3) In the long-term service process, affected by the environment, load, wear and other factors, the structure of auto parts has different degrees of aging, which needs to be updated and maintained regularly. Compared with materials such as metal and traditional plastics, the plant fiber/PLA composite material is more sensitive to service conditions, and the aging phenomenon is more obvious. If the service life of the components is exceeded, the use function may be lost, and even serious safety accidents may be caused. However, updating and maintaining in advance before reaching the limit of service life will also cause great waste of manpower and materials, and significantly increase the cost of use. The aging performance simulation and life prediction method provided by the invention can effectively support the structural design theory of green composite material parts, and provide a reference for formulating reasonable maintenance cycles and methods, thereby reducing economic costs.
附图说明Description of drawings
图1是本发明基于细观建模的植物纤维/聚乳酸复合材料湿热老化性能多尺度预测方法的工艺流程图。Fig. 1 is a process flow diagram of the multi-scale prediction method of the wet-heat aging performance of plant fiber/polylactic acid composite material based on mesoscopic modeling of the present invention.
具体实施方式Detailed ways
对本发明的基于细观建模的植物纤维/聚乳酸复合材料湿热老化性能多尺度预测方法作进一步详细说明:The multi-scale prediction method for wet-heat aging performance of plant fiber/polylactic acid composite material based on mesoscopic modeling of the present invention is further described in detail:
一种基于细观建模的植物纤维/聚乳酸复合材料湿热老化性能多尺度预测方法,如图1所示,包括以下步骤:A multi-scale prediction method of plant fiber/polylactic acid composites based on mesoscopic modeling for damp-heat aging properties, as shown in Figure 1, includes the following steps:
1)植物纤维/聚乳酸复合材料制备及老化试验:本发明选用黄麻纤维作为填充相,选用可降解聚乳酸作为基体材料,将黄麻纤维和聚乳酸颗粒在真空干燥箱中去除水分,用碱/硅烷偶联剂处理黄麻纤维,并按照1:9的比例(本实施例中所述比例仅仅为了说明模型建立,本发明适用于任何比例)将黄麻纤维和聚乳酸颗粒充分混合得到黄麻纤维/聚乳酸复合材料。通过双螺杆挤出机制备黄麻纤维/聚乳酸颗粒,经烘干处理后,再通过塑料注射成型机制备所需要的试样;同时采用上述方法分别制备纯黄麻纤维试样和聚乳酸试样,分别选取25℃、40℃、55℃水浴作为老化工况,对聚乳酸、黄麻纤维和复合材料进行湿热加速老化,老化完成后对试样干燥;1) plant fiber/polylactic acid composite material preparation and aging test: the present invention selects jute fiber as filling phase, selects degradable polylactic acid as matrix material, removes moisture from jute fiber and polylactic acid particles in a vacuum drying oven, and uses alkali/ The jute fiber is treated with a silane coupling agent, and the jute fiber and the polylactic acid particles are fully mixed according to the ratio of 1:9 (the ratio described in this example is only for illustrating the model establishment, and the present invention is applicable to any ratio) to obtain the jute fiber/polylactic acid particle. Lactic acid composite. The jute fiber/polylactic acid particles were prepared by a twin-screw extruder, and after drying, the required samples were prepared by a plastic injection molding machine; at the same time, pure jute fiber samples and polylactic acid samples were prepared by the above method, respectively. Select 25 ℃, 40 ℃, 55 ℃ water baths as aging conditions, respectively, carry out wet heat accelerated aging on polylactic acid, jute fiber and composite materials, and dry the samples after aging;
2)建立不同温度下的各老化材料吸水率随老化时间变化规律函数:通过电子秤称重法计算纯聚乳酸基体和复合材料的吸水率,通过定期监测黄麻纤维/聚乳酸复合材料、黄麻纤维和纯聚乳酸三者的平均吸水率,直到达到Fickian式中的平衡吸水率M∞,复合材料和纯聚乳酸材料中的吸水率M按照ASTM D5229标准来计算。2) Establish the law function of the water absorption rate of each aging material at different temperatures with aging time: calculate the water absorption rate of pure PLA matrix and composite materials by electronic scale weighing method, and regularly monitor jute fiber/polylactic acid composite material, jute fiber and the average water absorption of pure PLA until reaching the equilibrium water absorption M ∞ in Fickian formula, the water absorption M in the composite material and pure PLA material is calculated according to the ASTM D5229 standard.
纤维的吸水率通过下式获得The water absorption of the fiber is obtained by the following formula
ΔMc=ΔMf×Vf+ΔWm(1-Vf)ΔM c =ΔM f ×V f +ΔW m (1-V f )
通过实验分别得到黄麻纤维/聚乳酸复合材料、黄麻纤维和纯聚乳酸三者不同温度下的饱和吸水率,Through experiments, the saturated water absorption rates of jute fiber/PLA composite material, jute fiber and pure PLA at different temperatures were obtained, respectively.
表1不同温度下的各老化材料的饱和吸水率Table 1 Saturated water absorption of each aged material at different temperatures
整合试验数据分别得到不同老化温度下复合材料、聚乳酸基体和黄麻纤维吸水率随老化时间变化的规律;Integrate the test data to obtain the changing law of the water absorption of the composite material, PLA matrix and jute fiber with aging time under different aging temperatures;
3)建立不同温度下的各组分强度随老化时间变化规律函数:选取不同老化温度和老化时间的材料进行力学性能试验,对复合材料进行单纤维拔出试验、对纤维进行单纤维拉伸试验、对聚乳酸基体进行哑铃拉伸试验,分别得到三个组分达到饱和吸水率时的强度下降比例,如下表所示:3) Establish a function of the variation law of the strength of each component at different temperatures with aging time: select materials with different aging temperatures and aging times for mechanical property tests, single-fiber pull-out tests for composite materials, and single-fiber tensile tests for fibers , Carry out the dumbbell tensile test on the polylactic acid matrix, and obtain the strength reduction ratio of the three components when they reach the saturated water absorption rate, as shown in the following table:
表2不同温度下各组分强度变化情况(相比于初次强度下降百分比)Table 2 Changes in the strength of each component at different temperatures (percentage decrease compared to the initial strength)
通过对试验获得的数据整合分别得到基体-纤维界面的强度随老化温度、老化时间变化的规律,纤维强度随老化温度、老化时间变化的规律,聚乳酸基体强度随老化温度、老化时间变化的规律;By integrating the data obtained from the test, the law of the strength of the matrix-fiber interface changing with aging temperature and aging time, the law of fiber strength changing with aging temperature and aging time, and the changing law of polylactic acid matrix strength with aging temperature and aging time were obtained. ;
4)分别建立聚乳酸基体、纤维、复合材料各组分强度与吸水率、温度之间关系函数:对前面获得的数据进行整合,通过拟合数据分别得到植物纤维、聚乳酸基体和纤维-基体界面三者强度随吸水率和老化温度变化的函数4) Establish the relationship function between the strength of each component of the polylactic acid matrix, fiber and composite material, water absorption, and temperature: Integrate the data obtained before, and obtain the plant fiber, polylactic acid matrix and fiber-matrix by fitting the data. Interfacial triad strength as a function of water absorption and aging temperature
Sf(M,T)=14.82787+0.9456T-0.03143M-0.01163T2-0.02714M2-0.00991MTS f (M,T)=14.82787+0.9456T-0.03143M-0.01163T 2 -0.02714M 2 -0.00991MT
Sm(M,T)=42.16438+0.49302T-17.52793M-0.0065T2-24.68074M2+0.37515MTS m (M,T)=42.16438+0.49302T-17.52793M-0.0065T 2 -24.68074M 2 +0.37515MT
Si(M,T)=-1.23196+3.11211T+9.87753M-0.03748T2-1.82922M2-0.3378MTS i (M,T)=-1.23196+3.11211T+9.87753M-0.03748T 2 -1.82922M 2 -0.3378MT
5)复合材料细观RVE模型建立:从未老化植物纤维/聚乳酸复合材料中采样,进行X射线断层扫描,获得包含材料微观结构信息的断层灰度图像;采用Avizo9.0软件对植物纤维/聚乳酸复合材料的断层图像进行三维可视化分析。为了去除图像中杂质的影响,首先采用中值滤波技术对原始图像进行过滤简化处理;然后采用Avizo 9.0软件中的修复功能对天然纤维表面孔洞进行闭合清理,并对修复完成的植物纤维体素模型进行表面平滑处理;最后将重构的几何模型导入有限元分析软件Abaqus软件,建立复合材料的细观RVE模型;5) Establishment of the microscopic RVE model of the composite material: sampling from the unaged plant fiber/polylactic acid composite material, performing X-ray tomography to obtain a tomographic gray image containing the microstructure information of the material; using Avizo9.0 software to analyze the plant fiber/polylactic acid composite material. The tomographic images of PLA composites were analyzed for 3D visualization. In order to remove the influence of impurities in the image, the original image was filtered and simplified by median filtering technology; then, the repair function in Avizo 9.0 software was used to close and clean the pores on the surface of natural fibers, and the repaired plant fiber voxel model was The surface is smoothed; finally, the reconstructed geometric model is imported into the finite element analysis software Abaqus software to establish the meso RVE model of the composite material;
6)环境退化因子定义与引入:编写纤维、基体、界面模拟单元本构关系,对复合材料细观RVE模型进行属性定义。构建基于强度变化函数的各组分环境退化因子(D)函数6) Definition and introduction of environmental degradation factors: Compose the constitutive relationship of fiber, matrix, and interface simulation element, and define the properties of the meso-scale RVE model of composite materials. Constructing the environmental degradation factor (D) function of each component based on the intensity change function
分别得到纤维的退化因子函数The degradation factor functions of the fibers are obtained respectively
Df(M,T)=0.46337+0.02955T-0.00098m-0.00036T2-0.00085M2+0.00031MT,D f (M,T)=0.46337+0.02955T-0.00098m-0.00036T 2 -0.00085M 2 +0.00031MT,
基体的退化因子函数Degradation factor function of the matrix
Dm(M,T)=0.84329+0.00986T-0.35056m-0.00013T2-0.49361M2+0.00750MT,D m (M,T)=0.84329+0.00986T-0.35056m-0.00013T 2 -0.49361M 2 +0.00750MT,
纤维-基体界面的退化因子函数Degradation factor function of fiber-matrix interface
Di(M,T)=-0.02124+0.05366T+0.17030m-0.00065T2-0.03154M2-0.00582MT,D i (M,T)=-0.02124+0.05366T+0.17030m-0.00065T 2 -0.03154M 2 -0.00582MT,
并将退化因子引入对应本构,实现基于老化效应的各组分本构参数修正;The degradation factor is introduced into the corresponding constitutive to realize the correction of the constitutive parameters of each component based on the aging effect;
7)复合材料弹性性能计算:通过引入退化因子,对各组分弹性性能参数进行修正;在周期性边界条件下,对RVE模型施加线性不相关位移载荷并进行仿真,获得不同温度和吸水率条件下的宏观复合材料弹性性能;7) Calculation of elastic properties of composite materials: By introducing degradation factors, the elastic performance parameters of each component are corrected; under periodic boundary conditions, linear uncorrelated displacement loads are applied to the RVE model and simulated to obtain different temperature and water absorption conditions The elastic properties of macroscopic composites under
8)复合材料失效强度计算:由于在细观状态下研究尺度很小,所以在细观模型中将短植物纤维视为各向同性材料,将植物纤维表示为线弹性本构,纤维的失效准则采用最大应力准则;聚乳酸基体通常考虑为各向同性弹塑性材料,失效准则采用广义的Mises失效准则;对纤维-基体界面采用内聚力本构,采用二次应力准则;对失效参数进行修正,并进行数值模拟,得到不同吸水率和温度下的宏观复合材料失效强度;8) Calculation of the failure strength of composite materials: Since the research scale is very small in the mesoscopic state, the short plant fiber is regarded as an isotropic material in the mesoscopic model, and the plant fiber is expressed as a linear elastic constitutive, the failure criterion of the fiber The maximum stress criterion is adopted; the PLA matrix is usually considered as an isotropic elastic-plastic material, and the generalized Mises failure criterion is adopted for the failure criterion; the cohesive force constitutive and the quadratic stress criterion are adopted for the fiber-matrix interface; Numerical simulation was carried out to obtain the failure strength of macroscopic composites under different water absorption rates and temperatures;
9)宏观复合材料湿热老化性能预测:建立复合材料拉伸、三点弯曲试件宏观模型,其中,复合材料拉伸标准采用ISO 527-2,三点弯曲标准采用EN ISO14125,将基于细观RVE模型仿真得到的依赖吸水率变化的复合材料宏观弹性性能和失效强度,代入到宏观模型中,设置湿热环境进行有限元分析,对湿热老化工况下不同老化程度的植物纤维/聚乳酸复合材料力学性能进行预测,结果显示拉伸实验误差为7.86%,三点弯曲试验误差为8.13%,这表明本发明预测方法在植物纤维/聚乳酸复合材料湿热老化性能预测上具有一定的实际意义。9) Prediction of the wet-heat aging performance of macroscopic composite materials: establish a macroscopic model of composite tensile and three-point bending specimens. Among them, the tensile standard of composite materials adopts ISO 527-2, and the three-point bending standard adopts EN ISO14125, which will be based on the microscopic RVE The macro-elastic properties and failure strengths of composites that depend on changes in water absorption obtained from the model simulation are substituted into the macro-model, and the damp-heat environment is set for finite element analysis. The performance is predicted, and the results show that the tensile test error is 7.86%, and the three-point bending test error is 8.13%, which shows that the prediction method of the present invention has certain practical significance in predicting the wet heat aging performance of the plant fiber/polylactic acid composite material.
Claims (4)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210102395.0A CN114417634B (en) | 2022-01-27 | 2022-01-27 | Plant fiber/polylactic acid composite material wet heat aging performance multi-scale prediction method based on mesoscopic modeling |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210102395.0A CN114417634B (en) | 2022-01-27 | 2022-01-27 | Plant fiber/polylactic acid composite material wet heat aging performance multi-scale prediction method based on mesoscopic modeling |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114417634A true CN114417634A (en) | 2022-04-29 |
CN114417634B CN114417634B (en) | 2024-02-20 |
Family
ID=81280243
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210102395.0A Active CN114417634B (en) | 2022-01-27 | 2022-01-27 | Plant fiber/polylactic acid composite material wet heat aging performance multi-scale prediction method based on mesoscopic modeling |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114417634B (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190384878A1 (en) * | 2018-06-14 | 2019-12-19 | The United States Of America, As Represented By The Secretary Of The Navy | Fibrous Composite Failure Criteria with Material Degradation for Finite Element Solvers |
CN113420376A (en) * | 2021-06-17 | 2021-09-21 | 吉林大学 | Multi-scale-based impact-resistant mechanical property simulation method for carbon fiber composite material |
-
2022
- 2022-01-27 CN CN202210102395.0A patent/CN114417634B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190384878A1 (en) * | 2018-06-14 | 2019-12-19 | The United States Of America, As Represented By The Secretary Of The Navy | Fibrous Composite Failure Criteria with Material Degradation for Finite Element Solvers |
CN113420376A (en) * | 2021-06-17 | 2021-09-21 | 吉林大学 | Multi-scale-based impact-resistant mechanical property simulation method for carbon fiber composite material |
Non-Patent Citations (2)
Title |
---|
吕新颖;江龙;闫亮;王荣国;刘文博;: "碳纤维复合材料湿热性能研究进展", 玻璃钢/复合材料, no. 03, 28 May 2009 (2009-05-28) * |
胡丽娟;张少睿;李大永;苌群峰;彭颖红;: "细观参数对纤维增强金属基复合材料宏细观力学性能的影响", 上海交通大学学报, no. 03, 15 March 2008 (2008-03-15) * |
Also Published As
Publication number | Publication date |
---|---|
CN114417634B (en) | 2024-02-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chi et al. | Constitutive modeling of steel-polypropylene hybrid fiber reinforced concrete using a non-associated plasticity and its numerical implementation | |
CN104537259B (en) | A kind of Forecasting Methodology of the transverse isotropy fibrous material performance based on meso mechanical model | |
Ogierman et al. | A study on fiber orientation influence on the mechanical response of a short fiber composite structure | |
He et al. | Elastic properties and multi-scale design of long carbon fiber nonwoven reinforced plane-based isotropic composite | |
CN109858171B (en) | Prediction method for stress-strain response and strength of braided ceramic matrix composites | |
Gao et al. | A new probability model of residual strength of material based on interference theory | |
Wang et al. | Micromechanical modelling of the progressive failure in unidirectional composites reinforced with bamboo fibres | |
CN117804896A (en) | A method for analyzing the mesoscopic non-periodic structural strength of ceramic matrix composite materials | |
Zuccarello et al. | Environmental aging effects on high-performance biocomposites reinforced by sisal fibers | |
Wang et al. | Influence of moisture absorption on mechanical properties and damage mechanisms of three-dimensional six-directional braided composites under hydrostatic pressure | |
CN116386776A (en) | A method for predicting high and low temperature mechanical properties of flax fiber composites | |
CN114417634A (en) | A multi-scale prediction method for damp-heat aging properties of plant fiber/polylactic acid composites based on mesoscopic modeling | |
Califano et al. | Theoretical approach to the study of fatigue of composites under spectrum loading | |
Li et al. | Temperature-dependent strength modeling of fiber-reinforced composites considering critical properties evolution | |
Zhang et al. | Investigating the role of fibre-matrix interfacial degradation on the ageing process of carbon fibre-reinforced polymer under hydrothermal conditions | |
Califano et al. | Analysis of a phenomenological model for fatigue of composite materials | |
Kang et al. | Void stability process window and parametric optimization for filament‐wound composite riser | |
Huang et al. | Mathematical study of the mechanical properties of NOL rings with different winding process parameters | |
Despringre et al. | Micromechanical Fatigue Visco-Damage Model for Short Glass Fiber Reinforced Polyamide-66 | |
Despringre et al. | Multi-scale viscoelastic damage model of short glass fiber reinforced thermoplastics under fatigue loading | |
Zhang et al. | Experimental characterization and numerical analysis of 3D woven preforms torsional deformation | |
Wang et al. | Modeling and Simulation of Epoxy/Natural Fiber Composites | |
Al-Saadi et al. | Anexperimental investigation of moisture effect on fatigue behavior of composite materials | |
Sun et al. | Fatigue behavior and modeling of chopped carbon fiber reinforced sheet molding compound composites | |
Mohamed | Micromechanical Modelling of Random Short Fiber Reinforced Polymer Composites With Progressive Debonding Damage |
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 |