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 PDF

Info

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
Application number
CN201811013120.XA
Other languages
Chinese (zh)
Other versions
CN109446541B (en
Inventor
蒋建伟
邱浩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201811013120.XA priority Critical patent/CN109446541B/en
Publication of CN109446541A publication Critical patent/CN109446541A/en
Application granted granted Critical
Publication of CN109446541B publication Critical patent/CN109446541B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling 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)
  • Laser Beam Processing (AREA)
  • Processing Of Stones Or Stones Resemblance Materials (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

A kind of method of body diamond shape cutting finite element grid modeling
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.
CN201811013120.XA 2018-08-31 2018-08-31 Finite element mesh modeling method for projectile rhombus grooving Active CN109446541B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
叶帅奇: "锁环式同步器同步环热-结构耦合有限元分析", 《机械工程师》 *
姚兴佳等: "基于MSC.Patran风电机组轮毂的有限元分析", 《沈阳工业大学学报》 *
张健等: "中心孔对刻槽式MEFP成型及毁伤影响的数值模拟", 《沈阳理工大学学报》 *
王忠雷等: "三维有限元六面体网格几何自适应再生成方法", 《计算力学学报》 *
郝东升等: "斜齿轮精确接触分析有限元建模方法", 《大连理工大学学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
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
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
CN109977551A (en) A kind of lattice material modeling method based on NURBS
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