CN110334423B - Rapid calculation method for axial strain of hole core - Google Patents

Rapid calculation method for axial strain of hole core Download PDF

Info

Publication number
CN110334423B
CN110334423B CN201910557156.2A CN201910557156A CN110334423B CN 110334423 B CN110334423 B CN 110334423B CN 201910557156 A CN201910557156 A CN 201910557156A CN 110334423 B CN110334423 B CN 110334423B
Authority
CN
China
Prior art keywords
finite element
hole
straight beam
axial strain
preset
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
CN201910557156.2A
Other languages
Chinese (zh)
Other versions
CN110334423A (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.)
AVIC Aircraft Strength Research Institute
Original Assignee
AVIC Aircraft Strength Research Institute
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 AVIC Aircraft Strength Research Institute filed Critical AVIC Aircraft Strength Research Institute
Priority to CN201910557156.2A priority Critical patent/CN110334423B/en
Publication of CN110334423A publication Critical patent/CN110334423A/en
Application granted granted Critical
Publication of CN110334423B publication Critical patent/CN110334423B/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]
    • 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)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Complex Calculations (AREA)

Abstract

The application belongs to the technical field of structural strength, and particularly relates to a method for rapidly calculating axial strain of a hole core, which comprises the following steps: classifying the calculation models according to the structural shapes of the straight beam and the preset holes on the straight beam, and respectively establishing parameterized finite element models; determining the finite element node number and the unit number of a preset hole and a hole core thereof; establishing a series of finite element card data files corresponding to the hole cores of the preset holes according to the preset depth increment, and carrying out finite element calculation on the series of finite element card data files by adopting batch processing; extracting axial strain data from the finite element calculation result file by taking the unit number of the hole core as a criterion; and drawing an axial strain curve according to the axial strain data. The rapid calculation method for the axial strain of the hole core is rapid and practical, and has the advantages of clear physical concept, simple operation method and good implementation effect.

Description

Rapid calculation method for axial strain of hole core
Technical Field
The application belongs to the technical field of structural strength, and particularly relates to a rapid calculation method for axial strain of a hole core.
Background
In the technical field of engineering, the problem of axial strain calculation of a hole core is frequently encountered, the traditional calculation method is to build a fixed finite element model according to the structural form and the size of the hole core to perform first calculation, then build a finite element model again according to the change of the depth of the hole core to perform calculation, and the calculation result of each time must be read in a finite element result cloud picture display interface. This method is time consuming and laborious, requiring a repetitive series of calculations as described above if the core structure and dimensions change. In view of the above, a fast calculation method for axial strain of a hole core is proposed and realized.
Disclosure of Invention
In order to solve at least one of the above technical problems, the present application provides a method for rapidly calculating axial strain of a hole core.
The application discloses a method for rapidly calculating axial strain of a hole core, which comprises the following steps:
classifying calculation models according to the structural shapes of the straight beam and the preset holes on the straight beam, and respectively establishing parameterized finite element models;
step two, determining the finite element node number and the unit number of the preset hole and the hole core of the preset hole;
step three, establishing a series of finite element card data files corresponding to the hole cores of the preset holes according to preset depth increment, and carrying out finite element calculation on the series of finite element card data files by adopting batch processing;
step four, extracting axial strain data from the finite element calculation result file by taking the unit number of the hole core as a criterion;
and fifthly, drawing an axial strain curve according to the axial strain data.
According to at least one embodiment of the application, the straight beam has a given length L, a given width W and a given height H, and the constraint form is that the straight beam is fixedly supported at a position L/10 away from the left end of the straight beam, and is fixedly supported at a position L/10 away from the right end of the straight beam, and the load is that the two end faces of the straight beam are respectively stressed by pressure P; the preset hole is arranged at a position which is far from the right end fulcrum L/5 of the straight beam, and comprises a round hole, a round end hole and a rectangular hole.
According to at least one embodiment of the present application, in the step one, the step of establishing parameterized finite element models includes:
step 1.1, inputting parameters;
step 1.2, mesh dissection.
According to at least one embodiment of the present application, the parameter input step includes:
step 1.1.1, inputting the outline dimension of a preset hole structure;
step 1.1.2, inputting load data;
step 1.1.3, inputting physical property parameters;
step 1.1.4, inputting the center coordinates and the radius of the round hole core;
step 1.1.5, inputting the position and the radius of the round end hole core;
and 1.1.6, inputting the center coordinates and the length and width dimensions of the rectangular hole core.
According to at least one embodiment of the present application, the meshing step comprises:
step 1.2.1, dividing the front surface of the straight beam into three areas of a hole, a hole core and the rest part of the front surface of the straight beam;
step 1.2.2, plane mesh dissection;
step 1.2.3, stretching the formed body according to depth increment;
step 1.2.4, deleting the plane grid;
step 1.2.5, reordering finite element nodes and units;
step 1.2.6, merging the coincident points.
According to at least one embodiment of the present application, in the third step, the step of creating a series of finite element card data files corresponding to the predetermined Kong Kongxin according to the predetermined depth increment includes:
step 3.1, constraint and load are applied;
step 3.2, applying physical parameters on the units;
step 3.3, forming finite element card data;
and 3.4, forming a series of finite element card data by taking the serial number depth increment of the hole core as a criterion.
The application has at least the following beneficial technical effects:
the rapid calculation method for the axial strain of the hole core is rapid and practical, and has the advantages of clear physical concept, simple operation method and good implementation effect.
Drawings
FIG. 1 is a computational flow diagram of a method for rapidly calculating axial strain of a core of the present application;
2a, 2b, 2c in FIG. 2 are parameterized finite element models of rectangular holes of opinion for circular holes, round-ended holes, respectively;
FIG. 3 is a graph showing three axial strain curves of a core of the same area obtained by an embodiment of the method for rapidly calculating axial strain of a core of the present application;
fig. 4 is an axial strain curve of a circular core with different radii according to the depth of the core, obtained in another embodiment of the method for rapidly calculating the axial strain of the core.
Detailed Description
In order to make the purposes, technical solutions and advantages of the implementation of the present application more clear, the technical solutions in the embodiments of the present application will be described in more detail below with reference to the accompanying drawings in the embodiments of the present application. In the drawings, the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The described embodiments are some, but not all, of the embodiments of the present application. The embodiments described below by referring to the drawings are exemplary and intended for the purpose of explaining the present application and are not to be construed as limiting the present application. All other embodiments, which can be made by one of ordinary skill in the art based on the embodiments herein without making any inventive effort, are intended to be within the scope of the present application. Embodiments of the present application are described in detail below with reference to the accompanying drawings.
In the description of the present application, it should be understood that the terms "center," "longitudinal," "lateral," "front," "rear," "left," "right," "vertical," "horizontal," "top," "bottom," "inner," "outer," and the like indicate orientations or positional relationships based on the orientations or positional relationships illustrated in the drawings, merely to facilitate description of the present application and simplify the description, and do not indicate or imply that the device or element being referred to must have a specific orientation, be configured and operated in a specific orientation, and therefore should not be construed as limiting the scope of protection of the present application.
The method for rapidly calculating the axial strain of the core of the present application is described in further detail below with reference to fig. 1-4.
For the engineering straight beam with given length L, width W and height H, a hole is dug at a position which is far from a right end fulcrum L/5, the hole forms comprise a circular hole, a circular end hole and a rectangular hole, the constraint forms are fixedly supported at a position which is far from the left end L/10 of the straight beam, fixedly supported at a position which is far from the right end L/10, and the load is that the two end faces of the straight beam are respectively stressed by pressure P.
Because the structural form corresponding to the axial strain calculation of the hole core is relatively simple, the difference is mainly reflected in the difference of the size and the position and the difference of the shape of the hole. According to this feature, we can implement the series of designs by programming, freeing the designer from relatively repetitive manual means, enabling the calculation cycle to be greatly shortened and avoiding the drawbacks of manual error-prone. As long as the designer lifts the size of the structure and the change increment of the hole depth, the analysis work and the curve drawing are automatically completed by a computer.
Specifically, as shown in fig. 1, the application discloses a method for quickly calculating axial strain of a hole core, which may include the following steps:
classifying the calculation models according to the structural shapes of the straight beam and the preset holes on the straight beam, normalizing the calculation models, and respectively establishing parameterized finite element models. In this embodiment, a parameterized finite element model of rectangular holes can be obtained for the circular holes and the circular end holes shown in fig. 2a, 2b, and 2c, respectively.
Specifically, in this step, establishing the parameterized finite element models, respectively, may further include:
and 1.1, inputting parameters. This step may further comprise:
step 1.1.1, inputting the outline dimension of a preset hole structure;
step 1.1.2, inputting load data;
step 1.1.3, inputting physical property parameters;
step 1.1.4, inputting the center coordinates and the radius of the round hole core;
step 1.1.5, inputting the position and the radius of the round end hole core;
and 1.1.6, inputting the center coordinates and the length and width dimensions of the rectangular hole core.
Step 1.2, mesh dissection. This step may further comprise:
step 1.2.1, dividing the front surface of the straight beam into three areas of a hole, a hole core and the rest part of the front surface of the straight beam;
step 1.2.2, plane mesh dissection;
step 1.2.3, stretching the formed body according to depth increment;
step 1.2.4, deleting the plane grid;
step 1.2.5, reordering finite element nodes and units;
step 1.2.6, merging the coincident points.
Step two, determining the finite element node number and the unit number of the preset hole and the hole core of the preset hole.
And thirdly, establishing a series of finite element card data files corresponding to the hole cores of the preset holes according to the preset depth increment, and carrying out finite element calculation on the series of finite element card data files by adopting batch processing.
Further, the step of creating a series of finite element card data files corresponding to the predetermined Kong Kongxin according to the predetermined depth increment includes:
step 3.1, constraint and load are applied;
step 3.2, applying physical parameters on the units;
step 3.3, forming finite element card data;
and 3.4, forming a series of finite element card data by taking the predetermined depth increment of the hole core number as a criterion.
And step four, extracting axial strain data from the finite element calculation result file by taking the unit number of the hole core as a criterion. Specifically, the present embodiment automatically retrieves the axial strain value of the core through PCL language in a series of analysis result files.
And fifthly, drawing an axial strain curve according to the axial strain data. Specifically, the axial strain curve is drawn by using the ORIGIN commercial drawing software. It should be noted that, finally, axial strain curves corresponding to different depths of holes with different shapes as shown in fig. 3 can be obtained; it is also possible to obtain holes of the same shape with different diameters as shown in fig. 4, corresponding to axial strain curves at different depths.
The foregoing is merely specific embodiments of the present application, but the scope of the present application is not limited thereto, and any changes or substitutions easily conceivable by those skilled in the art within the technical scope of the present application should be covered in the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.

Claims (1)

1. The method for rapidly calculating the axial strain of the hole core is characterized by comprising the following steps of:
classifying calculation models according to the structural shapes of the straight beam and the preset holes on the straight beam, and respectively establishing parameterized finite element models;
step two, determining the finite element node number and the unit number of the preset hole and the hole core of the preset hole;
step three, establishing a series of finite element card data files corresponding to the hole cores of the preset holes according to preset depth increment, and carrying out finite element calculation on the series of finite element card data files by adopting batch processing;
step four, extracting axial strain data from the finite element calculation result file by taking the unit number of the hole core as a criterion;
drawing an axial strain curve according to the axial strain data, wherein the straight beam has a given length L, a given width W and a given height H, the constraint form is that the straight beam is fixedly supported at a position L/10 away from the left end of the straight beam, the straight beam is fixedly supported at a position L/10 away from the right end of the straight beam, and the load is that the two end faces of the straight beam are respectively stressed by pressure P; the preset hole is arranged at the position L/5 away from the fulcrum L at the right end of the straight beam, and comprises a round hole, a round end-shaped hole and a rectangular hole;
in the first step, the step of establishing parameterized finite element models respectively includes:
step 1.1, inputting parameters;
step 1.2, mesh dissection;
the parameter input step comprises the following steps:
step 1.1.1, inputting the outline dimension of a preset hole structure;
step 1.1.2, inputting load data;
step 1.1.3, inputting physical property parameters;
step 1.1.4, inputting the center coordinates and the radius of the round hole core;
step 1.1.5, inputting the position and the radius of the round end hole core;
step 1.1.6, inputting the center coordinates and the length and width dimensions of the rectangular hole core;
the mesh dissection step comprises the following steps:
step 1.2.1, dividing the front surface of the straight beam into three areas of a hole, a hole core and the rest part of the front surface of the straight beam;
step 1.2.2, plane mesh dissection;
step 1.2.3, stretching the formed body according to depth increment;
step 1.2.4, deleting the plane grid;
step 1.2.5, reordering finite element nodes and units;
step 1.2.6, merging the coincident points;
in the third step, the step of creating a series of finite element card data files corresponding to the predetermined Kong Kongxin according to the predetermined depth increment includes:
step 3.1, constraint and load are applied;
step 3.2, applying physical parameters on the units;
step 3.3, forming finite element card data;
and 3.4, forming a series of finite element card data by taking the predetermined depth increment of the hole core number as a criterion.
CN201910557156.2A 2019-06-25 2019-06-25 Rapid calculation method for axial strain of hole core Active CN110334423B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910557156.2A CN110334423B (en) 2019-06-25 2019-06-25 Rapid calculation method for axial strain of hole core

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910557156.2A CN110334423B (en) 2019-06-25 2019-06-25 Rapid calculation method for axial strain of hole core

Publications (2)

Publication Number Publication Date
CN110334423A CN110334423A (en) 2019-10-15
CN110334423B true CN110334423B (en) 2023-06-20

Family

ID=68142731

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910557156.2A Active CN110334423B (en) 2019-06-25 2019-06-25 Rapid calculation method for axial strain of hole core

Country Status (1)

Country Link
CN (1) CN110334423B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107577875A (en) * 2017-09-07 2018-01-12 东南大学 A kind of suture type battenboard parametric Finite Element Modeling Method
CN109325271A (en) * 2018-09-05 2019-02-12 北京理工大学 A kind of carrier rocket structural finite element model parametric modeling method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107577875A (en) * 2017-09-07 2018-01-12 东南大学 A kind of suture type battenboard parametric Finite Element Modeling Method
CN109325271A (en) * 2018-09-05 2019-02-12 北京理工大学 A kind of carrier rocket structural finite element model parametric modeling method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
圆孔密集分布模板的强度损失与应力集中;蔡中义等;《吉林大学学报(工学版)》;20031030(第04期);全文 *
带孔复合材料板和螺栓连接复合材料板的有限元计算分析;高宗战等;《飞机设计》;20031230(第04期);第24-28页 *
腹板开圆孔冷弯卷边槽钢轴压构件畸变屈曲承载力试验及计算方法;姚行友等;《工业建筑》;20160108(第04期);全文 *

Also Published As

Publication number Publication date
CN110334423A (en) 2019-10-15

Similar Documents

Publication Publication Date Title
CN100570639C (en) Use the reverse modeling method of grid data as feature
CN101388118B (en) Three-dimensional shoe last rebuilding method based on scanning and inputting
CN105900145A (en) Method for creating three dimensional lattice structures in computer-aided design models for additive manufacturing
JP2008186440A (en) Topology evolution optimization computing method for structural design
CN109388840B (en) Special programming computer for numerical simulation of metal forming processes with predefined load paths and corresponding grid adjustment schemes
CN110334423B (en) Rapid calculation method for axial strain of hole core
CN113792458B (en) Optimization method and device for finite element triangle mesh
CN109446749B (en) Design method and system suitable for ship jig, electronic terminal and storage medium
CN114861500A (en) Method and system for automatically generating tunnel structure finite element model based on three-dimensional point cloud
CN111415403A (en) Method for establishing virtual experimental digital pavement based on pavement scanning data and storage medium
US8477133B2 (en) Method and apparatus for generating three-dimensional finite element mesh
CN107330179B (en) Method for quickly creating three-dimensional process model with process auxiliary characteristics
CN113538695A (en) Method and device for quickly discretizing complex curved surface with arbitrary boundary and storage medium
CN109359360B (en) Structural stress processing method based on local features
JOSÉ Optimization of elastic structures using boundary elements and a topological-shape sensitivity formulation
CN115170734A (en) Three-dimensional revolution structure reconstruction method and system based on section slices
CN115592002A (en) Gas cap connecting column arrangement method, device, equipment and medium of stamping die
CN107832537B (en) Numerical simulation method for residual stress in complex form
EP2442279A1 (en) Evolvable control point grid layouts for the design optimization of physical bodies
JP4872583B2 (en) Shape model creation device
JP6308544B2 (en) Method for calculating contact surface pressure between two objects and computer program for calculating contact surface pressure between two objects
CN110059409B (en) BIM-based long-line type engineering longitudinal and transverse integration processing method
JP6469555B2 (en) Model design method
CN112509156B (en) UG-based tire pattern rapid editing system and method thereof
CN113779684A (en) Construction method of roadbed slope lattice beam construction model based on Revit

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