CN108595781A - 一种复合材料成型后纤维和基体的弹性参数识别方法 - Google Patents

一种复合材料成型后纤维和基体的弹性参数识别方法 Download PDF

Info

Publication number
CN108595781A
CN108595781A CN201810287409.4A CN201810287409A CN108595781A CN 108595781 A CN108595781 A CN 108595781A CN 201810287409 A CN201810287409 A CN 201810287409A CN 108595781 A CN108595781 A CN 108595781A
Authority
CN
China
Prior art keywords
matrix
fiber
parameters
formula
displacement
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.)
Pending
Application number
CN201810287409.4A
Other languages
English (en)
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.)
Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CN201810287409.4A priority Critical patent/CN108595781A/zh
Publication of CN108595781A publication Critical patent/CN108595781A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design 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)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种复合材料成型后纤维和基体的弹性参数识别方法,包括以下步骤:S1:建立复合材料的细观模型,通过有限元方法对细观模型进行静态结构分析,提取静态位移数据,计算静态位移场对纤维和基体弹性参数的相对灵敏度;S2:将试验和有限元模拟的位移残差的二范数作为目标函数,利用迭代优化算法最小化目标函数。本发明提出了考虑组分材料性能的复合材料细观模型,通过结构位移响应对组分材料参数的相对灵敏度分析,采用优化迭代算法识别成型后复合材料纤维和基体的弹性参数,具有十分重要的工程意义。

Description

一种复合材料成型后纤维和基体的弹性参数识别方法
技术领域
本发明涉及复合材料静力学性能参数计算技术领域,特别是涉及一种复合材料成型后纤维和基体的弹性参数识别方法。
背景技术
复合材料具有高比强、高比模量、可设计性强等优点,广泛应用于航空航天、土木工程、机械工程等领域。在不同工程领域中,对复合材料的宏观力学性能有不同的要求,因此需要复合材料准确的组分性能参数为材料结构设计提供基础信息。复合材料成型后的材料性能与其对应的单组分材料性能相差较大,其可以归因于实验室标准的制作工艺和工厂制造水平之间的差距。在材料制作过程中,需要对进行热处理、压力成型等工艺,这些工艺会导致纤维和基体性能发生改变,使制作出来的复合材料成品在性能确定时出现困难。
发明内容
发明目的:本发明的目的是提供一种复合材料成型后纤维和基体的弹性参数识别方法,能够有效识别成型后复合材料纤维和基体的弹性参数。
技术方案:为达到此目的,本发明采用以下技术方案:
本发明所述的复合材料成型后纤维和基体的弹性参数识别方法,包括以下步骤:
S1:建立复合材料的细观模型,通过有限元方法对细观模型进行静态结构分析,提取静态位移数据,计算静态位移场对纤维和基体弹性参数的相对灵敏度;
S2:将试验和有限元模拟的位移残差的二范数作为目标函数,利用迭代优化算法最小化目标函数。
进一步,所述步骤S1中,计算静态位移场对纤维和基体弹性参数的相对灵敏度具体包括以下步骤:
S1.1:根据先验信息确定材料参数初值,输入与实验相符的边界条件和载荷条件,根据经典混合理论得到细观组分材料参数与宏观材料性能的对应关系,根据经典混合理论和胡克定律,建立复合材料细观弹性本构,如式(1)和(2)所示:
σ=(kmDm+kfDf)·ε (1)
km+kf=1 (2)
其中,σ为材料应力,ε为材料应变,Dm为材料基体的弹性系数矩阵,Df为材料纤维的弹性系数矩阵,km为材料基体的体积分数;
S1.2:根据虚功原理确定的材料单元刚度矩阵以及单元刚度矩阵对材料参数的偏导,其中单元刚度矩阵的公式如式(3)所示:
其中,Ke为单元刚度矩阵,B为单元应变矩阵,Ve为细观模型单元的体积,D 为弹性系数矩阵;单元刚度矩阵对材料参数的偏导如式(4)所示:
其中,p为复合材料需要识别的弹性参数;
S1.3:推导复合材料整体刚度矩阵与各组分材料单元刚度矩阵的关系,并计算整体刚度矩阵对组分材料参数的偏导,其中材料整体刚度矩阵由纤维和基体的单元刚度矩阵叠加而成,如式(5)所示:
其中,K为材料整体刚度矩阵,G为单元和结构节点自由度的转换矩阵,为材料基体的刚度矩阵,为材料纤维的刚度矩阵;
整体刚度矩阵对组分材料参数的偏导如式(6)所示:
S1.4:根据离散有限元方程得到有限元求解的静态位移和灵敏度矩阵,其中有限元方程如式(7)所示:
K·u=P (7)
其中,u为节点位移向量,P为节点载荷向量;
有限元求解的静态位移u如式(8)所示:
静态位移对材料参数的偏导S如式(9)所示:
其中,S为位移对待识别组分参数的灵敏度矩阵;
S1.5:通过无量纲处理得到相对灵敏度矩阵,相对灵敏度矩阵的公式如式(10) 所示:
S'=u-1Sp (10)
其中,S'为相对灵敏度矩阵;
进一步,所述步骤S2具体包括以下步骤:
S2.1:将试验和有限元模拟的位移数据残差的二范数作为目标函数,如式(11) 所示:
F(p)=||u*-u||2 (11)
其中,F(p)为目标函数,u*表示实验的测量位移,u为有限元求解的静态位移;
S2.2:采用迭代优化算法,对目标函数进行优化,判定是否满足收敛准则:若满足,结束识别过程;若不满足,则更新材料参数重新计算。
有益效果:本发明公开了一种复合材料成型后纤维和基体的弹性参数识别方法,提出了考虑组分材料性能的复合材料细观模型,通过结构位移响应对组分材料参数的相对灵敏度分析,采用优化迭代算法识别成型后复合材料纤维和基体的弹性参数,具有十分重要的工程意义。
附图说明
图1为本发明具体实施方式中方法的流程图;
图2为本发明具体实施方式中的复合材料单胞实体模型;
图3为本发明具体实施方式中的复合材料单胞平面示意图;
图4为本发明具体实施方式中参数识别过程示意图。
具体实施方式
下面结合具体实施方式和附图对本发明的技术方案作进一步的介绍。
本具体实施方式公开了一种复合材料成型后纤维和基体的弹性参数识别方法,如图1所示,以纤维增强复合材料的单胞模型为例进行说明,其中纤维为横观各向同性材料,主方向为x方向,其中E11=220GPa,E22=13.8GPa,μ12=0.2,G12=9GPa,μ23=0.25;基体为各向同性材料,E=250GPa,μ=0.25。该方法包括以下步骤:
S1:建立复合材料的细观模型,单胞模型的示意图如图2、3所示,通过有限元方法对细观模型进行静态结构分析,提取静态位移数据,计算静态位移场对纤维和基体弹性参数的相对灵敏度。
S2:将试验和有限元模拟的位移残差的二范数作为目标函数,利用迭代优化算法最小化目标函数。
步骤S1中,计算静态位移场对纤维和基体弹性参数的相对灵敏度具体包括以下步骤:
S1.1:根据先验信息确定材料参数初值:纤维的材料初值为:E11=120GPa, E22=8GPa,μ12=0.3,G12=12GPa,μ23=0.47;基体的材料初值为:E=250GPa,μ=0.35。输入与实验相符的边界条件和载荷条件,给模型施加标准三点弯实验的约束,并在单胞模型的中间施加z向为100N的力,根据经典混合理论得到细观组分材料参数与宏观材料性能的对应关系,根据经典混合理论和胡克定律,建立复合材料细观弹性本构,如式(1)和(2)所示:
σ=(kmDm+kfDf)·ε (1)
km+kf=1 (2)
其中,σ为材料应力,ε为材料应变,Dm为材料基体的弹性系数矩阵,Df为材料纤维的弹性系数矩阵,km为材料基体的体积分数;Dm可通过基体的柔度矩阵 Sm取逆得到,Df可通过纤维的柔度矩阵Sf取逆得到,纤维的柔度矩阵Sf如式(3) 所示。
式(3)中,E11、E22、μ12、μ23、G12和G23均为材料弹性参数。
S1.2:根据虚功原理确定的材料单元刚度矩阵以及单元刚度矩阵对材料参数的偏导,其中单元刚度矩阵的公式如式(4)所示:
其中,Ke为单元刚度矩阵,B为单元应变矩阵,Ve为细观模型单元的体积,D 为弹性系数矩阵;单元刚度矩阵对材料参数的偏导如式(5)所示:
其中,p为复合材料需要识别的弹性参数;
S1.3:推导复合材料整体刚度矩阵与各组分材料单元刚度矩阵的关系,并计算整体刚度矩阵对组分材料参数的偏导,其中材料整体刚度矩阵由纤维和基体的单元刚度矩阵叠加而成,如式(6)所示:
其中,K为材料整体刚度矩阵,G为单元和结构节点自由度的转换矩阵,为材料基体的刚度矩阵,为材料纤维的刚度矩阵;
整体刚度矩阵对组分材料参数的偏导如式(7)所示:
S1.4:根据离散有限元方程得到有限元求解的静态位移和灵敏度矩阵,其中有限元方程如式(8)所示:
K·u=P (8)
其中,u为节点位移向量,P为节点载荷向量;
有限元求解的静态位移u如式(9)所示:
静态位移对材料参数的偏导S如式(10)所示:
其中,S为位移对待识别组分参数的灵敏度矩阵;
S1.5:通过无量纲处理得到相对灵敏度矩阵,相对灵敏度矩阵的公式如式(11) 所示:
S'=u-1Sp (11)
其中,S'为相对灵敏度矩阵;
步骤S2具体包括以下步骤:
S2.1:将试验和有限元模拟的位移数据残差的二范数作为目标函数,如式(12) 所示:
F(p)=||u*-u||2 (12)
其中,F(p)为目标函数,u*表示实验的测量位移,u为有限元求解的静态位移;
S2.2:采用迭代优化算法,对目标函数进行优化,判定是否满足收敛准则:若满足,结束识别过程;若不满足,则更新材料参数重新计算。复合材料纤维和基体各参数的识别过程示意图如图4所示。

Claims (3)

1.一种复合材料成型后纤维和基体的弹性参数识别方法,其特征在于:包括以下步骤:
S1:建立复合材料的细观模型,通过有限元方法对细观模型进行静态结构分析,提取静态位移数据,计算静态位移场对纤维和基体弹性参数的相对灵敏度;
S2:将试验和有限元模拟的位移残差的二范数作为目标函数,利用迭代优化算法最小化目标函数。
2.根据权利要求1所述的复合材料成型后纤维和基体的弹性参数识别方法,其特征在于:所述步骤S1中,计算静态位移场对纤维和基体弹性参数的相对灵敏度具体包括以下步骤:
S1.1:根据先验信息确定材料参数初值,输入与实验相符的边界条件和载荷条件,根据经典混合理论得到细观组分材料参数与宏观材料性能的对应关系,根据经典混合理论和胡克定律,建立复合材料细观弹性本构,如式(1)和(2)所示:
σ=(kmDm+kfDf)·ε (1)
km+kf=1 (2)
其中,σ为材料应力,ε为材料应变,Dm为材料基体的弹性系数矩阵,Df为材料纤维的弹性系数矩阵,km为材料基体的体积分数;
S1.2:根据虚功原理确定的材料单元刚度矩阵以及单元刚度矩阵对材料参数的偏导,其中单元刚度矩阵的公式如式(3)所示:
其中,Ke为单元刚度矩阵,B为单元应变矩阵,Ve为细观模型单元的体积,D为弹性系数矩阵;单元刚度矩阵对材料参数的偏导如式(4)所示:
其中,p为复合材料需要识别的弹性参数;
S1.3:推导复合材料整体刚度矩阵与各组分材料单元刚度矩阵的关系,并计算整体刚度矩阵对组分材料参数的偏导,其中材料整体刚度矩阵由纤维和基体的单元刚度矩阵叠加而成,如式(5)所示:
其中,K为材料整体刚度矩阵,G为单元和结构节点自由度的转换矩阵,为材料基体的刚度矩阵,为材料纤维的刚度矩阵;
整体刚度矩阵对组分材料参数的偏导如式(6)所示:
S1.4:根据离散有限元方程得到有限元求解的静态位移和灵敏度矩阵,其中有限元方程如式(7)所示:
K·u=P (7)
其中,u为节点位移向量,P为节点载荷向量;
有限元求解的静态位移u如式(8)所示:
静态位移对材料参数的偏导S如式(9)所示:
其中,S为位移对待识别组分参数的灵敏度矩阵;
S1.5:通过无量纲处理得到相对灵敏度矩阵,相对灵敏度矩阵的公式如式(10)所示:
S'=u-1Sp (10)
其中,S'为相对灵敏度矩阵。
3.根据权利要求1所述的复合材料成型后纤维和基体的弹性参数识别方法,其特征在于:所述步骤S2具体包括以下步骤:
S2.1:将试验和有限元模拟的位移数据残差的二范数作为目标函数,如式(11)所示:
F(p)=||u*-u||2 (11)
其中,F(p)为目标函数,u*表示实验的测量位移,u为有限元求解的静态位移;
S2.2:采用迭代优化算法,对目标函数进行优化,判定是否满足收敛准则:
若满足,结束识别过程;若不满足,则更新材料参数重新计算。
CN201810287409.4A 2018-03-30 2018-03-30 一种复合材料成型后纤维和基体的弹性参数识别方法 Pending CN108595781A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810287409.4A CN108595781A (zh) 2018-03-30 2018-03-30 一种复合材料成型后纤维和基体的弹性参数识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810287409.4A CN108595781A (zh) 2018-03-30 2018-03-30 一种复合材料成型后纤维和基体的弹性参数识别方法

Publications (1)

Publication Number Publication Date
CN108595781A true CN108595781A (zh) 2018-09-28

Family

ID=63625218

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810287409.4A Pending CN108595781A (zh) 2018-03-30 2018-03-30 一种复合材料成型后纤维和基体的弹性参数识别方法

Country Status (1)

Country Link
CN (1) CN108595781A (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110008530A (zh) * 2019-03-15 2019-07-12 东南大学 一种空间柔性复合材料分布式概率建模方法
CN110059368A (zh) * 2019-03-27 2019-07-26 东南大学 一种空间柔性复合材料的参数识别方法
CN111783347A (zh) * 2020-07-27 2020-10-16 内蒙古工业大学 一种复合纤维风力机叶片分区细观铺层参数优化方法
TWI758953B (zh) * 2020-11-16 2022-03-21 中國鋼鐵股份有限公司 耐火泥彈性模數的測量方法
CN116822282A (zh) * 2023-06-16 2023-09-29 南京航空航天大学 一种湿热环境下2.5维复合材料的动力学模型实现方法
CN117174217A (zh) * 2023-10-31 2023-12-05 中国矿业大学 一种满足热传导性能要求的复合材料纤维形状优化方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103871059A (zh) * 2014-03-13 2014-06-18 南京航空航天大学 纤维增强复合材料等效弹性参数的计算方法
US9081921B2 (en) * 2011-06-16 2015-07-14 Sumitomo Rubber Industries, Ltd. Method for simulating rubber compound
CN106066913A (zh) * 2016-05-31 2016-11-02 西北工业大学 复杂复合材料结构等效材料性能多尺度计算方法
CN107357992A (zh) * 2017-07-13 2017-11-17 东南大学 基于聚类分析的复合材料结构有限元模型修正方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9081921B2 (en) * 2011-06-16 2015-07-14 Sumitomo Rubber Industries, Ltd. Method for simulating rubber compound
CN103871059A (zh) * 2014-03-13 2014-06-18 南京航空航天大学 纤维增强复合材料等效弹性参数的计算方法
CN106066913A (zh) * 2016-05-31 2016-11-02 西北工业大学 复杂复合材料结构等效材料性能多尺度计算方法
CN107357992A (zh) * 2017-07-13 2017-11-17 东南大学 基于聚类分析的复合材料结构有限元模型修正方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
秦福溶 等: "基于灵敏度分析的复合材料组分参数识别方法", 《复合材料学报》 *
费庆国 等: "高温下编织复合材料热相关参数识别方法研究", 《力学学报》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110008530A (zh) * 2019-03-15 2019-07-12 东南大学 一种空间柔性复合材料分布式概率建模方法
CN110059368A (zh) * 2019-03-27 2019-07-26 东南大学 一种空间柔性复合材料的参数识别方法
CN110059368B (zh) * 2019-03-27 2021-02-09 东南大学 一种空间柔性复合材料的参数识别方法
CN111783347A (zh) * 2020-07-27 2020-10-16 内蒙古工业大学 一种复合纤维风力机叶片分区细观铺层参数优化方法
CN111783347B (zh) * 2020-07-27 2023-06-09 内蒙古工业大学 一种复合纤维风力机叶片分区细观铺层参数优化方法
TWI758953B (zh) * 2020-11-16 2022-03-21 中國鋼鐵股份有限公司 耐火泥彈性模數的測量方法
CN116822282A (zh) * 2023-06-16 2023-09-29 南京航空航天大学 一种湿热环境下2.5维复合材料的动力学模型实现方法
CN117174217A (zh) * 2023-10-31 2023-12-05 中国矿业大学 一种满足热传导性能要求的复合材料纤维形状优化方法
CN117174217B (zh) * 2023-10-31 2024-01-26 中国矿业大学 一种满足热传导性能要求的复合材料纤维形状优化方法

Similar Documents

Publication Publication Date Title
CN108595781A (zh) 一种复合材料成型后纤维和基体的弹性参数识别方法
CN110210103B (zh) 一种多相复合材料力学行为的多尺度模拟方法
Zhu et al. Buckling analysis of Euler–Bernoulli beams using Eringen’s two-phase nonlocal model
Zhang et al. A non-classical Mindlin plate finite element based on a modified couple stress theory
KR101383663B1 (ko) 가상 테스트에 근거한 파라미터화 재료 및 성능 특성
CN110389070B (zh) 确定高强度和高弹性点阵材料的方法
CN107273609A (zh) 一种基于Kriging模型齿轮传动可靠性评估方法
CN109271655B (zh) 一种基于非对称有限元算法的材料尺度效应分析方法
Sahmani et al. Nonlinear instability of hydrostatic pressurized microtubules surrounded by cytoplasm of a living cell including nonlocality and strain gradient microsize dependency
CN107330137A (zh) 一种金属材料切削本构模型参数的辨识方法
Sultanov Analysis of finite elasto-plastic strains. Medium kinematics and constitutive equations
CN105787167A (zh) 聚合物基复合材料有效应力松弛系数的预测方法
CN112749485A (zh) 一种晶格扰动模式下晶体材料理想强度的高通量计算方法
Pegios et al. Finite element static and stability analysis of gradient elastic beam structures
CN113065269A (zh) 复合材料多尺度并发模拟方法和系统
CN106777462A (zh) 电子组件振动应力的分析方法及系统
Salvadori Analytical integrations in 2D BEM elasticity
Shim et al. Integrated modeling framework for the interactions of plastic deformation, magnetic fields, and electrical circuits: Theory and applications to physics-informed real-time material monitoring
CN107563071A (zh) 一种综合环境应力的仿真方法
CN112989672B (zh) 适应于复杂应力变化的极小曲面梯度结构的构建方法
CN114388076A (zh) 基于深度学习的物体宏观物理属性预测方法及系统
Grbac et al. Hexahedral finite elements with enhanced fixed‐pole interpolation for linear static and vibration analysis of 3D micropolar continuum
CN102096737B (zh) 一种压电复合棒机-电耦合性能仿真模拟方法
Adly et al. Utilizing anisotropic Preisach-type models in the accurate simulation of magnetostriction
Vyas et al. Nonlinear geometric reduced order model for the response of a beam with a piezoelectric actuator

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20180928

RJ01 Rejection of invention patent application after publication