CN103714216A - 基于三维微观晶须构成的复合材料力学性能评估方法 - Google Patents

基于三维微观晶须构成的复合材料力学性能评估方法 Download PDF

Info

Publication number
CN103714216A
CN103714216A CN201310752049.8A CN201310752049A CN103714216A CN 103714216 A CN103714216 A CN 103714216A CN 201310752049 A CN201310752049 A CN 201310752049A CN 103714216 A CN103714216 A CN 103714216A
Authority
CN
China
Prior art keywords
whisker
stitch
dimensional
coordinate
matrix
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
CN201310752049.8A
Other languages
English (en)
Other versions
CN103714216B (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201310752049.8A priority Critical patent/CN103714216B/zh
Publication of CN103714216A publication Critical patent/CN103714216A/zh
Application granted granted Critical
Publication of CN103714216B publication Critical patent/CN103714216B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明涉及一种基于三维微观晶须构成的复合材料力学性能评估方法,属于复合材料力学领域。本方法采用欧拉法与蒙特卡洛法相结合的方法,实现了空间随机三维晶须微观构形的对增强复合材料力学性能的预估。该发明全面表征了晶须的微观结构,充分考虑了晶须数量、针脚长度与空间随机分布状况,具有晶须单元生成效率高、预估力学参数准确以及使用简单方便等优点。不但可以为四针状晶须增强聚合物的刚度分析、应力分布规律提供定性的研究,还可以推广应用到随机分布的晶须增强复合材料有效弹性性质、细观力学分析和细观损伤失效的数值计算,相比二维方法具有十分明显优势。

Description

基于三维微观晶须构成的复合材料力学性能评估方法
技术领域
本发明涉及一种基于三维微观晶须构成的复合材料力学性能评估方法,属于复合材料力学领域。
背景技术
氧化锌晶须增强复合材料作为一类结构和功能一体化的复合材料,具有优异的隐身吸波特性和结构承载能力,在国防、航天、航空和交通等领域具有广泛的应用前景。目前在考核氧化锌增强复合材料力学性能方面,主要依赖人工制备氧化锌增强复合材料,通过如拉伸试验机、霍普金森杆等设备进行力学性能的测试。但在研究优化氧化锌增强复合材料力学性能时,势必需要考察氧化锌晶须的数量与品质对整体氧化锌增强复合材料的力学性能影响,包括静态、准静态以及高应变率下的力学特性,这时就需要花费高昂的成本与时间去制备大量的实验样品,尤其是需要大量可变参数呈现梯度变化的实验样品,这些可变参数包括晶须针脚长度、针脚间的空间角度范围以及晶须的空间分布数量等。
应用有限元仿真方法开展复合材料力学性能的研究十分必要,但是借助有限元仿真方法,实现对氧化锌增强复合材料中微观嵌入单元属性的对其力学性能影响研究却面临着困难,因为晶须数量与构形特殊,在有限元建模时,现有商用有限元软件根本无法实现大量随机且可控晶须构形的生成与建模,因此也就无法考虑晶须微观结构对材料性能的影响。氧化锌晶须结构的特殊性以及随机分布的特性,使得通过有限元的手段研究多变量随机分布的氧化锌晶须复合材料的力学性能更加复杂,要进行定量的计算更加困难,需要大量的计算成本。
目前已有的研究成果中,仅有二维随机分布氧化锌晶须方法,如:氧化锌晶须应力传递机理的数值模拟研究(北京理工大学学报,2013年5月,第33卷第5期,页码450-453)。该方法将问题简化成二维平面模型,氧化锌晶须只有三个针脚,如图1所示,即为生成的两组长短数量各不相同的二维氧化锌晶须。但是,采用二维的氧化锌晶须模型对氧化锌复合材料的力学性能进行分析,其真实性与准确性显然有待商榷,而且二维模型无法实现三维力学参数的考察,因此,亟待解决三维随机氧化锌晶须的有限元生成方法,实现对三维氧化锌增强复合材料力学性能的预估。
发明内容
本发明的目的在于针对目前无法实现三维氧化锌晶须微观建模的问题,提供一种基于三维微观晶须构成的复合材料力学性能评估方法,能够在考虑晶须微观构形的情况下,基于三维微观晶须构成的复合材料有限元模型,实现对复合材料基本力学性能的评估与优化分析。
本发明的技术方案是:首先,通过刚体动力学中的欧拉法生成空间随机分布的三维晶须针脚的空间位置坐标,然后利用晶须针脚的空间坐标,利用晶须的中心坐标与n个针脚坐标,分别建立梁单元,再采用有限元方法中的嵌入单元方法,将晶须作为嵌入单元嵌入到基体材料中,进而实现复合材料的三维有限元建模。最后对有限元模型施加合适的边界与载荷,通过有限元力学分析方法,完成复合材料的力学性能评估。
对于m个、每个具有n个针脚的三维晶须,本发明的技术方案具体包括如下步骤:
步骤1,生成空间随机分布的三维晶须针脚的空间位置坐标。
步骤1.1,给定基体的长、宽、高(基体为长方体形状)以及针脚长度范围。根据蒙特卡洛方法生成一组(n+6)个随机数,将(n+6)个随机数分别作为其中一个晶须的中心节点在定坐标系下的三维坐标值、晶须所在动坐标系的三个欧拉角以及n个针脚长度。
所述晶须所在动坐标系为以晶须中心节点位置为原点,以其中一根针脚为Z轴,依据右手定则建立的坐标系。
所述定坐标系以基体的一个顶点为原点,以基体的长、宽、高作为三个坐标轴,依据右手定则建立的坐标系。
步骤1.2,根据步骤1.1生成的针脚长度及中心节点,采用刚体运动学理论和坐标变换,将n个针脚端点的动坐标转换到定坐标系下,得出晶须的中心节点与n个针脚端点在定坐标下的坐标值。
所述针脚定坐标的转换公式如下:
r=r0+A-1r'   (1)
其中r0为晶须中心节点在定坐标系下的位置矢量;A为坐标变换矩阵,具体为:
Figure BDA0000451259560000031
ψ,θ,
Figure BDA0000451259560000032
是三个随机欧拉角;r'为晶须针脚在动坐标系的位置矢量。
步骤1.3,判断定坐标系下晶须的n个节点坐标是否超出基体长宽高的范围,若超出则重复步骤1.1,否则执行步骤1.4。
步骤1.4,为了避免生成的氧化锌晶须针脚重叠,建立一个判断准则:要求任意两个晶须中心节点之间的距离要大于这两个晶须针脚的最大长度。若满足判断准则,则执行步骤1.5;若不满足,则判断新生成的晶须与已有晶须重叠,重复步骤1.1;
步骤1.5,存储晶须中心节点与各针脚顶点在定坐标系下的坐标值。循环步骤1.1至步骤1.5,直至生成m个晶须。
步骤2,建立连接关系:将每个晶须针脚中心节点分别与其针脚顶点作为一个关系并存储。根据现有商业CAD软件中关于线的定义,将晶须的每个针脚定义为一条线,并将编辑好的文件导入到CAD软件中,再输出商业有限元软件兼容的模型文件格式。
步骤3,将步骤2中的模型文件导入到商业有限元软件中,定义个针脚对应的线为梁单元;建立基体的有限元模型,并将晶须梁单元作为嵌入单元嵌入基体单元中,并分别给基体与晶须赋予对应的材料属性定义,最终实现晶须增强复合材料的有限元分析模型。
步骤4,根据需要分析的力学性能,对步骤3得到的有限元分析模型施加对应的边界条件与载荷,进行有限元仿真计算得出结果,最终完成考虑微观三维晶须的增强复合材料的力学性能评估。
有益效果
本方法采用欧拉法与蒙特卡洛法相结合的方法,实现了空间随机三维氧化锌晶须微观构形的对氧化锌增强复合材料力学性能的预估。该发明全面表征了晶须的微观结构,充分考虑了晶须数量、针脚长度与空间随机分布状况,具有晶须单元生成效率高、预估力学参数准确以及使用简单方便等优点。
本发明方法,不但可以为四针状晶须增强聚合物的刚度分析、应力分布规律提供定性的研究,还可以推广应用到随机分布的晶须增强复合材料有效弹性性质、细观力学分析和细观损伤失效的数值计算,相比二维方法具有十分明显优势。
附图说明
图1为随机生成的二维氧化锌晶须模型,(a)为15个随机短晶须,(b)为30个随机长晶须;
图2为本发明中晶须所在位置的空间动坐标系示意图;
图3为具体实施方式中含三维随机分布晶须的氧化锌增强复合材料的有限元模型;
图4为具体实施方式中考察氧化锌增强复合材料剪切力学性能的加载方式示意图;
图5为具体实施方式中氧化锌增强复合材料等效弹性模量评估值。
具体实施方式
1.根据设计要求给出基体的长度、宽度与高度,晶须针脚的长度范围(最大长度与最小长度)、晶须针脚的数量n以及晶须的数量m。以三维氧化锌晶须为例,四个针脚,n为4。
2.根据蒙特卡洛法,生成(4+6)个随机数,其中6个随机数作为晶须中心点坐标以及晶须所在动坐标系的三个欧拉角。
3.利用生成的4个随机数乘以晶须针脚的最大长度,确定晶须4个针脚长度;利用生成的3个随机数乘以基体边长,作为动坐标系的原点,同时也是晶须的中心节点坐标。坐标系的创建方法如下:
如图3所示,粗体的部分为一个四针状氧化锌的晶须,其中坐标系Oxyz为固定坐标系,O'x'y'z'为动坐标系,与晶须固连。晶须的中心点与动坐标系的原点重合,晶须的1号针脚沿着动坐标系O'x'y'z'的z'轴,2号针脚在O'x'y'平面内的投影沿着x'轴。
4.氧化锌晶须空间位置的随机性由动坐标系的随机分布而决定,利用生成的动坐标系的3个随机欧拉角,通过欧拉法坐标变换,确定动坐标系在定坐标系的相对位置。利用投影法求得到氧化锌晶须的4根针脚在定体坐标系下的节点坐标,其中r'的分量形式如公式(3):
Figure BDA0000451259560000051
其中,r1'、r2'、r3'和r4'分别是四个针脚的矢量分量,a、b、c和d分别为四根晶须针脚的随机长度。
5.利用公式(1)和公式(2)并结合公式(3),最终得到一个晶须中心节点与4个针脚顶点在定坐标系下的坐标值。
6.判断定坐标系下晶须的5个节点坐标是否超出基体范围,若超出则重复步骤2,否则继续执行;
7.在真实的结构中,氧化锌晶须不会发生重叠的状况。因此,为了避免生成的氧化锌晶须针脚重叠:要求每个晶须中心核之间的距离要大于针脚的最大长度。判断新生成的晶须是否与已有晶须重叠,若重叠则重复步骤2,否则继续执行;
8.判断晶须数量是否满足要求的数量,若不满足,则重复步骤2;否则继续执行;
9.每个氧化锌晶须有四根针脚,也就是说,每个晶须有四个连接关系。将每根针脚两端节点的节点编号存储在一个数组中,作为一个连接关系的存储。这样,每个晶须则需要一个4×2的二维数组来存储其连接关系。
10.将每个晶须的四个连接关系视为4条直线:根据Pro/e对直线的定义,编写.ibl文件,文件具体内容格式为:
Figure BDA0000451259560000061
这里,x1、y1、z1与x2、y2、z2分别是一条直线两个端点的坐标值,也就是每个针脚两个端点的坐标值。有多少个针脚变对应多少个“begincurve”将该文件导入到Pro/e中形成几何线,之后将几何导出为.igs文件。
11.将.igs文件导入到ABAQUS有限元软件中得到晶须的几何模型,对代表晶须针脚的几何线赋予梁单元属性,并划分单元,进而形成晶须的梁单元有限元模型。
12.建立基体的有限元模型,并将晶须梁单元作为嵌入单元嵌入基体单元中,并分别给基体与晶须赋予对应的材料属性定义,最终实现氧化锌晶须增强复合材料的有限元分析模型,如图3所示。
13.为了分析氧化锌晶须的剪切性能,可对氧化锌晶须增强复合材料的有限元分析模型赋予如图4所示的加载方式。若分析拉伸模量,可赋予模型一端面固定另一端面施加拉力的边界与载荷条件。图5为氧化锌增强复合材料等效弹性模量的有限元结果与理论值及实验值关系图。

Claims (2)

1.基于三维微观晶须构成的复合材料力学性能评估方法,其特征在于:对于m个、每个具有n个针脚的三维晶须,具体包括如下步骤:
步骤1,生成空间随机分布的三维晶须针脚的空间位置坐标;
步骤1.1,给定基体的长、宽、高以及针脚长度范围;根据蒙特卡洛方法生成一组(n+6)个随机数,将(n+6)个随机数分别作为其中一个晶须的中心节点在定坐标系下的三维坐标值、晶须所在动坐标系的三个欧拉角以及n个针脚长度;
所述晶须所在动坐标系为以晶须中心节点位置为原点,以其中一根针脚为Z轴,依据右手定则建立的坐标系;
步骤1.2,根据步骤1.1生成的针脚长度及中心节点,采用刚体运动学理论和坐标变换,将n个针脚端点的动坐标转换到定坐标系下,得出晶须的中心节点与n个针脚端点在定坐标下的坐标值;
所述针脚定坐标的转换公式如下:
r=r0+A-1r'   (1)
其中r0为晶须中心节点在定坐标系下的位置矢量;A为坐标变换矩阵,具体为:
Figure FDA0000451259550000011
ψ,θ,
Figure FDA0000451259550000012
是三个随机欧拉角;r'为晶须针脚在动坐标系的位置矢量;
步骤1.3,判断定坐标系下晶须的n个节点坐标是否超出基体长宽高的范围,若超出则重复步骤1.1,否则执行步骤1.4;
步骤1.4,建立判断准则:要求任意两个晶须中心节点之间的距离要大于这两个晶须针脚的最大长度;若满足判断准则,则执行步骤1.5;若不满足,则判断新生成的晶须与已有晶须重叠,重复步骤1.1;
步骤1.5,存储晶须中心节点与各针脚顶点在定坐标系下的坐标值;循环步骤1.1至步骤1.5,直至生成m个晶须;
步骤2,建立连接关系:将每个晶须针脚中心节点分别与其针脚顶点作为一个关系并存储;根据现有商业CAD软件中关于线的定义,将晶须的每个针脚定义为一条线,并将编辑好的文件导入到CAD软件中,再输出商业有限元软件兼容的模型文件格式;
步骤3,将步骤2中的模型文件导入到商业有限元软件中,定义个针脚对应的线为梁单元;建立基体的有限元模型,并将晶须梁单元作为嵌入单元嵌入基体单元中,并分别给基体与晶须赋予对应的材料属性定义,实现晶须增强复合材料的有限元分析模型;
步骤4,根据需要分析的力学性能,对步骤3得到的有限元分析模型施加对应的边界条件与载荷,进行有限元仿真计算得出结果,最终完成考虑微观三维晶须的增强复合材料的力学性能评估。
2.根据权利要求1所述的基于三维微观晶须构成的复合材料力学性能评估方法,其特征在于:所述定坐标系以基体的一个顶点为原点,以基体的长、宽、高作为三个坐标轴,依据右手定则建立的坐标系。
CN201310752049.8A 2013-12-31 2013-12-31 基于三维微观晶须构成的复合材料力学性能评估方法 Expired - Fee Related CN103714216B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310752049.8A CN103714216B (zh) 2013-12-31 2013-12-31 基于三维微观晶须构成的复合材料力学性能评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310752049.8A CN103714216B (zh) 2013-12-31 2013-12-31 基于三维微观晶须构成的复合材料力学性能评估方法

Publications (2)

Publication Number Publication Date
CN103714216A true CN103714216A (zh) 2014-04-09
CN103714216B CN103714216B (zh) 2016-06-29

Family

ID=50407187

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310752049.8A Expired - Fee Related CN103714216B (zh) 2013-12-31 2013-12-31 基于三维微观晶须构成的复合材料力学性能评估方法

Country Status (1)

Country Link
CN (1) CN103714216B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104573370A (zh) * 2015-01-16 2015-04-29 上海大学 一种基于蒙特卡洛模拟碳/碳复合材料织构的方法
CN106934109A (zh) * 2017-02-20 2017-07-07 南京航空航天大学 一种Ni3Al基合金代表体积单元模型的构建方法
CN108629147A (zh) * 2018-05-02 2018-10-09 西安电子科技大学 一种多晶体几何建模方法
CN109323830A (zh) * 2018-10-26 2019-02-12 珠海罗西尼表业有限公司 游丝刚度的确认方法
CN109507021A (zh) * 2018-09-30 2019-03-22 天津大学 一种复合材料力学性能的快速表征方法
CN109766669A (zh) * 2019-03-06 2019-05-17 四川大学 预测导电复合材料电阻及其响应的可视化数学模型方法
CN109829213A (zh) * 2019-01-21 2019-05-31 哈尔滨工业大学 一种非连续增强金属基复合材料力学性能设计与预测的方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107730434B (zh) * 2017-11-09 2020-12-01 哈尔滨工业大学 一种由两相复合材料细观数值模型生成点云模型的方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7050951B1 (en) * 2000-08-17 2006-05-23 Sandia Corporation Automatic detection of sweep-meshable volumes
CN101148780A (zh) * 2007-07-27 2008-03-26 北京理工大学 一种制备螺旋状氮化硅陶瓷晶须的方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7050951B1 (en) * 2000-08-17 2006-05-23 Sandia Corporation Automatic detection of sweep-meshable volumes
CN101148780A (zh) * 2007-07-27 2008-03-26 北京理工大学 一种制备螺旋状氮化硅陶瓷晶须的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
荣吉利等: "氧化锌晶须应力传递机理的数值模拟研究", 《北京理工大学学报》 *
邹敏等: "ZnO晶须/纳米TiO2复合粒子改善PP-R树脂性能的研究", 《功能材料》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104573370A (zh) * 2015-01-16 2015-04-29 上海大学 一种基于蒙特卡洛模拟碳/碳复合材料织构的方法
CN106934109A (zh) * 2017-02-20 2017-07-07 南京航空航天大学 一种Ni3Al基合金代表体积单元模型的构建方法
CN106934109B (zh) * 2017-02-20 2019-10-11 南京航空航天大学 一种Ni3Al基合金代表体积单元模型的构建方法
CN108629147A (zh) * 2018-05-02 2018-10-09 西安电子科技大学 一种多晶体几何建模方法
CN108629147B (zh) * 2018-05-02 2021-07-09 西安电子科技大学 一种多晶体几何建模方法
CN109507021A (zh) * 2018-09-30 2019-03-22 天津大学 一种复合材料力学性能的快速表征方法
CN109507021B (zh) * 2018-09-30 2021-03-30 天津大学 一种复合材料力学性能的快速表征方法
CN109323830A (zh) * 2018-10-26 2019-02-12 珠海罗西尼表业有限公司 游丝刚度的确认方法
CN109323830B (zh) * 2018-10-26 2020-05-15 珠海罗西尼表业有限公司 游丝刚度的确认方法
CN109829213A (zh) * 2019-01-21 2019-05-31 哈尔滨工业大学 一种非连续增强金属基复合材料力学性能设计与预测的方法
CN109766669A (zh) * 2019-03-06 2019-05-17 四川大学 预测导电复合材料电阻及其响应的可视化数学模型方法
CN109766669B (zh) * 2019-03-06 2022-09-27 四川大学 预测导电复合材料电阻及其响应的可视化数学模型方法

Also Published As

Publication number Publication date
CN103714216B (zh) 2016-06-29

Similar Documents

Publication Publication Date Title
CN103714216B (zh) 基于三维微观晶须构成的复合材料力学性能评估方法
Xie et al. The generation of arbitrary order curved meshes for 3D finite element analysis
Haga et al. A high-order unifying discontinuous formulation for the Navier-Stokes equations on 3D mixed grids
Borden et al. A phase-field description of dynamic brittle fracture
CN104038295B (zh) 一种基于机电耦合的变形阵列天线散射性能分析方法
Ayhan Simulation of three-dimensional fatigue crack propagation using enriched finite elements
Malekzadeh et al. Buckling analysis of functionally graded arbitrary straight-sided quadrilateral plates on elastic foundations
Sevilla et al. The use of hybrid meshes to improve the efficiency of a discontinuous Galerkin method for the solution of Maxwell’s equations
Munikrishna et al. Turbulent flow computations on a hybrid cartesian point distribution using meshless solver LSFD-U
Mustafa et al. Probabilistic micromechanical analysis of composite material stiffness properties for a wind turbine blade
Heo et al. Free vibration analysis of cracked plates using peridynamics
Mayeur et al. Numerical implementation of a crystal plasticity model with dislocation transport for high strain rate applications
Misra et al. Variational methods for continuum models of granular materials
Patni et al. Efficient modelling of beam-like structures with general non-prismatic, curved geometry
CN103984869A (zh) 一种预测复合材料热弹性有效属性和局部场的方法
Safdar et al. Stochastic compressive failure surface modelling for the unidirectional fibre reinforced composites under plainstress
Koohkan et al. A completely meshless analysis of cracks in isotropic functionally graded materials
Carrington et al. Three-dimensional ALE-FEM method for fluid flow in domains with moving boundaries part 1: algorithm description
Yang et al. Vertex-centered, high-order schemes for turbulent flows
Guan A moment quadrature method for uncertainty quantification of three-dimensional crack propagation via extremely few model runs
Rossikhin et al. Vibrations of Suspension Bridges: Fractional Derivative Model
Fotia et al. Quality index and improvement of the interfaces of general hybrid grids
Rossikhin et al. Vibrations of viscoelastic rods, Fractional operator models
Chen et al. Construction of n-sided polygonal spline element using area coordinates and B-net method
Genuit et al. An Implicit Discontinuous Galerkin Chimera Method for Unsteady Laminar Flow Problems with Multiple Bodies

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C53 Correction of patent of invention or patent application
CB03 Change of inventor or designer information

Inventor after: Rong Jili

Inventor after: Wang Dan

Inventor after: Xiang Dalin

Inventor after: Xu Tianfu

Inventor before: Rong Jili

Inventor before: Xiang Dalin

Inventor before: Wang Dan

Inventor before: Xu Tianfu

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: RONG JILI XIANG DALIN WANG DAN XU TIANFU TO: RONG JILI WANG DAN XIANG DALIN XU TIANFU

C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160629

Termination date: 20161231

CF01 Termination of patent right due to non-payment of annual fee