CN109902407A - A kind of body finite element grid parametric modeling method of high energy beam processing argyle design - Google Patents

A kind of body finite element grid parametric modeling method of high energy beam processing argyle design Download PDF

Info

Publication number
CN109902407A
CN109902407A CN201910168172.2A CN201910168172A CN109902407A CN 109902407 A CN109902407 A CN 109902407A CN 201910168172 A CN201910168172 A CN 201910168172A CN 109902407 A CN109902407 A CN 109902407A
Authority
CN
China
Prior art keywords
high energy
grid
energy beam
node
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
CN201910168172.2A
Other languages
Chinese (zh)
Other versions
CN109902407B (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 CN201910168172.2A priority Critical patent/CN109902407B/en
Publication of CN109902407A publication Critical patent/CN109902407A/en
Application granted granted Critical
Publication of CN109902407B publication Critical patent/CN109902407B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Welding Or Cutting Using Electron Beams (AREA)

Abstract

The present invention provides a kind of body finite element grid parametric modeling methods of high energy beam processing argyle design to complete the body grid dividing comprising the high energy beam line of rabbet joint using joint structure rhombohedron grid;Then polishing body grid front end face and rear end face;Then according to the high energy beam line of rabbet joint and the proportionate relationship of body radially and circumferentially, the unit number of high energy beam line of rabbet joint defect area is determined;Secondly modeling circumferentially is completed with the radial area grid for deleting high energy beam line of rabbet joint defect area respectively;The node and unit information of grid are finally exported according to the call format of finite element software.The present invention can be realized the body finite element grid modeling of the high energy beam processing argyle design of different geometrical characteristics, and method is easy and high-efficient.

Description

A kind of body finite element grid parametric modeling method of high energy beam processing argyle design
Technical field
The present invention relates to ammunition technology fields, and in particular to a kind of body finite element grid of high energy beam processing argyle design Parametric modeling method.
Background technique
In order to improve the power of warhead, using high energy beam pre-control technology, body can be made to obtain shape under explosion loading Shape and the more consistent pre-control fragmentation of mass ratio, to improve warhead power.High energy beam pre-control technology is using electron beam, laser Material is locally heated rapidly to melting state by the high energy beams such as beam or ion beam, and body cools down rapidly molten metal, Forming region unit control body is broken on body.High energy beam processing argyle design refers to the unit shape after high energy beam cutting Can shape be diamond shape, need be broken into body under explosion loading predefined size and the diamond shape fragmentation of shape calculates.It is logical It crosses and establishes the body finite element model progress fragmentation forming Analysis based on high energy beam processing argyle design, can intuitively obtain quick-fried The broken situation of the lower high energy beam processing body of Hong load.However, since the geometry of body high energy beam diamond shape cutting is more multiple It is miscellaneous, the region complicated composition of high energy beam processing, comprising melt zone, transition region and defect area, as shown in Figure 1.
However, utilizing business modeling software realization body diamond shape since the geometry of body diamond shape cutting is complex The foundation of the finite element grid of cutting is very time-consuming.The threedimensional model for initially setting up diamond shape cutting body, is then introduced into finite element Software carries out grid dividing, is generally divided using tetrahedral grid, and precision is not high;According to hexahedral mesh, need first to it In a cutting divided, then image copying is to other cuttings, but if the angle of each cutting is different, cannot be straight Connect mirror image, it is also necessary to repartition, complex steps.In addition, then needing to repartition finite element if you need to change cutting characteristic size Grid, a large amount of repetitive operation was not only cumbersome but also error-prone in modeling, seriously affected design efficiency.
Summary of the invention
In view of this, the present invention provides a kind of body finite element grid parametric modelings of high energy beam processing argyle design Method, can be realized the body finite element grid modeling of the high energy beam processing argyle design of different geometrical characteristics, method it is easy and It is high-efficient.
The technical solution adopted by the present invention is as follows:
A kind of body finite element grid parametric modeling method of high energy beam processing argyle design, comprising the following steps:
Step 1: completing the body grid dividing comprising the high energy beam line of rabbet joint using joint structure rhombohedron grid;
Step 2: polishing body grid front end face and rear end face;
Step 3: determining high energy beam line of rabbet joint defect area according to the high energy beam line of rabbet joint and the proportionate relationship of body radially and circumferentially Unit number;
Step 4: circumferentially completing modeling with the radial area grid for deleting high energy beam line of rabbet joint defect area respectively;
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 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 grid that high energy beam line of rabbet joint radial direction cutting is determined according to groove depth, body thickness and radial subdivision number Region;
Step 105 moves radially node and keeps slot bottom node consistent with groove depth;Moving radially node simultaneously lacks the line of rabbet joint It is consistent with defect area height to fall into area's upper layer node;
Step 106, using joint structure diamond-shaped element, obtain the rhombohedron grid comprising the high energy beam line of rabbet joint.
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.
Further, the step 3 method particularly includes:
Step 301 determines transition region territory element by circumferential cutting spacer units number and radial cutting unit number;
Step 302 determines melt zone territory element by circumferential cutting spacer units number and radial melt zone unit number;
Step 303 determines defect area territory element by circumferential cutting spacer units number and radial defect area unit number.
The utility model has the advantages that
1, the present invention can be realized the body finite element grid modeling of the high energy beam processing argyle design of different geometrical characteristics, Since processing graphic pattern is diamond shape, this method directly constructs rhombohedron grid and is divided, simple easily to realize;It is high to change Beam score line feature size, it is only necessary to which the geometrical characteristic parameter of change construction rhombohedron grid can model, and avoid change It while high energy beam score line feature size, re-establishes threedimensional model and divides the operation of finite element grid, avoid to analysis people Member brings a large amount of duplication of labour, solve high energy beam processing argyle design body geometry is complex, finite element net The cumbersome problem of the establishment process of lattice provides convenience to the body analysis of high energy beam processing argyle design, is high-efficient.
2, the present invention constructs rhombohedron grid by way of mobile node, passes through the geometrical characteristic essence of concept transfer The really size of the control high energy beam line of rabbet joint, 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 model area distribution schematic diagram processed based on high energy beam;
Fig. 2 is the main view of body rhombohedron grid of the present invention;
Fig. 3 is the partial cross-sectional view of body of the present invention;
Fig. 4 is diamond shape side length schematic diagram of the present invention;
Fig. 5 is node circumferential direction moving direction schematic diagram of the present invention;
Fig. 6 is the circumferential mobile front nodal point schematic diagram of the present invention;
Fig. 7 is the circumferential mobile posterior nodal point schematic diagram of the present invention;
Fig. 8 is the indexed sequential schematic diagram of eight joint structures of the invention, one unit;
Fig. 9 is node motion schematic diagram in front and rear end of the present invention;
Figure 10 is body rhombohedron grid dividing partial schematic diagram of the present invention;
Figure 11 is the rhombohedron grid dividing partial schematic diagram of the high energy beam line of rabbet joint of the present invention;
Figure 12 is left view local after body grid dividing of the present invention.
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 body finite element grid parametric modeling methods of high energy beam processing argyle design, utilize Joint structure rhombohedron grid carries out the body grid dividing comprising the high energy beam line of rabbet joint, and polishing body grid front end face is with after End face determines the unit of high energy beam line of rabbet joint defect area then according to the high energy beam line of rabbet joint and the proportionate relationship of body radially and circumferentially Number circumferentially completes modeling with the radial area grid for deleting high energy beam line of rabbet joint defect area respectively, such as Figure 10, Figure 11 and Figure 12 institute Show, specifically includes the following steps:
Step 1, the geometrical characteristic for concluding the diamond shape cutting of body high energy beam;
Cutting width, that is, line of rabbet joint width w, groove depth, that is, line of rabbet joint depth h, cutting interval a, cutting angle, θ, body length L, body inner radial R1, body thickness T and the line of rabbet joint intersect gained diamond shape along the element number N of cutting width.Defect area height h0, melt the ratio S that sector width accounts for cutting width0, as shown in Figure 2, Figure 3 and Figure 4.Take line of rabbet joint width w=1mm, line of rabbet joint depth h =5mm, cutting interval a=10mm, cutting angle, θ=93 °, body length L=50mm, body inner radial R1=20mm, bullet Body thickness T=10mm and the line of rabbet joint intersect gained diamond shape along the element number N=3 of cutting width.Defect area height h0=1mm, melts Change the ratio S that sector width accounts for cutting width0=1/3.
Step 2, the side length for calculating 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 calculates to obtain l=1.001mm 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, obtainNz=219.
Step 4 generates body node;
According to the radial subdivision number N of step 3r, circumferential subdivision numberAxial subdivision number NzBody is generated under cylindrical coordinate Internal node obtains the node of line by the endpoint of line, by obtaining the node in face to mid-side node, obtains the section of body by opposite node Point.Since cylinder feature has identical shaft section, shaft section Node distribution is determined according to mid-side node, then to each on shaft section Node using point the node of body is obtained to the algorithm that axis rotates, all nodes according to first z, thenThen the direction r stores, convenient 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, Wherein, i=1,2 ... ..., NrK=1,2 ... ..., Nz.The then moving distance of each nodeIt is calculated by formula (3), movement counterclockwise is as shown in figure 5, as shown in fig. 6, dotted line frame I is mobile front nodal point position, such as Fig. 7 institute Show, dotted line frame I ' is mobile posterior nodal point position.
Step 6, the grid for calculating radial cutting region;
Radial direction wants the grid number N ' of cuttingrN is calculated by formula (4)r'=5.
Step 7 moves radially node and 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, calculate to obtain Δ r=0.
Step 8 calculates defect area radial direction unit number;
Defect area radial direction unit number HrIt indicates, calculates to obtain H by formula (6)r=1.
Step 9 moves radially node and keeps defect area upper layer node consistent with defect area height;
Under cylindrical coordinate, r,The node index in tri- directions z is respectively i, j, k, to index i=Nr-N'r+Hr's Node is moved radially, and moving distance is calculated by formula (7), and defect area upper layer node is consistent with defect area height, calculates here Obtain Δ rh=0.
Step 10, by joint structure diamond-shaped element;
Eight joint structures, one unit, indexed sequential is as shown in figure 8, using cylindrical coordinate origin as node 1, on the direction r Point be node 2,Point in plane is node 3,Point on direction is node 4, and the point on the direction z is node 5, and rz is flat Point on face is node 6,Point in space is node 7,Point in plane is node 8.For according to first z, then Then the node of all nodes of the direction r storage, unit indexes indexnodeiIt can be calculated by formula (8).The diamond-shaped element of construction is pressed According to the first radial i.e. direction r, then circumferentially it isDirection, then axial is the storage of the direction z, facilitates subsequent cutting.
Step 11 generates 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. 9 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 12, 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 Node, the node of unit indexes indexnodeiIt can be calculated by formula (9).
Step 13 calculates 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 (10) ?
Step 14, setting transition region;
By circumferential cutting spacer units number and radial cutting unit number determine respectively cutting territory element circumferential index and Radial to index, all circumferential units for indexing and radially indexing of satisfaction are transition region in setting unit.
Step 15 calculates melt zone radial direction unit number and melts unit number on sector width;
Melt zone radial direction unit number H is calculated by formula (11)m, Hm=4.
Hm=Nr-Hr (11)
It is calculated by formula (12) and melts unit number N on sector widthm, Nm=1.
Nm=NS0 (12)
Step 16, setting melt zone;
Determine the circumferential index for melting territory element respectively by circumferential cutting spacer units number and radial melt zone unit number With radial index, all circumferential units for indexing and radially indexing of satisfaction are melt zone in setting unit.
Step 17 deletes defect area unit;
Determine the circumferential index of defect area unit respectively by circumferential cutting spacer units number and radial defect area unit number With radial index, all circumferential units for indexing and radially indexing of satisfaction are defect area in setting unit.
Step 18, cylindrical coordinates transform to rectangular co-ordinate;
Rectangular co-ordinate is obtained by formula (13) transformation, it can be 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 (5)

1. a kind of body finite element grid parametric modeling method of high energy beam processing argyle design, which is characterized in that including with Lower step:
Step 1: completing the body grid dividing comprising the high energy beam line of rabbet joint using joint structure rhombohedron grid;
Step 2: polishing body grid front end face and rear end face;
Step 3: determining the list of high energy beam line of rabbet joint defect area according to the high energy beam line of rabbet joint and the proportionate relationship of body radially and circumferentially First number;
Step 4: circumferentially completing modeling with the radial area grid for deleting high energy beam line of rabbet joint defect area respectively;
Step 4: exporting the node and unit information of grid according to the call format of finite element software.
2. the body finite element grid parametric modeling method of high energy beam processing argyle design as described in claim 1, special Sign is, 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 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 grid regions that high energy beam line of rabbet joint radial direction cutting is determined according to groove depth, body thickness and radial subdivision number Domain;
Step 105 moves radially node and keeps slot bottom node consistent with groove depth;Moving radially node makes line of rabbet joint defect area simultaneously Upper layer node is consistent with defect area height;
Step 106, using joint structure diamond-shaped element, obtain the rhombohedron grid comprising the high energy beam line of rabbet joint.
3. the body finite element grid parametric modeling method of high energy beam processing argyle design as claimed in claim 2, special Sign is, 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 body finite element grid parametric modeling method of high energy beam processing argyle design as claimed in claim 3, special Sign is, 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, 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.
5. the body finite element grid parametric modeling method of high energy beam processing argyle design as described in claim 1, special Sign is, the step 3 method particularly includes:
Step 301 determines transition region territory element by circumferential cutting spacer units number and radial cutting unit number;
Step 302 determines melt zone territory element by circumferential cutting spacer units number and radial melt zone unit number;
Step 303 determines defect area territory element by circumferential cutting spacer units number and radial defect area unit number.
CN201910168172.2A 2019-03-06 2019-03-06 Projectile finite element grid parametric modeling method for processing diamond patterns by high-energy beams Expired - Fee Related CN109902407B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910168172.2A CN109902407B (en) 2019-03-06 2019-03-06 Projectile finite element grid parametric modeling method for processing diamond patterns by high-energy beams

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910168172.2A CN109902407B (en) 2019-03-06 2019-03-06 Projectile finite element grid parametric modeling method for processing diamond patterns by high-energy beams

Publications (2)

Publication Number Publication Date
CN109902407A true CN109902407A (en) 2019-06-18
CN109902407B CN109902407B (en) 2020-10-27

Family

ID=66946506

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910168172.2A Expired - Fee Related CN109902407B (en) 2019-03-06 2019-03-06 Projectile finite element grid parametric modeling method for processing diamond patterns by high-energy beams

Country Status (1)

Country Link
CN (1) CN109902407B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110218781A1 (en) * 2010-03-08 2011-09-08 Livermore Software Technology Corporation Fem-bem coupling methods and systems for sliding contact interface
CN104484489A (en) * 2014-07-24 2015-04-01 江苏科技大学 Automatic generation method for quadrilateral finite element mesh of pitting corrosion damage cylindrical shell
CN105631117A (en) * 2015-12-28 2016-06-01 株洲时代新材料科技股份有限公司 Spline coupling anti-side-rolling torsion bar grid division method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110218781A1 (en) * 2010-03-08 2011-09-08 Livermore Software Technology Corporation Fem-bem coupling methods and systems for sliding contact interface
CN104484489A (en) * 2014-07-24 2015-04-01 江苏科技大学 Automatic generation method for quadrilateral finite element mesh of pitting corrosion damage cylindrical shell
CN105631117A (en) * 2015-12-28 2016-06-01 株洲时代新材料科技股份有限公司 Spline coupling anti-side-rolling torsion bar grid division method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈帅: "弹体改性区对弹丸过载性能和破片形成的影响", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 *

Also Published As

Publication number Publication date
CN109902407B (en) 2020-10-27

Similar Documents

Publication Publication Date Title
CN106513679B (en) A kind of laser gain material manufacturing method suitable for band large-scale metal part
CN108038266A (en) A kind of method for numerical simulation of selective laser repair process
CN106934826B (en) Rock slope structure refined modeling and block identification method
CN104765936A (en) Operation method of three-dimensional modeling expansion and lofting coordinate generation based on Autodesk inventor software
CN110188423B (en) Linear engineering structure rapid BIM modeling method based on finite element meshing
CN106650085B (en) Method for generating hollow model by entity model based on finite element mesh
KR20190044090A (en) Structural heat exchanger
CN112182921A (en) Prediction method for selective laser melting thermal-mechanical coupling behavior of high-performance alloy steel
CN104732028A (en) Method for generating powder laying and sintering supporting structure of aircraft structural part based on space lattice
WO2022036591A1 (en) Method and device for formulating printing process in additive manufacturing
CN111680338B (en) Modeling method of gradient porous structure based on STL file
CN115797568A (en) Modeling method and device based on three-dimensional GIS and BIM integration
Rosso et al. Optimization approaches in design for additive manufacturing
KR101462861B1 (en) System and Method for outputting hot working of curved board
CN109902407A (en) A kind of body finite element grid parametric modeling method of high energy beam processing argyle design
CN107134007A (en) The three-dimensional fast modeling method of traffic infrastructure
CN109948206A (en) A kind of plate finite element grid parametric modeling method of high energy beam processing argyle design
CN111080777A (en) Three-dimensional rapid modeling method for spacecraft thermal control product
CN116484668A (en) Electron beam additive manufacturing process simulation method
CN107977529B (en) Weldment modeling method based on three-dimensional software
CN110781555A (en) Structural member information symmetry method
US10962958B2 (en) Library of predefined shapes for additive manufacturing processes
Jin et al. Optimal design and modeling of 3D variable-density lattice structures
CN114472924A (en) Laser path planning method, system, computer and readable storage medium
CN117688813A (en) Fracture parameter calculation sub-model method for crack-containing structure under residual stress

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
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: 20201027