CN103559365B - A kind of parametric modeling method suitable for Generator Stator Winding Ends coil - Google Patents
A kind of parametric modeling method suitable for Generator Stator Winding Ends coil Download PDFInfo
- Publication number
- CN103559365B CN103559365B CN201310574323.7A CN201310574323A CN103559365B CN 103559365 B CN103559365 B CN 103559365B CN 201310574323 A CN201310574323 A CN 201310574323A CN 103559365 B CN103559365 B CN 103559365B
- Authority
- CN
- China
- Prior art keywords
- point
- coil
- ordinate
- polar
- rectangular
- 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.)
- Active
Links
Landscapes
- Windings For Motors And Generators (AREA)
- Cosmetics (AREA)
Abstract
The present invention relates to a kind of parametric modeling method suitable for Generator Stator Winding Ends coil, it is characterised in that step is:Step 1, input coil geometric parameter;Step 2, the cross sectional shape for defining coil;Step 3, the modeling of stator winding end coil launcher;Step 4, set up coil FEM model.It is an advantage of the invention that:Make stator winding end coil modeling parameters, and can simultaneously set up the FEM model of end coil, to be applied to the research of stator winding end vibration characteristics.
Description
Technical field
The present invention relates to a kind of parametric modeling method of Generator Stator Winding Ends coil, it is adaptable to which end coil is
Each model unit of involute expanded form, belongs to the dynamic characteristic calculation and improved skill of Generator Stator Winding Ends structure
Art field.
Background technology
Stator winding end coil carries the important function that induced-current is exported power network in generator.Its vibration feelings
Condition is always the problem being primarily upon during generator is researched and developed.Wherein stator winding end coil is due to the expansion using involute
Mode, the FEM model for hence setting up reasonable end coil has certain difficulty.
The modeling method of conventional end coil is to set up model using Three-dimensional Design Software, and its operating process is relatively
Complexity, and modification is difficult, while being not easily applicable to the finite element analysis of stator winding end structure.
The content of the invention
The technical problem to be solved in the present invention is to provide a kind of parametrization for Generator Stator Winding Ends coil and builds
Mould method, makes stator winding end coil modeling parameters, and can simultaneously set up the FEM model of end coil, so as to
It is applied to the research of stator winding end vibration characteristics.
In order to solve the above-mentioned technical problem, the technical scheme is that providing a kind of suitable for generator unit stator winding
The parametric modeling method of end coil, it is characterised in that step is:
Step 1, input coil geometric parameter, at least including coil section overall width ECBW, coil section total height ECBH,
Coil section copper cash partial width ECBCW, coil section copper cash Partial Height ECBCH, coil section hollow copper wire hole width
ECBCHW, coil section hollow copper string holes height ECBCHH, coil are to distance between center line TECRH unshakable in one's determination, coil base radius
TECRJ, coil straight line to the excessive chamfer radius TECRO of the conical surface, coil involute beginning chamfer radius TECR1, coil gradually
Burst at the seams chamfer radius TECR2 at end, coil goes out groove straight line portion length TECHZ, coil involute and starts preceding length of straigh line
TECBH1, coil involute terminate across angle TECBYY, P point of rear length of straigh line TECBH2, coil conical degree of conical surface CONTP, coil
The involute angle of spread and polar angle difference DEFAX, C point involute angle of spread DEFAY, wherein, P points were O points
The intersection point for gradually opening the tangent line with involute that terminate arc section after involute expansion of the coil under polar coordinates, O points exist for coil
The basic circle center of circle after involute expansion under polar coordinates, C points are the involute section after involute of the coil under polar coordinates launches
Terminal;
Step 2, the cross sectional shape for defining coil;
Step 3, the modeling of stator winding end coil launcher, including:
Step 3.1 solves coil key point coordinates:
The geometrical curve of coil is formed by connecting by multiple key points, and the coordinate of each key point is stored in aray variable KPCS
In (i, j), wherein i is coordinate components, i=1, and 2,3 represent the X-axis of key point, Y-axis and Z axis coordinate respectively, and j is that key point is compiled
Number, variable KPCS (m, n) of defining arrays, wherein m be coordinate components, m=1,2 respectively represent key point polar polar diameter R and
Polar angle φ, n are numbered for key point, and the formula for calculating crucial point coordinates at each section of coil is as follows:
Coil beginning straightway key point coordinates:
Starting point coordinate:KPCS (1,1)=0,
KPCS (2,1)=TECRH,
KPCS (3,1)=O;
Terminal point coordinate:KPCS (1,2)=0,
KPCS (2,2)=TECRH,
KPCS (3,2)=TECHZ;
The crucial point coordinates of coil cone angle transition circle section:
Starting point coordinate:KPCS (1,3)=O,
KPCS (2,3)=TECRH+TECRO-TECROcos (CONTP), KPCS (3,3)=TECRZ+TECROsin
(CONP);
Terminal point coordinate:KPCS (Isosorbide-5-Nitrae)=0,
KPCS (2,4)=TECRH+TECRO-TECROcos (CONP)+TECBH1sin (CONTP),
KPCS (3,4)=TECRZ+TECROsin (CONP)+TECBH1cos (CONTP);
The coordinate of the initial circular arc waypoint of involute:
1) A point coordinates, A points are the starting point of the initial arc section of involute:
A point polar diametersA point polar angles THTA=45;
A points rectangular co-ordinate X values XA=RAcos (THTA), A point rectangular co-ordinate Y value YA=RAsin (THTA), vertex of a cone Z sit
Mark OZ=KPCS (3,4)-RAcos (CONTP);
2) coordinate of fillet center of circle O1 is sought according to A point coordinates, O1 is the center of circle of the initial arc section of involute:
O1 points rectangular co-ordinate X values XO1=XA+TECR1cos (THTA+90), O1 point rectangular co-ordinate Y value YO1=YA+
TECR1·sin(THTA+90);
O1 point polar diametersD1 point polar angles
3) E points are coordinates, and E points are the point of contact of O1 points and basic circle tangent line:
E points rectangular co-ordinate X values XE=TECRJcos (THTE), E point rectangular co-ordinate Y value
YE=TECRJsin (THTE), wherein, THTE=THTO1+THTEOO1,
4) B point coordinates, B points are the starting point of involute section:
B points rectangular co-ordinate X values XB=XO1+TECR1cos (THTO1E), B point rectangular co-ordinate Y value YB=YO1+TECR1
Sin (THTO1E), wherein,
B point polar diametersB point polar angles
The involute angle of spread of B points:
5) angle of spread of the coil in sector:
TECBEA=TECBYY·sin(TECTP);
After step 3.2, zero point of launching polar coordinates zero point and involute overlap, each point coordinates is solved again:
1) A point coordinates:
A point rectangular co-ordinate X values XA '=RAcos (THTA '), A point rectangular co-ordinatesYValue
YA '=RAsin (THTA '), wherein, THTA '=THTA+ANGB-THTB;
2) O1 point coordinates:
O1 point rectangular co-ordinate X values XO1 '=XA '+TECR1cos (THTA '+90), O1 point rectangular co-ordinate Y values YO1 '=YA '+
TECR1·sin(THTA′+90);
O1 point polar angles
3) E point coordinates:
E point rectangular co-ordinate X values XE '=TECRJcos (THTE '),
E point rectangular co-ordinate Y values YE '=TECRJsin (THTE '), wherein,
THTE′=THTO1′+THTEOO1:
4) B point coordinates:
B point rectangular co-ordinate X values XB '=XO1 '+TECR1cos (THTO1E '),
B point rectangular co-ordinate Y values YB '=YO1 '+TECR1sin (THTO1E '),
Wherein,
B point polar anglesB point polar diameters
5) P point coordinates:
P point polar angle ANGP=DEFAX+THTP, wherein, THTP=TECBEA+THTA ';
6) C point coordinates:
C point rectangular co-ordinate X values
XC '=TECRJcos (ANGC)+TECRJANGC π/180sin (ANGC),
C point rectangular co-ordinate Y values
YC '=TECRJsin (ANGC)-TECRJANGC π/180cos (ANGC),
Wherein, ANGC=DEFAY;
C point polar diametersC point polar angles
7) O2 point coordinates, O2 points terminate the center of circle of arc section for involute:
O2 point rectangular co-ordinate X values:
XO2 '=TECRJcos (ANGC)+[TECRJANGC π/180SIN (ANGC)] sin (ANGC);
O2 point rectangular co-ordinate Y values:
YO2 '=TECRJsin (ANGC)-[TECRJANGC π/180COS (ANGC)] cos (ANGC);
8) D point coordinates, D points terminate the terminal of arc section for involute:
D point rectangular co-ordinate X values XD '=XO2 '+TECR2cos (THTP+90),
D point rectangular co-ordinate Y values YD '=XD ' tan (THTP),
D point polar diameters
Step 3.3, continuation solve each section of coordinate of key point of coil, and the key point number of the initial arc section of involute is
N1, the number of involute section key point is N2, and the key point number that involute terminates arc section is N3, then have:
1) cycle calculations obtain the coordinate of the initial each key point in arc section inside of involute respectively, wherein:
The N1 polar diameter of key point is RB ', and polar angle is THTB '-THTA ', the rectangular co-ordinate X of remaining k-th key point
Value, K=1,2 ... ..., (N1-1):
XTEMP=XO1 '+TECR1cos (TEMPTHTA+TEMPTHTABK/N1),
The rectangular co-ordinate Y value of k-th key point:
YTEMP=YO1 '+TECR1sin (TEMPTHTA+TEMPTHTABK/N1);
The polar diameter of k-th key point
Polar angle atan (YTEMP, XTEMP)-THTA ' of k-th key point, wherein:
2) cycle calculations obtain the coordinate of each key point of involute intersegmental part respectively, wherein:
The N2 polar diameter of key point is RC ', and polar angle is THTC '-THTA ', the rectangular co-ordinate X of remaining k-th key point
Value, K=1,2 ... ..., (N2-1):
XTEMP=TECRJcos (ANGB+TEMPBCK/N2)+
TECRJ (ANGB+TEMPBCK/N2) π/180sin (ANGB+TEMPBCK/N2) '
The rectangular co-ordinate Y value of k-th key point:
YTEMP=TECRJsin (ANGB+TEMPBCK/N2)-
TECRJ (ANGB+TEMPBCK/N2) π/180cos (ANGB+TEMPBCK/N2) '
The polar diameter of k-th key point
Polar angle atan (YTEMP, XTEMP)-THTA ' of k-th key point, wherein:
TEMPBC=ANGC-ANGB:
3) cycle calculations obtain the coordinate that involute terminates each key point in arc section inside respectively, wherein:
The N3 polar diameter of key point is RD ', and polar angle is THTD '-THTA ', the rectangular co-ordinate X of remaining k-th key point
Value, K=1,2 ... ..., (N3-1):
XTEMP=XO2 '+TECR2cos (TEMPTHTC-TEMPTHTCDK/N3);
The rectangular co-ordinate Y value of k-th key point:
YTEMP=YO2 '+TECR2sin (TEMPTHTC-TEMPTHTCDK/N3);
The polar diameter of k-th key point
Polar angle atan (YTEMP, XTEMP)-THTA ' of k-th key point, wherein:
If negative value, then it is modified to
TEMPTHTCD=TEMPTHTC-TEMPTHTD:
If negative value, then it is modified to
Step 3.4, the polar coordinates of each key point of coil are converted into cartesian coordinate using cycle calculations, wherein, k-th
KPCS (Isosorbide-5-Nitrae+K)=TECITEMPRsin (TEMPTHT) of key point,
KPCS (2,4+K)=TEMPRcos (TEMPTHT),
KPCS (3,4+K)=ICS (1, K) cos (CONTP)+OZ, wherein:
TEMPR=ICS (1, K) sin (CONTP), ICS (1, K) are the polar diameter of k-th key point;TEMPTHT=ICS (2,
K)/sin (CONTP), ICS (2, K) are the polar angle of k-th key point;
Step 3.5, set up key point:According to each crucial point coordinates that draws is solved in step 3.4, each key point is set up;
Step 3.6, the SPL for setting up coil:The key point set up in Connection Step 3.5, that is, generate the batten of coil
Curve;
Step 4, set up coil FEM model.
It is an advantage of the invention that:Make stator winding end coil modeling parameters, and can simultaneously set up end coil
FEM model, to be applied to the research of stator winding end vibration characteristics.
Brief description of the drawings
Fig. 1 is the stator winding end coil involute expanded form figure in the present invention.
Specific embodiment
To become apparent the present invention, hereby with preferred embodiment, the present invention will be further described.
The invention provides a kind of parametric modeling method suitable for Generator Stator Winding Ends coil, its step
For:
Step 1, input coil geometric parameter, including coil section overall width ECBW, coil section total height ECBH, coil
Section copper cash partial width ECBCW, coil section copper cash Partial Height ECBCH, coil section hollow copper wire hole width ECBCHW,
Coil section hollow copper string holes height ECBCHH, coil are to distance between center line TECRH unshakable in one's determination, coil base radius TECRJ, coil
Straight line terminates place to the excessive chamfer radius TECR0 of the conical surface, coil involute beginning chamfer radius TECR1, coil involute
Chamfer radius TECR2, coil go out groove straight line portion length TECHZ, coil involute and start preceding length of straigh line TECBH1, coil
The involute that involute terminates across angle TECBYY, P point of rear length of straigh line TECBH2, coil conical degree of conical surface CONTP, coil launches
The involute angle of spread DEFAY of difference DEFAX, the C point of angle and polar angle, wherein, P points were the coil of O points in polar coordinates
Under involute launch after the intersection point for gradually opening the tangent line and the involute that terminate arc section, O points be coil under polar coordinates gradually
The basic circle center of circle burst at the seams after launching, C points are the terminal of the involute section after involute of the coil under polar coordinates launches;
Step 2, the cross sectional shape for defining coil, in the present embodiment, the rectangular cross-section of coil, mainly by two kinds of materials
Constitute, section inner side is copper product, and outside is insulating materials.Pair cross-section is needed when setting up section carries out subdivision, and different
Positions of materials specifies corresponding material parameter.Mesh generation is carried out to coil section, the coil section that will be defined is saved as
.sect file for coil finite element modeling below when use.Idiographic flow is as follows:
Step 2.1, the square-section for setting up coil;
Step 2.2, by the different by two parts of section subdivision copper cash and insulation of coil method;
Step 2.3, specified mesh generation size;
Step 2.4, specified material parameter;
Step 2.5, division section grid;
Step 2.6, preservation .sect custom cross section files.
Step 3, the modeling of stator winding end coil launcher, including:
Step 3.1 solves coil key point coordinates:
Involute expanded form of Fig. 1 coils under polar coordinates, wherein, the O points position basic circle center of circle, at the beginning of A points are involute
The starting point of beginning arc section, B points position terminal is also the starting point of involute section, and O1 is the center of circle of the initial arc section of involute, C points
It is the terminal of involute section, the starting point of arc section is also terminated for involute, D points position involute terminates the terminal of arc section, O2 points
Terminate the center of circle of arc section for involute, E points are the point of contact of O1 points and basic circle tangent line, P points were that O points gradually open end arc section
Tangent line and involute intersection point.
The geometrical curve of coil is formed by connecting by multiple key points, and the coordinate of each key point is stored in aray variable KPCS
In (i, j), wherein i is coordinate components, i=1, and 2,3 represent the X-axis of key point, Y-axis and Z axis coordinate respectively, and j is that key point is compiled
Number, variable KPCS (m, n) of defining arrays, wherein m be coordinate components, m=1,2 respectively represent key point polar polar diameter R and
Polar angle φ, n are numbered for key point, and the formula for calculating crucial point coordinates at each section of coil is as follows:
Coil beginning straightway key point coordinates:
Starting point coordinate:KPCS (1,1)=0,
KPCS (2,1)=TECRH,
KPCS (3,1)=0;
Terminal point coordinate:KPCS (1,2)=0,
KPCS (2,2)=TECRH,
KPCS (3,2)=TECHZ;
The crucial point coordinates of coil cone angle transition circle section:
Starting point coordinate:KPCS (1,3)=0,
KPCS (2,3)=TECRH+TECRO-TECROcos (CONTP),
KPCS (3,3)=TECRZ+TECROsin (CONP);
Terminal point coordinate:KPCS (Isosorbide-5-Nitrae)=0,
KPCS (2,4)=TECRH+TECRO-TECROcos (CONP)+TECBH1sin (CONTP),
KPCS (3,4)=TECRZ+TECROsin (CONP)+TECBH1cos (CONTP);
The coordinate of the initial circular arc waypoint of involute:
1) A point coordinates, A points are the starting point of the initial arc section of involute:
A point polar diametersA point polar angles THTA=45;
A points rectangular co-ordinate X values XA=RAcos (THTA), A point rectangular co-ordinate Y value YA=RAsin (THTA), vertex of a cone Z sit
Mark OZ=KPCS (3,4)-RAcos (CONTP);
2) coordinate of fillet center of circle O1 is sought according to A point coordinates, O1 is the center of circle of the initial arc section of involute:
O1 points rectangular co-ordinate X values XO1=XA+TECR1cos (THTA+90), O1 point rectangular co-ordinate Y value YO1=YA+
TECR1·sin(THTA+90);
O1 point polar diametersO1 point polar angles
3) E point coordinates, E points are the point of contact of O1 points and basic circle tangent line:
E points rectangular co-ordinate X values AE=TECRJcos (THTE), E point rectangular co-ordinate Y value YE=TECRJsin (THTE),
Wherein, THTE=THTO1+THTEOO1,
4) B point coordinates, B points are the starting point of involute section:
B points rectangular co-ordinate X values XB=XO1+TECR1cos (THTO1E), B point rectangular co-ordinate Y value YB=YO1+TECR1
Sin (THTO1E), wherein,
B point polar diametersB point polar angles
The involute angle of spread of B points:
5) angle of spread of the coil in sector:
TECBEA=TECBYY·sin(TECTP);
After step 3.2, zero point of launching polar coordinates zero point and involute overlap, each point coordinates is solved again:
1) A point coordinates:
A point rectangular co-ordinate X values XA '=RAcos (THTA '), A point rectangular co-ordinate Y values
YA ' RAsin (THTA '), wherein, THTA '=THTA+ANGB-THTB;
2) O1 point coordinates:
O1 point rectangular co-ordinate X values XO1 '=XA '+TECR1cos (THTA '+90), O1 point rectangular co-ordinate Y values YO1 '=YA '+
TECR1·sin(THTA′+90);
O1 point polar angles
3) E point coordinates:
E point rectangular co-ordinate X values XE '=TECRJcos (THTE '),
E point rectangular co-ordinate Y values YE '=TECRJsin (THTE '), wherein,
THTE′=THTO1′+THTEOO1:
4) B point coordinates:
B point rectangular co-ordinate X values XB '=XO1 '+TECR1cos (THTO1E '),
B point rectangular co-ordinate Y values YB '=YO1 '+TECR1sin (THTO1E '),
Wherein,
B point polar anglesB point polar diameters
5) P point coordinates:
P point polar angle ANGP=DEFAX+THTP, wherein, THTP=TECBEA+THTA ';
6) C point coordinates:
C point rectangular co-ordinate X values
XC '=TECRJcos (ANGC)+TECRJANGC π/180sin (ANGC),
C point rectangular co-ordinate Y values
YC '=TECRJsin (ANGC)-TECRJANGC π/180cos (ANGC),
Wherein, ANGC=DEFAY;
C point polar diametersC point polar angles
7) O2 point coordinates, O2 points terminate the center of circle of arc section for involute:
O2 point rectangular co-ordinate X values:
XO2 '=TECRJcos (ANGC)+[TECRJANGC π/180SIN (ANGC)] sin (ANGC);
O2 point rectangular co-ordinate Y values:
YO2 '=TECRJsin (ANGC)-[TECRJANGC π/180COS (ANGC)] cos (ANGC);
8) D point coordinates, D points terminate the terminal of arc section for involute:
D point rectangular co-ordinate X values XD '=XO2 '+TECR2cos (THTP+90),
D point rectangular co-ordinate Y values YD '=XD ' tan (THTP),
D point polar diameters
Step 3.3, continuation solve each section of coordinate of key point of coil, and the key point number of the initial arc section of involute is
N1, the number of involute section key point is N2, and involute terminates the key point number of arc section for N3, N1=5, N2=20, N3=5,
Then have:
1) cycle calculations obtain the coordinate of the initial each key point in arc section inside of involute respectively, wherein:
The N1 polar diameter of key point is RB ', and polar angle is THTB '-THTA ', the rectangular co-ordinate X of remaining k-th key point
Value, K=1,2 ... ..., (N1-1):
XTEMP=XO1 '+TECR1cos (TEMPTHTA+TEMPTHTABK/N1),
The rectangular co-ordinate Y value of k-th key point:
YTEMP=YO1 '+TECR1sin (TEMPTHTA+TEMPTHTABK/N1);
The polar diameter of k-th key point
Polar angle atan (YTEMP, XTEMP)-THTA ' of k-th key point, wherein:
2) cycle calculations obtain the coordinate of each key point of involute intersegmental part respectively, wherein:
The N2 polar diameter of key point is RC ', and polar angle is THTC '-THTA ', the rectangular co-ordinate X of remaining k-th key point
Value, K=1,2 ... ..., (N2-1):
XTEMP=TECRJcos (ANGB+TEMPBCK/N2)+
TECRJ (ANGB+TEMPBCK/N2) π/180sin (ANGB+TEMPBCK/N2) '
The rectangular co-ordinate Y value of k-th key point:
YTEMP=TECRJsin (ANGB+TEMPBCK/N2)-
TECRJ (ANGB+TEMPBCK/N2) π/180cos (ANGB+TEMPBCK/ IV 2) '
The polar diameter of k-th key point
Polar angle atan (YTEMP, XTEMP)-THTA ' of k-th key point, wherein:
TEMPBC=ANGC-ANGB:
3) cycle calculations obtain the coordinate that involute terminates each key point in arc section inside respectively, wherein:
The N3 polar diameter of key point is RD ', and polar angle is THTD '-THTA ', the rectangular co-ordinate X of remaining k-th key point
Value, K=1,2 ... ..., (N3-1):
XTEMP=XO2 '+TECR2cos (TEMPTHTC-TEMPTHTCDK/N3);
The rectangular co-ordinate Y value of k-th key point:
YTEMP=YO2 '+TECR2sin (TEMPTHTC-TEMPTHTCDK/N3);
The polar diameter of k-th key point
Polar angle atan (YTEMP, XTEMP)-THTA ' of k-th key point, wherein:
If negative value, then it is modified to
TEMPTHTCD=TEMPTHTC-TEMPTHTD;
If negative value, then it is modified to
Step 3.4, the polar coordinates of each key point of coil are converted into cartesian coordinate using cycle calculations, wherein, k-th
KPCS (Isosorbide-5-Nitrae+K)=TECITEMPRsin (TEMPTHT) of key point,
KPCS (2,4+K)=TEMPRcos (TEMPTHT),
KPCS (3,4+K)=ICS (1, K) cos (CONTP)+OZ, wherein:
TEMPR=ICS (1, K) sin (CONTP), ICS (1, K) are the polar diameter of k-th key point;TEMPTHT=ICS (2,
K)/sin (CONTP), ICS (2, K) are the polar angle of k-th key point;
Step 3.5, set up key point:According to each crucial point coordinates that draws is solved in step 3.4, each key point is set up;
Step 3.6, the SPL for setting up coil:The key point set up in Connection Step 3.5, that is, generate the batten of coil
Curve;
Step 4, set up coil FEM model.
Claims (1)
1. a kind of parametric modeling method suitable for Generator Stator Winding Ends coil, it is characterised in that step is:
Step 1, input coil geometric parameter, at least including coil section overall width ECBW, coil section total height ECBH, coil
Section copper cash partial width ECBCW, coil section copper cash Partial Height ECBCH, coil section hollow copper wire hole width ECBCHW,
Coil section hollow copper string holes height ECBCHH, coil to core center linear distance TECRH, coil base radius TECRJ, coil
Straight line terminates place to the excessive chamfer radius TECR0 of the conical surface, coil involute beginning chamfer radius TECR1, coil involute
Chamfer radius TECR2, coil go out groove straight line portion length TECHZ, coil involute and start preceding length of straigh line TECBH1, coil
The involute that involute terminates across angle TECBYY, P point of rear length of straigh line TECBH2, coil conical degree of conical surface CONTP, coil launches
The involute angle of spread DEFAY of difference DEFAX, the C point of angle and polar angle, wherein, P points were the coil of O points in polar coordinates
Under involute launch after the intersection point for gradually opening the tangent line and the involute that terminate arc section, O points be coil under polar coordinates gradually
The basic circle center of circle burst at the seams after launching, C points are the terminal of the involute section after involute of the coil under polar coordinates launches;
Step 2, the cross sectional shape for defining coil;
Step 3, the modeling of stator winding end coil launcher, including:
Step 3.1 solves coil key point coordinates:
The geometrical curve of coil is formed by connecting by multiple key points, and the coordinate of each key point is stored in aray variable KPCS (i, j)
In, wherein i is coordinate components, i=1, and 2,3 represent the X-axis of key point, Y-axis and Z axis coordinate respectively, and j is key point numbering, fixed
Adopted aray variable KPCS (m, n), wherein m be coordinate components, m=1,2 respectively represent key point polar polar diameter R and polar angle
φ, n are numbered for key point, and the formula for calculating crucial point coordinates at each section of coil is as follows:
Coil beginning straightway key point coordinates:
Starting point coordinate:KPCS (1,1)=0,
KPCS (2,1)=TECRH,
KPCS (3,1)=0;
Terminal point coordinate:KPCS (1,2)=0,
KPCS (2,2)=TECRH,
KPCS (3,2)=TECHZ;
The crucial point coordinates of coil cone angle transition circle section:
Starting point coordinate:KPCS (1,3)=0,
KPCS (2,3)=TECRH+TECR0-TECR0cos (CONTP),
KPCS (3,3)=TECHZ+TECR0sin (CONTP);
Terminal point coordinate:KPCS (Isosorbide-5-Nitrae)=0,
KPCS (2,4)=TECRH+TECR0-TECR0cos (CONTP)+TECBH1sin (CONTP),
KPCS (3,4)=TECHZ+TECR0sin (CONTP)+TECBH1cos (CONTP);
The coordinate of the initial circular arc waypoint of involute:
1) A point coordinates, A points are the starting point of the initial arc section of involute:
A point polar diametersA point polar angles THTA=45;
A points rectangular co-ordinate X values XA=RAcos (THTA), A point rectangular co-ordinate Y value YA=RAsin (THTA), vertex of a cone Z coordinate
OZ=KPCS (3,4)-RAcos (CONTP);
2) coordinate of fillet center of circle O1 is sought according to A point coordinates, O1 is the center of circle of the initial arc section of involute:
O1 points rectangular co-ordinate X values XO1=XA+TECR1cos (THTA+90), O1 point rectangular co-ordinate Y value YO1=YA+TECR1
sin(THTA+90);
O1 point polar diametersO1 point polar angles
3) E point coordinates, E points are the point of contact of O1 points and basic circle tangent line:
E points rectangular co-ordinate X values XE=TECRJcos (THTE), E point rectangular co-ordinate Y value YE=TECRJsin (THTE), its
In, THTE=THTO1+THTEOO1,
4) B point coordinates, B points are the starting point of involute section:
B points rectangular co-ordinate X values XB=XO1+TECR1cos (THTO1E), B point rectangular co-ordinate Y value
YB=YO1+TECR1sin (THTO1E), wherein,
B point polar diametersB point polar angles
The involute angle of spread of B points:
5) angle of spread of the coil in sector:
TECBEA=TECBYYsin (CONTP);
After step 3.2, zero point of launching polar coordinates zero point and involute overlap, each point coordinates is solved again:
1) A point coordinates:
A points rectangular co-ordinate X values XA'=RAcos (THTA'), A point rectangular co-ordinate Y value
YA'=RAsin (THTA'), wherein, THTA'=THTA+ANGB-THTB;
2) O1 point coordinates:
O1 points rectangular co-ordinate X values XO1'=XA'+TECR1cos (THTA'+90), O1 point rectangular co-ordinate Y value
YO1'=YA'+TECR1sin (THTA'+90);
O1 point polar angles
3) E point coordinates:
E points rectangular co-ordinate X values XE'=TECRJcos (THTE'),
E points rectangular co-ordinate Y value YE'=TECRJsin (THTE'), wherein,
THTE'=THTO1'+THTEOO1;
4) B point coordinates:
B points rectangular co-ordinate X values XB'=XO1'+TECR1cos (THTO1E'),
B points rectangular co-ordinate Y value YB'=YO1'+TECR1sin (THTO1E'),
Wherein,
B point polar anglesB point polar diameters
5) P point coordinates:
P point polar angle ANGP=DEFAX+THTP, wherein, THTP=TECBEA+THTA';
6) C point coordinates:
C point rectangular co-ordinate X values
XC'=TECRJcos (ANGC)+TECRJANGC π/180sin (ANGC),
C point rectangular co-ordinate Y values
YC'=TECRJsin (ANGC)-TECRJANGC π/180cos (ANGC),
Wherein, ANGC=DEFAY;
C point polar diametersC point polar angles
7) O2 point coordinates, O2 points terminate the center of circle of arc section for involute:
O2 point rectangular co-ordinate X values:
XO2'=TECRJcos (ANGC)+[TECRJANGC π/180SIN (ANGC)] sin (ANGC);
O2 point rectangular co-ordinate Y values:
YO2'=TECRJsin (ANGC)-[TECRJANGC π/180COS (ANGC)] cos (ANGC);
8) D point coordinates, D points terminate the terminal of arc section for involute:
D points rectangular co-ordinate X values XD'=XO2'+TECR2cos (THTP+90),
D points rectangular co-ordinate Y value YD'=XD'tan (THTP),
D point polar diameters
Step 3.3, continuation solve each section of coordinate of key point of coil, and the key point number of the initial arc section of involute is N1, gradually
The number of open segment key point is N2, and the key point number that involute terminates arc section is N3, then have:
1) cycle calculations obtain the coordinate of the initial each key point in arc section inside of involute respectively, wherein:
The N1 polar diameter of key point is RB', and polar angle is THTB'-THTA', the rectangular co-ordinate X values of remaining k-th key point, K
=1,2 ... ..., (N1-1):
XTEMP=XO1'+TECR1cos (TEMPTHTA+TEMPTHTABK/N1),
The rectangular co-ordinate Y value of k-th key point:
YTEMP=YO1'+TECR1sin (TEMPTHTA+TEMPTHTABK/N1);
The polar diameter of k-th key point
Polar angle atan (YTEMP, XTEMP)-THTA' of k-th key point, wherein:
2) cycle calculations obtain the coordinate of each key point of involute intersegmental part respectively, wherein:
The N2 polar diameter of key point is RC', and polar angle is THTC'-THTA', the rectangular co-ordinate X values of remaining k-th key point, K
=1,2 ... ..., (N2-1):
XTEMP=TECRJcos (ANGB+TEMPBCK/N2)+
TECRJ·(ANGB+TEMPBC·K/N2)·π/180·sin(ANGB+TEMPBC·K/N2)
The rectangular co-ordinate Y value of k-th key point:
YTEMP=TECRJsin (ANGB+TEMPBCK/N2)-
TECRJ·(ANGB+TEMPBC·K/N2)·π/180·cos(ANGB+TEMPBC·K/N2)
The polar diameter of k-th key point
Polar angle atan (YTEMP, XTEMP)-THTA' of k-th key point, wherein:
TEMPBC=ANGC-ANGB;
3) cycle calculations obtain the coordinate that involute terminates each key point in arc section inside respectively, wherein:
The N3 polar diameter of key point is RD', and polar angle is THTD'-THTA', the rectangular co-ordinate X values of remaining k-th key point, K
=1,2 ... ..., (N3-1):
XTEMP=XO2'+TECR2cos (TEMPTHTC-TEMPTHTCDK/N3);
The rectangular co-ordinate Y value of k-th key point:
YTEMP=YO2'+TECR2sin (TEMPTHTC-TEMPTHTCDK/N3);
The polar diameter of k-th key point
Polar angle atan (YTEMP, XTEMP)-THTA' of k-th key point, wherein:
If negative value, then it is modified to
TEMPTHTCD=TEMPTHTC-TEMPTHTD;
If negative value, then it is modified to
Step 3.4, the polar coordinates of each key point of coil are converted into cartesian coordinate using cycle calculations, wherein, k-th is crucial
KPCS (Isosorbide-5-Nitrae+K)=TECRJTEMPRsin (TEMPTHT) of point,
KPCS(2,4+K)=TEMPR·cos(TEMPTHT),
KPCS (3,4+K)=ICS (1, K) cos (CONTP)+OZ, wherein:
TEMPR=ICS (1, K) sin (CONTP), ICS (1, K) are the polar diameter of k-th key point;TEMPTHT=ICS (2,
K)/sin (CONTP), ICS (2, K) are the polar angle of k-th key point;
Step 3.5, set up key point:According to each crucial point coordinates that draws is solved in step 3.4, each key point is set up;
Step 3.6, the SPL for setting up coil:The key point set up in Connection Step 3.5, that is, the batten for generating coil is bent
Line;
Step 4, set up coil FEM model.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310574323.7A CN103559365B (en) | 2013-11-15 | 2013-11-15 | A kind of parametric modeling method suitable for Generator Stator Winding Ends coil |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310574323.7A CN103559365B (en) | 2013-11-15 | 2013-11-15 | A kind of parametric modeling method suitable for Generator Stator Winding Ends coil |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103559365A CN103559365A (en) | 2014-02-05 |
CN103559365B true CN103559365B (en) | 2017-07-07 |
Family
ID=50013611
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310574323.7A Active CN103559365B (en) | 2013-11-15 | 2013-11-15 | A kind of parametric modeling method suitable for Generator Stator Winding Ends coil |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103559365B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105373675B (en) * | 2015-12-09 | 2018-09-14 | 河海大学 | A kind of SRM stator modes finite element modeling method |
CN108875104A (en) * | 2017-05-12 | 2018-11-23 | 上海电气集团股份有限公司 | A kind of binding-type generator stator end finite element modeling method |
CN109586448A (en) * | 2017-09-29 | 2019-04-05 | 比亚迪股份有限公司 | Conductor segment and stator module, motor with it |
CN108539883B (en) * | 2018-04-25 | 2024-02-02 | 上海法雷奥汽车电器系统有限公司 | Vehicle generator, method for determining height of stator holes of vehicle generator and stator |
CN111404333B (en) * | 2020-01-17 | 2021-07-02 | 华中科技大学 | Method and system for acquiring electromagnetic force waveform of motor end winding |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101976285A (en) * | 2010-10-26 | 2011-02-16 | 哈尔滨工业大学 | Parametric design method for turbonator key components |
CN103034756A (en) * | 2012-11-30 | 2013-04-10 | 上海电气电站设备有限公司 | Modularized rotor data generating method used in calculation of shaft system of power generator |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7844436B2 (en) * | 2008-12-18 | 2010-11-30 | Telefonaktiebolaget L M Ericsson (Publ) | Method and localization unit for detecting and locating load coils in a transmission line |
-
2013
- 2013-11-15 CN CN201310574323.7A patent/CN103559365B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101976285A (en) * | 2010-10-26 | 2011-02-16 | 哈尔滨工业大学 | Parametric design method for turbonator key components |
CN103034756A (en) * | 2012-11-30 | 2013-04-10 | 上海电气电站设备有限公司 | Modularized rotor data generating method used in calculation of shaft system of power generator |
Non-Patent Citations (2)
Title |
---|
CAD/CAE在大型汽轮发电机设计研发中的应用;黄磊,等.;《第三届中国CAE工程分析技术年会》;20070728;284-292 * |
定子端部线圈的设计计算方法与三维建模研究;胡建波,等.;《上海大中型电机》;20061231(第4期);14-18 * |
Also Published As
Publication number | Publication date |
---|---|
CN103559365A (en) | 2014-02-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103559365B (en) | A kind of parametric modeling method suitable for Generator Stator Winding Ends coil | |
CN103577654B (en) | A kind of finite element Precise modeling of large turbo-type generator stator bar | |
TW212257B (en) | ||
CN109255174B (en) | Magnetic coupling resonant wireless energy transmission coil simulation analysis method | |
JP2016525860A (en) | Windings of rotating electrical machines and methods for designing such windings | |
JP2013198398A (en) | Hollow cylindrical air-core winding | |
CN110571962A (en) | Flat wire motor stator | |
JP4809118B2 (en) | Method for manufacturing heating coil, heating coil for electromagnetic induction heating cooker, and electromagnetic induction heating cooker | |
JP4603810B2 (en) | Method for making two-layer lap winding | |
CN104599681A (en) | Method and device for processing audio file | |
CN206041666U (en) | Motor stator convenient to wire winding | |
CN103699752A (en) | Coupling method for processing moving boundary problems in electromagnetic field based on edge element method | |
CN105373675A (en) | SRM stator modal finite element modeling method | |
CN105391261A (en) | High-speed non-salient-pole electrically excited synchronous motor rotor in air gap magnetic field sine distribution and structural parameter determination method of rotor | |
CN109239544A (en) | A kind of oscillatory wave system air reactor design method | |
CN110321628B (en) | Design method of wearable device wireless charging system based on curved surface flexible coil | |
CN201707817U (en) | Electromagnetic induction demonstrator | |
CN103092130B (en) | Piston outer circle modeling method | |
US20110273051A1 (en) | Stator Core Winding Method for Motor and Structure Thereof | |
Serrano et al. | Mathematical description of PCB-adapted litz wire geometry for automated layout generation of WPT coils | |
CN115331927A (en) | Semiconductor structure, preparation method and semiconductor device | |
CN104272407B (en) | The flat copper winding that magnetic field is generated in electric transducer with high fill-factor | |
CN104022716A (en) | Method for controlling maximum torque current ratio of excitation-variation synchronous motor through coefficient fitting | |
JP2010055427A (en) | Simulation method for induction heating coil design support | |
CN114050046B (en) | Design method of wire gauge of wireless power transmission coil with double wire widths |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |