CN113420481A - Method for calculating material characteristics of back field magnet winding - Google Patents
Method for calculating material characteristics of back field magnet winding Download PDFInfo
- Publication number
- CN113420481A CN113420481A CN202110698446.6A CN202110698446A CN113420481A CN 113420481 A CN113420481 A CN 113420481A CN 202110698446 A CN202110698446 A CN 202110698446A CN 113420481 A CN113420481 A CN 113420481A
- Authority
- CN
- China
- Prior art keywords
- nodes
- conductor
- selecting
- insulator
- group
- 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
-
- 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)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method for calculating material characteristics of a back 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 calculation for the Bomb + modulus of elasticity ExAnd Poisson's ratio PRxy,PRxzThe boundary condition of (1); s04, calculating the elastic modulus E of the conductorxAnd Poisson's ratio PRxy,PRxz(ii) a S05, applying a calculation method for the elastic modulus E of the conductoryAnd Poisson's ratio PRyzThe boundary condition of (1); s06, calculating the elastic modulus of the conductorEyAnd Poisson's ratio PRyz(ii) a S07, applying a calculation method for the elastic modulus E of the conductorzThe boundary condition of (1); s08, calculating the elastic modulus E of the conductorz(ii) a S09, applying a shear modulus G for calculating the conductorxyThe boundary condition of (1); s10, calculating the shear modulus G of the conductorxy(ii) a The method for calculating the material characteristics of the back field magnet winding can quickly evaluate the mechanical properties and optimize the structure of the back field magnet part, and greatly improves the 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 field magnet winding.
Background
The superconducting cable is wound by a back field magnet winding formed by multistage stranded cables and through pipes, and the high-temperature superconducting material has the characteristics of no resistance and high critical current density, so that the superconducting cable becomes a technical breakthrough for power application. But the material property solution of the high-temperature superconducting material is still an unsolved problem.
Disclosure of Invention
In view of the problems mentioned in the background art, it is an object of the present invention 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 purpose of the invention is realized by the following technical scheme:
a method for calculating material characteristics of a back field magnet winding is characterized by comprising the following steps: 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 calculation method for the elastic modulus E of the conductorxAnd Poisson's ratio PRxy,PRxzThe boundary condition of (1);
s04, calculating the elastic modulus E of the conductorxAnd Poisson's ratio PRxy,PRxz;
S05, applying a calculation method for the elastic modulus E of the conductoryAnd Poisson's ratio PRyzThe boundary condition of (1);
s06, calculating the elastic modulus E of the conductoryAnd Poisson's ratio PRyz;
S07, applyingCalculation of the elastic modulus E of the conductorzThe boundary condition of (1);
s08, calculating the elastic modulus E of the conductorz;
S09, applying a shear modulus G for calculating the conductorxyThe boundary condition of (1);
s10, calculating the shear modulus G of the conductorxy;
S11, applying a shear modulus G for calculating the conductorxzThe boundary condition of (1);
s12, calculating the shear modulus G of the conductorxz;
S13, applying a shear modulus G for calculating the conductoryzThe boundary condition of (1);
s14, calculating the shear modulus G of the conductoryz。
Preferably, the S01 includes:
s011, defining variable parameters based on ANSYS parameterized 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;
and S012, creating key points of the conductor according to the dimension parameters of each component of the conductor under a Cartesian coordinate system and a cylindrical coordinate system, and then respectively creating a line model and a plane model of the section of the conductor by a bottom-up modeling method.
Preferably, the S02 includes:
s021, performing Boolean operation on the surface model of the conductor to enable the topological structures of all the sub-surfaces to meet the mapping mesh division requirement;
s022, respectively defining material characteristics such as elastic modulus, shear modulus, Poisson ratio and the like for the conductor armor and the insulating component;
s023, defining a MESH200 unit, and carrying out mapping meshing division on the surface model of the conductor to generate a two-dimensional surface MESH model;
and S024, creating a sweep line perpendicular to the conductor section, and generating a three-dimensional conductor finite element model through sweep meshing.
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, height and length of the insulator, and respectively applying freedom degree 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 the finite element model after the boundary condition is applied;
s042, extracting average deformation of conductor in x, y and z directionsAnd calculating the corresponding line strain SNx,SNy,SNz;
S043, selecting a group of nodes with x coordinate equal to the width of the insulator, and extracting the average stress S of the nodesx;
S044, calculating Sx/SNx,SNy/SNx,SNz/SNxAnd respectively assigned to the variable Ex,PRxy,PRxz。
Preferably, the S05 includes:
s051, deleting the load boundary condition in the S033;
s052, selecting a group of nodes with the y coordinate equal to the height of the insulator, and applying a pressure load of 1 MPa;
the S06 includes:
s061, solving the finite element model after the boundary condition is applied;
s062, extracting average deformation of conductorAnd calculating the corresponding line strain SNy,SNz;
S063, selecting a group of nodes with y coordinate equal to insulator height, and extracting average stress S of the nodesy;
S064, calculating Sy/SNy,SNy/SNzAnd subjecting it toAre respectively given to variable Ey,μyz。
Preferably, the S07 includes:
s071, deleting the load boundary condition in the S051;
s072, selecting a group of nodes with the z coordinate equal to the length of the insulator, and applying a pressure load of 1 MPa;
the S08 includes:
s081, solving the finite element model after the boundary condition is applied;
S083, selecting a group of nodes with the z coordinate equal to the length of the insulator, and extracting the average stress S of the nodesz;
S084, calculating Sz/SNzAnd assigning it to variable Ez。
Preferably, the S09 includes:
s091, deleting the boundary condition;
s092, selecting a set of nodes on the surface of the insulator where x is 0, and applying a displacement constraint UX is 0 and UY is 0;
s093, selectionxA set of nodes having coordinates equal to the width of the insulator and coupling the degrees of freedom of the set of nodes in the direction x;
s094, selecting a set of nodes of the insulator surface with z equal to 0, and applying displacement constraint UZ equal to 0;
s095, selecting a group of nodes with z coordinates equal to the length of the insulator, and coupling the z-direction freedom of the group of nodes;
s096, selecting a node on the surface of the insulator with z being 0, selecting a node with a node coordinate of (0,0,0) from the node group, and assigning the number of the node to the variable na, selecting a node with an x coordinate being equal to the width of the insulator from the node group, and assigning the number of the node to the variable nb, and assigning the number of any of the remaining nodes to the variable nc, and applying constraints to the nodes na, nb, nc;
s097, a set of nodes with x-coordinate equal to the width of the insulator is selected, and a load of 1N is applied to the nodes in the y-direction.
Preferably, the S10 includes:
s101, solving a finite element model applying boundary conditions in the S09;
s102, extracting shear stress tau in xy planexyAnd shear strain gammaxy;
S103, combining tauxy/γxyGiven variable Gxy;
The S11 includes:
s111, deleting boundary conditions;
s112, selecting a group of nodes on the surface of the insulator with x being 0, and applying displacement constraint UX being 0 and UZ being 0;
s113, selecting a group of nodes with x coordinates equal to the width of the insulator, and coupling the x-direction freedom degrees of the group of nodes;
s114, selecting a group of nodes on the surface of the insulator, where y is 0, and applying displacement constraint, UY is 0;
s115, selecting a group of nodes with y coordinates equal to the height of the insulator, and coupling the y-direction freedom degrees of the group of nodes;
s116, selecting a node group of the insulator surface with z being 0, selecting a node with a node coordinate of (0,0,0) from the node group, and assigning the number of the node to the variable na, selecting a node with an x coordinate being equal to the insulator width from the node group, and assigning the number of the node to the variable nb, and assigning the number of any of the remaining nodes to the variable nc, and applying constraints to the nodes na, nb, nc;
s117, a set of nodes with x-coordinate equal to the width of the insulator is selected, and a load of 1N is applied to the nodes in the z-direction.
Preferably, the S12 includes:
s121, solving the finite element model applying the boundary conditions in the S11;
s122, extracting shear stress tau in xy planexzAnd shear strain gammaxz;
S123, converting tauxz/γxzGiven variable Gxz;
The S13 includes:
s131, deleting the boundary condition in the S11;
s132, selecting a group of nodes on the surface of the insulator where y is 0, and applying displacement constraint UY is 0 and UZ is 0;
s133, selecting a group of nodes with x coordinates equal to the width of the insulator, and coupling the x-direction freedom degrees of the group of nodes;
s134, selecting a set of nodes on the insulator surface with x-0 coordinate 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 y-direction freedom degrees of the group of nodes;
s136, selecting a group of nodes on the surface of the insulator with x being 0, selecting a node with a node coordinate of (0,0,0) from the group of nodes, assigning the number of the node to the variable na, selecting a node with a (x, z) coordinate of (0,0) and a y coordinate equal to the height of the insulator from the group of nodes, assigning the number of the node to the variable nb, assigning the number of any of the remaining nodes to the variable nc, and applying constraints to the nodes na, nb, nc;
s137, a set of nodes having a y-coordinate equal to the insulator height is selected, and a load of 1N is applied to the nodes along z.
Preferably, the S14 includes:
s141, solving the finite element model applying the boundary conditions in the S13;
s142, extracting the shear stress tau in yz planeyzAnd shear strain gammayz;
S143, mixing tauyz/γyzGiven variable Gyz。
In summary, the invention mainly has the following beneficial effects:
the method for calculating the material characteristics of the back field magnet winding can quickly evaluate the mechanical properties and optimize the structure of a back field magnet part, and greatly improves the analysis and calculation efficiency of a back field magnet.
Drawings
FIG. 1 is a block flow diagram of the present invention;
FIG. 2 is a cross-sectional dimension block diagram of a back field magnet winding unit cell.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Example 1
Referring to fig. 1 and 2, a material property calculation method 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 calculation method for the elastic modulus E of the conductorxAnd Poisson's ratio PRxy,PRxzThe boundary condition of (1);
s04, calculating the elastic modulus E of the conductorxAnd Poisson's ratio PRxy,PRxz;
S05, applying a calculation method for the elastic modulus E of the conductoryAnd Poisson's ratio PRyzThe boundary condition of (1);
s06, calculating the elastic modulus E of the conductoryAnd Poisson's ratio PRyz;
S07, applying a calculation method for the elastic modulus E of the conductorzThe boundary condition of (1);
s08, calculating the elastic modulus E of the conductorz;
S09, applying a shear modulus G for calculating the conductorxyThe boundary condition of (1);
s10, calculating the shear modulus G of the conductorxy;
S11, applying a shear modulus G for calculating the conductorxzThe boundary condition of (1);
s12, calculating the shear modulus G of the conductorxz;
S13, applying a shear modulus G for calculating the conductoryzThe boundary condition of (1);
s14, calculating the shear modulus G of the conductoryz。
Wherein S01 includes:
s011, defining variable parameters based on ANSYS parameterized 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;
and S012, creating key points of the conductor according to the dimension parameters of each component of the conductor under a Cartesian coordinate system and a cylindrical coordinate system, and then respectively creating a line model and a plane model of the section of the conductor by a bottom-up modeling method.
Wherein S02 includes:
s021, performing Boolean operation on the surface model of the conductor to enable the topological structures of all the sub-surfaces to meet the mapping mesh division requirement;
s022, respectively defining material characteristics such as elastic modulus, shear modulus, Poisson' S ratio and the like for the conductor armor and the insulating member, wherein the material characteristics are defined as shown in a table 1;
table 1: material Properties of the conductor component (temperature: 4K)
S023, defining a MESH200 unit, and carrying out mapping meshing division on the surface model of the conductor to generate a two-dimensional surface MESH model;
and S024, creating a sweep line perpendicular to the conductor section, and generating a three-dimensional conductor finite element model through sweep meshing.
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, height and length of the insulator, and respectively applying freedom degree 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 the finite element model after the boundary condition is applied;
s042, extracting average deformation of conductor in x, y and z directionsAnd calculating the corresponding line strain SNx,SNy,SNz;
S043, selecting a group of nodes with x coordinate equal to the width of the insulator, and extracting the average stress S of the nodesx;
S044, calculating Sx/SNx,SNy/SNx,SNz/SNxAnd respectively assigned to the variable Ex,PRxy,PRxz。
Wherein S05 includes:
s051, deleting the load boundary condition in the S033;
s052, selectingyA set of nodes with coordinates equal to the height of the insulator and applying a pressure load of 1 MPa;
wherein S06 includes:
s061, solving the finite element model after the boundary condition is applied;
s062, extracting average deformation of conductorAnd calculate the corresponding line responseSN-changing devicey,SNz;
S063, selecting a group of nodes with y coordinate equal to insulator height, and extracting average stress S of the nodesy;
S064, calculating Sy/SNy,SNy/SNzAnd respectively assigned to the variable Ey,μyz。
Wherein S07 includes:
s071, deleting the load boundary condition in S051;
s072, selecting a group of nodes with the z coordinate equal to the length of the insulator, and applying a pressure load of 1 MPa;
wherein S08 includes:
s081, solving the finite element model after the boundary condition is applied;
S083, selecting a group of nodes with the z coordinate equal to the length of the insulator, and extracting the average stress S of the nodesz;
S084, calculating Sz/SNzAnd assigning it to variable Ez。
Wherein S09 includes:
s091, deleting the boundary condition;
s092, selecting a set of nodes on the surface of the insulator where x is 0, and applying a displacement constraint UX is 0 and UY is 0;
s093, selecting a group of nodes with x coordinates equal to the width of the insulator, and coupling the direction freedom degrees of the group of nodes x;
s094, selecting a set of nodes of the insulator surface with z equal to 0, and applying displacement constraint UZ equal to 0;
s095, selecting a group of nodes with z coordinates equal to the length of the insulator, and coupling the z-direction freedom of the group of nodes;
s096, selecting a node on the surface of the insulator with z being 0, selecting a node with a node coordinate of (0,0,0) from the node group, and assigning the number of the node to the variable na, selecting a node with an x coordinate being equal to the width of the insulator from the node group, and assigning the number of the node to the variable nb, and assigning the number of any of the remaining nodes to the variable nc, and applying constraints to the nodes na, nb, nc;
s097, a set of nodes with x-coordinate equal to the width of the insulator is selected, and a load of 1N is applied to the nodes in the y-direction.
Wherein S10 includes:
s101, solving a finite element model applying boundary conditions in the S09;
s102, extracting shear stress tau in xy planexyAnd shear strain gammaxy;
S103, combining tauxy/γxyGiven variable Gxy;
Wherein S11 includes:
s111, deleting boundary conditions;
s112, selecting a group of nodes on the surface of the insulator with x being 0, and applying displacement constraint UX being 0 and UZ being 0;
s113, selecting a group of nodes with x coordinates equal to the width of the insulator, and coupling the x-direction freedom degrees of the group of nodes;
s114, selecting a group of nodes on the surface of the insulator, where y is 0, and applying displacement constraint, UY is 0;
s115, selecting a group of nodes with y coordinates equal to the height of the insulator, and coupling the y-direction freedom degrees of the group of nodes;
s116, selecting a node group of the insulator surface with z being 0, selecting a node with a node coordinate of (0,0,0) from the node group, and assigning the number of the node to the variable na, selecting a node with an x coordinate being equal to the insulator width from the node group, and assigning the number of the node to the variable nb, and assigning the number of any of the remaining nodes to the variable nc, and applying constraints to the nodes na, nb, nc;
s117, a set of nodes with x-coordinate equal to the width of the insulator is selected, and a load of 1N is applied to the nodes in the z-direction.
Wherein S12 includes:
s121, solving the finite element model applying the boundary conditions in the S11;
s122, extracting shear stress tau in xy planexzAnd shear strain gammaxz;
S123, converting tauxz/γxzGiven variable Gxz;
Wherein S13 includes:
s131, deleting the boundary conditions in the S11;
s132, selecting a group of nodes on the surface of the insulator where y is 0, and applying displacement constraint UY is 0 and UZ is 0;
s133, selecting a group of nodes with x coordinates equal to the width of the insulator, and coupling the x-direction freedom degrees of the group of nodes;
s134, selecting a set of nodes on the insulator surface with x-0 coordinate 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 y-direction freedom degrees of the group of nodes;
s136, selecting a group of nodes on the surface of the insulator with x being 0, selecting a node with a node coordinate of (0,0,0) from the group of nodes, assigning the number of the node to the variable na, selecting a node with a (x, z) coordinate of (0,0) and a y coordinate equal to the height of the insulator from the group of nodes, assigning the number of the node to the variable nb, assigning the number of any of the remaining nodes to the variable nc, and applying constraints to the nodes na, nb, nc;
s137, selecting a group of nodes with y coordinates equal to the height of the insulator, and applying 1N load to the nodes along z, wherein the applied boundary conditions are shown in Table 2:
TABLE 2 loads and boundary conditions
Wherein S14 includes:
s141, solving the finite element model applying the boundary conditions in the S13;
s142, extracting the shear stress tau in yz planeyzAnd shear strain gammayz;
S143, mixing tauyz/γyzGiven variable GyzThe results of the calculations are shown in table 3:
TABLE 3 characteristics of back field magnet conductor winding materials
Referring to fig. 1 and 2, the method for calculating the material characteristics of the back field magnet winding can quickly perform mechanical property evaluation and structure optimization design of a back field magnet component, and greatly improves the analysis and calculation efficiency of the back field magnet.
Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that changes, modifications, substitutions and alterations can be made in these embodiments without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.
Claims (10)
1. A method for calculating material characteristics of a back field magnet winding is characterized by comprising the following steps: 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 calculation method for the elastic modulus E of the conductorxAnd Poisson's ratio PRxy,PRxzThe boundary condition of (1);
s04, calculating the elastic modulus E of the conductorxAnd Poisson's ratio PRxy,PRxz;
S05, applying a calculation method for the elastic modulus E of the conductoryAnd Poisson's ratio PRyzThe boundary condition of (1);
s06, calculating the elastic modulus E of the conductoryAnd Poisson's ratio PRyz;
S07, applying a calculation method for the elastic modulus E of the conductorzThe boundary condition of (1);
s08, calculating the elastic modulus E of the conductorz;
S09, applying a shear modulus G for calculating the conductorxyThe boundary condition of (1);
s10, calculating the shear modulus G of the conductorxy;
S11, applying a shear modulus G for calculating the conductorxzThe boundary condition of (1);
s12, calculating the shear modulus G of the conductorxz;
S13, applying a shear modulus G for calculating the conductoryzThe boundary condition of (1);
s14, calculating the shear modulus G of the conductoryz。
2. The method of calculating material properties of a back field magnet winding according to claim 1, wherein: the S01 includes:
s011, defining variable parameters based on ANSYS parameterized 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;
and S012, creating key points of the conductor according to the dimension parameters of each component of the conductor under a Cartesian coordinate system and a cylindrical coordinate system, and then respectively creating a line model and a plane model of the section of the conductor by a bottom-up modeling method.
3. The method of calculating material properties of a back field magnet winding according to claim 1, wherein: the S02 includes:
s021, performing Boolean operation on the surface model of the conductor to enable the topological structures of all the sub-surfaces to meet the mapping mesh division requirement;
s022, respectively defining material characteristics such as elastic modulus, shear modulus, Poisson ratio and the like for the conductor armor and the insulating component;
s023, defining a MESH200 unit, and carrying out mapping meshing division on the surface model of the conductor to generate a two-dimensional surface MESH model;
and S024, creating a sweep line perpendicular to the conductor section, and generating a three-dimensional conductor finite element model through sweep meshing.
4. The method of calculating material properties of a back field magnet winding according to claim 1, wherein: 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, height and length of the insulator, and respectively applying freedom degree 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 the finite element model after the boundary condition is applied;
s042, extracting average deformation of conductor in x, y and z directionsAnd calculating the corresponding line strain SNx,SNy,SNz;
S043, selecting a group of nodes with x coordinate equal to the width of the insulator, and extracting the average stress S of the nodesx;
S044, calculating Sx/SNx,SNy/SNx,SNz/SNxAnd respectively assigned to the variable Ex,PRxy,PRxz。
5. The method of calculating material properties of a back field magnet winding according to claim 4, wherein: the S05 includes:
s051, deleting the load boundary condition in the S033;
s052, selecting a group of nodes with the y coordinate equal to the height of the insulator, and applying a pressure load of 1 MPa;
the S06 includes:
s061, solving the finite element model after the boundary condition is applied;
s062, extracting average deformation of conductorAnd calculating the corresponding line strain SNy,SNz;
S063, selecting a group of nodes with y coordinate equal to insulator height, and extracting average stress S of the nodesy;
S064, calculating Sy/SNy,SNy/SNzAnd respectively assigned to the variable Ey,μyz。
6. The method of calculating material properties of a back field magnet winding according to claim 5, wherein: the S07 includes:
s071, deleting the load boundary condition in the S051;
s072, selecting a group of nodes with the z coordinate equal to the length of the insulator, and applying a pressure load of 1 MPa;
the S08 includes:
s081, solving the finite element model after the boundary condition is applied;
S083, selecting a group of nodes with the z coordinate equal to the length of the insulator, and extracting the average stress S of the nodesz;
S084, calculating Sz/SNzAnd assigning it to variable Ez。
7. The method of calculating material properties of a back field magnet winding according to claim 1, wherein: the S09 includes:
s091, deleting the boundary condition;
s092, selecting a set of nodes on the surface of the insulator where x is 0, and applying a displacement constraint UX is 0 and UY is 0;
s093, selecting a group of nodes with x coordinates equal to the width of the insulator, and coupling the direction freedom degrees of the group of nodes x;
s094, selecting a set of nodes of the insulator surface with z equal to 0, and applying displacement constraint UZ equal to 0;
s095, selecting a group of nodes with z coordinates equal to the length of the insulator, and coupling the z-direction freedom of the group of nodes;
s096, selecting a node on the surface of the insulator with z being 0, selecting a node with a node coordinate of (0,0,0) from the node group, and assigning the number of the node to the variable na, selecting a node with an x coordinate being equal to the width of the insulator from the node group, and assigning the number of the node to the variable nb, and assigning the number of any of the remaining nodes to the variable nc, and applying constraints to the nodes na, nb, nc;
s097, a set of nodes with x-coordinate equal to the width of the insulator is selected, and a load of 1N is applied to the nodes in the y-direction.
8. The method of calculating material properties of a back field magnet winding according to claim 7, wherein: the S10 includes:
s101, solving a finite element model applying boundary conditions in the S09;
s102, extracting shear stress tau in xy planexyAnd shear strain gammaxy;
S103, combining tauxy/γxyGiven variable Gxy;
The S11 includes:
s111, deleting boundary conditions;
s112, selecting a group of nodes on the surface of the insulator with x being 0, and applying displacement constraint UX being 0 and UZ being 0;
s113, selecting a group of nodes with x coordinates equal to the width of the insulator, and coupling the x-direction freedom degrees of the group of nodes;
s114, selecting a group of nodes on the surface of the insulator, where y is 0, and applying displacement constraint, UY is 0;
s115, selecting a group of nodes with y coordinates equal to the height of the insulator, and coupling the y-direction freedom degrees of the group of nodes;
s116, selecting a node group of the insulator surface with z being 0, selecting a node with a node coordinate of (0,0,0) from the node group, and assigning the number of the node to the variable na, selecting a node with an x coordinate being equal to the insulator width from the node group, and assigning the number of the node to the variable nb, and assigning the number of any of the remaining nodes to the variable nc, and applying constraints to the nodes na, nb, nc;
s117, a set of nodes with x-coordinate equal to the width of the insulator is selected, and a load of 1N is applied to the nodes in the z-direction.
9. The method of calculating material properties of a back field magnet winding according to claim 8, wherein: the S12 includes:
s121, solving the finite element model applying the boundary conditions in the S11;
s122, extracting shear stress tau in xy planexzAnd shear strain gammaxz;
S123, converting tauxz/γxzGiven variable Gxz;
The S13 includes:
s131, deleting the boundary condition in the S11;
s132, selecting a group of nodes on the surface of the insulator where y is 0, and applying displacement constraint UY is 0 and UZ is 0;
s133, selecting a group of nodes with x coordinates equal to the width of the insulator, and coupling the x-direction freedom degrees of the group of nodes;
s134, selecting a set of nodes on the insulator surface with x-0 coordinate 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 y-direction freedom degrees of the group of nodes;
s136, selecting a group of nodes on the surface of the insulator with x being 0, selecting a node with a node coordinate of (0,0,0) from the group of nodes, assigning the number of the node to the variable na, selecting a node with a (x, z) coordinate of (0,0) and a y coordinate equal to the height of the insulator from the group of nodes, assigning the number of the node to the variable nb, assigning the number of any of the remaining nodes to the variable nc, and applying constraints to the nodes na, nb, nc;
s137, a set of nodes having a y-coordinate equal to the insulator height is selected, and a load of 1N is applied to the nodes along z.
10. The method of calculating material properties of a back field magnet winding according to claim 9, wherein: the S14 includes:
s141, solving the finite element model applying the boundary conditions in the S13;
s142, extracting the shear stress tau in yz planeyzAnd shear strain gammayz;
S143, mixing tauyz/γyzGiven variable Gyz。
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 true CN113420481A (en) | 2021-09-21 |
CN113420481B 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 |
---|
IGA Y: "Finite element modelling of turbine generator stator end windings for vibration analysis", 《IET ELECTRIC POWER APPLICATIONS》 * |
何庆: "La2-xNdxCuO4+δ(0.1≤x≤1.2)体系中滞弹性弛豫与相变内耗研究", 《物理学报》 * |
孟坤: "考虑三维波动效应的浮承桩纵向振动理论研究", 《 中国博士学位论文全文数据库工程科技Ⅱ辑》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113420481B (en) | 2023-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ma et al. | A three-dimensional Monte Carlo model for electrically conductive polymer matrix composites filled with curved fibers | |
Górski et al. | Computation of mechanical properties of parts manufactured by fused deposition modeling using finite element method | |
US8056027B2 (en) | Characterizing thermomechanical properties of an organic substrate using three-dimensional finite element analysis | |
WO2007068674A1 (en) | Two-dimensional orthogonal wire harness representation | |
Garzon et al. | Improvements of explicit crack surface representation and update within the generalized finite element method with application to three‐dimensional crack coalescence | |
CN106960088B (en) | Three-dimensional process model automatic updating method facing process change | |
Gholaminezhad et al. | Multi-objective reliability-based robust design optimization of robot gripper mechanism with probabilistically uncertain parameters | |
US7725846B2 (en) | Conversion of a three-dimensional wire harness representative to a two-dimensional orthogonal wire harness representative | |
CN113420481A (en) | Method for calculating material characteristics of back field magnet winding | |
CN114399591A (en) | Laser point cloud data-based vectorization extraction method and device for power transmission line body | |
CN104268321B (en) | Satellite cable lightweight modeling method | |
CN107368563A (en) | The delet method and device of database data, electronic equipment, storage medium | |
CN114581871A (en) | Target object detection method and device | |
Zhao et al. | A divide-and-conquer linear scaling three-dimensional fragment method for large scale electronic structure calculations | |
CN113792420A (en) | Method for generating microstructure of unidirectional random fiber composite material by considering interface | |
Shen et al. | On embedding between 2D meshes of the same size | |
Arbelaez et al. | Cable deformation simulation and a hierarchical framework for Nb3Sn Rutherford cables | |
CN112906132A (en) | Method and device for generating aircraft harness component data | |
Mukhopadhyay et al. | Modified Hamiltonian chain: a graph theoretic approach to group technology | |
CN116476070A (en) | Method for adjusting scanning measurement path of large-scale barrel part local characteristic robot | |
Barzi et al. | A Model to Study Plastic Deformation in RRP ${\rm Nb} _ {3}{\rm Sn} $ Wires | |
Vogiatzis et al. | Multi-material topology optimization of structures infilled with conformal metamaterials | |
CN112950706B (en) | Mobile equipment positioning data processing method, device, equipment and storage medium | |
CN108197398A (en) | A kind of finite element method of D braided composites failure predicted based on space group P4 | |
CN109411025B (en) | Modeling method of mixed dislocation bifurcation atomic structure |
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 |