CN107239642A - 面向碳纤维增强复合材料钻削的多尺度数值仿真方法 - Google Patents

面向碳纤维增强复合材料钻削的多尺度数值仿真方法 Download PDF

Info

Publication number
CN107239642A
CN107239642A CN201710554685.8A CN201710554685A CN107239642A CN 107239642 A CN107239642 A CN 107239642A CN 201710554685 A CN201710554685 A CN 201710554685A CN 107239642 A CN107239642 A CN 107239642A
Authority
CN
China
Prior art keywords
msub
mrow
impaired loci
fiber
reinforced composite
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
Application number
CN201710554685.8A
Other languages
English (en)
Other versions
CN107239642B (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.)
Jiangsu Hangpu Guochuang Composite Material Research Institute Co ltd
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201710554685.8A priority Critical patent/CN107239642B/zh
Publication of CN107239642A publication Critical patent/CN107239642A/zh
Application granted granted Critical
Publication of CN107239642B publication Critical patent/CN107239642B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

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)
  • Manufacture Of Alloys Or Alloy Compounds (AREA)

Abstract

一种面向碳纤维增强复合材料钻削多尺度数值仿真方法,其特征是:首先,进行基体、纤维细观本构与宏观本构建模。其次,进行宏观损伤判断并输出损伤点的信息。第三,调用Python文件对损伤点进行细观建模分析。第四,利用细观建模分析的数据,计算出该点等效刚度并返回宏观模型。第五,进行多尺度数值仿真系统开发与模型评价。本发明充分的考虑了纤维和基体的区别,以及仿真的精度与效率。给出了碳纤维增强复合材料钻削多尺度数值仿真方法。

Description

面向碳纤维增强复合材料钻削的多尺度数值仿真方法
技术领域
本发明涉及一种机械加工技术,尤其是一种钻削仿真方法,具体的说就是一种碳纤维增强复合材料钻削多尺度数值仿真方法。
背景技术
航空航天产品中的CFRP通常需要与其它结构(复合材料、铝合金、钛合金等结构)进行铆接、螺接等机械连接,这就需要在航空航天产品装配过程中对CFRP构件进行大量的钻孔加工操作。由于CFRP是由纤维和基体组成的多相结构,具有非匀质性、各向异性、硬度高等特点,切削加工难度较大,钻削过程中极易出现入口处劈裂、出口处撕裂、毛边、层间分离、径向挤伤、微裂纹等损伤。从而降低机械连接质量,影响产品服役性能。国内外学者围绕CFRP切削力预测、工艺参数优化等方面,通过力学解析、数值模拟和实验观测展开了广泛的研究和探索随着航空航天制造工业对装配质量要求的进一步提高,对钻削损伤控制的要求越来越严格,因而CFRP切削机理、损伤的萌生和扩展规律等需要更进一步深入的研究。而有效的研究方法是对CFRP进行钻削模拟仿真研究。但由于复合材料钻削行为是一个囊括了几何信息、材料属性信息、应力应变场信息的动态变化的复杂系统,解析的方法很难提供令人满意的效果。并存在以下问题:
(1)采用单相均质材料等效CFRP,对于其中某一特定材料点而言,忽略了纤维和基体的区别,因此导致一些本该用纤维性质计算的单元被更弱的基体性质代替,从而提前失效;
(2)精度和效率难以达到让人满意的结果。
发明内容
本发明的目的是针对CFRP进行钻削模拟仿真中材料提前失效,精度效率不高的问题,发明一种面向碳纤维增强复合材料钻削多尺度数值仿真方法,它充分的考虑了纤维和基体的区别,以及仿真的精度与效率。给出了碳纤维增强复合材料钻削多尺度数值仿真方法。
本发明的技术方案是:
一种面向碳纤维增强复合材料钻削多尺度数值仿真方法,其特征是它包括以下步骤?
首先,进行基体、纤维细观本构与宏观本构建模;
其次,进行宏观损伤判断并输出损伤点的信息;
第三,调用Python文件对损伤点进行细观建模分析;
第四,利用细观建模分析的数据,计算出该点等效刚度并返回宏观模型;
具体步骤如下:
首先,进行基体、纤维等细观本构开发与建模;先通过对碳纤维增强复合材料的实验,采集并计算出碳纤维增强复合材料的在应力单元体中的X,Y,Z方向上的拉伸(压缩)极限应力和XY,XZ,YZ平面上的极限剪应力即(其中纤维方向与σx即应力单元体中X方向所受应力同向),以及X,Y,Z方向上的杨氏模量Ex,Ey,Ez、XY,XZ,YZ平面上的剪切模量、泊松比Gxy,Gxz,Gyzxyxzyz。利用上述数据和失效准则建立失效模型。通过更新的方式,把每一个积分点的相应数据代入失效模型进行失效判断。进而更新刚度矩阵,并代入应力更新函数。具体流程见说明书附图。
其次,进行宏观损伤判断并输出损伤点的信息;当碳纤维增强复合材料板受到钻削力作用时,通过Hashin失效准则来进行失效判断。若失效,则通过Abaqus的子程序vumat将失效点的信息输出到txt格式文件中。
第三,调用Python文件对损伤点进行细观建模分析;通过编写相应的Python程序,先读取上一步骤中输出的txt格式文件,获得损伤点的信息。然后,对损伤点进行再次建模,加大布种密度,将损伤点模型网格细化。而后,把读取到的损伤点信息以边界条件的形式施加到损伤点模型上。最后,启动分析。获得损伤点模型的数据,并输出到txt格式文件中。
第四,利用细观建模分析的数据,计算出该点等效刚度并返回宏观模型;由损伤点模型获得的应力数据,输入刚度退化模型中。刚度折减由以下公式求得:
其中λ为参数,通常取0.6。Ef、Em分别为纤维和基体的弹性模量,Vf、Vm分别为纤维和基体的体积含量。Df为碳纤维失效概率,可由以下公式求得:
其中σ0为对应测试跨距下尺度参数,σ为损伤点模型获得的应力数据,m为纤维的形状参数。σ0、m由实验测得。
通过上述步骤完成对损伤点的刚度折减。而后,将损伤点折减后的刚度返回到宏观模型中。在宏观模型的钻削过程中,同时进行细观与宏观分析并删除折减刚度为零的点,得到分析结果。
本发明的有益效果是:
本发明充分的考虑了纤维和基体的区别,以及仿真的精度与效率。给出了碳纤维增强复合材料钻削多尺度数值仿真方法。
本发明具有精度高,速度快的优点。
附图说明
图1是本发明的基体、纤维细观本构与宏观本构的建模流程。
图2是本发明的单层板钻削多尺度数值仿真结果图。
图3是本发明的单层板钻削多尺度数值仿真方法与普通仿真方法时间对比。
具体实施方式
下面结合实例对本发明进行进一步说明。
一种面向碳纤维增强复合材料钻削的多尺度数值仿真方法,它包括以下步骤:
首先,进行基体、纤维等细观本构与宏观本构建模。以T700碳纤维增强复合材料单层单向(0°)板为例。通过实验获得其纵向拉伸极限强度1830.7MPa,纵向极限压缩强度895.16MPa,横向拉伸极限31.31MPa,横向压缩极限124.54MPa,纵向杨氏模量99.22GPa,横向杨氏模量7.05GPa。将材料分为无限小应力单元体,根据经典材料力学理论计算出X,Y,Z方向上的拉伸(压缩)极限应力和XY,XZ,YZ平面上的极限剪应力。引入Hashin失效准则:
当纤维受拉时:
当RFT值大于1时,纤维拉伸失效。
当纤维受压时:
当RFC值大于1时,纤维压缩失效。
当基体受拉时:
当RMT大于1时,基体拉伸失效
当基体受压时:
当RMC大于1时,基体压缩失效。
通过更新的方式,把每一个积分点的相应数据代入上述准则,进行失效判断。进而更新刚度矩阵,并代入应力更新函数。具体流程见附图1。
其次,进行宏观损伤判断并输出损伤点的信息;当碳纤维增强复合材料板受到钻削力作用时,通过Hashin失效准则来进行失效判断。若失效,则通过Abaqus的子程序vumat将失效点的信息输出到txt格式文件中。
第三,调用Python文件对损伤点进行细观建模分析;通过编写相应的Python程序,先读取上一步骤中输出的txt格式文件,获得损伤点的信息。然后,对损伤点进行再次建模,加大布种密度,将损伤点模型网格细化。而后,把读取到的损伤点信息以边界条件的形式施加到损伤点模型上。最后,启动分析。获得损伤点模型的数据,并输出到txt格式文件中。
第四,利用细观建模分析的数据和T700碳纤维的尺度参数3.03GPa,形状参数6.54m,计算出该点等效刚度并返回宏观模型;由损伤点模型获得的应力数据,输入刚度退化模型中。刚度折减由以下公式求得:
计算碳纤维失效概率:
计算折减刚度:
通过上述步骤完成对损伤点的刚度折减。而后,将损伤点折减后的刚度返回到宏观模型中。在宏观模型的钻削过程中,同时进行细观与宏观分析并删除折减刚度为零的点,得到分析结果。
计算结果如图2所示。仿真时间对比如图3所示。
本发明未涉及部分与现有技术相同或可采用现有技术加以实现。

Claims (2)

1.一种面向碳纤维增强复合材料钻削的多尺度数值仿真方法,其特征是:
首先,对基体和纤维进行细观本构与宏观本构建模;
其次,进行宏观损伤判断并输出损伤点的信息;
第三,调用Python文件对损伤点进行细观建模分析;
第四,利用细观建模分析的数据,计算出损伤点等效刚度并返回宏观模型。
2.根据权利要求1所述的方法,其特征是它包括以下具体步骤:
首先,进行基体、纤维细观本构建模;先通过对碳纤维增强复合材料的实验,采集并计算出碳纤维增强复合材料的在应力单元体中的X,Y,Z方向上的拉伸或(压缩)极限应力和XY,XZ,YZ平面上的极限剪应力即以及X,Y,Z方向上的杨氏模量Ex,Ey,Ez、XY,XZ,YZ平面上的剪切模量、泊松比Gxy,Gxz,Gyzxyxzyz;其中纤维方向与σx即应力单元体中X方向所受应力同向,利用上述数据和失效准则建立失效模型;通过更新的方式,把每一个积分点的相应数据代入失效模型进行失效判断;进而更新刚度矩阵,并代入应力更新函数;
其次,进行宏观损伤判断并输出损伤点的信息;当碳纤维增强复合材料板受到钻削力作用时,通过上一步骤中准则来进行失效判断;若失效,则通过Abaqus的子程序vumat将失效点的信息输出到txt格式文件中;
第三,调用Python文件对损伤点进行细观建模分析;通过相应的Python程序,先读取上一步骤中输出的txt格式文件,获得损伤点的信息;然后,对损伤点进行再次建模,加大布种密度,将损伤点模型网格细化;而后,把读取到的损伤点信息以边界条件的形式施加到损伤点模型上;最后,启动分析,获得损伤点模型的数据,并输出到txt格式文件中;
第四,利用细观建模分析的数据,计算出损伤点等效刚度并返回宏观模型;由损伤点模型获得的应力数据,输入刚度退化模型中;刚度折减由以下公式求得:
<mrow> <msubsup> <mi>E</mi> <mn>1</mn> <mo>&amp;prime;</mo> </msubsup> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>E</mi> <mi>f</mi> </msub> <msub> <mi>V</mi> <mi>f</mi> </msub> <mo>+</mo> <msub> <mi>E</mi> <mi>m</mi> </msub> <msub> <mi>V</mi> <mi>m</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;lsqb;</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <msub> <mi>&amp;lambda;E</mi> <mi>f</mi> </msub> <msub> <mi>V</mi> <mi>f</mi> </msub> </mrow> <mrow> <msub> <mi>&amp;lambda;E</mi> <mi>f</mi> </msub> <msub> <mi>V</mi> <mi>f</mi> </msub> <mo>+</mo> <msub> <mi>E</mi> <mi>m</mi> </msub> <msub> <mi>V</mi> <mi>m</mi> </msub> </mrow> </mfrac> <mo>&amp;rsqb;</mo> <msub> <mi>D</mi> <mi>f</mi> </msub> </mrow>
其中λ为参数,通常取0.6。Ef、Em分别为纤维和基体的弹性模量,Vf、Vm分别为纤维和基体的体积含量;Df为碳纤维失效概率,由以下公式求得:
<mrow> <msub> <mi>D</mi> <mi>f</mi> </msub> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>E</mi> <mi>f</mi> </msub> <mi>&amp;sigma;</mi> </mrow> <mrow> <msub> <mi>E</mi> <mn>1</mn> </msub> <msub> <mi>&amp;sigma;</mi> <mn>0</mn> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mi>m</mi> </msup> <mo>)</mo> </mrow> </mrow>
其中σ0为对应测试跨距下尺度参数,σ为损伤点模型获得的应力数据,m为纤维的形状参数。σ0、m由实验测得;
通过上述步骤完成对损伤点的刚度折减;而后,将损伤点折减后的刚度返回到宏观模型中,在宏观模型的钻削过程中,同时进行细观与宏观分析并删除折减刚度为零的点,得到分析结果。
CN201710554685.8A 2017-07-10 2017-07-10 面向碳纤维增强复合材料钻削的多尺度数值仿真方法 Active CN107239642B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710554685.8A CN107239642B (zh) 2017-07-10 2017-07-10 面向碳纤维增强复合材料钻削的多尺度数值仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710554685.8A CN107239642B (zh) 2017-07-10 2017-07-10 面向碳纤维增强复合材料钻削的多尺度数值仿真方法

Publications (2)

Publication Number Publication Date
CN107239642A true CN107239642A (zh) 2017-10-10
CN107239642B CN107239642B (zh) 2020-08-25

Family

ID=59990949

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710554685.8A Active CN107239642B (zh) 2017-07-10 2017-07-10 面向碳纤维增强复合材料钻削的多尺度数值仿真方法

Country Status (1)

Country Link
CN (1) CN107239642B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108388763A (zh) * 2018-03-15 2018-08-10 北京航空航天大学 一种基于多尺度分析的定向纤维增强复合材料层合结构可靠性计算方法
CN109241650A (zh) * 2018-09-25 2019-01-18 南京航空航天大学 基于跨尺度仿真的碳纤维增强复合材料力学性能预测方法
CN113268897A (zh) * 2021-01-05 2021-08-17 北京航空航天大学 一种基于随机孔隙重构的多尺度建模与仿真方法
CN113420376A (zh) * 2021-06-17 2021-09-21 吉林大学 一种基于多尺度的碳纤维复合材料抗冲击力学性能仿真方法
CN113609688A (zh) * 2021-08-09 2021-11-05 大连理工大学 复合材料细观切削仿真中纤维剪断、弯断失效的精确判定方法
CN114741744A (zh) * 2022-05-13 2022-07-12 北京航空航天大学 一种适用于针刺复合材料的细观建模及多尺度分析方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100227153A1 (en) * 2008-12-22 2010-09-09 Florida State University Research Foundation Composite materials and methods for selective placement of nano-particulates within composites
CN105183990A (zh) * 2015-09-07 2015-12-23 哈尔滨工业大学 一种复合材料及其结构的多尺度响应分析方法
CN106066913A (zh) * 2016-05-31 2016-11-02 西北工业大学 复杂复合材料结构等效材料性能多尺度计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100227153A1 (en) * 2008-12-22 2010-09-09 Florida State University Research Foundation Composite materials and methods for selective placement of nano-particulates within composites
CN105183990A (zh) * 2015-09-07 2015-12-23 哈尔滨工业大学 一种复合材料及其结构的多尺度响应分析方法
CN106066913A (zh) * 2016-05-31 2016-11-02 西北工业大学 复杂复合材料结构等效材料性能多尺度计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
X.-Y. ZHOU 等: "Stochastic multi-scale finite element based reliability analysis for laminated composite structures", 《APPLIED MATHEMATICAL MODELLING》 *
李海斌 等: "复合材料宏细观结合的随机渐进失效分析", 《兵器材料科学与工程》 *
齐振超 等: "基于三维多相有限元的CFRP细观切削机理研究", 《机械工程学报》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108388763A (zh) * 2018-03-15 2018-08-10 北京航空航天大学 一种基于多尺度分析的定向纤维增强复合材料层合结构可靠性计算方法
CN109241650A (zh) * 2018-09-25 2019-01-18 南京航空航天大学 基于跨尺度仿真的碳纤维增强复合材料力学性能预测方法
CN109241650B (zh) * 2018-09-25 2023-06-23 南京航空航天大学 基于跨尺度仿真的碳纤维增强复合材料力学性能预测方法
CN113268897A (zh) * 2021-01-05 2021-08-17 北京航空航天大学 一种基于随机孔隙重构的多尺度建模与仿真方法
CN113268897B (zh) * 2021-01-05 2022-06-14 北京航空航天大学 一种基于随机孔隙重构的多尺度建模与仿真方法
CN113420376A (zh) * 2021-06-17 2021-09-21 吉林大学 一种基于多尺度的碳纤维复合材料抗冲击力学性能仿真方法
CN113420376B (zh) * 2021-06-17 2022-07-22 吉林大学 基于多尺度的碳纤维复合材料抗冲击力学性能仿真方法
CN113609688A (zh) * 2021-08-09 2021-11-05 大连理工大学 复合材料细观切削仿真中纤维剪断、弯断失效的精确判定方法
CN113609688B (zh) * 2021-08-09 2024-03-29 大连理工大学 复合材料细观切削仿真中纤维剪断、弯断失效的精确判定方法
CN114741744A (zh) * 2022-05-13 2022-07-12 北京航空航天大学 一种适用于针刺复合材料的细观建模及多尺度分析方法
CN114741744B (zh) * 2022-05-13 2024-07-09 北京航空航天大学 一种适用于针刺复合材料的细观建模及多尺度分析方法

Also Published As

Publication number Publication date
CN107239642B (zh) 2020-08-25

Similar Documents

Publication Publication Date Title
CN107239642A (zh) 面向碳纤维增强复合材料钻削的多尺度数值仿真方法
Caputo et al. Numerical study for the structural analysis of composite laminates subjected to low velocity impact
CN108427826B (zh) 缝合复合材料冲击损伤及剩余强度全程分析方法
CN113420376B (zh) 基于多尺度的碳纤维复合材料抗冲击力学性能仿真方法
CN111554362A (zh) 一种建立hbprc动态损伤本构模型的方法
Nguyen et al. A review of explicit finite element software for composite impact analysis
Arjangpay et al. Effects of structural characteristics of a bionic dragonfly wing on its low velocity impact resistance
CN108920734A (zh) Gfrp锚杆拉伸过程细观力学分析方法
CN114139300A (zh) 一种混杂纤维复合材料夹芯板抗冲击性能评价方法
Chang et al. Investigation on the mechanical properties of topologically optimized cellular structures for sandwiched morphing skins
EP4181007A2 (en) Systems and methods for semi-discrete modeling of progressive damage and failure in composite laminate materials
Mohammadi et al. A new FE modeling procedure to investigate the effects of toughening mechanisms on the fracture toughness of laminated composites
Szabo et al. Multi-scale modelling of woven carbon fibre reinforced epoxy
Alhussein et al. Simulation of the effect laminate sequence on delamination mode-I with elastic couplings
Zhicheng et al. Finite element analysis on the steel fiber-reinforced concrete beams: A systematic review
Shi et al. Short-term axial behavior of preloaded concrete columns strengthened with fiber reinforced polymer laminate
Este et al. A new approach to simulate interface damage in brittle matrix composites
Mohamed et al. Predicting the through-thickness enhancement of z-pinned composite laminates
Prasad et al. Numerical analysis of hybrid carbon fiber composite specimen and validation of results
Dayalen et al. Investigation on the influence of fiber loading on the tensile behavior of Glass/Jute hybrid composites using ANSYS 17.1
Raza et al. Development of a predictive damage methodology for hybrid wing body aircraft structures
Sabli et al. The effect of preloaded on delamination of composite laminate subjected to low energy impact
Bondarchuk et al. Investigation into residual stresses inherited during manufacturing of carbon-epoxy composites with delamination
Chernyakin et al. Probabilistic approach at delamination propagation problem in composite materials
Moldovan et al. Analysing the Performance of Hollow Core Slabs Strengthened with CFRP

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20211227

Address after: No. 29, Qinhuai District, Qinhuai District, Nanjing, Jiangsu

Patentee after: Nanjing University of Aeronautics and Astronautics Asset Management Co.,Ltd.

Address before: No. 29, Qinhuai District, Qinhuai District, Nanjing, Jiangsu

Patentee before: Nanjing University of Aeronautics and Astronautics

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220126

Address after: Room 703-29, block C, Kechuang headquarters building, No. 320, pubin Road, Nanjing, Jiangsu 211899

Patentee after: Jiangsu Hangpu Guochuang Composite Material Research Institute Co.,Ltd.

Address before: No. 29, Qinhuai District, Qinhuai District, Nanjing, Jiangsu

Patentee before: Nanjing University of Aeronautics and Astronautics Asset Management Co.,Ltd.