Next Article in Journal
A Study on the Improvement of Parallel Operation Characteristics of DC/DC Converter Using Improved Full-Duplex Communication Hardware
Previous Article in Journal
On Transparency and Accountability of Smart Assistants in Smart Cities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Processing Parameters on the Dynamic Characteristic of Material Extrusion Additive Manufacturing Plates

1
School of Mechanical Engineering and Automation, Northeastern University, Shenyang 110819, China
2
College of Information Science and Engineering, Northeastern University, Shenyang 110819, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(24), 5345; https://doi.org/10.3390/app9245345
Submission received: 29 October 2019 / Revised: 2 December 2019 / Accepted: 4 December 2019 / Published: 6 December 2019
(This article belongs to the Section Mechanical Engineering)

Abstract

:

Featured Application

With the development and deepening of this research, it will help improve the mechanical performance of ME components, and extend the applications of ME additive manufacturing technique.

Abstract

Material extrusion (ME), an additive manufacturing technique, can fabricate parts almost without geometrical limitations. With the growing application of ME parts, especially in actual working conditions, the dynamic characteristics are needed to be studied to accurately determine their reliability. This study provides an experimental validation of the theoretical model for predicting the dynamic characteristics of ME plates fabricated with three different key processing parameters, i.e., extrusion width, layer height and build direction. The model is set up based on the bidirectional beam function combination method, and a series of experimental tests are performed. It is found that different processing parameters result in the material properties of the samples to vary, thus leading to different dynamic characteristics. Through the comparison between predictions and measurements, it is shown that the influencing trend of the processing parameters is predicted precisely. The theoretical model gives reliable predictions in dynamic characteristics of ME plates. The natural frequency discrepancy is below 13.4%, and the predicted mode shapes are the same as the measured ones. This present work provides theoretical basis and technical support for further research in improving the dynamic performance of ME products, and helps extend the applications of this technique.

1. Introduction

Currently, increasing attention is being given to additive manufacturing techniques [1,2] compared to traditional processing techniques in order to promote competitiveness in the manufacturing market. Material extrusion (ME) technique, the term used here according to the standard of ISO/ASTM 52900:2015 to include the fused filament fabrication (FFF) or fused deposition modeling (FDM) technique, can fabricate products almost without geometrical limitations. It is an increasingly growing additive manufacturing technique in modern manufacturing industry, and the built end-use parts have been used in a wide range of applications, such as automotive, dentistry, electronic and aerospace industry [3,4]. In this technique, the raw material filament is first melted in an extrusion head. It is then extruded out from the nozzle, which moves in a scheduled mode. The molten filament then cools, solidifies and bonds with the extruded material layer beneath it. A three-dimensional model or part entity is finally completed by repeating this procedure until the last layer.
During the past few years, extensive research has been carried out in studying various quality characteristics of ME products, such as mechanical property, surface roughness and build quality Hwang [5] developed new composite filaments mixing the acrylonitrile butadiene styrene (ABS) thermoplastic with copper and iron particles. The thermo-mechanical properties of the ME samples made of these filaments were investigated. It was found that the tensile strength was decreased with the loading of metal particles, but the thermal conductivity was improved. Torres [6] studied the effects of key processing parameters (the layer height, infill density and postprocessing heat-treatment time at 100 °C) on the material properties of ME polylactic acid (PLA) parts tested in torsion. It was concluded that to improve the fabricated parts’ strength while preserving their ductility and reliability needs low levels of heat treatment. The layer height and infill density were important for optimizing strength, and ductility was significantly affected by infill and heat treatment. Tsouknidas [7] evaluated the shock mitigating properties of the polymeric structures fabricated by an ME device. The effect of layer height, infill pattern and density on the energy dissipation properties was examined. The results showed that the impact absorption capacity was influenced a lot by the infill density, whereas the effect of the other two factors was small. Li [8] studied the effects of layer height, deposition velocity and infill rate on ME parts’ mechanical properties, and summarized the quantitative relationship between these processing parameters and the parts’ tensile strength. The proposed theoretical model gave accurate predictions of the tensile strength, and it could reveal the work mechanism of these parameters. Rodríguez-Panes [9] compared the effect of layer height, infill density and build direction on the mechanical performance of PLA and ABS test samples, which included tensile yield stress, tensile strength, nominal strain at break, and modulus of elasticity. This work provided assistance in selecting proper filament materials for ME. Cuan-Urquizo et al. [10] thoroughly reviewed different experimental, analytical and computational studies for characterizing the mechanical properties of ME parts under different loading conditions. The applicability and limitations of these studies were summarized. The three approaches were recommended to be used together for accurately estimating the mechanical properties. The predictive tools, able to predict the mechanical performance, deformation and failures of ME parts, deserve further research. Durgun [11] investigated the selection of proper build directions and raster angles to achieve the optimum surface roughness and mechanical property of ME parts. The results showed that the direction has a more significant effect than the raster angle. The relationship between the surface roughness and mechanical properties was close. Kim [12] experimentally studied the effects of the gap between nozzle and substrate, filament feeding rate and printing speed on the ME parts’ surface roughness. The width, height and cross-sectional shapes of the extrudate were examined. The results showed that these factors had significant effect on the built parts’ surface roughness. Boschetto [13] developed a methodology able to improve the surface roughness of ME parts by computer numerical control machining. The proposed method was validated by a case study, achieving a great reduction of average roughness and a reliable uniformity of finished surfaces. Lalehpour [14] set up a reliable algorithm that could predict the surface roughness of ME parts with an accurate definition of the centerline taken into account. The method was then validated experimentally. With an ensemble learning algorithm considered, Li [15] proposed a data-driven predictive model for predicting the surface roughness of ME parts. Based on the condition monitoring data measured by multiple sensors of different types from an ME device, the predictive model is validated, and it was able to accurately predict the surface roughness of ME parts. Kestilä [16] combined atomic layer deposition (ALD) and ME technique to coat the aluminum oxide on the built small satellite propulsion parts. The benefits of the coating were studied in the context of in-space propulsion fluidics, with propellant flow properties taken into account. It was found that the parts’ surface roughness was improved, leading to improved flow properties. Sun [17] investigated the mechanisms controlling the bonding formation between the extruded filaments in ME process. The bonding quality was assessed in terms of measuring and analyzing the mesostructure and the degree of healing between the adjacent extrudate. The results showed that the build direction, liquefier temperature and convection coefficient had strong effects on the cooling temperature profile, as well as on the mesostructure and bonding strength of the built part. Gordeev [18] provided an efficient method to evaluate the quality of the parts built by ME. It was found that the wall permeability of the part was dependent on its geometric shape. The part’s quality was primarily influenced by the filament feeding rate, wall geometry and G-code-defined wall structure. Optimizing these parameters could improve the quality (decrease in the porosity) and sealing properties of ME parts. Based on short beam shear tests, Caminero [19] evaluated the effect of layer height and fiber content on the interlaminar bonding performance of the fiber reinforced nylon composites manufactured by ME technique. It was observed that the layer height of nylon influences the interlaminar shear performance slightly. Fiber reinforced samples showed higher interlaminar shear performance than unreinforced ones, but the level of increase was moderate with the increase of fiber content. Narahara [20] used the atmospheric pressure plasma to carry out the hydrophilic treatment to the ME printing layer. This method could significantly improve the rupture strength and surface free energy (up by 75%) of ME parts, and thereby enhance their inter-layer bonding strength. Besides, more and more attentions have been paid to ME parts’ accuracy levels and tolerances [21] to increase their industrial applications. Melenka [22] studied the dimensional accuracy of the parts built by an ME device under the same manufacturing condition. The analysis showed that the test specimens varied significantly from the nominal dimensions. Careful selection of processing parameters was suggested to optimize the dimensional accuracy. Boschetto [23] developed a design for manufacturing methodology able to improve the dimensional accuracy of ME parts. This approach, predicted the obtainable dimensional deviations in the model-design stage, provided compensations of the deviations generated in the fabrication stage. From the experiments, it was shown that the dimensional accuracy of the ME parts was improved significantly using this method. Turner [24] had done a literature review about the dimensional accuracy of ME parts. It was summarized that dimensional accuracy of the built parts was dependent upon the processing parameters and the properties of the feedstock filament properties. Particularly, the contours and layer height played key roles in determining the parts’ dimensional accuracy. Nuñez [25] analyzed the dimensional accuracy of ME parts made of ABS-plus material in order to establish the dimensional tolerance ranges. The effect of two processing parameters (the infill density and layer height) on the dimensional accuracy was experimentally studied. It was concluded that the optimum configurations for the fabrication of ME parts with ABS-plus was set up. The infill density was the most influencing factor.
With the increasingly wider application of ME products in actual working conditions, the study in terms of vibration characteristics requires urgent attention. However, very few investigations have been focused on this research area, and the built parts’ reliability can hardly be determined, especially in actual working conditions, such as machining, transportation, etc. Arivazhagan [26] used the DMA2980 equipment to conduct the frequency sweep test on ME samples from 1 to 100 Hz to determine their modulus, damping and viscosity values. The results illustrated that the strength of the part processed with the solid built style was higher than those fabricated with other built styles. Increasing temperature would increase their loss modulus, but decrease the viscosity values and storage modulus. Mohamed [27,28,29] used the same test method to investigate different processing parameters’ effect on the dynamic elastic modulus of samples, and pointed out that modulus values would decrease when the build direction, raster angle or air gap increased. In addition, the numerical values of maximum glass transition temperature and dynamic modulus were determined by optimizing the processing parameters. Domingo-Espin [30] studied the effect of processing and testing parameters on the dynamic response of ME samples. Dynamic mechanical analysis (DMA) was used to determine the storage and loss modulus under an oscillating load. The results showed that the number of contours was the most influential parameter. Processing parameters (nozzle diameter, air gap and number of contours) could control the storage modulus of the manufactured parts, and testing parameters (temperature, loading frequency, and amplitude) showed a great influence on their damping capacity. Mansour [31] performed a series of tests to investigate the dynamic behavior of ME parts made of polyethylene terephthalate glycol (PETG) and PETG reinforced with 20% carbon fibers. The results revealed that the loss factor and damping obtained from the cyclic compression and models tests dropped from 17.3% to 15.4% and 13.8% to 12.3%, respectively. Although the above studies were carried out under vibratory and cyclic conditions, the dynamics inherent characteristics of the samples were not taken into account, i.e., natural frequency and mode shape. Furthermore, the existing research mainly relies on experimental testing, and lacks theoretical basis and principle.
In this paper, a theoretical model was established and experimental research was performed to quantify the influence of different processing parameters on the dynamic inherent characteristics of ME plates. The sample preparations are introduced first, followed by the theoretical analysis and experimental tests. The samples’ material properties are then provided, and the effect of three processing parameters on the dynamic characteristics is quantified through the comparison between predictions and measurements.

2. Materials and Methods

To investigate the dynamic characteristics of ME plates theoretically and experimentally, the fabrication of the samples was first introduced in this section, followed by the theoretical modeling and experimental analysis.

2.1. Sample Preparation

To determine the single effect of the processing parameter on tensile property, three kinds of dense-structure samples were fabricated by an ME device (D-force V2) according to ISO 527-2-2012 (International standard, plastics determination of tensile properties Part 2: test conditions for molding and extrusion plastics), as the schematic shown in Figure 1a. The first kind could be divided into three types, and they had different extrusion widths (a certain width of the filament extruded from the nozzle), which were separately 0.4, 0.3 and 0.2 mm. These samples were represented by TE0.4_i, TE0.3_i and TE0.2_i. The second kind had another three types of samples built with different layer heights, which were 0.2, 0.15 and 0.1 mm. They were represented by TL0.2_i, TL0.15_i and TL0.1_i. The third kind of samples also had three different types, whose build directions were X, Z and Oblique, individually represented by TX0_i, TZ0_i and TO0_i. It is noted that only one processing parameter was set different from the original (or default) type of samples (i.e., TE0.4_i, TL0.15_i and TX0_i) for each other type based on fractional factorial design. Details of the processing parameter settings are provided in Table 1. To ensure the reliability and repeatability of the measurements, more than ten samples of each type were fabricated (thus, i > 10 ).
For the dynamic characteristic, another three kinds of dense-structure samples were fabricated with the same ME device, as the schematic shown in Figure 1b. They were also fabricated based on fractional factorial design, with the three processing parameters (i.e., extrusion width, layer height and build direction) individually set to the values shown in Table 1. These samples were separately represented by E0.4_j, E0.3_j, E0.2_j, L0.2_j, L0.15_j, L0.1_j, X0_j, Z0_j and O0_j, (j = 1–5). It is noted that TE0.4_i, E0.4_j, TL0.15_i, L0.15_j, TX0_i and X0_j were fabricated with the same processing parameters (the default values). All the samples mentioned above were made of polylactic acid (PLA [32], environmentally friendly material, having good workability and forming stability). Apart from only one single processing parameter set different for each type of samples, all the others were set the same, such as printing speed (moving speed of the printing head during the manufacturing process, set as 60 mm/s), extruder temperature (the temperature of the extruder to heat the filament to molten state, set as 200 °C), etc.

2.2. Theoretical Model

The 3D printing plate, composed of multiple layers of PLA fiber material, has the characteristic of orthogonal anisotropy. Figure 2 shows the geometric model of the rectangular plate, of which the length a was 130 mm, the width b was 50 mm and the total thickness h was 2.4 mm (the thickness of each single layer was the extrusion width of the extruded filament). Taking the central plane of the plate as the xoy plane, the three-dimensional coordinate system (O-xyz) was set up. The length, width and thickness directions are represented by x, y and z.
To model the plate, the following assumptions were made: (1) no delamination defect occurs inside the sample, the single layers are firmly bonded with each other, and no slip or relative displacement occurs between them. The sample can be regarded as an integral-structure plate, and there is no need to consider the interlayer coupling effect; (2) the bonding layer is thin, and no deformation occurs in itself, i.e., the deformation between the single-layer plates is continuous and (3) the plate is composed of multiple single-layer plates, but the total thickness is still consistent with the thin plate assumption. The modeling details are provided in the Appendix A.

2.3. Tensile Test

According to ISO 527-2-2012 standard, the tensile tests were carried out on the tensile test samples (TE0.4_i, TE0.3_i, TE0.2_i, TL0.2_i, TL0.15_i, TL0.1_i, TZ0_i, TX0_i and TO0_i, i > 10 ) using the testing machine (Shimadzu EHF-EV200k2-040, Figure 3). The loading rate during the test was set to 5 mm/min. Since the material was PLA (the mechanical strength is much lower than steel [32]), the clamping force at both ends was set to only 5 MPa to avoid breaking the samples, which ensured the accuracy in the measurements. It is noted that the results of the samples broken in the middle part during the tensile test were considered useful, and five typical testing results were chosen and summarized to represent the stress–strain relationship of the dynamic test samples, which were fabricated with the sample processing parameters as the tensile test ones.

2.4. Experimental Modal Analysis

To measure the dynamic characteristics of ME plates, modal tests were performed to the samples (E0.4_j, E0.3_j, E0.2_j, L0.2_j, L0.15_j, L0.1_j, X0_j, Z0_j and O0_j, j = 1–5). Figure 4 shows the testing system that mainly includes an impact hammer (PCB 086C01, with the sensitivity being 11.2 mV/N), a data acquisition card (NI USB 4431, with four analog input channels and one analog output channel) and a lightweight accelerometer (PCB 352C33, with the sensitivity being 103.5 mV/g). The plates were fixed with the clamping device and under the cantilever boundary conditions, with the length of the clamping area being 20 mm. The excitation point was about 10 mm away from the bottom, as shown in Figure 4. To ensure the accuracy of the measurements, the accelerometer was successively fixed to the measuring points on the plate with a large vibration response (the top and the middle positions of the sample, avoiding the node).
In the experiments, the impulse excitation was first applied to the samples by the hammer. The vibration responses of the samples due to the excitation were then measured by the accelerometer, with the sampling rate set as 2000 Hz. At the same time, the real-time excitation signals and response signals were collected through the data acquisition card. In order to obtain the mode shape of the samples, MISO (multiple-input single-output) method was used to identify the parameter, i.e., the accelerometer was fixed at a position with a large vibration response, and the other measuring positions were excited separately. It is noted that no filtering was implemented as no specific frequency band needed to be filtered out. For each sample, 10 groups of tests were completed which further ensures the accuracy of the measurements (350 groups of measured data obtained in total).

3. Results and Discussion

In this section, the material properties of the ME samples were provided, followed by the corresponding scanning electron microscope (SEM) analysis. The effect of three processing parameters (extrusion width, layer height and build direction) on the sample’s dynamic characteristics was then detailed.

3.1. Material Property of the Sample

The samples were weighted using the electrical analytical balance (FA2004), of which the precision was 0.1 mg. Together with the volume of the dynamic test plates, the density of each sample was determined. Based on Hooke’s law, the elastic modulus of the samples was obtained. According to the equation G   = E / 2 ( 1   +   υ ) , the ME plates’ shear elastic modulus could be determined, where υ is Poisson ratio, and υ = 0.36 [32]. Table 2 summarizes the average measured material properties of each type of dynamic test samples. Since E0.4, L0.15 and X0 were fabricated with the same processing parameters, they had the same properties.
It can be seen that the wider the extrusion width is, the larger the material property of the corresponding samples will be. This can be explained that the bonding area between adjacent extruded filaments will be decreased when the extrusion width decreases, as seen from the comparison between Figure 5a,b, leading to the deterioration of bonding strength. There would be more thermal cycles and more non-uniform thermal gradients, leading to more distortions, deformations and non-uniformities of the extruded filament (as shown from Figure 6a,b), which is a major reason for the bonded layers’ poor quality [27], and thus worse bonding quality. In addition, decreasing extrusion width would increase the bonding gap between the filaments of the cross section, generating more defects such as delamination, internal pores and void structures. Therefore, decreasing extrusion width will deteriorate the material property of the ME product.
Decreasing layer height can reduce the samples’ material property. When layer height is decreased, the number of layers needed to fabricate the sample will increase, resulting in the increase of bonding layers and thermal cycles, which further lead to more distortions and non-uniformities in the sample, as shown in Figure 6c. The bonding quality and strength between extruded filaments became worse. In addition, the defects like porosity, gap, separation, etc. would increase in the sample, as shown in Figure 5c. Therefore, the material property of the sample will be decreased when the layer height is reduced.
The samples built in the X direction had larger material properties than those built in other directions. As the build direction changed from X to Z, more layers were required to manufacture the part, resulting in more bonding layers, thermal cycles and distortions (Figure 6d, the distortion occurred due to high non-uniform temperature gradient between the bottom and the top layers during the printing). The bonding quality and strength would be reduced (see Figure 5d), and thus the material properties of the built part were deteriorated.

3.2. Effect of Extrusion Width on the Dynamic Property

Figure 7 presents the average measured frequency response curves of a set of ME plates fabricated with different extrusion widths separately being 0.4, 0.3 and 0.2 mm. It is noted that no filtering was implemented as no specific frequency band needed to be filtered out. By analyzing and identifying the peaks of the curves, the natural frequencies of each sample could be preliminarily determined. The results are detailed in Table 3. It can be seen that Sample E0.4 had the lowest values of resonant response (occurring due to the impulse excitation), meaning that the anti-vibration capability was the best compared with Sample E0.3 and E0.2. Therefore, the sample’s anti-vibration capability would be improved if the extrusion width was increased.
Table 4 summarizes the effect of different extrusion widths on the values of predicted and measured dynamic inherent characteristics. It can be seen from the comparison that the predictions agreed well with the measurements. The influencing trend of different extrusion width was predicted precisely. The discrepancy between predictions and measurements in natural frequency was in the range of 0.7–11%. The predicted mode shapes were the same as those of the measurements. The 1st order mode shape was the 1st order bending vibration, the 2nd order one was the 1st order torsional vibration and the 3rd order one was the 2nd order bending vibration. Therefore, the theoretical modeling based on the bidirectional beam function combination method could give reliable predictions on the inherent characteristics of ME plates. In addition, larger extrusion width mainly led to the increase in the material properties, i.e., elastic modulus, shear modulus and density, which would enhance the strain (or structural stiffness), potential and kinetic energy in the plates according to Equations (A12), (A14) and (A15). It eventually increased the natural frequency. For example, the values of material properties of Sample E0.4 were larger than those of E0.3 and E0.2 as shown in Table 2, the corresponding natural frequency was the highest, with the first three measured natural frequencies being 18 Hz, 154.9 Hz and 384.8 Hz. The corresponding predicted natural frequencies were 19.3 Hz, 170.1 Hz and 396.4 Hz. Since the structural characteristics were not changed with various extrusion widths, the mode shapes were kept the same.

3.3. Effect of Layer Height on the Dynamic Property

Figure 8 shows the effect of different layer height (0.1, 0.15 and 0.2 mm) on the frequency response of ME plates, with no filtering implemented. The results are detailed in Table 5. It could be seen that Sample L0.1 had the highest values of resonant response (occurring due to the impulse excitation), meaning that the anti-vibration capability was worse than Sample L0.15 and L0.2. Therefore, the sample’s anti-vibration capability would be improved if the layer height was increased.
The corresponding dynamic inherent characteristics are detailed in Table 6. It could be seen that the predicted natural frequencies and mode shapes agreed well with the measured data. The influencing trend of different layer heights was predicted precisely. The discrepancy of natural frequency was between 0.9% and 13.3%, the predicted mode shapes were the same as the measured ones. These confirmed the reliability of the theoretical modeling. Besides, different layer height caused the changes of elastic modulus, shear modulus and density of the samples, and eventually influenced the natural frequency. Since Sample L0.2 had the largest values of material property (shown in Table 2), the corresponding natural frequency was the highest. The first three measured natural frequencies were 20.1 Hz, 167.1 Hz and 392.3 Hz, with the predicted values being 21.3 Hz, 189.3 Hz and 434.3 Hz. Since different extrusion widths did not change the structural characteristics of the samples, the mode shapes were kept the same.

3.4. Effect of Build Direction on the Dynamic Property

Figure 9 provides the frequency response of the ME plates fabricated with different build directions, which were X, Z and Oblique directions. No filtering was implemented as no specific frequency band needed to be filtered out. Table 7 details the data. Since Sample X0 had the lowest values of resonant response, its anti-vibration capability was the best compared with Sample O0 and Z0. Therefore, the sample’s anti-vibration capability would be improved if it was built in the X direction.
Table 8 details the effect of different build directions on the inherent characteristics of ME plates. It could be seen that the influencing trend of different build directions was predicted precisely. The difference between predictions and measurements in natural frequency was from 2.3% to 13.4% and the predicted mode shapes were the same as the measurements. These further confirmed that the theoretical modeling was reliable in predicting the plate’s dynamic inherent characteristics. Furthermore, different build directions would change the samples’ elastic modulus, shear modulus and density, and further influenced the corresponding natural frequencies. For example, the material property values of Sample Z0 were the smallest among Z0, X0 and O0, as listed in Table 2, the corresponding natural frequency was the lowest, with the first three measured natural frequencies separately being 15.2 Hz, 130.3 Hz and 323 Hz. The predicted data was 14.1 Hz, 112.8 Hz and 290.9 Hz. Since different build directions did not change the structural characteristics of the samples, the mode shapes were kept the same.

4. Conclusions

Combined theory with experiment, this paper quantified the effect of extrusion width, layer height and build direction on the dynamic characteristics of material extrusion (ME) additive manufacturing plates under cantilever boundary conditions. The ME plate’s material properties will be deteriorated as the extrusion width decreases. Increasing layer height will improve the corresponding material property. When the build direction changes from X to Z, the material property of the sample will be reduced. Compared with the measurements, the theoretical model set up based on the bidirectional beam function combination method was validated, and it could give reliable predictions in the inherent characteristic of ME plates. The influencing trend of the processing parameters was predicted precisely. The natural frequency discrepancy was below 13.4%, and the predicted mode shapes were similar to the measured ones. In addition, different processing parameters, i.e., extrusion width, layer height and build direction, mainly led to the differences in samples’ elastic modulus, shear modulus and density. When these values were greater, the anti-vibration capability became better and the corresponding natural frequency was higher. As the structural characteristics were not changed with varying processing parameters, there was no change in the mode shape. Besides, the material used here was polylactic acid (PLA), which is just a type of common material. The modeling and testing methods mentioned in this present work could also be applied to other materials. This paper provided a theoretical basis and technical support for further research in improving the dynamic performance of ME products.

Author Contributions

Conceptualization, S.J.; data curation, Y.S. and Y.S.; formal analysis, S.J., Y.S. and Y.S.; funding acquisition, S.J.; investigation, Y.S., Y.S. and C.L.; methodology, S.J. and C.Z.; project administration, S.J. and M.Z.; resources, S.J., C.L. and C.Z.; software, Y.S.; supervision, S.J. and M.Z.; validation, S.J., Y.S. and Y.S.; writing—original draft preparation, Y.S. and Y.S.; writing—review and editing, S.J. and C.Z.

Funding

The work described was supported by National Natural Science Foundation of China (51705068) and the fundamental research funds for the central universities (N180703009).

Acknowledgments

The first author is grateful for the technical support from Hui Li, who is working in the same academic institute as the author.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

According to the small deflection theory of thin plates (Kirchoff hypothesis), the bidirectional beam function combination method [33] could be used to approximate the real mode shape of the plate under any boundary conditions (i.e., simply supported, fixed, free or any combinations), and further to analyze the inherent characteristics.
The plate’s mode shape can be expressed as
W p ( x , y ) = i = 1 n j = 1 n a i j X i ( x ) Y j ( y ) ,
where X i ( x ) is the i th order mode shape function meeting the boundary condition in x direction, Y j ( y ) is the jth order mode shape function meeting the boundary condition in y direction, a i j is an undetermined coefficient, which is utilized to adjust the combination of beam functions in different orders so as to make the mode shape function closer to the real mode shape of the rectangular plate.
The general expression of beam mode shapes is
Y ( x ) = A sin k x + B cos k x + C sh k x + D ch k x .
The x-axis direction is a free boundary condition, considering X i ( x ) as the nth order’s mode shape function of a free-free beam, the boundary condition needs to satisfy
When x = 0 and x = a,
E J d 2 Y d x 2 = 0 ,           d d x [ E J d 2 Y d x 2 ] = 0 .
Substitute Equation (A2) into Equation (A3), the expression of each order mode shape of the free-free beam is
X 1   =   1 X 2   =   3 ( 1 2 X / a ) , X i   =   ch k i X a   +   cos k i X a α i ( sh k i X a   +   sin k i X a ) ,
where k i is the coefficient related to beam frequency, k i 4 = ω 2 ( ρ A / E J ) ,   α i is the beam function coefficient and α i = ( ch k i a cos k i a ) / ( sh k i a sin k i a ) , ( i = 3, 4, 5…). The 1st order ( X 1 ) shows the movement of rigid body; the 2nd order ( X 2 ) shows the torsion of rigid body; starting from the 3rd order, the vibration mode is caused by the beam deformation.
In y-axis direction, the plate can be seen as a fixed-free beam, the boundary condition needs to satisfy
When y = 0,
Y = 0 ,     d Y d x = 0 .
When y = a,
E J d 2 Y d x 2 = 0 ,           d d x [ E J d 2 Y d x 2 ] = 0 .
Substitute Equation (A2) into Equations (A5) and (A6), the expression of each order’s mode shape function of a fixed-free beam is
Y i = ch k j y b cos k j y b α j ( sh k j y b sin k j y b ) ,
where k j is the same as k i (mentioned above) and α j = ( ch k j b cos k j b ) / ( sh k j b sin k j b ) , ( i = 3, 4, 5…). The frequency coefficient [33] is detailed in Table A1.
Substitute Equations (A4) and (A7) into Equation (A1), the mode shape function of the ME plate under cantilever boundary condition is given by
W p = X 1 · Y 1 · a 11 + X 1 · Y 2 · a 12 + X 1 · Y 3 · a 13 + X 2 · Y 1 · a 21 + X 2 · Y 2 · a 22 + X 2 · Y 3 · a 23 + X 3 · Y 1 · a 31 + X 3 · Y 2 · a 32 + X 3 · Y 3 · a 33 +
where {     X 1 = 1 X 2 = 3 ( 1 2 x / a ) X 3 = ch 4.73 x a + cos 4.73 x a ch k i a cos k i a sh k i a sin k i a ( sh 4.73 x a + sin 4.73 x a ) Y 1 = ch 1.875 y b cos 1.875 y b ch 1.875 b cos 1.875 b ch 1.875 b cos 1.875 b ( sh 1.875 y b sin 1.875 y b ) Y 2 = ch 4.694 y b cos 4.694 y b ch 4.694 b cos 4.694 b ch 4.694 b cos 4.694 b ( sh 4.694 y b sin 4.694 y b ) Y 3 = ch 4.694 y b cos 4.694 y b ch 7.855 b cos 7.855 b ch 7.855 b cos 7.855 b ( sh 7.855 y b sin 7.855 y b ) .
Table A1. The coefficient of beam frequency.
Table A1. The coefficient of beam frequency.
Boundary Conditions α i
x = 0   i = 1 i = 2 i = 3
FreeFree004.73
FixedFree1.8754.6947.855
Since the accurate theoretical results can hardly be determined by the analytical method, the Ritz method was adopted in this paper. It is an approximate solution method according to the energy variation principle, which could obtain the results by transforming functional extremum problem to the extremum problem of multivariate function. Due to the vibration excitation, the plate would be deformed, producing stress and strain. Therefore, there is deformation energy varying with time existing in the plate. The strain energy density can be expressed by
W p   =   1 2 ( σ X ε X + σ Y ε Y + τ xy γ xy ) .
Set the deformation function of central plane as w 0 ( x , y , t ) , the displacement function of any point in the plate can be given by
u ( x ,   y , z , t ) = z w 0 ( x , y , t ) x , v ( x ,   y , z , t ) = z w 0 ( x , y , t ) y .
According to geometric equation, the strain of any point in the plate can be obtained as
ε x = u x = z 2 w x 2 = z κ x , ε y = v y = z 2 w y 2 = z κ y , γ x y = u y + v x = 2 z 2 w x y = 2 z κ x y ,
where κ x , κ y and κ x y are the curvature and distortion of the central plane.
According to stress–strain relationship, the stress of any point in the plate could be expressed by
σ x = E 1 ν 2 ( ε x + ν ε y )   = E 1 ν 2 z ( 2 w x 2 + ν 2 w x 2 ) , σ y = E 1 ν 2 ( ε y + ν ε y )   = E 1 ν 2 z ( 2 w y 2 + ν 2 w x 2 ) ,   τ x y = G γ x y = 2 G   z 2 w x y ,
where ν is poisson’s ratio, E is elastic modulus and G is shear modulus.
Substituting Equations (A11) and (A12) into Equation (A10) gives
W p = 1 2 [ 1 E ( σ x 2 + σ y 2 2 ν σ x σ y ) + 1 G τ x y 2 ] = E z 2 2 ( 1 ν 2 ) [ ( 2 w x 2 ) 2 + ( 2 w y 2 ) 2 + 2 ( 1 ν ) ( 2 w x y ) 2 ] .
Therefore, the potential energy of the plate is given by
U = D 2 A { ( 2 w x 2 + 2 w y 2 ) 2 2 ( 1 υ ) [ ( 2 w x 2 ) ( 2 w y 2 ) ( 2 w x y ) 2 ] } d x d y .
The kinetic energy of the plate is given by
T = ρ h 2 A ( 2 W t 2 ) 2 d x d y .
Substitute the expression of displacement W ( x , y , t ) = W ( x , y ) sin ( w t + φ ) into Equations (A14) and (A15), and then into the variational formula based on Hamilton principle ( δ t 1 t 2 ( T U ) d t = 0 ), the time integral yields
δ ( U m a x T m a x ) = 0 .
U m a x = D 2 A { ( 2 w x 2 + 2 w y 2 ) 2 2 ( 1 υ ) [ ( 2 w x 2 ) ( 2 w y 2 ) ( 2 w x y ) 2 ] } d x d y .
T m a x = W 2 ρ h 2 A W 2 d x d y .
Substitute the mode shape function (Equation (A8)) into Equations (A17) and (A18), and let the partial derivative of U m a x T m a x in terms of the undetermined coefficient α i j equal 0, that is
d U m a x d α i j d T m a x d α i j = 0 .
The linear algebraic equation for α i j is thus obtained. Set the coefficient determinant of the equation to be 0, the circular frequency ω can be determined. According to f   =   ω / 2 π , the natural frequency of the ME plate can be obtained.

References

  1. Mohamed, O.A.; Masood, S.H.; Bhowmik, J.L. Optimization of fused deposition modeling process parameters: A review of current research and future prospects. Adv. Manuf. 2015, 3, 42–53. [Google Scholar] [CrossRef]
  2. Rezayat, H.; Zhou, W.; Siriruk, A.; Penumadu, D.; Babu, S.S. Structure-mechanical property relationship in fused deposition modelling. Mater. Sci. Technol. 2015, 31, 895–903. [Google Scholar] [CrossRef]
  3. Wu, P.; Wang, J.; Wang, X. A critical review of the use of 3-D printing in the construction industry. Autom. Constr. 2016, 68, 21–31. [Google Scholar] [CrossRef] [Green Version]
  4. Cao, Q.; Xie, H. Characterization for elastic constants of fused deposition modelling-fabricated materials based on the virtual fields method and digital image correlation. Acta Mech. Sin. 2017, 33, 1075–1083. [Google Scholar] [CrossRef]
  5. Hwang, S.; Reyes, E.I.; Kim, N.S.; Moon, K.S.; Rumpf, R.C. Thermo-mechanical characterization of metal/polymer composite filaments and printing parameter study for fused deposition modeling in the 3D printing process. J. Electron. Mater. 2015, 44, 771–777. [Google Scholar] [CrossRef]
  6. Torres, J.; Cotelo, J.; Karl, J.; Gordon, A.P. Mechanical property optimization of FDM PLA in shear with multiple objectives. JOM 2015, 67, 1183–1193. [Google Scholar] [CrossRef]
  7. Tsouknidas, A.; Pantazopoulos, M.; Katsoulis, I.; Fasnakis, D.; Maropoulos, S.; Michailidis, N. Impact absorption capacity of 3D-printed components fabricated by fused deposition modelling. Mater. Des. 2016, 102, 41–44. [Google Scholar] [CrossRef]
  8. Li, H.; Wang, T.; Sun, J.; Yu, Z. The effect of process parameters in fused deposition modelling on bonding degree and mechanical properties. Rapid Prototyp. J. 2018, 24, 80–92. [Google Scholar] [CrossRef]
  9. Rodríguez-Panes, A.; Claver, J.; Camacho, A. The influence of manufacturing parameters on the mechanical behaviour of pla and abs pieces manufactured by fdm: A comparative analysis. Materials 2018, 11, 1333. [Google Scholar] [CrossRef] [Green Version]
  10. Cuan-Urquizo, E.; Barocio, E.; Tejada-Ortigoza, V.; Pipes, R.B.; Rodriguez, C.A.; Roman-Flores, A. Characterization of the mechanical properties of FFF structures and materials: A review on the experimental, computational and theoretical approaches. Materials 2019, 12, 895. [Google Scholar] [CrossRef] [Green Version]
  11. Durgun, I.; Ertan, R. Experimental investigation of FDM process for improvement of mechanical properties and production cost. Rapid Prototyp. J. 2014, 20, 228–235. [Google Scholar] [CrossRef]
  12. Kim, M.K.; Lee, I.H.; Kim, H.-C. Effect of fabrication parameters on surface roughness of FDM parts. Int. J. Precis. Eng. Manuf. 2018, 19, 137–142. [Google Scholar] [CrossRef]
  13. Boschetto, A.; Bottini, L.; Veniali, F. Finishing of Fused Deposition Modeling parts by CNC machining. Robot. Comput. Integr. Manuf. 2016, 41, 92–101. [Google Scholar] [CrossRef]
  14. Lalehpour, A.; Barari, A. A more accurate analytical formulation of surface roughness in layer-based additive manufacturing to enhance the product’s precision. Int. J. Adv. Manuf. Technol. 2018, 96, 3793–3804. [Google Scholar] [CrossRef]
  15. Li, Z.; Zhang, Z.; Shi, J.; Wu, D. Prediction of surface roughness in extrusion-based additive manufacturing with machine learning. Robot. Comput. Integr. Manuf. 2019, 57, 488–495. [Google Scholar] [CrossRef]
  16. Kestilä, A.; Nordling, K.; Miikkulainen, V.; Kaipio, M.; Tikka, T.; Salmi, M.; Auer, A.; Leskelä, M.; Ritala, M. Towards space-grade 3D-printed, ALD-coated small satellite propulsion components for fluidics. Addit. Manuf. 2018, 22, 31–37. [Google Scholar] [CrossRef]
  17. Sun, Q.; Rizvi, G.M.; Bellehumeur, C.T.; Gu, P. Effect of processing conditions on the bonding quality of FDM polymer filaments. Rapid Prototyp. J. 2008, 14, 72–80. [Google Scholar] [CrossRef]
  18. Gordeev, E.G.; Galushko, A.S.; Ananikov, V.P. Improvement of quality of 3D printed objects by elimination of microscopic structural defects in fused deposition modeling. PLoS ONE 2018, 13, e0198370. [Google Scholar] [CrossRef]
  19. Caminero, M.A.; Chacón, J.M.; García-Moreno, I.; Reverte, J.M. Interlaminar bonding performance of 3D printed continuous fibre reinforced thermoplastic composites using fused deposition modelling. Polym. Test. 2018, 68, 415–423. [Google Scholar] [CrossRef]
  20. Narahara, H.; Shirahama, Y.; Koresawa, H. Improvement and evaluation of the interlaminar bonding strength of FDM parts by Atmospheric-Pressure Plasma. Procedia CIRP 2016, 42, 754–759. [Google Scholar] [CrossRef]
  21. Kretzschmar, N.; Chekurov, S.; Salmi, M.; Tuomi, J. Evaluating the readiness level of additively manufactured digital spare parts: An industrial perspective. Appl. Sci. 2018, 8, 1837. [Google Scholar] [CrossRef] [Green Version]
  22. Melenka, G.W.; Schofield, J.S.; Dawson, M.R.; Carey, J.P. Evaluation of dimensional accuracy and material properties of the MakerBot 3D desktop printer. Rapid Prototyp. J. 2015, 21, 618–627. [Google Scholar] [CrossRef]
  23. Boschetto, A.; Bottini, L. Design for manufacturing of surfaces to improve accuracy in Fused Deposition Modeling. Robot. Comput. Integr. Manuf. 2016, 37, 103–114. [Google Scholar] [CrossRef]
  24. Turner, B.N.; Gold, S.A. A review of melt extrusion additive manufacturing processes: II. Materials, dimensional accuracy, and surface roughness. Rapid Prototyp. J. 2015, 21, 250–261. [Google Scholar] [CrossRef]
  25. Nunez, P.J.; Rivas, A.; Garcia-Plaza, E.; Beamud, E.; Sanz-Lobera, A. Dimensional and surface texture characterization in Fused Deposition Modelling (FDM) with ABS plus. In Mesic Manufacturing Engineering Society International Conference 2015; Elsevier Science Bv: Amsterdam, The Netherlands, 2015; Volume 132, pp. 856–863. [Google Scholar]
  26. Arivazhagan, A.; Masood, S. Dynamic mechanical properties of ABS material processed by fused deposition modelling. Int. J. Eng. Res. Appl. 2012, 2, 2013–2014. [Google Scholar]
  27. Mohamed, O.A.; Masood, S.H.; Bhowmik, J.L.; Nikzad, M.; Azadmanjiri, J. Effect of process parameters on dynamic mechanical performance of FDM PC/ABS printed parts through design of experiment. J. Mater. Eng. Perform. 2016, 25, 2922–2935. [Google Scholar] [CrossRef]
  28. Mohamed, O.A.; Masood, S.H.; Bhowmik, J.L. Investigation of dynamic elastic deformation of parts processed by fused deposition modeling additive manufacturing. Adv. Prod. Eng. Manag. 2016, 11, 227–238. [Google Scholar] [CrossRef] [Green Version]
  29. 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. Mater. Manuf. Process. 2016, 31, 1983–1994. [Google Scholar] [CrossRef]
  30. Domingo-Espin, M.; Borros, S.; Agullo, N.; Garcia-Granada, A.A.; Reyes, G. Influence of building parameters on the dynamic mechanical properties of polycarbonate fused deposition modeling parts. 3D Print. Addit. Manuf. 2014, 1, 70–77. [Google Scholar] [CrossRef]
  31. Mansour, M.; Tsongas, K.; Tzetzis, D.; Antoniadis, A. Mechanical and dynamic behavior of fused filament fabrication 3D printed polyethylene terephthalate glycol reinforced with carbon fibers. Polym. Plast. Technol. Eng. 2018, 57, 1715–1725. [Google Scholar] [CrossRef]
  32. Farah, S.; Anderson, D.G.; Langer, R. Physical and mechanical properties of PLA, and their functions in widespread applications—A comprehensive review. Adv. Drug Deliv. Rev. 2016, 107, 367–392. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Berdichevsky, V.L. Theory of elastic plates and shells. In Variational Principles of Continuum Mechanics: II, Applications; Springer: Berlin/Heidelberg, Germany, 2009; pp. 589–714. [Google Scholar]
Figure 1. Two-dimensional drawing of the samples. (a) Tensile test sample and (b) dynamic test sample.
Figure 1. Two-dimensional drawing of the samples. (a) Tensile test sample and (b) dynamic test sample.
Applsci 09 05345 g001
Figure 2. Mechanism graph of the material extrusion (ME) plate.
Figure 2. Mechanism graph of the material extrusion (ME) plate.
Applsci 09 05345 g002
Figure 3. Tensile test machine (servo-hydraulic) and the sample.
Figure 3. Tensile test machine (servo-hydraulic) and the sample.
Applsci 09 05345 g003
Figure 4. Testing system for the determination of the dynamic characteristics of ME plates.
Figure 4. Testing system for the determination of the dynamic characteristics of ME plates.
Applsci 09 05345 g004
Figure 5. Scanning electron microscope (SEM) micrograph of the sample’s cross section. (a) Sample E0.4 (or L0.15, X0); (b) Sample E0.2; (c) Sample L0.1 and (d) Sample Z0.
Figure 5. Scanning electron microscope (SEM) micrograph of the sample’s cross section. (a) Sample E0.4 (or L0.15, X0); (b) Sample E0.2; (c) Sample L0.1 and (d) Sample Z0.
Applsci 09 05345 g005aApplsci 09 05345 g005b
Figure 6. Surface SEM observation of the sample. (a) Sample E0.4 (or L0.15, X0); (b) Sample E0.2; (c) Sample L0.1 and (d) Sample Z0.
Figure 6. Surface SEM observation of the sample. (a) Sample E0.4 (or L0.15, X0); (b) Sample E0.2; (c) Sample L0.1 and (d) Sample Z0.
Applsci 09 05345 g006
Figure 7. The average measured frequency response of the samples fabricated with different extrusion width.
Figure 7. The average measured frequency response of the samples fabricated with different extrusion width.
Applsci 09 05345 g007
Figure 8. The average measured frequency response of the samples fabricated with different layer height.
Figure 8. The average measured frequency response of the samples fabricated with different layer height.
Applsci 09 05345 g008
Figure 9. The average measured frequency response of the samples fabricated with different build directions.
Figure 9. The average measured frequency response of the samples fabricated with different build directions.
Applsci 09 05345 g009
Table 1. The key processing parameter settings of the samples.
Table 1. The key processing parameter settings of the samples.
Sample
(i > 10)
(j = 1–5)
Extrusion Width (mm)Layer Height (mm)Build DirectionLongitudinal
Section
Number of
Contours
TE0.4_i0.40.15X direction (0°)
Applsci 09 05345 i001
Applsci 09 05345 i0023
E0.4_j
TE0.3_i0.30.15 Applsci 09 05345 i0034
E0.3_j
TE0.2_i0.20.15 Applsci 09 05345 i0046
E0.2_j
TL0.2_i0.40.2X direction (0°) Applsci 09 05345 i0053
L0.2_j
TL0.15_i0.40.15 Applsci 09 05345 i0063
L0.15_j
TL0.1_i0.40.1 Applsci 09 05345 i0073
L0.1_j
TZ0_i0.40.15Z direction (90°)
Applsci 09 05345 i008
Applsci 09 05345 i0093
Z0_j
TX0_i0.40.15X direction (0°) Applsci 09 05345 i0103
X0_j
TO0_i0.40.15Oblique direction (45°)
Applsci 09 05345 i011
Applsci 09 05345 i0123
O0_j
Table 2. Average material properties of material extrusion (ME) samples.
Table 2. Average material properties of material extrusion (ME) samples.
OrderSampleAverage Elastic Modulus
(MPa)
Average Shear Modulus
(MPa)
Average Density
(g/cm³)
Poisson’s Ratio
1E0.4, L0.15 and X018.86.910.920.36
2E0.318.56.800.91
3E0.218.06.620.91
4L0.219.37.100.92
5L0.118.16.650.91
6Z017.16.290.90
7O017.66.470.93
Table 3. The resonant response values of the samples fabricated with different extrusion widths.
Table 3. The resonant response values of the samples fabricated with different extrusion widths.
SampleAverage Resonant Response (m/s2 N)
1st Order2nd Order3rd Order
E0.44.34.43.0
E0.35.05.03.6
E0.25.25.84.0
Table 4. The effect of different extrusion width on the samples’ inherent characteristics.
Table 4. The effect of different extrusion width on the samples’ inherent characteristics.
ItemSamplesOrder
1st2nd3rd
Average measured natural frequency A (Hz)E0.418.0154.9384.8
E0.317.1145.6371.1
E0.216.1134.2346.0
Average predicted natural frequency B (Hz)E0.419.3170.1396.4
E0.318.5161.6373.6
E0.217.0147.6334.4
Error (%)
| B A | / A
E0.47.29.83.0
E0.38.211.00.7
E0.25.610.03.4
Measured mode shape Applsci 09 05345 i013 Applsci 09 05345 i014 Applsci 09 05345 i015
Predicted mode shape Applsci 09 05345 i016 Applsci 09 05345 i017 Applsci 09 05345 i018
Table 5. The resonant response values of the samples fabricated with different layer heights.
Table 5. The resonant response values of the samples fabricated with different layer heights.
SampleAverage Resonant Response (m/s2 N)
1st Order2nd Order3rd Order
L0.14.85.05.1
L0.154.34.43.0
L0.24.03.52.7
Table 6. The effect of different layer height on samples’ inherent characteristics.
Table 6. The effect of different layer height on samples’ inherent characteristics.
ItemSamplesOrder
1st2nd3rd
Average measured natural frequency A (Hz)L0.116.6148.9375.1
L0.1518.0154.9384.8
L0.220.1167.1392.3
Average predicted natural frequency B (Hz)L0.117.0150.2342.7
L0.1519.3170.1396.4
L0.221.3189.3434.3
Error (%)
| B A | / A
L0.12.40.98.6
L0.157.29.83.0
L0.26.013.310.7
Measured mode shape Applsci 09 05345 i019 Applsci 09 05345 i020 Applsci 09 05345 i021
Predicted mode shape Applsci 09 05345 i022 Applsci 09 05345 i023 Applsci 09 05345 i024
Table 7. The resonant response values of the samples fabricated with different build directions.
Table 7. The resonant response values of the samples fabricated with different build directions.
SampleAverage Resonant Response (m/s2 N)
1st Order2nd Order3rd Order
X04.34.43.0
O04.84.93.9
Z04.87.75.3
Table 8. The effect of different build directions on the samples’ inherent characteristics.
Table 8. The effect of different build directions on the samples’ inherent characteristics.
ItemSamplesOrder
1st2nd3rd
Average measured natural frequency A (Hz)Z015.2130.3323
X017.8155.1387.6
O016.7145.5346.5
Average predicted natural frequency B (Hz)Z014.1112.8290.9
X019.3170.1396.4
O016.1132.2311.7
Error (%)
| B A | / A
Z07.213.49.9
X08.49.72.3
O03.69.110.0
Measured mode shape Applsci 09 05345 i025 Applsci 09 05345 i026 Applsci 09 05345 i027
Predicted mode shape Applsci 09 05345 i028 Applsci 09 05345 i029 Applsci 09 05345 i030

Share and Cite

MDPI and ACS Style

Jiang, S.; Shi, Y.; Siyajeu, Y.; Zhan, M.; Zhao, C.; Li, C. Effect of Processing Parameters on the Dynamic Characteristic of Material Extrusion Additive Manufacturing Plates. Appl. Sci. 2019, 9, 5345. https://doi.org/10.3390/app9245345

AMA Style

Jiang S, Shi Y, Siyajeu Y, Zhan M, Zhao C, Li C. Effect of Processing Parameters on the Dynamic Characteristic of Material Extrusion Additive Manufacturing Plates. Applied Sciences. 2019; 9(24):5345. https://doi.org/10.3390/app9245345

Chicago/Turabian Style

Jiang, Shijie, Yinfang Shi, Yannick Siyajeu, Ming Zhan, Chunyu Zhao, and Changyou Li. 2019. "Effect of Processing Parameters on the Dynamic Characteristic of Material Extrusion Additive Manufacturing Plates" Applied Sciences 9, no. 24: 5345. https://doi.org/10.3390/app9245345

APA Style

Jiang, S., Shi, Y., Siyajeu, Y., Zhan, M., Zhao, C., & Li, C. (2019). Effect of Processing Parameters on the Dynamic Characteristic of Material Extrusion Additive Manufacturing Plates. Applied Sciences, 9(24), 5345. https://doi.org/10.3390/app9245345

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop