CN109446541A - A kind of method of body diamond shape cutting finite element grid modeling - Google Patents
A kind of method of body diamond shape cutting finite element grid modeling Download PDFInfo
- Publication number
- CN109446541A CN109446541A CN201811013120.XA CN201811013120A CN109446541A CN 109446541 A CN109446541 A CN 109446541A CN 201811013120 A CN201811013120 A CN 201811013120A CN 109446541 A CN109446541 A CN 109446541A
- Authority
- CN
- China
- Prior art keywords
- node
- grid
- cutting
- diamond shape
- finite element
- 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]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Processing Of Stones Or Stones Resemblance Materials (AREA)
- Laser Beam Processing (AREA)
Abstract
The present invention provides a kind of methods of body diamond shape cutting finite element grid modeling, utilize joint structure rhombohedron grid, polishing body grid front end face and rear end face complete grid dividing, then circumferentially complete cutting with radial Delete Mesh respectively according to the side length, angle and depth of diamond shape cutting unit.The present invention can be realized the body diamond shape cutting finite element grid modeling of different geometrical characteristics, and method is easy and high-efficient.
Description
Technical field
The present invention relates to ammunition technology fields, and in particular to a kind of method of body diamond shape cutting finite element grid modeling.
Background technique
Body cutting be body inner and outer surfaces with machining or other methods process groove, utilize groove
The stress concentration region control body of formation is broken into the fragmentation of predefined size and shape by the grid system of design.Diamond shape cutting is
The mesh shape for referring to cutting is diamond shape.When carrying out body diamond shape cutting design, need to add the body of diamond shape cutting in detonation
Predefined size can be broken under carrying and the fragmentation of shape is calculated.By establish the finite element model of body diamond shape cutting into
Row fragmentation forming Analysis can intuitively obtain the broken situation of the diamond shape cutting body under explosion loading.However, due to body water chestnut
The geometry of shape cutting is complex, using business modeling software realization body diamond shape cutting finite element grid establish it is non-
It is often time-consuming.The threedimensional model for initially setting up diamond shape cutting body is then introduced into finite element software and carries out grid dividing, general to use
Tetrahedral grid divides, and precision is not high；According to hexahedral mesh, need first to divide one of cutting, then mirror
It, cannot direct mirror image as copying to other cuttings, but if the angle of each cutting is different, it is also necessary to repartition, walk
It is rapid cumbersome.In addition, then needing to repartition finite element grid if you need to change cutting characteristic size, a large amount of repetition is grasped in modeling
Work was not only cumbersome but also error-prone, seriously affected design efficiency.
Summary of the invention
In view of this, can be realized not the present invention provides a kind of method of body diamond shape cutting finite element grid modeling
Body diamond shape cutting finite element grid with geometrical characteristic models, and method is easy and high-efficient.
A specific embodiment of the invention is as follows:
A kind of method of body diamond shape cutting finite element grid modeling, comprising the following steps:
Step 1: carrying out body grid dividing using joint structure rhombohedron grid；
Step 2: polishing body grid front end face and rear end face；
Step 3: circumferentially being completed respectively with radial Delete Mesh according to the side length, angle and depth of diamond shape cutting unit
Cutting；
Step 4: exporting the node and unit information of grid according to the call format of finite element software.
Further, the step 1 method particularly includes:
Step 101, the subdivision number for calculating circumferential and axial under cylindrical coordinate, and set radial subdivision number；
Step 102 generates body node according to circumferential, axially and radially the subdivision number；
Step 103, circumferentially moved node, the line and axis after making former same axial node motion are in cutting angle
Half；
Step 104, the net region that radial cutting is determined according to groove depth, body thickness and radial subdivision number；
Step 105 moves radially node and keeps slot bottom node consistent with groove depth；
Step 106, using joint structure diamond-shaped element, obtain rhombohedron grid.
Further, the step 2 method particularly includes:
Step 201, the column of starting point two and two column node of end by duplication axial direction, then along axial positive and negative two
Direction movement obtains front end face node and rear end face node；
Step 202, by joint structure front end unit and backend unit.
Further, the moving method of the step 201 are as follows:
The first row of starting point and last column moving distance of end take L/2Nz, starting point secondary series and end
Secondary series take 3L/2Nz, it is ensured that the coordinate of the node axial direction after movement is all the same, and L is body length, NzFor axial subdivision number.
The utility model has the advantages that
1, the body diamond shape cutting finite element grid that the present invention can be realized different geometrical characteristics models, including interior cutting,
Outer cutting and inside and outside while cutting, since chamfered shape is diamond shape, this method directly constructs rhombohedron grid and is divided,
It is simple easily to realize；To change cutting characteristic size, it is only necessary to the geometrical characteristic parameter of change construction rhombohedron grid
Modeling re-establishes threedimensional model and divides the operation of finite element grid, avoid while avoiding change cutting characteristic size
The a large amount of duplication of labour is brought to analysis personnel, the geometry for solving body diamond shape cutting is complex, finite element grid
The cumbersome problem of establishment process, provide convenience to the analysis of the body of diamond shape cutting, be high-efficient.
2, the present invention constructs rhombohedron grid by way of mobile node, passes through the geometrical characteristic essence of concept transfer
Really control diamond shape groove dimensions, fast response time.
3, the present invention moves axially polishing end face by way of replicating endpoint node, constructs rectangular end face, end face one
Cause property is good, convenient for docking with other models.
4, the present invention limits the moving distance of endpoint node, will not so that the rectangular end face grid constructed is not too big
It is too small, therefore the calculating time of finite element step-length will not be influenced due to sizing grid.
Detailed description of the invention
Fig. 1 is the main view of body diamond shape cutting of the present invention；
Fig. 2 is body diamond shape cutting sectional view of the present invention；
Fig. 3 is diamond shape side length schematic diagram of the present invention；
Fig. 4 is node circumferential direction moving direction schematic diagram of the present invention；
Fig. 5 is the circumferential mobile front nodal point schematic diagram of the present invention；
Fig. 6 is the circumferential mobile posterior nodal point schematic diagram of the present invention；
Fig. 7 is the indexed sequential schematic diagram of 8 joint structures of the invention, one unit；
Fig. 8 is node motion schematic diagram in front and rear end of the present invention；
Fig. 9, Figure 10, Figure 11 are respectively the grid schematic diagram after the outer cutting of the present invention, interior cutting and inside and outside while cutting.
Specific embodiment
The present invention will now be described in detail with reference to the accompanying drawings and examples.
The present invention provides a kind of methods of body diamond shape cutting finite element grid modeling, utilize six face of joint structure diamond shape
Volume mesh, polishing body grid front end face and rear end face complete grid dividing, then according to the side length of diamond shape cutting unit, angle
Cutting circumferentially is completed with radial Delete Mesh respectively with depth, specifically includes the following steps:
Step 1: concluding the geometrical characteristic of body diamond shape cutting；
Cutting width w, groove depth h, cutting interval a, cutting angle, θ, body length L, body are concluded according to Fig. 1, Fig. 2
Inner radial R1, the number of grid N of body thickness T and diamond shape slot on cutting width.
Step 2: calculating the side length of line of rabbet joint intersection gained diamond shape；
As shown in figure 3, line of rabbet joint intersection gained diamond shape side length is indicated with l, i.e. rhombohedron side length of element, diamond shape side length l
It calculates and is calculated by formula (1),
Step 3: calculate cylindrical coordinate under r,The subdivision number in the direction z；
Taking body one end center of circle is cylindrical coordinates origin, and being directed toward the body other end center of circle by origin is z-axis positive direction.It is radial
Subdivision number is Nr, circumferential subdivision number isAxial subdivision number is Nz, radial subdivision number NrA definite value can be given as needed, pressed
Formula (2) calculates circumferential subdivision numberWith axial subdivision number Nz。
Step 4: generating body node；
According to the radial subdivision number N of step 3r, circumferential subdivision numberAxial subdivision number NzBody segment is generated under cylindrical coordinate
Point obtains the node of line by the endpoint of line, by obtaining the node in face to mid-side node, obtains the node of body by opposite node.Due to
Cylinder feature has identical shaft section, determines shaft section Node distribution according to mid-side node, then answers node each on shaft section
With point the body node of body is obtained to the algorithm that axis rotates, all nodes according to first z, thenThen the direction r stores, and facilitates subsequent
Step constructs diamond-shaped element.
Step 5: circumferential mobile node；
Circumferential movement is carried out to node, under cylindrical coordinate, r,Node index respectively i, j, the k in tri- directions z,
In, i=1,2 ... ..., Nr；J=1,2 ... ...,K=1,2 ... ..., Nz.The then moving distance of each nodeBy formula
(3) it calculates, movement counterclockwise is as shown in figure 4, as shown in figure 5, dotted line frame I is mobile front nodal point position, as shown in fig. 6, dotted line
Frame I ' is mobile posterior nodal point position.
Step 6: calculating the grid in radial cutting region；
Radial direction wants the grid number N ' of cuttingrIt is calculated by formula (4).
Step 7: moving radially node keeps slot bottom node consistent with groove depth；
Under cylindrical coordinate, r,The node index in tri- directions z is respectively i, j, k, to index i=Nr-Nr' node
It is moved radially, calculates moving distance by formula (5), keep trench bottom node consistent with groove depth.
Step 8: by joint structure diamond-shaped element；
8 joint structures, one unit, indexed sequential is as shown in fig. 7, using cylindrical coordinate origin as node 1, on the direction r
Point is node 2,Point in plane is node 3,Point on direction is node 4, and the point on the direction z is node 5, rz plane
On point be node 6,Point in space is node 7,Point in plane is node 8.For according to first z, thenThen
The node of all nodes of the direction r storage, unit indexes indexnodeiIt can be calculated by formula (6).The diamond-shaped element of construction is according to elder generation
Radial is the direction r, then is circumferentiallyDirection, then axial is the storage of the direction z, facilitates subsequent cutting.
Step 9: generating front end face and rear end face node；
Then the node of front end face is obtained by the two column node of starting point in the duplication direction z to the movement of z-axis negative direction, rear end
Then the node in face is obtained by the two column node of end in the duplication direction z to the movement of z-axis positive direction, the front end face node of generation
With rear end face node as shown in Fig. 8 dotted line frame, last column moving distance of first row and end for starting point is desirable
L/2Nz, the desirable 3L/2N of the secondary series of starting point secondary series and endz, it is ensured that the coordinate in the direction node z after movement is homogeneous
Together.By the node of duplication and the node being replicated according to first z, thenThen the direction r stores, and facilitates subsequent step structural unit.
Step 10: by joint structure front-end and back-end unit；
By 8 joint structures, one unit, indexed sequential as shown in fig. 7, for according to first z, thenThen the direction r stores
Body node, the node of unit indexes indexnodeiIt can be calculated by formula (7).
Step 11: calculating circumferential cutting interval grid number；
Circumferential cutting interval grid number is calculated by the number of grid of cutting width and single slot, can be calculated by formula (8).
Step 12: deleting cutting area grid；
By circumferential cutting interval grid number and radial cutting grid number determine respectively cutting area grid circumferential index and
Radial to index, grid such as Fig. 9 after cutting can be obtained in all circumferential units for indexing and radially indexing of satisfaction in Delete Mesh,
Shown in Figure 10 and Figure 11, respectively outer cutting, the grid schematic diagram after interior cutting and inside and outside while cutting.
Step 13: cylindrical coordinates transforms to rectangular co-ordinate；
By formula (9) transformation obtain rectangular co-ordinate, can according to the call format of different solvers, Formatting Output grid
Node and unit information.
In conclusion the above is merely preferred embodiments of the present invention, being not intended to limit the scope of the present invention.
All within the spirits and principles of the present invention, any modification, equivalent replacement, improvement and so on should be included in of the invention
Within protection scope.
Claims (4)
1. a kind of method of body diamond shape cutting finite element grid modeling, which comprises the following steps:
Step 1: carrying out body grid dividing using joint structure rhombohedron grid；
Step 2: polishing body grid front end face and rear end face；
Step 3: circumferentially completing cutting with radial Delete Mesh respectively according to the side length, angle and depth of diamond shape cutting unit；
Step 4: exporting the node and unit information of grid according to the call format of finite element software.
2. the method for body diamond shape cutting finite element grid modeling as described in claim 1, which is characterized in that the step 1
Method particularly includes:
Step 101, the subdivision number for calculating circumferential and axial under cylindrical coordinate, and set radial subdivision number；
Step 102 generates body node according to circumferential, axially and radially the subdivision number；
Step 103, circumferentially moved node, line and axis after making former same axial node motion are in the one of cutting angle
Half；
Step 104, the net region that radial cutting is determined according to groove depth, body thickness and radial subdivision number；
Step 105 moves radially node and keeps slot bottom node consistent with groove depth；
Step 106, using joint structure diamond-shaped element, obtain rhombohedron grid.
3. the method for body diamond shape cutting finite element grid modeling as claimed in claim 2, which is characterized in that the step 2
Method particularly includes:
Step 201, the column of starting point two and two column node of end by duplication axial direction, then along axial positive and negative both direction
Movement obtains front end face node and rear end face node；
Step 202, by joint structure front end unit and backend unit.
4. the method for body diamond shape cutting finite element grid modeling as claimed in claim 3, which is characterized in that the step
201 moving method are as follows:
The first row of starting point and last column moving distance of end take L/2Nz, the second of starting point secondary series and end
Column take 3L/2Nz, it is ensured that the coordinate of the node axial direction after movement is all the same, and L is body length, NzFor axial subdivision number.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811013120.XA CN109446541B (en) | 2018-08-31 | 2018-08-31 | Finite element mesh modeling method for projectile rhombus grooving |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811013120.XA CN109446541B (en) | 2018-08-31 | 2018-08-31 | Finite element mesh modeling method for projectile rhombus grooving |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109446541A true CN109446541A (en) | 2019-03-08 |
| CN109446541B CN109446541B (en) | 2022-08-23 |
Family
ID=65533152
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811013120.XA Active CN109446541B (en) | 2018-08-31 | 2018-08-31 | Finite element mesh modeling method for projectile rhombus grooving |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109446541B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111104716A (en) * | 2019-12-09 | 2020-05-05 | 北京航空航天大学 | Automatic generation method of groove type resistance reducing structure based on thermal diffusion facing to blade |
| CN113385901A (en) * | 2021-07-07 | 2021-09-14 | 国基检测有限公司 | Method for controlling notching and fragment breaking of electron beam projectile body |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017001890A1 (en) * | 2015-06-29 | 2017-01-05 | Bosch Car Multimedia Portugal, S.A. | Method for converting gerber data to finite element model for predicting printed circuit board warpage |
| CN108197418A (en) * | 2018-03-14 | 2018-06-22 | 上海理工大学 | A kind of hexahedron FEA Meshing Method for simulating thread fitting |
-
2018
- 2018-08-31 CN CN201811013120.XA patent/CN109446541B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017001890A1 (en) * | 2015-06-29 | 2017-01-05 | Bosch Car Multimedia Portugal, S.A. | Method for converting gerber data to finite element model for predicting printed circuit board warpage |
| CN108197418A (en) * | 2018-03-14 | 2018-06-22 | 上海理工大学 | A kind of hexahedron FEA Meshing Method for simulating thread fitting |
Non-Patent Citations (5)
| Title |
|---|
| 叶帅奇: "锁环式同步器同步环热-结构耦合有限元分析", 《机械工程师》 * |
| 姚兴佳等: "基于MSC.Patran风电机组轮毂的有限元分析", 《沈阳工业大学学报》 * |
| 张健等: "中心孔对刻槽式MEFP成型及毁伤影响的数值模拟", 《沈阳理工大学学报》 * |
| 王忠雷等: "三维有限元六面体网格几何自适应再生成方法", 《计算力学学报》 * |
| 郝东升等: "斜齿轮精确接触分析有限元建模方法", 《大连理工大学学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111104716A (en) * | 2019-12-09 | 2020-05-05 | 北京航空航天大学 | Automatic generation method of groove type resistance reducing structure based on thermal diffusion facing to blade |
| CN113385901A (en) * | 2021-07-07 | 2021-09-14 | 国基检测有限公司 | Method for controlling notching and fragment breaking of electron beam projectile body |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109446541B (en) | 2022-08-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2018059155A1 (en) | Method for constructing three-dimensional entity model having geometric error and computer readable storage medium | |
| CN103247041B (en) | A kind of dividing method of the cloud data of the many geometric properties based on local sampling | |
| JP4381743B2 (en) | Method and program for generating volume data from boundary representation data | |
| CN110704928B (en) | Method for converting BIM model into GIS model | |
| CN106934826B (en) | Rock slope structure refined modeling and block identification method | |
| CN106650085B (en) | Method for generating hollow model by entity model based on finite element mesh | |
| CN103413297A (en) | Cutting method based on integrated three-dimensional GIS model | |
| CN109446541A (en) | A kind of method of body diamond shape cutting finite element grid modeling | |
| Schön et al. | Octree-based indexing for 3D pointclouds within an Oracle Spatial DBMS | |
| CN108763827B (en) | Transmission tower finite element model establishing method and device | |
| CN111105503B (en) | Method for determining explosive-loading combustion surface of solid rocket engine | |
| CN111553008A (en) | Method for converting three-dimensional geological model of side slope into two-dimensional calculation slope | |
| CN111666645A (en) | Modeling method, system and medium for spiral bevel gear based on discrete point data | |
| CN111046571B (en) | Method for quickly determining symmetrical fuel charging combustion surface of solid rocket engine | |
| CN103279985B (en) | A kind of intelligent Modeling method of complicated landform structural system three-dimensional finite element model | |
| CN114510775A (en) | Complex model three-dimensional space curved grid division method | |
| CN107644139B (en) | Attribute mapping method from CAD model to CAE model | |
| CN101303770A (en) | Method for constructing triangle gridding as well as application thereof for geotechnical engineer modeling | |
| CN109902407A (en) | A kind of body finite element grid parametric modeling method of high energy beam processing argyle design | |
| CN109948206B (en) | Flat finite element grid parametric modeling method for processing diamond pattern by high energy beam | |
| CN115619961B (en) | Three-dimensional block cutting and searching method based on matrix expression | |
| Berggren | Shape calculus for fitted and unfitted discretizations: domain transformations vs. boundary-face dilations | |
| CN100464335C (en) | Transfer method from three-dimensional solid to standard technical layout in microelectromechanical system | |
| KR20000072426A (en) | System and method of unstructured tetrahedral mesh generation for a 3d object on semiconductor substrate | |
| Willenborg et al. | Semantic 3D city models serving as information hub for 3D field based simulations |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |