Inverse Determination of Johnson–Cook Parameters of Additively Produced Anisotropic Maraging Steel

In powder bed-based additive manufacturing (AM), complex geometries can be produced in a layer-wise approach. Results of material science experiments regarding material property identification, e.g., tensile strength, show interdependencies between the test load direction and the layer orientation. This goes hand-in-hand with the measured cutting force, changing with the relative angle between cutting direction and layer orientation in orthogonal cutting tests. However, due to the specific process characteristics, the layer orientation results in anisotropic material properties. Therefore, during machining, the material behaves depending on the buildup direction, which influences the cutting process. To predict this behavior, a simplified inverse approach is developed to determine the buildup direction-dependent parameters of a modified Johnson–Cook model for cutting simulation. To qualify these cutting models, mainly the cutting force and additionally the chip formation examined during orthogonal cuts are used. In the present paper, the influence of the laser-powder-bed-fusion (LPBF) process parameters on subtractive post-processing are shown. A good agreement between verification experiments and simulations is achieved.


Introduction
The LPBF process, as used for this research, utilizes lasers as a heat source to selectively melt powder in a thin layer. Repetition of this process over multiple layers can create complex three-dimensional structures. The layer-wise build-up material also leads to a rough surface that needs to be machined for most applications. For this machining, processes with a defined cutting edge are used in particular [1]. However, due to the selectively induced heat and large cooldown rates, a fine dendritic microstructure forms along the build-up direction [2]. This leads to an anisotropic material behavior, which also affects the cutting process [3].
The objective of machining simulation is to enable machining process optimization with a low expense of material-and time-intensive machining examinations. Suitable material models for machining simulation are often not available, especially for new materials. This also includes additively manufactured materials whose material properties are highly dependent on the manufacturing process. The aim of the work presented here is to model the influence of the build-up direction on chip formation in machining simulations and to provide adapted Johnson-Cook parameters for this purpose. The anisotropic behavior of specimens made from various steel and stainless steel alloys has been presented by numerous researchers. They all state a difference between build orientations. An almost linear decrease in ultimate tensile strength with an increase of build orientation ϕ is reported [4][5][6].
To take this influence into account in the machining simulation, inverse parameter identification can be used to determine the characteristic values of a material based on simulation calculations for a specific loading condition. In this way, the comparatively high strain rates and steep temperature gradients in machining processes can also be taken into account, which can only be mapped with great effort using experimental analyses. Simulation of machining processes is gaining more and more importance and is used both to get knowledge about basic machining mechanisms, as well as to optimize specific machining tasks [7,8]. The influence of a varied tool shape or of variable machining parameters on chip formation are exemplary subjects of intensive research work in recent years [9][10][11]. Furthermore, the machining temperature, residual stress states, and influences on microstructure formation were investigated by means of simulation [12][13][14]. In addition to the consideration of the machining process itself, the method of machining simulation has also been successively developed. The focus was, for example, on processes for remeshing and mesh refinement in the shear zone and approaches for the simulation of specific chip formation phenomena, such as shear chip formation [9,15]. However, the research topic was also the modeling of chip formation depending on the material as a function of its properties. Starting with differently conditioned low alloy steel [16], non-ferrous metals, such as titanium alloy Ti-6Al-4V [17] or nickel-based superalloys, such as Inconel 718 [18].
For the simulations to be conducted, the required material parameters, for instance Johnson-Cook parameters, usually are taken from the literature. There are mainly two approaches to determine the constants of the Johnson-Cook constitutive equation, which are used in the modeling of machining processes. The first approach is based on standardized experimental material tests, such as tensile and compression tests, or methods that reproduce the general conditions of machining processes, such as the split Hopkinson pressure bar (SHPB) test [19][20][21]. The searched individual Johnson-Cook constants are derived from the resulting flow curves using well-known mathematical methods [22]. A second approach, which is gaining in importance, uses the simulation of machining processes to get adapted material parameters [2,[23][24][25][26][27][28][29]. Here, the material parameters are modified within chip formation simulations starting from available initial values while observing the influence on previously defined evaluation criteria. In this way, for example, the Johnson-Cook parameters are varied repeatedly within an increasingly decreasing range until the simulated criteria deviate only minimally from the experimentally determined criteria. The required initial parameters are taken from the literature or determined on the basis of material tests. Based on this approach, the orientation-dependent properties of additively produced M300 steel and their influence on chip formation in the orthogonal cut are to be used for the simulation-based determination of the Johnson-Cook parameters and validated based on experimentally determined comparison criteria. Here, both the cutting force and chip shape criteria are integrated.

Material
The here described inverse parameter identification was conducted on the example of additively manufactured specimens made of grade 300 maraging steel (short M300), according to ASTM A646. Manufacturers characterize the martensitic hot work tool steel as a material with both high hardness and high ductility. These properties enable the material with a high wear resistance and thus allow applications in mechanically and thermally highly stressed components, e.g., forming tools and molds. Additive manufacturing is a suitable method to create complex internal contours, such as channels for cooling fluids, in components like these. Table 1 shows the chemical composition of the material M300. Due to the increasing use of M300, obtaining detailed knowledge of both material properties [30][31][32] and behavior [33,34] under various conditions is a focus of ongoing research. The Johnson-Cook constitutive equation [36] describing flow curves using few parameters in an easy-to-use manner is the basis of one of the most widely used material models. Equation (1) describes the von Mises flow stress with three terms considering different influences: the strain hardening, the strain rate, and the temperature.
where ε Equivalent plastic strain Dimensionless plastic strain rate for Influeence of strain rate on stress and choosing Equation (1) can be written in the generalized form: In the literature, Johnson-Cook parameters can be found for the most widely used materials. Usually, those values were determined experimentally with the known methods, e.g., based on tensile and compression tests, in some cases conducted under the thermal and mechanical conditions of machining processes. Unfortunately, it is challenging to find material models especially describing the behavior depending on the build-up orientation of additively manufactured parts. Table 2 shows the Johnson-Cook parameters of the above-mentioned maraging steel M300 used as material within this work. The data were found in the literature and has been elaborated from stress-strain curves for aged maraging steel [37].

Additive Manufacturing Process
The LPBF process was used within the described research work to create material specimens. It can be used to process a wide variety of metallic materials, such as aluminum alloys, steels, titanium alloys, and Ni-base super alloys. In theory, all weldable metallic materials are usable by means of the LPBF process [38]. The material specimens used here were produced on a RenAM 500Q LPBF system from Renishaw plc., Wotton-under-Edge, UK. It is equipped with four laser beam outputs, each with 500 W optical power, synchronously useable via 2-mirror galvanometers. The system has a build chamber with a built-up volume of 250 × 250 × 350 mm using a closed vacuum-assisted system for initial inertization and shielding gas recirculation (Renishaw plc., Wotton-under-Edge, UK). The specimens were produced with the inert gas nitrogen at a flow rate of 190 m 3 /h. The base plate was heated to 60 • C to reduce thermal influences on the support structures by gradually heating the build plate during the build-up. The layer height was set to 50 µm. The number of borders was set to two. The hatch distance was set to 75 µm. The hatch pattern was set to "stripes" with a length of 10 mm. The laser power was set to 150 W for each output for the borders and to 200 W for the hatching. The spot diameter of the machine was 80 µm. In this way, specimens with the shape of flat cuboids measuring 170 mm in length, 5 mm in thickness, and 25 mm in height were produced from the M300 powder. By orienting the specimens on the build platform, the build direction of the workpiece specimens was varied at angles of 0 and 90 degrees ( Figure 1). The stated fine dendritic microstructure, which forms along the build-up direction, can be seen in Figure 1b. It is known that this layer-wise structure anisotropic material behavior affects the cutting process [3]. Edge, UK. It is equipped with four laser beam outputs, each with 500 W optical power, synchronously useable via 2-mirror galvanometers. The system has a build chamber with a built-up volume of 250 × 250 × 350 mm using a closed vacuum-assisted system for initial inertization and shielding gas recirculation (Renishaw plc., Wotton-under-Edge, UK). The specimens were produced with the inert gas nitrogen at a flow rate of 190 m 3 /h. The base plate was heated to 60 °C to reduce thermal influences on the support structures by gradually heating the build plate during the build-up. The layer height was set to 50 μm. The number of borders was set to two. The hatch distance was set to 75 μm. The hatch pattern was set to "stripes" with a length of 10 mm. The laser power was set to 150 W for each output for the borders and to 200 W for the hatching. The spot diameter of the machine was 80 μm. In this way, specimens with the shape of flat cuboids measuring 170 mm in length, 5 mm in thickness, and 25 mm in height were produced from the M300 powder. By orienting the specimens on the build platform, the build direction of the workpiece specimens was varied at angles of 0 and 90 degrees ( Figure 1). The stated fine dendritic microstructure, which forms along the build-up direction, can be seen in Figure 1b. It is known that this layer-wise structure anisotropic material behavior affects the cutting process [3].

Machining Process
The specimens were used in experiments on orthogonal cutting on a test setup that enables a straight tool movement. The cuts were made both parallel (φ = 0 degrees) and normal (φ = 90 degrees) to the layer orientation regarding the long side of the specimen. Figure 2 shows the overall view of the test stand (left-hand side), as well as an exemplary specimen mounted on a dynamometer (Kistler Instrumente AG, Winterthur, Switzerland), an upstream run-in part, and the tool. The used tool insert is a parting and grooving tool type S229.0800.36 manufactured by Paul Horn GmbH, Tübingen, Germany. The inserts are made from a proprietary tungsten carbide compound called Ti 25. Table 3 lists essential tool characteristics. The experiments were conducted with a feed rate of f = 0.1 mm and a cutting speed of vc = 80 m/min.

Machining Process
The specimens were used in experiments on orthogonal cutting on a test setup that enables a straight tool movement. The cuts were made both parallel (ϕ = 0 degrees) and normal (ϕ = 90 degrees) to the layer orientation regarding the long side of the specimen. Figure 2 shows the overall view of the test stand (left-hand side), as well as an exemplary specimen mounted on a dynamometer (Kistler Instrumente AG, Winterthur, Switzerland), an upstream run-in part, and the tool. Edge, UK. It is equipped with four laser beam outputs, each with 500 W optical power, synchronously useable via 2-mirror galvanometers. The system has a build chamber with a built-up volume of 250 × 250 × 350 mm using a closed vacuum-assisted system for initial inertization and shielding gas recirculation (Renishaw plc., Wotton-under-Edge, UK). The specimens were produced with the inert gas nitrogen at a flow rate of 190 m 3 /h. The base plate was heated to 60 °C to reduce thermal influences on the support structures by gradually heating the build plate during the build-up. The layer height was set to 50 μm. The number of borders was set to two. The hatch distance was set to 75 μm. The hatch pattern was set to "stripes" with a length of 10 mm. The laser power was set to 150 W for each output for the borders and to 200 W for the hatching. The spot diameter of the machine was 80 μm. In this way, specimens with the shape of flat cuboids measuring 170 mm in length, 5 mm in thickness, and 25 mm in height were produced from the M300 powder. By orienting the specimens on the build platform, the build direction of the workpiece specimens was varied at angles of 0 and 90 degrees ( Figure 1). The stated fine dendritic microstructure, which forms along the build-up direction, can be seen in Figure 1b. It is known that this layer-wise structure anisotropic material behavior affects the cutting process [3].

Machining Process
The specimens were used in experiments on orthogonal cutting on a test setup that enables a straight tool movement. The cuts were made both parallel (φ = 0 degrees) and normal (φ = 90 degrees) to the layer orientation regarding the long side of the specimen. Figure 2 shows the overall view of the test stand (left-hand side), as well as an exemplary specimen mounted on a dynamometer (Kistler Instrumente AG, Winterthur, Switzerland), an upstream run-in part, and the tool. The used tool insert is a parting and grooving tool type S229.0800.36 manufactured by Paul Horn GmbH, Tübingen, Germany. The inserts are made from a proprietary tungsten carbide compound called Ti 25. Table 3 lists essential tool characteristics. The experiments were conducted with a feed rate of f = 0.1 mm and a cutting speed of vc = 80 m/min. The used tool insert is a parting and grooving tool type S229.0800.36 manufactured by Paul Horn GmbH, Tübingen, Germany. The inserts are made from a proprietary tungsten carbide compound called Ti 25. Table 3 lists essential tool characteristics. The experiments were conducted with a feed rate of f = 0.1 mm and a cutting speed of v c = 80 m/min. For statistical validation, 30 replicate tests were conducted for both material orientations. Samples of the resulting chips were examined microscopically, and both the chip thickness t and the lamella distance d were determined in order to correlate them with the builtup angle ϕ. Figure 3 shows both the resulting mean value of the cutting force F c of all individual tests, depending on the built-up orientation angle ϕ and sections of typical chips with a lamellar appearance. Mean values of the cutting forces were derived from the measured cutting force curves. The amplitude of the dynamic signal component is relatively low at less than 10% of the mean value. Due to the limited natural frequency of the dynamometer, it is not possible to derive any reliable quantified conclusions about the relationship between the dynamic component of the force curve and the lamella formation. The cutting force F c_0 = 1157.2 N in the case of ϕ = 0 • and F c_90 = 1069.7 N with ϕ = 90 • differ by 7.56%. This result corresponds with the above-mentioned context between tensile strength and built-up orientation ϕ. For statistical validation, 30 replicate tests were conducted for both material orientations. Samples of the resulting chips were examined microscopically, and both the chip thickness t and the lamella distance d were determined in order to correlate them with the builtup angle φ. Figure 3 shows both the resulting mean value of the cutting force Fc of all individual tests, depending on the built-up orientation angle φ and sections of typical chips with a lamellar appearance. Mean values of the cutting forces were derived from the measured cutting force curves. The amplitude of the dynamic signal component is relatively low at less than 10% of the mean value. Due to the limited natural frequency of the dynamometer, it is not possible to derive any reliable quantified conclusions about the relationship between the dynamic component of the force curve and the lamella formation. The cutting force Fc_0 = 1157.2 N in the case of φ = 0° and Fc_90 = 1069.7 N with φ = 90° differ by 7.56%. This result corresponds with the above-mentioned context between tensile strength and built-up orientation φ.

Simulation
To obtain adapted Johnson-Cook parameters in regard to different built-up directions of the additively manufactured material, inverse parameter identification was implemented by means of the finite-element method. The goal of this procedure is to compare the results of chip formation simulations and machining experiments based on predefined criteria. In the iterative approach chosen here, a defined number of parameter sets was varied within progressively narrower limits and used as input variables for the simulations. In doing so, parameter sets were identified for which the simulation results approximate the experimentally determined criteria as closely as possible. The limits of the next iteration step result from those two parameter sets with the best approximation. The starting parameter range of the first iteration step was created by means of a DoE on the basis of the Latin hyper-cube method. Due to the approach commonly described in the literature, the cutting force Fc was taken into account in the first simulation runs [39,40]. The two-dimensional simulation model was developed using DEFORM 3D™ V12.01 (V12.01, Scientific Forming Technologies Corp., Columbus, OH, USA). The shape of the previously mentioned turning tool has been taken into account by the clearance angle α = 12° and the rake angle γ = 15°, as well as the cutting-edge radius rb = 47.6 μm. The part's

Simulation
To obtain adapted Johnson-Cook parameters in regard to different built-up directions of the additively manufactured material, inverse parameter identification was implemented by means of the finite-element method. The goal of this procedure is to compare the results of chip formation simulations and machining experiments based on predefined criteria. In the iterative approach chosen here, a defined number of parameter sets was varied within progressively narrower limits and used as input variables for the simulations. In doing so, parameter sets were identified for which the simulation results approximate the experimentally determined criteria as closely as possible. The limits of the next iteration step result from those two parameter sets with the best approximation. The starting parameter range of the first iteration step was created by means of a DoE on the basis of the Latin hyper-cube method. Due to the approach commonly described in the literature, the cutting force F c was taken into account in the first simulation runs [39,40]. The two-dimensional simulation model was developed using DEFORM 3D™ V12.01 (V12.01, Scientific Forming Technologies Corp., Columbus, OH, USA). The shape of the previously mentioned turning tool has been taken into account by the clearance angle α = 12 • and the rake angle γ = 15 • , as well as the cutting-edge radius r b = 47.6 µm. The part's mesh is generated with local refining and a maximum mesh number of 10,000. Focusing on the Johnson-Cook parameters responsible for the strain hardening effect, sets with varied parameters A, B and n based on Table 2 were prepared in a design-of-experiment approach, which is based on the Latin hypercube method. The upper and lower limits of the parameter ranges shown in Table 4 were determined to be +20% and −20% of the values described in the literature. The parameters C and m remained unchanged during the simulation phase. In a first iteration, a number of 100 individual parameter sets were determined within the limits. To achieve an even cutting force progression, the initial simulations were conducted without any damage criterion resulting in the formation of flow chips. This simplification was also made with regard to the cutting force measured and subsequently averaged in the experimental investigations and the fact that the dynamic force component cannot be unambiguously quantified and thus used as a comparison criterion. The used simulation software prefers the generalized form of the Johnson-Cook material model. Figure 4 shows the 2D simulation model with the meshed part and tool, as well as the resulting chip shape using the start parameters.
Materials 2022, 14, x FOR PEER REVIEW 6 of 12 mesh is generated with local refining and a maximum mesh number of 10,000. Focusing on the Johnson-Cook parameters responsible for the strain hardening effect, sets with varied parameters A, B and n based on Table 2 were prepared in a design-of-experiment approach, which is based on the Latin hypercube method. The upper and lower limits of the parameter ranges shown in Table 4 were determined to be +20% and −20% of the values described in the literature. The parameters C and m remained unchanged during the simulation phase. In a first iteration, a number of 100 individual parameter sets were determined within the limits. To achieve an even cutting force progression, the initial simulations were conducted without any damage criterion resulting in the formation of flow chips. This simplification was also made with regard to the cutting force measured and subsequently averaged in the experimental investigations and the fact that the dynamic force component cannot be unambiguously quantified and thus used as a comparison criterion. The used simulation software prefers the generalized form of the Johnson-Cook material model. Figure 4 shows the 2D simulation model with the meshed part and tool, as well as the resulting chip shape using the start parameters.  Figure 5 shows the correlation between the cutting force value and the number of the individual parameter sets resulting from the first iteration. The experimentally determined mean values of the cutting forces based on the two built-up orientations, φ = 0° and 90°, are represented by a red (Fc_0 = 1157 N) and a blue (Fc_90 = 1067 N) line. Two parameter sets with the best cutting force approximation for each material orientation (set 0_1 and set 0_2, set 90_1 and 90_2) were determined.  Figure 5 shows the correlation between the cutting force value and the number of the individual parameter sets resulting from the first iteration. The experimentally determined mean values of the cutting forces based on the two built-up orientations, ϕ = 0 • and 90 • , are represented by a red (F c_0 = 1157 N) and a blue (F c_90 = 1067 N) line. Two parameter sets with the best cutting force approximation for each material orientation (set 0_1 and set 0_2, set 90_1 and 90_2) were determined. These sets are used as upper and lower limits for a second identification iteration to get closer to the experimentally gained cutting force values. The nearest approximation can be reached with the parameter sets #9 and #52 in the case of the orientation angle φ = 0° and with the sets #7 and #64 in the case of φ = 0°. These parameter sets (Table 5) were used as new limits for the second identification iteration based on a DOE with 25 runs each. Figures 6 and 7 show the calculated cutting forces depending on the individual sets. A renewed search for the parameter sets with the closest approximation to the experimentally gained cutting force values leads to sets #15 and #24 in the case of the orientation angle φ = 0° and sets #3 and #4 in the case of φ = 90°. The belonging Johnson-Cook parameters A, B and n are shown in Table 6.   These sets are used as upper and lower limits for a second identification iteration to get closer to the experimentally gained cutting force values. The nearest approximation can be reached with the parameter sets #9 and #52 in the case of the orientation angle ϕ = 0 • and with the sets #7 and #64 in the case of ϕ = 0 • . These parameter sets (Table 5) were used as new limits for the second identification iteration based on a DOE with 25 runs each. Figures 6 and 7 show the calculated cutting forces depending on the individual sets. A renewed search for the parameter sets with the closest approximation to the experimentally gained cutting force values leads to sets #15 and #24 in the case of the orientation angle ϕ = 0 • and sets #3 and #4 in the case of ϕ = 90 • . The belonging Johnson-Cook parameters A, B and n are shown in Table 6.  These sets are used as upper and lower limits for a second identification iteration to get closer to the experimentally gained cutting force values. The nearest approximation can be reached with the parameter sets #9 and #52 in the case of the orientation angle φ = 0° and with the sets #7 and #64 in the case of φ = 0°. These parameter sets (Table 5) were used as new limits for the second identification iteration based on a DOE with 25 runs each. Figures 6 and 7 show the calculated cutting forces depending on the individual sets. A renewed search for the parameter sets with the closest approximation to the experimentally gained cutting force values leads to sets #15 and #24 in the case of the orientation angle φ = 0° and sets #3 and #4 in the case of φ = 90°. The belonging Johnson-Cook parameters A, B and n are shown in Table 6.     To validate these final results and to gain information about the difference between them influencing the chip formation, the chip shape was taken into account as further criteria ( Figure 8). Consequently, both the mean chip thickness t and the mean distance d between the edge of two neighboring chip lamellas was determined on experimentally produced chips, as well as on chips resulting from simulations. φ = 0°, parameter set #15 φ = 0°, parameter set #24 φ = 90°, parameter set #17 φ = 90°, parameter set #3 Regarding this, the previously mentioned FE model was extended with a fracture model according to Cockroft and Latham [41] to enable the formation of lamella chips.  To validate these final results and to gain information about the difference between them influencing the chip formation, the chip shape was taken into account as further criteria ( Figure 8). Consequently, both the mean chip thickness t and the mean distance d between the edge of two neighboring chip lamellas was determined on experimentally produced chips, as well as on chips resulting from simulations.  To validate these final results and to gain information about the difference between them influencing the chip formation, the chip shape was taken into account as further criteria ( Figure 8). Consequently, both the mean chip thickness t and the mean distance d between the edge of two neighboring chip lamellas was determined on experimentally produced chips, as well as on chips resulting from simulations. φ = 0°, parameter set #15 φ = 0°, parameter set #24 φ = 90°, parameter set #17 φ = 90°, parameter set #3 Regarding this, the previously mentioned FE model was extended with a fracture model according to Cockroft and Latham [41] to enable the formation of lamella chips. Regarding this, the previously mentioned FE model was extended with a fracture model according to Cockroft and Latham [41] to enable the formation of lamella chips. Figure 8 shows the simulation results, where the chip formation depending on the chosen parameter set could be compared. It is noteworthy that chips resulting from the build-up angle ϕ = 0 • show a more regular succession of the individual lamellas than those with ϕ = 90 • . In addition, the chip shapes are little slimmer under condition ϕ = 0 • . In the case of the direct comparison of different sets at the same built-up angle ϕ = 0 • , it is noticeable that the lamellas originating from set #15 take a larger share of the total chip width than those from set #24.

Results and Discussion
To obtain information about the correlations between chip thickness t as well as lamella distance d and the parameter sets, several measurements were conducted on chips gained in both chip formation simulations and cutting experiments. The chip thickness was determined, regarding Figure 9, always at the thinnest point between the individual lamellas. The lamella distance was measured from edge to edge.
Materials 2022, 14, x FOR PEER REVIEW 9 of 12 Figure 8 shows the simulation results, where the chip formation depending on the chosen parameter set could be compared. It is noteworthy that chips resulting from the build-up angle φ = 0° show a more regular succession of the individual lamellas than those with φ = 90°. In addition, the chip shapes are little slimmer under condition φ = 0°. In the case of the direct comparison of different sets at the same built-up angle φ = 0°, it is noticeable that the lamellas originating from set #15 take a larger share of the total chip width than those from set #24.
To obtain information about the correlations between chip thickness t as well as lamella distance d and the parameter sets, several measurements were conducted on chips gained in both chip formation simulations and cutting experiments. The chip thickness was determined, regarding Figure 9, always at the thinnest point between the individual lamellas. The lamella distance was measured from edge to edge.   Table 7 shows the relative deviations between the individual criteria of the simulated chip shapes and those of the experimentally obtained ones. The variable cutting forces determined with the damage criterion are on average on a lower level than the cutting force with pure flow chip formation (without damage criterion). The maximum values always correspond to the cutting force with flow-chip formation when the individual lamellae reach their maximum thickness. Figure 11 illustrates this for the case of the built-up angle φ = 0° for the parameter sets #15 and #24.
Overall, it can be seen that the dependence of the criteria on the built-up angle can be mapped in the trend. However, it should also be mentioned that the inclusion of further comparison criteria, in particular the dynamic cutting force component but also, for example, the machining temperature, seems advisable for further minimization of the present deviations.   Table 7 shows the relative deviations between the individual criteria of the simulated chip shapes and those of the experimentally obtained ones.
Materials 2022, 14, x FOR PEER REVIEW 9 of 12 Figure 8 shows the simulation results, where the chip formation depending on the chosen parameter set could be compared. It is noteworthy that chips resulting from the build-up angle φ = 0° show a more regular succession of the individual lamellas than those with φ = 90°. In addition, the chip shapes are little slimmer under condition φ = 0°. In the case of the direct comparison of different sets at the same built-up angle φ = 0°, it is noticeable that the lamellas originating from set #15 take a larger share of the total chip width than those from set #24.
To obtain information about the correlations between chip thickness t as well as lamella distance d and the parameter sets, several measurements were conducted on chips gained in both chip formation simulations and cutting experiments. The chip thickness was determined, regarding Figure 9, always at the thinnest point between the individual lamellas. The lamella distance was measured from edge to edge.  Figure 10 shows that both the mean values of the simulated chip thickness and the lamella distance are below the mean values of the experimentally gained chips. Table 7 shows the relative deviations between the individual criteria of the simulated chip shapes and those of the experimentally obtained ones. The variable cutting forces determined with the damage criterion are on average on a lower level than the cutting force with pure flow chip formation (without damage criterion). The maximum values always correspond to the cutting force with flow-chip formation when the individual lamellae reach their maximum thickness. Figure 11 illustrates this for the case of the built-up angle φ = 0° for the parameter sets #15 and #24.
Overall, it can be seen that the dependence of the criteria on the built-up angle can be mapped in the trend. However, it should also be mentioned that the inclusion of further comparison criteria, in particular the dynamic cutting force component but also, for example, the machining temperature, seems advisable for further minimization of the present deviations.  The variable cutting forces determined with the damage criterion are on average on a lower level than the cutting force with pure flow chip formation (without damage criterion).
The maximum values always correspond to the cutting force with flow-chip formation when the individual lamellae reach their maximum thickness. Figure 11 illustrates this for the case of the built-up angle ϕ = 0 • for the parameter sets #15 and #24.

Summary and Conclusions
The objective of this paper was to present inverse parameter identification based on FE simulations to determine the Johnson-Cook parameters A, B, and n for maraging stee M300 used for additive manufacturing in the two built-up directions.

•
For a later comparison of experimentally and simulatively obtained criteria, marag ing steel M300 specimens were prepared with both build-up directions, φ = 0° and Overall, it can be seen that the dependence of the criteria on the built-up angle can be mapped in the trend. However, it should also be mentioned that the inclusion of further comparison criteria, in particular the dynamic cutting force component but also, for example, the machining temperature, seems advisable for further minimization of the present deviations.

Summary and Conclusions
The objective of this paper was to present inverse parameter identification based on FE simulations to determine the Johnson-Cook parameters A, B, and n for maraging steel M300 used for additive manufacturing in the two built-up directions.

•
For a later comparison of experimentally and simulatively obtained criteria, maraging steel M300 specimens were prepared with both build-up directions, ϕ = 0 • and 90 • , by using the LPBF Process. In machining experiments, both the cutting force Fc and the chip shape characteristic were determined. The cutting force depends on the built-up orientation and differs by 7.56%.

•
The required 2D finite element model was prepared using DEFORM 3D™ V12.01 and used to get sets of A, B and n Johnson-Cook parameters, which lead to the closest approximation to the experimental cutting forces. Start parameters to set up a DOE were taken from the literature.

•
As additional criteria, both the chip thickness and the lamella distance of experimentally and simulative gained chips were determined and compared. While the absolute differences between the criteria of the experimentally and simulatively obtained chip properties are comparatively large, the dependence of both the chip thickness and the lamella distance on the buildup direction can be confirmed in the trend. The closest results are 26.6% difference of the chip thickness between ϕ = 0 • and 90 • in the simulation and 30.6% in the experiment. For the case of the lamella distance, the closest approximation has been reached with 9.3% in the simulation compared with 7.3% in the experiment.

•
With reference to the simplified approach used in this work, further cutting experiments with a high-resolution and high-sensing force measurement system are planned. In this context, not only the influence of the built-up angle ϕ, but also the influence of individual layers in their interaction with lamellar chip formation is to be considered. In addition, temperature measurements in the area of the cutting zone are to be included. With additional criteria derived from this, further simulations are then to be performed with regard to a more accurate inverse parameter identification.