CN105677965A - Method for rapidly generating plate-type satellite-structure node merging grids - Google Patents
Method for rapidly generating plate-type satellite-structure node merging grids Download PDFInfo
- Publication number
- CN105677965A CN105677965A CN201610005236.3A CN201610005236A CN105677965A CN 105677965 A CN105677965 A CN 105677965A CN 201610005236 A CN201610005236 A CN 201610005236A CN 105677965 A CN105677965 A CN 105677965A
- Authority
- CN
- China
- Prior art keywords
- finite element
- limit
- max
- face
- interface
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
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
Abstract
The invention discloses a method for rapidly generating plate-type satellite-structure node merging grids. The method includes the following steps that 1, middle surfaces of a satellite-structure geometric model introduced into finite element modeling software are extracted; 2, to-be-processed middle surface sets and the number N of elements of the middle surface sets are obtained; 3, related operation is carried out on all surface pairs (A, B); 4, all the surfaces are sutured; 5, the geometric model of the structure is divided into finite element grids. According to the method, the geometry obtained in the mode that the middle surfaces of the satellite-structure three-dimensional solid model are extracted can be automatically processed, gaps between the faces are eliminated, redundant parts are deleted, and the high-quality node-merging finite element grids can be divided through the finally-generated geometric model. A large amount of human-computer interaction operation is avoided, required time is short, universality is high, the application range is wide, and the quite-high practical engineering value is achieved.
Description
Technical field:
The present invention relates to a kind of board-like satellite structure conode grid rapid generation, belong to the technical field of structural analysis of spacecraft.
Background technology:
Current satellite structure many employings honeycomb sandwich panel bears and transmits mechanical loading, and provides installation to support for instrument and equipment. In order to whether the design of verifying satellites structure meets requirement, engineering is generally adopted Finite Element Method it is carried out mechanics property analysis. The geometric model of satellite structure is generally adopted 3D solid and is modeled, but is generally adopted plate and shell structure unit when finite element modeling and carries out mechanical simulation. If dividing finite element grid by 3D solid in directly being extracted behind face, then owing to the gap between adjacent surface can cause the discontinuous of finite element grid; If adopting rigid coupling unit to carry out unit connection, then requiring that the grid comparison match on adjacent surface is beneficial to create rigid coupling unit on the one hand, structure can be introduced additional rigidity by creating of rigid coupling unit on the other hand; In addition, the establishment of rigid coupling unit also can increase the amount of calculation of finite element analysis, especially for the dynamic analysis situation such as common model analysis and frequency response analysis in satellite structure analysis. Therefore engineering wishes to generate the finite element grid of conode. The Major Difficulties that satellite structure conode grid quickly generates is: due to the existence of honeycomb interlayer plate thickness, in the geometric model obtained, there is gap, it is impossible to directly generate finite element grid between adjacent surface extracting from three-dimensional entity model behind face. In order to realize the division of conode grid, it is necessary to first opposite carries out extending or cutting, thus the gap eliminated between face and face or unnecessary part. Owing to satellite structure generally having tens structure plates, therefore rely only on the geometric manipulations function of business finite element software self and realize the generation of whole satellite structure conode finite element grid by skill higher for needs, substantial amounts of time and loaded down with trivial details work by man-machine interaction mode.
Summary of the invention:
1, purpose: in order to realize efficiently producing of satellite structure FEM (finite element) model, it is an object of the invention to provide a kind of board-like satellite structure conode grid rapid generation.The method obtains the middle face geometric model of satellite structure plate first with the middle face abstraction function of finite element software, then pass through faced by traversal obtains and need to carry out processing, secondly often group is automatically processed in the face of being sequentially carried out, thus obtaining the geometric model of whole satellite structure, the geometry again with finite element software sews up functional realiey by as a whole for the geometric model stitching of satellite structure. Carry out stress and strain model on this basis, obtain the high-quality conode finite element grid being simulated based on Shell Finite Element. Due to this method realize automatically process, therefore avoid substantial amounts of man-machine interactive operation, it is necessary to time few, highly versatile, applied widely.
2, technical scheme: the technical solution adopted for 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 comprises the following steps:
Step 1: face during the satellite structure geometric model importing finite element modeling software is extracted;
Step 2: obtain and need middle set to be processed and element number N thereof:
Step 3: carry out following operation in the face of (A, B) to all of:
Step 3.1: calculate the angle theta between the normal vector of A and BAB;
Step 3.2: if 10 ° of < θAB< 170 °, then perform step 3.3 to 3.10;
Step 3.3: find out and be there is a need in A to the B limit extended and number N thereofA;
Step 3.4: if NA> 0, then perform step 3.5 to 3.6; Otherwise perform step 3.7 to 3.10;
Step 3.5: calculating the limit in A needs to the B length δ extendedA;
Step 3.6: extended δ to the B limit extended by there is a need in AA;
Step 3.7: find out and be there is a need in B to the A limit extended and number N thereofB;
Step 3.8: if NB> 0, then perform step 3.7; Otherwise perform step 3.9;
Step 3.9: calculating the limit in B needs to the A length δ extendedB;
Step 3.10: extended δ to the A limit extended by there is a need in BB;
Step 3.11: utilize B to split A;
Step 3.12: delete the face that in A, area is less;
Step 3.13: utilize A to split B;
Step 3.14: delete the face that in B, area is less;
Step 4: all of face is sewed up;
Step 5: the geometric model of structure is divided finite element grid;
Wherein, the method in " during the satellite structure geometric model importing finite element modeling software is extracted face " described in step 1 is: existing finite element software generally all supports the middle face abstraction function of thin-walled parts, and satellite cellular Sandwich Plates is typically all in uniform thickness, therefore selects the geometric model of honeycomb sandwich construction plate in dummy satellite to carry out middle extraction either directly through the user interface of finite element software and can obtain corresponding middle face geometric model;
Wherein, the method of " obtain and need middle set to be processed and element number N thereof " described in step 2 is: existing finite element software generally all provides second development interface to make user can develop new function, and relevant function interface therefore can be utilized to obtain all middle set formed and its element number in model by secondary development. The second development interface function that in acquisition, face set is required can be obtained by the dependent help document of inquiry finite element software;
Wherein, " the angle theta between the normal vector of calculating A and B described in step 3.1AB" method be: assume the unit normal vector respectively n of A and BA=(nAx,nAy,nAz) and nB=(nBx,nBy,nBz), then between them, the cosine value of angle is cos θAB=nAxnBx+nAynBy+nAznBz, thus utilizing the antitrigonometric function in elementary algebra can obtain angle between two vectors is θAB=arccos (nAxnBx+nAynBy+nAznBz);
Wherein, " finding out and being there is a need in A to the B limit extended and number N thereof described in step 3.3A" method be: existing finite element software generally all provides second development interface to make user can develop new function, therefore can utilize relevant function interface by secondary development and in accordance with the following methods realize: 1) acquisition all of limit E of AA1,…,EAnWith B summit VB; 2) to each edge EAj(j=1 ..., n), it is assumed that its natural coordinates is expressed as PjS () (0≤s≤1), obtains N equablyj+ 1 some 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 analytical geometryjk=nBx(xjk-xB)+nBy(yjk-yB)+nBz(zjk-zB) calculate each some P successivelyjkTo the directed distance of B, if | djk| > dmax, then dmax=| djk|, d=djk; Wherein (xjk,yjk,zjk) for PjkCoordinate in rectangular coordinate system, (xB,yB,zB) for VBCoordinate, (nBx,nBy,nBz) for the unit normal vector of B; 4) abd=999999 is taken, to each some PjkIf, d djk> 0 and | djk| < abd, then abd=| djk|; 5) if abd < 0.04, then δmax=abd+0.01; Otherwise δmax=0.05; 6) NA=0, to each edge EAjIf: to k=0,1 ..., NjAll have | djk|≤δmax, then EAjIt is labeled as 1, NA=NA+ 1; Otherwise EAjIt is labeled as 0; The summit obtaining all of limit of A and B in this step all can be realized by the interface of software, and limit or curve generally can adopt natural coordinates to be indicated in finite element software, after its natural coordinates given, its coordinate representation under rectangular coordinate system can be obtained by interface. Second development interface function required in this step can be obtained by the dependent help document of inquiry finite element software;
Wherein, " calculating the limit in A and needing the length δ to B extension described in step 3.5A" method be: δA=δmax/sinθAB, wherein δmaxValue obtain in step 3.3, θABValue obtain in step 3.1;
Wherein, " will be there is a need in A to the B limit extended extension δ described in step 3.6A" method be: existing finite element software generally all provides second development interface to make user can develop new function, therefore can utilize relevant function interface by secondary development will A be there is a need to B extend limit extension δA. Some limits in one face carry out extend required second development interface function to be obtained by the dependent help document of inquiry finite element software;
Wherein, " finding out and being there is a need in B to the A limit extended and number N thereof described in step 3.7B" method be: existing finite element software generally all provides second development interface to make user can develop new function, therefore can utilize relevant function interface by secondary development and in accordance with the following methods realize: 1) acquisition all of limit E of BB1,…,EBnWith A summit VA; 2) to each edge EBj(j=1 ..., n) assume that its natural coordinates is expressed as PjS () (0≤s≤1), obtains N equablyj+ 1 some 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 analytical geometryjk=nAx(xjk-xA)+nAy(yjk-yA)+nAz(zjk-zA) calculate each some P successivelyjkTo the directed distance of A, if | djk| > dmax, then dmax=| djk|, d=djk; Wherein (xjk,yjk,zjk) for PjkCoordinate in rectangular coordinate system, (xA,yA,zA) for VACoordinate, (nAx,nAy,nAz) for the unit normal vector of A; 4) abd=999999 is taken, to each some PjkIf, d djk> 0 and | djk| < abd, then abd=| djk|;5) if abd < 0.04, then δmax=abd+0.01; Otherwise δmax=0.05; 6) NB=0, to each edge EBjIf: to k=1 ..., N all has | djk|≤δmax, then EBjIt is labeled as 1, NB=NB+ 1; Otherwise EBjIt is labeled as 0; The summit obtaining all of limit of B and A in this step all can be realized by the interface of software, and limit or curve generally can adopt natural coordinates to be indicated in finite element software, after its natural coordinates given, its coordinate representation under rectangular coordinate system can be obtained by interface. Second development interface function required in this step can be obtained by the dependent help document of inquiry finite element software;
Wherein, " calculating the limit in B and needing the length δ to A extension described in step 3.9B" method be: δB=δmax/sinθAB, wherein δmaxValue obtain in step 3.7, θABValue obtain in step 3.1; Although it should be noted that and this step and step 3.5 all have employed variable δmax, but in step 3.7, its value is altered.
Wherein, " will be there is a need in B to the A limit extended extension δ described in step 3.10B" method be: existing finite element software generally all provides second development interface to make user can develop new function, therefore can utilize relevant function interface by secondary development will B be there is a need to A extend limit extension δB. Some limits in one face carry out extend required second development interface function to be obtained by the dependent help document of inquiry finite element software;
Wherein, the method of " utilizing B to split A " described in step 3.11 is: existing finite element software generally all provides second development interface to make user can develop new function, and relevant function interface therefore can be utilized to realize utilizing B to split A by secondary development. Only need to specify divisional plane in second development interface function is B, and divided is A, and required second development interface function can be obtained by the dependent help document of inquiry finite element software;
Wherein, the method for " the deleting the face that in A, area is less " described in step 3.12 is: utilizes the second development interface of finite element software can obtain the area of each in A, so just can calculate the gross area of all in A by simple addition; If the ratio of the area in certain face and the gross area is less than given threshold gamma, then second development interface function is utilized to be deleted;
Wherein, the method of " utilizing A to split B " described in step 3.13 is: existing finite element software generally all provides second development interface to make user can develop new function, and relevant function interface therefore can be utilized to realize utilizing A to split B by secondary development; Only need to specify divisional plane in second development interface function is A, and divided is B, and required second development interface function can be obtained by the dependent help document of inquiry finite element software;
Wherein, the method for " the deleting the face that in B, area is less " described in step 3.14 is: utilizes the second development interface of finite element software can obtain the area of each in B, so just can calculate the gross area of all in B by simple addition; If the ratio of the area in certain face and the gross area is less than given threshold gamma, then second development interface function is utilized to be deleted;
Wherein, the method of " all of face is sewed up " described in step 4 is: existing finite element software is generally all supported to sew up as a whole by several adjacent faces, adjacent limit, two such face is changing to share limit, automatically will generate node during stress and strain model on shared limit, and the unit between adjacent surface shares node. Therefore sew up function either directly through the curved surface of finite element software, select the face of there is a need to grid division, then carry out curved surface and sew up;
Wherein, the method of " geometric model of structure is divided finite element grid " described in steps of 5 is: existing finite element software generally all supports grid division on curved surface, and user only needs designating unit type, size of mesh opening selection to need the curved surface of grid division can realize the automatic division of grid.Therefore either directly through the stress and strain model function of finite element software, and designating unit type is that in triangle or tetragon, given size of mesh opening preference pattern, all of curved surface can carry out stress and strain model. Owing to sewing up as a whole by all of curved surface in step 4, the grid therefore obtained in this step will be continuous print, and the grid between adjacent surface is conode.
By above each step, achieve board-like satellite structure conode grid to quickly generate, and in the generation process of grid, create new model based on the three-dimensional entity model of satellite structure without being specifically designed for mechanical analysis, the grid generated is continuous print, from without additionally introducing rigid coupling unit. Further, since adopt batch processing to realize, therefore avoid higher skill, substantial amounts of time and loaded down with trivial details work.
3, advantage and effect: the invention has the beneficial effects as follows: provide a kind of board-like satellite structure conode grid rapid generation, the geometry that in satellite structure three-dimensional entity model can being extracted by the method, face obtains automatically processes, gap between elimination face and face and delete unnecessary part, the geometric model ultimately generated can divide high-quality conode finite element grid. Due to this method realize automatically process, therefore avoid substantial amounts of man-machine interactive operation, it is necessary to time few, highly versatile, applied widely. Due to the finite element grid conode that this method generates, therefore need not create extra rigid coupling unit, it is to avoid introducing and extra amount of calculation of additional stiffness. Further, since board-like satellite structure is widely used, therefore the present invention has very strong engineering practical value in satellite structure finite element modeling.
Accompanying drawing illustrates a kind of board-like satellite structure conode grid rapid generation
Fig. 1. the method for the invention flow chart.
Fig. 2. face A is non-intersect with the plane at B place, face, and the minimum range of the point on the A of face and face B is more than δmax。
Fig. 3. face A is non-intersect with the plane at B place, face, and the minimum range of the point on the A of face and face B is less than δmax。
Fig. 4. face A intersects with the plane at B place, face.
Fig. 5. face B is non-intersect with the plane at A place, face, and the minimum range of the point on the B of face and face A is more than δmax。
Fig. 6. face B is non-intersect with the plane at A place, face, and the minimum range of the point on the B of face and face A is less than δmax。
Fig. 7. face B intersects with the plane at A place, face.
Fig. 8. extended surface A needs the limit extended.
Fig. 9. extended surface B needs the limit extended.
Figure 10. utilize face B divisional plane A.
Figure 11. delete the face that in A, area is only small.
Figure 12. delete the face that in B, area is only small.
In figure, symbol description is as follows:
I, j be two enumerators respectively;
nA,nBThe respectively normal vector of face A and face B;
θABFor nA,nBBetween angle;
N needs the number of middle to be processed in model;
NA,NBFace A and face B respectively needs the number on the limit extended;
δmaxFor the distance threshold calculated or give;
δANeed to the B length extended for the limit in A;
δBNeed to the A length extended for the limit in B;
Detailed description of the invention
Below in conjunction with drawings and Examples, the present invention being elaborated, wherein finite element software adopts MSC.SimXpert2012.
A kind of board-like satellite structure conode grid rapid generation of the present invention, as shown in Figure 1, it comprises the following steps:
Step 1: the CAD 3D of satellite structure is designed a model and imports in MSC.SimXpert software with parasolid neutral file format, then utilize the middle face abstraction function of this software 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, wherein in the user interface ejected, Bodytype option selects Constantthickness, DeleteSolid option is chosen, and will automatically be deleted by original three-dimensional entity model in so extracting behind face;
Step 2: obtain and need middle set to be processed and element number N:MSC.SimXpert software thereof to provide second development interface function getPartList to make user can obtain all of parts list in model in secondary development program, owing to model having deleted original three-dimensional entity model, therefore the set that this parts list is elementary composition is needs middle geometry to be processed, and the element number in set is N;
Step 3: carry out following operation in the face of (A, B) to all of:
Step 3.1: calculate the angle theta between the normal vector of A and BAB: utilize the second development interface function normal that MSC.SimXpert software provides to obtain the unit normal vector respectively n of A and B respectivelyA=(nAx,nAy,nAz) and nB=(nBx,nBy,nBz), then between them, the cosine value of angle is cos θAB=nAxnBx+nAynBy+nAznBz, thus utilizing the antitrigonometric function in elementary algebra can obtain angle between two vectors is θAB=arccos (nAxnBx+nAynBy+nAznBz);
Step 3.2: if 10 ° of < θAB< 170 °, then perform step 3.3 to 3.10;
Step 3.3: find out and be there is a need in A to the B limit extended and number N thereofA: MSC.SimXpert software provides second development interface function getCurves to obtain all of curve list in a face, getVerts obtains the vertex list in part, getCoordinateAtParameter obtains its rectangular coordinate according to the natural coordinates that curve is given, these interface functions therefore can be utilized by secondary development and realize in accordance with the following methods: 1) obtain all of limit E of AA1,…,EAnWith B summit VB; 2) to each edge EAj(j=1 ..., n), it is assumed that its natural coordinates is expressed as PjS () (0≤s≤1), obtains N equablyj+ 1 some Pjk=Pj(k/Nj) (k=0,1 ..., Nj), wherein Nj=2000; 3) d is takenmax=-999999, d=0, then according to the formula d in elementary analytical geometryjk=nBx(xjk-xB)+nBy(yjk-yB)+nBz(zjk-zB) calculate each some P successivelyjkTo the directed distance of B, if | djk| > dmax, then dmax=| djk|, d=djk; Wherein (xjk,yjk,zjk) for PjkCoordinate in rectangular coordinate system, (xB,yB,zB) for VBCoordinate, (nBx,nBy,nBz) for the unit normal vector of B; 4) abd=999999 is taken, to each some PjkIf, d djk> 0 and | djk| < abd, then abd=| djk|; 5) if abd < 0.04, then δmax=abd+0.01; Otherwise δmax=0.05; 6) NA=0, to each edge EAjIf: to k=0,1 ..., NjAll have | djk|≤δmax, then EAjIt is labeled as 1, NA=NA+ 1; Otherwise EAjIt is labeled as 0; Three kinds of possible cases of this step are respectively as shown in Fig. 2,3,4, and wherein in Fig. 2, face A and face B is non-intersect, but from the distance of limit nearest for face B and face B more than δ on the A of facemax, therefore NA=0; In Fig. 3, face A and face B be non-intersect and the A upper distance from limit nearest for face B and face B in face is less than δmax, therefore NA> 0; In Fig. 4, face A intersects with face B, and so there is no limit needs to extend, NA=0;
Step 3.4: if NA> 0, then perform step 3.5 to 3.6; Otherwise perform step 3.7 to 3.10;
Step 3.5: calculating the limit in A needs to the B length δ extendedA=δmax/sinθAB, wherein δmaxValue obtain in step 3.3, θABValue obtain in step 3.1;
Step 3.6: extended δ to the B limit extended by there is a need in AA, as shown in Figure 5: MSC.SimXpert software provides second development interface ExtendSurfaces to make user by length given for surface extending, this function interface therefore can be utilized to be extended δ by there is a need in A to the B limit extended by secondary developmentA.When calling this function, the value of method parameter is the value of ByLength, m_entities_to_extend parameter is the curve needing to extend, and obtains in step 3.3, and the value of m_length parameter is δA;
Step 3.7: find out and be there is a need in B to the A limit extended and number N thereofB: MSC.SimXpert software provides second development interface function getCurves to obtain all of curve list in a face, getVerts obtains the vertex list in part, getCoordinateAtParameter obtains its rectangular coordinate according to the natural coordinates that curve is given, these interface functions therefore can be utilized by secondary development and realize in accordance with the following methods: 1) obtain all of limit E of BB1,…,EBnWith A summit VA; 2) to each edge EBj(j=1 ..., n) assume that its natural coordinates is expressed as PjS () (0≤s≤1), obtains N equablyj+ 1 some Pjk=Pj(k/Nj) (k=0,1 ..., Nj) wherein Nj=2000; 3) d is takenmax=-999999, d=0, then according to the formula d in elementary analytical geometryjk=nAx(xjk-xA)+nAy(yjk-yA)+nAz(zjk-zA) calculate each some P successivelyjkTo the directed distance of A, if | djk| > dmax, then dmax=| djk|, d=djk; Wherein (xjk,yjk,zjk) for PjkCoordinate in rectangular coordinate system, (xA,yA,zA) for VACoordinate, (nAx,nAy,nAz) for the unit normal vector of A; 4) abd=999999 is taken, to each some PjkIf, d djk> 0 and | djk| < abd, then abd=| djk|; 5) if abd < 0.04, then δmax=abd+0.01; Otherwise δmax=0.05; 6) NB=0, to each edge EBjIf: to k=1 ..., N all has | djk|≤δmax, then EBjIt is labeled as 1, NB=NB+ 1; Otherwise EBjIt is labeled as 0; Three kinds of possible cases of this step are respectively as shown in Fig. 6,7,8, and wherein in Fig. 6, face B and face A is non-intersect, but from the distance of limit nearest for face A and face A more than δ on the B of facemax, therefore NB=0; In Fig. 7, face B and face A be non-intersect and the B upper distance from limit nearest for face A and face A in face is less than δmax, therefore NB> 0; In Fig. 8, face B intersects with face A, and so there is no limit needs to extend, NB=0;
Step 3.8: if NB> 0, then perform step 3.7; Otherwise perform step 3.9;
Step 3.9: calculating the limit in B needs to the A length δ extendedB=δmax/sinθAB, wherein δmaxValue obtain in step 3.7, θABValue obtain in step 3.1; Although it should be noted that and this step and step 3.5 all have employed variable δmax, but in step 3.7, its value is altered.
Step 3.10: extended δ to the A limit extended by there is a need in BB, as shown in Figure 9: MSC.SimXpert software provides second development interface ExtendSurfaces to make user by length given for surface extending, this function interface therefore can be utilized to be extended δ by there is a need in B to the A limit extended by secondary developmentB. When calling this function, the value of method parameter is the value of ByLength, m_entities_to_extend parameter is the curve needing to extend, and obtains in step 3.7, and the value of m_length parameter is δB;
Step 3.11: utilize B to split A, as shown in Figure 10: MSC.SimXpert software provides second development interface TrimSurfaces to make user can split curved surface, therefore can utilize this function interface to realize utilizing B to split A by secondary development. Only need to specify divisional plane in second development interface function is B, and divided is A;
Step 3.12: delete the face that area is less in A, as shown in figure 11: utilize the second development interface function DeleteObject of MSC.SimXpert software to delete in A area less than the face of the gross area 5%;
Step 3.13: utilize A to split B, as shown in figure 11: MSC.SimXpert software provides second development interface TrimSurfaces to make user can split curved surface, therefore can utilize this function interface to realize utilizing A to split B by secondary development. Only need to specify divisional plane in second development interface function is A, and divided is B;
Step 3.14: delete the face that area is less in B, as shown in figure 12: utilize the second development interface function DeleteObject of MSC.SimXpert software to delete in B area less than the face of the gross area 5%;
Step 4: all of face is sewed up: the curved surface that in MSC.SimXpert software, this step uses is sewed up the menu of function and is followed successively by Geometry-> Stitch, wherein in the user interface ejected, Allowdisjoint, Treatasmanifold and ForceStitch option all chooses and selects all of curved surface can realize the purpose sewed up;
Step 5: the geometric model of structure is divided finite element grid: in MSC.SimXpert software, the menu of the stress and strain model function that this step uses is followed successively by Meshing-> Surface, wherein in the user interface ejected, Meshtype selects Mixed, Meshmethod selects Autodecided, and given size of mesh opening 0.02m also selects all of curved surface can obtain the finite element grid of conode.
By above each step, achieve board-like satellite structure conode grid based on Three Dimensional Design Model and quickly generate. Grid owing to generating is continuous print, therefore without additionally introducing rigid coupling unit. Further, since adopt batch processing to realize, therefore avoid higher skill, substantial amounts of time and loaded down with trivial details work.
Claims (10)
1. a board-like satellite structure conode grid rapid generation, it is characterised in that: it comprises the following steps:
Step 1: face during the satellite structure geometric model importing finite element modeling software is extracted;
Step 2: obtain and need middle set to be processed and element number N thereof:
Step 3: carry out following operation in the face of (A, B) to all of:
Step 3.1: calculate the angle theta between the normal vector of A and BAB;
Step 3.2: if 10 ° of < θAB< 170 °, then perform step 3.3 to 3.10;
Step 3.3: find out and be there is a need in A to the B limit extended and number N thereofA;
Step 3.4: if NA> 0, then perform step 3.5 to 3.6; Otherwise perform step 3.7 to 3.10;
Step 3.5: calculating the limit in A needs to the B length δ extendedA;
Step 3.6: extended δ to the B limit extended by there is a need in AA;
Step 3.7: find out and be there is a need in B to the A limit extended and number N thereofB;
Step 3.8: if NB> 0, then perform step 3.7; Otherwise perform step 3.9;
Step 3.9: calculating the limit in B needs to the A length δ extendedB;
Step 3.10: extended δ to the A limit extended by there is a need in BB;
Step 3.11: utilize B to split A;
Step 3.12: delete the face that in A, area is less;
Step 3.13: utilize A to split B;
Step 3.14: delete the face that in B, area is less;
Step 4: all of face is sewed up;
Step 5: the geometric model of structure is divided finite element grid.
2. one according to claim 1 board-like satellite structure conode grid rapid generation, it is characterised in that: " finding out and being there is a need in A to the B limit extended and number N thereof described in step 3.3A" method be: existing finite element software generally all provides second development interface to enable a user to the new function of exploitation, and realizes in accordance with the following methods by secondary development hence with relevant function interface: 1) the acquisition all of limit E of AA1,…,EAnWith B summit VB;2) to each edge EAj(j=1 ..., n), it is assumed that its natural coordinates is expressed as PjS () (0≤s≤1), obtains N equablyj+ 1 some 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 analytical geometryjk=nBx(xjk-xB)+nBy(yjk-yB)+nBz(zjk-zB) calculate each some P successivelyjkTo the directed distance of B, if | djk| > dmax, then dmax=| djk|, d=djk; Wherein (xjk,yjk,zjk) for PjkCoordinate in rectangular coordinate system, (xB,yB,zB) for VBCoordinate, (nBx,nBy,nBz) for the unit normal vector of B; 4) abd=999999 is taken, to each some PjkIf, d djk> 0 and | djk| < abd, then abd=| djk|; 5) if abd < 0.04, then δmax=abd+0.01; Otherwise δmax=0.05; 6) NA=0, to each edge EAjIf: to k=0,1 ..., NjAll have | djk|≤δmax, then EAjIt is labeled as 1, NA=NA+ 1; Otherwise EAjIt is labeled as 0; The summit of the acquisition all of limit of A and B realizes each through the interface of software in this step, and limit and curve generally can adopt natural coordinates to be indicated in finite element software, after its natural coordinates given, obtain its coordinate representation under rectangular coordinate system by interface; Second development interface function required in this step obtains by inquiring about the dependent help document of finite element software.
3. one according to claim 1 board-like satellite structure conode grid rapid generation, it is characterised in that: " calculating the limit in A and needing the length δ to B extension described in step 3.5A" method be: δA=δmax/sinθAB, wherein δmaxValue obtain in step 3.3, θABValue obtain in step 3.1.
4. one according to claim 1 board-like satellite structure conode grid rapid generation, it is characterised in that: " will be there is a need in A to the B limit extended extension δ described in step 3.6A" method be: existing finite element software generally all provides second development interface to enable a user to the new function of exploitation, will be there is a need to the B limit extended extension δ by secondary development hence with relevant function interface in AA; Each bar limit in one face carries out extend required second development interface function obtain by inquiring about the dependent help document of finite element software.
5. one according to claim 1 board-like satellite structure conode grid rapid generation, it is characterised in that: " finding out and being there is a need in B to the A limit extended and number N thereof described in step 3.7B" method be: existing finite element software generally all provides second development interface to enable a user to the new function of exploitation, and realizes in accordance with the following methods by secondary development hence with relevant function interface: 1) the acquisition all of limit E of BB1,…,EBnWith A summit VA; 2) to each edge EBj(j=1 ..., n) assume that its natural coordinates is expressed as PjS () (0≤s≤1), obtains N equablyj+ 1 some 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 analytical geometryjk=nAx(xjk-xA)+nAy(yjk-yA)+nAz(zjk-zA) calculate each some P successivelyjkTo the directed distance of A, if | djk| > dmax, then dmax=| djk|, d=djk; Wherein (xjk,yjk,zjk) for PjkCoordinate in rectangular coordinate system, (xA,yA,zA) for VACoordinate, (nAx,nAy,nAz) for the unit normal vector of A; 4) abd=999999 is taken, to each some PjkIf, d djk> 0 and | djk| < abd, then abd=| djk|; 5) if abd < 0.04, then δmax=abd+0.01; Otherwise δmax=0.05; 6) NB=0, to each edge EBjIf: to k=1 ..., N all has | djk|≤δmax, then EBjIt is labeled as 1, NB=NB+ 1;Otherwise EBjIt is labeled as 0; The summit of the acquisition all of limit of B and A realizes each through the interface of software in this step, and limit or curve generally can adopt natural coordinates to be indicated in finite element software, after its natural coordinates given, obtain its coordinate representation under rectangular coordinate system by interface; Second development interface function required in this step obtains by inquiring about the dependent help document of finite element software.
6. one according to claim 1 board-like satellite structure conode grid rapid generation, it is characterised in that: " calculating the limit in B and needing the length δ to A extension described in step 3.9B" method be: δB=δmax/sinθAB, wherein δmaxValue obtain in step 3.7, θABValue obtain in step 3.1; Although it should be noted that and this step and step 3.5 all have employed variable δmax, but in step 3.7, its value is altered.
7. one according to claim 1 board-like satellite structure conode grid rapid generation, it is characterised in that: " will be there is a need in B to the A limit extended extension δ described in step 3.10B" method be: existing finite element software generally all provides second development interface to enable a user to the new function of exploitation, will be there is a need to the A limit extended extension δ by secondary development hence with relevant function interface in BB; Each bar limit in one face carries out extend required second development interface function obtain by inquiring about the dependent help document of finite element software.
8. one according to claim 1 board-like satellite structure conode grid rapid generation, it is characterized in that: the method for " utilizing B to split A " described in step 3.11 is: existing finite element software generally all provides second development interface to enable a user to the function that exploitation is new, realize utilizing B to split A by secondary development hence with relevant function interface, only need to specify divisional plane in second development interface function is B, divided is A, and required second development interface function obtains by inquiring about the dependent help document of finite element software.
9. one according to claim 1 board-like satellite structure conode grid rapid generation, it is characterized in that: the method for " deleting the face that in A, area is less " described in step 3.12 is: utilize the second development interface of finite element software to obtain the area of each in A, so just calculated the gross area of all in A by simple addition; If the ratio of the area in certain face and the gross area is less than given threshold gamma, then second development interface function is utilized to be deleted.
10. one according to claim 1 board-like satellite structure conode grid rapid generation, it is characterized in that: the method for " deleting the face that in B, area is less " described in step 3.14 is: utilize the second development interface of finite element software can obtain the area of each in B, so just calculated the gross area of all in B by simple addition; If the ratio of the area in certain face and the gross area is less than given threshold gamma, then second development interface function is utilized to be deleted.
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 true CN105677965A (en) | 2016-06-15 |
| CN105677965B 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) |
Cited By (4)
| 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 |
| CN110263426A (en) * | 2019-06-18 | 2019-09-20 | 珠海格力电器股份有限公司 | A kind of middle face abstracting method, device and the storage medium of thin-wall construction model |
| CN110727989A (en) * | 2018-06-28 | 2020-01-24 | 中车大同电力机车有限公司 | Structural fatigue strength analysis method, device and computer readable storage medium |
| CN115797600A (en) * | 2023-02-08 | 2023-03-14 | 安世亚太科技股份有限公司 | Three-dimensional geometric model envelope 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 |
|---|
| 付志方等: "多工况线性结构稳健拓扑优化设计", 《力学学报》 * |
| 王炳杰等: "基于三维最小二乘方法的空间直线度误差评定", 《北京航空航天大学学报》 * |
Cited By (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 |
| CN110727989A (en) * | 2018-06-28 | 2020-01-24 | 中车大同电力机车有限公司 | Structural fatigue strength analysis method, device and computer readable storage medium |
| CN110727989B (en) * | 2018-06-28 | 2023-08-04 | 中车大同电力机车有限公司 | Structural fatigue strength analysis method, device and computer readable storage medium |
| CN110263426A (en) * | 2019-06-18 | 2019-09-20 | 珠海格力电器股份有限公司 | A kind of middle face abstracting method, device and the storage medium of thin-wall construction model |
| CN115797600A (en) * | 2023-02-08 | 2023-03-14 | 安世亚太科技股份有限公司 | Three-dimensional geometric model envelope grid generation system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105677965B (en) | 2018-09-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105677965A (en) | Method for rapidly generating plate-type satellite-structure node merging grids | |
| CN110188423B (en) | Linear engineering structure rapid BIM modeling method based on finite element meshing | |
| CN106600700B (en) | Three-dimensional model data processing system | |
| CN102693281B (en) | AUTOCAD-based method for generating auxiliary line in PDMS | |
| CN109858161A (en) | A kind of Abaqus Meshing Method based on Midas modeling and Matlab conversion | |
| JPH0816629A (en) | Mesh generating method/device for analysis | |
| CN101908215B (en) | Spatial data fusing method | |
| CN113486429B (en) | Automatic construction method of space intersection structure based on interpolation algorithm | |
| CN113656875A (en) | BIM model lightweight realization method | |
| CN114266093A (en) | Standardized BIM model lightweight realization method | |
| CN101794337A (en) | CAD (Computer-Aided Design)/CAE (Computer-Aided Engineering) system integration method based on layered architecture | |
| CN107895065A (en) | A kind of plant design method based on 3-D technology | |
| CN105701275A (en) | Method for rapidly generating beam plate composite satellite structure common node mesh | |
| CN109933588B (en) | Method and system for converting dwg data into gdb data | |
| CN110096540A (en) | Surveying and mapping data conversion method, equipment, storage medium and device | |
| CN101055601B (en) | Quick inserting pattern frame system and method | |
| CN114022649A (en) | GPU-CPU (graphics processing unit-central processing unit) cooperative raster data rapid coordinate conversion method and system | |
| CN108073776B (en) | Complex river network trunk and branch intersection grid drawing and river center continent grid processing method | |
| US9152741B2 (en) | Three-dimensional shape generation method, program, and recording medium | |
| CN105205263B (en) | A kind of power equipment threedimensional model fast construction equipment and its application | |
| CN108090962B (en) | Method for quickly generating PLAAXIS geometric model based on ANSYS data file | |
| JP7333096B2 (en) | Color information conversion program, color information conversion method, and color information conversion information processing apparatus | |
| Yuan et al. | Accuracy analysis on the automatic registration of multi-source remote sensing images based on the software of ERDAS imagine | |
| JP2012128609A (en) | Drawing creation support method and apparatus | |
| Li et al. | Research on parallel solution of GRAPES Helmholtz equation |
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 |