2019, Vol. 40
缩短产品的生产周期是工业领域提升市场竞争力的主要考虑因素之一, 而且关注的焦点已经从传统加工技术转变为快速生产制造技术, 如3D打印技术[1].3D打印技术是一种新兴发展的通过将材料层层累加堆积成型的制造技术.其中, 熔融沉积成型(fused deposition modeling, FDM)技术具有成本低廉、原材料范围广、环境污染小及后处理简单等优势, 是目前应用最为广泛的3D打印技术之一.FDM成型原理简单, 利用热熔喷头将材料丝加热至熔融状态, 喷头按照计算机辅助设计(CAD)的模型轮廓和填充轨迹进行运动, 与此同时将熔融材料挤出, 以与之前挤出的材料相粘结并凝固, 最终逐层累加从而堆积成实体的模型或零件[2].
由于材料逐层累加的制造工艺所导致的气孔、夹杂、无层间压力等缺陷, 致使FDM 3D打印零件的机械性能很难与传统加工工件相媲美, 严重阻碍了FDM技术的发展.近几年, 学者们的研究重点集中在FDM产品的材料力学性能方面[3-7], 但是随着3D打印技术的发展, 其产品越来越多地应用于机械加工、交通运输等实际工况中, 仅仅进行材料力学特性分析已难以准确判定零件的可靠性, 因此需要对制品进行更为实际的振动特性分析, 然而目前却很少有学者进行此方面的研究.Arivazhagan[2]利用DMA 2980设备对FDM 3D打印样件进行了0~100 Hz动力学特性参数的扫频试验, 进而分析了温度与储能模量和黏度值的关系.Mohamed等[8-9]采用同样的方法实验研究了不同过程参数对样件动态弹性模量的影响规律, 并指出动态弹性模量随着打印间隙、光栅角、打印方向和路径宽度的增加而减小.此外, 通过过程参数的优化, 确定了最大动态模量和玻璃过渡温度的数值.然而, 上述研究都是基于实验测试开展的, 缺乏理论基础，并且只考虑了动态载荷下的弹性模量等参数的变化情况, 仍属于材料力学性能的研究范畴, 而并非振动特性方面的研究.
本文首次针对不同打印方向的FDM 3D打印薄板的振动特性进行了理论建模和实验研究，通过将固有特性的理论结果与实验结果相对比, 验证了FDM 3D打印薄板固有特性理论模型的正确性.
1 FDM薄板固有特性理论分析 1.1 理论建模本文由六层PLA纤维材料组合而成的3D打印薄板具有正交各向异性特点，因此在研究3D打印薄板的理论模型时, 首先作了如下假设:①层与层之间粘结牢固, 无滑移, 无相对位移, 不考虑层间耦合效应; ②薄层之间粘结良好, 可作为一个整体结构板, 并且粘结层很薄, 其本身不发生变形, 即各单层板之间变形连续; ③薄板虽由多层单层板叠合而成, 但其总厚度仍符合薄板假定.
采用双向梁函数组合法[10]对3D打印薄板进行动力学建模与分析, 薄板的振动解析是运用梁函数与另一方向满足边界条件的乘积型振型函数对建立的振动基本微分方程进行求解, 在满足边界条件情况下, 列出振型方程, 获得固有频率, 从而得到振型.
3D打印薄板的动力学建模可根据薄板理论假设(即Kirchhoff假设), 采用薄板小挠度理论研究.以板变形前的中面为xoy平面, 建立空间固定的直角坐标系o-xyz, 长度方向为x, 宽度方向为y, 厚度方向为z.矩形薄板的几何模型如图 1所示, 其中, 长度a为130 mm, 宽度b为50 mm, 厚度h为2.4 mm(每层厚度均为0.4 mm).
|
图 1 3D打印薄板三维模型 Fig.1 3D model of the 3D printed sheet |
薄板自由振动的基本微分方程为
|
(1) |
式中:ρ是密度, kg/m3; D是板的弯曲刚度, N·m, 
本文的FDM 3D打印薄板边界条件为固定与自由组合情况下的矩形板, 采用双向梁函数组合法对其进行振动特性分析, 从而描述薄板样件振动的真实振型.将3D打印薄板的振型设为
|
(2) |
式中:aij是待定系数, 其作用是调整不同阶次梁函数的组合以实现振型函数表达式逼近FDM 3D打印薄板的真实振型的目的；Xi(x)是满足x方向边界条件的第i阶振型函数表达式; Yj(y)是满足y方向边界条件的第j阶振型函数表达式.
梁振型函数的一般表达式为
|
(3) |
沿x轴方向为自由边界条件, Xi(x)可取自由-自由单向板第n阶振型函数, 需要满足的边界条件:
当x=0和x=a时，
|
(4) |
将式(3)代入边界条件(4)中, 可得边界条件为自由-自由梁的振型函数如下:
|
(5) |
式中:ki是与梁频率有关的系数, ki4=ω2(ρA/EJ)；αi为梁频率系数, αi=(chkia-coskia)/(shkia-sinkia)(i=3, 4, 5, …)；第一阶振型X1和第二阶振型X2分别表示刚体的移动和扭转, 第三阶振型Xi(i≥3)才是由梁变形引起的振型.
梁在Y方向的边界条件为固定-自由, 其表达式为
当y=0时，
|
(6) |
当y=a时,
|
(7) |
将式(3)代入边界条件(6)和(7)中, 可得边界条件为固定-自由梁的振型函数如下:
|
(8) |
梁频率系数如表 1所示.
|
表 1 梁频率系数 Table 1 Coefficient of beam frequency |
将式(5)和(8)代入式(2)中, 便可以得到FDM悬臂板的振型函数:
|
(9) |
Ritz(瑞利-里兹)法是根据能量变分原理, 将泛函极值问题转化为多元函数极值问题来近似求解工程结构或系统固有特性参数的一种方法.FDM 3D打印薄板在振动过程中会发生变形, 产生应力、应变以及随时间发生变化的板内变形能.
薄板势能为
|
(10) |
薄板动能为
|
(11) |
位移表达式为
|
(12) |
将式(12)代入式(10), (11)中, 再代入由Hamilton原理得到的变分式
|
(13) |
|
(14) |
|
(15) |
将振型函数(9)代入式(14), (15)中, 且Umax-Tmax对待定系数αij的偏导数为0, 即
|
(16) |
得到关于αij的线性代数方程式, 令方程组的系数行列式为0, 可求出圆频率ω, 根据圆频率与固有频率换算公式
本文利用D-FORCE V2型桌面级3D打印机制备了外形结构如图 2所示的试件(Z, X和45°方向分别打印), 其长度为150 mm, 测试宽度为50 mm, 厚度为2.4 mm.
|
图 2 试件的二维示意图 Fig.2 Two-dimensional drawing of the specimen |
试验试件共有三种类型:①Z方向打印试件(Z01和Z02);②X方向打印试件(X01和X02);③ 45°方向打印试件(S01和S02).需要注意的是, 除了打印方向不同外, 其他如打印层厚度、喷嘴直径、打印速度、路径宽度等加工过程参数的设置都是相同的.
根据ASTM D638标准, 利用拉伸试验机完成了试件的静力拉伸实验, 进而确定FDM薄板的弹性模量和剪切模量.通过电子分析天平, 分析获取了样件的密度.样件的材料属性如表 2所示，泊松比为0.3.
|
|
表 2 FDM 3D打印薄板样件材料属性 Table 2 Material properties of the FDM 3D printing plate |
为获取FDM试件的固有特性参数, 搭建了如图 3所示的测试系统.测试系统主要包括模态力锤(PCB 086C01)、数据采集卡(NI USB 4431)和加速度传感器(B&K 4517)等设备.其中, 选用的B&K 4517轻质加速度传感器仅重0.6 g, 因此可忽略其对薄板振动特性的影响.
|
图 3 测试系统示意图 Fig.3 Schematics of the test system |
利用夹具固定保持FDM薄板试件的悬臂状态, 其夹持区长度为20 mm, 激励点距离底边约10 mm, 加速度传感器分别固定于薄板振动响应较大的测点处(样件顶部和中部, 避开节点), 以获取准确的测试结果.
实验过程中, 首先用力锤对不同方向打印制成的试件(Z方向、X方向和45°方向)施加脉冲激励, 同时利用加速度传感器测试由脉冲激励所引起的FDM薄板的动态响应, 并利用数据采集卡对输入的激励信号和输出的响应信号进行实时采集.为获得薄板的模态振型, 采用单点拾振法完成对试件的参数识别, 即将加速度传感器固定在振动较大的测点处, 分别对其余各测点进行激励.针对每个试件, 分别完成了10组实验, 以进一步确保测试结果的准确性.
图 4为一组不同打印方向的FDM薄板试件的频率响应曲线, 通过分析辨识该曲线的峰值数据即可初步确定薄板试件的各阶固有频率, 具体如表 3所示.
|
图 4 FDM薄板的频率响应函数 Fig.4 Frequency response function of the FDM sheet |
|
|
表 3 FDM薄板样件的固有频率 Table 3 Natural frequency of the FDM sheet samples |
为了验证本文所提出的悬臂边界下FDM 3D打印薄板固有特性计算方法的正确性, 将理论计算结果与实验结果进行对比分析, 结果如表 4所示.其中误差为|C-1|/A.
|
|
表 4 FDM 3D打印薄板前3阶固有特性的理论和实验结果 Table 4 Theoretical and experimental results of the first three order inherent characteristics of the FDM 3D printed sheets |
通过与实验结果进行对比验证可知, 基于梁函数组合法的3D打印薄板计算振型和测试结果一致, 且固有频率的计算结果与实验结果吻合度较好, 误差仅为5 % ~9.7 %, 验证了理论模型的正确性, 即利用本文所提出的理论计算方法可以较为准确地分析和预测出悬臂边界条件下FDM 3D打印薄板的固有特性参数(固有频率和模态振型).而误差产生的原因主要来自理论建模和实验测试两个方面, 如模型忽略了各层之间的应力应变情况、薄板纤维排列不规则、残余应力、实验过程中边界条件不一致(夹具夹紧程度不同)、测试仪器灵敏度漂移等问题.
4 结论1) 基于小挠度薄板理论对FDM 3D打印薄板在悬臂状态下进行了理论建模, 计算了3D打印薄板的固有频率及振型, 该方法同样适用于同类型多层无滑移的复合薄板的固有特性分析.
2) 给出了计算3D打印薄板固有特性的具体流程.主要包括输入薄板的几何参数和不同打印方向的薄板的输入参数、基于双向梁函数组合法来表示振型函数以及通过Ritz法对该类型薄板的固有特性进行求解等3个关键步骤.
3) 以FDM 3D打印薄板为研究对象, 将固有特性的理论分析结果与实验测试结果进行了对比, 误差范围为5 % ~9.7 %, 吻合度较好, 验证了理论计算方法的正确性.
4) 本文的研究有助于提高FDM产品的动力学性能, 并且为今后的研究提供重要的参考和技术支持.
| [1] |
Panda S K, Padhee S, Sood A K, et al.
Optimization of fused deposition modeling(FDM) process parameters using bacterial foraging technique[J]. Intelligent Information Management, 2009, 1(2): 89–97.
|
| [2] |
Arivazhagan A, Masood S H.
Dynamic mechanical properties of ABS material processed by fused deposition modeling[J]. International Journal of Engineering Research and Applications, 2014, 2(3): 2009–2014.
|
| [3] |
Brensons I, Polukoshko S, Silins A, et al.
FDM prototype experimental research of processing parameter optimization to achieve higher tensile stress[J]. Solid State Phenomena, 2015, 220/221: 767–773.
DOI:10.4028/www.scientific.net/SSP.220-221 |
| [4] |
Mohamed O A, Masood S H, Bhowmik J L.
Experimental investigations of process parameters influence on rheological behavior and dynamic mechanical properties of FDM manufactured parts[J]. Advanced Manufacturing Processes, 2016, 31(15): 1983–1994.
|
| [5] |
Torres J, Cotelo J, Karl J, et al.
Mechanical property optimization of FDM PLA in shear with multiple objectives[J]. Journal of the Minerals Metals & Materials Society, 2015, 67(5): 1183–1193.
|
| [6] |
Rayegani F, Onwubolu G C.
Fused deposition modeling(FDM) process parameter prediction and optimization using group method for data handling(GMDH) and differential evolution(DE)[J]. International Journal of Advanced Manufacturing Technology, 2014, 73(1/2/3/4): 509–519.
|
| [7] |
Sood A K, Ohdar R K, Mahapatra S S.
Experimental investigation and empirical modeling of FDM process for compressive strength improvement[J]. Journal of Advanced Research, 2012, 3(1): 81–90.
DOI:10.1016/j.jare.2011.05.001 |
| [8] |
Mohamed O A, Masood S H, Bhowmik J L.
Investigation of dynamic elastic deformation of parts processed by fused deposition modeling additive manufacturing[J]. Advances in Production Engineering & Management, 2016, 13(3): 227–238.
|
| [9] |
Mohamed O A, Masood S H, Bhowmik J L, et al.
Effect of process parameters on dynamic mechanical performance of FDM PC/ABS printed parts through design of experiment[J]. Journal of Materials Engineering & Performance, 2016, 25(7): 2922–2935.
|
| [10] |
曹志远.
板壳振动理论[M]. 北京: 中国铁道出版社, 1989.
( Cao Zhi-yuan. Vibration theory of plate and shell[M]. Beijing: China Railway Publishing House, 1989. ) |