CN105787226B - The parameterized model of four side finite element mesh models is rebuild - Google Patents
The parameterized model of four side finite element mesh models is rebuild Download PDFInfo
- Publication number
- CN105787226B CN105787226B CN201610307722.0A CN201610307722A CN105787226B CN 105787226 B CN105787226 B CN 105787226B CN 201610307722 A CN201610307722 A CN 201610307722A CN 105787226 B CN105787226 B CN 105787226B
- Authority
- CN
- China
- Prior art keywords
- quadrangle
- model
- boundary
- point
- domain
- 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.)
- Expired - Fee Related
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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention provides a kind of parameterized reconstruction algorithm of four side finite element mesh models, based on the network of quadrilaterals lattice model that finite element analysis software obtains, using the element number of network of quadrilaterals lattice model and node coordinate as the topological relation between input information reconstruction quadrilateral mesh unit, and the merging of unit grid specification is carried out according to obtained topology information, the union operations such as block domain self-adopt combination, subdomain after reducing block domain quantity and being merged, by retaining the characteristic point in the subzone boundaries for merging and obtaining and carrying out line to characteristic point, obtain the final subdomain of a certain number of quadrangles, the building of body parameter model is carried out according to the boundary of obtained final subdomain, adjustment and optimization, finally obtain body parameterized model.The algorithm can carry out parameterized reconstruction to the isostructural network of quadrilaterals lattice model of porous plate, and obtained body parameterized model is suitable for equal geometrical analysis, and have and be concise in expression, the advantages such as universe fairing.
Description
Technical field
The present invention relates to a kind of finite element model parameterized reconstruction algorithms, and in particular to a kind of four side finite element mesh models
Parameterized reconstruction algorithm.
Background technique
In current product design process, CAD is led over due to the development of CAE (computer-aided engineering) technology and (is calculated
Machine Computer Aided Design), therefore the Seamless integration- of the two is a bottleneck problem, about 70% runing time of CAE software is used for data
Exchange, and the generation phase workload of discrete grid block accounts for 80% or so of entire CAD/CAE process.
For the Seamless integration- for solving the problems, such as CAE and CAD, there are many analytic approach for application, medium geometry analysis method
Have many advantages, such as that strong flexibility, accuracy are high, have convergence, boundary condition setting succinct, and based on body parameterized model
Etc. geometric analysis methods have and solve the problems, such as the potentiality of CAE Yu CAD Seamless integration-.Body parameterized model is applied to wait geometry point
There is presently no too many researchs for analysis, but have been subjected to scholars and have more and more paid close attention to.
Construct parameter model is the first step of the geometric analysis methods such as above-mentioned.For zero deficiency entity, bivariate is created
(corresponding to areal model) or ternary (corresponding to threedimensional model) tensor body, makes itself and plane square or cube establish one
One mapping relations, the parameterized model obtained under both of these case are referred to as body parameterized model；For plurigenus entity, then will
It is decomposed into multiple zero deficiency subdomains and is handled, same available body parameterized model.
If using NUBRS batten during the creation of above-mentioned body parameterized model as mapping function, create
Model is batten body parameterized model.Using NUBRS batten as mapping function, the body parameterized model that may make has table
Up to succinct, universe fairing, it is convenient for and the advantages such as CAD is integrated, very big effect can be played in equal geometrical analysis and figure modelling field.
Currently, limited element analysis technique is the main stream approach that the field CAE carries out simulation analysis to CAD product, this method is utilized
Obtained grid model can be used for the building of above-mentioned body parameterized model.In the creation for carrying out grid model using this method
When, it can encounter and need the case where modeling to perforated plate construction, since a large amount of dense holes are deposited in these perforated plate constructions
The grid for dividing limited element analysis technique often exceeds the calculating capacity of computer, to be unfavorable for subsequent calculating
It carries out.
Summary of the invention
To solve the above problems, present invention employs following technical solutions：
The parameterized model algorithm for reconstructing of a kind of four side finite element mesh model, for being fitted to network of quadrilaterals lattice model
Body parameterized model for equal geometrical analysis is rebuild, the model that wherein network of quadrilaterals lattice model will have been created by finite element software
It carries out quadrangular mesh partition and obtains, by the number and node coordinate of quadrilateral mesh units all in network of quadrilaterals lattice model
The progress of parameterized model algorithm for reconstructing is used for as input information, which to be characterised by comprising following steps：
Step S1 rebuilds the topological relation between quadrilateral mesh unit according to input information, including cell abutment is closed
System, unit Rankine-Hugoniot relations；
Step S2 carries out the merging of grid cell specification according to topological relation obtained in step S1 and block domain is adaptively closed
And obtain final subdomain and final subzone boundaries；
Step S3 is inserted by carrying out control point in boundary Control point reverse and domain to final subdomain obtained in step S2
Value, obtains two-dimensional parameter domain model, the control point in two-dimensional parameter domain model is adjusted and is optimized, and is suitable for etc. several
The body parameterized model of what analysis；
Wherein, step S2 specifically comprises the following steps：
Step S2.1, grid cell specification merges, in conjunction with UV growth algorithm and relative edge's growth algorithm, to quadrilateral mesh
Unit carries out specification merging, obtains a certain number of basic blocks domains；
Step S2.2 carries out block domain self-adopt combination, including completely shared side block to the basic blocks domain that step S2.1 is obtained
Domain merges, most of shared side block domain merges and isolated island block domain merges, and obtains a certain number of subdomains；
Step S2.3, point characterized by borderline all angle points of the subdomain obtained in the step S2.2, keeping characteristics point,
It is remained unchanged according to outer boundary characteristic point, the rule that inner boundary characteristic point is adjusted to boundary Limit position adjusts characteristic point
It is whole；Characteristic point adjusted is connected with straight line, generates new quadrangle domain, obtains final subdomain and corresponding final subdomain
Boundary.
Further, in the parameterized model algorithm for reconstructing of four side provided by the invention finite element mesh model, step S1
Include the following steps：
Step S1.1 extracts essential information by analyzing the element number and node coordinate of quadrilateral mesh unit；
Step S1.2 traverses all quadrilateral mesh units, analyzes the shared letter of node and side of quadrilateral mesh unit
Breath, obtains all public side informations and identifies the syntople between quadrilateral mesh unit；
Step S1.3 selectes a quadrilateral mesh unit as starting, respectively from root in two mutually perpendicular directions
Search for obtain the Rankine-Hugoniot relations of quadrilateral mesh unit according to adjacency information, both direction be respectively defined as U to V to；
Step S1.4 classifies all quadrilateral mesh units, will be located at the quadrilateral mesh unit on boundary line
It is set as boundary quadrangle, remaining is internal quadrangle；Boundary quadrangle containing boundary angle point is set to boundary angle point four
Side shape, remaining is the general quadrangle in boundary；
Quadrangle and boundary by internal quadrangle midpoint adjacent quadrangular number not equal to side adjacent quadrangular number
Quadrangle definition of the adjacent quadrangular number not less than side adjacent quadrangular number is characterized quadrangle in quadrangle, by surrounding
The node definition that quadrangle is characterized quadrangle is characterized node, and marks feature quadrangle and characteristic node；Wherein, phase is put
Adjacent definition is two quadrangle one and only one common nodes, the adjacent definition in side be two quadrangles have and only there are two
Common node, putting two adjacent or adjacent side quadrangles is adjacent quadrangular.
Further, the step S1 of the parameterized model algorithm for reconstructing of four side provided by the invention finite element mesh model
In, UV growth algorithm includes the following steps：
According to unit Rankine-Hugoniot relations in step S1, starting quadrangle is selected, respectively along the U of starting quadrangle to, V to progress
Merge, until both direction all encounters feature quadrangle, merging terminates.
Further, the step S1 of the parameterized model algorithm for reconstructing of four side provided by the invention finite element mesh model
In, relative edge's growth algorithm includes the following steps：
During unit grid specification merges, when the U of quadrangle is to, V to when can not judge, using the quadrangle as
Quadrangle is originated, according to cell abutment relationship in step S1, two adjacent edges of selected starting quadrangle have been used as initial line, increase
The long side adjacent quadrangular opposite with initial line is played, until two sides belong to feature quadrangle on growing direction, merging terminates.
Further, the step S2.3 of the parameterized model algorithm for reconstructing of four side provided by the invention finite element mesh model
In, the Connection Step of characteristic point is：
Step S2.3.1 traverses all inner boundary characteristic point and outer boundary characteristic point, will be in recently with straight line
Edge feature point and outer boundary characteristic point connect；
Step S2.3.2 traverses all inner boundary characteristic points, will be apart from nearest inner boundary characteristic point connection with straight line
Get up, obtain a certain number of sides N shape, the numerical value of N is 4 to 20；
Step S2.3.3 generates straight line inside the shape of the side N according to the method that polygon is divided into quadrangle, generates new son
Domain, making subdomain is quadrangle domain.
Further, the parameterized model algorithm for reconstructing of four side provided by the invention finite element mesh model also has as follows
Technical characteristic：
In step s3, the method for boundary Control point reverse is：To the side in final subzone boundaries obtained in step S2
Boundary's point carries out B-spline curves fitting, obtains initial control point, then carries out node insertion to curve, rises rank, splicing, makes itself and phase
The knot vector of adjacent subzone boundaries curve unifies and same order, to obtain ultimate bound control point；
The method of control point interpolation is in domain：It is theoretical according to Coons curved surface, it sets four and surrounds into closing quadrangle domain
Spline curve P (u, 0), P (u, 1), P (0, v), P (1, v), u, v ∈ [0,1], corresponding control vertex areFour angle points are P (0,0), P (0,1), P (1,0), P (1,1), following institute
Show.Obtaining control point in four side curved surface domains is：
Invention action and effect
The parameterized model algorithm for reconstructing of four side finite element mesh models of offer according to the present invention, with finite element analysis
It is adaptive to grid model progress unit grid specification merging, block domain based on the network of quadrilaterals lattice model that software obtains
The union operations such as merging reduce block domain quantity, obtain a certain number of subdomains；By retaining the characteristic point and right in subzone boundaries
Characteristic point carries out line, obtains a certain number of final subdomains, carries out body parameter model according to the boundary of obtained final subdomain
Building, adjustment and optimization, finally obtain the body parameterized model for being suitable for equal geometrical analysis, the model have be concise in expression,
Universe fairing is convenient for and the advantages such as CAD is integrated, can play very big effect in equal geometrical analysis and figure modelling field.
Detailed description of the invention
Fig. 1 is the general flow chart of the parameterized model algorithm for reconstructing of four side provided by the invention finite element mesh model；
Fig. 2 is the common edge structural schematic diagram of grid cell；
Fig. 3 be UV to divide schematic diagram；
Fig. 4 is the schematic diagram that unit grid specification merges；
Fig. 5 is the preferential model of growth schematic diagram of quadrangle of UV growth algorithm；
Fig. 6 is the flow chart of UV growth algorithm；
Fig. 7 is that wrong schematic diagram is exchanged in the direction UV；
Fig. 8 is that completely shared side block domain merges schematic diagram；
Fig. 9 is the flow chart that completely shared side block domain merges；
Figure 10 is that most of shared side block domain merges schematic diagram；
Figure 11 (A), Figure 11 (B) are that isolated island block domain merges schematic diagram；
Figure 12 is characterized an adjustment schematic diagram；
Figure 13 is U and V to ruled surface schematic diagram；
Figure 14 is quadrangle domain control point schematic diagram；
Figure 15 is that body parameterizes grid model quality optimization flow chart；
Figure 16 is the nine orifice plate grid model instance graphs divided in ANSYS；
Figure 17 is the unit information table ELIST of nine orifice plate grid models derived from ANSYS；
Figure 18 is the informational table of nodes NLIST of nine orifice plate grid models derived from ANSYS；
Figure 19 is result schematic diagram of the nine orifice plate grid models after the merging of grid cell specification；
Figure 20 is the merging final result schematic diagram of nine orifice plate grid models；
Figure 21 is characteristic point schematic diagram after the adjustment of nine orifice plate grid models；
Figure 22 is the final subzone boundaries schematic diagram of nine orifice plate grid models；
Figure 23 is the running example figure of L shape Slab with holes；
Figure 24 is the running example figure of foot type Slab with holes；
Figure 25 is the running example figure of random nine hole Slab.
Specific embodiment
Illustrate a specific embodiment of the invention with reference to the accompanying drawings and embodiments.
【Embodiment one】
Figure 16 is the nine orifice plate grid model instance graphs divided in ANSYS.The model selected in embodiment one is rule point
Nine hole Slab of cloth, the model are one piece of square plate, and rule distributed nine through holes thereon, and as shown in figure 16, utilization is limited
Meta analysis software ANSYS analyzes the model, and obtained model is grid model, and grid cell therein is four sides
Shape.
Figure 17 is that nine orifice plate grid model unit information table ELIST, Figure 18 derived from ANSYS are nine holes derived from ANSYS
Expanded metals lattice model informational table of nodes NLIST.After obtaining network of quadrilaterals lattice model using finite element software, by finite element analysis software
The element number and node coordinate of obtained grid model are as input information, and as shown in figure 17, in NLIST, first is classified as
Node serial number, the x for being followed successively by node below, y, z coordinate；As shown in figure 18, in ELSIT, first is classified as element number, most
Four it is classified as node serial number afterwards.
Fig. 1 is the general flow chart of the parameterized model algorithm for reconstructing of four side provided by the invention finite element mesh model.
As shown in Figure 1, following the steps below after obtaining input information, body parameterized model is rebuild.
Step S1, grid model pretreatment.According to above-mentioned input information, the topological relation of reconstruction model is obtained relevant
Topology information, specific steps are as described in step S1.1~1.4.
Step S1.1, essential information are extracted, and the element number ELSIT and node coordinate of quadrilateral mesh unit are analyzed
NLIST extracts essential information；
Step S1.2, syntople extracts, the syntople by the shared information of node and side, between recognition unit.
Two kinds of syntoples between definition unit here：
Side is adjacent, and two units have and only there are two common nodes；
Point is adjacent, two unit one and only one common nodes.
Fig. 2 is the common edge structural schematic diagram of grid cell.As shown in Fig. 2, the adjacent unit in any two side must be deposited
In a common edge, unit 1 and unit 2 are respectively indicated with firstquad and secondquad, Node [0] and Node [1] are two
The common node of a unit.All grid cells are traversed, finally obtain all public side informations, and store spare.
Step S1.3, Rankine-Hugoniot relations extract, and select some quadrilateral mesh unit as start unit, adjacent using unit
Relationship searches for obtain the Rankine-Hugoniot relations of unit respectively according to syntople from both direction.
Fig. 3 be UV to divide schematic diagram.As shown in figure 3, introducing the direction UV storage rule in Rankine-Hugoniot relations extraction, appoint
Quadrangles at most adjacent there are four side around meaning quadrangle, defining 0 is the direction u+, and 2 be the direction u-, and 1 is the direction v+, and 3 be v-
Direction can be by adjacent cells around center cell rule compositor using the relationship of adjacent cells mass center and center cell mass center, will
Relevant information storage is spare.
Step S1.4, characteristic information mark.
Classify first to all quadrangles.In all units, become boundary quadrangle on boundary line, remaining
For internal quadrangle；In all boundary quadrangles, the angle point on four vertex boundary Zhong You is boundary angle point quadrangle, remaining is
The general quadrangle in boundary.
Then, special element present in network of quadrilaterals lattice model and node definition are characterized unit and characteristic node,
And it is labeled.Feature unit and characteristic node are defined as follows：
It, should if point adjacent quadrangular number is not equal to the quadrangle of side adjacent quadrangular number for internal quadrangle
Quadrangle is characterized quadrangle；It, should if adjacent quadrangular number is not less than side adjacent quadrangular number for boundary quadrangle
Boundary quadrangle is also characterized quadrangle.It, should if the quadrangle around a certain node is characterized quadrangle in all nodes
Node is characterized node.
Step S2, body parameterize domain boundary and generate, and according to topology information obtained in step S1, carry out a series of merging behaviour
Make, specific steps are as described in step S2.1~S 2.3.
Step S2.1, grid cell specification merge.
Fig. 4 is the schematic diagram that grid cell specification merges.As shown in figure 4, a selected boundary angle point quadrangle is to merge
Quadrangle is originated, is merged using UV growth algorithm and relative edge's growth algorithm.
Fig. 5 is the preferential model of growth schematic diagram of quadrangle of UV growth algorithm.UV growth algorithm is excellent with diagonal line model
First model of growth makes to originate quadrangle respectively along the merging of UV both direction, and diagonal line model of growth is required to merge each time and be obtained
Result be still quadrangle domain, as shown in figure 5, U is to a step is increased, then V to each unit require in V to increasing a step；
Similarly, V is to one step of every growth, U to all units be also required to increase a step upwards in U, to guarantee the figure increased
It is still quadrangle.
Fig. 6 is the flow chart of UV growth algorithm.As shown in fig. 6, defining index [0] [0] is to increase starting quadrangle,
UAddFlag is that U marks for V to growth to label, vAddFlag is increased, and initial value is true, is stopped if being set to false
Increase.U first judges whether current quadrangle is characterized quadrangle to increasing a step；If then uAddFlag is set to
Otherwise false and U continue to increase to stopping increasing；Then V judges whether current quadrangle is characterized four sides to a step is increased
Shape；If then vAddFlag is set to false and V is to stopping increasing, otherwise continue to increase.The growth label of UV both direction is set
For false, then merge termination, amalgamation result is stored spare.
Fig. 7 is that wrong schematic diagram is exchanged in the direction UV.Under some cases, circular arc boundary as shown in Figure 7, the UV of quadrangle
To can occur angulation change with circular arc, to the last direction complete exchange, UV has been to being unable to judge accurately at this time, if after
It is continuous that algorithm is merged using UV, it will lead to merging mistake, therefore selection relative edge merging algorithm continues to merge at this time.According to rise
The unit number of beginning quadrangle has been initial line by any two adjacent edges of quadrangle are originated according to relative edge's storage rule, is increased
Side adjacent quadrangular on this long two adjacent edges, judges whether the quadrangle on growing direction is characterized quadrangle, if then
Stop increasing, otherwise continues to increase, algorithm flow is with UV growth algorithm, and details are not described herein.
Figure 19 is result schematic diagram of the nine orifice plate grid models after the merging of grid cell specification.As shown in figure 19, pass through
Grid cell specification merges, and grid model originally merges into several basic blocks domains.
Step S2.2, block domain self-adopt combination.Merge by unit grid specification, element number is reduced to original 1/10
~1/30, a certain number of basic blocks domains are obtained, but the quantity in these basic blocks domains is still larger, has needed to continue to merge.By
In obtained basic blocks domain be quadrangle domain, therefore the four edges circle information in available each basic blocks domain, according to these
Boundary information carries out being faster than self-adopt combination, and it is described that steps are as follows.
Shared the merging of side block domain completely first.Merging rule is, if any two basic blocks Yu You a line circle is complete
Full weight is closed, then is merged the two basic blocks domains.
Fig. 8 is that completely shared side block domain merges schematic diagram.As shown in figure 8, the unit grid specification by step S2.1 is closed
And basic blocks domain 1~4 is obtained, and according to the rule that completely shared side block domain merges, domain 1 and domain 3 merge, and domain 2 and domain 4 merge, and two
The secondary obtained block domain that merges merges again, this final four block domains merge into one piece.
Fig. 9 is the flow chart that completely shared side block domain merges.As shown in figure 9, the process that completely shared side block domain merges is,
All basic blocks domains are traversed, searches with the presence or absence of the boundary being completely coincident between two block domains, if it exists then merges two pieces of domains,
Next piece of domain is continued to search, until the boundary being completely coincident all is not present in all pieces of domains.
Then most of shared side block domain is carried out to merge.Merging rule is all pieces of domains to be traversed, if the boundary in two pieces of domains has
Shared angle point, and boundary degree of share is more than a certain threshold value, then new block domain is merged into this two pieces of domains.Wherein, threshold value can root
It is set according to actual conditions.
Figure 10 is that most of shared side block domain merges schematic diagram.As shown in Figure 10, angle point is shared on the side in domain 1 and domain 2, and altogether
Enjoying has been more than threshold value, therefore domain 1 and domain 2 are merged；Similarly, domain 3 and domain 4 also meet the condition of merging, by domain 3 and domain 4
Merge.
Isolated island block domain is finally carried out to merge.By above-mentioned merging, block domain quantity is remarkably decreased, but still has some orphans in model
Vertical grid cell or block domain there are problems that not being merged or be completely included by other block domains, these grid cells and block domain
Referred to as isolated island block domain.The merging algorithmic rule in isolated island block domain is to identify isolated unit or completely by comprising block domain, be incorporated into
In adjacent maximum subdomain.
Figure 11 (A), Figure 11 (B) are that isolated island block domain merges schematic diagram.As shown in Figure 11 (A), isolated unit 1 and isolated unit
2 include by bulk domain 3 and bulk domain 4 respectively, and domain 5 includes by bulk domain 6, therefore is respectively merged into isolated unit 1 greatly
In block domain 3, isolated unit 2 is merged into bulk domain 4, and domain 5 is merged into bulk domain 6, amalgamation result such as Figure 11 (B).
Figure 20 is the merging final result schematic diagram of nine orifice plate grid models.As shown in figure 20, model is adaptive by block domain
It should merge, obtain a certain number of subdomains.
Step S2.3 generates final subzone boundaries.By above-mentioned merging, model partition is several quadrangle subdomains, but side
Boundary's shape is complex, is unfavorable for expressing.The point characterized by borderline all angle points of these subdomains, retains these characteristic points
And it is adjusted.
Figure 12 is characterized an adjustment schematic diagram.The rule of characteristic point adjustment is that outer boundary characteristic point remains unchanged, inner boundary
Characteristic point is adjusted to boundary Limit position.As shown in figure 12, side's point is characterized a little, needs to be adjusted, and dot is general side
Boundary's point, without adjustment.
Figure 21 is characteristic point schematic diagram after the adjustment of nine orifice plate grid models, and as shown in figure 21, keeping characteristics point simultaneously carries out
Adjustment, is adjusted rear characteristic point.
Characteristic point adjustment after the completion of, then by characteristic point connect quadrangularly domain.Connection Step such as step S2.3.1~
Described in S2.3.3.
Step S2.3.1 traverses all inner boundary characteristic point and outer boundary characteristic point, will be in recently with straight line
Edge feature point and outer boundary characteristic point connect；
Step S2.3.2 traverses all inner boundary characteristic points, will be apart from nearest inner boundary characteristic point connection with straight line
Get up, obtain a certain number of sides N shape, the numerical value of N is 4 to 20；
Step S2.3.3, according to polygon be divided into the method for quadrangle in above-mentioned steps obtained in inside the shape of the side N it is raw
It is in line, generates new subdomain, these subdomains are final subdomain, and are quadrangle domain.
Figure 22 is the final subzone boundaries schematic diagram of nine orifice plate grid models.As shown in figure 22, to characteristic point adjusted
Line operation is carried out, the final subdomain of several quadrangles is obtained.
Step S3, two-dimensional parameter domain model generates and optimization, and process is as described in step S3.1~3.2.
Step S3.1 controls point interpolation in the generation of two-dimensional parameter model, including boundary curve control point reverse and domain.It is right
Final subzone boundaries point obtained in step S2 carries out B-spline curves fitting, obtains initial control point, then saves to curve
Point insertion rises rank, splicing, keeps the knot vector of itself and adjacent subzone boundaries curve unified and same order, thus obtains ultimate bound
Control point.By carrying out controlling point interpolation in domain to boundary control point, initial internal control point is obtained.
Figure 13 is U and V to ruled surface schematic diagram, and Figure 14 is quadrangle domain control point schematic diagram.As shown in figure 13, in initial
Portion control point is generated using Coons surface interpolation theory, sets four spline curve P (u, 0) surrounded into closing quadrangle domain,
P (u, 1), P (0, v), P (1, v), u, v ∈ [0,1], corresponding control vertex areFour angle points are P (0,0), P (0,1), P (1,0), P (1,1), following institute
Show.Obtaining control point in four side curved surface domains is：
Step S3.2, model adjustment and optimization.By the process that above-mentioned model generates, obtained output result is body parameter
Change model internal control point, but these control points have the defects of distortion and selfing, influence subsequent analysis, it is therefore desirable to control
Point carries out reasonably optimizing.
Figure 15 is that body parameterizes mesh quality optimized flow chart.Initialization according to process as shown in figure 15, to obtaining
Each point control of model is denoted as T using the method divided and rule by the control point in grid modelH, time numerical value α=H is set,
The integration constraint equation at each control point is optimized respectively, every numerical value for solving primary then α subtracts 1 and continues to solve,
Until α=1, terminates after carrying out last time optimization, finally obtain and export body parameterized model.
【Embodiment two】
In the present embodiment two, the explanation with identical algorithms and step in embodiment one is omitted.
Figure 23 is the running example figure of L shape Slab with holes, which is one piece of L shaped plate, and rule distributed four thereon
Through hole, as shown in figure 23, the grid model obtained from finite element analysis software is by merging, characteristic point retains and adjusts, is special
A sign point line, obtains the final subdomain of several quadrangles, carries out the generations of two-dimensional parameter model, adjustment and excellent according to these final subdomains
Change, can be obtained and export the body parameterized model of geometrical analysis such as being suitable for.
【Embodiment three】
Figure 24 is the running example figure of foot type Slab with holes, which is one block of foot type plate, is distributed with thereon one big
Through hole and three small through holes, as shown in figure 24, the grid model obtained from finite element analysis software is by merging, characteristic point
Retain and adjust, characteristic point line, obtain the final subdomain of several quadrangles, carries out two-dimensional parameter model according to these final subdomains
It generates, adjust and optimizes, can be obtained and export the body parameterized model for being suitable for equal geometrical analysis.
【Example IV】
Figure 25 is the running example figure of random nine hole Slab, which is one piece of square plate, is divided at random thereon
Be furnished with nine through holes, as shown in figure 25, the grid model that obtains from finite element analysis software is by merging, characteristic point retain with
Adjustment, characteristic point line, obtain the final subdomain of several quadrangles, according to these final subdomains carry out the generations of two-dimensional parameter model,
Adjustment and optimization, can be obtained and export the body parameterized model of geometrical analysis such as being suitable for.
Embodiment action and effect
The parameterized model algorithm for reconstructing of the four side finite element mesh models provided according to the present invention, in example 1,
The body parameterized model based on four side finite element mesh models is carried out to nine hole Slabs to rebuild, as a result as shown in figure 22, warp
To cross merging step can effectively reduce subdomain quantity, retain using characteristic point and adjust, characteristic point line, can obtain
To the final subdomain of a certain number of quadrangles, two-dimensional parameter model is carried out to the boundary of these subdomains and generates, optimize and adjusts,
It can be obtained and export body parameterized model, can be applied to subsequent Deng geometrical analysis.Similarly, in embodiment two, implementation
In example three and example IV, algorithm provided by the invention can trouble-free operation, reduce the quantity of unit grid and obtain final
Subdomain is finally generated by two-dimensional parameter model, optimization and the operation of adjustment aspect can be obtained and exports body parametrization mould
Type, these body parameterized models, which all have, to be concise in expression, universe fairing, is convenient for and advantage that CAD is integrated, can be used in subsequent
Equal geometrical analysis.
The parameterized model that above embodiments are merely to illustrate four side provided by the invention finite element mesh model, which is rebuild, to be calculated
The embodiment and operation result of method, but the present invention is not merely defined in and carries out in the above-described embodiments, for other types
Perforated plate construction and grid model in the less general structure of grid cell quantity, algorithm provided by the invention is also the same
Effectively.
Claims (6)
1. a kind of parameterized model algorithm for reconstructing of four side finite element mesh models, applicable for carrying out to network of quadrilaterals lattice model
It is rebuild in the body parameterized model of equal geometrical analysis, wherein the mould that the quadrilateral mesh model will have been created by finite element software
Type carries out quadrangular mesh partition and obtains, by the number and section of quadrilateral mesh units all in the network of quadrilaterals lattice model
Point coordinate is used for the progress of the parameterized model algorithm for reconstructing as input information, which is characterised by comprising as follows
Step：
Step S1 rebuilds the topological relation between quadrilateral mesh unit according to the input information, including cell abutment is closed
System, unit Rankine-Hugoniot relations；
Step S2, the topological relation carries out the merging of grid cell specification according to obtained in step S1 and block domain is adaptively closed
And obtain final subdomain and final subzone boundaries；
Step S3 is inserted by carrying out control point in boundary Control point reverse and domain to the final subdomain obtained in step S2
Value, obtains two-dimensional parameter domain model, the control point in the two-dimensional parameter domain model is adjusted and is optimized, is suitable for
The body parameterized model of equal geometrical analysis；
Wherein, step S2 specifically comprises the following steps：
Step S2.1, grid cell specification merges, in conjunction with UV growth algorithm and relative edge's growth algorithm, to the quadrilateral mesh
Unit carries out specification merging, obtains a certain number of basic blocks domains；
Step S2.2 carries out block domain self-adopt combination, including completely shared side block to the basic blocks domain that step S2.1 is obtained
Domain merges, most of shared side block domain merges and isolated island block domain merges, and obtains a certain number of subdomains；
Step S2.3, point characterized by borderline all angle points of the subdomain obtained in the step S2.2, retains the spy
Point is levied, is remained unchanged according to outer boundary characteristic point, inner boundary characteristic point is adjusted to the rule of boundary Limit position to the feature
Point is adjusted；The characteristic point adjusted is connected with straight line, generates four new side domains, obtains final subdomain and corresponding
Final subzone boundaries.
2. the parameterized model algorithm for reconstructing of four side according to claim 1 finite element mesh model, which is characterized in that institute
Step S1 is stated to include the following steps：
Step S1.1, by analyzing the information table of the element number of quadrilateral mesh unit and the information of the node coordinate
Table extracts the essential information including the element number and the node coordinate；
Step S1.2 traverses all quadrilateral mesh units, analyzes the node and side shared information of quadrilateral mesh unit, obtains
To all public side informations and identify the syntople between quadrilateral mesh unit；
Step S1.3 selectes a quadrilateral mesh unit as starting, respectively according to institute from two mutually perpendicular directions
State adjacency information and search for obtain the Rankine-Hugoniot relations of quadrilateral mesh unit, described two directions be respectively defined as U to V to；
Step S1.4 classifies all quadrilateral mesh units, will be located at the quadrilateral mesh unit on boundary line and sets
For boundary quadrangle, remaining is internal quadrangle；The boundary quadrangle containing boundary angle point is set to boundary angle point four
Side shape, remaining is the general quadrangle in boundary；
Will the internal quadrangle midpoint adjacent quadrangular number not equal to the quadrangle of side adjacent quadrangular number and described
Quadrangle definition of the adjacent quadrangular number not less than side adjacent quadrangular number is characterized quadrangle in the quadrangle of boundary, will be all
The node definition that the quadrangle enclosed is characterized quadrangle is characterized node, and marks the feature quadrangle and the feature section
Point；Wherein, putting adjacent definition is two quadrangle one and only one common nodes, and the adjacent definition in side is two quadrangles
Have and only there are two common nodes, putting two adjacent or adjacent side quadrangles is adjacent quadrangular.
3. the parameterized model algorithm for reconstructing of four side according to claim 2 finite element mesh model, which is characterized in that
In step S2.1, the UV growth algorithm includes the following steps：
According to unit Rankine-Hugoniot relations described in step S2.1, starting quadrangle is selected, respectively along the U of the starting quadrangle
To, the V to merging, until both direction all encounters feature quadrangle, merging terminates.
4. the parameterized model algorithm for reconstructing of four side according to claim 2 finite element mesh model, which is characterized in that
In step S2.1, relative edge's growth algorithm includes the following steps：
During the unit grid specification merges, when the U of quadrangle is to, V to when can not judge, using the quadrangle as
Quadrangle is originated, according to cell abutment relationship described in step S1, two adjacent edges for selecting the starting quadrangle are used as
Initial line increases the side adjacent quadrangular opposite with described initial line, until two sides belong to feature quadrangle on growing direction,
Merging terminates.
5. the parameterized model algorithm for reconstructing of four side according to claim 1 finite element mesh model, which is characterized in that
In step S2.3, the Connection Step of the characteristic point is：
Step S2.3.1 traverses all inner boundary characteristic point and outer boundary characteristic point, will be apart from nearest inner boundary with straight line
Characteristic point and outer boundary characteristic point connect；
Step S2.3.2 traverses all inner boundary characteristic points, will be connected apart from nearest inner boundary characteristic point with straight line,
A certain number of sides N shape is obtained, the numerical value of N is 4 to 20；
Step S2.3.3 generates straight line inside the shape of the side N according to the method that polygon is divided into quadrangle, generates new son
Domain, making the subdomain is four side domains.
6. the parameterized model algorithm for reconstructing of four side according to claim 1 finite element mesh model, it is characterised in that：
In step s3, the method for the boundary Control point reverse is：To in final subzone boundaries obtained in the step S2
Boundary point carry out B-spline curves fitting, obtain initial control point, then node insertion is carried out to curve, rises rank, splicing, make it
The unified and same order with the knot vector of adjacent subzone boundaries curve, to obtain ultimate bound control point；
The method of control point interpolation is in the domain：It is theoretical according to Coons curved surface, set four samples surrounded into four side domains of closing
Curve P (u, 0), P (u, 1), P (0, v), P (1, v), u, v ∈ [0,1], P (1,1), corresponding control vertex areFour angle points are P (0,0), P (0,1), P (1,0), and P (1,1) is as follows
Shown in figure, obtaining control point in four side curved surface domains is：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610307722.0A CN105787226B (en) | 2016-05-11 | 2016-05-11 | The parameterized model of four side finite element mesh models is rebuild |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610307722.0A CN105787226B (en) | 2016-05-11 | 2016-05-11 | The parameterized model of four side finite element mesh models is rebuild |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105787226A CN105787226A (en) | 2016-07-20 |
| CN105787226B true CN105787226B (en) | 2018-11-20 |
Family
ID=56401283
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610307722.0A Expired - Fee Related CN105787226B (en) | 2016-05-11 | 2016-05-11 | The parameterized model of four side finite element mesh models is rebuild |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105787226B (en) |
Families Citing this family (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106384384B (en) * | 2016-09-18 | 2020-05-05 | 上海理工大学 | Shape optimization method of three-dimensional product model |
| CN106910216B (en) * | 2017-03-03 | 2018-08-03 | 广东工业大学 | A kind of grid surface based on ship hull plate point cloud data determines method and device |
| CN107085865B (en) * | 2017-05-12 | 2020-10-16 | 杭州电子科技大学 | Quadrilateral segmentation method applied to finite element analysis |
| CN107451371B (en) * | 2017-08-11 | 2021-01-26 | 四川大学 | Method for improving calculation precision of triangular unit |
| CN107767452A (en) * | 2017-10-10 | 2018-03-06 | 上海理工大学 | The AMF general file generation methods of heterogeneous solid parameterized model |
| CN109033601B (en) * | 2018-07-18 | 2023-03-17 | 雅马哈发动机（厦门）信息系统有限公司 | CAE design method for pre-analyzing boundary conditions |
| CN110188423B (en) * | 2019-05-16 | 2022-08-09 | 广西交通设计集团有限公司 | Linear engineering structure rapid BIM modeling method based on finite element meshing |
| CN111984746B (en) * | 2019-05-21 | 2024-04-05 | 北京京东乾石科技有限公司 | Method and device for generating simulation map based on grid map |
| CN112329282B (en) * | 2019-07-31 | 2023-12-12 | 比亚迪股份有限公司 | Method, device and storage medium for classifying hole perimeter grids |
| CN111161129B (en) * | 2019-11-25 | 2021-05-25 | 佛山欧神诺云商科技有限公司 | Three-dimensional interaction design method and system for two-dimensional image |
| CN111177837B (en) * | 2019-12-30 | 2023-08-22 | 北京禹数技术有限公司 | Data processing method and device based on three-dimensional building model and electronic equipment |
| CN111223185B (en) * | 2019-12-30 | 2023-10-27 | 苏州数设科技有限公司 | Grid cell information display method of three-dimensional model, device and electronic equipment thereof |
| CN111177968B (en) * | 2020-01-02 | 2022-01-18 | 大连理工大学 | Method for quickly inserting two-dimensional zero-thickness cohesive force unit into target area |
| CN111104764B (en) * | 2020-01-21 | 2020-08-25 | 湖南科技大学 | Structured grid optimization division method for alternating current motor rotor conducting bar thermal analysis model |
| CN111581776B (en) * | 2020-04-07 | 2023-09-26 | 华南理工大学 | Iso-geometric analysis method based on geometric reconstruction model |
| CN111797555B (en) * | 2020-07-01 | 2022-06-24 | 上海理工大学 | Geometric reconstruction method based on finite element model |
| CN112182489B (en) * | 2020-10-14 | 2023-11-24 | 杭州电子科技大学 | Two-dimensional high-order grid generation method based on partial differential equation solution |
| CN114627265A (en) * | 2021-08-09 | 2022-06-14 | 上海理工大学 | Volume parameterization modeling method based on curve quadrilateral subdivision |
| CN113792458B (en) * | 2021-09-09 | 2023-11-14 | 中国航天科工集团第二研究院 | Optimization method and device for finite element triangle mesh |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101655992A (en) * | 2009-09-15 | 2010-02-24 | 中国科学院软件研究所 | Reconstruction method of triangular mesh model |
| CN102034272A (en) * | 2010-09-29 | 2011-04-27 | 浙江大学 | Generating method of individualized maxillofacial soft tissue hexahedral mesh |
| CN102306396A (en) * | 2011-09-15 | 2012-01-04 | 山东大学 | Three-dimensional entity model surface finite element mesh automatic generation method |
| CN104157013A (en) * | 2014-08-06 | 2014-11-19 | 杭州新迪数字工程系统有限公司 | Reestablishing method for smooth curved surface on arbitrary quadrilateral grid |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110310101A1 (en) * | 2010-06-22 | 2011-12-22 | Schlumberger Technology Corporation | Pillar grid conversion |
-
2016
- 2016-05-11 CN CN201610307722.0A patent/CN105787226B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101655992A (en) * | 2009-09-15 | 2010-02-24 | 中国科学院软件研究所 | Reconstruction method of triangular mesh model |
| CN102034272A (en) * | 2010-09-29 | 2011-04-27 | 浙江大学 | Generating method of individualized maxillofacial soft tissue hexahedral mesh |
| CN102306396A (en) * | 2011-09-15 | 2012-01-04 | 山东大学 | Three-dimensional entity model surface finite element mesh automatic generation method |
| CN104157013A (en) * | 2014-08-06 | 2014-11-19 | 杭州新迪数字工程系统有限公司 | Reestablishing method for smooth curved surface on arbitrary quadrilateral grid |
Non-Patent Citations (2)
| Title |
|---|
| "Determination of nonlinear behavior of a ball joint system by model updating";Mohammad RezaDavoodi etc.;《Journal of Constructional Steel Research》;20120831;第71卷;Pages 52-62 * |
| "复杂机械零件的六面体有限元网格生成方法";曹华军等;《机械工程学报》;20140831;第50卷(第15期);第114-118页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105787226A (en) | 2016-07-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105787226B (en) | The parameterized model of four side finite element mesh models is rebuild | |
| CN109977507B (en) | Efficient generation method of porous scaffold with complex morphology | |
| CN109377561A (en) | A kind of digital-to-analogue surface mesh generation method based on conformal geometry | |
| CN109461209B (en) | Novel structural grid generation method | |
| CN106981097B (en) | A kind of T spline surface approximating method based on subregion Local Fairing weight factor | |
| CN103970949B (en) | A kind of rapid shaping intermediate cam patch model by side layered approach | |
| Mencl | A graph–based approach to surface reconstruction | |
| CN102436217B (en) | Method for reconstructing web processing drive geometry of slot characteristic of aircraft structure | |
| CN106023297A (en) | Texture dynamic organization method for fine three-dimensional model | |
| CN101853526A (en) | Density self-adapting non-uniform point cloud simplifying treatment method | |
| CN108389263A (en) | The IGES surface grids rapid generations calculated are solved towards Element BEM | |
| CN115481497A (en) | Volume parameterization modeling method based on feature framework | |
| CN114239356A (en) | Conformal lattice material design method based on finite element grid | |
| CN117708904A (en) | Mixed unit grid generation method oriented to semiconductor numerical simulation | |
| CN113792458B (en) | Optimization method and device for finite element triangle mesh | |
| CN116958385B (en) | Material texture dynamic updating method supporting mass monomer models, storage medium and equipment | |
| CN110349262A (en) | A kind of tunnel geological forecast three-dimensional modeling method and system based on BIM | |
| CN113378279A (en) | Building modeling method, building modeling system, electronic equipment and storage medium | |
| CN103617291B (en) | Equivalent characterization method for reservoir stratum cause unit interface | |
| CN103886635A (en) | Self-adaption LOD model establishing method based on face clustering | |
| CN106780747B (en) | A kind of method that Fast Segmentation CFD calculates grid | |
| CN117333630A (en) | Local implicit incremental updating and automatic reconstructing method for three-dimensional geological model | |
| CN106097446A (en) | A kind of three dimensional geological curved surface method for reconstructing based on extensive geological data | |
| CN114119628B (en) | Point cloud accurate segmentation method based on feature template | |
| CN108986212A (en) | Three-dimensional virtual terrain L OD model generation method based on crack elimination |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181120 Termination date: 20210511 |