Digital Implementation of LCC Resonant Converters for X-ray Generator with Optimal Trajectory Startup Control

: High voltage LCC resonant converters have been widely used in X-ray imaging systems in automobile nondestructive testing (NDT) applications. Low ripple voltage waveforms with fast-rising time under no-overshoot response are required for safety in such applications. The optimal state trajectory control (OTC) based on the state plane model is one of the most effective control methods to optimize transient response. Dynamic variations of the resonant voltages/currents are described as corresponding trajectories on the state plane. The transient relations can be determined by evaluating the geometric relationships of the trajectories. However, the LCC resonant converter has more state variables, resulting in more complex calculations that make the state trajectory control challenging. Furthermore, the startup duration is the most demanding process of the state trajectory control. In this paper, a digital implementation based on a hybrid controller built in a ﬁeld-programmable gate array (FPGA) is proposed for LCC resonant converters with optimal trajectory startup control. A coordinated linear compensator is employed to control the switching frequency during steady-state conditions, hence eliminating the steady-state error. The experimental results were conducted on a 140-kV/42-kW LCC resonant converter for an X-ray generator. It achieves a short rising time of output voltage with no additional current or voltage stress in the resonant tank during startup compared to the conventional digital implementation control. This article proposes a state trajectory startup control and real-time implementation for ISOS Type LCC resonant converters with a capacitive output ﬁlter in a 140 kV/42 X-ray generator. The state trajectory startup control for LCC converter is optimized and conducted using FPGA to achieve low ripple voltage waveforms with fast-rising time under no-overshoot response. Furthermore, the state trajectory-based startup is successfully applied to an industry-level X-ray generator with a Power/ Voltage Level of 42 kw/140 kV for the ﬁrst time. The rest of the paper is organized as follows.


Introduction
Effective detection and identification must be carried out to ensure the product quality of the automobile hub. Nondestructive testing technology (NDT) plays an essential role in improving product quality and production efficiency. At the same time, it is one of the most effective testing methods in testing technology. X-ray imaging technology is widely used in nondestructive testing because of its intuitive and accurate test results [1]. Due to the high quality of X-ray imaging, it can accurately provide real information about automobile hub defects. Therefore, the digital image of the wheel hub can be obtained through the X-ray imaging system. The X-ray generator is a device that supplies electric power to the X-ray tube and permits the selection of X-ray energy (voltage), X-ray quantity (current), and exposure time. The generator has three main interrelated electric circuits to serve three main functions: (1) The power supply for heating the cathode filament and evaporating electrons is provided by the filament circuit; (2) the high-voltage circuits are used to increase the speed of electrons from the cathode to the anode to produce X-rays; (3) the timer circuit (exposure timer) controls the duration of X-ray. Figure 1 shows a block diagram of a typical X-ray generator [2]. The X-ray tube application presented in this paper utilizes an Input Single Output Series (ISOS) LCC resonant converter [3], as illustrated in Figure 2, which features a multilevel series arrangement with diode rectified output capacitors [4], increasing flexibility and scalability. converter (PRC), and LLC resonant converters [10], state trajectory analysis and control have been applied in addressing unpredictable dynamics, burst mode for light load [11], soft startup, and short circuit protection [12]. However, more state variables make it more complicated for LCC resonant converters to present their relationships in a 2-D plane coordinate system. Therefore, it is very challenging to implement the state trajectory control. Furthermore, the startup duration is the toughest process of the state trajectory control [13].   soft startup, and short circuit protection [12]. However, more state variables make it more complicated for LCC resonant converters to present their relationships in a 2-D plane coordinate system. Therefore, it is very challenging to implement the state trajectory control. Furthermore, the startup duration is the toughest process of the state trajectory control [13].   Figure 2. The topology of LCC resonant converter. Researchers have proposed various modeling methods to operate resonant converters properly. Fundamental harmonic approximation (FHA) is widely utilized, but the ignorance of harmonic components reduces the accuracy of circuit analysis [5]. Extended describing function (EDF) modeling methods [6] or time-domain analysis methods [7] are adopted to solve such problems. A state trajectory analysis was proposed to achieve higher accuracy of resonant tank behavior control [8]. The resonant process was analyzed in the state plane, and the analytical formula can be obtained based on geometric relationships of the trajectory [9]. For series resonant converter (SRC), parallel resonant converter (PRC), and LLC resonant converters [10], state trajectory analysis and control have been applied in addressing unpredictable dynamics, burst mode for light load [11], soft startup, and short circuit protection [12]. However, more state variables make it more complicated for LCC resonant converters to present their relationships in a 2-D plane coordinate system. Therefore, it is very challenging to implement the state trajectory control. Furthermore, the startup duration is the toughest process of the state trajectory control [13]. This article proposes a state trajectory startup control and real-time implementation for ISOS Type LCC resonant converters with a capacitive output filter in a 140 kV/42 kW X-ray generator. The state trajectory startup control for LCC converter is optimized and conducted using FPGA to achieve low ripple voltage waveforms with fast-rising time under no-overshoot response. Furthermore, the state trajectory-based startup is successfully applied to an industry-level X-ray generator with a Power/ Voltage Level of 42 kw/140 kV for the first time. The rest of the paper is organized as follows. Section 2 introduces the state plane analysis method and state trajectory, including six operation conditions in detail, and describes the optimal trajectory control method for the soft startup process. Section 3 introduces the implementation of the state trajectory control in FPGA. Finally, the experimental results on a 140 kV/42 kW HV power supply used for the X-ray generator are presented to verify the proposed model and control strategy.

Analysis of State Trajectory-Based Startup for LCC Converter
The control mechanism for the LCC converter's state trajectory-based startup is depicted in Figure 3. The figure illustrates the topology of the LCC converter and the sampling and control processes used by the control system. The state trajectory-based startup for the LCC converter is divided into three stages [14].
This article proposes a state trajectory startup control and real-time implem for ISOS Type LCC resonant converters with a capacitive output filter in a 140k X-ray generator. The state trajectory startup control for LCC converter is optimi conducted using FPGA to achieve low ripple voltage waveforms with fast-ris under no-overshoot response. Furthermore, the state trajectory-based startup cessfully applied to an industry-level X-ray generator with a Power/ Voltage 42kw/140kV for the first time. The rest of the paper is organized as follows. Sect troduces the state plane analysis method and state trajectory, including six o conditions in detail, and describes the optimal trajectory control method for startup process. Section 3 introduces the implementation of the state trajectory c FPGA. Finally, the experimental results on a 140kV/42kW HV power supply use X-ray generator are presented to verify the proposed model and control strategy

Analysis of State Trajectory-Based Startup for LCC Converter
The control mechanism for the LCC converter's state trajectory-based startu picted in Figure 3. The figure illustrates the topology of the LCC converter and pling and control processes used by the control system. The state trajectory-based for the LCC converter is divided into three stages [14]. In Stage 1, the resonant current of the LCC converter rises to -imaxN in switching cycle. ±imaxN is the limit value of the resonant current in the state tr control, as shown in Figure 4. At the end of converter startup, the frequency val pensated by PI control is added to the frequency value predicted by state traject trol [15]. Vin and Vo are sampled for the calculation of state trajectory control scribed in Figure 4, the LCC converter works in the 3rd operating mode with S1 In the 3rd operating mode, the state trajectory is arc AB, and the state trajectory e is given as:  In Stage 1, the resonant current of the LCC converter rises to −i maxN in the first switching cycle. ±i maxN is the limit value of the resonant current in the state trajectory control, as shown in Figure 4. At the end of converter startup, the frequency value compensated by PI control is added to the frequency value predicted by state trajectory control [15]. V in and V o are sampled for the calculation of state trajectory control. As described in Figure 4, the LCC converter works in the 3rd operating mode with S1 and S4. In the 3rd operating mode, the state trajectory is arc AB, and the state trajectory equation is given as:  Then, the LCC converter operates in the 6th operating mode with S2 a state trajectory is arc BCD. The state trajectory of the 6th operating mode where, iLrN(t5) and VCprN(t5) are the initial values of iLrN(t) and VCprN(t) in the mode.
The arc radius of the state trajectories in the 3rd operating mode and ating mode are defined as ρT0 and ρT1. The central angle of arc AB is defined central angle of arc BCD is defined as θ1+θ2. The process of solving the cen be obtained as:  Figure 4. ρT0 and ρT1 can be expressed as: Then, the LCC converter operates in the 6th operating mode with S2 and S3, and the state trajectory is arc BCD. The state trajectory of the 6th operating mode is given as: where, i LrN (t 5 ) and V CprN (t 5 ) are the initial values of i LrN (t) and V CprN (t) in the 6th operating mode. The arc radius of the state trajectories in the 3rd operating mode and the 6th operating mode are defined as ρ T0 and ρ T1 . The central angle of arc AB is defined as θ 0 , and the central angle of arc BCD is defined as θ 1 + θ 2 . The process of solving the central angle can be obtained as: According to Figure 4. ρ T0 and ρ T1 can be expressed as: Therefore, the conduction time T 0 of S 1 and S 4 and T 1 of S 2 and S 3 can be derived as: In stage 2, the maximum values of resonant current are kept constant [16]. The energy in the resonator increases rapidly and is transmitted to the load. The output voltage of LCC converter rises quickly. The power in the resonator is low, and the resonant current value reaches the maximum at the switching time [17].
As shown in Figure 5, the state trajectory of the whole cycle is defined as 6 operating modes. In the 1st operating mode, the center coordinates of the state trajectory arc are (1, 0), and the radius is R [18]. In the 2nd operating mode, the state trajectory is an ellipse. In the 3rd operating mode, the state trajectory equation is the same as in stage 1. The trajectory equation of Mode I-III is given as: where, Z 0 is the double-element resonance impedance [19], Z 1 is the triple-element resonance impedance, i LrN (t 0 ) and V CprN (t 0 ) are the initial values of i LrN (t) and V CprN (t) in the 1st operating mode, i LrN (t 1 ) and V CprN (t 1 ) are the initial values of i LrN (t) and V CprN (t) in the 2nd operating mode.
orld Electr. Veh. J. 2022, 13, x FOR PEER REVIEW where, Z0 is the double-element resonance impedance [19], Z1 is the triple nance impedance, iLrN(t0) and VCprN(t0) are the initial values of iLrN(t) and VC operating mode, iLrN(t1) and VCprN(t1) are the initial values of iLrN(t) and VCpr operating mode. Similarly, according to the geometric symmetry of the trajectory, equation of Mode IV-VI can be obtained as:  Similarly, according to the geometric symmetry of the trajectory, the trajectory equation of Mode IV-VI can be obtained as: The coordinates of points A, B, C, D are substituted into the state trajectory equation, and the equation can be rewritten as: R is the radius of state trajectory arc in the 3rd operating mode, r is the radius of state trajectory arc in the 4th operating mode.
Through the analysis of the state trajectory in stage 2, and (12) can be simplified as: The parameters A 1 , B 1 , and C 1 can be expressed as (14).
The equation for solving the central angle and centrifugal angle of the state trajectory curve and the equation of the operating frequency can be obtained as: where, the base values of the circuit parament are listed in Table 1. By normalizing all voltages with the voltage factor V in and all currents with the current factor V in /Z 0 , the following expressions can be derived as [20]: In stage 3, the resonant current has reached the maximum value before switching. The energy in the resonator is high. The maximum resonant current set value is gradually reduced, and PI control is added.

Base Value
Similarly, according to the state trajectory in Figure 6, the coordinates of points E, F, G, H are substituted into the state trajectory equation, and the equation can be rewritten as: The parameters a1, b1, and c1 can be expressed by (19)  The solutions of Equation (17) can be simplified as: The parameters a 1 , b 1 , and c 1 can be expressed by (19) The equation for solving the center angle and centrifugal angle of the state trajectory curve and the equation of the operating frequency is:

Implementation of State Trajectory Control Based on FPGA
In stage 1, T0 and T1 are calculated by DSP, and then the data of T0 and T1 are sent to FPGA before the LCC converter starts. As shown in Figure 7 Figure 7 shows the calculation process of stage 2 in FPGA. As shown in Figure 7, A1, B1, and C1 can be calculated in the 1 st , 2 nd, and 3 rd calculation of multiplication correspondingly. θ0 can be calculated in the calculation of arcsine. θ1 can be calculated in the 3 rd calculation of arc cosine. θ2 is the difference between the 2nd and 1st calculation of arc cosine. Figure 8 shows the calculation process of stage 3 in FPGA. As shown in Figure 8, a1 can be calculated in the 1 st calculation of multiplication, and c1 is the sum of the 4 th and 6 th calculations. θ3+θ5 can be calculated when the 3 rd calculation of arc cosine is completed. θ4 can be calculated in the 4 th calculation of arc cosine.
Mult 10 Vo/Vin  In FPGA, division and square root operation modules consume many logic resources. The logic resources of FPGA are seriously wasted when the state trajectory control program is parallel computing. Therefore, the state trajectory control program adopts serial calculation. FPGA has only one division, square root, and inverse trigonometric function module. The red arrow shows the operation sequence of the program. Figure 7 shows the calculation process of stage 2 in FPGA. As shown in Figure 7, A 1 , B 1, and C 1 can be calculated in the 1st, 2nd, and 3rd calculation of multiplication correspondingly. θ 0 can be calculated in the calculation of arcsine. θ 1 can be calculated in the 3rd calculation of arc cosine. θ 2 is the difference between the 2nd and 1st calculation of arc cosine. Figure 8 shows the calculation process of stage 3 in FPGA. As shown in Figure 8, a 1 can be calculated in the 1st calculation of multiplication, and c 1 is the sum of the 4th and 6th calculations. θ 3 + θ 5 can be calculated when the 3rd calculation of arc cosine is completed. θ 4 can be calculated in the 4th calculation of arc cosine.

Div1
Vin Mult1 2×Cp/(n1)×(1/Cr+1/Cp) Mult 10 The data in FPGA is stored in 25-bit fixed-point register. The 25th bit of the register is used to store the sign of data. The 17th to 24th bits are used to store integer data. The 1st to 16th bits of the register are used to store decimal data. FPGA can not directly process decimal, so the CORDIC(Coordinate Rotation Digital Computer) algorithm is used in the division and square root operation module. The look-up table is used to operate the inverse trigonometric function.
The calculation speed of the state trajectory control in FPGA is significantly faster than that in DSP. The calculation time of state trajectory control in DSP28335 needs more than 20 us. But it just takes 200 clock cycles to complete all calculations of state trajectory control in FPGA, and the calculation time is less than 2us. In FPGA, each division and square root operation requires 17 clock cycles, and Multiplication and inverse trigonometric function operations require 4 clock cycles.

Experiment Results
In order to verify the correctness of the proposed control algorithm, the 140 kV/42 kW for the X-ray generator prototype platform was built up, as shown in Figure 9. The experiment parameters are listed in Table 2.

Load (RL) Transformer Ratio Resonant inductor (Lr) Series Resonant capacitor (Cr) Parallel Resonant capacitor (Cp)
Switching frequency (fsw) Figure 9. Prototype of the 140kV/42kW. Figure 10a, the exposure time of the And the output voltage Vo is 100 kV, imax is set as 200A, th the set value of 200A during the first cycle, and no extra injected into the resonant circuit. A fast rise time is achie  As shown in Figure 10a, the exposure time of the X-ray generator is set to 10 ms. And the output voltage V o is 100 kV, i max is set as 200 A, the resonant tank current achieves the set value of 200 A during the first cycle, and no extra voltage or current overshoot is injected into the resonant circuit. A fast rise time is achieved during the startup process, where (10-90%) of V o is 350 µs. At the same time, for conventional PFM control, the experiment waveforms at startup are presented in Figure 10b A fast rise time is achieved during the startup process. Based on the same circuit parameters and classical PI control, the rise time of the output voltage of the generator is 1.182 ms. where (10-90%) of Vo is 350 μs. At the same time, for conventional PFM control, the experiment waveforms at startup are presented in Figure 10b A fast rise time is achieved during the startup process. Based on the same circuit parameters and classical PI control, the rise time of the output voltage of the generator is 1.182 ms.

As shown in
(a) (b)  The switching process of different kV levels is presented between 140 kV, 120 kV, and 80 kV to achieve a faster transient time and keep the resonant tank's voltage or current from overshoot. The control strategy combines the benefits of OTC and PI. It has an excellent dynamic characteristic with a rise time of less than 100μs and a fall time of less than 150 μs under the resistive load, which benefits from the resonant tank being precisely controlled. Therefore, the validity of the proposed OTC is verified.  where (10-90%) of Vo is 350 μs. At the same time, for conventional PFM control, the experiment waveforms at startup are presented in Figure 10b A fast rise time is achieved during the startup process. Based on the same circuit parameters and classical PI control, the rise time of the output voltage of the generator is 1.182 ms.
(a) (b)  The switching process of different kV levels is presented between 140 kV, 120 kV, and 80 kV to achieve a faster transient time and keep the resonant tank's voltage or current from overshoot. The control strategy combines the benefits of OTC and PI. It has an excellent dynamic characteristic with a rise time of less than 100μs and a fall time of less than 150 μs under the resistive load, which benefits from the resonant tank being precisely controlled. Therefore, the validity of the proposed OTC is verified.

Conclusions
In this paper, a novel topology of ISOS Type LCC-SPRC is presented. The starting strategy of the LCC converter based on state trajectory control is analyzed in detail, and a state trajectory control based on FPGA for LCC resonant converter is proposed. Finally, the proposed approach has been implemented in FPGA and DSP. The experimental results on a 140 kV/42 kW X-ray imaging system generator show the model's good dynamics and the control scheme's implementation ability and significantly reduce the rising time of output voltage.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.