CN114721276A - Transmission chain system collaborative modeling and multi-physical field analysis method - Google Patents

Transmission chain system collaborative modeling and multi-physical field analysis method Download PDF

Info

Publication number
CN114721276A
CN114721276A CN202210644031.5A CN202210644031A CN114721276A CN 114721276 A CN114721276 A CN 114721276A CN 202210644031 A CN202210644031 A CN 202210644031A CN 114721276 A CN114721276 A CN 114721276A
Authority
CN
China
Prior art keywords
chain system
transmission chain
model
control
electromagnetic
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
CN202210644031.5A
Other languages
Chinese (zh)
Other versions
CN114721276B (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.)
Hunan University
Original Assignee
Hunan 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 Hunan University filed Critical Hunan University
Priority to CN202210644031.5A priority Critical patent/CN114721276B/en
Publication of CN114721276A publication Critical patent/CN114721276A/en
Application granted granted Critical
Publication of CN114721276B publication Critical patent/CN114721276B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
    • 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
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

Compared with the traditional transmission chain multi-physical field analysis modeling method, the main advantages are that the influence of the transmission chain converter side on the whole electromagnetism, temperature and stress field of the transmission chain system is considered, through the establishment of a drive chain system control algorithm and a drive external circuit mathematical model, the electromagnetic loss of the generator under the PWM waveform is calculated, and based on the electromagnetic loss, performing electromagnetic-temperature field iterative calculation on the transmission chain system, taking the obtained temperature distribution result of the transmission chain system as the initial condition of a stress field after the precision requirement is met, calculating the thermal stress distribution of the transmission chain system, the invention fully considers the influence of the converter side on the transmission chain system, carries out cooperative modeling and multi-physical field analysis on the traditional chain system (converter-generator), and the modeling process and result thereof are more in line with the actual structure and the operation condition of the transmission chain system.

Description

Transmission chain system collaborative modeling and multi-physical field analysis method
Technical Field
The invention relates to the technical field of transmission chains, in particular to a transmission chain system collaborative modeling and multi-physical field analysis method.
Background
The deep and offshore wind generating set transmission chain system has the characteristics of large capacity, severe operating environment, huge maintenance cost and the like, and the high reliability of the transmission chain system becomes the key point for designing the deep and offshore large transmission chain system; the operation environment is considered to be severe, and the test cost is high; after the design is finished, high-precision and high-reliability modeling analysis is carried out on the transmission chain system to become a key link for the overall optimal design and manufacture of the transmission chain system; therefore, aiming at the real-time running condition of the system, the influence of the control system and the strong coupling relation among various physical fields are considered, the modeling analysis of the system is of great significance to the design and manufacture of high-reliability and high-efficiency transmission chain systems in China, and the basis for realizing the double-carbon goal in China is laid.
The deep and open sea oriented transmission chain system mainly comprises a converter-generator, wherein the converter plays an important role in modeling analysis of the whole transmission chain system, and different control algorithms of the converter can generate great influence on the overall electromagnetic performance of the motor, so that the indexes of the thermal balance state and the mechanical performance of the material of the motor are changed. In order to ensure the reliability of an offshore wind power transmission chain system, people develop modeling analysis and optimization design on a generator part, however, accurate modeling analysis is the basis for realizing optimization design, and people develop a large amount of modeling analysis on the generator part, so that the reliability of the transmission chain system is verified to a certain extent, but the following defects still exist:
(1) the real-time running state of the transmission chain system is complex, the influence of a converter side in the transmission chain system on the coupling relation of the whole multi-physical field is not considered, or a control system of a generator in the transmission chain system is simplified, and the generator control system is considered as an ideal three-phase sine system during modeling.
Secondly, in the prior art, when analyzing the reliability of the drive train system, only the generator side in the drive train system is usually considered, which only reveals the multi-physical field coupling relationship of the generator side in the drive train system, however, in the actual operation process, the drive train system operates as a whole, and all the physical fields between all the components are mutually coupled and interacted. Considering only the physical field performance on the motor side, it is difficult to evaluate and detect the reliability of the entire drive train system.
The side part of the converter in the traditional chain system comprises a driving external circuit and a control strategy, and the driving external circuit and the control strategy are used for driving the converter to work by outputting PWM pulse waveforms with different amplitudes and frequencies to the motor, compared with the three-phase sine wave shape under the ideal condition, the whole electromagnetic loss of the transmission chain is increased, the loss of the power element and the motor side is increased, the temperature is increased, and the mechanical stress is further increased, therefore, when the reliability of the drive chain system is verified, the whole structure of the drive chain needs to be modeled and analyzed, the traditional method only considers the multi-physical field distribution of the motor in the drive chain system, the influence of the converter side is ignored, the obtained model is often difficult to represent the performance of the whole system of the transmission chain, and when the multi-physical field analysis is carried out on the model, compared with the real result, the result obtained by the electromagnetic-temperature iterative calculation is larger in error, and the stress distribution is larger in error than that in the actual operation state.
Disclosure of Invention
The invention aims to solve the technical problem of providing a method for collaborative modeling and multi-physical field analysis of a transmission chain system.
In order to solve the technical problem, the invention adopts the following technical scheme:
a method for collaborative modeling and multi-physical field analysis of a drive chain system comprises the following steps;
s1, establishing a drive chain system control model, wherein the drive chain system control model comprises a control algorithm model and a control external circuit model;
s2, after receiving a driving signal provided by a control algorithm in a control model, a driving circuit power tube generates a periodic PWM control signal, an electromagnetic finite element model of the generator is established on the basis of considering the PWM control signal, namely a control equation generated by the control model in the simultaneous S1 is established, and the electromagnetic parameters of the motor are solved by establishing a magnetic vector equation similar to a steady field;
wherein, the electromagnetic finite element model in the cell is as follows:
Figure 562993DEST_PATH_IMAGE001
；（1）
in the formula
Figure 346142DEST_PATH_IMAGE002
Is in the x-axis direction;
Figure 809746DEST_PATH_IMAGE003
is in the y-axis direction;
Figure 641436DEST_PATH_IMAGE004
is an initial value;
Figure 5421DEST_PATH_IMAGE005
is the reluctance ratio;
Figure 693892DEST_PATH_IMAGE006
is a magnetic vector, wherein the magnetic vector has only a component in the direction of the z-axis;
Figure 376283DEST_PATH_IMAGE007
is the source current density;
Figure 542823DEST_PATH_IMAGE008
is the tangential component of the magnetic field strength;
Figure 26894DEST_PATH_IMAGE009
is a first type boundary;
Figure 886265DEST_PATH_IMAGE010
is a second type boundary;
s3, calculating the electromagnetic parameters of the drive chain system under the condition of PWM control waveform, establishing a temperature field model of the drive chain system, calculating a temperature distribution result by taking electromagnetic loss as a heat source in the temperature field model of the drive chain system, updating the electromagnetic parameters of the drive chain system, and performing iterative calculation until the temperature distribution error is less than 5%;
s4, taking the stable temperature distribution result as the initial condition of the stress field differential equation model of the transmission chain system, and calculating the thermal stress distribution of the transmission chain system;
s5: and outputting the structural model results of the calculated electromagnetic loss, temperature distribution diagram and stress distribution diagram, then comparing the results with the design indexes of the transmission chain system, and finally checking the reliability of the transmission chain.
Further, the control algorithm model is as follows;
Figure 760680DEST_PATH_IMAGE011
；（2）
Figure 966796DEST_PATH_IMAGE012
；（3）
in the formula (I), the compound is shown in the specification,
Figure 39794DEST_PATH_IMAGE013
、
Figure 70067DEST_PATH_IMAGE014
being the d-q axis component of the stator voltage,
Figure 494095DEST_PATH_IMAGE015
、
Figure 506831DEST_PATH_IMAGE016
is the d-q component of the stator current,
Figure 699915DEST_PATH_IMAGE017
is a resistance of the stator, and is,
Figure 104352DEST_PATH_IMAGE018
is the electrical angular velocity of the beam of light,
Figure 15676DEST_PATH_IMAGE019
、
Figure 530971DEST_PATH_IMAGE020
is the inductive component of the d-q axis,
Figure 47403DEST_PATH_IMAGE021
is a permanent magnet flux linkage, and is provided with a permanent magnet,
Figure 186523DEST_PATH_IMAGE022
is time.
Further, the external control circuit model is as follows;
Figure 522826DEST_PATH_IMAGE023
；（4）
Figure 904129DEST_PATH_IMAGE024
；（5）
in the formula (I), the compound is shown in the specification,
Figure 275067DEST_PATH_IMAGE025
is an electromotive force of the x-phase,
Figure 286886DEST_PATH_IMAGE026
is the voltage drop of N with respect to point O of reference,
Figure 844906DEST_PATH_IMAGE027
is the current of the x-phase,
Figure 731697DEST_PATH_IMAGE028
in order to be the load current,
Figure 222721DEST_PATH_IMAGE029
in order to be an equivalent resistance, the resistance,
Figure 405441DEST_PATH_IMAGE030
in order to be an equivalent inductance,
Figure 981916DEST_PATH_IMAGE031
in order to be an equivalent capacitance,
Figure 642704DEST_PATH_IMAGE032
is a direct-current voltage, and the voltage is,
Figure 988235DEST_PATH_IMAGE033
are respectively paired
Figure 341856DEST_PATH_IMAGE034
The three phases are calculated and,
Figure 140048DEST_PATH_IMAGE035
is time;
and (5) synthesizing the control models (2) to (5) to obtain the control model of the transmission chain system.
Further, in the S2;
Figure 371571DEST_PATH_IMAGE036
；（6）
in the formula (I), the compound is shown in the specification,
Figure 571608DEST_PATH_IMAGE037
、
Figure 361710DEST_PATH_IMAGE038
is the d-q component of the stator current;
the magnetic vector position distribution in the electromagnetic field of the drive chain system is influenced by a control algorithm and a driving external circuit in a control model, and the electromagnetic field performance under the control of the PWM waveform is considered in the formula (6).
Further, in S3, the temperature field model in the cell is:
Figure 381618DEST_PATH_IMAGE039
； (7)
in the formula (I), the compound is shown in the specification,
Figure 915368DEST_PATH_IMAGE040
is the temperature of the boundary of the motor,
Figure 969912DEST_PATH_IMAGE041
is the temperature of the fluid near the motor boundary,
Figure 665335DEST_PATH_IMAGE042
and
Figure 438119DEST_PATH_IMAGE043
is an electric motor
Figure 713243DEST_PATH_IMAGE044
And
Figure 911706DEST_PATH_IMAGE045
the thermal conductivity of the material in the direction,
Figure 246873DEST_PATH_IMAGE046
is the sum of the heat flux densities in the motor,
Figure 241373DEST_PATH_IMAGE047
in order to obtain a heat transfer coefficient,
Figure 116925DEST_PATH_IMAGE048
to pass through the boundary x1The density of the heat flow is such that,
Figure 146061DEST_PATH_IMAGE049
in order to obtain the magnetic induction intensity,
Figure 183288DEST_PATH_IMAGE050
as a matter of time, the time is,
Figure 930664DEST_PATH_IMAGE051
is in the direction of the x-axis,
Figure 642530DEST_PATH_IMAGE052
in the y-axis direction.
Further, in S4, the differential equation model of the stress field in the cell is:
Figure 526172DEST_PATH_IMAGE053
；（8）
in the formula (I), the compound is shown in the specification,
Figure 937562DEST_PATH_IMAGE054
、
Figure 703393DEST_PATH_IMAGE055
is composed of
Figure 124010DEST_PATH_IMAGE056
、
Figure 862159DEST_PATH_IMAGE057
The direction of the positive stress is the direction of the positive stress,
Figure 506767DEST_PATH_IMAGE058
in order to be able to apply shear forces in the respective planes,
Figure 963156DEST_PATH_IMAGE059
、
Figure 718622DEST_PATH_IMAGE060
in order to be a surface force,
Figure 75392DEST_PATH_IMAGE061
in the direction of the x-axis,
Figure 218797DEST_PATH_IMAGE062
in the y-axis direction.
Compared with the prior art, the invention has the following beneficial effects:
1. compared with the traditional transmission chain multi-physical field analysis modeling method, the method has the main advantages that the influence of the transmission chain converter side on the whole electromagnetic, temperature and stress fields of the transmission chain system is considered, the electromagnetic loss of the generator under the PWM waveform is calculated through the establishment of a transmission chain system control algorithm and a driving external circuit mathematical model, the transmission chain system is subjected to electromagnetic-temperature field iterative calculation on the basis of the electromagnetic loss and the temperature loss, and after the precision requirement is met, the obtained transmission chain system temperature distribution result is used as the initial condition of the stress field to calculate the thermal stress distribution. The invention fully considers the influence of the converter side on the transmission chain system, carries out cooperative modeling and multi-physical field analysis on the traditional system (converter-generator), and the modeling process and the result thereof are more in line with the actual structure and the operation condition of the transmission chain system.
2. The invention establishes a cooperative computing and analyzing platform of the transmission chain system, and fully considers the influence of a control module of the transmission chain system on the overall performance through the steps of building a control external circuit, realizing a control algorithm and the like; secondly, a multi-physical-field calculation model of the transmission chain system electromagnetism-temperature-stress is established through an electromagnetic model interface, temperature rise distribution caused by the electromagnetic loss of the transmission chain and stress distribution of the transmission chain system are calculated under the PWM power supply waveform, the method comprehensively considers the coupling relation among all parts of the transmission chain system and all physical fields in all the parts, the multi-physical-field model of the transmission chain system is accurately established, and the transmission chain system performance under different working conditions and running states is accurately calculated.
Detailed Description
To make the technical problems, technical solutions and advantages to be solved by the present invention clearer. In the following description, characteristic details such as specific configurations and components are provided only to help the embodiments of the present invention be fully understood. Thus, it will be apparent to those skilled in the art that various changes and modifications may be made to the embodiments described herein without departing from the scope and spirit of the invention. In addition, descriptions of well-known functions and constructions are omitted for clarity and conciseness.
It should be appreciated that reference throughout this specification to "one embodiment" or "an embodiment" means that a particular feature, structure or characteristic described in connection with the embodiment is included in at least one embodiment of the present invention. Thus, the appearances of the phrases "in one embodiment" or "in an embodiment" in various places throughout this specification are not necessarily all referring to the same embodiment. Furthermore, the particular features, structures, or characteristics may be combined in any suitable manner in one or more embodiments.
In various embodiments of the present invention, it should be understood that the sequence numbers of the following processes do not mean the execution sequence, and the execution sequence of each process should be determined by its function and inherent logic, and should not constitute any limitation to the implementation process of the embodiments of the present invention.
It should be understood that the term "and/or" herein is only one kind of association relationship describing the association object, and means that there may be three kinds of relationships, for example, a and/or B, and may mean: a exists alone, A and B exist simultaneously, and B exists alone. In addition, the character "/" herein generally indicates that the former and latter related objects are in an "or" relationship.
In the embodiments provided herein, it should be understood that "B corresponding to a" means that B is associated with a from which B can be determined. It should also be understood that determining B from a does not mean determining B from a alone, but may be determined from a and/or other information.
The invention provides a method for collaborative modeling and multi-physical field analysis of a transmission chain system, which comprises the following steps;
s1, establishing a drive chain system control model, wherein the drive chain system control model comprises a control algorithm model and a control external circuit model, and the control algorithm model is as follows;
Figure 428062DEST_PATH_IMAGE011
；（2）
Figure 721640DEST_PATH_IMAGE012
；（3）
in the formula (I), the compound is shown in the specification,
Figure 434381DEST_PATH_IMAGE013
、
Figure 922256DEST_PATH_IMAGE014
being the d-q axis component of the stator voltage,
Figure 353238DEST_PATH_IMAGE015
、
Figure 450507DEST_PATH_IMAGE016
is the d-q component of the stator current,
Figure 17754DEST_PATH_IMAGE017
is a resistance of the stator, and is,
Figure 175066DEST_PATH_IMAGE018
is the electrical angular velocity of the beam of light,
Figure 827764DEST_PATH_IMAGE019
、
Figure 728724DEST_PATH_IMAGE020
is the inductive component of the d-q axis,
Figure 416058DEST_PATH_IMAGE021
is a permanent magnet flux linkage, and is characterized in that,
Figure 514244DEST_PATH_IMAGE022
is time;
the external circuit model is controlled as follows;
Figure 654239DEST_PATH_IMAGE023
；（4）
Figure 93310DEST_PATH_IMAGE024
；（5）
in the formula (I), the compound is shown in the specification,
Figure 697467DEST_PATH_IMAGE025
is an electromotive force of the x-phase,
Figure 196581DEST_PATH_IMAGE026
is the voltage drop of N with respect to point of reference O,
Figure 823872DEST_PATH_IMAGE027
is the current of the x-phase,
Figure 568099DEST_PATH_IMAGE028
in order to be the load current,
Figure 230025DEST_PATH_IMAGE029
in order to be an equivalent resistance, the resistance,
Figure 900041DEST_PATH_IMAGE030
in order to be an equivalent inductance,
Figure 749048DEST_PATH_IMAGE031
in order to be an equivalent capacitance,
Figure 795501DEST_PATH_IMAGE032
is a direct-current voltage, and the voltage is,
Figure 311933DEST_PATH_IMAGE033
are respectively paired
Figure 152850DEST_PATH_IMAGE034
The three phases are calculated and,
Figure 489154DEST_PATH_IMAGE035
is time;
synthesizing the control models (2) - (5) to obtain a control model of the transmission chain system;
s2, after receiving a driving signal provided by a control algorithm in a control model, a driving circuit power tube generates a periodic PWM control signal, an electromagnetic finite element model of the generator is established on the basis of considering the PWM control signal, namely a control equation generated by the control model in the simultaneous S1 is established, and the electromagnetic parameters of the motor are solved by establishing a magnetic vector equation similar to a steady field;
wherein, the electromagnetic finite element model in the cell is as follows:
Figure 572254DEST_PATH_IMAGE001
；（1）
in the formula
Figure 943192DEST_PATH_IMAGE002
Is in the x-axis direction;
Figure 955011DEST_PATH_IMAGE003
is in the y-axis direction;
Figure 247452DEST_PATH_IMAGE004
is an initial value;
Figure 166866DEST_PATH_IMAGE005
is the reluctance ratio;
Figure 657890DEST_PATH_IMAGE006
is a magnetic vector, wherein the magnetic vector has only a component in the direction of the z-axis;
Figure 840610DEST_PATH_IMAGE007
is the source current density;
Figure 151506DEST_PATH_IMAGE008
is the tangential component of the magnetic field strength;
Figure 579338DEST_PATH_IMAGE009
is a first type boundary;
Figure 924869DEST_PATH_IMAGE010
is a second type boundary;
wherein:
Figure 481752DEST_PATH_IMAGE036
；（6）
in the formula (I), the compound is shown in the specification,
Figure 76682DEST_PATH_IMAGE037
、
Figure 806740DEST_PATH_IMAGE038
is the d-q component of the stator current;
the magnetic vector position distribution in the electromagnetic field of the drive chain system is influenced by a control algorithm and a drive external circuit in the control model, and the electromagnetic field performance under the control of the PWM waveform is considered in the formula (6);
s3, under the condition of PWM control waveform, calculating the electromagnetic parameters of the drive chain system, establishing a temperature field model of the drive chain system, taking the electromagnetic loss as a heat source in the temperature field model of the drive chain system, calculating the temperature distribution result, updating the electromagnetic parameters of the drive chain system, and performing iterative calculation until the temperature distribution error is less than 5%;
wherein the temperature field model within the cell is:
Figure 6777DEST_PATH_IMAGE039
； (7)
in the formula (I), the compound is shown in the specification,
Figure 531300DEST_PATH_IMAGE040
is the temperature of the boundary of the motor,
Figure 816788DEST_PATH_IMAGE041
is the temperature of the fluid near the motor boundary,
Figure 679150DEST_PATH_IMAGE042
and
Figure 733694DEST_PATH_IMAGE043
is an electric motor
Figure 429118DEST_PATH_IMAGE044
And
Figure 405164DEST_PATH_IMAGE045
the thermal conductivity of the material in the direction,
Figure 8183DEST_PATH_IMAGE046
is the sum of the heat flux densities in the motor,
Figure 713971DEST_PATH_IMAGE047
in order to obtain a heat transfer coefficient,
Figure 845875DEST_PATH_IMAGE048
to pass through the boundary x1The density of the heat flow is such that,
Figure 341841DEST_PATH_IMAGE049
in order to obtain the magnetic induction intensity,
Figure 420655DEST_PATH_IMAGE050
in the form of a time, the time,
Figure 184212DEST_PATH_IMAGE051
is in the direction of the x-axis,
Figure 487018DEST_PATH_IMAGE052
is in the y-axis direction;
s4, taking the stable temperature distribution result as the initial condition of the stress field differential equation model of the transmission chain system, and calculating the thermal stress distribution of the transmission chain system;
wherein the differential equation model of the stress field in the cell is:
Figure 968814DEST_PATH_IMAGE053
；（8）
in the formula (I), the compound is shown in the specification,
Figure 382478DEST_PATH_IMAGE054
、
Figure 266121DEST_PATH_IMAGE055
is composed of
Figure 739827DEST_PATH_IMAGE056
、
Figure 708920DEST_PATH_IMAGE057
The direction of the positive stress is the direction of the positive stress,
Figure 955969DEST_PATH_IMAGE058
in order to provide a shear force for each plane,
Figure 490855DEST_PATH_IMAGE059
、
Figure 135463DEST_PATH_IMAGE060
in order to be a surface force,
Figure 591852DEST_PATH_IMAGE061
in the direction of the x-axis,
Figure 347319DEST_PATH_IMAGE062
is in the y-axis direction;
s5: and outputting the structural model results of the calculated electromagnetic loss, temperature distribution diagram and stress distribution diagram, then comparing the results with the design indexes of the transmission chain system, and finally checking the reliability of the transmission chain.
Therefore, the influence of the transmission chain converter side on the whole electromagnetic, temperature and stress fields of the transmission chain system is considered, the electromagnetic loss of the generator under the PWM waveform is calculated through establishing a transmission chain system control algorithm and a driving external circuit mathematical model, the transmission chain system is subjected to electromagnetic-temperature field iterative calculation on the basis of the electromagnetic loss, and after the precision requirement is met, the obtained transmission chain system temperature distribution result is used as the initial condition of the stress field to calculate the thermal stress distribution;
in addition, the influence of the converter side on the transmission chain system is fully considered, the traditional chain system (converter-generator) is subjected to collaborative modeling and multi-physical field analysis, and the modeling process and the result of the method are more consistent with the actual structure and the operation condition of the transmission chain system.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the technical solutions of the present invention have been described in detail with reference to the foregoing embodiments, it should be understood by those skilled in the art that the technical solutions described in the foregoing embodiments can be modified or some technical features can be replaced equally; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims (6)

1. A method for collaborative modeling and multi-physics field analysis of a transmission chain system is characterized by comprising the following steps: comprises the following steps;
s1, establishing a drive chain system control model, wherein the drive chain system control model comprises a control algorithm model and a control external circuit model;
s2, after receiving a driving signal provided by a control algorithm in a control model, a driving circuit power tube generates a periodic PWM control signal, an electromagnetic finite element model of the generator is established on the basis of considering the PWM control signal, namely a control equation generated by the control model in the simultaneous S1 is established, and the electromagnetic parameters of the motor are solved by establishing a magnetic vector equation similar to a steady field;
wherein, the electromagnetic finite element model in the cell is as follows:
Figure 215413DEST_PATH_IMAGE001
；（1）
in the formula
Figure 312682DEST_PATH_IMAGE002
Is in the x-axis direction;
Figure 942246DEST_PATH_IMAGE003
is in the y-axis direction;
Figure 833979DEST_PATH_IMAGE004
is an initial value;
Figure 486677DEST_PATH_IMAGE005
is the reluctance ratio;
Figure 387637DEST_PATH_IMAGE006
is a magnetic vector position whichThe middle magnetic vector position only has a component in the direction of the z axis;
Figure 809391DEST_PATH_IMAGE007
is the source current density;
Figure 639069DEST_PATH_IMAGE008
is the tangential component of the magnetic field strength;
Figure 779063DEST_PATH_IMAGE009
is a first type boundary;
Figure 749293DEST_PATH_IMAGE010
is a second type boundary;
s3, calculating the electromagnetic parameters of the drive chain system under the condition of PWM control waveform, establishing a temperature field model of the drive chain system, calculating a temperature distribution result by taking electromagnetic loss as a heat source in the temperature field model of the drive chain system, updating the electromagnetic parameters of the drive chain system, and performing iterative calculation until the temperature distribution error is less than 5%;
s4, taking the stable temperature distribution result as the initial condition of the stress field differential equation model of the transmission chain system, and calculating the thermal stress distribution of the transmission chain system;
s5: and outputting the structural model results of the calculated electromagnetic loss, temperature distribution diagram and stress distribution diagram, then comparing the results with the design indexes of the transmission chain system, and finally checking the reliability of the transmission chain.
2. The method of claim 1, wherein the method comprises: the control algorithm model is as follows;
Figure 25554DEST_PATH_IMAGE011
；（2）
Figure 790248DEST_PATH_IMAGE012
；（3）
in the formula (I), the compound is shown in the specification,
Figure 151959DEST_PATH_IMAGE013
、
Figure 660300DEST_PATH_IMAGE014
being the d-q axis component of the stator voltage,
Figure 56647DEST_PATH_IMAGE015
、
Figure 225198DEST_PATH_IMAGE016
is the d-q component of the stator current,
Figure 870943DEST_PATH_IMAGE017
as the resistance of the stator,
Figure 651817DEST_PATH_IMAGE018
is the electrical angular velocity of the beam of light,
Figure 168249DEST_PATH_IMAGE019
、
Figure 9166DEST_PATH_IMAGE020
is the inductive component of the d-q axis,
Figure 611049DEST_PATH_IMAGE021
is a permanent magnet flux linkage, and is provided with a permanent magnet,
Figure 195614DEST_PATH_IMAGE022
is time.
3. The method of claim 2, wherein the method comprises: the external control circuit model is as follows;
Figure 68017DEST_PATH_IMAGE023
；（4）
Figure 79835DEST_PATH_IMAGE024
；（5）
in the formula (I), the compound is shown in the specification,
Figure 903435DEST_PATH_IMAGE025
is an electromotive force of the x-phase,
Figure 291691DEST_PATH_IMAGE026
is the voltage drop of N with respect to point O of reference,
Figure 48294DEST_PATH_IMAGE027
is the current of the x-phase,
Figure 231014DEST_PATH_IMAGE028
in order to be the load current,
Figure 807489DEST_PATH_IMAGE029
in order to be an equivalent resistance, the resistance,
Figure 238251DEST_PATH_IMAGE030
in order to be an equivalent inductance,
Figure 318202DEST_PATH_IMAGE031
in order to be an equivalent capacitance,
Figure 937403DEST_PATH_IMAGE032
is a direct-current voltage, and is,
Figure 470015DEST_PATH_IMAGE033
are respectively paired
Figure 200074DEST_PATH_IMAGE034
The three phases are calculated and,
Figure 134532DEST_PATH_IMAGE035
is time;
and (5) synthesizing the control models (2) to (5) to obtain the control model of the transmission chain system.
4. The method of claim 3, wherein the method comprises: in said S2;
Figure 924633DEST_PATH_IMAGE036
；（6）
in the formula (I), the compound is shown in the specification,
Figure 210121DEST_PATH_IMAGE037
、
Figure 979756DEST_PATH_IMAGE038
is the d-q component of the stator current;
the magnetic vector position distribution in the electromagnetic field of the drive chain system is influenced by a control algorithm and a driving external circuit in a control model, and the electromagnetic field performance under the control of the PWM waveform is considered in the formula (6).
5. The method of claim 1, wherein the method comprises: in S3, the temperature field model in the cell is:
Figure 34300DEST_PATH_IMAGE039
； (7)
in the formula (I), the compound is shown in the specification,
Figure 995303DEST_PATH_IMAGE040
is the temperature of the boundary of the motor,
Figure 502507DEST_PATH_IMAGE041
is the temperature of the fluid near the motor boundary,
Figure 574368DEST_PATH_IMAGE042
and
Figure 748998DEST_PATH_IMAGE043
is an electric motor
Figure 615323DEST_PATH_IMAGE044
And
Figure 875403DEST_PATH_IMAGE045
the thermal conductivity of the material in the direction,
Figure 983911DEST_PATH_IMAGE046
is the sum of the heat flux densities in the motor,
Figure 13047DEST_PATH_IMAGE047
in order to obtain a heat transfer coefficient,
Figure 847010DEST_PATH_IMAGE048
to pass through the boundary x1The density of the heat flow is such that,
Figure 594387DEST_PATH_IMAGE049
in order to obtain the magnetic induction intensity,
Figure 8050DEST_PATH_IMAGE050
as a matter of time, the time is,
Figure 891693DEST_PATH_IMAGE051
in the direction of the x-axis,
Figure 365399DEST_PATH_IMAGE052
in the y-axis direction.
6. The method of claim 1, wherein the method comprises: in S4, the differential equation model of the stress field in the cell is:
Figure 835957DEST_PATH_IMAGE053
；（8）
in the formula (I), the compound is shown in the specification,
Figure 787733DEST_PATH_IMAGE054
、
Figure 791461DEST_PATH_IMAGE055
is composed of
Figure 170490DEST_PATH_IMAGE056
、
Figure 626879DEST_PATH_IMAGE057
The direction of the positive stress is the direction of the positive stress,
Figure 382345DEST_PATH_IMAGE058
in order to be able to apply shear forces in the respective planes,
Figure 240580DEST_PATH_IMAGE059
、
Figure 56089DEST_PATH_IMAGE060
in order to be a surface force,
Figure 226871DEST_PATH_IMAGE061
in the direction of the x-axis,
Figure 520449DEST_PATH_IMAGE062
in the y-axis direction.
CN202210644031.5A 2022-06-09 2022-06-09 Transmission chain system collaborative modeling and multi-physical field analysis method Active CN114721276B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210644031.5A CN114721276B (en) 2022-06-09 2022-06-09 Transmission chain system collaborative modeling and multi-physical field analysis method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210644031.5A CN114721276B (en) 2022-06-09 2022-06-09 Transmission chain system collaborative modeling and multi-physical field analysis method

Publications (2)

Publication Number Publication Date
CN114721276A true CN114721276A (en) 2022-07-08
CN114721276B CN114721276B (en) 2022-08-09

Family

ID=82232670

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210644031.5A Active CN114721276B (en) 2022-06-09 2022-06-09 Transmission chain system collaborative modeling and multi-physical field analysis method

Country Status (1)

Country Link
CN (1) CN114721276B (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103678835A (en) * 2014-01-15 2014-03-26 三峡大学 Modeling method of motor in electromagnetic field-flow field-temperature field coupling calculation
CN103745070A (en) * 2014-01-28 2014-04-23 中国科学院电工研究所 Modeling and simulating method for mechanical transient characteristics of transmission chain of wind generating set
US20140175796A1 (en) * 2012-12-21 2014-06-26 Envision Energy (Denmark) Aps Wind turbine having a hts generator with a plurality of phases
CN106777459A (en) * 2016-11-10 2017-05-31 北京交通大学 The computational methods in generator shaft radial rotor temperature field
CN109063337A (en) * 2018-08-03 2018-12-21 内蒙古工业大学 A kind of more coupling magnetic field with electric circuit simulation methods of permanent-magnetic wind driven generator
US20200394345A1 (en) * 2017-09-11 2020-12-17 Romax Technology Limited Driveline Designer
CN112800685A (en) * 2021-03-26 2021-05-14 湖南大学 High-precision motor multi-physical field coupling simulation calculation method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140175796A1 (en) * 2012-12-21 2014-06-26 Envision Energy (Denmark) Aps Wind turbine having a hts generator with a plurality of phases
CN103678835A (en) * 2014-01-15 2014-03-26 三峡大学 Modeling method of motor in electromagnetic field-flow field-temperature field coupling calculation
CN103745070A (en) * 2014-01-28 2014-04-23 中国科学院电工研究所 Modeling and simulating method for mechanical transient characteristics of transmission chain of wind generating set
CN106777459A (en) * 2016-11-10 2017-05-31 北京交通大学 The computational methods in generator shaft radial rotor temperature field
US20200394345A1 (en) * 2017-09-11 2020-12-17 Romax Technology Limited Driveline Designer
CN109063337A (en) * 2018-08-03 2018-12-21 内蒙古工业大学 A kind of more coupling magnetic field with electric circuit simulation methods of permanent-magnetic wind driven generator
CN112800685A (en) * 2021-03-26 2021-05-14 湖南大学 High-precision motor multi-physical field coupling simulation calculation method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
DARIA KEPSU 等: "Interdisciplinary Design of a High-Speed Drivetrain for a Kinetic Compressor in a High-Temperature Heat Pump", 《IEEE》 *
孟祥玉 等: "电动汽车用感应电机电磁热耦合仿真研究", 《农业装备与车辆工程》 *
张松: "航空高力能永磁电机多物理场耦合优化设计与分析", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 *

Also Published As

Publication number Publication date
CN114721276B (en) 2022-08-09

Similar Documents

Publication Publication Date Title
CN1157845C (en) Non-synchronous motor parameter identification method
CN110768299B (en) Sequence impedance modeling and stability analysis method of load virtual synchronous machine
CN106357184A (en) Temperature compensation method for output torque of permanent magnet synchronous motor for vehicle based on neural network
CN114172412B (en) Parameter-free model prediction current control method for double three-phase permanent magnet motor
Hu et al. Loss minimization control strategy for linear induction machine in urban transit considering normal force
WO2022134661A1 (en) Method for selecting magnetization state of adjustable-flux permanent magnet synchronous motor in case of optimal control of full-speed domain efficiency and online control method
CN113285481B (en) Grid-connected converter inductance parameter online estimation method, prediction control method and system
Nasirian et al. High-fidelity magnetic characterization and analytical model development for switched reluctance machines
Yan et al. Torque estimation and control of PMSM based on deep learning
Sun et al. Efficient feedback linearization control for an IPMSM of EVs based on improved firefly algorithm
Wang et al. Low-switching-loss finite control set model predictive current control for IMs considering rotor-related inductance mismatch
CN109066784A (en) A kind of micro-capacitance sensor stability control method based on bifurcation theory
CN109802385B (en) Impedance modeling method of voltage source inverter
Yi et al. Loss minimization DTC electric motor drive system based on adaptive ANN strategy
Borisevich Numerical method for power losses minimization of vector-controlled induction motor
CN114721276B (en) Transmission chain system collaborative modeling and multi-physical field analysis method
CN111707938B (en) Motor simulator based on finite element reverse lookup table model
CN115642845B (en) Electromechanical actuating system multi-software joint simulation method based on model predictive control
WO2022257258A1 (en) Predicted current increment control method suitable for high-speed zone operation of permanent magnet synchronous motor
Zhu The key technologies for powertrain system of intelligent vehicles based on switched reluctance motors
CN114244216A (en) Permanent magnet synchronous motor parameter identification method, device and system
Jiang et al. Multi-objective Optimization Considering PET's Vibration Suppression of Dual Active Bridge Converter Based on BP-NSGA-II
Du et al. Research on torque observer of permanent magnet synchronous motor based on model integration
CN115051611B (en) Open winding motor simulator based on power electronic converter and control system thereof
Bizhani et al. A Multi-Input Multi-Output Model Predictive Direct Torque Control for Dual Mechanical Port Machine

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