CN105677965B - A kind of board-like satellite structure conode grid rapid generation - Google Patents
A kind of board-like satellite structure conode grid rapid generation Download PDFInfo
- Publication number
- CN105677965B CN105677965B CN201610005236.3A CN201610005236A CN105677965B CN 105677965 B CN105677965 B CN 105677965B CN 201610005236 A CN201610005236 A CN 201610005236A CN 105677965 B CN105677965 B CN 105677965B
- Authority
- CN
- China
- Prior art keywords
- face
- finite element
- max
- interface
- function
- 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Complex Calculations (AREA)
Abstract
A kind of board-like satellite structure conode grid rapid generation, it includes the following steps：Step 1：Face in satellite structure geometrical model extraction to importing finite element modeling software；Step 2：Obtain middle face set to be treated and its element number N；Step 3：To all related operation is carried out in face of (A, B)；Step 4：All faces are sutured；Step 5：Finite element grid is divided to the geometrical model of structure.The obtained geometry in face automatically processes in capable of being extracted to satellite structure three-dimensional entity model by this method, the gap between elimination face and face and deletes extra part, and the geometrical model ultimately generated can divide the conode finite element grid of high quality.A large amount of man-machine interactive operation is avoided, the time needed is few, versatile, applied widely, has very strong engineering practical value.
Description
Technical field：
The present invention relates to a kind of board-like satellite structure conode grid rapid generations, belong to structural analysis of spacecraft
Technical field.
Background technology：
Satellite structure mostly uses honeycomb sandwich panel and bears and transmit mechanical loading at present, and provides installation branch for instrument and equipment
Support.Whether met the requirements for the design of verifying satellites structure, mechanical property is generally carried out to it using finite element method in engineering
It can analysis.The geometrical model of satellite structure is generally modeled using 3D solid, but plate is generally used in finite element modeling
Shell structure unit carries out mechanical simulation.If dividing finite element grid in directly being extracted by 3D solid behind face, due to adjacent
Gap between face can lead to the discontinuous of finite element grid；If carrying out unit connection, a side using rigid coupling unit
Face requires the grid comparison match on adjacent surface in favor of creating rigid coupling unit, the establishment of another aspect rigid coupling unit
Additional rigidity can be introduced to structure；In addition to this, the establishment of rigid coupling unit can also increase the calculation amount of finite element analysis,
The dynamic analyses situation such as common model analysis and frequency response analysis in being analyzed especially for satellite structure.Therefore engineering
In wish to generate conode finite element grid.The Major Difficulties that satellite structure conode grid quickly generates are：Due to bee
The presence of nest interlayer plate thickness, from the geometrical model obtained behind face in three-dimensional entity model extraction, between existing between adjacent surface
Gap can not directly generate finite element grid.In order to realize the division of conode grid, opposite first is needed to be extended or cut out
It cuts, to eliminate gap or the extra part between face and face.Due to generally there is tens structure plates in satellite structure,
It relies only on the geometric manipulations function of business finite element software itself and realizes that entire satellite structure saves altogether by man-machine interaction mode
The generation of point finite element grid will need higher skill, a large amount of time and loaded down with trivial details work.
Invention content：
1, purpose：In order to realize efficiently producing for satellite structure finite element model, the object of the present invention is to provide a kind of plates
Formula satellite structure conode grid rapid generation.This method is defended first with the middle face abstraction function of finite element software
Then the middle face geometrical model of star structural slab is faced by what traversal acquisition was handled, secondly to every group in face of successively
It is automatically processed, to obtain the geometrical model of entire satellite structure, the geometry of finite element software is utilized to suture function again
Realize that by the geometrical model suture of satellite structure be an entirety.Mesh generation is carried out on this basis, obtains being based on plate shell list
The high quality conode finite element grid that member is simulated.Since this method realization automatically processes, a large amount of people is avoided
Machine interactive operation, the time needed is few, versatile, applied widely.
2, technical solution：The technical solution adopted by the present invention to solve the technical problems is：
A kind of board-like satellite structure conode grid rapid generation of the present invention, it is characterised in that it includes following step
Suddenly：
Step 1：Face in satellite structure geometrical model extraction to importing finite element modeling software；
Step 2：Obtain middle face set to be treated and its element number N：
Step 3：To all following operation is carried out in face of (A, B)：
Step 3.1：Calculate the angle theta between the normal vector of A and BAB；
Step 3.2：If 10 ° of ＜ θAB170 ° of ＜, thens follow the steps 3.3 to 3.10；
Step 3.3：Find out in A the side in need extended to B and its number NA；
Step 3.4：If NA＞ 0, thens follow the steps 3.5 to 3.6；It is no to then follow the steps 3.7 to 3.10；
Step 3.5：The side calculated in A needs the length δ extended to BA；
Step 3.6：By in A it is in need to B extend side extend δA；
Step 3.7：Find out in B the side in need extended to A and its number NB；
Step 3.8：If NB＞ 0, thens follow the steps 3.7；It is no to then follow the steps 3.9；
Step 3.9：The side calculated in B needs the length δ extended to AB；
Step 3.10：By in B it is in need to A extend side extend δB；
Step 3.11：Divide A using B；
Step 3.12：Delete the face that area is smaller in A；
Step 3.13：Divide B using A；
Step 3.14：Delete the face that area is smaller in B；
Step 4：All faces are sutured；
Step 5：Finite element grid is divided to the geometrical model of structure；
Wherein, " face in the satellite structure geometrical model extraction to importing finite element modeling software " in step 1
Method be：Existing finite element software generally all supports the middle face abstraction function of thin-walled parts, and satellite cellular sandwich plate
Structure is typically in homogeneous thickness, therefore directly selects honeycomb interlayer in dummy satellite by the user interface of finite element software
Face extraction can be obtained corresponding middle face geometrical model during the geometrical model of structural slab carries out；
Wherein, the method for described " the obtaining middle face set to be treated and its element number N " is in step 2：It is existing
Some finite element softwares generally all provide second development interface and user are allow to develop new function, therefore can utilize correlation
Function interface the set and its element number that all middle faces in model are formed are obtained by secondary development.Face collection in acquisition
Closing required second development interface function can be obtained by inquiring the dependent help document of finite element software；
Wherein, " angle theta between the normal vector of A and B is calculated described in step 3.1AB" method be：Assuming that A and B
Unit normal vector be respectively nA=(nAx,nAy,nAz) and nB=(nBx,nBy,nBz), then the cosine value of angle is cos between them
θAB=nAxnBx+nAynBy+nAznBz, to using the antitrigonometric function in elementary algebra can obtain two vectors between angle be
θAB=arccos (nAxnBx+nAynBy+nAznBz)；
Wherein, described in step 3.3 " find out in A the side in need extended to B and its number NA" method be：
Existing finite element software generally all provides second development interface and user is allow to develop new function, therefore can utilize phase
The function interface of pass is by secondary development and realizes in accordance with the following methods：1) all side E of A are obtainedA1,…,EAnWith one of B
Vertex VB；2) to each edge EAj(j=1 ..., n), it is assumed that its natural coordinates is expressed as Pj(s) (0≤s≤1), equably obtains Nj
+ 1 point Pjk=Pj(k/Nj) (k=0,1 ..., Nj), wherein Nj≥1000；3) d is takenmax=-999999, d=0, then according to first
Formula d in equal analytic geometryjk=nBx(xjk-xB)+nBy(yjk-yB)+nBz(zjk-zB) each point P is calculated successivelyjkTo having for B
To distance, if | djk| ＞ dmax, then dmax=| djk|, d=djk；Wherein (xjk,yjk,zjk) it is PjkSeat in rectangular coordinate system
Mark, (xB,yB,zB) it is VBCoordinate, (nBx,nBy,nBz) be B unit normal vector；4) abd=999999 is taken, to each point Pjk,
If ddjk＞ 0 and | djk| ＜ abd, then abd=| djk|；If 5) abd ＜ 0.04, δmax=abd+0.01；Otherwise δmax=
0.05；6)NA=0, to each edge EAj：If to k=0,1 ..., NjHave | djk|≤δmax, then EAjLabeled as 1, NA=NA+1；It is no
Then EAjLabeled as 0；Obtaining the vertex on all sides and B A in this step can be realized by the interface of software, and had
Side or curve can be generally indicated using natural coordinates in limit meta software, can be obtained by interface after giving its natural coordinates
Its coordinate representation under rectangular coordinate system.Required second development interface function can be by inquiring finite element in this step
The dependent help document of software obtains；
Wherein, described in step 3.5 " calculate A in side need the length δ extended to BA" method be：δA=
δmax/sinθAB, wherein δmaxValue obtained in step 3.3, θABValue obtained in step 3.1；
Wherein, described in step 3.6 " by A it is in need to B extend side extend δA" method be：It is existing
Finite element software generally all provides second development interface and user is allow to develop new function, therefore can utilize relevant letter
Number interfaces by secondary development by A institute it is in need to B extend side extension δA.Several sides in one face are extended
Required second development interface function can be obtained by inquiring the dependent help document of finite element software；
Wherein, described in step 3.7 " find out in B the side in need extended to A and its number NB" method be：
Existing finite element software generally all provides second development interface and user is allow to develop new function, therefore can utilize phase
The function interface of pass is by secondary development and realizes in accordance with the following methods：1) all side E of B are obtainedB1,…,EBnWith one of A
Vertex VA；2) to each edge EBj(j=1 ..., n) assume that its natural coordinates is expressed as Pj(s) (0≤s≤1), equably obtains Nj+
1 point Pjk=Pj(k/Nj) (k=0,1 ..., Nj) wherein Nj≥1000；3) d is takenmax=-999999, d=0, then according to first
Formula d in equal analytic geometryjk=nAx(xjk-xA)+nAy(yjk-yA)+nAz(zjk-zA) each point P is calculated successivelyjkTo having for A
To distance, if | djk| ＞ dmax, then dmax=| djk|, d=djk；Wherein (xjk,yjk,zjk) it is PjkSeat in rectangular coordinate system
Mark, (xA,yA,zA) it is VACoordinate, (nAx,nAy,nAz) be A unit normal vector；4) abd=999999 is taken, to each point Pjk,
If ddjk＞ 0 and | djk| ＜ abd, then abd=| djk|；If 5) abd ＜ 0.04, δmax=abd+0.01；Otherwise δmax=
0.05；6)NB=0, to each edge EBj：If to k=1 ..., N has | djk|≤δmax, then EBjLabeled as 1, NB=NB+1；Otherwise
EBjLabeled as 0；Obtaining the vertex on all sides and A B in this step can be realized by the interface of software, and limited
Side or curve can be generally indicated using natural coordinates in meta software, it can be obtained by interface after giving its natural coordinates
Coordinate representation under rectangular coordinate system.Required second development interface function can be soft by inquiring finite element in this step
The dependent help document of part obtains；
Wherein, described in step 3.9 " calculate B in side need the length δ extended to AB" method be：δB=
δmax/sinθAB, wherein δmaxValue obtained in step 3.7, θABValue obtained in step 3.1；It should be noted that
Although being all made of variable δ in this step and step 3.5max, but its value has changed in step 3.7.
Wherein, described in step 3.10 " by B it is in need to A extend side extend δB" method be：It is existing
Finite element software second development interface is generally all provided user is allow to develop new function, therefore can utilize relevant
Function interface by secondary development by B it is in need to A extend side extend δB.Several sides in one face are prolonged
Stretching required second development interface function can be obtained by inquiring the dependent help document of finite element software；
Wherein, the method for " dividing A using B " described in step 3.11 is：Existing finite element software generally all carries
Allow user to develop new function for second development interface, therefore secondary development can be passed through using relevant function interface
It realizes and divides A using B.It is B that divisional plane need to be only specified in second development interface function, and the face that is divided is A, required
Second development interface function can be obtained by inquiring the dependent help document of finite element software；
Wherein, the method for " deleting the face that area is smaller in A " described in step 3.12 is：Utilize finite element software
Second development interface can obtain the area in each face in A, can calculate the total of all faces in A by simple addition in this way
Area；If the ratio between the area in some face and the gross area are less than given threshold gamma, using second development interface function by its
It deletes；
Wherein, the method for " dividing B using A " described in step 3.13 is：Existing finite element software generally all carries
Allow user to develop new function for second development interface, therefore secondary development can be passed through using relevant function interface
It realizes and divides B using A；It is A that divisional plane need to be only specified in second development interface function, and the face that is divided is B, required
Second development interface function can be obtained by inquiring the dependent help document of finite element software；
Wherein, the method for " deleting the face that area is smaller in B " described in step 3.14 is：Utilize finite element software
Second development interface can obtain the area in each face in B, can calculate the total of all faces in B by simple addition in this way
Area；If the ratio between the area in some face and the gross area are less than given threshold gamma, using second development interface function by its
It deletes；
Wherein, the method for described " being sutured to all faces " is in step 4：Existing finite element software is general
All support that by several adjacent face sutures be an entirety, two such face it is adjacent when being changing to shared, grid is drawn
Timesharing will generate node, and the unit shared node between adjacent surface on shared side automatically.Therefore directly by limited
The curved surface of meta software sutures function, the face of selection institute grid division in need, then progress curved surface suture；
Wherein, the method for described " the dividing finite element grid to the geometrical model of structure " is in steps of 5：It is existing to have
Limit meta software generally all supports the grid division on curved surface, and user only needs designating unit type, size of mesh opening and selects to need
Want the curved surface of grid division that the automatic division of grid can be realized.Therefore directly by the mesh generation function of finite element software,
And designating unit type can carry out grid for curved surface all in triangle or quadrangle, given size of mesh opening and preference pattern
It divides.It is an entirety due to suturing all curved surfaces in step 4, the grid obtained in this step will be continuous
, the grid between adjacent surface is conode.
By above each step, realizes board-like satellite structure conode grid and quickly generate, and in the generation of grid
In the process based on the three-dimensional entity model of satellite structure without creating new model, the grid of generation specifically for mechanical analysis
To be continuous, without additionally introducing rigid coupling unit.In addition, due to being realized using batch processing, avoid higher
Skill, a large amount of time and loaded down with trivial details work.
3, advantage and effect：The beneficial effects of the invention are as follows：It is quick to provide a kind of board-like satellite structure conode grid
Generation method, the geometry that face obtains in being extracted to satellite structure three-dimensional entity model by this method automatically process,
Gap between elimination face and face and extra part is deleted, the geometrical model ultimately generated can divide the total section of high quality
Point finite element grid.Since this method realization automatically processes, a large amount of man-machine interactive operation is avoided, the time needed
It is few, it is versatile, it is applied widely.Due to the finite element grid conode that this method generates, need not create additional
Rigid coupling unit, avoid the introducing of additional stiffness and additional calculation amount.In addition, since the application of board-like satellite structure is wide
It is general, therefore the present invention has very strong engineering practical value in satellite structure finite element modeling.
Illustrate a kind of board-like satellite structure conode grid rapid generation
Fig. 1 the method for the invention flow charts.
Fig. 2 face A and the plane where the B of face are non-intersecting, and the minimum range of the point on the A of face and face B is more than δmax。
Fig. 3 face A and the plane where the B of face are non-intersecting, and the minimum range of the point on the A of face and face B is no more than δmax。
Fig. 4 face A intersect with the plane where the B of face.
Fig. 5 face B and the plane where the A of face are non-intersecting, and the minimum range of the point on the B of face and face A is more than δmax。
Fig. 6 face B and the plane where the A of face are non-intersecting, and the minimum range of the point on the B of face and face A is no more than δmax。
Fig. 7 face B intersect with the plane where the A of face.
The side extended is needed in Fig. 8 extended surfaces A.
The side extended is needed in Fig. 9 extended surfaces B.
Figure 10 utilize face B divisional planes A.
Figure 11 delete the face of area very little in A.
Figure 12 delete the face of area very little in B.
Symbol description is as follows in figure：
I, j are respectively two counters；
nA,nBThe respectively normal vector of face A and face B；
θABFor nA,nBBetween angle；
N is the number in middle face to be treated in model；
NA,NBThe number on the side extended is needed in respectively face A and face B；
δmaxFor the distance threshold for calculating or giving；
δAThe length for needing to extend to B for the side in A；
δBThe length for needing to extend to A for the side in B；
Specific implementation mode
It elaborates with reference to the accompanying drawings and examples to the present invention, wherein finite element software uses MSC.SimXpert
2012。
A kind of board-like satellite structure conode grid rapid generation of the present invention, as shown in Figure 1, it includes following step
Suddenly：
Step 1：The CAD 3D of satellite structure is designed a model and is imported with parasolid neutral file formats
In MSC.SimXpert softwares, then utilize the software middle face abstraction function to selectively honeycomb panel extract in
Face；The menu of the middle face abstraction function used in this step is followed successively by Geometry->Mid-Surface->Automatic,
In in the user interface of pop-up Body type options selection Constant thickness, Delete Solid options choose,
Original three-dimensional entity model will be deleted automatically behind face in extracting in this way；
Step 2：Obtain middle face set to be treated and its element number N：MSC.SimXpert softwares provide secondary open
Hair interface function getPartList allows user to obtain parts list all in model in secondary development program, by
Original three-dimensional entity model is deleted in model, therefore the set of parts list element composition is to need to handle
Middle face geometry, the element number in set is N；
Step 3：To all following operation is carried out in face of (A, B)：
Step 3.1：Calculate the angle theta between the normal vector of A and BAB:Secondary the opening of MSC.SimXpert softwares offer is provided
It is respectively n that hair interface function normal obtains the unit normal vector of A and B respectivelyA=(nAx,nAy,nAz) and nB=(nBx,nBy,
nBz), then the cosine value of angle is cos θ between themAB=nAxnBx+nAynBy+nAznBz, to utilize anti-three in elementary algebra
Angle is θ between angle function can obtain two vectorsAB=arccos (nAxnBx+nAynBy+nAznBz)；
Step 3.2：If 10 ° of ＜ θAB170 ° of ＜, thens follow the steps 3.3 to 3.10；
Step 3.3：Find out in A the side in need extended to B and its number NA：MSC.SimXpert softwares provide secondary
Development interface function getCurves obtains curve list all in a face, and getVerts obtains the vertex list in part,
GetCoordinateAtParameter obtains its rectangular co-ordinate according to the natural coordinates that curve gives, therefore can utilize these
Interface function is by secondary development and realizes in accordance with the following methods：1) all side E of A are obtainedA1,…,EAnWith a vertex of B
VB；2) to each edge EAj(j=1 ..., n), it is assumed that its natural coordinates is expressed as Pj(s) (0≤s≤1), equably obtains Nj+ 1
Point Pjk=Pj(k/Nj) (k=0,1 ..., Nj), wherein Nj=2000；3) d is takenmax=-999999, d=0, then according to non -zero solution
Analyse the formula d in geometryjk=nBx(xjk-xB)+nBy(yjk-yB)+nBz(zjk-zB) each point P is calculated successivelyjkTo B it is oriented away from
From, if | djk| ＞ dmax, then dmax=| djk|, d=djk；Wherein (xjk,yjk,zjk) it is PjkCoordinate in rectangular coordinate system,
(xB,yB,zB) it is VBCoordinate, (nBx,nBy,nBz) be B unit normal vector；4) abd=999999 is taken, to each point PjkIf
d·djk＞ 0 and | djk| ＜ abd, then abd=| djk|；If 5) abd ＜ 0.04, δmax=abd+0.01；Otherwise δmax=0.05；
6)NA=0, to each edge EAj：If to k=0,1 ..., NjHave | djk|≤δmax, then EAjLabeled as 1, NA=NA+1；Otherwise EAj
Labeled as 0；Three kinds of possible cases of this step are respectively such as Fig. 2, and shown in 3,4, face A and face B are non-intersecting in wherein Fig. 2, but on the A of face
The side nearest from face B is more than δ at a distance from the B of facemax, therefore NA=0；Face A is non-intersecting with face B in Fig. 3 and face A on it is nearest from face B
Side at a distance from the B of face be no more than δmax, therefore NA＞ 0；Face A intersects with face B in Fig. 4, therefore no side needs to extend, NA=0；
Step 3.4：If NA＞ 0, thens follow the steps 3.5 to 3.6；It is no to then follow the steps 3.7 to 3.10；
Step 3.5：The side calculated in A needs the length δ extended to BA=δmax/sinθAB, wherein δmaxValue in step
It is obtained in rapid 3.3, θABValue obtained in step 3.1；
Step 3.6：By in A it is in need to B extend side extend δA, as shown in Figure 5：MSC.SimXpert softwares provide
Second development interface ExtendSurfaces allows user by the given length of surface extending, therefore can utilize this function
Interface by secondary development by A it is in need to B extend side extend δA.When calling this function, the value of method parameters is
The value of By Length, m_entities_to_extend parameters is the curve for needing to extend, and is obtained in step 3.3, m_
The value of length parameters is δA；
Step 3.7：Find out in B the side in need extended to A and its number NB：MSC.SimXpert softwares provide secondary
Development interface function getCurves obtains curve list all in a face, and getVerts obtains the vertex list in part,
GetCoordinateAtParameter obtains its rectangular co-ordinate according to the natural coordinates that curve gives, therefore can utilize these
Interface function is by secondary development and realizes in accordance with the following methods：1) all side E of B are obtainedB1,…,EBnWith a vertex of A
VA；2) to each edge EBj(j=1 ..., n) assume that its natural coordinates is expressed as Pj(s) (0≤s≤1), equably obtains Nj+ 1
Point Pjk=Pj(k/Nj) (k=0,1 ..., Nj) wherein Nj=2000；3) d is takenmax=-999999, d=0, then according to non -zero solution
Analyse the formula d in geometryjk=nAx(xjk-xA)+nAy(yjk-yA)+nAz(zjk-zA) each point P is calculated successivelyjkTo A it is oriented away from
From, if | djk| ＞ dmax, then dmax=| djk|, d=djk；Wherein (xjk,yjk,zjk) it is PjkCoordinate in rectangular coordinate system,
(xA,yA,zA) it is VACoordinate, (nAx,nAy,nAz) be A unit normal vector；4) abd=999999 is taken, to each point PjkIf
d·djk＞ 0 and | djk| ＜ abd, then abd=| djk|；If 5) abd ＜ 0.04, δmax=abd+0.01；Otherwise δmax=0.05；
6)NB=0, to each edge EBj：If to k=1 ..., N has | djk|≤δmax, then EBjLabeled as 1, NB=NB+1；Otherwise EBjLabel
It is 0；Three kinds of possible cases of this step are respectively such as Fig. 6, and shown in 7,8, face B and face A are non-intersecting in wherein Fig. 6, but from face on the B of face
Side nearest A is more than δ at a distance from the A of facemax, therefore NB=0；Face B and face A be non-intersecting in Fig. 7 and face B on the side nearest from face A
With at a distance from the A of face be no more than δmax, therefore NB＞ 0；Face B intersects with face A in Fig. 8, therefore no side needs to extend, NB=0；
Step 3.8：If NB＞ 0, thens follow the steps 3.7；It is no to then follow the steps 3.9；
Step 3.9：The side calculated in B needs the length δ extended to AB=δmax/sinθAB, wherein δmaxValue in step
It is obtained in rapid 3.7, θABValue obtained in step 3.1；Although it should be noted that being adopted in this step and step 3.5
With variable δmax, but its value has changed in step 3.7.
Step 3.10：By in B it is in need to A extend side extend δB, as shown in Figure 9：MSC.SimXpert softwares carry
Allow user by the given length of surface extending for second development interface ExtendSurfaces, therefore this letter can be utilized
Number interfaces by secondary development by B institute it is in need to A extend side extension δB.When calling this function, the value of method parameters
Value for By Length, m_entities_to_extend parameters is the curve for needing to extend, and is obtained in step 3.7,
The value of m_length parameters is δB；
Step 3.11：Divide A using B, as shown in Figure 10：MSC.SimXpert softwares provide second development interface
TrimSurfaces allows user to divide curved surface, therefore this function interface can be utilized to be realized by secondary development and utilize B
Divide A.It is B that divisional plane need to be only specified in second development interface function, and the face that is divided is A；
Step 3.12：The face that area is smaller in A is deleted, as shown in figure 11：Secondary using MSC.SimXpert softwares is opened
It sends out interface function DeleteObject and deletes the face that area in A is less than the gross area 5%；
Step 3.13：Divide B using A, as shown in figure 11：MSC.SimXpert softwares provide second development interface
TrimSurfaces allows user to divide curved surface, therefore this function interface can be utilized to be realized by secondary development and utilize A
Divide B.It is A that divisional plane need to be only specified in second development interface function, and the face that is divided is B；
Step 3.14：The face that area is smaller in B is deleted, as shown in figure 12：Secondary using MSC.SimXpert softwares is opened
It sends out interface function DeleteObject and deletes the face that area in B is less than the gross area 5%；
Step 4：All faces are sutured:The curved surface suture function that this step uses in MSC.SimXpert softwares
Menu is followed successively by Geometry->Stitch, wherein Allow disjoint, Treat as in the user interface of pop-up
The purpose that manifold and Force Stitch options are chosen and select all curved surfaces that suture can be realized；
Step 5：Finite element grid is divided to the geometrical model of structure：The net that this step uses in MSC.SimXpert softwares
The menu of lattice partition functionality is followed successively by Meshing->Surface, wherein Mesh type are selected in the user interface of pop-up
Mixed, Mesh method select Auto decided, given size of mesh opening 0.02m simultaneously all curved surfaces to be selected to can be obtained altogether
The finite element grid of node.
By above each step, board-like satellite structure conode grid is realized based on Three Dimensional Design Model and is quickly generated.
Since the grid of generation is continuous, there is no need to additionally introduce rigid coupling unit.In addition, due to being realized using batch processing,
Therefore higher skill, a large amount of time and loaded down with trivial details work are avoided.
Claims (9)
1. a kind of board-like satellite structure conode grid rapid generation, it is characterised in that：It includes the following steps：
Step 1：Face in satellite structure geometrical model extraction to importing finite element modeling software；
Step 2：Obtain middle face set to be treated and its element number N：
Step 3：All are faced, it is assumed that carry out following operation for (A, B)：
Step 3.1：Calculate the angle theta between the normal vector of A and BAB；
Step 3.2：If 10 ° of ＜ θAB170 ° of ＜, thens follow the steps 3.3 to 3.10；
Step 3.3：Find out in A the side in need extended to B and its number NA；
Step 3.4：If NA＞ 0, thens follow the steps 3.5 to 3.6；It is no to then follow the steps 3.7 to 3.10；
Step 3.5：The side calculated in A needs the length δ extended to BA；
Step 3.6：By in A it is in need to B extend side extend δA；
Step 3.7：Find out in B the side in need extended to A and its number NB；
Step 3.8：If NB＞ 0, thens follow the steps 3.7；It is no to then follow the steps 3.9；
Step 3.9：The side calculated in B needs the length δ extended to AB；
Step 3.10：By in B it is in need to A extend side extend δB；
Step 3.11：Divide A using B；
Step 3.12：Delete the face that area is smaller in A；
Step 3.13：Divide B using A；
Step 3.14：Delete the face that area is smaller in B；
Step 4：All faces are sutured；
Step 5：Finite element grid is divided to the geometrical model of structure；
The method of described " dividing A using B " is：Existing finite element software all provides second development interface and enables a user to out
The function of sending out new, therefore realized by secondary development using relevant function interface and divide A using B, in second development interface letter
It is B that divisional plane need to be only specified in number, and the face that is divided is A, and required second development interface function is by inquiring finite element
The dependent help document of software obtains；
The method of described " dividing B using A " is：Existing finite element software generally all provides second development interface and makes user
New function is developed, therefore is realized by secondary development using relevant function interface and divides B using A；In second development interface
It is A that divisional plane need to be only specified in function, and the face that is divided is B, and required second development interface function is limited by inquiring
The dependent help document of meta software obtains.
2. a kind of board-like satellite structure conode grid rapid generation according to claim 1, it is characterised in that：
Described in step 3.3 " find out in A the side in need extended to B and its number NA" method be：Existing finite element is soft
Part all provides second development interface and enables a user to develop new function, therefore passes through secondary development using relevant function interface
And it realizes in accordance with the following methods：1) all side E of A are obtainedA1,…,EAnWith a vertex V of BB；2) to each edge EAj(j=
1 ..., n), it is assumed that its natural coordinates is expressed as Pj(s) (0≤s≤1),
Equably obtain Nj+ 1 point, Pjk=Pj(k/Nj) (k=0,1 ..., Nj),
Wherein Nj≥1000；3) d is takenmax=-999999, d=0,
Then according to the formula d in elementary analytic geometryjk=nBx(xjk-xB)+nBy(yjk-yB)+nBz(zjk-zB) calculate successively often
A point PjkTo the directed distance of B, if | djk| ＞ dmax, then dmax=| djk|, d=djk；
Wherein (xjk,yjk,zjk) it is PjkCoordinate in rectangular coordinate system, (xB,yB,zB) it is VBCoordinate, (nBx,nBy,nBz) be
The unit normal vector of B；4) abd=999999 is taken, to each point PjkIf ddjk＞ 0 and | djk| ＜ abd, then abd=| djk|；
If 5) abd ＜ 0.04,
Then δmax=abd+0.01；Otherwise δmax=0.05；6)NA=0, to each edge EAj：If to k=0,1 ..., NjHave | djk|
≤δmax, then EAjLabeled as 1, NA=NA+1；Otherwise EAjLabeled as 0；The vertex for obtaining all sides and B A in this step is logical
The interface for crossing software is realized, and side and curve can be indicated using natural coordinates in finite element software, gives its nature
Its coordinate representation under rectangular coordinate system is obtained by interface after coordinate；Required second development interface function in this step
Dependent help document by inquiring finite element software obtains.
3. a kind of board-like satellite structure conode grid rapid generation according to claim 1, it is characterised in that：
Described in step 3.5 " calculate A in side need the length δ extended to BA" method be：δA=δmax/sinθAB, wherein δmax
Value obtained in step 3.3, θABValue obtained in step 3.1.
4. a kind of board-like satellite structure conode grid rapid generation according to claim 1, it is characterised in that：
Described in step 3.6 " by A it is in need to B extend side extend δA" method be：Existing finite element software all carries
It enables a user to develop new function for second development interface, therefore will be in A by secondary development using relevant function interface
It is in need to B extend side extend δA；Each side in one face is carried out extending required second development interface function
Dependent help document by inquiring finite element software obtains.
5. a kind of board-like satellite structure conode grid rapid generation according to claim 1, it is characterised in that：
Described in step 3.7 " find out in B the side in need extended to A and its number NB" method be：Existing finite element is soft
Part all provides second development interface and enables a user to develop new function, therefore passes through secondary development using relevant function interface
And it realizes in accordance with the following methods：1) all side E of B are obtainedB1,…,EBnWith a vertex V of AA；2) to each edge EBj(j=
1 ..., n) assume that its natural coordinates is expressed as Pj(s) (0≤s≤1),
Equably obtain Nj+ 1 point, Pjk=Pj(k/Nj) (k=0,1 ..., Nj),
Wherein Nj≥1000；3) d is takenmax=-999999, d=0,
Then according to the formula d in elementary analytic geometryjk=nAx(xjk-xA)+nAy(yjk-yA)+nAz(zjk-zA) calculate successively often
A point PjkTo the directed distance of A, if | djk| ＞ dmax, then dmax=| djk|, d=djk；
Wherein (xjk,yjk,zjk) it is PjkCoordinate in rectangular coordinate system, (xA,yA,zA) it is VACoordinate, (nAx,nAy,nAz) be
The unit normal vector of A；4) abd=999999 is taken, to each point PjkIf ddjk＞ 0 and | djk| ＜ abd, then abd=| djk|；
If 5) abd ＜ 0.04,
Then δmax=abd+0.01；Otherwise δmax=0.05；6)NB=0, to each edge EBj：If to k=0,1 ..., NjHave | djk|
≤δmax, then EBjLabeled as 1, NB=NB+1；Otherwise EBjLabeled as 0；The vertex for obtaining all sides and A B in this step is logical
The interface for crossing software is realized, and side or curve can be indicated using natural coordinates in finite element software, gives its nature
Its coordinate representation under rectangular coordinate system is obtained by interface after coordinate；Required second development interface function in this step
Dependent help document by inquiring finite element software obtains.
6. a kind of board-like satellite structure conode grid rapid generation according to claim 1, it is characterised in that：
Described in step 3.9 " calculate B in side need the length δ extended to AB" method be：δB=δmax/sinθAB, wherein δmax
Value obtained in step 3.7, θABValue obtained in step 3.1；Although should be noted that this step and step
Variable δ has been all made of in rapid 3.5max, but its value has changed in step 3.7.
7. a kind of board-like satellite structure conode grid rapid generation according to claim 1, it is characterised in that：
Described in step 3.10 " by B it is in need to A extend side extend δB" method be：Existing finite element software all carries
It enables a user to develop new function for second development interface, therefore will be in B by secondary development using relevant function interface
It is in need to A extend side extend δB；Each side in one face is carried out extending required second development interface function
Dependent help document by inquiring finite element software obtains.
8. a kind of board-like satellite structure conode grid rapid generation according to claim 1, it is characterised in that：
The method of " deleting the smaller face of area in A " described in step 3.12 is：It is obtained using the second development interface of finite element software
To the area in each face in A, the gross area in all faces in A is just calculated by simple addition in this way；If the area in some face
It is less than given threshold gamma with the ratio between the gross area, then is deleted using second development interface function.
9. a kind of board-like satellite structure conode grid rapid generation according to claim 1, it is characterised in that：
The method of " deleting the smaller face of area in B " described in step 3.14 is：Utilize the second development interface energy of finite element software
The area in each face in B is obtained, just calculates the gross area in all faces in B by simple addition in this way；If the face in some face
The ratio between product and the gross area are less than given threshold gamma, then are deleted using second development interface function.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610005236.3A CN105677965B (en) | 2016-01-04 | 2016-01-04 | A kind of board-like satellite structure conode grid rapid generation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610005236.3A CN105677965B (en) | 2016-01-04 | 2016-01-04 | A kind of board-like satellite structure conode grid rapid generation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105677965A CN105677965A (en) | 2016-06-15 |
| CN105677965B true CN105677965B (en) | 2018-09-18 |
Family
ID=56298946
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610005236.3A Expired - Fee Related CN105677965B (en) | 2016-01-04 | 2016-01-04 | A kind of board-like satellite structure conode grid rapid generation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105677965B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107066676A (en) * | 2017-01-19 | 2017-08-18 | 北京航空航天大学 | A kind of finite element automation modeling method based on satellite plate and shell structure |
| CN110727989B (en) * | 2018-06-28 | 2023-08-04 | 中车大同电力机车有限公司 | Structural fatigue strength analysis method, device and computer readable storage medium |
| CN110263426B (en) * | 2019-06-18 | 2020-12-04 | 珠海格力电器股份有限公司 | Medium surface extraction method and device of thin-wall structure model and storage medium |
| CN115828662B (en) * | 2022-11-04 | 2024-10-01 | 广州中望龙腾软件股份有限公司 | Finite element simulation solving method, system, equipment and medium based on grid division |
| CN115797600B (en) * | 2023-02-08 | 2023-06-16 | 安世亚太科技股份有限公司 | Three-dimensional geometric model wrapping grid generation system |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1386200A (en) * | 2000-06-29 | 2002-12-18 | 目标储油层公司 | Feature modeling method in a finite element model |
| CN101882168A (en) * | 2009-05-06 | 2010-11-10 | 利弗莫尔软件技术公司 | Be used at finite element analysis simulating beam and the surperficial method and system that contacts |
| CN102306396A (en) * | 2011-09-15 | 2012-01-04 | 山东大学 | Three-dimensional entity model surface finite element mesh automatic generation method |
| US8180605B1 (en) * | 2005-11-23 | 2012-05-15 | Livermore Software Technology Corporation | Methods and systems for creating a smooth contact-impact interface in finite element analysis |
| CN104408267A (en) * | 2014-12-17 | 2015-03-11 | 湖北工业大学 | Finite element modeling method for special-shaped spatial structure |
-
2016
- 2016-01-04 CN CN201610005236.3A patent/CN105677965B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1386200A (en) * | 2000-06-29 | 2002-12-18 | 目标储油层公司 | Feature modeling method in a finite element model |
| US8180605B1 (en) * | 2005-11-23 | 2012-05-15 | Livermore Software Technology Corporation | Methods and systems for creating a smooth contact-impact interface in finite element analysis |
| CN101882168A (en) * | 2009-05-06 | 2010-11-10 | 利弗莫尔软件技术公司 | Be used at finite element analysis simulating beam and the surperficial method and system that contacts |
| CN102306396A (en) * | 2011-09-15 | 2012-01-04 | 山东大学 | Three-dimensional entity model surface finite element mesh automatic generation method |
| CN104408267A (en) * | 2014-12-17 | 2015-03-11 | 湖北工业大学 | Finite element modeling method for special-shaped spatial structure |
Non-Patent Citations (2)
| Title |
|---|
| 基于三维最小二乘方法的空间直线度误差评定;王炳杰等;《北京航空航天大学学报》;20141031;第40卷(第10期);第1477-1480页 * |
| 多工况线性结构稳健拓扑优化设计;付志方等;《力学学报》;20150731;第47卷(第4期);第642-650页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105677965A (en) | 2016-06-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105677965B (en) | A kind of board-like satellite structure conode grid rapid generation | |
| Natarajan et al. | Integrating strong and weak discontinuities without integration subcells and example applications in an XFEM/GFEM framework | |
| Hoffman et al. | Decomposition plans for geometric constraint problems, part II: New algorithms | |
| CN100595768C (en) | System and method for mesh and body hybrid modeling using 3d scan data | |
| Tang et al. | A hybrid geometric modeling method for lattice structures fabricated by additive manufacturing | |
| Lee et al. | On error estimation and adaptive refinement for element free Galerkin method: Part I: stress recovery and a posteriori error estimation | |
| CN108876923A (en) | A kind of three-dimensional pore space Scale Model method for reconstructing based on the micro- CT image of rock | |
| Liu et al. | Topology optimization for hybrid additive-subtractive manufacturing | |
| CN104331933B (en) | A kind of fabrication orientation self adaptation fast selecting method | |
| US10926469B2 (en) | Three-dimensional laminating and fabricating system, three-dimensional laminating and fabricating method, laminating and fabricating control apparatus and method of controlling the same, and control program | |
| JPH0816629A (en) | Mesh generating method/device for analysis | |
| Huang et al. | A high performance multi-bit-width booth vector systolic accelerator for NAS optimized deep learning neural networks | |
| CN104123411B (en) | A kind of succinct derivation method of powertrain system of automobile rotary inertia synthesis | |
| Larwood et al. | Domain decomposition approach for parallel unstructured mesh generation | |
| Kumar et al. | Geometry based triangulation of multiple trimmed NURBS surfaces | |
| CN105701275B (en) | A kind of beam slab combined type satellite structure conode grid rapid generation | |
| Merkley et al. | Methods and applications of generalized sheet insertion for hexahedral meshing | |
| CN111027274B (en) | Three-dimensional chip layout method | |
| CN108389154B (en) | Method for implementing cutting area segmentation technology for parallel drawing in GPU | |
| Yang et al. | Feature extraction and volume decomposition for orthogonal layered manufacturing | |
| Kallis et al. | Incremental 2D Delaunay triangulation core implementation on FPGA for surface reconstruction via high-level synthesis | |
| Matsunaga et al. | Timing-constrained area minimization algorithm for parallel prefix adders | |
| CN101625708A (en) | Visual structure analysis system | |
| Dong et al. | The hybrid stress model for Mindlin-Reissner plates based on a stress optimization condition | |
| Ranaivomiarana et al. | Optimal orthotropy and density distribution of two-dimensional structures |
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: 20180918 Termination date: 20200104 |