CN109829213A - 一种非连续增强金属基复合材料力学性能设计与预测的方法 - Google Patents
一种非连续增强金属基复合材料力学性能设计与预测的方法 Download PDFInfo
- Publication number
- CN109829213A CN109829213A CN201910054500.6A CN201910054500A CN109829213A CN 109829213 A CN109829213 A CN 109829213A CN 201910054500 A CN201910054500 A CN 201910054500A CN 109829213 A CN109829213 A CN 109829213A
- Authority
- CN
- China
- Prior art keywords
- reinforcement
- volume
- model
- polyhedron
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 239000002131 composite material Substances 0.000 title claims abstract description 43
- 239000002184 metal Substances 0.000 title claims abstract description 13
- 238000013461 design Methods 0.000 title abstract description 11
- 230000002787 reinforcement Effects 0.000 claims abstract description 115
- 239000002245 particle Substances 0.000 claims abstract description 52
- 239000011159 matrix material Substances 0.000 claims abstract description 39
- 230000008569 process Effects 0.000 claims abstract description 21
- 239000000463 material Substances 0.000 claims abstract description 20
- 238000009826 distribution Methods 0.000 claims abstract description 16
- 239000013078 crystal Substances 0.000 claims abstract description 9
- 238000001514 detection method Methods 0.000 claims description 15
- 239000011156 metal matrix composite Substances 0.000 claims description 12
- 230000006399 behavior Effects 0.000 claims description 8
- 238000004088 simulation Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000003014 reinforcing effect Effects 0.000 claims description 4
- 238000005482 strain hardening Methods 0.000 claims description 4
- 239000010953 base metal Substances 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000009864 tensile test Methods 0.000 claims description 3
- 239000008187 granular material Substances 0.000 claims description 2
- 238000012821 model calculation Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 8
- -1 stub Substances 0.000 abstract 1
- 238000011160 research Methods 0.000 description 6
- 239000000956 alloy Substances 0.000 description 3
- 229910045601 alloy Inorganic materials 0.000 description 3
- 239000011208 reinforced composite material Substances 0.000 description 3
- 230000005476 size effect Effects 0.000 description 3
- PCTMTFRHKVHKIS-BMFZQQSSSA-N (1s,3r,4e,6e,8e,10e,12e,14e,16e,18s,19r,20r,21s,25r,27r,30r,31r,33s,35r,37s,38r)-3-[(2r,3s,4s,5s,6r)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-19,25,27,30,31,33,35,37-octahydroxy-18,20,21-trimethyl-23-oxo-22,39-dioxabicyclo[33.3.1]nonatriaconta-4,6,8,10 Chemical compound C1C=C2C[C@@H](OS(O)(=O)=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2.O[C@H]1[C@@H](N)[C@H](O)[C@@H](C)O[C@H]1O[C@H]1/C=C/C=C/C=C/C=C/C=C/C=C/C=C/[C@H](C)[C@@H](O)[C@@H](C)[C@H](C)OC(=O)C[C@H](O)C[C@H](O)CC[C@@H](O)[C@H](O)C[C@H](O)C[C@](O)(C[C@H](O)[C@H]2C(O)=O)O[C@H]2C1 PCTMTFRHKVHKIS-BMFZQQSSSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000011165 3D composite Substances 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013329 compounding Methods 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Body Structure For Vehicles (AREA)
Abstract
一种非连续增强金属基复合材料力学性能设计与预测方法，本发明涉及一种新材料设计领域的技术，具体是非连续增强金属基复合材料中增强体形状、含量、尺寸、取向的设计、复合材料结构建模及其变形行为、损伤行为和力学性能的预测方法。其操作流程包括：基于颗粒、短棒、晶须状增强体随机分布状态，构建非连续增强复合材料的三维几何模型；对三维模型进行网格划分；通过公式计算来修正金属基体和增强体的强度；将各组分的力学性能赋予模型；对模型施加边界条件及载荷；通过仿真技术计算复合材料的力学性能。该方法具有操作简便、适用复合材料体系广、精度高等特点。
Description
技术领域
本发明涉及一种新材料设计领域的技术，具体是一种基于增强体形状、含量、尺寸、取向的非连续增强金属基复合材料结构生成及其力学性能、变形行为、损伤行为的预测方法。
背景技术
近些年，随着工程上对结构材料要求的提高，具有良好综合性能的金属基复合材料已成为国际上研究热点之一。非连续增强金属基复合材料是由增强体与金属基体复合而成，其中增强体的形状、含量、尺寸和取向设计及其与力学性能的相关性，是研发高性能金属基复合材料的关键。由于非连续增强金属基复合材料的变形过程中，塑性变形首先发生于增强体周围的应力集中区，导致基体发生不均匀的塑性变形。因此外加应力小于基体合金的屈服强度时，复合材料的应力～应变曲线就偏离线弹性区，这将给复合材料的屈服强度的预测增加难度。此外，增强体的类型、形状、含量、尺寸和取向分布状态等因素，也会影响复合材料的变形及损伤行为。因此传统金属与合金的设计方法已经不适用于非连续增强金属基复合材料的设计与性能预测。近几年有限元数值仿真技术越来越多地应用于复合材料研究中，可以较准确地预测复合材料增强体特征与力学性能的关系，已经成为了一种复合材料组分设计的重要方法。开发一种可广泛适用于构建多种微观结构的复合材料有限元建模技术，将极大地推动复合材料优化设计，满足高性能复合材料设计研发的需求。
发明内容
本发明针对现有三维复合材料建模技术上的不足，提出一种非连续增强复合材料模型构建成及其力学性能、变形行为及损伤过程的预测方法。本发明充分考虑增强体的形状、含量、尺寸、取向、分布等结构特征，生成三维代表性体积单元模型。对生成的三维代表性体积单元模型施加拉伸边界条件及载荷，通过有限元模拟仿真技术，预测复合材料的力学性能、变形行为及损伤行为。所采取的技术方案为：
一种非连续增强金属基复合材料力学性能设计与预测的方法，所述方法包括：
步骤一、选取增强体类型、尺寸及体积分数；所述增强体类型包括：颗粒、短棒和晶须；其中，颗粒可以选择为球形、六面体、二十面体或二十六面体，增强体的尺寸为颗粒直径；短棒、晶须为二十六面体，增强体的尺寸为短棒、晶须的端面直径和长径比。本发明支持多种尺寸、多种形状的增强体混合使用。
步骤二、基于增强体的分布状态构建复合材料的三维代表性体积单元模型；其中，颗粒、短棒和晶须的分布状态分别为：颗粒为随机分布；短棒和晶须为完全随机分布和定向分布；
步骤三、对复合材料的三维模型进行网格划分；
步骤四、通过模型计算来修正金属基体和增强体的强度，获得修正后的材料属性；将修正后的材料属性赋予复合材料有限元模型；其中，金属基体强度采用Taylor等效塑性应变模型进行强度修正；增强体的断裂强度采用Griffith断裂模型进行修正；
步骤五、为有限元模型施加边界条件，选取固定的O点(0，0，0)和加载点RF(100，100，100)，对面x＝0上的所有节点约束施加约束使之与O点具有相同的x坐标，即节点始终保持在x＝0平面上；面y＝0上的所有节点始终在y＝0平面上；面z＝0上的所有节点始终在z＝0平面上；对面x＝100上的所有节点约束施加约束使之与RF点具有相同的x坐标；面y＝100上的所有节点与RF点具有相同的y坐标；面z＝100上的所有节点与RF点具有相同的z坐标；对加载点RF施加x方向的位移载荷，最大载荷为4～8；即4～8％的应变量，载荷的应变速率为4×10-4s-1；
步骤六：对复合材料的拉伸试验过程进行数值仿真，获得复合材料在室温下的应力应变曲线并确定材料的变形、断裂行为。
进一步地，步骤一中所述颗粒的形状为球体或多面体，并且增强体的体积分数为0％～30％；其中，所述多面体包括六面体、二十面体和二十六面体；
进一步地，步骤一中采用颗粒直径表示颗粒尺寸；采用短棒、晶须的端面直径和长径比表示短棒、晶须尺寸。
进一步地，步骤二所述随机分布的具体建模过程包括：
第一步、计算单个增强体的体积：其中，颗粒的体积采用球体体积公式计算单个增强体的体积，短棒、晶须的体积采用圆柱体体积公式计算单个增强体的体积；若增强体形状为多面体，则对多面体进行等比例缩放使其体积等于Vparticle/Vwhisker，其中，Vparticle和Vwhisker分别表示为球体体积和圆柱体体积；
第二步、随机旋转所述增强体，分别绕Y轴Z轴旋转，旋转角度分别为θrotate和旋转操作由旋转矩阵计算表示：
第三步、将增强体移动到随机位置(xrandom,yrandom,zrandom)；
第四步、进行两个增强体的碰撞检测，避免增强体相交；
第五步、检测增强体的总体积是否达到体积分数，若满足则建模过程结束，否则再次添加增强体。
进一步地，第一步中：
所述球体体积公式为：
其中，dparticle表示颗粒的直径；
短棒、晶须的体积采用圆柱体体积公式计算单个增强体的体积，所述圆柱体体积公式包括：
其中，dsection表示短棒、晶须的截面直径；a表示短棒、晶须的长度；
进一步地，第四步所述两个增强体的碰撞检测的碰撞条件为：
若两个增强体均为球体，则碰撞检测需满足以下条件：
dparticle 2>＝(x1-x2)2+(y1-y2)2+(z1-z2)2
若两个增强体为球体和多面体，则多面体的顶点不可在球体内，碰撞检测需满足以下条件：
其中，xvex，yvex和zvex表示多面体的顶点坐标；xsphere，ysphere和zsphere表示球体的球心坐标；
若两个增强体均为多面体，则其中一个多面体的顶点不可在另一个多面体内；为方便判断顶点与多面体的关系，将多面体被划分为若干4面体，若顶点不在所有四面体内则定点不在多面体内；顶点(点O)是否在4面体(体ABCD)内需以下判据：
若D0、D1、D2、D3、D4同正或同负则顶点在四面体内，否则点在四面体外。
进一步地，步骤二所述定向分布的具体建模过程包括：
第1步、计算单个增强体的体积：
颗粒的体积采用球体体积公式计算单个增强体的体积，所述球体体积公式包括：
其中，dparticle表示颗粒的直径；
短棒、晶须的体积采用圆柱体体积公式计算单个增强体的体积，所述圆柱体体积公式包括：
其中，dsection表示短棒、晶须的直径；a表示短棒、晶须的长度；
若增强体形状为多面体，则对多面体进行等比例缩放使其体积等于Vparticle/Vwhisker；
第2步、随机旋转所述增强体，其中，所述短棒、晶须平行于向量则定向旋转矩阵具体为：
其中，θ为为晶须沿Y轴旋转的角度；
第3步、将增强体移动到随机位置(xrandom,yrandom,zrandom)；
第4步、进行两个增强体的碰撞检测，避免增强体相交；
第5步、检测增强体的总体积是否达到体积分数，若满足则建模过程结束，否则再次添加增强体。
进一步地，步骤二中所述三维代表性体积单元模型的总尺寸为100×100×100μm。
进一步地，步骤三中所述对复合材料的三维模型进行网格划分的过程为：采用有限元软件组装增强体；基体为一个实体，然后对所述实体进行三维网格的划分三维；网格类型为4节点4面体网格，划分的网格总数为800000～850000个。
进一步地，步骤四中：
所述Taylor等效塑性应变模型表示为：
其中，σ表示为基体的流变应力；σm为基体强度，αm为基体热膨胀系数，b为柏氏矢量，d为增强体尺寸，c为增强体体积分数，Em为基体弹性模量，α＊为与基体金属相关的系数，σm0为基体屈服强度，n为加工硬化系数；
所述Griffith断裂模型表示为：
其中，表示为断裂强度；K为增强体材料系数。两种理论都与增强体的尺寸有关，可用与考虑增强体的尺寸效应。
本发明有益效果：
本发明提出的一种非连续增强金属基复合材料力学性能设计与预测的方法可用于复合材料相关的细观力学研究，是材料科学研究领域中重要研究方向，对现实中的复合材料体系设计有借鉴意义。本发明相比于当前建模仿真技术，有以下几个优点：
1)增强体的空间位置通过随机数来指定，这种建模方式能够使细观结构更加贴近于实际材料。
2)增强体类型多样，模型能够有更广泛的适用性。
3)本发明中的建模方式能够建立随机分布及定向分布的短棒、晶须模型。
4)本发明提供了多种形状增强体混合使用的碰撞算法，能够更好的消除增强体重叠的现象。
5)本发明充分考虑了复合材料的增强体尺寸效应。
6)本发明构建的模型尺寸较大，数值模拟结果更加的接近于实际测试结果。
附图说明
图1为本发明具体实施方案的工艺流程图。
图2为本发明的模型构建流程图。
图3为实施例一建立的挤压态10％、20％二十六面体SiCp/6061Al三维几何模型。
图4为实施例一的模拟应力应变曲线。
图5为实施例二建立的10％二十面体、二十六面体SiCp/6061Al，5:1SiCrod/6061Al和10:1SiCw/6061Al三维几何模型。
图6为实施例二的模拟应力应变曲线。
图7为实施例二的基体等效塑性应变云图。
具体实施方式
下面结合具体实施例对本发明做进一步说明，但本发明不受实施例的限制。
实施例1：
一种非连续增强金属基复合材料力学性能设计与预测的方法，所述方法包括：
步骤一、选取增强体形状类型SiC颗粒，形状为二十六面体，颗粒的平均直径为10μm，体积分数分别选取10％及20％。
步骤二、基于增强体的分布状态构建复合材料的三维代表性体积单元模型。颗粒的分布状态为完全随机分布。
步骤三、采用有限元软件ABAQUS组装增强体、基体为一个实体，并对该实体划分三维网格，网格类型为4节点4面体网格，划分的网格总数约为825000个。
步骤四、通过公式计算来修正金属基体和增强体的强度。其中金属基体强度经Taylor等效塑性应变理论进行强度修正，基体热膨胀系数为23.6×10-61/K，伯格斯矢量长度为基体弹性模量为70GPa，α＊为3/7，基体屈服强度为270MPa，加工硬化指数为0.7。SiC颗粒的断裂强度经Griffith断裂理论修正后为2GPa；
步骤五、为有限元模型施加边界条件，选取固定的O点(0,0,0)和加载点RF(100,100,100)，对面x＝0上的所有节点约束施加约束使之与O点具有相同的x坐标，即节点始终保持在x＝0平面上；面y＝0上的所有节点始终在y＝0平面上；面z＝0上的所有节点始终在z＝0平面上。对面x＝100上的所有节点约束施加约束使之与RF点具有相同的x坐标；面y＝100上的所有节点与RF点具有相同的y坐标；面z＝100上的所有节点与RF点具有相同的z坐标。对加载点RF施加x方向的位移载荷，最大载荷为6％的应变量，载荷的应变速率为4×10- 4s-1。
步骤六：对复合材料的拉伸试验过程进行数值仿真，获得复合材料在室温下的应力应变曲线。
其中，本实施例中步骤一颗粒形状选择三种：二十面体颗粒，二十六面体短棒和二十六面体晶须。其中，颗粒的平均直径为10μm，短棒的端面平均直径为5.85μm、长径比为5:1，短棒的端面平均直径为4.64μm、长径比为10:1。增强体的体积分数均为10％。步骤五最大载荷为4％的应变量。
另外，步骤二所述随机分布的具体建模过程包括：
第一步、计算单个增强体的体积：其中，颗粒的体积采用球体体积公式计算单个增强体的体积，短棒、晶须的体积采用圆柱体体积公式计算单个增强体的体积；若增强体形状为多面体，则对多面体进行等比例缩放使其体积等于Vparticle/Vwhisker，其中，Vparticle和Vwhisker分别表示为球体体积和圆柱体体积；
第二步、随机旋转所述增强体，分别绕Y轴Z轴旋转，旋转角度分别为θrotate和旋转操作由旋转矩阵计算表示：
第三步、将增强体移动到随机位置(xrandom,yrandom,zrandom)；
第四步、进行两个增强体的碰撞检测，避免增强体相交；
第五步、检测增强体的总体积是否达到体积分数，若满足则建模过程结束，否则再次添加增强体。
第一步中：
所述球体体积公式为：
其中，dparticle表示颗粒的体积；
短棒、晶须的体积采用圆柱体体积公式计算单个增强体的体积，所述圆柱体体积公式包括：
其中，dsection表示短棒、晶须的截面直径；a表示短棒、晶须的长度；
第四步所述两个增强体的碰撞检测的碰撞条件为：
若两个增强体均为球体，则碰撞检测需满足以下条件：
dparticle 2>＝(x1-x2)2+(y1-y2)2+(z1-z2)2
若两个增强体为球体和多面体，则多面体的顶点不可在球体内，碰撞检测需满足以下条件：
其中，xvex，yvex和zvex表示多面体的顶点坐标；xsphere，ysphere和zsphere表示球体的球心坐标；
若两个增强体均为多面体，则其中一个多面体的顶点不可在另一个多面体内；为方便判断顶点与多面体的关系，将多面体被划分为若干4面体，若顶点不在所有四面体内则定点不在多面体内；顶点(点O)是否在4面体(体ABCD)内需以下判据：
若D0、D1、D2、D3、D4同正或同负则顶点在四面体内，否则点在四面体外。
步骤二所述定向分布的具体建模过程包括：
第1步、计算单个增强体的体积：
颗粒的体积采用球体体积公式计算单个增强体的体积，所述球体体积公式包括：
其中，dparticle表示颗粒的直径；
短棒、晶须的体积采用圆柱体体积公式计算单个增强体的体积，所述圆柱体体积公式包括：
其中，dsection表示短棒、晶须的截面直径；a表示短棒、晶须的长度；
若增强体形状为多面体，则对多面体进行等比例缩放使其体积等于Vparticle/Vwhisker；
第2步、随机旋转所述增强体，其中，所述短棒、晶须平行于向量则定向旋转矩阵具体为：
其中，θ为晶须沿Y轴旋转的角度；
第3步、将增强体移动到随机位置(xrandom,yrandom,zrandom)；
第4步、进行两个增强体的碰撞检测，避免增强体相交；
第5步、检测增强体的总体积是否达到体积分数，若满足则建模过程结束，否则再次添加增强体。
步骤二中所述三维代表性体积单元模型的总尺寸为100μm×100μm×100μm。
步骤三中所述对复合材料的三维模型进行网格划分的过程为：采用有限元软件组装增强体；基体为一个实体，然后对所述实体进行三维网格的划分三维；网格类型为4节点4面体网格，划分的网格总数为800000～850000个。
步骤四中：
所述Taylor等效塑性应变模型表示为：
其中，σ表示为基体的流变应力；σm为基体强度，αm为基体热膨胀系数，b为柏氏矢量，d为增强体尺寸，c为增强体体积分数，Em为基体弹性模量，α＊为与基体金属相关的系数，σm0为基体屈服强度，n为加工硬化系数；
所述Griffith断裂模型表示为：
其中，表示为断裂强度；K为增强体材料系数。两种理论都与增强体的尺寸有关，可用与考虑增强体的尺寸效应。
采用以下方式验证本发明的有益效果：
方式一：一种非连续增强金属基复合材料力学性能设计与预测的方法，通过比对10％、20％SiCp/6061Al复合材料的模拟应力应变曲线及实验数据(图4)，可以看出数值结果与实验结果使吻合的。该方法获得的仿真结果稳定可靠。
方式二：一种非连续增强金属基复合材料力学性能设计与预测的方法，通过比对10％SiCp、SiCrod、SiCw/6061Al复合材料的模拟应力应变曲线及实验数据(图6)，可以看出小变形范围内，强度变化为SiCw/6061Al>SiCrod/6061Al>SiCp/6061Al。但是SiCw/6061Al复合材料的曲线下降最早，这意味着复合材料的延伸率变化趋势为SiCw/6061Al<SiCrod/6061Al<SiCp/6061Al。从等效塑性应变云图(图7，应变载荷为4％)可以观察到，增强体开裂的个数为SiCw/6061Al>SiCrod/6061Al>SiCp/6061Al，这是由于载荷更多传递到短棒、晶须上，造成增强体承担的载荷过早达到断裂强度。而塑性应变分布显示，SiCp/6061Al的平均应变值比SiCrod/6061Al和SiCw/6061Al，由于短棒、晶须增强复合材料在小变形范围内的强度更大，可以预见其塑性应变也应该更大。但是由于增强体开裂过早造成了应力释放，因此在载荷为4％时，短棒、晶须增强复合材料的合金基体中的应变反而小。
虽然本发明已以较佳的实施例公开如上，但其并非用以限定本发明，任何熟悉此技术的人，在不脱离本发明的精神和范围内，都可以做各种改动和修饰，因此本发明的保护范围应该以权利要求书所界定的为准。
Claims (10)
1.一种非连续增强金属基复合材料力学性能设计与预测的方法，其特征在于，所述方法包括：
步骤一、选取增强体类型、尺寸及体积分数；所述增强体类型包括：颗粒、短棒和晶须；
步骤二、基于增强体的分布状态构建复合材料的三维代表性体积单元模型；其中，颗粒、短棒和晶须的分布状态分别为：颗粒为随机分布；短棒和晶须为完全随机分布和定向分布；
步骤三、对复合材料的三维模型进行网格划分；
步骤四、通过模型计算来修正金属基体和增强体的强度，获得修正后的材料属性；将修正后的材料属性赋予复合材料有限元模型；其中，金属基体强度采用Taylor等效塑性应变模型进行强度修正；增强体的断裂强度采用Griffith断裂模型进行修正；
步骤五、为有限元模型施加边界条件，选取固定的O点(0，0，0)和加载点RF(100，100，100)，对面x＝0上的所有节点约束施加约束使之与O点具有相同的x坐标，即节点始终保持在x＝0平面上；面y＝0上的所有节点始终在y＝0平面上；面z＝0上的所有节点始终在z＝0平面上；对面x＝100上的所有节点约束施加约束使之与RF点具有相同的x坐标；面y＝100上的所有节点与RF点具有相同的y坐标；面z＝100上的所有节点与RF点具有相同的z坐标；对加载点RF施加x方向的位移载荷，最大载荷为4～8；即4～8％的应变量，载荷的应变速率为4×10-4s-1；
步骤六：对复合材料的拉伸试验过程进行数值仿真，获得复合材料在室温下的应力应变曲线并确定材料的变形、断裂行为。
2.根据权利要求1所述方法，其特征在于，步骤一中所述颗粒的形状为球体或多面体，并且增强体的体积分数为0％～30％；其中，所述多面体包括六面体、二十面体和二十六面体。
3.根据权利要求1所述方法，其特征在于，步骤一中采用颗粒直径表示颗粒尺寸；采用短棒、晶须的端面直径和长径比表示短棒、晶须尺寸。
4.根据权利要求1或2所述方法，其特征在于，步骤二所述随机分布的具体建模过程包括：
第一步、计算单个增强体的体积：其中，颗粒的体积采用球体体积公式计算单个增强体的体积，短棒、晶须的体积采用圆柱体体积公式计算单个增强体的体积；若增强体形状为多面体，则对多面体进行等比例缩放使其体积等于Vparticle/Vwhisker，其中，Vparticle和Vwhisker分别表示为球体体积和圆柱体体积；
第二步、随机旋转所述增强体，分别绕Y轴Z轴旋转，旋转角度分别为θrotate和旋转操作由旋转矩阵计算表示：
第三步、将增强体移动到随机位置(xrandom,yrandom,zrandom)；
第四步、进行两个增强体的碰撞检测，避免增强体相交；
第五步、检测增强体的总体积是否达到体积分数，若满足则建模过程结束，否则再次添加增强体。
5.根据权利要求4所述方法，其特征在于，第一步中：
所述球体体积公式为：
其中，dparticle表示颗粒的直径；
短棒、晶须的体积采用圆柱体体积公式计算单个增强体的体积，所述圆柱体体积公式包括：
其中，dsection表示短棒、晶须的截面直径；a表示短棒、晶须的长度。
6.根据权利要求4所述方法，其特征在于，第四步所述两个增强体的碰撞检测的碰撞条件为：若两个增强体均为球体，则碰撞检测需满足以下条件：
dparticle 2>＝(x1-x2)2+(y1-y2)2+(z1-z2)2
若两个增强体为球体和多面体，则多面体的顶点不可在球体内，碰撞检测需满足以下条件：
其中，xvex，yvex和zvex表示多面体的顶点坐标；xsphere，ysphere和zsphere表示球体的球心坐标；
若两个增强体均为多面体，则其中一个多面体的顶点不可在另一个多面体内；为方便判断顶点与多面体的关系，将多面体被划分为若干4面体，若顶点不在所有四面体内则定点不在多面体内；顶点是否在4面体内需以下判据：
若D0、D1、D2、D3、D4同正或同负则顶点在四面体内，否则点在四面体外。
7.根据权利要求1所述方法，其特征在于，步骤二所述定向分布的具体建模过程包括：
第1步、计算单个增强体的体积：
颗粒的体积采用球体体积公式计算单个增强体的体积，所述球体体积公式包括：
其中，dparticle表示颗粒的直径；
短棒、晶须的体积采用圆柱体体积公式计算单个增强体的体积，所述圆柱体体积公式包括：
其中，dsection表示短棒、晶须的截面直径；a表示短棒、晶须的长度；
若增强体形状为多面体，则对多面体进行等比例缩放使其体积等于Vparticle/Vwhisker；
第2步、随机旋转所述增强体，其中，所述短棒、晶须平行于向量则定向旋转矩阵具体为：
其中，θ为晶须沿Y轴旋转的角度；
第3步、将增强体移动到随机位置(xrandom,yrandom,zrandom)；
第4步、进行两个增强体的碰撞检测，避免增强体相交；
第5步、检测增强体的总体积是否达到体积分数，若满足则建模过程结束，否则再次添加增强体。
8.根据权利要求1所述方法，其特征在于，步骤二中所述三维代表性体积单元模型的总尺寸为100μm×100μm×100μm。
9.根据权利要求1所述方法，其特征在于，步骤三中所述对复合材料的三维模型进行网格划分的过程为：采用有限元软件组装增强体；基体为一个实体，然后对所述实体进行三维网格的划分三维；网格类型为4节点4面体网格，划分的网格总数为800000～850000个。
10.根据权利要求1所述方法，其特征在于，步骤四中：
所述Taylor等效塑性应变模型表示为：
其中，σ表示为基体的流变应力；σm为基体强度，αm为基体热膨胀系数，b为柏氏矢量，d为增强体尺寸，c为增强体体积分数，Em为基体弹性模量，α＊为与基体金属相关的系数，σm0为基体屈服强度，n为加工硬化系数；
所述Griffith断裂模型表示为：
其中，表示为断裂强度；K为增强体材料系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910054500.6A CN109829213B (zh) | 2019-01-21 | 2019-01-21 | 一种非连续增强金属基复合材料力学性能设计与预测的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910054500.6A CN109829213B (zh) | 2019-01-21 | 2019-01-21 | 一种非连续增强金属基复合材料力学性能设计与预测的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109829213A true CN109829213A (zh) | 2019-05-31 |
| CN109829213B CN109829213B (zh) | 2021-04-30 |
Family
ID=66860410
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910054500.6A Active CN109829213B (zh) | 2019-01-21 | 2019-01-21 | 一种非连续增强金属基复合材料力学性能设计与预测的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109829213B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111931397A (zh) * | 2020-07-01 | 2020-11-13 | 中国科学院金属研究所 | 网状增强金属基复合材料的数值模拟方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103604729A (zh) * | 2013-11-19 | 2014-02-26 | 西安电子科技大学 | 一种颗粒随机分布的复合材料宏观有效性质的预测方法 |
| CN103714216A (zh) * | 2013-12-31 | 2014-04-09 | 北京理工大学 | 基于三维微观晶须构成的复合材料力学性能评估方法 |
| US20150057988A1 (en) * | 2013-08-26 | 2015-02-26 | The University Of Akron | Method for Identifying Stochastic Information of Heterogeneous Materials |
| CN104899393A (zh) * | 2015-06-19 | 2015-09-09 | 山东大学 | 一种离散相增强复合材料细观结构的建模方法 |
| CN108090963A (zh) * | 2017-12-25 | 2018-05-29 | 大连理工大学 | 一种纤维增强复合材料在低温下热残余应力的数值计算方法 |
-
2019
- 2019-01-21 CN CN201910054500.6A patent/CN109829213B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150057988A1 (en) * | 2013-08-26 | 2015-02-26 | The University Of Akron | Method for Identifying Stochastic Information of Heterogeneous Materials |
| CN103604729A (zh) * | 2013-11-19 | 2014-02-26 | 西安电子科技大学 | 一种颗粒随机分布的复合材料宏观有效性质的预测方法 |
| CN103714216A (zh) * | 2013-12-31 | 2014-04-09 | 北京理工大学 | 基于三维微观晶须构成的复合材料力学性能评估方法 |
| CN104899393A (zh) * | 2015-06-19 | 2015-09-09 | 山东大学 | 一种离散相增强复合材料细观结构的建模方法 |
| CN108090963A (zh) * | 2017-12-25 | 2018-05-29 | 大连理工大学 | 一种纤维增强复合材料在低温下热残余应力的数值计算方法 |
Non-Patent Citations (2)
| Title |
|---|
| 赵明明: "颗粒增强镁基复合材料力学性能的有限元数值分析", 《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》 * |
| 霍金星等: "ABAQUS二次开发在Csf/Mg复合材料微观几何建模中的应用", 《特种铸造及有色合金》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111931397A (zh) * | 2020-07-01 | 2020-11-13 | 中国科学院金属研究所 | 网状增强金属基复合材料的数值模拟方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109829213B (zh) | 2021-04-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Li et al. | Architecture design of periodic truss-lattice cells for additive manufacturing | |
| Song et al. | Octet-truss cellular materials for improved mechanical properties and specific energy absorption | |
| CN114462124B (zh) | 一种混凝土三维多相细观模型的建立与数值模拟方法 | |
| WO2019010859A1 (zh) | 一种高致密度离散颗粒多相体系的建模方法 | |
| Zheng et al. | Mesoscopic models and numerical simulations of the temperature field and hydration degree in early-age concrete | |
| Wang et al. | Compression behavior of strut-reinforced hierarchical lattice—Experiment and simulation | |
| WO2021139287A1 (zh) | 负泊松比结构的非线性本构关系分析方法、系统及装置 | |
| Cao et al. | Prediction on dynamic modulus of asphalt concrete with random aggregate modeling methods and virtual physics engine | |
| JP2019114265A (ja) | トポロジー最適化による部品の設計 | |
| CN110162849B (zh) | 一种混杂纤维混凝土的建模方法 | |
| CN106650018B (zh) | 一种用于vcfem分析的高体积分数rve模型生成方法 | |
| CN114186434B (zh) | 基于逐级划分的沥青混合料三维细观结构模型构建方法 | |
| CN114636360B (zh) | 五模冲击隐身复合点阵环状结构及其参数优化方法 | |
| Jun et al. | Experimental and numerical simulation of bird-strike performance of lattice-material-infilled curved plate | |
| CN112507581A (zh) | 可避免畸形计算单元的岩石细观数值模型生成方法及系统 | |
| CN109829213B (zh) | 一种非连续增强金属基复合材料力学性能设计与预测的方法 | |
| Zheng et al. | Mesoscale models and uniaxial tensile numerical simulations of concrete considering material heterogeneity and spatial correlation | |
| CN113139320A (zh) | 构建颗粒增强复合材料三维微观构型的方法 | |
| Peng et al. | 3D multiscale tri-level finite element analysis of aluminum matrix composites with Nano&Micro hybrid inclusions | |
| Bahmani et al. | Three-dimensional micromechanical assessment of bio-inspired composites with non-uniformly dispersed inclusions | |
| Xu et al. | 3D mesostructure generation of fully-graded concrete based on hierarchical point cloud and aggregate coarsening | |
| Murotani et al. | Adaptive finite elements using hierarchical mesh and its application to crack propagation analysis | |
| Chao et al. | Monte carlo based strategy to simulate the microstructure evolution of the short-fiber reinforced metal matrix composites | |
| Chen et al. | Multiscale topology optimization of gradient lattice structure based on volume parametric modeling | |
| Takeoka et al. | Non-parametric design of free-form shells with curved boundaries and specified reaction forces |
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 |