CN107679341A - A kind of barrel configuration parametric Finite Element Modeling Method - Google Patents
A kind of barrel configuration parametric Finite Element Modeling Method Download PDFInfo
- Publication number
- CN107679341A CN107679341A CN201711034851.8A CN201711034851A CN107679341A CN 107679341 A CN107679341 A CN 107679341A CN 201711034851 A CN201711034851 A CN 201711034851A CN 107679341 A CN107679341 A CN 107679341A
- Authority
- CN
- China
- Prior art keywords
- barrel
- node
- finite element
- python
- fem
- 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
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000013515 script Methods 0.000 claims abstract description 8
- 238000010586 diagram Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 210000000038 chest Anatomy 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000000843 powder Substances 0.000 description 2
- 241000282470 Canis latrans Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000008188 pellet Substances 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
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]
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)
- Image Generation (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a kind of barrel configuration parametric Finite Element Modeling Method, comprise the following steps：Straight cylinder barrel 3-D geometric model is established by finite element software；Finite element grid is divided to geometrical model；The execution journal File that finite element CAE software automatically generates, file name suffix .rpy is changed to Python shell script file .py, forms Python scripts, and barrel configuration parameter is replaced using variable, to realize parametric modeling；Export the Inp files of straight cylinder barrel FEM model；The node coordinate along each layer unit of barrel axis included in Inp files is read using Python, forms barrel FEM mesh；Inp is updated, and is imported into commercial finite element CAE software, that is, generates barrel FEM mesh.
Description
Technical field
The present invention relates to a kind of computer simulation technique, particularly a kind of barrel configuration parametric Finite Element Modeling Method.
Background technology
Barrel is to provide the basic components of gun interior ballistics process, bears high powder pressure.Pellet injection mistake
Cheng Zhong, powder burning moment produces a large amount of high temperature, gases at high pressure promote bullet to be moved along barrel axis direction, in the process,
With Pipe rifling mutually extruding occurs for bullet, and bearing band material undergoes the process such as high transient state large deformation, friction.Bullet is in chamber inner movement
It is a nonlinearity problem.First motion process in Projectile is studied, has important meaning for research Projectile Motion rule
Justice.With the development of science and technology, finite element method, which is increasingly becoming, studies such problem main stream approach.In FEM model,
The construction and grid precision of barrel finite element grid are particularly important to the accuracy of simulation result.
The interior thorax of rifling barrel is complicated, and stand under load is severe., be to Pipe rifling geometry during using finite element method
Entity discreteness.General barrel has more riflings, and the cross section geometry of rifling is smaller relative to barrel bore, and around barrel
Axis spiral rotating, influence the foundation of barrel finite element grid and the division of good cell.
In finite element analysis technology field, more than the 60% of pre-treatment modeling shared whole analysis time, routinely build
Mould technological means is the barrel geometrical model established according to design drawing in 3 d modeling software after simplifying, and is then introduced into preceding
Handle in software, carry out mesh generation, finally the grid model after dividing well is imported into finite element software and carries out emulation meter
Calculate.The mutual importing of each data between software, do not only result in it is cumbersome, when and easily causing loss of data, and spending
Between it is longer, grid model precision is poor.In addition, when preceding processing carries out barrel mesh generation, due to the complexity of barrel configuration,
The mesh quality of division is not fine.
The content of the invention
The present invention provides a kind of barrel configuration parametric Finite Element Modeling Method, can be realized to barrel by secondary development
Coyote hole, forcing cone structure and guide part and barrel exterior contour carry out quick, accurate parametersization modeling.
The technical scheme for realizing the object of the invention is：A kind of barrel configuration parametric Finite Element Modeling Method, including it is following
Step：
Step 1, straight cylinder barrel 3-D geometric model is established by finite element software；
Step 2, finite element grid is divided to geometrical model；
Step 3, the execution journal File that finite element CAE software automatically generates, by file name suffix .rpy
Python shell script file .py are changed to, form Python scripts, and barrel configuration parameter is replaced using variable, with
Realize parametric modeling；
Step 4, the Inp files of straight cylinder barrel FEM model are exported；
Step 5, the node along each layer unit of barrel axis included in Inp files is read using Python to sit
Mark, form barrel FEM mesh；
Step 6, Inp is updated, and is imported into commercial finite element CAE software, that is, generates barrel FEM mesh.
Using the above method, the detailed process of step 5 is：
Step 5.1, the z directions coordinate using the front/rear end node of each layer unit of barrel is identical, navigates to rifling and rises
Beginning portion's unit, the end face node in breech direction is will be close to as source face node, other end node is as whole face node；
Step 5.2, using source face node as fixed point, whole face node is as skew node, source face node to whole face node
For barrel axial direction；
Step 5.3, Arbitrary Rotation is carried out according to barrel degree of twining and rifling form to a node (x, y) to offset to obtain newly
Coordinate (x', y')
Step 5.4, next layer unit node offset in the unit basis of last layer repeat more than operate both can be with complete
The skew of barrel node is helped, forms barrel FEM mesh.
This method quickly and efficiently can accurately establish the FEM mesh of barrel configuration, and can be to barrel knot
Structure carries out parametric modeling, to realize the foundation of the FEM mesh of the barrel to any different bores.
The present invention is described further with reference to Figure of description.
Brief description of the drawings
Fig. 1 is the method flow schematic diagram of the present invention.
Fig. 2 is Pipe rifling structural representation.
Fig. 3 is straight cylinder barrel grid FEM model schematic diagram.
Fig. 4 is elementary layer identification schematic diagram.
Fig. 5 is cell node skew principle schematic.
Fig. 6 is node calculations of offset schematic diagram.
Fig. 7 is certain 76mm gun barrel FEM mesh schematic diagram.
Embodiment
With reference to Fig. 1 to Fig. 6, a kind of barrel configuration parametric modeling method, comprise the following steps：
Step 1, according to thorax structure actual size in barrel, in ABAQUS softwares, according to center roundlet as shown in Figure 2
(center axis of small circle is negative line width) method, establishes Artillery rifle section.And barrel section is stretched, form straight cylinder
Barrel.
Step 2, finite element grid is divided to straight cylinder gun barrel, as shown in Figure 3.
Step 3, during straight cylinder barrel grid is generated, ABAQUS softwares automatically generate an execution journal file
(.rpy) file, almost each step CAE operations can find corresponding sentence in this document, file name suffix are changed to
Python shell scripts file (.py), Python scripts are formed, and barrel configuration parameter is replaced using variable, with side
Just parametric modeling is carried out.
Step 4, Inp file of the generation comprising unit and node；
Step 5, the node along each layer unit of barrel axis included in Inp files is read using Python to sit
Mark, form barrel FEM mesh.
Such as：As shown in figure 4, it is identical using the z directions coordinate of the front/rear end node of each layer unit of barrel, navigate to
Rifling initial part unit (assuming that the 3rd layer unit is rifling initial part unit), will be close to the end face node in breech direction as source
Face node, other end node is as whole face node.After the source face node and the whole face node that identify unit, source face node is made
For fixed point, whole face node is as skew node.As shown in figure 5, for example, in the unit, node 100,101,102,103 is
Source face node, node 200,201,202,203 are whole face node, and source face node to whole face node is barrel axial direction.Node
During skew, by taking node 203 as an example, node 203 is one of the unit end face node, and another unit source face node it
One, this node is entered into line displacement according to barrel degree of twining and rifling form, rotated according to upper 1: 203 (x, y) of same circle any
After angle, the new coordinate 203 (x', y') after being offset by formula (1), that is, the node skew of the layer unit, such as Fig. 6 are completed
It is shown.The node of next layer unit is offset repeated in the unit basis of last layer more than operate and can both complete whole body pipe node
Skew, form barrel FEM mesh.Building for barrel FEM mesh is carried out in finite element unit nodal basis
Mould, the model accuracy of foundation is higher, and models and spend the period, relative to traditional modeling method, will not produce the geometry of model
Distortion and loss of data, it can provide precision higher barrel FEM mesh to be related to the correlation of barrel FEM calculation.
Step 6, amended node coordinate is imported into new Inp files, that is, forms barrel FEM mesh.
Fig. 7 show the 76mm gun barrel FEM meshes generated using this method.
Claims (2)
1. a kind of barrel configuration parametric Finite Element Modeling Method, it is characterised in that comprise the following steps：
Step 1, straight cylinder barrel 3-D geometric model is established by finite element software；
Step 2, finite element grid is divided to geometrical model；
Step 3, the execution journal File that finite element CAE software automatically generates, file name suffix .rpy is changed to
Python shell script file .py, Python scripts are formed, and barrel configuration parameter is replaced using variable, to realize
Parametric modeling；
Step 4, the Inp files of straight cylinder barrel FEM model are exported；
Step 5, the node coordinate along each layer unit of barrel axis included in Inp files, shape are read using Python
Into barrel FEM mesh；
Step 6, Inp is updated, and is imported into commercial finite element CAE software, that is, generates barrel FEM mesh.
2. according to the method for claim 1, it is characterised in that the detailed process of step 5 is：
Step 5.1, the z directions coordinate using the front/rear end node of each layer unit of barrel is identical, navigates to rifling initial part
Unit, the end face node in breech direction is will be close to as source face node, other end node is as whole face node；
Step 5.2, using source face node as fixed point, for whole face node as skew node, source face node to whole face node is body
Pipe axial direction；
Step 5.3, Arbitrary Rotation is carried out according to barrel degree of twining and rifling form to a node (x, y) to offset to obtain new coordinate
(x',y')
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mi>x</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
<mo>-</mo>
<mi>y</mi>
<mi> </mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mi>y</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mi>x</mi>
<mi> </mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
<mo>+</mo>
<mi>y</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
Step 5.4, next layer unit node offset in the unit basis of last layer repeat more than operation can both complete entirely
Barrel node is offset, and forms barrel FEM mesh.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711034851.8A CN107679341B (en) | 2017-10-30 | 2017-10-30 | Finite element parametric modeling method for barrel structure |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711034851.8A CN107679341B (en) | 2017-10-30 | 2017-10-30 | Finite element parametric modeling method for barrel structure |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107679341A true CN107679341A (en) | 2018-02-09 |
| CN107679341B CN107679341B (en) | 2021-03-02 |
Family
ID=61142397
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711034851.8A Expired - Fee Related CN107679341B (en) | 2017-10-30 | 2017-10-30 | Finite element parametric modeling method for barrel structure |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107679341B (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108959786A (en) * | 2018-07-12 | 2018-12-07 | 华东交通大学 | A kind of parametric modeling method of imitative mother-of-pearl laminated composite plate |
| CN109858161A (en) * | 2019-02-01 | 2019-06-07 | 东北大学 | A kind of Abaqus Meshing Method based on Midas modeling and Matlab conversion |
| CN111125956A (en) * | 2019-12-23 | 2020-05-08 | 奇瑞汽车股份有限公司 | Automobile finite element model assembling method |
| CN111736530A (en) * | 2020-06-19 | 2020-10-02 | 山东大学 | Method and system for simulating tool wear morphology in machining process |
| CN112949128A (en) * | 2021-03-03 | 2021-06-11 | 中铁十八局集团有限公司 | Finite element analysis method for outer tank model of LNG storage tank |
| CN113190905A (en) * | 2021-05-13 | 2021-07-30 | 中国电子工程设计院有限公司 | Building model analysis method and device and storage medium |
| CN113792386A (en) * | 2021-09-27 | 2021-12-14 | 中国第一汽车股份有限公司 | Finite element automatic modeling method for suspension system parts |
| CN114970196A (en) * | 2022-06-14 | 2022-08-30 | 北京科技大学 | Method and device for synchronously and rapidly generating model and structured grid of metallurgical container |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5279446A (en) * | 1991-01-11 | 1994-01-18 | The Cornelius Company | Beverage cooling system |
| CN104281726A (en) * | 2014-05-27 | 2015-01-14 | 北京宇航系统工程研究所 | Unfolding finite element modeling method for complex cylinder structural plane |
| CN106372288A (en) * | 2016-08-26 | 2017-02-01 | 伊顿上飞（上海）航空管路制造有限公司 | Pre-processing module for development and automatic interaction analysis of aerial hydraulic pipelines |
| CN107145651A (en) * | 2017-04-24 | 2017-09-08 | 西南交通大学 | The three-dimensional infinite element border fast modeling methods of ABAQUS based on INP files |
| CN107229786A (en) * | 2017-05-24 | 2017-10-03 | 上海理工大学 | Model reconstruction method based on FEM post-processing |
| CN107256304A (en) * | 2017-06-02 | 2017-10-17 | 中北大学 | A kind of gun bullet emissive porwer reliability estimation method |
-
2017
- 2017-10-30 CN CN201711034851.8A patent/CN107679341B/en not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5279446A (en) * | 1991-01-11 | 1994-01-18 | The Cornelius Company | Beverage cooling system |
| CN104281726A (en) * | 2014-05-27 | 2015-01-14 | 北京宇航系统工程研究所 | Unfolding finite element modeling method for complex cylinder structural plane |
| CN106372288A (en) * | 2016-08-26 | 2017-02-01 | 伊顿上飞（上海）航空管路制造有限公司 | Pre-processing module for development and automatic interaction analysis of aerial hydraulic pipelines |
| CN107145651A (en) * | 2017-04-24 | 2017-09-08 | 西南交通大学 | The three-dimensional infinite element border fast modeling methods of ABAQUS based on INP files |
| CN107229786A (en) * | 2017-05-24 | 2017-10-03 | 上海理工大学 | Model reconstruction method based on FEM post-processing |
| CN107256304A (en) * | 2017-06-02 | 2017-10-17 | 中北大学 | A kind of gun bullet emissive porwer reliability estimation method |
Non-Patent Citations (1)
| Title |
|---|
| 丁传俊 等: "身管内膛参数化模型及其磨损有限元模型的生成方法", 《兵工学报》 * |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108959786A (en) * | 2018-07-12 | 2018-12-07 | 华东交通大学 | A kind of parametric modeling method of imitative mother-of-pearl laminated composite plate |
| CN108959786B (en) * | 2018-07-12 | 2022-09-16 | 华东交通大学 | Parametric modeling method of nacre-like layered composite board |
| CN109858161A (en) * | 2019-02-01 | 2019-06-07 | 东北大学 | A kind of Abaqus Meshing Method based on Midas modeling and Matlab conversion |
| CN111125956A (en) * | 2019-12-23 | 2020-05-08 | 奇瑞汽车股份有限公司 | Automobile finite element model assembling method |
| CN111125956B (en) * | 2019-12-23 | 2023-08-22 | 奇瑞汽车股份有限公司 | Automobile finite element model assembling method |
| CN111736530A (en) * | 2020-06-19 | 2020-10-02 | 山东大学 | Method and system for simulating tool wear morphology in machining process |
| CN111736530B (en) * | 2020-06-19 | 2021-10-26 | 山东大学 | Method and system for simulating tool wear morphology in machining process |
| CN112949128A (en) * | 2021-03-03 | 2021-06-11 | 中铁十八局集团有限公司 | Finite element analysis method for outer tank model of LNG storage tank |
| CN113190905A (en) * | 2021-05-13 | 2021-07-30 | 中国电子工程设计院有限公司 | Building model analysis method and device and storage medium |
| CN113190905B (en) * | 2021-05-13 | 2024-05-10 | 中国电子工程设计院有限公司 | Building model analysis method, device and storage medium |
| CN113792386A (en) * | 2021-09-27 | 2021-12-14 | 中国第一汽车股份有限公司 | Finite element automatic modeling method for suspension system parts |
| CN114970196A (en) * | 2022-06-14 | 2022-08-30 | 北京科技大学 | Method and device for synchronously and rapidly generating model and structured grid of metallurgical container |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107679341B (en) | 2021-03-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107679341A (en) | A kind of barrel configuration parametric Finite Element Modeling Method | |
| Wang et al. | Numerical simulation of explosive welding using the material point method | |
| CN103366052A (en) | Method for analyzing hot pneumostatic elasticity of hypersonic aerocraft | |
| Yu et al. | Aerodynamic analysis of projectile in gun system firing process | |
| Pengfei et al. | Solid rocket motor propellant grain burnback simulation based on fast minimum distance function calculation and improved marching tetrahedron method | |
| Shen et al. | Implementation of density-based implicit LU-SGS solver in the framework of OpenFOAM | |
| Zhao et al. | Optimization of the aerodynamic configuration of a tubular projectile based on blind kriging | |
| Yang et al. | An efficient discrete velocity method with inner iteration for steady flows in all flow regimes | |
| Li et al. | Ground effects on the hypervelocity jet flow and the stability of projectile | |
| Zhang | Research progress of hypersonic inlet inverse design based on curved shock compression system | |
| CN104834795B (en) | Band connection structure nonlinear contact with friction simulated behavior method and system | |
| Li et al. | Analysis of ground distance effect on hypervelocity jet flow and stability of projectile in half-space | |
| CN110276142A (en) | The emulation mode of plasma density Conduce Disciplinarian mathematical model under high pressure conditions | |
| Yu et al. | Numerical simulation and analysis of the muzzle flow during the revolving barrel gun firing | |
| Karimian et al. | Immersed boundary method for the solution of 2d inviscid compressible flow using finite volume approach on moving cartesian grid | |
| Liu et al. | RSAL-iMFS: A framework of randomized stacking with active learning for incremental multi-fidelity surrogate modeling | |
| Cook et al. | Applications of the CESE method in LS-DYNA | |
| Wang et al. | Application of Three-Dimensional Meshless Method in Muzzle Flow Field of Projectile with Large Displacement. | |
| Hu et al. | A thermal compensated model for predicting the internal combustion of energetic materials in a combustion chamber | |
| Yunpeng et al. | Theories and methods for designing hypersonic high-enthalpy flow nozzles | |
| Phoenix et al. | Mach 5–3.5 morphing waverider accuracy and aerodynamic performance evaluation | |
| CN109408908A (en) | Bearing band based on two phase flow inner trajectory gets into numerical simulation method | |
| CN112084723A (en) | Method and device for simulating one-step stretching process of optical fiber preform | |
| Alfaro et al. | Recent advances in the meshless simulation of aluminium extrusion and other related forming processes | |
| CN115048730B (en) | Axial symmetry supersonic velocity spray pipe optimization design method and device based on grid displacement |
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 |
Granted publication date: 20210302 |