The invention relates generally to the measurement and control of laser cladding process. In particular, the invention relates to the temperature profile control of direct metal deposition.
Direct Metal Deposition (DMD) is a material additive manufacturing technology utilizing a precisely controlled laser beam to melt powders onto a substrate to form products. DMD with a closed loop control system has been successfully applied in complicated part prototyping, repairs and surface modifications [1].
DMD is a multi-parameter process where laser power, traverse speed and powder feed rate are considered the most dominant parameters that determine the dimensional accuracy and mechanical properties of products. Other secondary important parameters include laser beam size, delivery and shielding gases, nozzle design, bead overlap, z increment, tool path design, and powder qualities. Any disturbance from the controlling parameters, environment, and pool itself (surface tension, flow-ability), may shift the process away from its stable point and result in defects in the produced parts.
Existing sensing and modeling efforts have been focused on cladding tracks and molten pools. Monitoring cladding tracks can directly provide dimensional information regarding depositions [8]. However, monitoring cladding tracks introduces inherent process delays which must be compensated for in the controller. On the other hand, sensing molten pools can provide online process information, which could enable real time process control without process delays [1].
Optical intensity [1] and infrared images [10] of molten pools have been successfully employed to control the cladding process. Pool temperature measurement and transient mathematical modeling of the process have been reported by Han et al [6, 7]. Processing infrared images leads to complex calculations, and is therefore slower than either optical intensity or temperature measurements.
A dynamic model of the process is essential for advanced model based closed-loop controller designs. Several theoretical and numerical models have been studied to give the insight of the process [3-7]. However, because of limitations, complexities and extensive numerical operations of the simulations, these models are not practical for in-process control. Experimental-based modeling using system identification has been reported to identify the nonlinear input-output dynamic relationship between traverse velocity and deposition bead height [8]. However, significant deviations existed between the actual data and the model outputs.
To overcome the difficulties of the system modeling, a fuzzy logic controller was implemented where only the fuzzy knowledge of the process was needed [9]. Mazumder et al proposed a closed-loop controlled DMD system, in which three photo-detectors were used to monitor the molten pool height [1, 2]. A control unit, where an OR logic function was operated on the three signals from photo-detectors, was used to trigger off the laser when the detected pool height was above the pre-set limits. This closed-loop control system proved to be successful in controlling the dimensional accuracy of the produced parts. POM Group Inc. in Auburn Hills has commercialized the system and installed the system on three different continents.
While various methods have been developed to monitor and control the laser cladding process, such methods can, nevertheless, be the subject of certain improvements. In this regard, conventional measurement and controlling methods for laser cladding are not sufficiently efficient and robust for large scale production. Thus, it would be advantageous to provide robust, reliable and efficient methods for direct metal deposition for commercial production.
This invention improves upon existing process-control methodologies by providing a model predictive system that controls the molten pool temperature during DMD process. The preferred embodiment includes a two-color pyrometer used to measure the molten pool temperature, and a real-time controller implementing a generalized predictive control algorithm with constraints.
The dynamics describing the relationship between the laser power and the molten pool temperature are used to design the generalized predictive controller. A Kalman filter is used to estimate the states. A reference temperature profile including a sine wave and three step changes demonstrated that the predictive controller successfully stabilizes the DMD process. More particularly, the approach improves the temperature profile during the deposition process to improve end-product microstructure and/or dimensional accuracy.
According to the invention, a method of controlling a DMD process comprising the steps of identifying temperature dynamics associated with the melt pool, and generating excitation signals to control the laser as a function of the temperature dynamics using a generalized predictive control algorithm with input constraints.
The step of identifying temperature dynamics associated with the melt pool is carried out with a two-color pyrometer that senses in regions of the spectrum different from that used by the laser used to form the melt pool, which may be a diode laser, a fiber laser, or a CO2 laser. In the preferred embodiment, the two-color pyrometer senses in bands at 1.3 μm and 1.64 μm.
The excitation signals may comprise random amplitudes or random durations in a predetermined range. The step of identifying temperature dynamics associated with the melt pool may comprise model order selections, step response comparisons and residual analysis among different models structures.
The generalized predictive control algorithm may use space-state models, including space-state models that can be scaled into multiple-input and multiple-output systems to implement other control parameters such as the pool geometry and plume plasma radiation so as to control product dimensions or compositions. The generalized predictive control algorithm may further use a dual active-set method with modifications.
A two-color pyrometer 102 is connected by fiber 104 to a collecting lens to monitor the molten pool temperatures. Two-color detection was chosen for its accurate temperature measurement. A dSPACE 1104 controller was used as the real time controller to implement the generalized predictive control algorithm. The measured molten pool temperature was relayed to the controller. The function of the controller is to compare the molten pool temperature to the reference values and calculate the optimal output of the laser power.
The selection of the model structures and the excitations is critical to obtain an accurate dynamic model. The characterization of the input-output pair determines the maximum accuracy that can be achieved by a model independent structure. For a linear system, a pseudo-random binary signal train is normally used to excite the system. The system dynamic model can be obtained by a least-square algorithm. For a nonlinear system, the excitation signals need to cover the entire plant's operating range because the nonlinear models seldom extrapolate accurately. A rich spectrum of excitation amplitudes and frequencies is thus desirable.
The amplitudes of the excitations should be changed around a desired working point. The range of the amplitude reflects the operating range where the model parameters are valid. The frequency components of the excitations determine if a frequency response is correct. Low frequency signals have long pulse durations, which give the correct steady state response. High frequency signals, on the other hand, have shorter pulse durations, which give the transient response [11]. Therefore, the best excitation signal is a series of pulses of random amplitudes and widths.
a shows a voltage train that was applied to the control port of the laser. The voltage values and the voltage periods were randomly changed. The voltage amplitudes were random variables with Gaussian probability density functions. The mean value is 1.7 volts and the standard deviation is 0.2 volts. The pulse durations randomly changed between 10 milliseconds to 5 seconds. The randomly generated voltage levels passed through a saturation gate with a lower saturation value of 1.5 volts and an upper saturation value of 1.9 volts
H13 tool steel powder was deposited on the low carbon steel to form a single track. The powder flow rate was 12 grams per second. The shielding gas was Argon (25 psi) and delivery gas was Helium (20 psi). The traverse speed was 14.4 inches per minute. The beam size on the substrate was 1.0 mm.
In
where z−1 is the single sampling interval delay operator. The filter has a static gain of 1. The low pass filter should be able to filter out the high frequency noise, but still capture the transient response of the dynamics. In order to capture a 300 ms transient response, a 3 dB bandwidth of the filter should be greater than 3.3 Hz. Therefore, β was chosen to be 0.8, which corresponds to a 3 dB bandwidth of 3.5 Hz. The frequency response of the filter is shown in
In order to get the system dynamic model, two portion signals were used, as shown in
The model was identified using four different model structures, state space model, Box Jenkins model, output error model and auto-regressive with moving average with external inputs (ARMAX) model. Comparing the four step responses in
State-space model has the form
X(k+1)=AX(k)+Bu(k)+Ke(k) (2.1)
y(k)=CX(k)+Du(k)+e(k) (2.2)
where X is the state vector, y denotes the process output to be controlled and u denotes the process input (controller output). A, B, and C are the matrices defining the state-space model.
The identified model has matrix values of:
From
The identified model output was compared to the measured data, as shown in
Predictive control is a multi-step approach, combining feed forward and feedback control design [3]. Feed forward is represented by predictions based on a mathematical model and is the dominant component of control actions. Feedback from measured output serves as compensation for some bounded model inaccuracies and low external disturbance. The design consists of local minimization of quadratic criterion, in which the predictions of future outputs are involved. The predictions are determined from the model describing the system dynamics. At each time step, predictions and minimization of the quadratic criterion are repeated to give the next optimal control.
Generalized Predictive Control Algorithm with Input Constraints
From equation (2.1-2.2), the N step prediction {circumflex over (X)}(k+N) and ŷ(k+N) can be expressed as:
{circumflex over (X)}(k+N)=ANX(k)+AN−1Bu(k)+ . . . +Bu(k+N−1) (3.1)
{circumflex over (y)}(k+N)=CANX(k)+CAN−1Bu(k)+ . . . +CBu(k+N−1) (3.2)
The cost function to minimize is:
The cost function is expressed in step k, over indicated horizons. N is the optimization horizon, N0 is the initial insensitive horizon, and Nu is the control horizon.
One of the major advantages of generalized predictive control is its ability to take systematic account of constraints, as they can easily be incorporated into the optimization (Equation 4). The DMD system considered here only has an input constraint that constrains the laser power since the model is valid only within a certain laser input power range. Assuming that the input of the plant after prediction horizon Nu is the same as at step Nu (u(k+N
The minimization of equation 4 with constraints Equation 5 is known as a quadratic programming (QP) problem. The algorithm solving this problem is based on Goldfarb and Idnani's dual active-set method [12] with modifications from [13].
The T-filter is a low pass filter that improves prediction accuracy in the high frequency range by reducing the transference of high frequency noise. A T-filter can also improve the high frequency range sensitivity by reducing the input sensitivity to high frequency noise. Equation 1 is the form of the T filter that was used in the control system design.
The noises in the state space model (Equation 2.1 and 2.2) are assumed to be white, mutually independent and normally distributed □ (mean, covariance) with zero mean and known positive definite covariance. The state estimate of {circumflex over (X)}(k+1,k) can be expressed as
{circumflex over (X)}(k+1,k)=A{circumflex over (X)}(k,k−1)+AK(k)(y(k)−C{circumflex over (X)}(k,k−1))+Bu(k) (6)
where K(k) is the Kalman filter gain.
The generalized predictive control algorithm with input constraints was first simulated in Matlab-Simulink environment using the model identified in the previous section. Then the control algorithm was implemented in dSPACE real time controller. In view of the strong noise from the molten pool temperature measurement, an extra 20-point moving average filter was used to filter out the noise. A temperature profile including a sine wave and three step changes was used as the tracking reference. In order to test the large range input controllability using the identified model, the upper limit and lower limit of the voltage applied to the laser was softened to +0.4V and −0.5V, respectively. The reference temperature ranged from −200° C. to +200° C. The control results are shown in
The results showed that the controller can successfully track the reference temperature by adjusting the voltage supplied to the laser power controller. Compared to the on-off controller, a generalized predictive controller can provide smooth tracking of the references. It would be difficult for an on-off controller to get the desired temperature profile.
This application claims priority from U.S. Provisional Patent Application Ser. No. 60/866,150, filed Nov. 16, 2006, the entire content of which is incorporated herein by reference.
Number | Date | Country | |
---|---|---|---|
60866150 | Nov 2006 | US |