Next Article in Journal
Irradiation Damage Independent Deuterium Retention in WMoTaNbV
Previous Article in Journal
Feasibility Analysis of Machining Cobalt-Chromium Alloy (Stellite-6) Using TiN Coated Binary Inserts
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Injection Molding Process Simulation of Polycaprolactone Sticks for Further 3D Printing of Medical Implants

Department of Mechanical Engineering Fundamentals, Faculty of Mechanical Engineering and Computer Science, University of Bielsko-Biala, 43-309 Bielsko-Biala, Poland
*
Author to whom correspondence should be addressed.
Materials 2022, 15(20), 7295; https://doi.org/10.3390/ma15207295
Submission received: 30 September 2022 / Revised: 13 October 2022 / Accepted: 17 October 2022 / Published: 18 October 2022

Abstract

:
The aim of the present study was a simulation of the injection molding process of polycaprolactone filament sticks for further 3D printing of osteochondral implants. Polycaprolactone data are not available in the data banks of popular injection molding simulation programs. Therefore, thermal and rheological data from the literature were imported to the material database of Solidworks Plastics software to simulate the injection molding process of filament sticks. The influence of several injection molding parameters including melt temperature, injection time, and injection pressure on the geometry of filament stick (final part) was investigated. Based on the results of the performed simulation and analyses, it was possible to improve the injection process parameters. The accuracy of simulation predictions, based on the literature data, demonstrates the potential of using simulation as a tool to develop polycaprolactone parts for future implants and to optimize the injection molding process.

1. Introduction

Facial injuries are a growing clinical problem not only for patients, but also for laryngologist and plastic surgeons. The main cause of injuries is various types of accidents, as well as diseases, congenital deformities, or cancer of supporting tissues such as bone and cartilage. Treatment of osteochondral defects remains a great challenge for many reconstructive surgical procedures.
Three-dimensional (3D) printing is revolutionizing and enhancing the medical industry as it allows for the personalization of the treatment to match a patient’s individual needs. FDM is a 3D printing techniques and offers great flexibility in the handling and processing of materials. In this method, implants are created through the layer-by-layer deposition of a thermoplastic polymer using a 3D printer. The polymeric filament is loaded into the 3D printer, heated above its melting point, softened, and extruded through the nozzle. Usually, the polymers used for biomedical applications are very expensive. Filament materials are generally supplied in spools that are not entirely used during the production of small implants for cartilage or bone replacement [1]. In order to produce filament without such a consumption of biomedical materials, the idea of small sticks that can be joined together was developed by the authors. Moreover, this will provide us an opportunity to put together differently modified sticks and produce personalized implants using a commercially available 3D printer. For their production, injection molding can be applied. The results of the injection molding simulation process presented in this paper will also help us to set the production parameters for the polymer modified with various additives such as hydroxyapatite, bioglass, graphene, drugs, or antibiotics. Therefore, soon, instead of buying entire filament spool, it will be possible to buy a set of differently modified sticks and organize them depending on the patient’s requirements, using a 3D printer. We will be able to produce various scaffolds without having to change the spool as it is in the traditional FDM process. The sticks developed during our project are designed to be used in commercially available printers without any hardware or software modification.
Injection molding (IM) is a method used to obtain products by injecting molten polymers into the mold cavity with high pressure and then cooling and solidifying them. When the specimen has cooled down, the mold opens and the solid part is ejected, before the mold closes, starting the cycle again [2]. IM is able to manufacture products with precise dimension steadiness, low manufacture cycles, and low costs. Compared with other techniques, injection molding will be found to feature very small tolerances. Therefore, it is also widely used for the production of medical parts. Having fabricated a few sets of filament sticks, it will be possible to join them and obtain filament that can be used in standard or mid-range FDM 3D printer [1]. The flow of molten thermoplastic material in the runners and cavity during the injection molding process is a complex phenomenon that depends on various factors. The viscosity of polymer melt, polymer flow rate, melt temperature, and wall thickness all contribute to the behavior of polymer melt flow [3]. Moreover, the processing parameters also have a great impact on the final properties of the molded product, i.e., filament stick. The changes in processing conditions generate changes in molecular orientation, residual stress, molecular degradation, and the degree of crystallinity [3]. The control of the injection molding machine and the design of the mold are a great challenge for the manufacturer. The lack of knowledge about the material and process conditions leads to molding of parts with some critical defects such as shrinkage and warpage [4]. The type of polymer and its structure affect the shrinkage; however, most of those factors arise from process conditions, such as the injection time, packing time, melt temperature, mold surface temperature, and packing pressure. Warpage can be a result of non-uniform packing pressures in the mold cavities, non-uniform molding temperatures or cooling rates, and non-uniform wall thicknesses [4,5]. When it comes to product quality in injection molding, the control of process parameters is fundamental to create a stable manufacturing process and high quality parts. In the industries, the trial and error method is normally practiced to achieve a desirable product; however, it is a very time-consuming method. Therefore, nowadays, computational simulation of injection molding processes has been widely used. It is a powerful tool that can be used in the early stages of part development to identify and correct mistakes in injection molds and to avoid the risks of costly re-engineering. The commercial CAE packages such as Autodesk Moldflow, Moldex3D, and Solidworks Plastics are used to eliminate the traditional trial and error approach and can provide useful assistance for the optimization of process parameters [6]. However, there are not many types of research for injection molding simulation of the biomedical parts and the material data for biomedical polymers are seldom available from data banks of simulation tools [4].
Polycaprolactone (PCL) is a semi-crystalline aliphatic polyester with high biocompatibility, approved by the FDA [7,8]. Polycaprolactone is non-toxic and biodegradable; the good solubility of PCL, as well as its low melting point (59–64 °C), long-term degradation properties (>24 months to lose total mass), and exceptional blend-compatibility, have stimulated extensive research into its potential application in the biomedical field [9]. PCL displays a high viscosity in the melt phase. Therefore, elevated temperatures (above 140 °C) are required to process it [10]. PCL has been used for various biomedical applications such as controlled drug delivery system, tissue-engineering scaffolds, and wound healing [11]. As we mentioned, PCL has a longer degradation time than other polymers from the group of aliphatic polyesters, thus biomaterials made using this polymer can be implanted in areas subjected to increased load [12]. The use of PCL, which has better mechanical parameters, to make scaffolds for the treatment of osteochondral defects allows for maintaining sufficient space for extracellular matrix formation, not only until the moment of producing new tissue, but until it obtains the required mechanical parameters [13]. Unfortunately, PCL data are not available in the data banks of popular injection molding simulation programs. Therefore, in this study, we have create polycaprolactone material data from the literature and use Solidworks Plastic software to simulate the process of injection molding in order to evaluate the influence of operational conditions and mold design on the efficiency of an injection process used to produce PCL filament stick for further 3D printing of facial implants. Thermal and rheological data from the literature were imported to the material database of Solidworks Plastics software for simulating the injection molding process of filament sticks. The influence of several injection molding parameters including melt temperature, injection time, packing time, and packing pressure on the geometry of final part (polymer sticks that will be used for 3D printing) was investigated.

2. Materials and Methods

Solidworks Plastics was used for the process of polymer injection simulation and material model building. The material model was created based on PCL polymer CAPA® 6500 with a molecular weight of 50 kDa. All thermal and rheological data necessary for stimulation were collected from the literature and imported to the material database Solidworks Plastics software. CAD Siemens NX was used to build the CAD model of the filament sticks.
In order to verify the simulation process, a Babyplast 6/10P injection molding machine was used to make PCL stick filaments. Chrome-molybdenum 1.2316 X38CrMo16 cold work tool steel was used to make the mold. The size of the mold was 75 mm × 75 mm × 76 mm (height) and the dimensions of the molded part were 60.6 mm (along the sticks) × 57.8 mm (across the sticks) × 23 mm (height).

2.1. Material Data Collection

The viscosity of polymers is a very important property in determining the rheological nature of plastics. In order to accurately predict the viscosity of PCL polymer with respect to shear rate, the cross-WLF viscosity model was used. The cross-WLF viscosity model describes the temperature, shear rate, and pressure dependency of the viscosity for thermoplastic materials. The viscosity with the corresponding shear rate is determined through the following formulas:
η = η 0 1 + ( η 0 γ ˙ τ * ) 1 n
η 0 = D 1 exp [ A 1 ( T T ref ) A 2 + ( T T ref ) ]
where η is the melt viscosity, Pa∙s; η0 is the zero viscosity, Pa∙s; γ ˙ is the shear rate, s−1; τ is the critical shear stress, Pa; n is the power law index, dimensionless; T is the temperature, K; and D1, A1, A2, and Tref are the parameters of the cross-WLF equation.
Tref = D2 + D3·P
D2 and D3 are model coefficients determined from experimental measurements for a specific polymer. In most cases, the effect of pressure on viscosity is neglected; therefore, D3 is considered to be zero, meaning that the pressure influence is ignored. The cross-WLF viscosity model use the materials’ constants (data-fitted constants) along with predetermined values of temperature and shear rate to calculate the viscosity at a required point. The rheological parameters values were selected on the basis of the literature [14,15] (Table 1) and imported to the SolidWorks program to generate the cross-WLF viscosity model (Figure 1a). The thermal characteristics of the polycaprolactone were obtained from the literature [16,17,18] and are reported in Table 1. The information contains the following: transition temperature, melt flow rate (MFR), thermal conductivity, and melting point. Physical and mechanical properties of polycaprolactone such as solid density, tensile modulus, yield stress, and Poisson’s ratio were also obtained from the literature and are reported in Table 1 [18,19,20].
The simulation of the injection molding process requires knowledge of the time and temperature dependence of the material density at the pressures that occur during the process [21]. The quantity of shrinkage in the injection molded part is a direct function of PVT behavior [22]. Therefore, the pressure–volume–temperature (PVT) relationship is the key rule to characterize the change in density with respect to pressure and temperature [23]. For PVT, we used the Tait model:
υ ( T ,   p ) = υ 0 ( T ) [ 1 C l n ( 1 + p B ( T ) ) ] + υ t ( T , p )
where C = 0.0894
υ 0 ( T ) = { b 1 m + b 2 m ( T b 5 )     i f   T > T t   b 1 s + b 2 s ( T b 5 )     i f   T < T t
B ( T ) = { b 3 m e x p [ b 4 m ( T b 5 ) ]     i f   T > T t   b 3 s e x p [ b 4 s ( T b 5 ) ]     i f   T < T t
υ t ( T , p ) = {   0                                           i f   T > T t   b 7 e x p [ b 8 ( T b 5 ) b 9 p ]   i f   T < T t
T t ( p ) = b 5 + b 6 · p
where v0 is the specific volume at zero pressure; b1 (b1m, b1s) and b2 (b2m, b2s) are the parameters to describe the dependence of v0 on pressure and temperature; B is the sensitivity to pressure and a function of temperature alone with two material constants b3 (b3m, b3s) and b4 (b4m, b4s); vt is the specific volume decrease due to crystallization; Tt is the transition temperature; b5 and b6 are parameters that describe the change of Tt with pressure; and b7, b8, and b9 are particular parameters of semi-crystalline polymers that describe the form of the state transition [23,24].
All parameters are presented in Table 2 and were generated using SolidWorks Plastic based on plot (Figure 1b).

2.2. Injection Molding Simulation

Simulations were carried out using SolidWorks Plastics 2021 software 29.3.0.0059, Dassault Systemes, Pars, France. The complete project of filament mold based on designed part was executed in the CAD Siemens NX (computer-aided design). An image of the filament model and the CAD model of injection mold are shown in Figure 2. Based on the full CAD assembly of the mold, the simplified CAD model was built. Additional parts such as the inset, plugs, and cooling channels were deleted or transformed to one solid geometry. Figure 3a,b shows the upper and the lower side of the simplified CAD model. It was used to isolate the solid geometry of the injection part. The model of the part was isolated from the simplified mold model as a result of the use of Boolean operations (geometry of solids intersection) (Figure 3c). Figure 3d shows the simplified injection part model with only one side (four sticks). The aim of this simplification was to obtain a shorter calculation time. This operation was possible because the part is symmetrical. All models were prepared in CAD software.
Based on the CAD model of the injection part, the mesh model was developed. The mesh subdivides domains of the simulation into discrete cells (in this case, the runner system and cavity area). The solid mesh procedure subdivides the volume of each domain into polyhedral-shaped elements. These solid elements have a tetrahedral shape. The software creates a surface mesh too, but the surface elements (faces) represent the boundaries of the closed solid polyhedral elements. A solid mesh is suitable for models with a complex geometry. The mesh model was prepared using injection molding simulation software. Figure 3e shows the mesh model of the injection part. Table 3 contains the mesh parameters. Figure 3f shows the runner system and sticks area of the injection part. The mesh of sticks is more precise (higher number of tetrahedral elements). This means the analysis results will be more accurate. Figure 3g presents the discretization of stick volume.
Figure 4 shows the longitudinal section through the mesh model of stick. Figure 4b,c presents the joint regions of sticks. The areas are very small (0.8–1.2 mm diameter). Discretization in these areas is much more precise. Additionally, two external boundary layers were manually added at the outside of sticks for better accuracy of calculations.
Additionally, two external boundary layers were manually added at the outside of the sticks. The nominal thickness of a single layer was set in the software at 0.04 mm. As we see in the mesh model, the layers thickness is not regular. The true thickness of the layers is determined by the geometry of the CAD model (especially narrow areas). Their shape and size are also influenced by the other discretization elements. By adding external boundary layers, it will be possible to obtain more accurate results (better accuracy of calculations); for example, better quality of the temperature distribution result on the workpiece outer surface. Figure 5a,b presents the front and back views of sticks. Figure 5d,e presents the cross section through the mesh model of a stick. Two external boundary layers are visible. Figure 5c,f shows the distribution runners.
In the next step, the process parameters based on our previous works were established [1,25]. The process parameters (mold temperature, melt temperature, and injection pressure) are presented in Table 4.

3. Results

Figure 6a shows fifty percent of filling time (about 0.6 s). Until this time, all runner systems were filled. Liquid material was starting to flow to the sticks. Figure 6b presents the total fill time (about 1.2 s). After that time, the entire volume of the part was filled. The results of the fill time analysis indicate that the entire capacity of the sticks was filled evenly.
Figure 6c shows the pressure at the end of fill. Pressure falls evenly from the maximum values closer to the injection point to zero at the end sticks. Figure 6d shows the pressure at the end of packing. The greatest pressure drop is visible after material flows from distributing runners to sticks (through small gates). It is important to remember the PIR (pressure intensification ratio). The value of injection pressure from the analysis could be higher than the real machine pressure during the real injection process. The PIR value is the relationship between the injection pressure and the hydraulic pressure (oil pressure in the piston of the injection molding machine). In practice, the pressure at the machine could be much lower than the injection pressure. A typical range of PRI value is between 7 and 15. Therefore, during the injection process, pressure should be set experimentally depending on the visual quality assessment of the part.
Figure 7a shows the temperature at the flow front. The highest value of temperature (170 °C) is visible near the injection point. This value is assumed as the flow temperature in the injection process. Temperature drops further from the injection point. This is because of the heat exchange between liquid plastic and mold at an ambient temperature (22 °C). The even temperature distribution in the process is noticeable. Each stick is filled at the same temperature. It is important for the high quality of the finished part and to maintain consistent properties of each stick.
Figure 7b shows the temperature at the end of the fill. The temperature distribution is irregular. Some regions of the part have reached the solidification temperature. The remaining places are still hot (plastic is in a liquid state). After this stage, material flow is still possible. Figure 7c shows the temperature at the end of the packing. All boundary areas of sticks reached an ambient temperature (22 °C). Most of runner’s areas did so too, except for a few regions (wider gates of main and distribution runners), where the temperature is higher. However, the entire outside boundary surface is in a solid state. Figure 7d shows the temperature at the end of the packing inside of the main runner. Only the outer walls are in contact with the cold surfaces of the mold. Inside the part, heat transfer is slow, so the core is still very hot (above 120 °C).
Figure 8 shows the temperature at the end of the packing inside of the part. The temperature inside the sticks is still high (about 70–90 °C). After this stage, there is an ejection from the mold. The interior of the part cools in the open air.
Figure 9a shows the volumetric shrinkage of the part at the end of the fill. The results of the analysis indicate that the shrinkage at the end of fill in most areas is large. Especially at the end of sticks, it is about 10 percent. During the filling process, the material cools down. Under these conditions (temperature drop, reaching the solid form), volumetric shrinkage takes place. It must be compensated by additional material from the piston of the injection molding machine (pack stage). Figure 9b presents volumetric shrinkage at the end of packing. At the sticks, shrinkage does not exceed 2 percent.
Figure 10 shows places where defects, such as air bubbles and empty cavities, may occur. The most risky regions are located at both ends of the sticks (internal and external joint areas).
In order to verify the simulation process, a Babyplast 6/10P injection molding machine was used to make PCL filaments in the form of sticks using the mold and parameters used previously for simulation. The shape of the mold was properly reproduced and good quality sticks were obtained (Figure 11).

4. Discussion

Achieving high-quality parts is a crucial task for the injection molding of sticks for further 3D printing of implants. With computer simulation, potential defects, such as air trap, weld line positions, and sink mark, emerging during the injection molding process can be predicted, thus mistakes in the design can be modified without excessive expenses of time and money [4]. With thermoplastics injection molding simulation, the material data are generally available from the simulation tools. Users could easily choose thermoplastics material data to simulate entire phases of the thermoplastics injection molding process [26]. However, polycaprolactone data are not available in the data banks of common injection molding simulation programs. Therefore, in this study, we created polycaprolactone data from the literature. The basic material properties are listed in Table 1. The selected materials properties were implemented in Solidworks Plastic software to simulate the process of injection molding. We designed the mold and evaluated the influence of the injection time, pressure, and temperature on the geometry of the part. Our studies demonstrated that the selected parameters had a good effect on the quality of the part (sticks). The mold filling process was correct. The material temperature was almost the same in the entire sticks area. The entire capacity of sticks was filled evenly and shrinkage did not exceed 2 percent. In our previous work, we started to develop the idea of filament stick production by injection molding; however, the parameters of the process were selected using the traditional trial and error method [1]. Then, computer-aided engineering was used to verify the part and mold construction [25]. The simulation of the injection process showed the faults and imperfections of the former injection mold and sticks’ design. The sticks were arranged in series; therefore, the flow path of the material to those more distant from the injection site was much longer and the last stick was filled with a lower pressure than the others in the series. The problem of temperature and pressure drop was caused by the narrow gates near the sticks. The negative impact of their small cross section was also visible in the simulation of the injection time. After passing through them, the material propagated much slower. The changes introduced in the newly designed sticks and the design of a new mold as a result of the simulation presented in this work allowed us to obtain a new quality of sticks.

5. Conclusions

Our studies indicate that SolidWorks software can be used as a simulation tool for investigating the stick quality. The parameters selected during the simulation were applied in a real process using the injection molding machine. The good results of the simulation were confirmed by the good quality of the obtained stick. Comparisons of the results of simulation and experimental methods show that SolidWorks offers similar results. Our study has shown the most significant parameters of the injection molding process, such as the injection time, pressure, and temperature. By combining the experimental results and simulation methods, it is possible to obtain a qualitative prediction of the morphological characteristics of the polymer obtained with various additives. Therefore, the promising results of these preliminary studies encourage us to continue working on the modification of PCL. Different additives (bioactive powders, antibacterial particles, and drugs) will be used in the future to study the effect of the mentioned process parameters on the modified PCL polymer. Moreover, other simulation software will be used to perform a verification of the obtained results.

6. Patents

Application number: P.428429; Filing date: 31 December 2018. Patent number: Pat.240243. Date of granting the right: 12 September 2021. Authors: Izabella Rajzer, Adam Jabłoński, Anna Kurowska, Jerzy Kopeć, and Marcin Sidzina. “A method of producing implants from bioresorbable thermoplastic polymer composites, especially in the form of three-dimensional scaffolds, intended for the reconstruction of cartilage and bone tissue defects, using 3D printing”.

Author Contributions

Conceptualization, I.R. and K.F.; methodology, K.F. and I.R.; validation, K.F. and J.J.; formal analysis, I.R. and K.F.; investigation, K.F., J.J., A.K. and P.S.; resources, I.R.; data curation, K.F. and I.R.; writing—original draft preparation, I.R. and K.F.; writing—review and editing, I.R., K.F., A.K., J.J. and P.S.; visualization, K.F.; supervision, I.R.; project administration, I.R.; funding acquisition, I.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the National Science Centre, Poland, in the frame of project “3D and 4D printing of stimuli-responsive and functionally graded biomaterials for osteochondral defects regeneration”, grant number: 2020/39/I/ST5/00569 (OPUS-LAP) and by The National Centre for Research and Development, Poland, grant number TANGO-IV-A/0002/2019.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Rajzer, I.; Kurowska, A.; Jabłoński, A.; Kwiatkowski, R.; Piekarczyk, W.; Hajduga, M.B.; Kopeć, J.; Sidzina, M.; Menaszek, E. Scaffolds modified with graphene as future implants for nasal cartilage. J. Mater. Sci. 2020, 55, 4030–4042. [Google Scholar] [CrossRef]
  2. Fuenmayor, E.; O’Donnell, C.; Gately, N.; Doran, P.; Devine, D.M.; Lyons, J.G.; McConville, C.; Major, I. Mass-customization of oral tablets via the combination of 3D printing and injection molding. Int. J. Pharm. 2019, 569, 118611. [Google Scholar] [CrossRef] [PubMed]
  3. Yizong, T.; Ariff, Z.M.; Khalil, A.M. Influence of processing parameters on injection molded polystyrene using Taguchi method as design of experiment. Procedia Eng. 2017, 184, 350–359. [Google Scholar] [CrossRef]
  4. Heidari, B.S.; Oliaei, E.; Shayesteh, H.; Davachi, S.M.; Hejazi, I.; Seyfi, J.; Bahrami, M.; Rashedi, H. Simulation of mechanical behavior and optimization of simulated injection molding process for PLA based antibacterial composite and nanocomposite bone screws using central composite design. J. Mech. Behav. Biomed. Mater. 2017, 65, 160–176. [Google Scholar] [CrossRef] [PubMed]
  5. Surace, R.; Bellantone, V.; Trotta, G.; Fassi, I. Replicating capability investigation of micro features in injection moulding process. J. Manuf. Process. 2017, 28, 351–361. [Google Scholar] [CrossRef]
  6. Tosello, G.; Costa, F.S. High precision validation of micro injection molding process simulations. J. Manuf. Process. 2019, 48, 236–248. [Google Scholar] [CrossRef]
  7. Rajzer, I.; Rom, M.; Menaszek, E.; Fabia, J.; Kwiatkowski, R. Conductive polyaniline patterns on electrospun polycaprolactone/hydroxyapatite scaffolds for bone tissue engineering. Materials 2021, 14, 4837. [Google Scholar] [CrossRef] [PubMed]
  8. Hajduga, M.B.; Bobiński, R.; Dutka, M.; Ulman-Włodarz, I.; Bujok, J.; Pająk, C.; Ćwiertnia, M.; Kurowska, A.; Dziadek, M.; Rajzer, I. Analysis of the antibacterial properties of polycaprolactone modified with graphene, bioglass and zinc doped bioglass. Acta. Bioeng. Biomech. 2021, 23, 131–138. [Google Scholar] [CrossRef] [PubMed]
  9. Woodruff, M.A.; Hutmacher, D.W. The return of a forgotten polymer—Polycaprolactone in the 21st century. Prog. Polym. Sci. 2010, 35, 1217–1256. [Google Scholar] [CrossRef] [Green Version]
  10. Kelly, C.A.; Murphy, S.H.; Leeke, G.A.; Howdle, S.M.; Shakesheff, K.M.; Jenkins, M.J. Rheological studies of polycaprolactone in supercritical CO2. Eur. Polym. J. 2013, 49, 464–470. [Google Scholar] [CrossRef]
  11. Domalik-Pyzik, P.; Morawska-Chochół, A.; Chłopek, J.; Rajzer, I.; Wrona, A.; Menaszek, E.; Ambroziak, M. Polylactide/polycaprolactone asymmetric membranes for guided bone regeneration. E-Polymers 2016, 16, 351–358. [Google Scholar] [CrossRef]
  12. Rajzer, I.; Menaszek, E.; Kwiatkowski, R.; Planell, J.A.; Castano, O. Electrospun gelatin/poly(ε-caprolactone) fibrous scaffold modified with calcium phosphate for bone tissue engineering. Mater. Sci. Eng. C 2014, 44, 183–190. [Google Scholar] [CrossRef] [PubMed]
  13. Koch, F.; Thaden, O.; Conrad, S.; Tröndle, K.; Finkenzeller, G.; Zengerle, R.; Kartmann, S.; Zimmermann, S.; Koltay, P. Mechanical properties of polycaprolactone (PCL) scaffolds for hybrid 3D-bioprinting with alginate-gelatin hydrogel. J. Mech. Behav. Biomed. Mater. 2022, 130, 105219. [Google Scholar] [CrossRef] [PubMed]
  14. Noroozi, N.; Thomson, J.A.; Noroozi, N.; Schafer, L.L.; Hatzikiriakos, S.G. Viscoelastic behaviour and flow instabilities of biodegradable polycaprolactone polyesters. Rheol. Acta. 2012, 51, 179–192. [Google Scholar] [CrossRef]
  15. Schaible, S.; Bernet, C.; Ledergeber, P.; Balmer, T.; Brändli, C. Basic study on the evaluation of thermoplastic polymers as hot-melt adhesives for mixed-substrate joining. Open J. Appl. Sci. 2016, 6, 579–592. [Google Scholar] [CrossRef] [Green Version]
  16. Wunderlich, B. Thermal Analysis of Polymeric Materials; Springer: Berlin/Heidelberg, Germany, 2005; p. 381. [Google Scholar]
  17. Gumede, T.P.; Luyt, A.S.; Hassan, M.K.; Pérez-Camargo, R.A.; Tercjak, A.; Müller, A.J. Morphology, nucleation, and isothermal crystallization kinetics of poly(ε-caprolactone) mixed with a polycarbonate/MWCNTs masterbatch. Polymers 2017, 9, 709. [Google Scholar] [CrossRef] [Green Version]
  18. Properties & Processing of CAPA® Thermoplastics. Available online: http://www.rapstrap.com/TDS-CAPA6500.pdf (accessed on 26 September 2022).
  19. Density of Plastics: Technical Properties. Available online: https://omnexus.specialchem.com/polymer-properties/properties/density (accessed on 26 September 2022).
  20. Lu, L.; Zhang, Q.; Wootton, D.M.; Chiou, R.; Li, D.; Lu, B.; Lelkes, P.I.; Zhou, J. Mechanical study of polycaprolactone-hydroxyapatite porous scaffolds created by porogen-based solid freeform fabrication method. J. Appl. Biomater. Funct Mater. 2014, 12, 145–154. [Google Scholar] [CrossRef] [PubMed]
  21. Pionteck, J. Determination of pressure dependence of polymer phase transitions by pVT analysis. Polymers 2018, 10, 578. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Oliaei, E.; Heidari, B.S.; Davachi, S.M.; Bahrami, M.; Davoodi, S.; Hejaz, I.; Seyfi, J. Warpage and shrinkage optimization of injection-molded plastic spoon parts for biodegradable polymers using Taguchi, ANOVA and artificial neural network methods. J. Mater. Sci. Technol. 2016, 32, 710–720. [Google Scholar] [CrossRef]
  23. Wang, X.; Li, H.; Gu, J.; Li, Z.; Ruan, S.; Shen, C.; Wang, M. Pressure analysis of dynamic injection molding and process parameter optimization for reducing warpage of injection molded products. Polymers 2017, 9, 85. [Google Scholar] [CrossRef]
  24. Wang, J.; Hopmann, C.; Schmitz, M.; Hohlweck, T.; Wipperfürth, J. Modeling of pvT behavior of semi-crystalline polymer based on the two-domain Tait equation of state for injection molding. Mater. Des. 2019, 183, 108149. [Google Scholar] [CrossRef]
  25. Formas, K.; Rajzer, I.; Wróbel, I. Symulacja procesów wtrysku sztyftów używanych w druku 3D. Tworzywa Sztuczne w Przemyśle 2020, 5, 58–64. [Google Scholar]
  26. Tran, N.T.; Gehde, M. Creating material data for thermoset injection molding simulation process. Polym. Test. 2019, 73, 284–292. [Google Scholar] [CrossRef]
Figure 1. (a) Viscosity–shear rate plot (based on cross model data); (b) temperature-specific volume plot (based on Tait equation data).
Figure 1. (a) Viscosity–shear rate plot (based on cross model data); (b) temperature-specific volume plot (based on Tait equation data).
Materials 15 07295 g001
Figure 2. (a) The part model; (b) injection mold CAD data translucency view; (c) injection mold: major CAD data—upper side; (d) injection mold: major CAD data—lower side.
Figure 2. (a) The part model; (b) injection mold CAD data translucency view; (c) injection mold: major CAD data—upper side; (d) injection mold: major CAD data—lower side.
Materials 15 07295 g002
Figure 3. Injection mold, simplified CAD data: (a) lower side and (b) upper side. Injection part: (c) fully CAD data and (d) simplified CAD data—one side. 3D mesh simulation model: (e) ready to analyze (injection location is visible); (f) view from above; and (g) stick region.
Figure 3. Injection mold, simplified CAD data: (a) lower side and (b) upper side. Injection part: (c) fully CAD data and (d) simplified CAD data—one side. 3D mesh simulation model: (e) ready to analyze (injection location is visible); (f) view from above; and (g) stick region.
Materials 15 07295 g003
Figure 4. Three-dimensional (3D) mesh model of stick: (a) clipping section; (b) external joint region (boundary layers are visible); and (c) internal joint region (boundary layers are visible).
Figure 4. Three-dimensional (3D) mesh model of stick: (a) clipping section; (b) external joint region (boundary layers are visible); and (c) internal joint region (boundary layers are visible).
Materials 15 07295 g004
Figure 5. Three-dimensional (3D) mesh model of a stick: (a) front view of sticks; (b) back view of stick; and (d,e) clipping section (boundary layers are visible). Sticks and distribution runners: (c) view from above and (f) bottom view.
Figure 5. Three-dimensional (3D) mesh model of a stick: (a) front view of sticks; (b) back view of stick; and (d,e) clipping section (boundary layers are visible). Sticks and distribution runners: (c) view from above and (f) bottom view.
Materials 15 07295 g005
Figure 6. (a) Fifty percent of the fill time; (b) total fill time; (c) pressure at the end of fill; and (d) pressure at the end of packing.
Figure 6. (a) Fifty percent of the fill time; (b) total fill time; (c) pressure at the end of fill; and (d) pressure at the end of packing.
Materials 15 07295 g006
Figure 7. (a) Flow front temperature; (b) temperature at the end of fill—boundary surface; (c) temperature at the end of packing—boundary surface; and (d) temperature at the end of packing—clipping section through the runner.
Figure 7. (a) Flow front temperature; (b) temperature at the end of fill—boundary surface; (c) temperature at the end of packing—boundary surface; and (d) temperature at the end of packing—clipping section through the runner.
Materials 15 07295 g007
Figure 8. (a) Temperature at the end of packing—clipping section through the stick; (b) external joint areas; and (c) internal joint areas.
Figure 8. (a) Temperature at the end of packing—clipping section through the stick; (b) external joint areas; and (c) internal joint areas.
Materials 15 07295 g008
Figure 9. (a) Volumetric shrinkage at the end of fill; (b) volumetric shrinkage at the end of packing.
Figure 9. (a) Volumetric shrinkage at the end of fill; (b) volumetric shrinkage at the end of packing.
Materials 15 07295 g009
Figure 10. Air traps.
Figure 10. Air traps.
Materials 15 07295 g010
Figure 11. PCL stick produced using parameters developed during the simulation: (a) stick at the end of the injection molding process; (b) produced part; and (c) PCL filament sticks.
Figure 11. PCL stick produced using parameters developed during the simulation: (a) stick at the end of the injection molding process; (b) produced part; and (c) PCL filament sticks.
Materials 15 07295 g011
Table 1. Rheological, thermal, and mechanical properties of polycaprolactone needed to build the model, collected from the literature.
Table 1. Rheological, thermal, and mechanical properties of polycaprolactone needed to build the model, collected from the literature.
ParameterValueLiterature
Rheological properties. Viscosity model: Cross-WLF
D2 [K]243.15[10,14,15]
n0.1547
D1 [Pa∙s]2.21915 × 1012
τ* [Pa]429,551
A129.20632
A2 [K]51.6
Physical Properties
Solid density, kg/m31140[19]
Thermal Properties
Glass transition, °C−60[16,17,18]
Specific heat, [J/kgK]2100
Melt flow rate (MFR),
g/10 min (190 °C/2.16 kg)
28
Thermal conductivity, [W/m·K] (20 °C)0.1
Melting point [°C]60–62
Mechanical properties
Yield stress, σy, [MPa]17.5[20]
Tensile modulus [MPa]430[18]
Poisson’s ratio0.44[20]
Table 2. Pressure–volume–temperature properties [21].
Table 2. Pressure–volume–temperature properties [21].
P–V–T (Pressure–Volume–Temperature) Properties
13-Parameter Modified Tait Equation
ParameterValue
b1m [m3/kg]9.463 × 10−4
b2m [m3/(kg-K)]6.514 × 10−7
b3m [Pa]2.030180 × 108
b4m [1/K]4.521 × 10−3
b1s [m3/kg]9.045 × 10−4
b2s [m3/(kg-K)]6.291 × 10−7
b3s [Pa]2.52593 × 108
b4s [1/K]7.368 × 10−3
b5 [K]336.77
b6 [K/Pa]1.497 × 10−8
b7 [m3/kg]4.181 × 10−5
b8 [1/K]4.042 × 10−1
b9 [1/Pa]7.393 × 10−9
Table 3. Discretization details.
Table 3. Discretization details.
Mesh Details
solver typeAMG (algebraic multigrid)
mesh typesolid mesh
mesh volume, [mm3]955
mesh typetetrahedral hybrid
total elements of surface mesh59,756
total nodes of surface mesh29,857
total elements of solid mesh394,602
total nodes of solid mesh119,344
analyze typefill and pack
Table 4. Manufacturing parameters.
Table 4. Manufacturing parameters.
Process Details
Process ParameterValue
ambient temperature, [°C]22
melt temperature, [°C]170
mold temperature, [°C]22
filling time, [sec] 1
injection pressure at the nozzle (filling), [MPa]130
packing time, [sec]2.5
injection pressure at the nozzle (packing), [MPa]20
cooling time, [sec]25
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Formas, K.; Kurowska, A.; Janusz, J.; Szczygieł, P.; Rajzer, I. Injection Molding Process Simulation of Polycaprolactone Sticks for Further 3D Printing of Medical Implants. Materials 2022, 15, 7295. https://doi.org/10.3390/ma15207295

AMA Style

Formas K, Kurowska A, Janusz J, Szczygieł P, Rajzer I. Injection Molding Process Simulation of Polycaprolactone Sticks for Further 3D Printing of Medical Implants. Materials. 2022; 15(20):7295. https://doi.org/10.3390/ma15207295

Chicago/Turabian Style

Formas, Krzysztof, Anna Kurowska, Jarosław Janusz, Piotr Szczygieł, and Izabella Rajzer. 2022. "Injection Molding Process Simulation of Polycaprolactone Sticks for Further 3D Printing of Medical Implants" Materials 15, no. 20: 7295. https://doi.org/10.3390/ma15207295

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