CN107679341B - Finite element parametric modeling method for barrel structure - Google Patents

Finite element parametric modeling method for barrel structure Download PDF

Info

Publication number
CN107679341B
CN107679341B CN201711034851.8A CN201711034851A CN107679341B CN 107679341 B CN107679341 B CN 107679341B CN 201711034851 A CN201711034851 A CN 201711034851A CN 107679341 B CN107679341 B CN 107679341B
Authority
CN
China
Prior art keywords
finite element
barrel
node
model
file
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.)
Active
Application number
CN201711034851.8A
Other languages
Chinese (zh)
Other versions
CN107679341A (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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201711034851.8A priority Critical patent/CN107679341B/en
Publication of CN107679341A publication Critical patent/CN107679341A/en
Application granted granted Critical
Publication of CN107679341B publication Critical patent/CN107679341B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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]

Abstract

The invention provides a finite element parametric modeling method for a barrel structure, which comprises the following steps: establishing a three-dimensional geometric model of the straight cylindrical barrel through finite element software; dividing a finite element mesh into the geometric model; an execution log file automatically generated by finite element CAE software is changed into a Python script program file rpy to form a Python script language, and the parameters of the barrel structure are replaced by variables to realize parametric modeling; exporting an Inp file of a finite element model of the straight cylindrical barrel; reading the node coordinates of each layer of unit along the barrel axis contained in the Inp file by using a Python language to form a barrel finite element mesh model; and updating Inp, and importing the Inp into commercial finite element CAE software, namely generating a barrel finite element mesh model.

Description

Finite element parametric modeling method for barrel structure
Technical Field
The invention relates to a computer simulation technology, in particular to a barrel structure finite element parametric modeling method.
Background
The barrel is the basic component that provides the ballistic process within the gun, withstanding extremely high powder gas pressures. During the process of launching the projectile, a large amount of high-temperature and high-pressure gas is generated instantly when gunpowder burns to push the projectile to move along the axis direction of the barrel, in the process, the projectile and the rifling of the barrel are mutually extruded, and the belt material undergoes high transient large deformation, friction and other processes. The movement of the projectile in the bore is a highly non-linear problem. The research on the initial motion process in the projectile chamber has important significance for researching the motion rule of the projectile. With the development of science and technology, finite element methods gradually become mainstream methods for studying such problems. In the finite element model, the structure and the grid precision of the barrel finite element grid are particularly important for the accuracy of the simulation calculation result.
The inner bore structure of the rifling barrel is complex and is loaded badly. When using finite element analysis, the barrel rifling geometry is discretized. The general barrel is provided with a plurality of rifling lines, the geometrical size of the cross section of each rifling line is smaller than the caliber of the barrel, and the rifling lines rotate spirally around the axis of the barrel to influence the establishment of a finite element mesh of the barrel and the division of excellent units.
In the technical field of finite element analysis, the whole analysis time occupied by pretreatment modeling is more than 60%, and the conventional modeling technical means is to establish a simplified barrel geometric model in three-dimensional modeling software according to a design drawing, then introduce the simplified barrel geometric model into the pretreatment software for grid division, and finally introduce the divided grid model into finite element software for simulation calculation. The mutual introduction of data among various software not only causes complex operation, but also easily causes data loss, and takes longer time, and the precision of the grid model is poorer. Furthermore, when the pre-processing is performed for barrel mesh division, the quality of the divided mesh is not very good due to the complexity of the barrel structure.
Disclosure of Invention
The invention provides a finite element parametric modeling method for a barrel structure, which can realize rapid and accurate parametric modeling of a barrel medicine chamber, a slope chamber structure, a guide part and the outer contour of the barrel through secondary development.
The technical scheme for realizing the purpose of the invention is as follows: a finite element parametric modeling method for a barrel structure comprises the following steps:
step 1, establishing a three-dimensional geometric model of a straight cylindrical barrel through finite element software;
step 2, dividing a finite element mesh into the geometric model;
step 3, automatically generating an execution log file by finite element CAE software, changing a filename suffix of rpy into a Python script program file of py, forming a Python script language, and replacing the structure parameters of the body pipe by using variables to realize parametric modeling;
step 4, exporting an Inp file of the finite element model of the straight cylindrical barrel;
step 5, reading the node coordinates of each layer of unit along the barrel axis contained in the Inp file by using Python language to form a barrel finite element mesh model;
and 6, updating Inp, and importing the Inp into commercial finite element CAE software, namely generating a barrel finite element mesh model.
By adopting the method, the specific process of the step 5 is as follows:
step 5.1, positioning to a rifling initial part unit by using the same z-direction coordinates of front and rear end surface nodes of each layer of unit of the barrel, taking an end surface node close to the gun tail direction as a source surface node, and taking the other end surface node as a terminal surface node;
step 5.2, taking the source surface node as a fixed point, taking the terminal surface node as an offset node, and taking the source surface node to the terminal surface node as the barrel axial direction;
step 5.3, carrying out any angle rotation deviation on one node (x, y) according to the barrel winding degree and the rifling form to obtain a new coordinate (x ', y')
Figure BDA0001450223690000021
And 5.4, the node offset of the next layer of unit is performed on the basis of the unit of the previous layer, and the operations are repeated, so that the node offset of the whole body pipe can be completed, and a body pipe finite element mesh model is formed.
The method can quickly, effectively and accurately establish the finite element mesh model of the barrel structure, and can carry out parametric modeling on the barrel structure so as to establish the finite element mesh model of the barrel with any different calibers.
The invention is further described below with reference to the accompanying drawings.
Drawings
FIG. 1 is a schematic flow diagram of the process of the present invention.
Figure 2 is a schematic view of a barrel rifling configuration.
FIG. 3 is a schematic diagram of a straight cylindrical barrel mesh finite element model.
FIG. 4 is a schematic diagram of cell layer identification.
Fig. 5 is a schematic diagram of the cell node offset principle.
Fig. 6 is a schematic diagram of node offset calculation.
FIG. 7 is a schematic diagram of a finite element mesh model of a 76mm barrel.
Detailed Description
With reference to fig. 1 to 6, a method for parametrizing and modeling a barrel structure includes the following steps:
step 1, according to the actual size of the barrel bore structure, in ABAQUS software, according to the method of a central small circle (the diameter of the central small circle is the width of a female line) as shown in figure 2, establishing the rifling section of the gun barrel. And the section of the barrel is stretched to form the straight cylindrical barrel.
And 2, dividing a finite element mesh for the straight cylindrical artillery barrel, as shown in figure 3.
And 3, in the process of generating the straight cylindrical barrel grid, automatically generating an execution log file (. rpy) by ABAQUS software, finding a corresponding statement in the file by almost each step of CAE operation, changing the suffix of the file name into a Python script program file (.py) to form a Python script language, and replacing the barrel structure parameter by a variable to facilitate parametric modeling.
Step 4, generating an Inp file containing units and nodes;
and 5, reading the node coordinates of each layer of unit along the axis of the barrel contained in the Inp file by using a Python language to form a barrel finite element mesh model.
For example: as shown in fig. 4, the same z-direction coordinates of the front and rear end face nodes of each layer unit of the barrel are used to position the barrel to the rifling start unit (assuming that the 3 rd layer unit is the rifling start unit), the end face node near the tail direction is used as the source face node, and the other end face node is used as the end face node. After the source node and the terminal node of the unit are identified, the source node is used as a fixed point, and the terminal node is used as an offset node. As shown in fig. 5, for example, in this unit, nodes 100, 101, 102, 103 are source nodes, nodes 200, 201, 202, 203 are terminal nodes, and the source to terminal nodes are barrel axial directions. When the node is shifted, taking the node 203 as an example, the node 203 is one of the terminal surface nodes of the unit and is also one of the source surface nodes of another unit, the node is shifted according to the barrel winding degree and the rifling form, and after rotating by any angle according to a point 203(x, y) on the same circle, the new coordinates 203(x ', y') after the shift can be obtained by the formula (1), that is, the node shift of the layer unit is completed, as shown in fig. 6. Node deviation of the next layer of units is carried out on the basis of the units of the previous layer, and the operations are repeated, so that node deviation of the whole body pipe can be completed, and a body pipe finite element mesh model is formed. The method has the advantages that the barrel finite element mesh model is modeled on the basis of the finite element unit nodes, the accuracy of the established model is high, time period is consumed for modeling, geometric distortion and data loss of the model cannot be generated compared with the traditional modeling method, and the barrel finite element mesh model with higher accuracy can be provided for correlation related to barrel finite element calculation.
Figure BDA0001450223690000041
And 6, importing the modified node coordinates into a new Inp file to form a barrel finite element mesh model.
FIG. 7 shows a 76mm barrel finite element mesh model generated using this method.

Claims (1)

1. A finite element parametric modeling method for a barrel structure is characterized by comprising the following steps:
step 1, establishing a three-dimensional geometric model of a straight cylindrical barrel through finite element software;
step 2, dividing a finite element mesh into the geometric model;
step 3, automatically generating an execution log file by finite element CAE software, changing a filename suffix of rpy into a Python script program file of py, forming a Python script language, and replacing the structure parameters of the body pipe by using variables to realize parametric modeling;
step 4, exporting an Inp file of the finite element model of the straight cylindrical barrel;
step 5, reading the node coordinates of each layer of unit along the barrel axis contained in the Inp file by using Python language to form a barrel finite element mesh model;
step 6, updating the Inp file, and importing the Inp file into commercial finite element CAE software, namely generating a barrel finite element mesh model;
the specific process of the step 5 is as follows:
step 5.1, positioning to a rifling initial part unit by using the same z-direction coordinates of front and rear end surface nodes of each layer of unit of the barrel, taking an end surface node close to the gun tail direction as a source surface node, and taking the other end surface node as a terminal surface node;
step 5.2, taking the source surface node as a fixed point, taking the terminal surface node as an offset node, and taking the source surface node to the terminal surface node as the barrel axial direction;
step 5.3, carrying out any angle rotation deviation on one node (x, y) according to the barrel winding degree and the rifling form to obtain a new coordinate (x ', y')
x ′ = x c o s ( θ ) - y s i n ( θ ) y ′ = x s i n ( θ ) + y c o s ( θ ) ;
Figure DEST_PATH_IMAGE001
θ is the angle of rotation of (x, y);
and 5.4, repeating the steps 5.1 to 5.3 on the basis of the node offset of the unit of the next layer and the unit of the previous layer, so that the node offset of the whole body pipe can be completed, and a body pipe finite element mesh model is formed.
CN201711034851.8A 2017-10-30 2017-10-30 Finite element parametric modeling method for barrel structure Active CN107679341B (en)

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 CN107679341A (en) 2018-02-09
CN107679341B true CN107679341B (en) 2021-03-02

Family

ID=61142397

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711034851.8A Active 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)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
CN111125956B (en) * 2019-12-23 2023-08-22 奇瑞汽车股份有限公司 Automobile finite element model assembling method
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
CN113792386B (en) * 2021-09-27 2022-08-26 中国第一汽车股份有限公司 Finite element automatic modeling method for suspension system parts
CN114970196B (en) * 2022-06-14 2022-12-09 北京科技大学 Method and device for synchronously and quickly generating model and structured grid of metallurgical container

Citations (6)

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

Patent Citations (6)

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

* Cited by examiner, † Cited by third party
Title
身管内膛参数化模型及其磨损有限元模型的生成方法;丁传俊 等;《兵工学报》;20161231;第37卷(第12期);第2212-2219页 *

Also Published As

Publication number Publication date
CN107679341A (en) 2018-02-09

Similar Documents

Publication Publication Date Title
CN107679341B (en) Finite element parametric modeling method for barrel structure
CN109446601B (en) Uncertain optimization method for projectile initial disturbance
Ding et al. A mesh generation method for worn gun barrel and its application in projectile-barrel interaction analysis
CN111105503B (en) Method for determining explosive-loading combustion surface of solid rocket engine
CN110345392A (en) The localization method and device of oil pipeline leak source
CN109697309B (en) Quick acquisition method for high-speed impact extrusion resistance of projectile
Short et al. Steady detonation propagation in a circular arc: a detonation shock dynamics model
CN108459993B (en) Complex high-dimensional system optimization method based on rapid peak-tracking sampling
CN112182955A (en) Construction method and application of shield general type segment point selection model
CN112784458B (en) Bullet and gun coupled finite element parameterized modeling method
CN111046571A (en) Method for quickly determining symmetrical explosive-loading combustion surface of solid rocket engine
CN107256304B (en) Method for evaluating reliability of bullet head emission intensity
Zhao et al. Optimization of the aerodynamic configuration of a tubular projectile based on blind kriging
CN115619930A (en) Combustion surface estimation method based on parameterized volume Boolean operation
CN110276142A (en) The emulation mode of plasma density Conduce Disciplinarian mathematical model under high pressure conditions
CN111931404A (en) Bore complex surface contact collision response prediction method based on self-optimization CNN
Leonhardt et al. Analysis of the effect of bore centerline on projectile exit conditions in small arms
CN109446541A (en) A kind of method of body diamond shape cutting finite element grid modeling
Valko et al. Simulation of an Air Shock Wave Using the Equations of State for the Jones–Wilkins–Lee Detonation Products
Zhou et al. Research on parametric modeling optimization of trajectory correction projectile based on fvm
Milukov et al. Efficiency improvement of technological preparation of power equipment manufacturing
Wadibhasme Implementation of dynamicMesh in OpenFOAM
Song et al. Combined size and shape optimization of structures with DOE, RSM and GA
Xie et al. Multi-Stage Multidisciplinary Design Optimization Method for Enhancing Complete Artillery Internal Ballistic Firing Performance.
Ren et al. Research on Simulation System and Simulation Technology of the Guided Ammunition

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