CN108197360A - The automatic dividing system of turbine rotor grid and method - Google Patents

The automatic dividing system of turbine rotor grid and method Download PDF

Info

Publication number
CN108197360A
CN108197360A CN201711405395.3A CN201711405395A CN108197360A CN 108197360 A CN108197360 A CN 108197360A CN 201711405395 A CN201711405395 A CN 201711405395A CN 108197360 A CN108197360 A CN 108197360A
Authority
CN
China
Prior art keywords
seed
characteristic
aspect ratio
size
edge
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
CN201711405395.3A
Other languages
Chinese (zh)
Other versions
CN108197360B (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.)
Harbin Turbine Co Ltd
Hadian Power Equipment National Engineering Research Center Co Ltd
Original Assignee
Harbin Turbine Co Ltd
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 Harbin Turbine Co Ltd filed Critical Harbin Turbine Co Ltd
Priority to CN201711405395.3A priority Critical patent/CN108197360B/en
Publication of CN108197360A publication Critical patent/CN108197360A/en
Application granted granted Critical
Publication of CN108197360B publication Critical patent/CN108197360B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Image Processing (AREA)

Abstract

The automatic dividing system of turbine rotor grid and method, are related to mesh generation service system and method.In order to solve can not to carry out ensureing the mesh generation of finite element simulation precision for different size turbine rotors at present.The present invention includes collecting the characteristic information rotator model feature acquisition module for importing rotator model；The Eigenvalues analysis module numbered and classified according to characteristic coordinates to characteristic edge for each characteristic edge and characteristic node；Most suitable fabric mode will be found in preposition feature when postposition aspect ratio input logic library, and curvature control is carried out to seed and is responsible for the characteristic edge fabric module of characteristic edge fabric；To preposition aspect ratio R1And postposition aspect ratio R2Carry out the logical base of logic judgment；It is characterized model partition and divides hexahedral mesh and the mesh generation module limited the geometric eccentricity factor；The node number that stroke seed units are applied is stored to the gridding information memory module in cell data library.The present invention divides automatically suitable for turbine rotor grid.

Description

The automatic dividing system of turbine rotor grid and method
Technical field
The present invention relates to mesh generation service system and methods, and in particular to the automatic dividing system of turbine rotor grid and Division methods.
Background technology
Finite Element has become the important means of Steam Turbine turbine rotor Strength co-mputation now, and mesh generation is done Restrictive condition for FEM calculation precision is also increasingly taken seriously, and to improve rotor strength computational accuracy, ensures steam turbine Operational safety, it is necessary to divide suitable grid.Turbine rotor mesh generation is still in the manual division stage at home at present, turns Sub- automatic grid dividing system still in blank stage, main reasons is that：The complicated and form of turbine rotor is not Unified, functional structure is more, can not establish the mesh generation rule of uniformity；The large scale of turbine rotor and small size phase More than poor an order of magnitude, the fabric rule of uniformity can not be formed；Turbine blade root groove and balancing drum root are small sizes Intensively, while it is also the restriction site of position, i.e. intensity that stress concentration often occurs, to ensure that precision need to be to carrying out net herein Lattice encrypted work, but number of grid can be caused excessive, occupancy computing resource is more, leads to the lance of computational accuracy and computing resource Shield；Since the mesh generation thinking for calculating analysis personnel is different, lack relevant industries standard, accuracy and computational efficiency are irregular not Together, it is difficult to realize and improve working efficiency and widespread adoption.
Invention content
The present invention in order to solve at present can not in the case of less calculation amount for different size turbine rotors carry out Ensure the mesh generation of finite element simulation precision.And then propose a kind of automatic dividing system of turbine rotor grid and method.
The automatic dividing system of turbine rotor grid, including：
Rotator model feature acquisition module collects the characteristic information for importing rotator model, is sat including each endpoint of rotator model The characteristic informations such as mark, each feature edge lengths, angle；
Eigenvalues analysis module divides characteristic edge for each characteristic edge and characteristic node number and according to characteristic coordinates Class, while obtain the length information of every characteristic edge, and by the length n of every characteristic edge respectively with preposition feature edge lengths n2、 Postposition feature edge lengths n3It does and compares, respectively obtain preposition aspect ratioPostposition aspect ratio
Characteristic edge fabric module, by preposition aspect ratio R1And postposition aspect ratio R2Most suitable cloth is found in input logic library Kind mode, and curvature control is carried out to seed, curvature control is carried out by the maximum deviation factor and the minimum dimension factor, limitation is most Big discrepancy factor；Until the seed of cloth meet it is above-mentioned required, complete the work of characteristic edge fabric；
Logical base is to preposition aspect ratio R1And postposition aspect ratio R2Carry out logic judgment；Logic it is as follows：
To preposition aspect ratio R1Judged,
If R1>=5, seed is deviated to common node direction, and grid full-size isMinimum dimension is
If 3≤R1＜ 5, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If 1≤R1＜ 3, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If R1＜ 1, seed are set as without offset, size of mesh opening
To postposition aspect ratio R2Judged,
If R2>=5, seed is deviated to common node direction, and grid full-size isMinimum dimension is
If 3≤R2＜ 5, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If 1≤R2＜ 3, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If R2＜ 1, seed are set as without offset, size of mesh opening
Mesh generation module according to the seed position arranged in characteristic edge, model is characterized by way of scanning and is drawn Graduation divides hexahedral mesh, and the geometric eccentricity factor is limited；And rotor physical model is separated into based on mesh generation Limited a unit；
Each node coordinate is stored ready-portioned seed to number of nodes according to node number sequence by gridding information memory module According in library, ready-portioned unit is stored according to the node number that unit number sequence is applied unit into cell data library, it will Node database is integrated with cell data library, forms the inp files that general finite meta software can be directly read.
Turbine rotor grid Automated Partition Method, includes the following steps：
Step 1 establishes rotator model using modeling software；
Step 2, tagsort and Eigenvalues analysis：
Step 2.1 collects the characteristic information for importing rotator model, including each extreme coordinates of rotator model, each feature length of side The characteristic informations such as degree, angle；
Step 2.2 is numbered for each characteristic edge and characteristic node, and characteristic edge refers to the side information of rotator model, sat according to feature Mark classifies to characteristic edge, if a plurality of characteristic edge shares a common node, is divided into same class, is that every characteristic edge is distributed Preposition feature and postposition feature, there are two characteristic nodes for a characteristic edge, and preposition node is defined as with after according to number size Put node, except characteristic edge in itself in addition to other sides for sharing preposition nodes be known as preposition feature, except characteristic edge is outer in itself, other are shared The side of postposition node is known as postposition feature；Obtain the length information of every characteristic edge, and by the length n of every characteristic edge respectively with Preposition feature edge lengths n2, postposition feature edge lengths n3It does and compares, respectively obtain preposition aspect ratioPostposition aspect ratio
Step 2.3, by preposition aspect ratio R1And postposition aspect ratio R2Most suitable fabric mode is found in input logic library, I.e. logically library requires tentatively to arrange seed；The logic of logical base is as follows：
To preposition aspect ratio R1Judged,
If R1>=5, seed is deviated to common node direction, and grid full-size isMinimum dimension is
If 3≤R1＜ 5, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If 1≤R1＜ 3, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If R1＜ 1, seed are set as without offset, size of mesh opening
To postposition aspect ratio R2Judged,
If R2>=5, seed is deviated to common node direction, and grid full-size isMinimum dimension is
If 3≤R2＜ 5, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If 1≤R2＜ 3, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If R2＜ 1, seed are set as without offset, size of mesh opening
And curvature control is carried out to seed, curvature control, limitation are carried out by the maximum deviation factor and the minimum dimension factor The maximum deviation factor and the minimum dimension factor；
Until the seed of cloth meet it is above-mentioned required, complete the work of characteristic edge fabric；
Step 3, characteristic model mesh generation：
According to the seed position arranged in characteristic edge, model partition is characterized by way of scanning and divides hexahedron net Lattice, and the geometric eccentricity factor is limited；
If mesh generation fails, return to step 2 re-starts tagsort and Eigenvalues analysis, adjusts fabric sequence, directly It is completed to mesh generation and inerrancy unit exists, if mesh generation work can not still be completed by recycling ten times, reexamine feature Model re-establishes and imports rotator model；
Rotor physical model is separated by limited a unit based on mesh generation；
Each node coordinate is stored ready-portioned seed into node database according to node number sequence by step 4, will be drawn In the node number storage to cell data library that the unit divided is applied unit according to unit number sequence, by node database with Cell data library is integrated, and forms the file that general finite meta software can be directly read.
Further, the limitation maximum deviation factor described in step 2 is 0.1.
Further, the minimum dimension factor described in step 2 is set as 0.1.
Further, the limitation geometric eccentricity factor described in step 3 is less than 0.2.
The present invention has the following effects that：
For different types of turbine rotor, the present invention believes automatically by the feature of the turbine rotor model to input Breath is identified, and characteristic information is screened and classified, and effect and the pass of each characteristic edge are determined by fuzzy diagnosis mode System, and selects best fabric mode in fabric typelib, by the rational hexahedral element of the discrete forming shape of rotor, so as to fulfill The automatic partition functionality of turbine rotor grid, while can realize the work(of the automatic grid division of different types of turbine rotor Can, it almost can realize mesh generation for more than 90% turbine rotor model, net will not be caused by the present invention Lattice quantity is excessive, it is more to occupy computing resource, and has good computational accuracy for various sizes of turbine rotor.As long as In the case of the model of foundation is rational, by the present invention to the mesh generation of turbine rotor, output result can ensure limited First simulation result precision can calculate for steam turbine limited strength member and provide technical support, and the rotor for meeting Steam Turbine is strong Spend design requirement.
The system of the present invention reserves a large amount of interfaces convenient for secondary development simultaneously, using the present invention by a continuous entity from Limited a unit and node are dissipated into, Strength co-mputation can be carried out with Finite Element for turbine rotor, technical support is provided.
Description of the drawings
Fig. 1 is present system mesh generation effect total figure；
Fig. 2 is present system to intensity key point mesh generation effect detail view；
Fig. 3 is logic judgment flow chart；
Fig. 4 is Grid Edge and the departure degree schematic diagram on realistic model side.
Specific embodiment
Present invention invention establishes a set of automatic dividing system of turbine rotor grid, to the rotator model in import system Carry out auto grid plot work, and reserve a large amount of second development interfaces, be reached for turbine rotor with Finite Element into Row intensity and dynamic analysis provide the purpose of technical preparation, and engineering foundation is provided for turbine rotor safety accounting.
Specific embodiment one：
The automatic dividing system of turbine rotor grid, including：
Rotator model feature acquisition module, it is automatic to collect the characteristic information for importing rotator model, it is respectively held including rotator model The characteristic informations such as point coordinates, each feature edge lengths, angle；
Eigenvalues analysis module divides characteristic edge for each characteristic edge and characteristic node number and according to characteristic coordinates Class, while obtain the length information of every characteristic edge, and by the length n of every characteristic edge respectively with preposition feature edge lengths n2、 Postposition feature edge lengths n3It does and compares, respectively obtain preposition aspect ratioPostposition aspect ratio
Characteristic edge fabric module, by preposition aspect ratio R1And postposition aspect ratio R2Most suitable cloth is found in input logic library Kind mode, and curvature control is carried out to seed, curvature control is carried out by the maximum deviation factor and the minimum dimension factor, limitation is most Big discrepancy factor；Until the seed of cloth meet it is above-mentioned required, complete the work of characteristic edge fabric；
Logical base is to preposition aspect ratio R1And postposition aspect ratio R2Carry out logic judgment；Logic it is as follows：
To preposition aspect ratio R1Judged,
If R1>=5, seed is deviated to common node direction, and grid full-size isMinimum dimension is
If 3≤R1＜ 5, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If 1≤R1＜ 3, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If R1＜ 1, seed are set as without offset, size of mesh opening
To postposition aspect ratio R2Judged,
If R2>=5, seed is deviated to common node direction, and grid full-size isMinimum dimension is
If 3≤R2＜ 5, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If 1≤R2＜ 3, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If R2＜ 1, seed are set as without offset, size of mesh opening
Mesh generation module according to the seed position arranged in characteristic edge, model is characterized by way of scanning and is drawn Graduation divides hexahedral mesh, and the geometric eccentricity factor is limited；And rotor physical model is separated into based on mesh generation Limited a unit；
Each node coordinate is stored ready-portioned seed to number of nodes according to node number sequence by gridding information memory module According in library, ready-portioned unit is stored according to the node number that unit number sequence is applied unit into cell data library, it will Node database is integrated automatically with cell data library, forms the inp files that general finite meta software can be directly read.
For different types of turbine rotor, the present invention believes automatically by the feature of the turbine rotor model to input Breath is identified, and characteristic information is screened and classified, and effect and the pass of each characteristic edge are determined by fuzzy diagnosis mode System, and selects best fabric mode in fabric typelib, by the rational hexahedral element of the discrete forming shape of rotor, so as to fulfill The automatic partition functionality of turbine rotor grid, while can realize the work(of the automatic grid division of different types of turbine rotor Can, almost it can realize mesh generation for more than 90% turbine rotor model.As long as the model in foundation is reasonable In the case of, by the present invention to the mesh generation of turbine rotor, output result can ensure FEM Numerical Simulation precision, Ensure the intensity of key point simultaneously.Present system mesh generation effect total figure is as shown in Figure 1；Present system is to intensity key It is as shown in Figure 2 that dot grid divides effect detail view.
The system of the present invention reserves a large amount of interfaces convenient for secondary development simultaneously, using the present invention by a continuous entity from Limited a unit and node are dissipated into, Strength co-mputation can be carried out with Finite Element for turbine rotor, technical support is provided, Meet the rotor strength design requirement of Steam Turbine.
Specific embodiment two：
Turbine rotor grid Automated Partition Method, includes the following steps：
Step 1 establishes rotator model using modeling softwares such as UG；
Rotator model of different nature can be established according to rotor configurations, for T-shaped blade root and double-T shaped blade root, two can be established Axisymmetric model is tieed up, and for palm type blade root, three-dimensional circular symmetry model can be established；Established rotator model is stored as Stp forms, and be directed into the automatic dividing system of turbine rotor grid；
Step 2, tagsort and Eigenvalues analysis：
Step 2.1, the automatic characteristic information for collecting importing rotator model, including each extreme coordinates of rotator model, each feature The characteristic informations such as edge lengths, angle；
Step 2.2 is numbered for each characteristic edge and characteristic node, and the side information that characteristic edge refers to rotator model (including length, is sat Mark and angle information), classify according to characteristic coordinates to characteristic edge, if a plurality of characteristic edge shares a common node, divide It is that every characteristic edge distributes preposition feature and postposition feature, a characteristic edge is there are two characteristic nodes, according to volume for same class Number size is defined as preposition node and postposition node, except characteristic edge in itself in addition to other sides for sharing preposition nodes be known as preposition spy Sign, except sides of characteristic edge other outer shared postposition nodes in itself are known as postposition feature；The length information of every characteristic edge is obtained, and By the length n of every characteristic edge respectively with preposition feature edge lengths n2, postposition feature edge lengths n3It does and compares, respectively obtain preposition spy Sign ratioPostposition aspect ratioA practical common node can have up to three sides, be directed to preposition node and Speech just has a characteristic edge and two preposition characteristic edges, is required for calculating preposition aspect ratio for every preposition characteristic edge at this time, And find fabric mode；Postposition node and postposition characteristic edge are handled using same way；
Step 2.3, by preposition aspect ratio R1And postposition aspect ratio R2Most suitable fabric mode is found in input logic library, I.e. logically library requires tentatively to arrange seed；As shown in figure 3, the logic of logical base is as follows：
To preposition aspect ratio R1Judged,
If R1>=5, seed is deviated to common node direction, and grid full-size isMinimum dimension is
If 3≤R1＜ 5, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If 1≤R1＜ 3, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If R1＜ 1, seed are set as without offset, size of mesh opening
To postposition aspect ratio R2Judged,
If R2>=5, seed is deviated to common node direction, and grid full-size isMinimum dimension is
If 3≤R2＜ 5, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If 1≤R2＜ 3, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If R2＜ 1, seed are set as without offset, size of mesh opening
Seed is finite element term, refers to the dummy node established for grid division arranged in characteristic edge；
And curvature control is carried out to seed, curvature control is finite element term, passes through the maximum deviation factor and minimum dimension The factor carries out curvature control, limits the maximum deviation factor and the minimum dimension factor (minimum feature resolution)；Grid Edge and reality The departure degree on border model side, as shown in figure 4, wherein h is string deviation, L is chord length, and h/L is discrepancy factor；
Until the seed of cloth meet it is above-mentioned required, complete the work of characteristic edge fabric；
Step 3, characteristic model mesh generation：
According to the seed position arranged in characteristic edge, model partition is characterized by way of scanning and divides hexahedron net Lattice, and the geometric eccentricity factor is limited；
Scan as mechanical term, represent to generate quadrilateral surface grid first on rotor section, by this surface grids along Rotor rotates in a circumferential direction, you can forms hexahedral mesh；The geometric eccentricity factor be finite element term, represent the element sides marked off with The deviation situation on model side；
If mesh generation fails, failure cause is analyzed, return to step 2 re-starts tagsort and characteristic value point Analysis adjusts fabric sequence, until mesh generation is completed and inerrancy unit exists, if mesh generation can not still be completed by recycling ten times Work, then prompt " please reexamine rotator model " dialog box, and reexamine characteristic model, re-establish and import rotor mould Type；
Rotor physical model is separated by limited a unit based on mesh generation；
Each node coordinate is stored ready-portioned seed into node database according to node number sequence by step 4, will be drawn In the node number storage to cell data library that the unit divided is applied unit according to unit number sequence, by node database with Cell data library is integrated automatically, forms the inp files that general finite meta software can be directly read.
Specific embodiment three：
The limitation maximum deviation factor described in present embodiment step 2 is 0.1.
Other steps and parameter are identical with embodiment two.
Specific embodiment four：
In present embodiment step 2, prevent from causing inessential feature grid division sharp unit occur, the minimum Size factor is set as 0.1.
Other steps and parameter are identical with specific embodiment two or three.
Specific embodiment five：
The limitation geometric eccentricity factor described in present embodiment step 3 is less than 0.2.
Other steps and one of parameter and specific embodiment two to four are identical.

Claims (5)

1. the automatic dividing system of turbine rotor grid, which is characterized in that including：
Rotator model feature acquisition module collects the characteristic information for importing rotator model, including each extreme coordinates of rotator model, respectively The characteristic informations such as feature edge lengths, angle；
Eigenvalues analysis module classifies to characteristic edge for each characteristic edge and characteristic node number and according to characteristic coordinates, together When obtain the length information of every characteristic edge, and by the length n of every characteristic edge respectively with preposition feature edge lengths n2, postposition it is special Levy edge lengths n3It does and compares, respectively obtain preposition aspect ratioPostposition aspect ratio
Characteristic edge fabric module, by preposition aspect ratio R1And postposition aspect ratio R2Most suitable fabric side is found in input logic library Formula, and curvature control is carried out to seed, curvature control is carried out by the maximum deviation factor and the minimum dimension factor, is limited maximum inclined From the factor；Until the seed of cloth meet it is above-mentioned required, complete the work of characteristic edge fabric；
Logical base is to preposition aspect ratio R1And postposition aspect ratio R2Carry out logic judgment；Logic it is as follows：
To preposition aspect ratio R1Judged,
If R1>=5, seed is deviated to common node direction, and grid full-size isMinimum dimension is
If 3≤R1＜ 5, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If 1≤R1＜ 3, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If R1＜ 1, seed are set as without offset, size of mesh opening
To postposition aspect ratio R2Judged,
If R2>=5, seed is deviated to common node direction, and grid full-size isMinimum dimension is
If 3≤R2＜ 5, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If 1≤R2＜ 3, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If R2＜ 1, seed are set as without offset, size of mesh opening
Mesh generation module according to the seed position arranged in characteristic edge, model partition is characterized by way of scanning and is drawn Divide hexahedral mesh, and the geometric eccentricity factor is limited；And rotor physical model is separated into based on mesh generation limited A unit；
Each node coordinate is stored ready-portioned seed to node database according to node number sequence by gridding information memory module In, ready-portioned unit is stored according to the node number that unit number sequence is applied unit into cell data library, by node Database is integrated with cell data library, forms the inp files that general finite meta software can be directly read.
2. turbine rotor grid Automated Partition Method, which is characterized in that include the following steps：
Step 1 establishes rotator model using modeling software；
Step 2, tagsort and Eigenvalues analysis：
Step 2.1 collects the characteristic information for importing rotator model, including each extreme coordinates of rotator model, each feature edge lengths, angle The characteristic informations such as degree；
Step 2.2 is numbered for each characteristic edge and characteristic node, and characteristic edge refers to the side information of rotator model, according to characteristic coordinates pair Characteristic edge is classified, if a plurality of characteristic edge shares a common node, is divided into same class, is that every characteristic edge distribution is preposition Feature and postposition feature, there are two characteristic nodes for a characteristic edge, and preposition node and postposition section are defined as according to number size Point, except characteristic edge in itself in addition to other sides for sharing preposition nodes be known as preposition feature, except characteristic edge other outer shared postpositions in itself The side of node is known as postposition feature；Obtain the length information of every characteristic edge, and by the length n of every characteristic edge respectively with it is preposition Feature edge lengths n2, postposition feature edge lengths n3It does and compares, respectively obtain preposition aspect ratioPostposition aspect ratio
Step 2.3, by preposition aspect ratio R1And postposition aspect ratio R2Most suitable fabric mode is found in input logic library, that is, is pressed It requires tentatively to arrange seed according to logical base；The logic of logical base is as follows：
To preposition aspect ratio R1Judged,
If R1>=5, seed is deviated to common node direction, and grid full-size isMinimum dimension is
If 3≤R1＜ 5, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If 1≤R1＜ 3, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If R1＜ 1, seed are set as without offset, size of mesh opening
To postposition aspect ratio R2Judged,
If R2>=5, seed is deviated to common node direction, and grid full-size isMinimum dimension is
If 3≤R2＜ 5, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If 1≤R2＜ 3, seed are deviated to common node direction, and grid full-size isMinimum dimension is
If R2＜ 1, seed are set as without offset, size of mesh opening
And curvature control is carried out to seed, curvature control is carried out by the maximum deviation factor and the minimum dimension factor, limitation is maximum Discrepancy factor and the minimum dimension factor；
Until the seed of cloth meet it is above-mentioned required, complete the work of characteristic edge fabric；
Step 3, characteristic model mesh generation：
According to the seed position arranged in characteristic edge, model partition is characterized by way of scanning and divides hexahedral mesh, And the geometric eccentricity factor is limited；
If mesh generation fails, return to step 2 re-starts tagsort and Eigenvalues analysis, adjustment fabric sequence, until net Lattice divide completion and inerrancy unit exists, if mesh generation work can not still be completed by recycling ten times, reexamine characteristic model, It re-establishes and imports rotator model；
Rotor physical model is separated by limited a unit based on mesh generation；
Each node coordinate is stored ready-portioned seed into node database according to node number sequence by step 4, will be divided The unit node number storage to cell data library of being applied unit according to unit number sequence in, by node database and unit Database is integrated, and forms the file that general finite meta software can be directly read.
3. turbine rotor grid Automated Partition Method according to claim 2, which is characterized in that described in step 2 It is 0.1 to limit the maximum deviation factor.
4. turbine rotor grid Automated Partition Method according to claim 3, which is characterized in that described in step 2 The minimum dimension factor is set as 0.1.
5. the turbine rotor grid Automated Partition Method according to one of claim 2 to 4, which is characterized in that in step 3 The limitation geometric eccentricity factor is less than 0.2.
CN201711405395.3A 2017-12-22 2017-12-22 Automatic division system and method for steam turbine rotor grids Active CN108197360B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711405395.3A CN108197360B (en) 2017-12-22 2017-12-22 Automatic division system and method for steam turbine rotor grids

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711405395.3A CN108197360B (en) 2017-12-22 2017-12-22 Automatic division system and method for steam turbine rotor grids

Publications (2)

Publication Number Publication Date
CN108197360A true CN108197360A (en) 2018-06-22
CN108197360B CN108197360B (en) 2021-07-13

Family

ID=62583279

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711405395.3A Active CN108197360B (en) 2017-12-22 2017-12-22 Automatic division system and method for steam turbine rotor grids

Country Status (1)

Country Link
CN (1) CN108197360B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112487575A (en) * 2020-11-25 2021-03-12 南京航空航天大学 Method for inspecting connection node of high-pressure rotor and low-pressure rotor of double-rotor system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5844564A (en) * 1996-04-12 1998-12-01 Institute Francais Du Petrole Method for generating a 3D-grid pattern matching the geometry of a body in order to achieve a model representative of this body
CN102306396A (en) * 2011-09-15 2012-01-04 山东大学 Three-dimensional entity model surface finite element mesh automatic generation method
CN103236087A (en) * 2013-04-25 2013-08-07 河海大学 Triangular prism-shaped geological model construction method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5844564A (en) * 1996-04-12 1998-12-01 Institute Francais Du Petrole Method for generating a 3D-grid pattern matching the geometry of a body in order to achieve a model representative of this body
CN102306396A (en) * 2011-09-15 2012-01-04 山东大学 Three-dimensional entity model surface finite element mesh automatic generation method
CN103236087A (en) * 2013-04-25 2013-08-07 河海大学 Triangular prism-shaped geological model construction method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
KULESZA Z 等: "Rigid finite element model of a cracked rotor", 《JOURNAL OF SOUND & VIBRATION》 *
王璋奇 等: "基于参数化建模的转子有限元剖分", 《热能动力工程》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112487575A (en) * 2020-11-25 2021-03-12 南京航空航天大学 Method for inspecting connection node of high-pressure rotor and low-pressure rotor of double-rotor system
CN112487575B (en) * 2020-11-25 2023-06-06 南京航空航天大学 Method for checking connection node of high-voltage rotor and low-voltage rotor of double-rotor system

Also Published As

Publication number Publication date
CN108197360B (en) 2021-07-13

Similar Documents

Publication Publication Date Title
CN104331621B (en) A kind of wind-resources computational methods
CN102306396B (en) Three-dimensional entity model surface finite element mesh automatic generation method
CN103268572B (en) A kind of microcosmic structure method of ten million multikilowatt large-scale wind electricity base wind measurement network
CN101257149B (en) Method for dividing aerial reflecting plane graticule based on structure electromagnetic coupling
CN103871102B (en) A kind of road three-dimensional fine modeling method based on elevational point and road profile face
CN107045526A (en) A kind of pattern recognition method of electronics architectural working drawing
CN110134752B (en) Three-dimensional large-scene modeling data processing method and device
CN110030972A (en) Tunnel Overbreak & Underbreak detection method based on ExcelVBA
JP4322615B2 (en) Component information processing apparatus, component information processing method, and component information processing program provided with three-dimensional layout adjustment CAD system
CN103745459A (en) Detection method of an unstructured point cloud feature point and extraction method thereof
CN103838907A (en) Curved surface cutting trajectory obtaining method based on STL model
CN107944203A (en) A kind of visual architectural design method of wind-velocity streamline
KR101190882B1 (en) Apparatus and method of cloud computing for online visualization of noise map at high speed
CN107464285A (en) The Meshing Method and device of a kind of threedimensional model
CN115758938B (en) Adhesive boundary flow field numerical simulation-oriented boundary layer grid generation method
CN109711680A (en) River System Hydropower planing method based on GIS digitizing technique
CN105894553B (en) A kind of Street Space form layout method based on grid selection
CN102982567B (en) The deformable body collision detection elimination method that a kind of Corpus--based Method is analyzed
CN104537254B (en) A kind of drafting method that becomes more meticulous based on social statistics data
CN108197360A (en) The automatic dividing system of turbine rotor grid and method
CN107560539A (en) For the method and system measured in complex pattern structure
CN104732589B (en) Quick hybrid grid generation method
CN110334450A (en) A kind of multi-blocked structure grid generate in object plane erroneous projection restorative procedure
CN113987610A (en) Matching method of different-resolution clothing curved surface meshes based on mesh mapping
CN107276807B (en) Hierarchical network community tree pruning method based on community dynamic compactness

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
TR01 Transfer of patent right

Effective date of registration: 20221226

Address after: 150000 building 3, high tech production base, Nangang District, Harbin City, Heilongjiang Province

Patentee after: HARBIN TURBINE Co.,Ltd.

Patentee after: HADIAN POWER EQUIPMENT NATIONAL ENGINEERING RESEARCH CENTER CO.,LTD.

Address before: 150046 No. three power road 345, Xiangfang District, Heilongjiang, Harbin

Patentee before: HARBIN TURBINE Co.,Ltd.

TR01 Transfer of patent right