Impact crushing experiment and numerical evaluation method for composite laminates
-
摘要: 纤维增强复合材料层板是电路板基板的重要组成部分, 预测其破碎过程是废旧电路板破碎机设计的重要依据.首先, 基于实验测得复合材料层板的冲击能量, 对比面内和面外两种冲击破碎方式.然后, 基于Hashin损伤理论建立复合材料层板冲击破碎预测模型, 预测的冲击消耗能量与冲击速度与实验结果吻合度高, 表明Hashin损伤理论可以应用于玻璃纤维复合材料层板冲击破碎问题.通过复合材料层板的面内和面外冲击损伤过程对比分析, 发现面内冲击破碎可以避免层间开裂造成的额外能耗, 所以沿平面方向的冲击破碎效果要好于沿厚度方向的破碎效果.当前冲击损伤预测模型可进一步应用于预测电路板基板的破碎过程, 也可用于指导电路板破碎机传动系统设计.Abstract: Fiber reinforced laminates are commonly used in the manufacture of circuit board substrates, impact crushing process prediction of the fiber reinforced lamitates is an important basis for the design of a crushing machine of used circuit board. Firstly, impact energy corresponding to in plane and out-of-plane impact is obtained by pendulum impact tests. Then, the Hashin damage theory as well as the cohesive zone method is used to simulate the impact crushing process of glass fiber reinforced laminates. The numerical results of impact energy and velocity changes and are in good agreement with the experiments, which indicates that the present model can be applied to the low-speed impact problem of glass fiber composite laminates. According to the impact damage process analysis, it is found that in-plane impact crushing can avoid additional energy consumption caused by interface delamination, so that the in-plane impact crushing effect is better than the out-plane impact crushing effect. The present prediction model can be further applied to predict the crushing process of circuit board substrates and design the transmission system of the impact crushing machine.
-
Keywords:
- composite laminates /
- impact /
- crushing process /
- damage criterion /
- numerical simulation[1]
-
玻璃纤维增强复合材料因其强度高、质量轻,具有减震性、抗疲劳性、耐化学品腐蚀性等优点,广泛应用于汽车工业、航空航天、电子通信等各种领域中[1].作为各工业基础产品的印制电路板主要由铜箔基板、导电图形以及电子元器件等部分组成.铜箔基板在电路板制造中使用最多的是环氧树脂玻璃纤维增强复合材料[2].全世界每年都有大量的电子产品需要报废,而对于废旧印制电路板的处理也日益受到重视.在废旧电路板中,不仅金属组分具有极高的利用价值,而作为电路板基板的玻璃纤维复合材料也能作为填料或混凝土添加剂,其颗粒作为黏弹性材料时,具有更高的温度敏感[3-4],能使混凝土更坚实耐用;其粉末也可用来制备高密度聚乙烯复合材料、活性炭、吸附剂等[5-6].因此废旧印制电路板基板的处理及资源化能够产生较高的经济、社会和环境效益.
目前大多数学者研究玻璃纤维复合材料面内冲击损伤机理是为找到增强其剩余强度的方法. Xin, et al.[7]研究了以磨细玻璃纤维为填料对玻璃/环氧复合材料层合板冲击和冲击后压缩性能的影响,研究表明磨碎玻璃纤维填充复合材料的冲击性能和冲击后压缩(CAI)性能的得到改善,承载性能进一步提高. Mars[8]等研究了玻璃纤维增强聚酰胺作为基体材料的低速冲击性能,研究了玻璃纤维体积分数和附带冲击能对复合材料圆板冲击响应的影响.研究表明冲击速度对复合材料板的冲击响应有较大影响. Kevin, et al.[9]研究了相同表面密度和材料组成的二维和三维编织玻璃/环氧复合材料的低速冲击损伤,结果表明三维编织复合材料的分层长度和开口均小于相同能量下的二维复合材料.
恰恰相反,在基于废旧电路板基板破碎回收的研究中,需要找到便于废旧电路板冲击破碎的方式,从而可以提高其回收再利用的效率.文中将利用冲击实验和复合材料损伤理论研究玻璃纤维复合材料破碎过程,对比分析面内与面外玻璃纤维复合材料层板两个方向的冲击破碎效果,确定合理的冲击破碎方式,为进一步应用于实际电路板破碎预测和破碎机器的设计提供理论依据.
1 冲击试验
本次实验采用的是钢研纳克摆锤式仪器化NI300冲击试验机,用规定高度的摆锤对处于冲击试验机砧座上的冲击试样进行一次性打击,测量试样弯折时的冲击吸收能量Ak,即冲断试样消耗的能量.
摆锤是以冲击试验机摆轴为中心,在处于预仰角150°位置释放摆锤,以一定速度冲击试验机砧座上的冲击试样,冲击过程要保证打击点与打击中心重合.两种方向不同的标准试样见图 1(a)、图 1(c),试样尺寸均为10 mm×10 mm×55 mm.在进行试验前进行几次空打试验,确定仪器的空打能损在标称能量的0.5%以内.按GB/T 229—2007的试验标准在室温下进行实验,其传感器型号为2011R2013-5v-80,摆锤的初始落脚为150°,初始速度为5.5 m/s.
按实验操作规范将两组试样分别进行15次冲击,得出冲击消耗能量如表 1所列,其中沿平面方向冲击方式称为面内冲击,沿厚度方向的冲击方式称为面外冲击,要说明的是该测试结果包含了冲击锤与试样之间的能量、试样与底座之间的摩擦功、支架轴承摩擦功、空气阻力功等外界因素;同时表 1列出了15组冲击前后的能量消耗及速度变化,进一步对比沿平面方向及厚度方向冲击破碎效果.传感器测得沿平面方向冲击过程的冲击吸收功保持在3.5~4.0 J之间,沿厚度方向冲击消耗能量值在9.3~10.3 J之间,冲击后两种情况下冲击锤速度分别为面内冲击5.454~5.463 m/s、厚度方向冲击5.404~5.414 m/s.
表 1 冲击过程的各项参数Table 1. Table of parameters of impact process冲击能量是破碎过程中重要的依据,能量的损耗值代表玻璃纤维复合材料的抗冲击性能的强度.冲击完成后的试样见图 1(b)和图 1(d).平面方向开口槽的试样冲击脆断情况如图 1(b)所示,其断口平滑,所消耗的冲击能量较小;厚度方向开口槽的试样断裂如图 1(d)所示,该试样所需的冲击能量更大,是平面方向开口槽试样的3~4倍,且断口粗糙、参差不齐,有粘连现象,即试样在冲击处未完全断裂.
采用体视显微镜拍摄试样断口形貌如图 2(a)所示,沿平面方向冲击过程中绝大部分的冲击力由纤维层承载,断口平滑且无明显分层现象;由图 2(b)可知,沿厚度方向冲击的过程中由于层间胶体的阻碍作用,冲击力沿胶体传递到纤维层,因此纤维层之间出现明显的分层现象.
2 复合材料层板材料属性的获取
2.1 测试方法
玻璃纤维增强复合材料基本参数测试实验采用MTS-E45电子万能试验机和材料测试及生物力学光测设备(Digital Image Correlation,DIC).
2.2 实验过程
拉伸试验:采用常规的板状试样进行拉伸试验时,在试样的夹持部位易产生破坏和滑移.因此,参考ASTM D 3039-08[10],采用哑铃型试样(如图 3(a)所示)开展拉伸实验,玻璃环氧纤维复合材料板拉伸速率为1 mm/s.
面内剪切试验:Iosipescu剪切法是目前公认的测试多向层合板面内剪切强度的试验方法,其原理是通过专用夹具对带有双V形槽层压板试样(如图 3(b)所示)两边施加一对大小相等、方向相反的集中载荷,在试样工作区形成均匀的面内剪切应力场[11].为准确测得玻璃纤维层板的剪切参数,设计制作了(如图 3(b)所示)的双V型剪切夹具.试样的相关尺寸设计参考ASTM D5379标准[12],拉伸速率为0.5 mm/s.在试样的夹持部位制作与夹具相对应的螺纹孔,是为了避免试样在拉伸过程中出现滑移,从而保证实验的准确性.
拉伸剪切试验:采用双切口拉伸剪切法测定复合材料层合板的层间剪切强度,是用于测定复合材料层合板层间剪切强度的有效试验方法.双切口试样的设计(如图 3(c)所示)使剪切面上尽可能得到纯剪切应力,并使应力尽可能均匀分布,同时避免了压缩剪切法试验夹具结构的烦琐对试验结果的影响[13].拉伸速率为0.5 mm/s.
基于以上拉伸实验、剪切实验的结果,进行合理的转化运算后得出玻璃纤维复合材料的相关材料属性,如表 2所列.
表 2 玻璃纤维复合材料属性Table 2. Table of material properties of glass fiber composites3 复合材料层板冲击破碎模拟
3.1 Hashin损伤准则
玻璃纤维复合材料破坏准则采用Hashin损伤准则,Hashin损伤起始准则考虑4种失效模式[14]:
纤维拉伸失效:
(1) 纤维压缩失效:
(2) 基体拉伸失效:
(3) (4) 式中:系数a代表剪切应力σ12对纤维拉伸失效判据的影响系数(设定α=0)[15].
Hashin损伤模型是基于能量的渐进损伤模型,参数d表示损伤变量[16]:
(5) 式中:δeq指在有限元分析中考虑特征有限元长度来评估的等效位移;δeq0指材料刚发生破坏时的等效位移;δeqf指材料破坏之后的等效位移;这种方法是采用应力-位移而不是用应力-应变关系来描述损伤演化规律,有效减少了模型的网格依赖性[17].
3.2 Cohesive界面脱黏模型
采用Cohesive单元模拟玻璃纤维层板的分层失效过程,采用线弹性-线性软化本构模型[18].分层之前,界面是线弹性的,应力包括3部分:
(6) 开始分层阶段,界面是线性软化阶段,即当δ0≤δ≤δmax时,为材料软化区域.此时有:
(7) 式中:tn为界面正应力;ts和tt为剪切应力;Knn、Kss、Ktt分别为相应的界面刚度;δ为位移;D是损伤系数,当D=0时,表示材料没有屈服或刚开始屈服,当D=1时,表示材料破坏,失去承载能力.
采用二次应变准则用于判定混合型载荷下的损伤产生:
(8) 损伤演化以Benzeggagh–Kenan断裂准则为依据[19],表达式如下:
(9) 式中:GC的为混合型分层扩展的临界断裂能释放率;GnC、GsC、GtC分别为界面在法向和2个剪切方向破坏所需要的断裂能量;η为界面性能参数,取值范围一般为1~2.
3.3 冲击过程模拟
将仪器化冲击试验机摆锤的冲头部分进行建模,其各项参数如表 3所列.
表 3 摆锤冲头的材料参数Table 3. Table of material parameters for pendulum punch纤维复合材料采用Hashin准则,其参数由实验测得并参考文献[20],如表 4所列.层合板层间材料参数如表 5所列[21].
表 4 纤维材料参数一览Table 4. Table of fiber material parameters表 5 Cohesive的材料参数Table 5. Table of cohesive parameters在有限元仿真中,如图 4(a)所示,将边界条件设为固定底座,同时在锤头与试样、试样与底座之间采取摩擦系数为0.2的通用接触.建立参考点1将冲头约束为刚体,同时对于刚体的运动方向进行约束,只保证其在Y轴上的自由度.
冲击锤的初速度与实验相等,设为5 500 mm/s.如图 4(b)所示,网格类型为SC8R,采用缩减积分,单元数量为43 980.采用适用于瞬时动态问题的ABAQUS/Explicit显示动力学分析方法,将冲击时间设为0.005 s.
ABAQUS通过绘制伪应变能(ALLAE)和内能(ALLIE)来评价沙漏模式对计算结果的影响[22].为验证计算结果的可靠性,输出冲击过程中的伪应变能/内能比ALLAE/ALLIE,如图 5所示,其值保持在3%以内,表明可接受有限元模拟分析结果.
4 冲击实验与仿真结果对比
4.1 对比冲击能量
冲击过程中冲击锤头的动能变化曲线如图 6所示,沿厚度方向冲击能量消耗明显大于沿平面方向冲击的能量损耗值,与实验现象一致;沿平面方向冲击时,冲击锤头动能由最初的300 J逐渐减少至296.24 J,冲击能量变化值为3.76 J,对比分析,发现仿真与实验的误差平均值为2.21%;沿厚度方向冲击试样时,摆锤锤头的动能由300 J下降至290.794 J,试样的冲击吸收功为9.206 J,对比分析误差保持在6.4%.
4.2 对比冲击速度
冲击过程中冲击锤速度的变化曲线如图 7所示.冲击开始后,冲击锤受到试样的阻碍作用,速度急剧减小;随后试样出现损伤,对其阻碍作用减弱,且试样在冲击过程中是自由状态,因此速度下降趋势变缓,直至冲击结束.冲击锤在0.005 s的时间历程中,沿平面方向冲击的摆锤速度值由最初的5.5 m/s降低为5.465 m/s,沿厚度方向冲击的锤头速度降至5.415 m/s.对比分析实验与仿真结果,数值模拟的摆锤速度变化结果与实验高度吻合.
在整个冲击过程中,根据能量守恒原理:
(10) 式中:m为冲击锤的质量;v0为冲击试样前的初速度;v1为冲击完成时的末速度;h为冲击过程中摆锤质心的高度差;W为冲击吸收功.
将2种冲击实验中所得的15组冲击吸收功取平均值,得出沿平面方向冲击的能量为3.843 J,沿厚度方向冲击的能量损耗为9.84 J.由上述公式可求出沿平面方向冲击后的末速度平均值分别为5.458 m/s,与模拟所得摆锤的速度差距0.007 m/s,计算两者误差为0.128%;沿厚度方向冲击的实验平均速度值为5.408 m/s,两者误差为0.129%.这进一步证明了当前数值模拟模型的正确性.
4.3 两种冲击破碎方式的机理分析
受到沿平面的横向冲击时,有限元仿真输出的纤维复合材料层板的各种损伤参量云图如图 8所示,由损伤云图可知,最先出现的是基体损伤,纤维损伤则随后才出现,且基体损伤程度远大于纤维损伤;在锤头接触端出现纤维与基体压缩损伤,在U型缺口处出现纤维与基体拉伸损伤;在冲击过程中纤维拉伸损伤大于纤维压缩损伤,基体拉伸损伤与压缩损伤程度相当.实验中冲击损伤主要是由于缺口处纤维和基体的共同断裂,在远离缺口试样末端并未出现明显损伤.
纤维复合材料层板在受到沿厚度方向的冲击损伤情况如图 9所示,与平面内冲击破坏不同,面外损伤区域呈层条状,此时纤维与基体几乎同时出现损伤,且损伤程度相当;在冲击过程中纤维拉伸损伤远大于纤维压缩损伤,基体拉伸损伤与压缩损伤程度相当;如图 10所示,由于各纤维层间由胶体相连接,冲击断裂过程中会伴随层间开裂,在不同纤维层间出现多条裂纹,这与实验中的分层现象相同;同时由于在冲击过程中出现的层板分层现象,会造成额外的能量消耗.
5 结论
1)通过设计不同方向的冲击试样进行对比试验,得出玻璃纤维复合材料沿平面方向的抗冲击性能弱于沿厚度方向的抗冲击性能.
2)在有限元仿真中,基于Hashin损伤理论建立复合材料层板冲击破碎预测模型, 理论预测的冲击消耗能量、冲击速度与实验结果吻合度高,表明Hashin损伤理论可以应用于玻璃纤维复合材料层板冲击破碎问题.
3)通过玻璃纤维复合材料层板的面内冲击损伤过程分析,发现纤维损伤、基体损伤相继出现;在锤头接触端出现压缩损伤,在U型缺口处出现拉伸损伤;纤维拉伸损伤远大于纤维压缩损伤,基体拉伸损伤与压缩损伤程度相当.在面外冲击损伤中,纤维与基体几乎同时出现损伤,且损伤程度相当;在冲击过程中纤维拉伸损伤远大于纤维压缩损伤,基体拉伸损伤与压缩损伤程度相当.
4)在面外冲击过程中会出现复合材料层板的分层现象,会造成额外的能量消耗,所以推荐用面内冲击的方式进行冲击破碎.
-
表 1 冲击过程的各项参数
Table 1 Table of parameters of impact process
表 2 玻璃纤维复合材料属性
Table 2 Table of material properties of glass fiber composites
表 3 摆锤冲头的材料参数
Table 3 Table of material parameters for pendulum punch
表 4 纤维材料参数一览
Table 4 Table of fiber material parameters
表 5 Cohesive的材料参数
Table 5 Table of cohesive parameters
-
[1] ZHANG Y H, LIU S L, XIE H H, et al. Current status on leaching precious metals from waste printed circuit boards[J]. Procedia Environmental Sciences, 2012, 16(4):560-568. http://cn.bing.com/academic/profile?id=0bf6d7946bdecee158cddcdce94c4674&encoded=0&v=paper_preview&mkt=zh-cn
[2] HICKS C, DIETMAR R, EUGSTER M. The recycling and disposal of electrical and electronic waste in China-legislative and market responses[J]. Environmental Impact Assessment Review, 2005, 25(5):459-471. doi: 10.1016/j.eiar.2005.04.007
[3] 赵春虎.废旧印制电路板中非金属的热解处理及金的回收技术研究[D].广州: 华南理工大学, 2017. http://kns.cnki.net/KCMS/detail/detail.aspx?dbcode=CDFD&dbname=CDFD&filename=1017734288.nh [4] DANGTUNGEE R, SOMCHUA S, SIENGCHIN S. Recycling glass fiber/epoxy resin of waste printed circuit boards:morphology and Mechanical properties[J]. Mechanics of Composite Materials, 2012, 48(3):325-330. doi: 10.1007/s11029-012-9279-1
[5] ZHOU Y H, QIU K Q. A new technology for recycling materials from waste printed circuit boards[J]. Journal of Hazardous Materials, 2010, 175(1/2/3): 823-828. http://cn.bing.com/academic/profile?id=1cff8ff4a3e453b2f75c4d160171b83c&encoded=0&v=paper_preview&mkt=zh-cn
[6] SARAVANAKUMAR K, SUBRAMANIAN H, ARUMUGAM V, et al. Influence of milled glass fillers on the impact and compression after impact behavior of glass/epoxy composite laminates[J]. Polymer Testing, 2019, 75:133-141. doi: 10.1016/j.polymertesting.2019.02.007
[7] XIN H H, LIU Y Q, MOSALLAM A S, et al. Evaluation on material behaviors of pultruded glass fiber reinforced polymer (GFRP) laminates[J]. Composite Structures, 2017, 182: 283-300. doi: 10.1016/j.compstruct.2017.09.006
[8] MARS J, CHEBBI E, WALI M, et al. Numerical and experimental investigations of low velocity impact on glass fiber-reinforced polyamide[J]. Composites Part B: Engineering, 2018, 146: 116-123. doi: 10.1016/j.compositesb.2018.04.012
[9] KEVIN R, PATRICK X L, LAWRENCE E, et al. Mechanisms and characterization of impact damage in 2D and 3D woven fiber-reinforced composites[J]. Composites Part A: Applied Science and Manufacturing, 2017, 101:432-443. doi: 10.1016/j.compositesa.2017.07.004
[10] ASTM D 3039-08. Standard test method for tensile properties of fiber reinforced metal matrix Composites[Z]. American Society for Testing and Materials, 2008.
[11] 徐琪.复合材料面内剪切性能测试方法的研究[J].玻璃纤维, 2012(3):6-10. doi: 10.3969/j.issn.1005-6262.2012.03.002 [12] ASTM D5379/D5379M-12. Standard test method for shear properties of composite materials by the V-notched beam method[Z]. American Society for Testing and Materials, 2012.
[13] 王瑞, 陈海霞, 郭兴峰, 等.层合板复合材料的层间剪切强度评价方法及其改进研究[J].玻璃钢/复合材料, 2004(3):8-11. doi: 10.3969/j.issn.1003-0999.2004.03.002 [14] 吴振, 陈健.基于Hashin准则的复合材料层合结构低速冲击研究[J].沈阳航空航天大学学报, 2017, 34(5):12-20. doi: 10.3969/j.issn.2095-1248.2017.05.002 [15] LONG S, YAO X, ZHANG X. Delamination prediction in composite laminates under low-velocity impact[J]. Composite Structures, 2015, 132:290-298. doi: 10.1016/j.compstruct.2015.05.037
[16] LINDE P, BOER H D. Modelling of inter-rivet buckling of hybrid composites[J]. Composite Structures, 2006, 73(2):221-228. doi: 10.1016/j.compstruct.2005.11.062
[17] FAGGIANI A, FALZON B G. Predicting low-velocity impact damage on a stiffened composite panel[J]. Composites Part A: Applied Science and Manufacturing, 2010, 41(6):740-749. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=c7b17252d522457ee11e215643ee8f0b
[18] DUARTE A P C, DíAZ SáEZ A, SILVESTRE N. Comparative study between XFEM and hashin damage criterion applied to failure of composites[J]. Thin-Walled Structures, 2015, 115:277-288. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=6606edeee94a09c519f6d11f47a37c69
[19] LI Z, GHOSH S, GETINET N. Micromechanical modeling and characterization of damage evolution in glass fiber epoxy matrix composites[J]. Mechanics of Materials, 2016, 99:37-52. doi: 10.1016/j.mechmat.2016.05.006
[20] 拓宏亮, 马晓平, 卢智先.基于连续介质损伤力学的复合材料层合板低速冲击损伤模型[J].复合材料学报, 2018, 35(7):202-212. http://d.old.wanfangdata.com.cn/Periodical/fhclxb201807026 [21] 陈海立.冲破式复合材料易碎盖破坏机理研究[D].南京: 南京航空航天大学, 2013. http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=D565394 [22] 庄茁.基于ABAQUS的有限元分析和应用[M].北京:清华大学出版社, 2009:258. -
期刊类型引用(7)
1. 李悦,谢贤,李加文,张守逊. 铜硫浮选分离技术与药剂研究进展. 化工矿物与加工. 2024(01): 50-59 . 百度学术
2. 李超,凌石生. 内蒙古某铜矿低碱度铜硫分离浮选试验研究. 矿冶. 2024(01): 39-45 . 百度学术
3. 周立波,曹飞,曹进成,吕良. 河南某含金银硫化铜矿选矿试验研究. 黄金科学技术. 2024(02): 377-386 . 百度学术
4. 汪国辉,何海洋,汪浩翔,蔡锦鹏,赖浩,申培伦,刘殿文. 铜铁硫化矿浮选分离捕收剂的研究进展. 中国有色金属学报. 2024(12): 4090-4110 . 百度学术
5. 田长顺,饶运章,苏港,黄涛. 磁黄铁矿对黄铁矿粉尘爆炸特性参数的影响. 有色金属科学与工程. 2023(05): 692-699 . 本站查看
6. 马妮娅,林圆圆,杨会兵. 龙桥矿业铜硫分离工艺优化研究. 现代矿业. 2022(02): 143-146+149 . 百度学术
7. 钟春晖,冯博,严华山,周鹰,刘刚,陈渊淦,宁湘菡,汪惠惠. 三种有机抑制剂在辉钼矿与滑石浮选分离中的作用. 中国有色金属学报. 2022(12): 3843-3852 . 百度学术
其他类型引用(2)