CN113420481A - Method for calculating material characteristics of back field magnet winding - Google Patents

Method for calculating material characteristics of back field magnet winding Download PDF

Info

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
Application number
CN202110698446.6A
Other languages
Chinese (zh)
Other versions
CN113420481B (en
Inventor
汪献伟
何庆
李秀莲
蒋琴仙
谢文涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangsu University of Technology
Original Assignee
Jiangsu University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jiangsu University of Technology filed Critical Jiangsu University of Technology
Priority to CN202110698446.6A priority Critical patent/CN113420481B/en
Publication of CN113420481A publication Critical patent/CN113420481A/en
Application granted granted Critical
Publication of CN113420481B publication Critical patent/CN113420481B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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 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

Method for calculating material characteristics of back field magnet winding
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 directions
Figure BDA0003129465520000031
And 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 conductor
Figure BDA0003129465520000032
And 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;
s082, extracting average deformation of conductor
Figure BDA0003129465520000041
And calculating the corresponding line strain SNz；
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)
Figure BDA0003129465520000091
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 directions
Figure BDA0003129465520000092
And 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 conductor
Figure BDA0003129465520000101
And 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;
s082, extracting average deformation of conductor
Figure BDA0003129465520000102
And calculating the corresponding line strain SNz；
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
Figure BDA0003129465520000131
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
Figure BDA0003129465520000141
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 directions
Figure FDA0003129465510000021
And 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 conductor
Figure FDA0003129465510000031
And 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;
s082, extracting average deformation of conductor
Figure FDA0003129465510000032
And calculating the corresponding line strain SNz；
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。
CN202110698446.6A 2021-06-23 2021-06-23 Material characteristic calculation method for back field magnet winding Active CN113420481B (en)

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)

* Cited by examiner, † Cited by third party
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

Patent Citations (9)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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
Górski et al. Computation of mechanical properties of parts manufactured by fused deposition modeling using finite element method
EP1960234A1 (en) Two-dimensional orthogonal wire harness representation
Chai et al. Dense matrix inversion of linear complexity for integral-equation-based large-scale 3-D capacitance extraction
US20090312960A1 (en) Characterizing Thermomechanical Properties of an Organic Substrate Using Three-Dimensional Finite Element Analysis
US20220130082A1 (en) Method for rapid reconstruction of woven composite material microstructure based on topological features
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
Frei et al. Complex approaches for the calculation of EMC problems of large systems
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
CN113792420A (en) Method for generating microstructure of unidirectional random fiber composite material by considering interface
Arbelaez et al. Cable deformation simulation and a hierarchical framework for Nb3Sn Rutherford cables
Mukhopadhyay et al. Modified Hamiltonian chain: a graph theoretic approach to group technology
CN116502595A (en) Random walking-based capacitance extraction method and device, equipment and medium
Barzi et al. A Model to Study Plastic Deformation in RRP ${\rm Nb} _ {3}{\rm Sn} $ Wires
CN116228158A (en) Big data platform plan management method and related device based on constructional engineering
CN108197398A (en) A kind of finite element method of D braided composites failure predicted based on space group P4
CN113778024A (en) Method and device for scheduling equipment, electronic equipment and storage medium
Pan et al. A fast multipole method for embedded structure in a stratified medium
CN104809280A (en) Cable interference detection method based on geometrical characteristics
Han et al. A novel field-line-circuit hybrid algorithm for transient responses prediction of transmission lines based on FDTD method
Kocer An automated incremental finite element study of Hertzian cone crack growth
Harten Method of artificial compression. I. shocks and contact discontinuities
CN115687492B (en) Automobile intelligent service platform and method based on multi-party service connection

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