CN113420481B - Material characteristic calculation method for back field magnet winding - Google Patents
Material characteristic calculation method for back field magnet winding Download PDFInfo
- Publication number
- CN113420481B CN113420481B CN202110698446.6A CN202110698446A CN113420481B CN 113420481 B CN113420481 B CN 113420481B CN 202110698446 A CN202110698446 A CN 202110698446A CN 113420481 B CN113420481 B CN 113420481B
- Authority
- CN
- China
- Prior art keywords
- nodes
- conductor
- selecting
- applying
- insulator
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/16—Cables, cable trees or wire harnesses
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention discloses a material characteristic calculation method of a back surface field magnet winding, which has the technical scheme that: the method comprises the following steps: s01, establishing a parameterized conductor section model; s02, creating a three-dimensional finite element model of the conductor; s03, applying a method for calculating the elastic modulus E of the conductor x And poisson ratio PR xy ,PR xz Boundary conditions of (2); s04, calculating the elastic modulus E of the conductor x And poisson ratio PR xy ,PR xz The method comprises the steps of carrying out a first treatment on the surface of the S05, applying the elastic modulus E for calculating the conductor y And poisson ratio PR yz Boundary conditions of (2); s06, calculating the elastic modulus E of the conductor y And poisson ratio PR yz The method comprises the steps of carrying out a first treatment on the surface of the S07, applying the elastic modulus E for calculating the conductor z Boundary conditions of (2); s08, calculating elastic modulus E of the conductor z The method comprises the steps of carrying out a first treatment on the surface of the S09, applying a method for calculating the shear modulus G of the conductor xy Boundary conditions of (2); s10, calculating the shear modulus G of the conductor xy The method comprises the steps of carrying out a first treatment on the surface of the The material characteristic calculation method of the back field magnet winding can rapidly develop mechanical property evaluation and structure optimization design of the back field magnet component, and greatly improve analysis and calculation efficiency of the back field magnet.
Description
Technical Field
The invention relates to a material characteristic calculation method, in particular to a material characteristic calculation method of a back surface field magnet winding.
Background
The superconducting cable passes through a back field magnet winding formed by multi-stage cable penetrating and winding, and the high-temperature superconducting material has the characteristics of no resistance and high critical current density, so that the back field magnet winding becomes a technical break for power application. However, the material property solving of the high-temperature superconducting material is still an unsolved problem at present.
Disclosure of Invention
In view of the problems mentioned in the background art, an object of the present invention is to provide a method for calculating material characteristics of a back field magnet winding, so as to solve the problems mentioned in the background art.
The technical aim of the invention is realized by the following technical scheme:
a material characteristic calculation method of a back field magnet winding is characterized by comprising the following steps of: the method comprises the following steps:
s01, establishing a parameterized conductor section model;
s02, creating a three-dimensional finite element model of the conductor;
s03, applying a method for calculating the elastic modulus E of the conductor x And poisson ratio PR xy ,PR xz Boundary conditions of (2);
s04, calculating the elastic modulus E of the conductor x And poisson ratio PR xy ,PR xz ;
S05, applying the elastic modulus E for calculating the conductor y And poisson ratio PR yz Boundary conditions of (2);
s06, calculating the elastic modulus E of the conductor y And poisson ratio PR yz ;
S07, applying the elastic modulus E for calculating the conductor z Boundary conditions of (2);
s08, calculating elastic modulus E of the conductor z ;
S09, applying a method for calculating the shear modulus G of the conductor xy Boundary conditions of (2);
s10, calculating the shear modulus G of the conductor xy ;
S11, applying the method for calculatingConductor shear modulus G xz Boundary conditions of (2);
s12, calculating the shear modulus G of the conductor xz ;
S13, applying a method for calculating the shear modulus G of the conductor yz Boundary conditions of (2);
s14, calculating the shear modulus G of the conductor yz 。
Preferably, the S01 includes:
s011, defining variable parameters based on ANSYS parameterization language to describe the thickness of the conductor insulating layer, the width, the height, the thickness and the chamfer radius of the armor, and preparing for creating the conductor insulating layer and the armor;
s012, establishing key points of the conductor according to the dimensional parameters of each component of the conductor in a Cartesian coordinate system and a cylindrical coordinate system, and then respectively establishing a line model and a surface model of the conductor section by a bottom-up modeling method.
Preferably, the S02 includes:
s021, carrying out Boolean operation on the surface model of the conductor to enable the topological structures of all the sub-surfaces to meet the requirement of mapping grid division;
s022, respectively defining material characteristics such as elastic modulus, shear modulus, poisson ratio and the like for the conductor armor and the insulating part;
s023, defining a MESH200 unit, and carrying out mapping grid division on the surface model of the conductor to generate a two-dimensional surface grid model;
s024, creating a sweep line perpendicular to the section of the conductor, and generating a three-dimensional conductor finite element model through sweep grid division.
Preferably, the S03 includes:
s031, selecting three groups of nodes satisfying x=0, y=0, and z=0, respectively, and applying a fixed constraint;
s032, respectively selecting three groups of nodes with x, y and z coordinates equal to the width, the height and the length of the insulator, and respectively applying degree-of-freedom coupling constraint;
s033, selecting a group of nodes with x coordinates equal to the width of the insulator, and applying a pressure load of 1 MPa;
the S04 includes:
s041, solving a finite element model after boundary conditions are applied;
s042 extracting average deformation of the conductor in x, y and z directionsAnd calculates the corresponding line strain SN x ,SN y ,SN z ;
S043, selecting a group of nodes with x coordinates equal to the width of the insulator, and extracting the average stress S of the nodes x ;
S044, calculate S x /SN x ,SN y /SN x ,SN z /SN x And assign them to variables E x ,PR xy ,PR xz 。
Preferably, the S05 includes:
s051, deleting the load boundary conditions in the S033;
s052, selecting a group of nodes with y coordinates equal to the height of the insulator, and applying a pressure load of 1 MPa;
the S06 includes:
s061, solving a finite element model after applying boundary conditions;
s062 extracting average deformation of conductorAnd calculates the corresponding line strain SN y ,SN z ;
S063, selecting a group of nodes with y coordinates equal to the height of the insulator, and extracting the average stress S of the nodes y ;
S064, calculate S y /SN y ,SN y /SN z And assign them to variables E y ,μ yz 。
Preferably, the S07 includes:
s071 deleting the load boundary condition in S051;
s072, selecting a group of nodes with z coordinates equal to the length of the insulator, and applying a pressure load of 1 MPa;
the S08 includes:
s081, solving a finite element model after boundary conditions are applied;
s082, extracting average deformation of conductorAnd calculates the corresponding line strain SN z ;
S083, selecting a group of nodes with z coordinates equal to the length of the insulator, and extracting average stress S of the nodes z ;
S084, calculate S z /SN z And assign it to variable E z 。
Preferably, the S09 includes:
s091, deleting boundary conditions;
s092, selecting a set of nodes of the insulator surface with x=0, and applying a displacement constraint ux=0, uy=0;
s093, select x A set of nodes having coordinates equal to the width of the insulator and coupling the degrees of freedom of the direction of the set of nodes x;
s094 selecting a set of nodes of the insulator surface with z=0 and applying a displacement constraint uz=0;
s095, selecting a group of nodes with z coordinates equal to the length of the insulator, and coupling the degrees of freedom of the direction of the group of nodes z;
s096, selecting a node of the insulator surface with z=0, selecting a node of which the node coordinates are (0, 0) from the node group, assigning the number of the node to the variable na, selecting a node of which the node coordinates are (y=0, z=0) from the node group, and the x-coordinate is equal to the insulator width, assigning the number of the node to the variable nb, assigning the number of any of the remaining nodes to the variable nc, and imposing a constraint on the nodes na, nb, nc;
s097, selecting a set of nodes with x-coordinates equal to the insulator width, and applying a load of 1N to the nodes in the y-direction.
Preferably, the S10 includes:
s101, solving a finite element model applying boundary conditions in S09;
s102, extracting shear stress tau in xy plane xy Shear strain gamma xy ;
S103, τ is set xy /γ xy Assigned to variable G xy ;
The step S11 includes:
s111, deleting boundary conditions;
s112, selecting a set of nodes of the insulator surface with x=0, and applying a displacement constraint ux=0, uz=0;
s113, selecting a group of nodes with x coordinates equal to the width of the insulator, and coupling the degrees of freedom of the group of nodes in the x direction;
s114, selecting a set of nodes of the insulator surface with y=0, and applying a displacement constraint uy=0;
s115, selecting a group of nodes with y coordinates equal to the height of the insulator, and coupling the degrees of freedom of the group of nodes in the y direction;
s116, selecting a group of nodes on the insulator surface with z=0, selecting a node with a node coordinate of (0, 0) from the group of nodes, assigning the number of the node to a variable na, selecting a node with a node coordinate of (y=0, z=0) from the group of nodes, and an x coordinate equal to the width of the insulator, assigning the number of the node to a variable nb, assigning the number of any of the remaining nodes to a variable nc, and applying a constraint to the nodes na, nb, nc;
s117, selecting a group of nodes with x coordinates equal to the width of the insulator, and applying a load of 1N to the nodes along the z direction.
Preferably, the S12 includes:
s121, solving the finite element model applying the boundary condition in the S11;
s122, extracting shear stress tau in xy plane xz Shear strain gamma xz ;
S123, will τ xz /γ xz Assigned to variable G xz ;
The step S13 includes:
s131, deleting the boundary conditions in the S11;
s132, selecting a set of nodes of the insulator surface with y=0, and applying a displacement constraint uy=0, uz=0;
s133, selecting a group of nodes with x coordinates equal to the width of the insulator, and coupling the degrees of freedom of the group of nodes in the x direction;
s134, selecting a set of nodes of the insulator surface with x=0 coordinates equal to zero, and applying a displacement constraint ux=0;
s135, selecting a group of nodes with y coordinates equal to the height of the insulator, and coupling the degrees of freedom of the group of nodes in the y direction;
s136, selecting a group of nodes on the insulator surface with x=0, selecting a node with a node coordinate of (0, 0) from the node group, assigning the number of the node to a variable na, selecting a node with a (x, z) coordinate of (0, 0) from the node group and a y coordinate equal to the height of the insulator, assigning the number of the node to a variable nb, assigning the number of any other node to a variable nc, and applying constraint to the nodes na, nb, nc;
s137, selecting a group of nodes with y coordinates equal to the height of the insulator, and applying a load of 1N to the nodes along z.
Preferably, the S14 includes:
s141, solving the finite element model applying the boundary condition in the S13;
s142, extracting shear stress tau in yz plane yz Shear strain gamma yz ;
S143, will τ yz /γ yz Assigned to variable G yz 。
In summary, the invention has the following advantages:
the material characteristic calculation method of the back field magnet winding can rapidly develop mechanical property evaluation and structure optimization design of the back field magnet part, greatly improve analysis and calculation efficiency of the back field magnet, firstly establish relative position coordinates of key points of a conductor section under a Cartesian coordinate system and a cylindrical coordinate system respectively by using analysis software according to the section size of the single-turn conductor of the magnet, sequentially generate a conductor section contour line, a conductor surface model and a three-dimensional solid model from bottom to top by the key points, then divide a grid for the conductor, add constraint equations and loads and solve the grids, finally extract stress strain of the conductor and reversely calculate the material properties of the conductor by combining the constitutive equations of the conductor, and realize full-parametric driving of the single-turn conductor modeling of the magnet.
Drawings
FIG. 1 is a flow chart of the present invention;
fig. 2 is a structural diagram of the cross-sectional dimensions of a back field magnet winding unit cell.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Example 1
Referring to fig. 1 and 2, a method for calculating material characteristics of a back field magnet winding includes the steps of:
s01, establishing a parameterized conductor section model;
s02, creating a three-dimensional finite element model of the conductor;
s03, applying a method for calculating the elastic modulus E of the conductor x And poisson ratio PR xy ,PR xz Boundary conditions of (2);
s04, calculating the elastic modulus E of the conductor x And poisson ratio PR xy ,PR xz ;
S05, applying the elastic modulus E for calculating the conductor y And poisson ratio PR yz Boundary conditions of (2);
s06, calculating the elastic modulus E of the conductor y And poisson ratio PR yz ;
S07, applying the elastic modulus E for calculating the conductor z Boundary conditions of (2);
s08, calculating elastic modulus E of the conductor z ;
S09, applying a method for calculating the shear modulus G of the conductor xy Boundary conditions of (2);
s10, calculating the shear modulus G of the conductor xy ;
S11, applying a method for calculating the shear modulus G of the conductor xz Boundary conditions of (2);
s12, calculating the shear modulus G of the conductor xz ;
S13, applying a method for calculating the shear modulus G of the conductor yz Boundary conditions of (2);
s14, calculating the shear modulus G of the conductor yz 。
Wherein S01 includes:
s011, defining variable parameters based on ANSYS parameterization language to describe the thickness of the conductor insulating layer, the width, the height, the thickness and the chamfer radius of the armor, and preparing for creating the conductor insulating layer and the armor;
s012, establishing key points of the conductor according to the dimensional parameters of each component of the conductor in a Cartesian coordinate system and a cylindrical coordinate system, and then respectively establishing a line model and a surface model of the conductor section by a bottom-up modeling method.
Wherein S02 includes:
s021, carrying out Boolean operation on the surface model of the conductor to enable the topological structures of all the sub-surfaces to meet the requirement of mapping grid division;
s022, respectively defining material characteristics such as elastic modulus, shear modulus, poisson ratio and the like for the conductor armor and the insulating part, wherein the definition is shown in table 1;
table 1: material properties of the conductor Assembly (temperature: 4K)
S023, defining a MESH200 unit, and carrying out mapping grid division on the surface model of the conductor to generate a two-dimensional surface grid model;
s024, creating a sweep line perpendicular to the section of the conductor, and generating a three-dimensional conductor finite element model through sweep grid division.
Wherein S03 includes:
s031, selecting three groups of nodes satisfying x=0, y=0, and z=0, respectively, and applying a fixed constraint;
s032, respectively selecting three groups of nodes with x, y and z coordinates equal to the width, the height and the length of the insulator, and respectively applying degree-of-freedom coupling constraint;
s033, selecting a group of nodes with x coordinates equal to the width of the insulator, and applying a pressure load of 1 MPa;
s04 includes:
s041, solving a finite element model after boundary conditions are applied;
s042 extracting average deformation of the conductor in x, y and z directionsAnd calculates the corresponding line strain SN x ,SN y ,SN z ;
S043, selecting a group of nodes with x coordinates equal to the width of the insulator, and extracting the average stress S of the nodes x ;
S044, calculate S x /SN x ,SN y /SN x ,SN z /SN x And assign them to variables E x ,PR xy ,PR xz 。
Wherein S05 includes:
s051, deleting the load boundary conditions in S033;
s052, selection y A set of nodes with coordinates equal to the insulator height and applying a pressure load of 1 MPa;
wherein S06 includes:
s061, solving a finite element model after applying boundary conditions;
s062 extracting average deformation of conductorAnd calculates the corresponding line strain SN y ,SN z ;
S063, selecting a group of nodes with y coordinates equal to the height of the insulator, and extracting the average stress S of the nodes y ;
S064, calculate S y /SN y ,SN y /SN z And assign them to variables E y ,μ yz 。
Wherein S07 includes:
s071, deleting the load boundary condition in S051;
s072, selecting a group of nodes with z coordinates equal to the length of the insulator, and applying a pressure load of 1 MPa;
wherein S08 includes:
s081, solving a finite element model after boundary conditions are applied;
s082, extracting average deformation of conductorAnd calculates the corresponding line strain SN z ;
S083, selecting a group of nodes with z coordinates equal to the length of the insulator, and extracting average stress S of the nodes z ;
S084, calculate S z /SN z And assign it to variable E z 。
Wherein S09 includes:
s091, deleting boundary conditions;
s092, selecting a set of nodes of the insulator surface with x=0, and applying a displacement constraint ux=0, uy=0;
s093, selecting a group of nodes with x coordinates equal to the width of the insulator, and coupling the direction degrees of freedom of the group of nodes x;
s094 selecting a set of nodes of the insulator surface with z=0 and applying a displacement constraint uz=0;
s095, selecting a group of nodes with z coordinates equal to the length of the insulator, and coupling the degrees of freedom of the direction of the group of nodes z;
s096, selecting a node of the insulator surface with z=0, selecting a node of which the node coordinates are (0, 0) from the node group, assigning the number of the node to the variable na, selecting a node of which the node coordinates are (y=0, z=0) from the node group, and the x-coordinate is equal to the insulator width, assigning the number of the node to the variable nb, assigning the number of any of the remaining nodes to the variable nc, and imposing a constraint on the nodes na, nb, nc;
s097, selecting a set of nodes with x-coordinates equal to the insulator width, and applying a load of 1N to the nodes in the y-direction.
Wherein S10 includes:
s101, solving a finite element model applying boundary conditions in S09;
s102, extracting shear stress tau in xy plane xy Shear strain gamma xy ;
S103, τ is set xy /γ xy Assigned to variable G xy ;
Wherein S11 includes:
s111, deleting boundary conditions;
s112, selecting a set of nodes of the insulator surface with x=0, and applying a displacement constraint ux=0, uz=0;
s113, selecting a group of nodes with x coordinates equal to the width of the insulator, and coupling the degrees of freedom of the group of nodes in the x direction;
s114, selecting a set of nodes of the insulator surface with y=0, and applying a displacement constraint uy=0;
s115, selecting a group of nodes with y coordinates equal to the height of the insulator, and coupling the degrees of freedom of the group of nodes in the y direction;
s116, selecting a group of nodes on the insulator surface with z=0, selecting a node with a node coordinate of (0, 0) from the group of nodes, assigning the number of the node to a variable na, selecting a node with a node coordinate of (y=0, z=0) from the group of nodes, and an x coordinate equal to the width of the insulator, assigning the number of the node to a variable nb, assigning the number of any of the remaining nodes to a variable nc, and applying a constraint to the nodes na, nb, nc;
s117, selecting a group of nodes with x coordinates equal to the width of the insulator, and applying a load of 1N to the nodes along the z direction.
Wherein S12 includes:
s121, solving a finite element model applying boundary conditions in S11;
s122, extracting shear stress tau in xy plane xz Shear strain gamma xz ;
S123, will τ xz /γ xz Assigned to variable G xz ;
Wherein S13 includes:
s131, deleting the boundary conditions in S11;
s132, selecting a set of nodes of the insulator surface with y=0, and applying a displacement constraint uy=0, uz=0;
s133, selecting a group of nodes with x coordinates equal to the width of the insulator, and coupling the degrees of freedom of the group of nodes in the x direction;
s134, selecting a set of nodes of the insulator surface with x=0 coordinates equal to zero, and applying a displacement constraint ux=0;
s135, selecting a group of nodes with y coordinates equal to the height of the insulator, and coupling the degrees of freedom of the group of nodes in the y direction;
s136, selecting a group of nodes on the insulator surface with x=0, selecting a node with a node coordinate of (0, 0) from the node group, assigning the number of the node to a variable na, selecting a node with a (x, z) coordinate of (0, 0) from the node group and a y coordinate equal to the height of the insulator, assigning the number of the node to a variable nb, assigning the number of any other node to a variable nc, and applying constraint to the nodes na, nb, nc;
s137, selecting a group of nodes with y coordinates equal to the height of the insulator, and applying a load of 1N to the nodes along z, wherein the applied boundary conditions are as shown in table 2:
TABLE 2 load and boundary conditions
Wherein S14 includes:
s141, solving a finite element model applying boundary conditions in S13;
s142, extracting shear stress tau in yz plane yz Shear strain gamma yz ;
S143, will τ yz /γ yz Assigned to variable G yz The results of the calculations are shown in Table 3:
TABLE 3 Back field magnet conductor winding Material Properties
Referring to fig. 1 and 2, the material characteristic calculation method of the back field magnet winding can rapidly develop mechanical property evaluation and structure optimization design of the back field magnet component, greatly improve analysis and calculation efficiency of the back field magnet, firstly, according to the section size of a single-turn conductor of the magnet, relative position coordinates of key points of the section of the conductor are created under a Cartesian coordinate system and a cylindrical coordinate system respectively by using analysis software, a conductor section contour line, a conductor surface model and a three-dimensional entity model are sequentially generated from bottom to top by the key points, then grids are divided on the conductor, constraint equations and loads are added and solved, finally stress strain of the conductor is extracted, and the material properties of the conductor are reversely calculated by combining the constitutive equations of the conductor, so that full-parametric driving of single-turn conductor modeling of the magnet is realized.
Although embodiments of the present invention have been shown and described, it will be understood by those skilled in the art that various changes, modifications, substitutions and alterations can be made therein without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.
Claims (9)
1. A material characteristic calculation method of a back field magnet winding is characterized by comprising the following steps of: the method comprises the following steps:
s01, establishing a parameterized conductor section model;
s02, creating a three-dimensional finite element model of the conductor;
s03, applying a method for calculating the elastic modulus E of the conductor x And poisson ratio PR xy ,PR xz Boundary conditions of (2);
s04, calculating the elastic modulus E of the conductor x And poisson ratio PR xy ,PR xz ;
S05, applying the elastic modulus E for calculating the conductor y And poisson ratio PR yz Boundary conditions of (2);
s06, calculating the elastic modulus E of the conductor y And poisson ratio PR yz ;
S07, applying the elastic modulus E for calculating the conductor z Boundary conditions of (2);
s08, calculating elastic modulus E of the conductor z ;
S09, applying a method for calculating the shear modulus G of the conductor xy Boundary conditions of (2);
s10, calculating the shear modulus G of the conductor xy ;
S11, applying a method for calculating the shear modulus G of the conductor xz Boundary conditions of (2);
s12, calculating the shear modulus G of the conductor xz ;
S13, applying a method for calculating the shear modulus G of the conductor yz Boundary conditions of (2);
s14, calculating the shear modulus G of the conductor yz ;
The step S03 comprises the following steps:
s031, selecting three groups of nodes satisfying x=0, y=0, and z=0, respectively, and applying a fixed constraint;
s032, respectively selecting three groups of nodes with x, y and z coordinates equal to the width, the height and the length of the insulator, and respectively applying degree-of-freedom coupling constraint;
s033, selecting a group of nodes with x coordinates equal to the width of the insulator, and applying a pressure load of 1 MPa;
the S04 includes:
s041, solving a finite element model after boundary conditions are applied;
s042 extracting average deformation of the conductor in x, y and z directionsAnd calculates the corresponding line strain SN x ,SN y ,SN z ;
S043, selecting a group of nodes with x coordinates equal to the width of the insulator, and extracting the average stress S of the nodes x ;
S044, calculate S x /SN x ,SN y /SN x ,SN z /SN x And assign them to variables E x ,PR xy ,PR xz 。
2. The method for calculating the material characteristics of the back field magnet winding according to claim 1, wherein: the S01 includes:
s011, defining variable parameters based on ANSYS parameterization language to describe the thickness of the conductor insulating layer, the width, the height, the thickness and the chamfer radius of the armor, and preparing for creating the conductor insulating layer and the armor;
s012, establishing key points of the conductor according to the dimensional parameters of each component of the conductor in a Cartesian coordinate system and a cylindrical coordinate system, and then respectively establishing a line model and a surface model of the conductor section by a bottom-up modeling method.
3. The method for calculating the material characteristics of the back field magnet winding according to claim 1, wherein: the S02 includes:
s021, carrying out Boolean operation on the surface model of the conductor to enable the topological structures of all the sub-surfaces to meet the requirement of mapping grid division;
s022, respectively defining material characteristics of elastic modulus, shear modulus and Poisson ratio for the conductor armor and the insulating part;
s023, defining a MESH200 unit, and carrying out mapping grid division on the surface model of the conductor to generate a two-dimensional surface grid model;
s024, creating a sweep line perpendicular to the section of the conductor, and generating a three-dimensional conductor finite element model through sweep grid division.
4. The method for calculating the material characteristics of the back field magnet winding according to claim 1, wherein: the step S05 includes:
s051, deleting the load boundary conditions in the S033;
s052, selecting a group of nodes with y coordinates equal to the height of the insulator, and applying a pressure load of 1 MPa;
the S06 includes:
s061, solving a finite element model after applying boundary conditions;
s062 extracting average deformation of conductorAnd calculates the corresponding line strain SN y ,SN z ;
S063, selecting a group of nodes with y coordinates equal to the height of the insulator, and extracting the average stress S of the nodes y ;
S064, calculate S y /SN y ,SN y /SN z And assign them to variables E y ,μ yz 。
5. The method for calculating the material characteristics of the back field magnet winding according to claim 4, wherein: the S07 includes:
s071 deleting the load boundary condition in S051;
s072, selecting a group of nodes with z coordinates equal to the length of the insulator, and applying a pressure load of 1 MPa;
the S08 includes:
s081, solving a finite element model after boundary conditions are applied;
s082, extracting average deformation of conductorAnd calculates the corresponding line strain SN z ;
S083, selecting a group of nodes with z coordinates equal to the length of the insulator, and extracting average stress S of the nodes z ;
S084, calculate S z /SN z And assign it to variable E z 。
6. The method for calculating the material characteristics of the back field magnet winding according to claim 1, wherein: the S09 includes:
s091, deleting boundary conditions;
s092, selecting a set of nodes of the insulator surface with x=0, and applying a displacement constraint ux=0, uy=0;
s093, select x A set of nodes having coordinates equal to the width of the insulator and coupling the degrees of freedom of the direction of the set of nodes x;
s094 selecting a set of nodes of the insulator surface with z=0 and applying a displacement constraint uz=0;
s095, selecting a group of nodes with z coordinates equal to the length of the insulator, and coupling the degrees of freedom of the direction of the group of nodes z;
s096, selecting a node on the insulator surface with z=0, selecting a node with a node coordinate of (0, 0) from the node group of S094, assigning the number of the node to a variable na, selecting a node with a node coordinate of (y=0, z=0) from the node group and x-coordinate equal to the insulator width, assigning the number of the node to a variable nb, assigning the number of any of the remaining nodes to a variable nc, and imposing a constraint on the nodes na, nb, nc;
s097, selecting a set of nodes with x-coordinates equal to the insulator width, and applying a load of 1N to the nodes in the y-direction.
7. The method for calculating the material characteristics of the back field magnet winding according to claim 6, wherein: the S10 includes:
s101, solving a finite element model applying boundary conditions in S09;
s102, extracting shear stress tau in xy plane xy Shear strain gamma xy ;
S103, τ is set xy /γ xy Assigned to variable G xy ;
The step S11 includes:
s111, deleting boundary conditions;
s112, selecting a set of nodes of the insulator surface with x=0, and applying a displacement constraint ux=0, uz=0;
s113, selecting a group of nodes with x coordinates equal to the width of the insulator, and coupling the degrees of freedom of the group of nodes in the x direction;
s114, selecting a set of nodes of the insulator surface with y=0, and applying a displacement constraint uy=0;
s115, selecting a group of nodes with y coordinates equal to the height of the insulator, and coupling the degrees of freedom of the group of nodes in the y direction;
s116, selecting a node with a node coordinate of (0, 0) from a node group of the insulator surface with z=0, assigning the number of the node to a variable na, selecting a node with a coordinate of (y=0, z=0) and an x-coordinate equal to the insulator width from the node group, assigning the number of the node to a variable nb, assigning the number of any other node to a variable nc, and applying a constraint to the nodes na, nb, nc;
s117, selecting a group of nodes with x coordinates equal to the width of the insulator, and applying a load of 1N to the nodes along the z direction.
8. The method for calculating the material characteristics of the back field magnet winding according to claim 7, wherein: the S12 includes:
s121, solving the finite element model applying the boundary condition in the S11;
s122, extracting shear stress tau in xy plane xz Shear strain gamma xz ;
S123, will τ xz /γ xz Assigned to variable G xz ;
The step S13 includes:
s131, deleting the boundary conditions in the S11;
s132, selecting a set of nodes of the insulator surface with y=0, and applying a displacement constraint uy=0, uz=0;
s133, selecting a group of nodes with x coordinates equal to the width of the insulator, and coupling the degrees of freedom of the group of nodes in the x direction;
s134, selecting a set of nodes of the insulator surface with x=0 coordinates equal to zero, and applying a displacement constraint ux=0;
s135, selecting a group of nodes with y coordinates equal to the height of the insulator, and coupling the degrees of freedom of the group of nodes in the y direction;
s136, selecting a node with a node coordinate of (0, 0) from a node group of the insulator surface with x=0, assigning the number of the node to a variable na, selecting a node with a (x, z) coordinate of (0, 0) from the node group and a y coordinate equal to the height of the insulator, assigning the number of the node to a variable nb, assigning the number of any other node to a variable nc, and applying a constraint to the nodes na, nb, nc;
s137, selecting a group of nodes with y coordinates equal to the height of the insulator, and applying a load of 1N to the nodes along z.
9. The method for calculating the material characteristics of the back field magnet winding according to claim 8, wherein: the S14 includes:
s141, solving the finite element model applying the boundary condition in the S13;
s142, extracting shear stress tau in yz plane yz Shear strain gamma yz ;
S143, will τ yz /γ yz Assigned to variable G yz 。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110698446.6A CN113420481B (en) | 2021-06-23 | 2021-06-23 | Material characteristic calculation method for back field magnet winding |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110698446.6A CN113420481B (en) | 2021-06-23 | 2021-06-23 | Material characteristic calculation method for back field magnet winding |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113420481A CN113420481A (en) | 2021-09-21 |
CN113420481B true CN113420481B (en) | 2023-07-11 |
Family
ID=77716335
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110698446.6A Active CN113420481B (en) | 2021-06-23 | 2021-06-23 | Material characteristic calculation method for back field magnet winding |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113420481B (en) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001104522A (en) * | 1999-10-04 | 2001-04-17 | Mizuno Corp | Design aiding system for golf shaft made of fiber reinforced plastic |
CN105468826A (en) * | 2015-11-18 | 2016-04-06 | 北京玻钢院复合材料有限公司 | Design method of composite material |
CN106650141A (en) * | 2016-12-29 | 2017-05-10 | 北京航空航天大学 | Uncertainty analysis method for predicting performance of periodic material |
CN109583052A (en) * | 2018-11-15 | 2019-04-05 | 西南交通大学 | Fiber reinforced ceramic matrix composites micro-structure Accurate Model and modulus calculation method |
CN110162846A (en) * | 2019-04-30 | 2019-08-23 | 西南交通大学 | The design method of composite rigidity under a kind of elasticity separated rail fastening system iron chair |
CN110502793A (en) * | 2019-07-23 | 2019-11-26 | 南京航空航天大学 | A kind of unidirectional ceramic matric composite off-axis tension testpieces optimum design method |
CN111079334A (en) * | 2019-12-16 | 2020-04-28 | 北京理工大学 | Mesomechanics-based short fiber composite material effective elastic modulus prediction method |
CN112016233A (en) * | 2020-08-31 | 2020-12-01 | 江苏骠马智能工业设计研究有限公司 | Dynamics optimization simulation analysis method for driving mechanism of rail-mounted inspection robot |
CN112541289A (en) * | 2020-12-15 | 2021-03-23 | 西南交通大学 | Sound-vibration equivalent modeling method for rail vehicle body aluminum profile structure |
-
2021
- 2021-06-23 CN CN202110698446.6A patent/CN113420481B/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001104522A (en) * | 1999-10-04 | 2001-04-17 | Mizuno Corp | Design aiding system for golf shaft made of fiber reinforced plastic |
CN105468826A (en) * | 2015-11-18 | 2016-04-06 | 北京玻钢院复合材料有限公司 | Design method of composite material |
CN106650141A (en) * | 2016-12-29 | 2017-05-10 | 北京航空航天大学 | Uncertainty analysis method for predicting performance of periodic material |
CN109583052A (en) * | 2018-11-15 | 2019-04-05 | 西南交通大学 | Fiber reinforced ceramic matrix composites micro-structure Accurate Model and modulus calculation method |
CN110162846A (en) * | 2019-04-30 | 2019-08-23 | 西南交通大学 | The design method of composite rigidity under a kind of elasticity separated rail fastening system iron chair |
CN110502793A (en) * | 2019-07-23 | 2019-11-26 | 南京航空航天大学 | A kind of unidirectional ceramic matric composite off-axis tension testpieces optimum design method |
CN111079334A (en) * | 2019-12-16 | 2020-04-28 | 北京理工大学 | Mesomechanics-based short fiber composite material effective elastic modulus prediction method |
CN112016233A (en) * | 2020-08-31 | 2020-12-01 | 江苏骠马智能工业设计研究有限公司 | Dynamics optimization simulation analysis method for driving mechanism of rail-mounted inspection robot |
CN112541289A (en) * | 2020-12-15 | 2021-03-23 | 西南交通大学 | Sound-vibration equivalent modeling method for rail vehicle body aluminum profile structure |
Non-Patent Citations (3)
Title |
---|
Finite element modelling of turbine generator stator end windings for vibration analysis;Iga Y;《IET Electric Power Applications》;75-81 * |
La2-xNdxCuO4+δ(0.1≤x≤1.2)体系中滞弹性弛豫与相变内耗研究;何庆;《物理学报》;第59卷(第11期);7979-7985 * |
考虑三维波动效应的浮承桩纵向振动理论研究;孟坤;《 中国博士学位论文全文数据库工程科技Ⅱ辑》(第04期);C038-10 * |
Also Published As
Publication number | Publication date |
---|---|
CN113420481A (en) | 2021-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Górski et al. | Computation of mechanical properties of parts manufactured by fused deposition modeling using finite element method | |
CN105528503B (en) | A kind of large-scale component dynamic optimization design method based on STRUCTURE DECOMPOSITION | |
Garzon et al. | Improvements of explicit crack surface representation and update within the generalized finite element method with application to three‐dimensional crack coalescence | |
US4858150A (en) | Shape modeling system | |
CN106960088B (en) | Three-dimensional process model automatic updating method facing process change | |
US20090312960A1 (en) | Characterizing Thermomechanical Properties of an Organic Substrate Using Three-Dimensional Finite Element Analysis | |
CN107680168B (en) | Grid simplifying method based on plane fitting in three-dimensional reconstruction | |
US20220130082A1 (en) | Method for rapid reconstruction of woven composite material microstructure based on topological features | |
CN1991846A (en) | Method for creating three-dimensional finite element mesh | |
US8345955B2 (en) | Characterizing thermomechanical properties of an organic substrate using finite element analysis | |
CN114719962B (en) | Mechanical vibration digital twin model construction method based on machine learning | |
CN111898332B (en) | Frequency domain simulation adaptive frequency point extraction and calculation method for very large scale integrated circuit | |
Brant et al. | A novel electrochemical micro additive manufacturing method of overhanging metal parts without reliance on support structures | |
CN113420481B (en) | Material characteristic calculation method for back field magnet winding | |
CN111027140A (en) | Airplane standard part model rapid reconstruction method based on multi-view point cloud data | |
Persova et al. | Improving the computational efficiency of solving multisource 3-D airborne electromagnetic problems in complex geological media | |
Herbst et al. | Improved transmission line parameter calculation through TCAD process modeling for superconductor integrated circuit interconnects | |
CN116502595A (en) | Random walking-based capacitance extraction method and device, equipment and medium | |
CN113094960B (en) | Method for quickly generating electromagnetic model of multi-core aviation connector based on moment method | |
CN108137128B (en) | For determining the method and system of connecting element manufacture size | |
CN105894575B (en) | The three-dimensional modeling method and device of road | |
CN109411025B (en) | Modeling method of mixed dislocation bifurcation atomic structure | |
CN114299243A (en) | Point cloud feature enhancement method and device based on multi-scale fusion | |
CN116476070B (en) | Method for adjusting scanning measurement path of large-scale barrel part local characteristic robot | |
Cho et al. | Enhancing buckling performance of perforated composite laminates by manipulating fiber direction using a genetic algorithm |
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 |