基于Mori⁃Tanaka方法的短纤维复合材料开孔板分析

doi:10.19491/j.issn.1001-9278.2020.07.010

中国塑料 ›› 2020, Vol. 34 ›› Issue (7): 61-68.DOI: 10.19491/j.issn.1001-9278.2020.07.010

基于Mori⁃Tanaka方法的短纤维复合材料开孔板分析

董治男¹^,², 徐峰祥¹^,²(), 郭巍¹^,², 华林¹^,², 苏建军³

^1.武汉理工大学现代汽车零部件技术湖北省重点实验室，武汉 430070
^2.武汉理工大学汽车零部件技术湖北省协同创新中心，武汉 430070
^3.湖北齐星汽车车身股份有限公司，湖北随州 441300

收稿日期:2019-12-02 出版日期:2020-07-26 发布日期:2020-07-26
基金资助:
国家重点研发计划(2018YFB1106700);国家自然科学基金(51975438);高等学校学科创新引智计划(B17034)

Analysis of Opening Plate Injection⁃Molded by Short⁃Fiber⁃Reinforced Thermoplastic Composites Based on Mori⁃Tanaka Method

Zhinan DONG¹^,², Fengxiang XU¹^,²(), Wei GUO¹^,², Lin HUA¹^,², Jianjun SU³

^1.Hubei Key Laboratory of Advanced Technology for Automotive Components, Wuhan University of Technology, Wuhan 430070, China
^2.Hubei Collaborative Innovation Center for Automotive Components Technology, Wuhan University of Technology, Wuhan 430070, China
^3.Hubei Qixing Automobile Body Co, Ltd, Suizhou 441300, China

Received:2019-12-02 Online:2020-07-26 Published:2020-07-26
Contact: Fengxiang XU E-mail:xufx@whut.edu.cn

摘要/Abstract

摘要：

基于Mori?Tanaka均匀化方法，对短纤维复合材料开孔板力学性能进行研究。首先，使用Mori?Tanaka方法预测单向和二维2种分布形式的短纤维复合材料应力?应变曲线曲线。其次，通过试验方法分析不同区域纤维取向对开孔板力学性能的影响。最后，使用多尺度方法研究不同注塑位置和不同开孔方式对开孔板力学性能的影响，该部分在宏观有限元分析中耦合了注塑模拟的纤维方向分布以及Mori?Tanaka均匀化方法得到的复合材料参数。结果表明，短边注塑开孔板综合力学性能更好；孔周围纤维分布是影响应力集中的主要因素；整体的纤维取向分布影响开孔板的刚度；单向纤维材料的角度对材料性能的影响不是单调的。

关键词: 开孔板, 多尺度, 短纤维, 方向, 应力集中

Abstract:

This paper reported an investigation on the mechanical properties of opening plate injection?molded by short?fiber?reinforced thermoplastic composites by means of the Mori?Tanaka homogenization method. The Mori?Tanaka method was adopted to predict the stress?strain curves of the composites in both unidirectional and two?dimensional distributions, and then the effect of fiber orientation on the mechanical properties of the opening plate were analyzed by experimental methods. Finally, a multi?scale method was used to study the effects of different injection positions and different opening methods on the mechanical properties of the opening plate. This study was coupled with the fiber direction distribution of the injection molding simulation and the material parameters obtained by the Mori?Tanaka homogenization method in the macroscopic finite element analysis. The results indicated that the opening plate molded by short?side injection exhibited better comprehensive mechanical properties, and the fiber distribution around the hole was the main factor affecting the stress concentration. The overall fiber orientation distribution influenced the stiffness of the opening plate, and however there is a complicated influence from the angle of unidirectional fiber?reinforced composites on the material properties.

Key words: shell with a circle hole, multi?scale, short fiber, orientation, stress concentration

中图分类号:

TQ 320.66⁺2

董治男, 徐峰祥, 郭巍, 华林, 苏建军. 基于Mori⁃Tanaka方法的短纤维复合材料开孔板分析[J]. 中国塑料, 2020, 34(7): 61-68.

Zhinan DONG, Fengxiang XU, Wei GUO, Lin HUA, Jianjun SU. Analysis of Opening Plate Injection⁃Molded by Short⁃Fiber⁃Reinforced Thermoplastic Composites Based on Mori⁃Tanaka Method[J]. China Plastics, 2020, 34(7): 61-68.

图/表 22

图1 用基于角度θ和φ的坐标系定义方向矢量p

Fig.1 Axes de?nition based on the angles of θ and φ, de?ning the orientation vector of p

图2 多尺度方法流程图

Fig.2 Flow chart of multi?scale method

图3 不同区域纤维取向影响分析模型

Fig.3 Analysis model of fiber orientation influence in different regions

图4 注塑模型

Fig.4 Injection model

图5 开孔板有限元模型

Fig.5 Finite element model of the shell with a circle hole

样品	E/ MPa	泊松比	密度/g?cm^?3	屈服应力/ MPa	硬化系数	硬化模量/ MPa
PA66	3 100	0.35	1.14	25	325	150
玻璃纤维	76 000	0.22	2.54	-	-	-

表1 组分的参数

Tab.1 Properties of the composite constituents

图6 单向纤维材料应力?应变曲线

Fig.6 Stress?strain curve of unidirectional fiber material

图7 二维纤维材料应力?应变曲线

Fig.7 Stress?strain curve of 2D fiber material

试验号	区域A角度/°	区域B角度/°
1	0	0
2	15	15
3	30	30
4	45	45
5	60	60
6	75	75
7	90	90

表2 不同区域单向纤维试验

Tab.2 Unidirectional fiber tests in different regions

试验号	区域A张量第一特征值	区域B张量第一特征值
1	0.1	0.1
2	0.2	0.2
3	0.3	0.3
4	0.4	0.4
5	0.5	0.5
6	0.6	0.6
7	0.7	0.7
8	0.8	0.8
9	0.9	0.9

表3 不同区域二维纤维试验

Tab.3 Two?dimensional fiber tests in different regions

图8 角度对应力集中系数的影响

Fig.8 Effect of the angle on stress concentration factor

图9 方向张量第一特征值对应力集中系数的影响

Fig.9 Effect of the first?order components on the stress concentration factors

图10 角度对刚度的影响

Fig.10 Effect of the angle on the stiffness

图11 方向张量第一特征值对刚度的影响

Fig.11 Effect of the first eigenvalue of the directional tensor on the stiffness

图12 短边注塑方案纤维取向图

Fig.12 Fiber orientation pattern for long side injection

图13 长边注塑方案纤维取向图

Fig.13 Fiber orientation pattern for long side injection

图14 纤维方向张量的映射结果

Fig.14 Mapping results of fiber direction tensors

图15 短边注塑方案路径图

Fig.15 Short side injection path diagram

图16 长边注塑方案路径图

Fig.16 Long side injection path diagram

图17 短边注塑方案有限元结果

Fig.17 Finite element results of short side injection

图18 长边注塑方案有限元结果

Fig.18 Finite element results of long side injection

项目	应力集中系数	刚度/MPa
短边注塑机械开孔	2.203	5 166.4
短边注塑开孔板	2.00	5 160.9
长边注塑机械开孔	1.99	4 873.4
长边注塑开孔板	2.16	4 870.7

表4 应力集中系数和刚度

Tab. 4 Stress concentration factor and stiffness

参考文献 15

1	YU S Y, MAI Y W, LAUKE B. Science and Engineering of Short Fibre Reinforced Polymer Composites[M]. 2009.
2	BABU K P, MOHITE P M, UPADHYAY C. Development of an RVE and Its Stiffness Predictions Based on Mathematical Homogenization Theory for Short Fibre Composites[J]. International Journal of Solids and Structures，2018, 130⁃131(1):80⁃104.
3	SHOKRIEH M M, MOSHREFZADEH⁃SANI H. An optimized Representative Volume Element to Predict The Stiffness of Aligned Short Fiber Composites[J]. Journal of Composite Materials，2016, 50(23):3 301⁃3 310.
4	TIAN W, QI L, SU C, et al. Numerical Evaluation on Mechanical Properties of Short⁃fiber⁃reinforced Metal Matrix Composites: Two⁃step Mean⁃field Homogenization Procedure[J]. Composite Structures ，2016, 139(4):96⁃103.
5	陈玉丽,马勇,潘飞,等. 多尺度复合材料力学研究进展[J].固体力学学报,2018,39(01):1⁃68.
	CHEN Y L, MA Y, PAN F, et al. Research Progress of Multi⁃scale Composite Mechanics [J]. Journal of Solid Mechanics, 2018, 39 (1): 1⁃68.
6	BERBENNI S, CAPOLUNGO L A. Mori⁃Tanaka Homogenization Scheme for Non⁃linear Elasto⁃viscoplastic Heterogeneous Materials Based on Translated Fields: An Affine Extension[J]. Comptes Rendus Mécanique, 2015, 343(2):95⁃106.
7	GONG S, LI Z, ZHAO Y Y. An Extended Mori⁃Tanaka Model for The Elastic Moduli of Porous Materials of Finite Size[J]. Acta Materialia, 2011, 59(17):6 820⁃6 830.
8	LURIE S, SOLYAEV Y, SHRAMKO K. Comparison Between The Mori⁃Tanaka and Generalized Self⁃consistent Methods in The Framework of Anti⁃plane Strain Inclusion Problem in Strain Gradient Elasticity[J].Mechanics of Materials,2018, 122(7):133⁃144.
9	胡学川, 李又兵. 基于Moldflow的汽车中央控制面板中玻纤取向的分析[J]. 塑料工业, 2014,42(10):53⁃55.
	HU X C, LI Y B. Analysis of Glass Fiber Orientation in Automobile Central Control Panel Based on Moldflow [J]. Plastics Industry, 2014,42 (10): 53⁃55.
10	黄达勇, 丁智平, 荣继刚, 等. 短玻纤增强工程塑料纤维取向张量影响因素研究[J]. 湖南工业大学学报, 2016, 30(3):12⁃18.
	HUANG D Y, DING Z P, RONG J G, et al. Study on Influencing Factors of Orientation Tensor of Short Glass Fiber Reinforced Engineering Plastics [J]. Journal of Hunan University of Technology, 2016, 30 (3): 12⁃18
11	ADVANI, SURESH G. The Use of Tensors to Describe and Predict Fiber Orientation in Short Fiber Composites[J]. Journal of Rheology, 1987, 31(8):751⁃784.
12	GRUBER B, HUFENBACH W, KROLL L, et al. Stress Concentration Analysis of Fibre⁃reinforced Multilayered Composites with Pin⁃loaded Holes[J]. Composites Science and Technology, 2007, 67(7/8):1 439⁃1 450.
13	NGUYEN D H, WANG H. Multi⁃scale Analyses of Three Dimensional Woven Composite 3D Shell with a Cut Out Circle[J]. Applied Composite Materials, 2019, 26(1): 339⁃356.
14	PIERARD O, GONZALEZ C, SEGGURADO J, et al. Micromechanics of Elasto⁃plastic Materials Reinforced with Ellipsoidal Inclusions[J]. International Journal of Solids and Structures, 2007, 44(21):6 945⁃6 962.
15	TUCKER III C L, LIANG E. Stiffness Predictions for Unidirectional Short⁃fiber Composites: Review and Evaluation[J]. Composites Science and Technology, 1999, 59(5):655⁃671.

[1]	杨惠娣. 聚氯乙烯及其稳定剂现状与发展趋势[J]. 中国塑料, 2019, 33(4): 111-119.
[2]	黄薇. 电动车箱盖纤维填充取向与浇口位置关系分析[J]. 中国塑料, 2016, 30(09): 54-57 .
[3]	刘喜山谷晓昱赵斌许晓光孙军张胜. 硅酸盐纳米短纤维增强EVA复合材料的阻燃和力学性能[J]. 中国塑料, 2012, 26(07): 34-38 .
[4]	刘先锋洪剑城卜海涛张雄伟吴世见. 脆性材料拉伸试样结构与应力集中关系研究[J]. 中国塑料, 2009, 23(09): 80-84 .
[5]	郑津洋秦永泉施建峰李翔徐平郭伟灿王非. 聚乙烯管道热熔接头卷边对拉伸性能的影响[J]. 中国塑料, 2009, 23(07): 77-80 .
[6]	邓可;卜海涛;洪剑城;刘先锋;吴世见. 脆性材料哑铃形拉伸试样的应力分布及断裂位置研究[J]. 中国塑料, 2009, 23(02): 88-92 .
[7]	祝铁丽;盛占辉;刘永云;王敏杰. 图像处理在注塑制品模内收缩研究中的应用[J]. 中国塑料, 2008, 22(11): 52-55 .

基于Mori⁃Tanaka方法的短纤维复合材料开孔板分析

Analysis of Opening Plate Injection⁃Molded by Short⁃Fiber⁃Reinforced Thermoplastic Composites Based on Mori⁃Tanaka Method

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表 22

参考文献 15

相关文章 7

编辑推荐

Metrics

本文评价

试验号	区域A张量第一特征值	区域B张量第一特征值
1	0.1	0.1
2	0.2	0.2
3	0.3	0.3
4	0.4	0.4
5	0.5	0.5
6	0.6	0.6
7	0.7	0.7
8	0.8	0.8
9	0.9	0.9

试验号	区域A张量第一特征值	区域B张量第一特征值
1	0.1	0.1
2	0.2	0.2
3	0.3	0.3
4	0.4	0.4
5	0.5	0.5
6	0.6	0.6
7	0.7	0.7
8	0.8	0.8
9	0.9	0.9

试验号	区域A张量第一特征值	区域B张量第一特征值
1	0.1	0.1
2	0.2	0.2
3	0.3	0.3
4	0.4	0.4
5	0.5	0.5
6	0.6	0.6
7	0.7	0.7
8	0.8	0.8
9	0.9	0.9

试验号	区域A张量第一特征值	区域B张量第一特征值
1	0.1	0.1
2	0.2	0.2
3	0.3	0.3
4	0.4	0.4
5	0.5	0.5
6	0.6	0.6
7	0.7	0.7
8	0.8	0.8
9	0.9	0.9

试验号	区域A张量第一特征值	区域B张量第一特征值
1	0.1	0.1
2	0.2	0.2
3	0.3	0.3
4	0.4	0.4
5	0.5	0.5
6	0.6	0.6
7	0.7	0.7
8	0.8	0.8
9	0.9	0.9

试验号	区域A张量第一特征值	区域B张量第一特征值
1	0.1	0.1
2	0.2	0.2
3	0.3	0.3
4	0.4	0.4
5	0.5	0.5
6	0.6	0.6
7	0.7	0.7
8	0.8	0.8
9	0.9	0.9