CN114552657A - Vehicle network system dynamic stability analysis method considering phase-locked loop angular frequency change - Google Patents

Vehicle network system dynamic stability analysis method considering phase-locked loop angular frequency change Download PDF

Info

Publication number
CN114552657A
CN114552657A CN202210317322.3A CN202210317322A CN114552657A CN 114552657 A CN114552657 A CN 114552657A CN 202210317322 A CN202210317322 A CN 202210317322A CN 114552657 A CN114552657 A CN 114552657A
Authority
CN
China
Prior art keywords
network system
angular frequency
phase
locked loop
vehicle network
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210317322.3A
Other languages
Chinese (zh)
Other versions
CN114552657B (en
Inventor
刘志刚
喻文倩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN202210317322.3A priority Critical patent/CN114552657B/en
Publication of CN114552657A publication Critical patent/CN114552657A/en
Application granted granted Critical
Publication of CN114552657B publication Critical patent/CN114552657B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
    • H02M7/00—Conversion of ac power input into dc power output; Conversion of dc power input into ac power output
    • H02M7/02—Conversion of ac power input into dc power output without possibility of reversal
    • H02M7/04—Conversion of ac power input into dc power output without possibility of reversal by static converters
    • H02M7/12—Conversion of ac power input into dc power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode
    • H02M7/145—Conversion of ac power input into dc power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a thyratron or thyristor type requiring extinguishing means
    • H02M7/155—Conversion of ac power input into dc power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a thyratron or thyristor type requiring extinguishing means using semiconductor devices only
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/40—Arrangements for reducing harmonics

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Electric Propulsion And Braking For Vehicles (AREA)

Abstract

The invention discloses a dynamic stability analysis method of a vehicle network system considering the angular frequency change of a phase-locked loop, which specifically comprises the following steps: firstly, a state equation of a circuit part of the vehicle network system is established through the simplified equivalent circuit. Then, dynamic equation description is carried out on the angular frequency in the phase-locked loop and the time delay in the controller, and a vehicle network system time domain state equation model considering the dynamic angular frequency is obtained; secondly, solving a balance point of the system according to the established SSA model; furthermore, the influence of circuit parameters and control parameters on the stability of the car network system is analyzed by solving the Jacobian matrix and by a root locus method; then, drawing a stable region of the parameter by a design-oriented analysis method; and finally, verifying the accuracy of the theoretical analysis result through simulation and a hardware-in-loop platform. The method improves the precision of the model, can reflect the influence of the circuit and the control parameter on the stability of the car network system, and has quick calculation and short time consumption.

Description

Vehicle network system dynamic stability analysis method considering phase-locked loop angular frequency change
Technical Field
The invention belongs to the field of high-speed railway time domain state equation modeling, and particularly relates to a dynamic stability analysis method for a vehicle network system by considering phase-locked loop angular frequency change.
Background
In recent years, four-quadrant voltage source converter motor train units with an AC-DC-AC structure have been widely adopted, and the main reason for the LFO phenomenon of the train network system is generally considered to be parameter mismatching between a traction network and the motor train units. In order to avoid the LFO phenomenon and ensure the reliable operation of the system, it is necessary to research the influence of the parameters on the stability of the vehicle network system.
At present, most documents mainly analyze the stability of the car network system by establishing a small signal impedance model 2 in a dq domain. Although this method is not difficult in each step of calculation, the frequency coupling process is omitted in the modeling process, so that the analysis result is only effective in the low frequency band, but the analysis result in the middle and high frequency bands has errors. Document 3 adopts a modeling mode based on harmonic linearization, considers the coupling process of partial frequencies, but ignores wp±nw1The frequency coupling term of (1). In document 4, an α β impedance model is derived based on a dq impedance model, and the critical value changes of various parameters when the converter switches between a rectification mode and an inversion mode are analyzed. The nature of the impedance model is to calculate the input and output impedance of the vehicle network system, and the operation process is complicated. Document 5 establishes an HSS model of a single-phase rectifier, and analyzes the dynamic characteristics of a harmonic transfer function using characteristic values. For the HSS modeling process, the higher the harmonic order considered, the more accurate the resulting model will be. Meanwhile, it is worth noting that due to the asymmetry of the single-phase rectifier, compared with a three-phase rectifier, the impedance modeling process of the single-phase rectifier is more difficult, and the calculation is more complicated.
In summary, the following steps:
1) most of modeling modes of the vehicle network system adopt impedance models in a frequency domain, and the modeling process is complex. Equation models in the time domain are less studied, and are described by an accurate mathematical model which takes the dynamic change of the phase-locked loop into consideration.
2) When stability analysis is performed, the traditional frequency domain analysis method is complex and takes a long time. And the characteristic value and design-oriented analysis method can more quickly and clearly show the result of the stability analysis.
Therefore, the time domain state space modeling method is introduced into the stability analysis of the vehicle network system, and meanwhile, the dynamic angular frequency in the phase-locked loop is considered, so that an accurate mathematical model is established, the influence of parameters on the system can be judged quickly and accurately, and the stable range of the parameters is provided.
Reference documents:
[1]：LIU Zhigang，ZHANG Guinan，LIAO Yicheng.Stability research ofhigh-speed railway EMUs andtraction network cascade system considering impedance matching[J].IEEE Transactions on Industry Applications，2016，52(5)：4315-4326.
[2]: zhang Guinan, Zhang Bo, jin, et al. study of inter-harmonic propagation mechanism of traction drive systems taking into account network pressure fluctuations [ J ] electric automation equipment, 2021, 41 (2): 186-192.
[3]：Zhang Han,Liu Zhigang,Wu Siqi,et al.Input impedance modeling and verification of single-phase voltage source converters based on harmonic linearization[J].IEEE Trans-actions on Power Electronics,2018,34(9):8544-8554.
[4]：HARNEFORS Lennart，YEPES Alejandro G.，Vidal Ana，et al.Passivity-based controller design ofgrid-connected VSCs forprevention ofelectrical resonance instability[J].IEEE Transactions on Industrial Electronics，2014，62(2)：702-710.
[5]：WANG Xiongfei，BLAABJERG Frede，LOH Poh Chiang.Passivity-based stability analysis and damping injection for multi-paralleled VSCs with LCL filters[J].IEEE Transactions on Power Electronics，2017，32(11)：8922-8935.
Disclosure of Invention
The invention aims to establish a time domain state space equation aiming at the low-frequency oscillation phenomenon of a vehicle network system, consider the dynamic angular frequency of a phase-locked loop, analyze the influence of a circuit and control parameters on the stability of the vehicle network system through a root track and a design-oriented method, draw the stable range of various parameters and provide a theoretical basis for the stable operation of the system. Therefore, the dynamic stability analysis method of the vehicle network system considering the change of the angular frequency of the phase-locked loop is provided.
The invention discloses a dynamic stability analysis method of a vehicle network system considering the angular frequency change of a phase-locked loop, which comprises the following steps of:
step 1: listing state equations of the circuit part according to the simplified equivalent circuit model; and constructing a vehicle network system time domain state space equation model considering the dynamic angular frequency by considering the change of the angular frequency in the phase-locked loop and the time delay part in the control and neglecting the secondary coupling frequency.
S11: state equation of the circuit part:
Figure BDA0003570230220000021
Figure BDA0003570230220000022
Figure BDA0003570230220000023
in the formula, edAnd eqIs the voltage of PCC on the dq axis, idAnd iqIs the grid side current in the dq coordinate system. u. ofdcRepresenting the output dc voltage. OmegaPLLIs the output angular frequency of the PLL.
Figure BDA0003570230220000024
And
Figure BDA0003570230220000025
is the output of the current inner loop controller. R is used for the equivalent resistance and inductance of the vehicle-mounted sidenAnd LnAnd (4) showing. CdRepresenting the capacitance on the DC side, RdRepresenting the equivalent resistance on the dc side.
S12: PCC node voltage enQ-axis voltage component, e, decoupled at dq and coordinate-transformedqThe expression of (c) is:
Figure BDA0003570230220000031
wherein L issAnd RsRespectively representing the equivalent inductance and the equivalent resistance, u, of the traction network sidesD、usQRespectively representing the supply voltage usThe components of the D axis and the Q axis in the DQ coordinate system on the traction net side; theta represents an included angle between a DQ coordinate system on the traction network side and a DQ coordinate system on the single-phase rectifier side, and the dynamic equation expression of theta/dt is omega-omegaPLL；
S13: the angular frequency dynamics in the phase locked loop are described as:
Figure BDA0003570230220000032
in the formula, kupllAnd kipllRespectively representing the proportional gain and the integral coefficient of a PI controller in a phase-locked loop; t represents time, ωkRepresenting the dynamic angular frequency in the phase locked loop.
S14: the time delay part, after being simplified by using three-order Pade approximation, is written into an expression of a state equation as follows:
Figure BDA0003570230220000033
yx＝Cxxx+Dxux
wherein x isx＝[x1,x2,x3,x4,x5,x6]TIs a new auxiliary state variable that is,
Figure BDA0003570230220000034
representing the time derivative of the auxiliary state variable,
Figure BDA0003570230220000035
is an output variable, ux＝[dd,dq]TIs an input variable, Ax，Bx，CxAnd DxIs a parameter matrix related to time delay, and the specific expression is：
Figure BDA0003570230220000036
Figure BDA0003570230220000037
Figure BDA0003570230220000038
Where τ denotes a time constant and I denotes an identity matrix.
S15: and synthesizing the expression to obtain a time domain state space equation of the single-phase vehicle network system as follows:
Figure BDA0003570230220000039
wherein X ═ id,iq,udc,mid,miq,mdc,θ,ωPLL,xx]TWhich represents the state variable of the device,
Figure BDA0003570230220000041
represents an input variable, u ═ usD,usQ,Udcref]TRepresenting the control variable. m isid、miqAnd mdcAre all auxiliary state variables in the controller, usDAnd usQRespectively representing the supply voltage usThe D-axis and Q-axis components in the DQ coordinate system on the trailing wire side.
Step 2: and according to the established time domain state space equation model, solving the balance point of the system by using a Newton iteration method.
And step 3: and solving the Jacobian matrix of the car network system, and analyzing the influence of the circuit and the control acceptance number on the stability of the car network system according to the root track.
And 4, step 4: and drawing a stable region of the parameter by a design-oriented analysis method on the basis of the established model.
And 5: and verifying the theoretical analysis result on a Simulink simulation and hardware-in-loop platform.
The beneficial technical effects of the invention are as follows:
1. according to the method, the dynamic angular frequency and the time delay part of the phase-locked loop are brought into the modeling process of the time domain state space model, so that the model precision is improved, and the influence of the circuit and the control parameters on the stability of the vehicle network system can be reflected.
2. The modeling process and the analyzing process are both simple, signals do not need to be converted from a time domain to a frequency domain in the modeling process, and the characteristic value and design-oriented analyzing method is rapid in calculation and short in time consumption.
The invention can well analyze the influence of the circuit and control parameters on the low-frequency oscillation phenomenon, reveal the generation reason and mechanism of the phenomenon and be beneficial to making corresponding inhibition measures aiming at the phenomenon by considering the dynamic angular frequency and time delay of the phase-locked loop to establish a time domain state space vehicle network system model after the simulation and hardware verification.
Drawings
Fig. 1 is a model of a vehicle network system of the present invention.
Fig. 2 is a control block diagram of the circuit of the present invention.
FIG. 3 shows a network side inductor LsSampling time TsDC side capacitor CdProportional gain k of sum current loopipRoot trace map of changes.
FIG. 4 shows the network side inductor LsSampling time TsDC side capacitor CdProportional gain k of sum current loopipAnd (5) a low-frequency oscillation simulation result graph.
FIG. 5 shows the network side inductor LsSampling time TsDC side capacitor CdProportional gain k of sum current loopipAnd (5) a low-frequency oscillation semi-physical experiment result graph.
FIG. 6 shows inductance values L for different net side inductancessWhen T issWhen it is 5e-5s, kip_CdOn a planeThe boundary is stabilized.
FIG. 7 shows the sampling time T for different timessWhen C is presentdWhen equal to 4.5mF, Ls_kipStable boundary on the plane.
FIG. 8 shows the capacitance C for different DC-side valuesdWhen L is presentsWhen 2mH, kip_TsA stable boundary on the face.
FIG. 9 shows the inner loop proportional gain k for different currentsipWhen C is presentdWhen the average grain size is 4.5mF, Ts_LsStable boundary on the plane.
Detailed Description
The invention is described in further detail below with reference to the figures and the specific embodiments.
The example takes a simplified single-phase vehicle network system as an example. Fig. 1 is a model of a single-phase vehicle network system, which is composed of a simplified traction network equivalent circuit and a vehicle-mounted single-phase rectifier equivalent circuit. Wherein u issRepresenting equivalent power supply on the side of the traction network, LsAnd RsRespectively representing the equivalent inductance and the equivalent resistance of the traction network side. The voltage at PCC is denoted as en。inRepresenting the current on the net side. R is used for the equivalent resistance and inductance of the vehicle-mounted sidenAnd LnAnd (4) showing. CdRepresenting the capacitance on the DC side, RdRepresenting the equivalent resistance on the dc side.
FIG. 2 is a control block diagram of a circuit, where UdcrefAnd eqrefRespectively representing the reference value of the DC side voltage and the voltage e at PCCnReference value of q-axis voltage after dq decomposition. I isdrefAnd IqrefRepresenting the d-axis and q-axis current reference values of the net side current after dq decomposition, respectively. k is a radical ofupAnd kuiRespectively representing the proportional gain and the integral coefficient of the voltage outer loop PI controller. k is a radical ofipAnd kiiRespectively representing the proportional gain and the integral coefficient in the current inner loop PI controller. k is a radical ofupllAnd kipllRespectively representing the proportional gain and the integral coefficient of the PI controller in the phase locked loop. The parameters of the vehicle network system are shown in table 1.
TABLE 1 vehicle network System model parameters
Figure BDA0003570230220000051
The root trajectory, the simulation result, the semi-physical experiment result, and the design-oriented analysis result of this embodiment are composed of the following steps.
1) By varying the net side inductance L in the built modelsSampling time TsDC side capacitor CdProportional gain k of sum current loopipThe method comprises the steps of (1) obtaining a Jacobian matrix by carrying out linearization processing on a model at a balance point, further obtaining a characteristic value of a system, and when a real part of a certain characteristic value crosses a virtual axis by observing the change of the characteristic value of the system, the system is unstable;
FIGS. 3(a) - (d) are respectively the network side inductors LsSampling time TsDC side capacitor CdProportional gain k of sum current loopipThe root trajectory change map. Current network side inductance LsAt 4.2mH, the real part of the maximum eigenvalue crosses the imaginary axis, and the system will be unstable, and the frequency of the system oscillation is 2.1 Hz. In the same way, it can be observed that when the sampling time T issIs 2e-4s, and a DC side capacitor CdIs 8mF and the current loop proportional gain kipAt 1.3, the real part of the eigenvalue crosses to the right half plane from the left half plane, which means that the system is unstable and the oscillation frequency is 3.1Hz, 2.1Hz and 2.2 Hz.
2) By changing the inductance L at the network side in Matlab \ sillink simulationsSampling time TsDC side capacitor CdProportional gain k of sum current loopipThe oscillation waveform of the system can be observed;
FIGS. 4(a) - (d) are respectively the network side inductances LsEqual to 4.3mH, sample time TsThe capacitance C on the direct current side is increased to 2e-4sdIncrease to 8mF and current loop proportional gain kipThe simulation waveform is reduced to 1.4, and the system has low-frequency oscillation phenomena with frequencies of 2.1Hz, 3Hz, 2Hz and 2.3Hz respectively. Can seeIt is shown that the critical value of the parameter and the oscillation frequency result are very close compared to the theoretical analysis result.
3) Performing a test on a semi-physical platform, and observing the influence of a circuit and control on the stability of a vehicle network system;
in order to further verify the theoretical analysis result, a semi-physical experiment platform is set up. A model of a vehicle network system is divided into two parts, namely a controlled object (a main circuit) and a control algorithm, the two parts are respectively downloaded into a hardware-in-loop real-time simulation system and a rapid prototype control system by using StarSim software and are connected through a physical input and output module to form a closed-loop test circuit.
FIGS. 5(a) - (d) are respectively the network side inductances LsIncrease to 4.5mH, sample time TsThe capacitance C on the direct current side is increased to 2.2e-4sdIncrease to 8mF and current loop proportional gain kipWhen the voltage is reduced to 1.3, the waveforms of the ac side voltage current and the dc side voltage are extracted from the oscilloscope. The experimental result on the semi-physical platform is similar to the Matlab \ Simulink simulation result.
Specific parameter critical values and oscillation frequencies between theoretical analysis results, Matlab \ Simulink simulation results and experimental results on a semi-physical platform are shown in tables 2 and 3.
TABLE 2 parameter thresholds
Figure BDA0003570230220000061
TABLE 3 Low frequency oscillation frequency
Figure BDA0003570230220000062
Figure BDA0003570230220000071
It can be seen that the theoretical analysis result, the Matlab \ Simulink simulation result and the semi-physical platform experiment result are basically consistent.
4) Drawing the stable region of the parameter by adopting a design-oriented analysis method
First, inductance L is measured from the network sidesSampling time TsDC side capacitor CdProportional gain k of sum current loopipSelecting 3 of them to form a parameter plane, and then determining that one variable of the 3 is temporarily used as a fixed value, one is used as an independent variable, and the other parameter is used as a numerical value of the ordinate, and calculating the stable critical value by using the above eigenvalue analysis method, and recording the stable critical value. And respectively solving corresponding critical values of the ordinate parameters by changing the values of the independent variable parameters. And finally, drawing the recorded critical points into a stable boundary curve through least square fitting. The results are shown in FIGS. 6-9.
As can be seen from fig. 6, when the sampling time is constant, the larger the net side inductance is, the smaller the stable region is. Fig. 7 illustrates that the larger the sampling time, the smaller the stable region when the dc-side capacitance is constant. It can be observed from fig. 8 that for different dc side capacitances, the larger the capacitance, the smaller the stable region when the grid side inductance parameter is constant. Fig. 9 shows that when the proportional gain of the current loop is changed and the capacitance parameter is not changed, the larger the gain is, the larger the stable region is.

Claims (2)

1. A dynamic stability analysis method of a vehicle network system considering phase-locked loop angular frequency change is characterized by comprising the following steps:
step 1: listing state equations of the circuit part according to the simplified equivalent circuit model; by considering the change of angular frequency in a phase-locked loop and a time delay part in control and neglecting secondary coupling frequency, a vehicle network system time domain state space equation model considering dynamic angular frequency is constructed;
step 2: according to the established time domain state space equation model, a Newton iteration method is used for solving the balance point of the system;
and step 3: solving a Jacobian matrix of the car network system, and analyzing the influence of the circuit and the control acceptance number on the stability of the car network system according to the root track;
and 4, step 4: on the basis of the established model, drawing a stable region of the parameter by a design-oriented analysis method;
and 5: and verifying the theoretical analysis result on a Simulink simulation and hardware-in-loop platform.
2. The method for analyzing the dynamic stability of the vehicle network system in consideration of the angular frequency change of the phase-locked loop according to claim 1, wherein the step 1 specifically comprises:
s11: state equation of the circuit part:
Figure FDA0003570230210000011
Figure FDA0003570230210000012
Figure FDA0003570230210000013
in the formula, edAnd eqIs the voltage of PCC on the dq axis, idAnd iqIs the grid side current in dq coordinate system, udcIndicating the output DC voltage, ωPLLIs the output angular frequency of the PLL and,
Figure FDA0003570230210000014
and
Figure FDA0003570230210000015
is the output of the current inner loop controller; r is used for the equivalent resistance and inductance of the vehicle-mounted sidenAnd LnIs represented by CdCapacitance representing the DC side, RdRepresents the equivalent resistance on the direct current side;
s12: PCC node voltage enQ-axis voltage component, e, decoupled at dq and coordinate-transformedqThe expression of (a) is:
Figure FDA0003570230210000016
wherein L issAnd RsRespectively representing the equivalent inductance and the equivalent resistance, u, of the traction network sidesD、usQRespectively representing the supply voltage usThe components of the D axis and the Q axis in the DQ coordinate system on the traction net side; theta represents an included angle between a DQ coordinate system on the traction network side and a DQ coordinate system on the single-phase rectifier side, and the dynamic equation expression of theta/dt is omega-omegaPLL；
S13: the angular frequency dynamics in the phase locked loop are described as:
Figure FDA0003570230210000017
in the formula, kupllAnd kipllRespectively representing the proportional gain and the integral coefficient of a PI controller in a phase-locked loop; t represents time, ωkRepresenting the dynamic angular frequency in the phase-locked loop;
s14: the time delay part, after being simplified by using three-order Pade approximation, is written into an expression of a state equation as follows:
Figure FDA0003570230210000021
yx＝Cxxx+Dxux
wherein x isx＝[x1,x2,x3,x4,x5,x6]TIs a new auxiliary state variable that is,
Figure FDA0003570230210000022
representing the time derivative of the auxiliary state variable,
Figure FDA0003570230210000023
is an output variable, ux＝[dd,dq]TIs an input variable, Ax，Bx，CxAnd DxIs a parameter matrix related to time delay, and the specific expression is as follows:
Figure FDA0003570230210000024
Figure FDA0003570230210000025
Figure FDA0003570230210000026
wherein τ represents a time constant, and I represents an identity matrix;
s15: and synthesizing the expression to obtain a time domain state space equation of the single-phase vehicle network system as follows:
Figure FDA0003570230210000027
wherein X ═ id,iq,udc,mid,miq,mdc,θ,ωPLL,xx]TWhich represents the state variable of the device,
Figure FDA0003570230210000028
represents an input variable, u ═ usD,usQ,Udcref]TRepresents a control variable; m isid、miqAnd mdcAre all auxiliary state variables in the controller, usDAnd usQRespectively representing the supply voltage usThe D-axis and Q-axis components in the DQ coordinate system on the trailing wire side.
CN202210317322.3A 2022-03-29 2022-03-29 Dynamic stability analysis method for vehicle network system considering phase-locked loop angular frequency change Active CN114552657B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210317322.3A CN114552657B (en) 2022-03-29 2022-03-29 Dynamic stability analysis method for vehicle network system considering phase-locked loop angular frequency change

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210317322.3A CN114552657B (en) 2022-03-29 2022-03-29 Dynamic stability analysis method for vehicle network system considering phase-locked loop angular frequency change

Publications (2)

Publication Number Publication Date
CN114552657A true CN114552657A (en) 2022-05-27
CN114552657B CN114552657B (en) 2023-05-05

Family

ID=81664978

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210317322.3A Active CN114552657B (en) 2022-03-29 2022-03-29 Dynamic stability analysis method for vehicle network system considering phase-locked loop angular frequency change

Country Status (1)

Country Link
CN (1) CN114552657B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115102191A (en) * 2022-07-05 2022-09-23 兰州交通大学 Stability analysis method for photovoltaic access traction power supply system
CN116914757A (en) * 2023-07-20 2023-10-20 国网上海市电力公司 Method and device for determining maximum estimated attraction domain of flexible interconnection system of power distribution area

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012163355A1 (en) * 2011-05-31 2012-12-06 Vestas Wind Systems A/S System and method for generating an inertial response to a change in the voltage of an electricial grid
CN105006834A (en) * 2015-06-10 2015-10-28 合肥工业大学 Optimal virtual inertia control method based on virtual synchronous generator
CN108988384A (en) * 2018-07-24 2018-12-11 西安工业大学 Grid-connected current DC component suppressing method based on fractional order PIR
CN110020510A (en) * 2019-06-06 2019-07-16 西南交通大学 A kind of method for analyzing stability of the EMU net side Pulse rectifier based on Floquet theory
CN110138241A (en) * 2019-06-06 2019-08-16 西南交通大学 A kind of method for analyzing stability of the single-phase electricity potential source based on Jacobi's theory

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012163355A1 (en) * 2011-05-31 2012-12-06 Vestas Wind Systems A/S System and method for generating an inertial response to a change in the voltage of an electricial grid
CN105006834A (en) * 2015-06-10 2015-10-28 合肥工业大学 Optimal virtual inertia control method based on virtual synchronous generator
CN108988384A (en) * 2018-07-24 2018-12-11 西安工业大学 Grid-connected current DC component suppressing method based on fractional order PIR
CN110020510A (en) * 2019-06-06 2019-07-16 西南交通大学 A kind of method for analyzing stability of the EMU net side Pulse rectifier based on Floquet theory
CN110138241A (en) * 2019-06-06 2019-08-16 西南交通大学 A kind of method for analyzing stability of the single-phase electricity potential source based on Jacobi's theory

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115102191A (en) * 2022-07-05 2022-09-23 兰州交通大学 Stability analysis method for photovoltaic access traction power supply system
CN115102191B (en) * 2022-07-05 2024-04-02 兰州交通大学 Stability analysis method for photovoltaic access traction power supply system
CN116914757A (en) * 2023-07-20 2023-10-20 国网上海市电力公司 Method and device for determining maximum estimated attraction domain of flexible interconnection system of power distribution area
CN116914757B (en) * 2023-07-20 2024-05-07 国网上海市电力公司 Method and device for determining maximum estimated attraction domain of flexible interconnection system of power distribution area

Also Published As

Publication number Publication date
CN114552657B (en) 2023-05-05

Similar Documents

Publication Publication Date Title
CN109446653B (en) Modeling method of mixed MMC small-signal model of mixed MMC state space model
CN114552657A (en) Vehicle network system dynamic stability analysis method considering phase-locked loop angular frequency change
CN110739678A (en) Control method for series virtual impedance of grid-connected converters
WO2022226709A1 (en) Power control-based method for analyzing stability of permanent magnet synchronous fan accessing weak grid
CN112993998B (en) System harmonic stability judgment method based on LCC closed-loop impedance model
CN110224431B (en) Control method for reducing influence of phase-locked loop in grid-connected inverter system
CN111884218B (en) Stability evaluation method and system for double-fed VSC power transmission system
CN109782173A (en) Asynchronous machine excitation mutual inductance curve measurement system and its measurement method
CN112039123A (en) Control method for grid-connected inverter without alternating voltage sensor
Xue et al. A comprehensive study on impedance models of grid-tied voltage-source converters
CN114935690A (en) Disturbance injection and impedance measurement method and system suitable for new energy grid-connected system impedance measurement
CN112952901B (en) Distributed stability analysis method for multi-fan grid-connected system
Zhong et al. CDSC-based adaptive impedance measurement method for grid-tied inverter system under adverse grid voltage conditions
CN104993487A (en) SAPF control method for voltage source type non-linear load harmonic compensation
CN107294096B (en) The determination unit function coefficient method of AC harmonic filter is adjusted in a kind of consideration more
CN110912182B (en) Design method of virtual resistance parameters of parallel impedance adapter
CN111308207A (en) Dq impedance measuring method for single-phase alternating current system
CN112564131B (en) LCC-HVDC system power oscillation damping suppression method and system
CN115236404A (en) Self-measuring method for port impedance of grid-connected inverter
CN115378025A (en) Wind power converter grid-connected control method, device and medium without phase-locked loop
CN110095654B (en) Power grid inductance detection method
CN113189532A (en) Method and device for online correction of harmonic measurement error of capacitor voltage transformer
CN111106618A (en) Harmonic analysis method and device for new energy power generation equipment to access power system
CN113922390B (en) Vehicle network low-frequency oscillation suppression method based on improved phase-locked loop
CN110190627B (en) Grid-connected inverter power grid adaptability control algorithm based on auxiliary link

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant