CN110837706A - Parameterization calculation method for GVPI stator core mode - Google Patents
Parameterization calculation method for GVPI stator core mode Download PDFInfo
- Publication number
- CN110837706A CN110837706A CN201911050835.7A CN201911050835A CN110837706A CN 110837706 A CN110837706 A CN 110837706A CN 201911050835 A CN201911050835 A CN 201911050835A CN 110837706 A CN110837706 A CN 110837706A
- Authority
- CN
- China
- Prior art keywords
- core
- stator
- stator core
- elastic modulus
- gvpi
- 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
Links
Images
Classifications
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Iron Core Of Rotating Electric Machines (AREA)
Abstract
The invention discloses a parametric calculation method for a GVPI stator core mode. When the stator core mode adopting the GVPI technology is calculated, the improvement of the overall rigidity of the core by the stator winding behind the GVPI is simulated by establishing the stator winding with actual mass and low elastic modulus, so that the accuracy of the mode calculation of the stator core adopting the GVPI technology is improved. The method has the advantages that the calculation of multiple GVPI iron core modes can be realized through secondary development, the operation is simple, the standardization of the calculation method and the flow is realized, the human errors caused by different analysts are avoided, the calculation time is obviously shortened, and the efficiency is improved.
Description
Technical Field
The invention relates to a parameterization calculation method for a GVPI stator core mode, which is suitable for a two-pole or four-pole generator adopting a GVPI technology and belongs to the technical field of motors.
Background
When the rotor of the generator is connected with exciting current, the rotor becomes an electromagnet, and most of a rotating magnetic field generated by the electromagnet forms a closed loop through the rotor, an air gap and a stator core. The stator core is subjected to an electromagnetic excitation force of a certain frequency, and if the natural frequency of the core is close to the electromagnetic excitation frequency, destructive resonance is caused, so that this must be avoided in the structural design.
In the whole vacuum pressure impregnation (GVPI) process of the stator core and the winding, the stator core and the winding are crosslinked into a whole by impregnating resin in a vacuum and pressurized state, and the whole rigidity of the stator core is higher than that of a core which does not adopt the GVPI technology. In the conventional iron core mode calculation, only the rigidity contribution of an iron core yoke part is generally considered, the iron core is simply equivalent to a ring, and the corresponding mode is calculated. The method has the following defects: 1) the influence of iron core rigidity improvement brought by the GVPI process cannot be considered; 2) the stiffness contribution of the core teeth is slight. Therefore, the results obtained by the conventional core mode calculation method are lower than the actual values.
Disclosure of Invention
The purpose of the invention is: and the influence of the GVPI rear stator winding on the rigidity of the iron core is considered, so that a more accurate iron core mode calculation value is obtained.
In order to achieve the above object, the technical solution of the present invention is to provide a parametric calculation method for GVPI stator core mode, wherein a stator core is formed by stacking silicon steel sheets, and the method is characterized by comprising the following steps:
Establishing a simplified stator core two-dimensional model according to the design size of a generator stator, wherein the simplified stator core two-dimensional model consists of an iron core yoke part, an iron core tooth part and a stator winding;
the iron core yoke part and the tooth part are made of orthotropic materials;
the stator winding is an isotropic material having a low elastic modulus, and the stator winding is regarded as being formed by compounding copper wires with resin, and the elastic modulus of the stator winding is determined by the following formula:
wherein E is the elastic modulus of the stator winding, u is the volume ratio of the gelatinized coil bars in the stator winding, E1Is the elastic modulus of the gelled bar, E2Is the elastic modulus of the resin;
the density of the stator windings is determined by the following equation:
ρ＝M/(W*D*N*L)
in the formula, rho is the density of the stator winding, M is the total mass of parts in the stator core slot, W is the width of the core slot, D is the depth of the core slot, N is the number of the core slots, and L is the length of the core;
and extracting the Nth mode as a calculated value of the stator core mode according to N-2.5 POLE-4 and N-2.5 POLE-3, wherein POLE represents the number of generator stages.
Preferably, in step 2, the axial elastic modulus of the core yoke and the tooth part is far smaller than the elastic modulus in the radial direction and the tangential direction; the density of the iron core yoke part and the density of the tooth part are determined by multiplying the density of the silicon steel sheets by the iron core stacking coefficient.
Preferably, in step 2, the modulus of elasticity E of the gelled bar1Obtained through experiments or calculation; elastic modulus E of the resin2Obtained by experiments.
Preferably, in step 4, for the circular stator core geometric model, the first 7-order mode vibration modes are respectively ellipse, three-lobe, respiration, four-lobe and four-lobe; according to the characteristics of magnetic pulling force generated by a rotor, the two-pole synchronous generator extracts 1 st and 2 nd order modes of a stator core geometric model; the four-pole synchronous generator extracts the 6 th and 7 th order modes of the stator core geometric model.
The secondary development is carried out according to the calculation method and the steps, so that the parametric analysis of the GVPI stator core modal calculation can be realized. For the stator core adopting the GVPI technology, the analysis can be completed only by inputting parameters in a dialog box. The parameters to be input specifically include: the core length, the core outer diameter, the core inner diameter, the core slot depth, the core slot width, the core slot number, the rotor pole number, the stacking coefficient and the volume ratio of the gelled wire bars in the stator winding.
The implementation of the invention has the following beneficial effects: when the stator core mode adopting the GVPI technology is calculated, the improvement of the overall rigidity of the core by the stator winding behind the GVPI is simulated by establishing the stator winding with actual mass and low elastic modulus, so that the accuracy of the mode calculation of the stator core adopting the GVPI technology is improved.
The method has the advantages that the calculation of multiple GVPI iron core modes can be realized through secondary development, the operation is simple, the standardization of the calculation method and the flow is realized, the human errors caused by different analysts are avoided, the calculation time is obviously shortened, and the efficiency is improved.
Drawings
FIG. 1 is a schematic diagram of a three-dimensional model of a stator core of a generator using GVPI technology;
fig. 2 is a two-dimensional simplified model of a stator core, in which 1 denotes a core yoke, 2 denotes a core tooth, and 3 denotes a stator winding;
fig. 3 is a deformation diagram of a stator core under the action of magnetic tension, wherein 4 represents elliptical frequency doubling vibration, 5 represents magnetic tension, 6 represents air gap flux density, 7 represents a rotor, and 8 represents the core;
fig. 4(a) and 4(b) are diagrams of the modal shape of the stator core obtained by finite element method calculation;
fig. 5 is a diagram of a modal shape of a stator core obtained by an experiment.
Detailed Description
The invention will be further illustrated with reference to the following specific examples. It should be understood that these examples are for illustrative purposes only and are not intended to limit the scope of the present invention. Further, it should be understood that various changes or modifications of the present invention may be made by those skilled in the art after reading the teaching of the present invention, and such equivalents may fall within the scope of the present invention as defined in the appended claims.
The invention provides a parameterization calculation method for a GVPI stator core mode, which comprises the following steps of:
According to the design size of the generator stator, a simplified stator core two-dimensional model is established, and the simplified model comprises three parts, namely a core yoke part, a core tooth part and a stator winding, as shown in fig. 2. The simplified stator core two-dimensional model is established, and 5 size parameters of the core outer diameter, the core inner diameter, the core groove depth, the core groove width and the core groove number are needed.
The yoke part and the tooth part of the iron core are made of orthogonal anisotropic materials, because the stator iron core is formed by laminating silicon steel sheets, and the axial elastic modulus of the stator iron core is far smaller than the radial elastic modulus and the tangential elastic modulus; the density is determined by multiplying the density of the silicon steel sheets by the iron core stacking coefficient.
The stator winding is defined as isotropic material with lower elastic modulus because after GVPI, resin material is filled in the insulation layer of the stator winding and the gaps between the winding and the slots, the stator winding is crosslinked with the iron core into a whole, the whole stator winding can be regarded as being formed by compounding copper wires and resin, and the elastic modulus can be determined by the following formula:
wherein E is the elastic modulus of the stator winding; u is the volume ratio of the gelled wire rod in the stator winding, and in the embodiment, the volume of the copper wire is 52%; e1The elastic modulus of the gelled wire rod can be obtained through tests or calculation, and is 76000MPa in the embodiment; e2The elastic modulus of the resin was obtained by a tensile test, and the elastic modulus of the resin used in this example was 1230 MPa. In this way,in this embodiment, the calculated elastic modulus of the stator winding is 2518 Mpa.
The density of the stator winding is determined according to the principle of mass equivalence, and can be specifically determined by the following formula:
ρ＝M/(W*D*N*L)
where ρ is the density of the stator winding, M is the total mass of the stator core in-slot components, W is the core slot width, D is the core slot depth, N is the number of core slots, and L is the core length.
And 3, calculating the iron core mode by adopting a finite element method. In this embodiment, the established calculation model is modal solved by using ANSYS.
And 4, extracting the Nth-order mode as a calculated value of the stator core mode according to the formulas N-2.5 POLE-4 and N-2.5 POLE-3, wherein POLE represents the number of generator steps. For the ring model, the first 7 orders of modal vibration modes are respectively ellipse, three-valve, respiration, four-valve and four-valve; according to the characteristics of magnetic pulling force generated by a rotor, the two-pole synchronous generator should pay attention to the natural frequency of the elliptical vibration mode, so that the 1 st and 2 nd order modes need to be extracted; the quadrupole synchronous generator should pay attention to the natural frequency of four-lobe vibration mode, so that the 6 th and 7 th order modes need to be extracted
The calculation object of the present embodiment is a two-pole synchronous generator. The magnetic pulling force generated by the rotor reaches the maximum value at the center of the magnetic poles and is the minimum value between the magnetic poles as shown in fig. 3, so that the stator core is forced to present 4-node elliptical vibration, and the 1 st and 2 nd order elliptical mode shapes of the stator core should be concerned to avoid the resonance of the core.
The invention is further illustrated by the following specific data:
the following parameters are input into a parameter input dialog of the GVPI core mode calculation program:
the length of the iron core is 3400 mm;
the outer diameter of the iron core is 2600 mm;
the inner diameter of the iron core is 1170 mm;
the depth of the iron core groove is 245.4 mm;
the width of the iron core slot is 25.6 mm;
the number of the iron core slots is 66;
the number of rotor poles is 2;
the stacking coefficient is 0.95;
the volume ratio of the gelled bars in the stator winding is 0.52.
After the calculation is completed, the core mode results output by the calculation program are shown in fig. 4(a) and 4(b), and it can be seen from the figure that the mode frequency of the stator core is 201.059Hz/203.398Hz, and the mode shape is an elliptical shape; the actual elliptical modal frequency of the stator core in this embodiment is measured by experiments to be 199.4Hz, and as shown in fig. 5, the stator core modal method provided by the present invention has a relatively high calculation accuracy.
Claims (4)
1. A parameterization calculation method for a GVPI stator core mode is characterized by comprising the following steps of:
step 1, establishing a simplified stator core geometric model
Establishing a simplified stator core two-dimensional model according to the design size of a generator stator, wherein the simplified stator core two-dimensional model consists of an iron core yoke part, an iron core tooth part and a stator winding;
step 2, defining material parameters of a stator core geometric model, wherein:
the iron core yoke part and the tooth part are made of orthotropic materials;
the stator winding is an isotropic material having a low elastic modulus, and the stator winding is regarded as being formed by compounding copper wires with resin, and the elastic modulus of the stator winding is determined by the following formula:
wherein E is the elastic modulus of the stator winding, u is the volume ratio of the gelatinized coil bars in the stator winding, E1Is the elastic modulus of the gelled bar, E2Is the elastic modulus of the resin;
the density of the stator windings is determined by the following equation:
ρ＝M/(W*D*N*L)
in the formula, rho is the density of the stator winding, M is the total mass of parts in the stator core slot, W is the width of the core slot, D is the depth of the core slot, N is the number of the core slots, and L is the length of the core;
step 3, calculating the mode of the stator core geometric model by adopting a finite element method;
step 4, extracting the modal shape:
and extracting the Nth mode as a calculated value of the stator core mode according to N-2.5 POLE-4 and N-2.5 POLE-3, wherein POLE represents the number of generator stages.
2. The parametric calculation method for the GVPI stator core mode shape according to claim 1, wherein in step 2, the axial elastic modulus of the core yoke and the tooth part is far smaller than the radial elastic modulus and the tangential elastic modulus; the density of the iron core yoke part and the density of the tooth part are determined by multiplying the density of the silicon steel sheets by the iron core stacking coefficient.
3. The parametric calculation method for the mode of the GVPI stator core as claimed in claim 1, wherein in step 2, the elastic modulus E of the colloidizing bar1Obtained through experiments or calculation; elastic modulus E of the resin2Obtained by experiments.
4. The parametric calculation method for the mode of the GVPI stator core according to claim 1, wherein in step 4, for the circular stator core geometric model, the first 7-order mode vibration modes are respectively ellipse, three-lobe, breathing, four-lobe and four-lobe; according to the characteristics of magnetic pulling force generated by a rotor, the two-pole synchronous generator extracts 1 st and 2 nd order modes of a stator core geometric model; the four-pole synchronous generator extracts the 6 th and 7 th order modes of the stator core geometric model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911050835.7A CN110837706B (en) | 2019-10-31 | 2019-10-31 | Parameterized calculation method for GVPI stator core mode |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911050835.7A CN110837706B (en) | 2019-10-31 | 2019-10-31 | Parameterized calculation method for GVPI stator core mode |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110837706A true CN110837706A (en) | 2020-02-25 |
| CN110837706B CN110837706B (en) | 2023-06-02 |
Family
ID=69575913
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911050835.7A Active CN110837706B (en) | 2019-10-31 | 2019-10-31 | Parameterized calculation method for GVPI stator core mode |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110837706B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113310520A (en) * | 2021-05-28 | 2021-08-27 | 包头市威丰稀土电磁材料股份有限公司 | Method for evaluating quality of ultrathin silicon steel core for saturated reactor of direct-current converter valve |
| CN114547946A (en) * | 2022-03-03 | 2022-05-27 | 哈尔滨理工大学 | Method for calculating orthogonal anisotropic material parameters of motor |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105205285A (en) * | 2015-10-20 | 2015-12-30 | 河海大学 | BSRM stator mode finite element modeling method considering influence of windings |
| CN108181944A (en) * | 2018-01-18 | 2018-06-19 | 广东电网有限责任公司电力科学研究院 | A kind of generator windings end portion vibration control method and system |
-
2019
- 2019-10-31 CN CN201911050835.7A patent/CN110837706B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105205285A (en) * | 2015-10-20 | 2015-12-30 | 河海大学 | BSRM stator mode finite element modeling method considering influence of windings |
| CN108181944A (en) * | 2018-01-18 | 2018-06-19 | 广东电网有限责任公司电力科学研究院 | A kind of generator windings end portion vibration control method and system |
Non-Patent Citations (4)
| Title |
|---|
| V.KRISHNAN: "《Analysis and experimental validation of performance for a switched reluctance motor with pre-fabricated soft magnetic core》" * |
| 张艺华;闫兵;张胜杰;唐琴;吴磊平;: "基于ANSYS的车用交流发电机定子建模方法分析" * |
| 张金华;夏黎;彭涵;陈永庚;: "GVPI定子线棒稳态热应力分析" * |
| 贺岩松;张浩;赵勤;孔祥杰;: "基于物理模型爪极发电机定子系统模态分析" * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113310520A (en) * | 2021-05-28 | 2021-08-27 | 包头市威丰稀土电磁材料股份有限公司 | Method for evaluating quality of ultrathin silicon steel core for saturated reactor of direct-current converter valve |
| CN114547946A (en) * | 2022-03-03 | 2022-05-27 | 哈尔滨理工大学 | Method for calculating orthogonal anisotropic material parameters of motor |
| CN114547946B (en) * | 2022-03-03 | 2022-07-22 | 哈尔滨理工大学 | Method for calculating orthogonal anisotropic material parameters of motor |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110837706B (en) | 2023-06-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Li et al. | Comparison of 12/8 and 6/4 switched reluctance motor: Noise and vibration aspects | |
| JP5646119B1 (en) | Permanent magnet type rotating electric machine | |
| Liu et al. | Nonlinear equivalent magnetic network of a linear permanent magnet vernier machine with end effect consideration | |
| CN110837706B (en) | Parameterized calculation method for GVPI stator core mode | |
| Verez et al. | Influence of slots and rotor poles combinations on noise and vibrations of magnetic origins inu'-core flux-switching permanent magnet machines | |
| CN113489387B (en) | Method for weakening electromagnetic vibration and noise of permanent magnet synchronous motor with specific frequency | |
| Frias et al. | Rotor and stator shape optimization of a synchronous machine to reduce iron losses and acoustic noise | |
| Kim et al. | Core loss analysis of permanent magnet synchronous generator with slotless stator | |
| Kahourzade et al. | Investigation of emerging magnetic materials for application in axial-flux PM machines | |
| Boroujeni et al. | Design, prototyping, and analysis of a novel tubular permanent-magnet linear machine | |
| CN113765324B (en) | Asymmetric modular permanent magnet auxiliary synchronous reluctance motor and vibration noise optimization method thereof | |
| CN105205285A (en) | BSRM stator mode finite element modeling method considering influence of windings | |
| Hasan et al. | Effect of pole shaping on cogging torque, torque ripple and vibrational performance in consequent pole TFM | |
| Verez et al. | Impact of pole and slot combination on noise and vibrations of flux-switching PM machines | |
| Estenlund et al. | Dovetail Design for Direct Cooled Rotor: Design and Manufacturing | |
| Ferkova | Comparison between 2D and 3D modelling of induction machine using finite element method | |
| Cao et al. | Magnetic field analytical solution and electromagnetic force calculation of coreless stator axial-flux permanent-magnet machines | |
| CN106383971B (en) | Improved motor stator core vibration analysis model caused by magnetostriction | |
| Fu et al. | A novel 2-D simplified model for investigating the rotor shape of homopolar inductor alternator | |
| Lim et al. | Embedded stator end-windings in soft magnetic composite and laminated surface PM machines | |
| Podhajecki et al. | Numerical and experimental vibration analysis of BLDC motor | |
| Zarate et al. | A study of the effects of the magnetization tolerance in the manufacturing process on permanent magnet synchronous machines | |
| Verez et al. | Natural frequencies analytical modeling of small industrial radial flux permanent magnet motors | |
| Ferraris et al. | Design optimization of bonded PM BLDC motors with reference to the cogging torque amplitude | |
| Isfanuti et al. | Outer ferrite-PM-rotor BLAC motor characterization: FEM-assisted optimal design and preliminary experiments |
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 |