SYSTEM AND METHOD FOR DETERMINING BEAM POWER LEVEL ALONG AN ADDITIVE DEPOSITION PATH

Abstract

A laser power control method for additive manufacturing includes a pre-processing component and an intra-processing component. The pre-processing component creates in the system controller a machine code expressed additive-path-with-geometry-metadata that includes a path description. The path description represents the path of the beam source on the build and includes a geometry index for the build. The intra-processing component calculates required power for the beam at intervals and events on the beam path based upon the additive-path-with-geometry-metadata and calculations of the energy balance at the melt pool and the total energy for each point P(s) on the path.

Claims

1. A method of controlling beam power during an additive deposition process that is controlled by a system controller and that uses a beam source to project energy to a point P(s) along an additive path that will be traveled to form a build, the build having a geometry and being formed from deposited material added to a structure, the method comprising: a pre-processing component and an intra-processing component; the pre-processing component comprising creating in the system controller a machine code expressed additive-path-with-geometry-metadata, the additive-path-with-geometry-metadata including: a) a path description, the path description representing the path of the beam source through space during the additive process; and b) a geometry index determined for the point P(s) based upon a geometric description representing the geometry of the build during the additive process, the geometry index comprising predicted thermal characteristics as to a determined melt pool, hot zone and bulk portion for the point P(s); the intra-processing component comprising at a selected time: 1) calculating an energy balance at the melt pool (EB.sub.meltpool) for the point P(s) on the additive path, the energy balance calculation being based upon the following calculations: a) a calculation of energy radiated from the melt pool (H.sub.rad); b) a calculation of energy conducted from the melt pool to the hot zone (H.sub.cond), the calculation of H.sub.cond being based upon the additive-path-with-geometry-metadata and c) a calculation of energy lost due to convection at the melt pool (H.sub.conv); 2) calculating total energy (H.sub.total) needed at the point P(s) on the additive path according to the following equation H.sub.total=H.sub.deposited material+H.sub.remelt+EB.sub.melt pool wherein H.sub.deposited material represents energy required to melt the deposited material and H.sub.remelt it represents energy required to remelt existing material; 3) calculating optimum beam source power Q.sub.source for the point P(s) according to the following equation Q.sub.source=min(Q.sub.max, H.sub.total/αΔt) wherein Q.sub.max represent maximum laser power, α represents a beam absorption coefficient and Δt represents a calculation interval between laser power computations; the selected time being the occurrence of a defined event during the additive deposition process or the beam source arriving at a selected point on the additive path during the additive deposition process; and modulating the energy projected by the beam source based upon the calculated optimum beam source power Q.sub.source.

2. The method of claim 1 wherein the geometry index is determined by a geometry index calculation, the geometry index calculation including: calculating a geometry index hot zone volume (V.sub.hot) for the point P(s) on the additive path according to the formula V.sub.hot=G∩Z(r.sub.hot, P(s), {circumflex over (t)}(s)) wherein G represents deposition geometry, Z represents a hemisphere of a certain radius (r) and having a circular surface centered at a point P(s) on the additive path with a normal direction (ň), r.sub.hot represents the radius of the hot zone and {circumflex over (t)}(s) represents tool axis direction at a distance s along the additive path; and calculating a geometry index hot zone area (A.sub.hot) for the point P(s) on the additive path based upon the calculated geometry index hot zone volume (V.sub.hot) according to the formula A.sub.hot=V.sub.hot/A.sub.max/V.sub.max wherein A.sub.max represents the maximum area of the hot zone and V.sub.max represents the maximum volume of the hot zone.

3. The method of claim 2 wherein calculating EB.sub.melt pool is performed according to the formula EB.sub.melt pool=H.sub.rad+H+H.sub.conv

4. The method of claim 1 wherein the geometry index is: a) based upon a determined geometry case, the geometry case being based upon the geometric description and a point-along-path-data element, the point-along-path-data element being based upon the path description and the selected time; and b) not based upon a calculated hot zone radius.

5. The method of claim 1 wherein the geometry index is: a) based upon a pre-determined value for a path type description for the point P(s) that is based upon the path description; and b) not based upon a calculated hot zone radius.

6. A system for fabricating a part on a substrate using a deposition beam source that emits a beam that follows an additive path composed of a plurality of points, the system comprising: a substrate support for supporting the substrate and manipulating it through space; a metal stock delivery system; a controllable beam source capable of emitting a beam onto the substrate and adapted to form a melt pool thereon; and a system controller including machine code containing the additive-path-with-geometry-metadata of claim 1 and controlling power to the beam source, the system controller programmed to regulate energy of the produced beam for any point on the additive path in accordance with the machine code and the required power predicted by the intra-processing component of the method of claim 1.

7. A system for fabricating a part on a substrate using a deposition beam source that emits a beam that follows an additive path composed of a plurality of points, the system comprising: a substrate support for supporting the substrate and manipulating it through space; a metal stock delivery system; a controllable beam source capable of emitting a beam onto the substrate and adapted to form a melt pool thereon; and a system controller including machine code containing the additive-path-with-geometry-metadata of claim 2 and controlling power to the beam source, the system controller programmed to regulate energy of the produced beam for any point on the additive path in accordance with the machine code and the required power predicted by the intra-processing component of the method of claim 2.

8. A system for fabricating a part on a substrate using a deposition beam source that emits a beam that follows an additive path composed of a plurality of points, the system comprising: a substrate support for supporting the substrate and manipulating it through space; a metal stock delivery system; a controllable beam source capable of emitting a beam onto the substrate and adapted to form a melt pool thereon; and a system controller including machine code containing the additive-path- with-geometry-metadata of claim 4 and controlling power to the beam source, the system controller programmed to regulate energy of the produced beam for any point on the additive path in accordance with the machine code and the required power predicted by the intra-processing component of the method of claim 4.

9. A system for fabricating a part on a substrate using a deposition beam source that emits a beam that follows an additive path composed of a plurality of points, the system comprising: a substrate support for supporting the substrate and manipulating it through space; a metal stock delivery system; a controllable beam source capable of emitting a beam onto the substrate and adapted to form a melt pool thereon; and a system controller including machine code containing the additive-path-with-geometry-metadata of claim 5 and controlling power to the beam source, the system controller programmed to regulate energy of the produced beam for any point on the additive path in accordance with the machine code and the required power predicted by the intra-processing component of the method of claim 5.

10. A method of controlling beam power during an additive deposition process that is controlled by a system controller and that employs a beam source using a power schedule calculation that calculates optimum beam power for any point P(s) along an additive path that will be traveled to form a build, the build having a geometry and being formed from deposited material added to a structure, the method comprising: a pre-processing component and an intra-processing component; during the pre-processing component creating in the system controller a machine code expressed additive-path-with-geometry-metadata. that includes a path description, the path description: a) representing the path of the beam source through space during the additive process; and b) including the geometry index calculated for the point P(s) based upon a geometric description representing the geometry of the build during the additive process, the geometry index comprising predicted thermal characteristics as to a determined melt pool, hot zone and bulk portion for the point P(s); during the intra-processing component calculating at a selected time: a) the temperature of the build at point P(s) during the additive process, the calculation of the temperature of the build including a calculation of hot zone temperature based upon additive-path-with-geometry-metadata; an energy loss of the build at point P(s) during the additive process, the calculation of energy loss being based upon the calculation of temperature and including a calculation of energy conducted from the melt pool to the hot zone; an enthalpy of the build at point P(s) during the additive process; total energy needed at the point P(s) based upon the calculated energy loss and enthalpy of the build at point P(s) on the additive path; an optimum beam source power based upon the calculation of total. energy needed at point P(s) during the additive process; the selected time being the occurrence of a defined event during the additive deposition process or the passage of a selected time interval during the additive deposition process; and modulating the power of the beam emitted by the beam source at point P(s) of the additive path in accordance with the optimum beam source power calculated for point P(s).

11. The method of claim 10 wherein the geometry index is determined by an geometry index calculation, the geometry index calculation including: calculating a geometry index hot zone volume V.sub.hot) for the point P(s) on the additive path according to the formula V.sub.hot=G∩Z (r.sub.hot, P(s), {circumflex over (t)}(s) wherein G represents deposition geometry, Z represents a hemisphere of a certain radius (r) and having a circular surface centered at a point P(s) on the additive path with a normal direction (ň), r h.sub.01 represents the radius of the hot zone and {circumflex over (t)}(s) represents tool axis direction at a distance s along the additive path; and calculating a geometry index hot zone area (A.sub.hot) for the point P(s) on the additive path based upon the calculated geometry index hot zone volume (V.sub.hot) according to the formula A.sub.hot=V.sub.hotA.sub.max/V.sub.max wherein A.sub.max represents the maximum area of the hot zone and V.sub.max represents the maximum volume of the hot zone.

12. The method of claim 10 wherein the geometry index is: a) based upon a determined. geometry case, the geometry case being based upon the geometric description and a point-along-path-data element, the point-along-path-data element being based upon the path description and the selected time; and b) not based upon a calculated hot zone radius.

13. The method of claim 10 wherein the geometry index is: a) based upon a pre-determined value for a path type description for the point P(s) that is based upon the path description; and b) not based upon a calculated hot zone radius.

14. A system for fabricating a part on a substrate using a deposition beam source that emits a beam that follows an additive path composed of a plurality of points, the system comprising: a substrate support for supporting the substrate and manipulating it through space; a metal stock delivery system; a controllable beam source capable of emitting a beam onto the substrate and adapted to form a melt pool thereon; and a system controller including machine code containing the additive-path-with-geometry-metadata of claim 10 and controlling power to the beam source, the system controller programmed to regulate energy of the produced beam for any point on the additive path in accordance with the machine code and the required power predicted by the intra-processing component of the method of claim 10.

15. A system for fabricating a part on a substrate using a deposition beam source that emits a beam that follows an additive path composed of a plurality of points, the system comprising: a substrate support for supporting the substrate and manipulating it through space; a metal stock delivery system; a controllable beam source capable of emitting a beam onto the substrate and adapted to form a melt pool thereon; and a system controller including machine code containing the additive-path-with-geometry-metadata of claim 11 and controlling power to the beam source, the system controller programmed to regulate energy of the produced beam for any point on the additive path in accordance with the machine code and the required power predicted by the intra-processing component of the method of claim 11.

16. A system for fabricating a part on a substrate using a deposition beam source that emits a beam that follows an additive path composed of a plurality of points, the system comprising: a substrate support for supporting the substrate and manipulating it through space; a metal stock delivery system; a controllable beam source capable of emitting a beam onto the substrate and adapted to form a melt pool thereon; and a system controller including machine code containing the additive-path-with-geometry-metadata of claim 12 and controlling power to the beam source, the system controller programmed to regulate energy of the produced beam for any point on the additive path in accordance with the machine code and the required power predicted by the intra-processing component of the method of claim 12.

17. A system for fabricating a part on a substrate using a deposition beam source that emits a beam that follows an additive path composed of a plurality of points, the system comprising: a substrate support for supporting the substrate and manipulating it through space; a metal stock delivery system; a controllable beam source capable of emitting a beam onto the substrate and adapted to form a melt pool thereon; and a system controller including machine code containing the additive-path-with-geometry-metadata of claim 13 and controlling power to the beam source, the system controller programmed to regulate energy of the produced beam for any point on the additive path in accordance with the machine code and the required power predicted by the intra-processing component of the method of claim 13.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0020] FIG. 1 is a diagram showing the hemisphere oriented technique used to estimate hot zone shape and connectivity of the bulk of the workpiece structure.

[0021] FIG. 2A depicts a simplified deposition path utilized when building a wall-shaped structure using additive laser deposition. FIG. 2B depicts juxtaposed graphs showing predicted needed laser power and temperatures of the hot zone, deposit and machine environment over time for an additive metal deposition process used to build the wall structure of FIG. 2A when determined according to the prior art predictive method. FIG. 2C depicts juxtaposed graphs showing predicted needed laser power and temperatures of the hot zone, deposit and machine environment over time for an additive metal deposition process used to build the wall structure of FIG. 2A when determined according to the present invention method. The predicted needed laser power is calculated at various intervals along the additive path process according to a time and event based calculation technique.

[0022] FIG. 3A depicts a simplified deposition path utilized when building a cylindrical-shaped structure using additive laser deposition. FIG. 3B depicts juxtaposed graphs showing predicted needed laser power and temperatures of the hot zone, deposit and machine environment over time for an additive metal deposition process used to build the cylindrical structure of FIG. 3A when determined according to the prior art predictive method. FIG. 3C depicts juxtaposed graphs showing predicted needed laser power and temperatures of the hot zone, deposit and machine environment over time for an additive metal deposition process used to build the cylindrical structure of FIG. 3A when determined according to the present invention method. The predicted needed laser power is calculated at various intervals along the additive path process according to a time and event based calculation technique.

[0023] FIG. 4 is a schematic diagram of a preferred embodiment system for fabricating a part on a substrate using a deposition beam source that follows an additive path and that is controlled in accordance with the two-part laser power prediction method of the present invention.

[0024] FIG. 5 is a flow chart of an embodiment of the overall, two-component method for controlling beam power according to the present invention predictive technique.

[0025] FIGS. 6A-7C depict three alternative methods of determining the additive-path-with-geometry-metadata as part of the first (pre-processing) component of the overall method for controlling beam power according to the present invention predictive technique.

DETAILED DESCRIPTION

[0026] The present invention is directed to a laser power prediction method. The method utilizes the following items to create an improved thermodynamic model of the work piece to compute needed laser power at intervals during the additive process. Like the method of the U.S. Pat. No. 9,573,224, the method described herein does so using: a) an additive path describing the path of the laser through space; b) a geometric representation of the geometry that the additive path is intended to create; c) a description of the thermodynamic characteristics of the manufacturing environment; and d) computed thermophysical characteristics of the materials involved. However, the method does so in a two-part process.

[0027] The first processing step is a pre-processing step that assigns connectivity information to the path for later use. This connectivity information is encoded in the path plan data (namely, the additive-path-with-geometry-metadata). Of note, the determined model of the thermodynamic characteristics of the build is used to predict an appropriate input laser power at definable intervals along the laser path. It accomplishes this by using the path and geometric representation of the part being produced to create a geometry index (a/k/a an idealized geometry) that allows for tenable calculations. However, whereas in the previous method the geometry index was an intermediary value derived at the computer-aided manufacturing (“CAM”) level (i.e., it was continuously being sampled from CAD data), the geometry index is now provided in the machine code from the CAM solution. The geometry index therefore resides as part of the additive manufacturing system instructions (e.g., g-code for a CNC machine) instead of being part of the CAM post-processor action to compute output power. This allows for the system controller to update the output power as needed, rather than expressly follow a predetermined power schedule. The end result is a laser power that is more responsive to the actual activities of the additive manufacturing system. The power update rate is not limitedly tied to the verbosity of the expression in the additive manufacturing system instructions, but rather is a result of the calculations on the laser power prediction control system, which use cues generated from events on the additive manufacturing system to create accurate timing. In particular, when used with a feedback sensor to estimate hot zone enthalpy, corrections to input laser power level can be made much more quickly than with the prior art method. The faster correcting of input laser power level results in a more stable build.

[0028] The second processing step is an on-machine process that computes laser powers as they are required. These computations are performed at event-driven and time driven intervals. In this respect the prediction system receives laser “on” and “off” events from the system executing the additive path. A laser power is computed at a laser “on” event after simulating the cooling process that occurred since the laser “off” event. In terms of time-driven events, laser power is recomputed at a set rate during the period the laser is active.

[0029] A workpiece undergoing additive laser deposition includes various regions that respond to input beam energy. These regions are utilized to create an idealized geometry for a thermodynamic model for the article being built. The workpiece subject to laser impact comprises:

[0030] a. A molten puddle (melt pool);

[0031] b. A region of hot material in direct contact with the molten puddle (hot zone); and

[0032] c. A mass representing the substrate material and part geometry not represented by the other two regions (bulk).

[0033] The laser power input is calculated by performing an energy balance calculation at intervals along the laser's path with the following considerations:

[0034] 1. Conduction from the melt pool to the hot zone;

[0035] 2. Conduction from the hot zone to the bulk of the particular workpiece;

[0036] 3. Conduction from the bulk of the part to the damping system of the CNC machine;

[0037] 4. Radiation from the melt pool;

[0038] 5. Radiation from the hot zone;

[0039] 6. Convection from the melt pool;

[0040] 7. Convection from the hot zone;

[0041] 8. Convection from the bulk structure;

[0042] 9. Process activity at the current point (i.e. is the laser supposed to be on now or not?);

[0043] 10. Mass of the hot zone; and

[0044] 11. Area for conduction between the hot zone and bulk structure.

[0045] Relevant to describing and demonstrating the inventive method are the following elements and symbols, which have the meanings indicated.

[0046] Material properties

[0047] Specific heat, G

[0048] Thermal conductivity, k

[0049] Density, ρ

[0050] Laser absorption coefficient, α

[0051] Liquidus temperature, T.sub.l

[0052] Emissivity, ϵ

[0053] heat of fusion, H.sub.f

[0054] Machine parameters

[0055] Maximum laser power, Q.sub.max

[0056] Ambient temperature, ∞

[0057] Machine temperature, T.sub.mach

[0058] Area for part.fwdarw.machine heat transfer, A.sub.mach

[0059] Powder supply temperature, T.sub.p

[0060] Convection coefficient for the machine environment, h

[0061] Deposition parameters

[0062] Layer thickness, d.sub.l

[0063] track width, w

[0064] Remelt thickness as a fraction of layer thickness, r

[0065] Scan speed, f

[0066] Substrate mass, m.sub.sub

[0067] Target superheat, T.sub.super

[0068] Laser Predictive Power Settings

[0069] Target hot zone radius, r.sub.hot

[0070] Target melt pool temperature, T.sub.melt=T.sub.l+T.sub.super

[0071] Target hot zone temperature, T.sub.target

[0072] Laser power prediction interval, i.sub.p

[0073] Nomenclature

[0074] Let s describe the length along the deposition path.

[0075] Let P describe the deposition path such that P(s) is the Cartesian coordinate of a point s along the path.

[0076] Let {circumflex over (t)}(s) describe the tool axis direction at a distance of s along the path.

[0077] Let G describe the deposition geometry.

[0078] Let Z (r, p, {circumflex over (n)}) describe a. hemisphere of radius r with the circular surface centered at point p with a normal direction of {circumflex over (n)}.

[0079] An exemplary calculation method thus proceeds as follows. First, before beginning the formal calculations representing the in-process thermodynamic environment, initial values are input into the data structures.

[0080] The initial energy content of the build is expressed in Equation 1:

H.sub.build=C.sub.pT.sub.ambm.sub.sub  (1)

[0081] The initial build mass is assumed to be the equivalent to the substrate, as expressed in Equation 2:

m.sub.build=m.sub.sub  (2)

[0082] The initial temperature of the build is expressed as follows:

T.sub.build=T.sub.∞  (3)

[0083] The initial temperature of the hot zone temperature is expressed in Equation 4:

T.sub.hot=T.sub.build  (4)

[0084] The distance along the path is initialized at 0 in Equation 5:

s=0  (5)

[0085] It is useful to pre-compute items that will remain constant throughout the computation,

[0086] Maximum volume of the hot zone is calculated as such:

V.sub.max=2/3r.sup.3.sub.hot.sup.π  (6)

[0087] The maximum hot zone area is calculated in Equation 7 as:

A.sub.max=2/r.sup.2.sub.hot.sup.π(7)

[0088] The surface area of the melt pool can be expressed as:

A.sub.melt1/4w.sup.2π  (8)

[0089] While mass per unit length of the additive path is expressed in Equation 9 as:

m.sub.ul=wdip  (9)

[0090] The time interval needed to travel one inspection interval along the additive path is computed via the following equation:

Δt=i.sub.p/f  (10)

[0091] With the above initial pre-calculation steps undertaken, one can proceed to calculate needed laser power along a given deposition path according to the following method.

[0092] A. Calculate the geometry index for the current position.

V.sub.hot=G∩Z(r.sub.hot, P(s), {circumflex over (t)}(s))  (11)

A.sub.hot=V.sub.hot/A.sub.max/V.sub.max  (12)

[0093] Equation 12 can be simplified using the geometry index to:

A.sub.hot=A.sub.max*g.sub.idx

[0094] B. Compute the rate of heat transfer from the melt pool. Note that equations 13-15 refer to a time period, Δt. In the on-machine predictive system discussed herein, Δt is defined as the elapsed time since the last event triggering the prediction system.

[0095] 1. Energy transferred via conduction from the melt pool:

H.sub.cond=kA.sub.hot/r.sub.hot(T.sub.melt-T.sub.hot)Δt  (13)

[0096] 2. Energy lost due to convection from the melt pool:

H.sub.conv=hA.sub.melt(T.sub.melt-T.sub.hot)Δt  (14)

[0097] 3. Energy released by the melt pool via radiation:

H.sub.rod=ϵσA.sub.melt(T.sup.4.sub.melt-T.sup.4.sub.hot)Δt  (15)

[0098] C. Compute an amount of energy needed to bring incoming material and remelted material up to the desired melt pool temperature. Note that the following equations refer to a prediction interval i.sub.p. In the disclosed invention i.sub.p is defined as the distance traveled since the last event triggering the prediction system.

[0099] 1. Energy required to melt the powder:

H.sub.powder=m.sub.uli.sub.p(T.sub.melt-T.sub.p)C.sub.p+H.sub.f)  (16)

[0100] 2. Energy required to remelt the existing material:

H.sub.conv=hA.sub.melt(T.sub.melt-T.sub.hot)Δt  (17)

[0101] 3 By summing the energy losses at the melt pool and the energy needed to melt material one can calculate a total input energy at this step used to estimate a viable laser power as follows:

H.sub.total=H.sub.powder+H.sub.remeltH.sub.rad+H.sub.cond+.sub.conv  (18)

[0102] 4. Calculate a viable laser power:

Q.sub.laser=min(Q.sub.max, H.sub.total/αΔt)  (19)

[0103] The expression in Equation 19 is constrained by the output capacity of the laser system employed.

[0104] D. The mass of the build is now updated to include the newly deposited material.

m.sub.build=m.sub.build+m.sub.uli.sub.p  (20)

[0105] E. Calculate the energy losses at the build.

[0106] 1If accurate model data is not readily available at computation time, an effective radius and area of the build for the purposes of heat transfer calculations can be estimated via Equations 21 and 22.

[00001]$\begin{array}{cc}{r}_{\mathrm{build}}=3\sqrt[3]{3m_{\mathrm{build}}\text{/}{\rho }^{4\pi }}& \left(21\right)\\ A_{\mathrm{build}}=4r_{\mathrm{build}}^2\pi & \left(22\right)\end{array}$

[0107] Next significant losses to the environment from the hot zone, the build, and the machine are estimated in Equations 23, 24 and 25.

[0108] 2 Energy radiated from the hot zone:

H.sub.hot=∈σA.sub.hot(T.sup.4.sub.hot-T.sup.4.sub.∞)Δt  (23)

[0109] 3. Energy conducted from the build to the machine:

H.sub.mach=(kA.sup.machine/r.sub.build)(T.sub.build-T.sub.mach)Δt  (24)

[0110] 4. Energy lost due to convection to the machine environment:

H.sub.env=hA.sub.build(T.sub.build−T.sub.∞)Δt  (25)

[0111] F. Update enthalpy (internal energy) of the build.

H.sub.build=H.sub.build+Q.sub.laserΔt−H.sub.mach−H.sub.hot−H.sub.env  (26)

[0112] G. Update (recompute) the build temperature.

T.sub.build/=H.sub.build/C.sub.p  (27)

[0113] H. The hot zone temperature estimate is then recomputed in Equation 28. The tuning parameter is used to scale how quickly the hot zone temperature estimate scales with increasing build temperature. In practice, a value of 0.5 was sufficient to provide reasonable results.

T.sub.hot=ƒT.sub.build+(1-ƒ)T.sub.target  (28)

[0114] I. In Equation 29 the position along the deposition path is updated. Note, however, that in using the on-machine predictive system of the present invention, the position along the path s, is determined by the motion system, The predictive system is event-driven rather than position-driven. The additive manufacturing system sends a signal to the predictive software for when the additive process begins and ends. It is the timing of these signals that drives the computation of Δt and i.sub.p.

s=s+i.sub.p  (29)

[0115] J. If the path is not complete, return to step 1.

[0116] FIG. 1 shows the “hemisphere” hot zone size estimation of various points on a build. Here, the hot zone estimation is referred to herein as a “geometry index.” In the method of the prior art patent, the geometry index was used wholly within the CAM system as part of the preparatory step to embed output power into the additive manufacturing system instructions. In contrast, in the method disclosed herein, the pre-computation i.e,, preparatory) steps end with the geometry index. The geometry index is a normalized, composite term, defined at any point along the additive path as the ratio of the ratio of the volume of the intersection of the ideal hot zone geometry with the part geometry, divided by the ideal hot zone volume. Thus, a value of 1 would indicate the case of a semi-infinite plate, while a smaller value such as 0.1 would suggest a less connected condition, such as a thin wall or rib. From a human standpoint, it is a more natural thing to consider when working with the prediction system. Essentially, it says how well the particular hot zone region under consideration is connected to the bulk structure in terms proportion. So a value of 0 means it is not connected at all (i.e. some fictitious hot zone alone in its own universe), a 1 means it is fully connected (i.e. the hot zone lies on an infinite flat plate). This does not change the model mathematically, but it simplifies its execution and makes it more intuitive. The geometry index is computed as:

[0117] g.sub.idx=GZ(r,p,{circumflex over (n)})/V.sub.max

[0118] FIG. 2A depicts a simplified deposition path utilized when building a wall-shaped structure using additive laser deposition. FIG. 2B depicts juxtaposed graphs showing predicted needed laser power and temperatures of the hot zone, deposit and machine environment over time for an additive metal deposition process used to build the wall structure of FIG. 2A when determined according to the prior art predictive method. The periodicity of the predicted laser power shown in FIG. 2B is due to the laser reaching the corners of the wall structure, which reduces the area available for heat transfer away from the melt pool, resulting in less required laser power to maintain the desired energy balance. The predicted needed laser power is calculated at various intervals along the additive path process according to the prior art calculation technique. The computer controlling the laser power can be programmed to adjust laser power along the deposition path in accordance with the predicted levels of needed energy.

[0119] Of note, the present invention laser power prediction method can be compared to the prior art method. FIG. 2C depicts juxtaposed graphs showing improved predicted needed laser power and temperatures of the hot zone, deposit and machine environment over time for an additive metal deposition process used to build the wall structure of FIG. 2A when determined according to the present invention method. The predicted needed laser power is calculated at various intervals along the additive path process according to a time and event based calculation technique. By comparing the graph of FIG. 2B to the graph of FIG. 2C one can see that for the simplified structure of FIG. 2A, the initial predicted power levels are similar. Again, note that the predicted laser power requirement has a periodic component superpositioned with an exponential decay component. The periodic component is due to the hack-and-forth motion caused by the shape of the desired geometry. The exponential decay component is a result of the part and workpiece heading towards some steady state temperature. Note, however, as shown in FIG. 2C, when event time-based interruptions are taken into account, the predicted power levels and temperatures of the system components drop at certain points.

[0120] FIG. 3A depicts a simplified deposition path utilized when building a cylindrical-shaped structure using additive laser deposition. The helical path for creating a cylinder as seen in FIG. 3A is as close to a steady state process as metal additive manufacturing will ever get. This is one of the reasons that a hollow cylindrical boss is a common demonstration structure for new metal additive manufacturing systems. FIG. 3B depicts juxtaposed graphs showing predicted needed laser power and temperatures of the hot zone, deposit and machine environment over time for an additive metal deposition process used to build the cylinder structure of FIG. 3A when determined according to the prior art predictive method. The predicted power shown in top graph of FIG. 3B is very stable compared to the case described in FIG. 2B. This is because the helical path shown in FIG. 3A is much closer to a steady state than the raster path used in the rectangular solid of FIG. 2A.

[0121] FIG. 3C depicts juxtaposed graphs showing predicted needed laser power and temperatures of the hot zone, deposit and machine environment over time for an additive metal deposition process used to build the cylinder structure of FIG. 3A when determined according to the present invention method. The predicted needed laser power is calculated at various intervals along the additive path process according to a time and event based calculation technique. By comparing the graph of FIG. 3C to the graph of FIG. 3B one can see that for the cylindrical structure of FIG. 3A, the predicted power levels are similar at the beginning of the process, but deviate significantly when event interruptions are factored in. With the current prediction method the needed power input is computed to stabilize discretely and not gradually.

[0122] Using the calculations above, the laser energy delivered during an additive path deposition process can be controlled on a time basis for each point along the path by regulating power to the laser at each point in accordance with the required power predicted by the above calculation methods. An embodiment of the inventive method can be programmed into a computer-controlled laser metal deposition system to improve the fabrication of components of almost any geometry that can be produced from a computer database. This system for fabricating a part on a substrate using a deposition beam source that follows an additive path is depicted in FIG. 4.

[0123] In the depicted system, the pre-processing steps leading to the calculation of the additive-path-with-geometry-metadata 40 (discussed in more detail infra) are shown in the left side panel of FIG. 4. Additive-path-with-geometry-metadata 40 is input into system controller 50 for use in controlling the additive manufacturing process during operation. The system includes a substrate support (CNC system or table) 55 for supporting the workpiece (substrate) 54 and manipulating it through space, a material feeder 51 (a metal stock delivery system such as a powdered metal injector/sprayer or wire feeder) and a laser 52, both of which are typically found in a processing head 53. Laser 52, system controller 50 and processing head 53 work together to emit a beam (dotted line impacting workpiece) onto workpiece 54 and form a melt pool thereon and control power to the beam source, System controller 50 is programmed to regulate energy of the produced beam for any point on the additive path in accordance with the required power predicted by power prediction service component 56.

[0124] System controller 50 regulates power to laser 52 for each point on the additive path as a function of preprocessing data 40 and the other input data shown in FIGS. 4 and 5. As shown in the right panel of FIG. 4, the system calculates H.sub.total 60 using enthalpy (H.sub.build) and temperature (T.sub.build) data 58, machine environment data (data regarding ambient temperature, initial temperature of the build and other environmental conditions) 57 plus input from the power prediction service 56. That determined H.sub.total 60 is used to compute laser power (Q.sub.laser) 46, which computation is fed back into the power prediction service 56 to input into the system controller 50. With that input along with the additive path and geometry metadata 40, the system controller sends control signals to the material feeder 51, laser 52 and CNC system 55. These inputs discussed hereafter are used to produce overall additive path and power data inputs for the system controller. As the deposition process is additive and the formula recursive, computed laser power 46 is also a data element that factors into enthalpy and temperature data 58.

[0125] As noted, the method described herein has certain components in common with the method of U.S. Pat. No. 9,573,224. In the prior art method the important geometry factor is computed according to that patent's equations nos. 11-12 and depicted in the method schematically shown in that patent's FIGS. 7A-7C. The geometry factor computing process was part of the in-process deposition process. In contrast, with the instant invention, there are now two separate, but interrelated processes and the geometry index determination component of the method, though still important to predicting laser power, is separated from the in process computations.

[0126] In one embodiment the invention is directed to a method of controlling beam power during an additive deposition process. FIG. 5 is a flow chart showing the intra-process component of beam power controlling method according to the present invention predictive technique, As shown in FIG. 5, system controller 50 uses event timing input 70 to determine a computation event trigger 71 which leads to a determination of time since the last power computation 72. That determination 72 along with the determined geometric index 20 for the point on the build in system controller 50 is used to compute a predicted laser power 46 for the point P(s) on the build. Computed laser power 46 and the additive-path-with-geometry-metadata 40 are input into system controller 50. System controller 50 then sends a power modulation signal to material feeder 51 and laser 52 of processing head 53. As can be appreciated, timing inputs of 70, 71 and 72 provide an indication to system controller as to where the additive process is timing-wise and thus a more accurate prediction of needed laser power. In other word, the additive-path-with-geometry-metadata plus the event selection inputs in system controller 50 is used by as a basis to determine when, as opposed to simply where, to do a beam power modulation calculation.

[0127] As shown in the figures, the deposition process is controlled by a system controller 50, The deposition process uses a beam source (laser) 52 to project energy to a point P(s) along an additive path that will be traveled to form a build. The build has a geometry and is formed from deposited material added to a structure. The point P(s) along the additive path has associated with it a geometry index (idealized geometry) that comprises a melt pool, a hot zone and a bulk portion. The method includes a pre-processing component and an intra-processing component.

[0128] The pre-processing component involves creating in the system controller a machine code expressed additive-path-with-geometry-metadata that includes a path description. The path description represents the path of the beam source through space during the additive process and includes the idealized geometry calculated for the point P(s) based upon a geometric description representing the geometry of the build during the additive process.

[0129] The intra-processing component comprises several calculations. The first calculation calculates an energy balance at the melt pool (EB.sub.melt pool) for the point P(s) on the additive path. The energy balance calculation is based upon the following calculations: a) a calculation of energy radiated from the melt pool (H.sub.rod); b) a calculation of energy conducted from the melt pool to the hot zone (H.sub.cond), the calculation of H.sub.cond being based upon the additive-path-with-geometry metadata; and c) a calculation of energy lost due to convection at the melt pool (H.sub.conv).

[0130] The second calculation calculates total energy (H.sub.total) needed at the point P(s) on the additive path according to the following equation

[0131] H.sub.total=H.sub.deposited material+H.sub.remelt+EB.sub.melt pool

[0132] wherein H.sub.deposited material represents energy required to melt the deposited material and H.sub.remelt represents energy required to remelt existing material.

[0133] The third calculation calculates optimum beam source power Q.sub.source for the point P(s) according to the following equation

[0134] Q.sub.source=min(Q.sub.max, H.sub.total/αΔt)

[0135] wherein Q.sub.max represent maximum laser power, α represents a beam absorption coefficient and Δt represents a calculation interval (the time between recomputing laser power)

[0136] The energy projected by the beam source is modulated based upon the calculated optimum beam source power Q.sub.source.

[0137] The additive-path-with-geometry-metadata can be created in alternative ways. These alternative methods are depicted in FIGS. 6A through 6C. The first method (referred to as the “detailed method”) is shown in FIG. 6A. This figure describes the creation of the geometry index, which is basically the volume fraction of the hemispheres seen in FIG. 1 that intersect the part. The creation of a geometry index is part of the original laser power prediction method. However, in the present invention, what was an intermediary value in the original predictive method is now the final number of interest to be included with the path plan data in lieu of the laser power. According to FIG. 6A, path planning data 30 and inspection interval data 11 (distance along the path between recomputing the geometry index) are used as inputs to determine a point along the path data 12, Hot zone radius data 13 and point along the path data 12 are used to calculate hot zone hemisphere data 14. That hot zone hemisphere data 14 along with the input geometry 10 of the build are used to determine an intersect hot zone hemisphere with input geometry value 15. The determined intersect hot zone hemisphere with input geometry value 15 yields geometric description 25. In the case of FIG. 6A the geometric description is labeled as 25A and is a determined intersection to hot zone hemisphere ratio. That ratio yields a geometry index (a/k/a geometry factor) g.sub.idx 20. The geometry index 20 and the path planning data 30 are combined to create “additive-path-with-geometry-metadata” 40, which will be expressed in machine code for configuration in the system controller. Due to the reliance on boolean operations to compute the intersection, the metadata creation can be computationally expensive. Thus, applicant has devised the two alternate levels of approximation in FIGS. 6B and 6C.

[0138] In a second method (referred to as the “geometry case method”) shown in FIG. 6B path planning data 30 and inspection interval data 11 are used as inputs to determine a point along path data 12. That point along path data 12 along with the input geometry 10 of the build are used to determine a geometric description (in this case a geometry case referable to the particular structure; edge, corner, infill, overhang) 25B for the subject point along the path, which yields a geometry index (a/k/a geometry factor) g.sub.idx 20. The geometry index 20 and the path planning data 30 are combined to create “additive-path-with-geometry-metadata” 40, which will be expressed in machine code for configuration in the system controller. With the method of FIG. 6B, the geometry index value can be chosen categorically based on a case-based approach. In this respect, a corner may have set value and a point on a long edge would have another value. This is tenable via a signed distance function for the model, which is more expensive than a single boolean operation, but only needs to he done once. After computation, one can sample the field generated via the signed distance function (i.e. a big voxel domain in which each voxel contains the nearest distance to the surface of the part) to determine the case.

[0139] In a third method (referred to as the “path case method”) shown in FIG. 6C, path planning data 30 is used to determine a geometric description in the form of a. path type 25C on the build (e.g., perimeter, infill) for the subject point under processing. That path type 25C information yields a geometry index (a/k/a geometry factor) g.sub.idx 20 for the point. The geometry index 20 and the path planning data 30 are combined to create “additive-path-with-geometry-metadata” 40, which will be expressed in machine code for configuration in the system controller, The method of FIG. 6C is the simplest geometry index estimation. This method involves skipping analyzing the model in its entirety. A value is assigned based on the type of path in the path plan. For example, it is already known whether the path is on the perimeter or in the interior based on metadata already within the path plan because this information is needed for path sorting algorithms in other areas of my CAM tools. Assigning geometry index values based on this method is essentially “free” in a computational sense. In testing, the method proved to be acceptable, but less reactive to changes in geometry. On the other hand, the method. of FIG. 6A will be the most reactive to what is actually being built, However, if the part being built is small, faster estimates are more than acceptable. For example, the method of FIG. 6C can be used to create a multitude of small parts on a plate for some characterization of a material.

[0140] The invention is further directed to a method of controlling beam power during an additive deposition process that is controlled by a system controller and that employs a beam source. The method uses a power schedule calculation that calculates optimum beam power for any point P(s) along an additive path that will be traveled to form a build. The build has a geometry and is formed from deposited material added to a structure. The point P(s) along the additive path has associated with it an idealized geometry comprising a melt pool, hot zone and bulk portion. The method comprises a pre-processing component and an intra-processing component.

[0141] The pre-processing component of the method involves creating in the system controller a machine code expressed additive-path-with-geometry-metadata that includes a path description. The path description: a) represents the path of the beam source through space during the additive process; and b) includes the idealized geometry calculated for the point P(s) based upon a geometric description representing the geometry of the build during the additive process.

[0142] The intra-processing component involves:

[0143] calculating the temperature of the build at the point P(s) during the additive process, the calculation of the temperature of the build includes a calculation of hot zone temperature based upon additive-path-with-geometry-metadata;

[0144] calculating an energy loss of the build at point P(s) during the additive process, the calculation of energy loss being based upon the calculation of temperature and including a calculation of energy conducted from the melt pool to the hot zone;

[0145] calculating an enthalpy of the build at point P(s) in time during the additive process; calculating total energy needed at the point P(s) based upon the calculated energy loss and enthalpy of the build at point P(s) on the additive path;

[0146] calculating an optimum beam source power based upon the calculation of total energy needed at point P(s) during the additive process; and

[0147] modulating the power of the beam emitted by the beam source at point P(s) of the additive path in accordance with the optimum beam source power calculated for point P(s).

[0148] The invention is also directed to a system for fabricating a part on a substrate using a deposition beam source that emits a beam that follows an additive path composed of a plurality of points. The system comprises a substrate support for supporting the substrate and manipulating it through space. The system further includes a metal stock delivery system and a controllable beam source capable of emitting a beam onto the substrate and adapted to form a melt pool thereon. The system also includes a system controller that has machine code containing the additive-path-with-geometry-metadata described above and that controls power to the beam source. The system controller is programmed to regulate energy of the produced beam for any point on the additive path in accordance with the machine code and the required power predicted by the intra-processing component via any of the methods described above.

[0149] The live computation from a runtime-derived timing allows for the double checking of results via feedback. In this respect sensors such as spectrometers, welding cameras (for welding applications) or a 3D scanning system could be integrated into the system and method, A spectrometer, for instance, can be used to detect the emission spectra from the melt pool and hot zone regions. A neural network could be added to predict model errors or process problems using a spectrometer and camera.

[0150] The process described herein can also be used in connection with other traditional welding techniques, such as tungsten inert gas (“TIG”) welding, gas metal arc welding (“GMAW”), plasma transferred arc (“PTA”) welding and electron beam (“EB”) welding. While the embodiments of the method and system of the present invention have been described herein, numerous modifications, alterations and changes to the described embodiments are possible without departing from the scope of the invention. The embodiments described herein are not intended to be limiting.