CN111814382B - Wavefront recognition method for propagation of non-planar wave in multicellular material - Google Patents

Wavefront recognition method for propagation of non-planar wave in multicellular material Download PDF

Info

Publication number
CN111814382B
CN111814382B CN202010717279.0A CN202010717279A CN111814382B CN 111814382 B CN111814382 B CN 111814382B CN 202010717279 A CN202010717279 A CN 202010717279A CN 111814382 B CN111814382 B CN 111814382B
Authority
CN
China
Prior art keywords
strain
multicellular
octahedral
wavefront
tensor
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
CN202010717279.0A
Other languages
Chinese (zh)
Other versions
CN111814382A (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.)
General Engineering Research Institute China Academy of Engineering Physics
Original Assignee
General Engineering Research Institute China Academy of Engineering Physics
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 General Engineering Research Institute China Academy of Engineering Physics filed Critical General Engineering Research Institute China Academy of Engineering Physics
Priority to CN202010717279.0A priority Critical patent/CN111814382B/en
Publication of CN111814382A publication Critical patent/CN111814382A/en
Application granted granted Critical
Publication of CN111814382B publication Critical patent/CN111814382B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a wavefront recognition method for non-planar wave propagation in a multicellular material, which is based on a grid node displacement field generated by numerical simulation software, adopts a local strain field calculation method based on a least square discrete local deformation gradient, adopts a continuous medium mechanical calculation method to convert Green strain tensor into engineering strain tensor, and determines a non-planar wavefront position by an octahedral equivalent strain calculation method, thereby solving the problem of recognition of the non-planar wave in the multicellular material, providing a more accurate calculation method for recognizing the non-planar wave wavefront under the general application environment that a foam macrostructure presents a complex structure or a load is in a non-single direction, and simultaneously providing a powerful tool for recognizing propagation rules of the non-planar wavefront in the cellular material and laying a research foundation. The invention can be applied to basic scientific research and engineering design of anti-explosion and impact-resistant protection designs such as aerospace, military protection, automobile manufacturing, product packaging and the like.

Description

Wavefront recognition method for propagation of non-planar wave in multicellular material
Technical Field
The invention relates to the technical field of multi-cell material shock wave front recognition, in particular to a wave front recognition method for non-plane waves propagating in multi-cell materials under explosion or impact conditions.
Background
The multicellular material is used as an impact protection buffer material, is widely applied to aerospace, military protection, automobile manufacture, product packaging and the like, and the propagation rule of the impact wave in the multicellular material has basic scientific research significance for antiknock impact protection design and important guiding significance for military and civil protection, and the identification technology of the impact wave front in discontinuous media such as the multicellular material is an important subject for research in the field.
At present, in the technical field of multi-cell material shock wave front recognition, the method is divided into two routes. The first route is an experimental means, as shown in fig. 1, a high-speed image of a multicellular material test piece under dynamic compression is shot through high-speed photography, then analysis is carried out through DIC digital correlation technique, a strain distribution field on the surface of the test piece is obtained, and the position of a wave front is determined according to the strain jump position; the second route is numerical simulation and theoretical calculation, a strain distribution field in a test piece cannot be obtained based on an experimental means, displacement field data at a unit node under an impact condition is obtained by adopting a multi-cell material microscopic numerical simulation method, a Green strain tensor at the node is obtained by adopting a local strain field calculation method based on a least square discrete local deformation gradient, and the strain is converted into strain in an impact direction through continuous medium mechanics, so that a strain jump position is identified, and a wavefront position is determined.
Of the two routes in the technical field of multi-cell material shock wave wavefront recognition, the first experimental route is limited by the fact that high-speed photography can only shoot the dynamic deformation of the outer surface of a test piece, and the internal deformation and the form of a shock wave wavefront cannot be obtained; the second numerical simulation+theoretical calculation route can only identify plane shock wave fronts in the multicellular material, and non-plane waves cannot be identified or interface identification is inaccurate because only positive strain components are considered and shear strain components are not involved. Based on the above, the present invention is proposed to solve the problem of recognizing the non-planar wave inside the multicellular material.
Disclosure of Invention
The present invention aims to solve the above problems by providing a wavefront recognition method for propagation of non-planar waves in multicellular material.
To achieve the above object, the present disclosure provides a wavefront recognition method of propagation of a non-planar wave in a multicellular material, including:
establishing a numerical simulation model of the multicellular material by a microscopic numerical simulation method, and carrying out dynamic numerical simulation under explosion or impact conditions;
obtaining a local strain field of the multicellular material by a local strain field calculation method;
converting the Green strain tensor into an engineering strain tensor by adopting a continuous medium mechanical calculation method;
the octahedral equivalent strain field formed by each node of the material is obtained by an octahedral equivalent strain calculation method, and the non-planar wavefront position is determined by software.
Alternatively, the multicellular material is modeled as a finite element model, pi at a designated macroscopic region of the material 0 The N nucleation points are distributed according to a uniform probability distribution, and the distance delta between the nucleation points is controlled to be not smaller than a given distance delta min
Generating cell edge and cell surface geometric information through MATLAB software, and performing grid division on cell surfaces to obtain a numerical simulation model of the multicellular material in finite element software ABAQUS.
Alternatively, for a multicellular material of regular mesostructure, a finite element model is constructed periodically using a finite element method.
Alternatively, for multicellular materials of random mesostructure, a finite element model is constructed using 3D-Voronoi random mesostructure.
Alternatively, in the multicellular material, discrete local deformation gradients based on least squares are used to describe discrete displacement fields resulting from multicellular material, resulting in local strain fields for the multicellular material.
Alternatively, the positive strain component and the shear strain component at the material node form an engineering strain tensor, and the Green strain tensor is solved by a continuous medium mechanical method.
Alternatively, the engineering strain tensor is mapped from 2 th order to 0 th order, and the octahedral equivalent strain calculation method is adopted to perform 0 th order characterization on the octahedral equivalent strain, so that an octahedral equivalent strain field is formed.
Optionally, according to the octahedral equivalent strain field formed by each node of the material, interpolation is performed on the three-dimensional discrete octahedral equivalent strain field through MATLAB software to obtain a three-dimensional continuous strain field, so that the purpose of identifying the non-plane wave front is achieved.
The invention has the beneficial effects that:
the invention relates to a wavefront recognition method for spreading non-planar waves in a multicellular material, which is based on a microscopic numerical simulation method and a local strain field calculation method, converts Green strain tensor into engineering strain tensor by adopting a continuous medium mechanical calculation method, and determines the recognition method of the equivalent strain field jump position and the non-planar wavefront position by an octahedral equivalent strain calculation method, thereby solving the recognition problem of the non-planar waves in the multicellular material, providing a more accurate calculation method for recognizing the non-planar wave wavefront under the general application environment that a foam macrostructure presents a complex structure or a load is non-unidirectional, and providing a powerful tool for recognizing the spreading rule of the non-planar wave wavefront in the cellular material, and laying a research foundation. The invention can be applied to basic scientific research and engineering design of anti-explosion and impact-resistant protection designs such as aerospace, military protection, automobile manufacturing, product packaging and the like.
Additional aspects and advantages of the invention will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention.
Drawings
The accompanying drawings are included to provide a further understanding of the disclosure, and are incorporated in and constitute a part of this specification, illustrate the disclosure and together with the description serve to explain, but do not limit the disclosure. In the drawings:
FIG. 1 is a schematic diagram of the microstructure geometry of an aluminum foam according to the example;
FIG. 2 is a numerical model of the microstructure of aluminum foam of the example;
FIG. 3 is a cross-sectional view of a numerical simulated equivalent stress of an aluminum foam microstructure according to the example;
FIG. 4 is a mid-plane cross-sectional view of the ε 11 positive strain field under the aluminum foam microstructure local strain method described in the examples;
FIG. 5 is a graph of ε under the local strain of the microstructure of aluminum foam described in the examples eff A mid-plane cross-sectional view of an octahedral equivalent strain field.
Detailed Description
Specific embodiments of the present disclosure are described in detail below with reference to the accompanying drawings. It should be understood that the detailed description and specific examples, while indicating and illustrating the disclosure, are not intended to limit the disclosure.
The invention relates to a wavefront recognition method for spreading non-planar waves in a multicellular material, which is based on a grid node displacement field generated by numerical simulation software, adopts a local strain field calculation method based on a least square discrete local deformation gradient to obtain a Green strain tensor at a node, converts the Green strain tensor into an engineering strain tensor through continuous medium mechanics, adopts an octahedral equivalent strain calculation method to determine an equivalent strain field jump position and a non-planar wavefront position, thereby solving the problem of recognizing the non-planar waves in the multicellular material, providing a more accurate calculation method for recognizing the non-planar wave wavefront under the general application environment that a foam macrostructure presents a complex structure or a load is non-unidirectional, and simultaneously providing a powerful tool for recognizing the propagation rule of the non-planar wavefront in the cellular material and laying a research foundation. The invention can be applied to basic scientific research and engineering design of anti-explosion and impact-resistant protection designs such as aerospace, military protection, automobile manufacturing, product packaging and the like.
The specific implementation process is as follows:
1. mesoscopic numerical simulation method
The multicellular material is divided into a regular type and a random type microstructure, and for the multicellular material with the regular type microstructure, a finite element method can be utilized to directly and periodically build a model. For multicellular materials with irregular microscopic structures, such as foam materials, a finite element model is built by adopting a 3D-Voronoi random microscopic structure, and pi is firstly formed in a designated macroscopic area of the material 0 The N nucleation points are distributed according to a uniform probability distribution, and the distance delta between the nucleation points is controlled to be not smaller than a given distance delta min
δ min =(1-k)·δ 0 (1)
(1) Wherein k is the irregularity in the foam material characterizing the microstructure; delta 0 Is the minimum distance of any nucleation point piece of a decahedron Kelvin cell edge topology structure:
(2) Wherein V is cell Is the volume of a unit cell.
After the nucleation points are generated under the conditions, generating cell edges and cell surface geometric information through a Voronoi function in MATLAB. The cell surface is subjected to grid division by adopting a Hypermesh shell unit, and the foam cell wall thickness is set to be t in finite element software ABAQUS. And obtaining a numerical simulation model of the foam microstructure, and carrying out dynamic numerical simulation under explosion or impact conditions.
2. Local strain field calculation method
In the process of obtaining the non-planar shock wave of the microstructure of the foam material, not only the positive strain component at the node point is considered, but also the shear strain component at the node point is considered. The two components form an engineering strain tensor and are solved by the Green strain tensor through a continuous medium mechanical method.
Characterization of local deformation according to Green strain tensor E in continuous medium mechanics:
(3) Where F is the deformation gradient, I is the unit tensor, and superscript T denotes the tensor. In the multicellular material, discrete local deformation gradient based on the least square method is adopted to describe discrete displacement fields brought by multicellular property, so that local strain fields of the multicellular material are obtained.
Defining the node configuration of the foam material before deformation as omega 0 The deformed node is configured to be the current configuration omega 1 . Node i and its neighboring node j are in configuration Ω 0 And configuration omega 1 The relative displacement vectors in (a) are U respectively ij and uij
U ij =X i -X j (4)
u ij =x i -x j (5)
wherein ,Xi and Xj Is the configuration omega of node i and adjacent node j 0 Position vector, x i and xj Is the configuration omega of node i and adjacent node j 1 Is included in the position vector. Based on the assumption of local uniform deformation in continuous medium mechanics, the existence of a unique deformation gradient F between a node i and an adjacent node j will form a configuration omega 0 Mapping to configuration Ω 1 The method comprises the following steps:
u ij =F·U ij
(6)
because of the discrete type and the complexity of the microstructure of the multicellular material, the optimal deformation gradient F of the node i is obtained by the least square method of mapping the deformation gradients of all the configurations of the node i and other nodes in a certain range of the field of the node i i
If the equation and hence the determinant of the second term is 0, then it is assumed that node i is in a zero strain state at this time. Wherein the domain radius R of node i is set to 0.5R, R being the average radius of the cell.
3. Method for converting Green strain tensor E into engineering strain tensor E1
The Green strain tensor E obtained by the formula (7) and the formula (3) is:
wherein ,est Is a component of the Green strain tensor E, and subscripts s and t take 1,2,3, respectively.
For engineering strain tensor E 1
wherein ,εst Is a component of the Green strain tensor E, and subscripts s and t take 1,2,3, respectively. For engineering strain tensor E 1 The mid-diagonal upper component ε ss The component of the Green strain tensor E can be expressed as:
engineering strain tensor E 1 The mid-off-diagonal upper component ε st The component of the Green strain tensor E can be expressed as:
4. octahedral equivalent strain calculation method
E obtained above 1 Representing the strain state at the node, but its attribute is the engineering strain tensor, which is highly dependent on the type and direction of the load when mapped to the foam microstructure. For example in the simplest case: the foam microstructure is a regular hexahedron with normal lines of each surface pointing to the x, y or z axis, and the load is in a single direction and along the x axisIf the foam microstructure is considered to be macroscopically uniform and continuous, the shock wave is a plane wave in the foam microstructure, in which case ε is obtained from formula (10) 11 The positive strain field, through strain jump position identification, can confirm the shock wave wavefront position. However, if the foam microstructure exhibits a non-regular hexahedral structure or the load is not unidirectional, the shock wave in the foam microstructure is a non-planar wave and the positive strain field determined by equation (10) cannot accurately identify the shock wave wavefront position from the strain relief position due to the presence of shear strain. The invention provides a method for adopting octahedral equivalent strain, which considers the influence of positive strain and shear strain and tenses engineering strain E 1 Complete consideration is made. Octahedral equivalent strain ε eff The method comprises the following steps:
5. non-plane wave front recognition method
According to the octahedral equivalent strain field formed by each node obtained through calculation, the ScatteredInterpolant function in MATLAB software is used for interpolating the three-dimensional discrete octahedral equivalent strain field to obtain a three-dimensional continuous strain field, so that the purpose of identifying the non-planar wave front is achieved.
Examples:
the macroscopic foam has a region size of 80mm×20mm, 600 nucleation points distributed according to uniform probability are randomly scattered in the region, the irregularity k of the microstructure is set to 0.2, and the minimum given distance delta between the nucleation points can be obtained according to the formula (1) and the formula (2) min About 3.3mm, the microstructure geometry of the aluminum foam is shown in figure 1. Grid division is carried out by adopting an S4R shell unit, the cell wall thickness of a foam material is set to be 0.045mm, an aluminum alloy is adopted as a foam material base material, an elastic ideal shaping material model is adopted, and the density is 2.77g/cm 3 Young's modulus of 69GPa, poisson's ratio of 0.3 and yield strength of 170MPa. In order to generate non-planar shock wave, a spherical shell is adopted at the impact end, a flat plate is adopted at the supporting end, and an explosion shock wave with characteristic time of 50 mu s and amplitude of 20MPa acts on the inner part of the spherical shellThe numerical model of the concave surface and the microstructure of the aluminum foam is shown in figure 2. The middle section of the equivalent stress distribution obtained after the numerical calculation is shown in fig. 3, and according to the principle that the front and back stress, strain and speed of the wavefront can generate abrupt transition, the abrupt transition interface cannot be identified, namely the position of the wavefront cannot be directly identified. After the displacement field information is output, the epsilon under the condition of positive strain is only considered by adopting a local strain method, namely the formulas (3) - (10) 11 A positive strain field, wherein the plane cross section is shown in figure 4; epsilon under the conditions of positive strain and shear strain can be considered simultaneously by adopting a local strain method and an octahedral equivalent strain calculation method, namely the formulas (3) - (12) eff The positive strain field, wherein the plane section is shown in figure 5, can more accurately identify the jump interface, thereby more accurately acquiring the position and morphological characteristics of the wave front.
The invention is based on the existing microscopic numerical simulation technology and local strain field calculation method, adopts a continuous medium mechanical calculation method to convert Green strain tensor into engineering strain tensor, and provides an identification method for determining the equivalent strain field jump position and the non-planar wavefront position through an octahedral equivalent strain calculation method, thereby solving the identification problem of non-planar waves in a multicellular material, providing a more accurate calculation method for identifying non-planar wave fronts in general application environments such as complex structures of foam macrostructures or non-single load directions, and providing a powerful tool for recognizing propagation rules of non-planar wave fronts in cellular materials. The invention can be applied to basic scientific research and engineering design of anti-explosion and impact-resistant protection designs such as aerospace, military protection, automobile manufacturing, product packaging and the like.
In addition, the specific features described in the foregoing embodiments may be combined in any suitable manner, and in order to avoid unnecessary repetition, the present disclosure does not further describe various possible combinations.
Moreover, any combination between the various embodiments of the present disclosure is possible as long as it does not depart from the spirit of the present disclosure, which should also be construed as the disclosure of the present disclosure.

Claims (3)

1. A method for wavefront recognition of a non-planar wave propagating in a multicellular material, comprising:
establishing a numerical simulation model of a multicellular material by a microscopic numerical simulation method, and carrying out dynamic numerical simulation under explosion or impact conditions, wherein the multicellular material is divided into a regular microscopic structure and a random microscopic structure; for the multicellular material with an irregular microscopic structure, a finite element model is constructed by adopting a 3D-Voronoi random microscopic structure, and the method specifically comprises the following steps:
in the designated macro-area of the materialΠ 0 Distributing N nucleation points according to uniform probability distribution and controlling the distance between the nucleation pointsδNot less than a given distanceδ min
(1)
(1) in the formula ,kcharacterizing irregularities of the microstructure in the multicellular material that is of an irregular microstructure;δ 0 is the minimum distance of any nucleation point piece of a decahedron Kelvin cell edge topology structure:
(2)
(2) in the formula ,V cell is the volume of a unit cell;
generating nucleation points under the above conditions, generating cell edge and cell surface geometric information through a Voronoi function in MATLAB, meshing cell surfaces by adopting Hypermesh shell units, setting the cell wall thickness of a multicellular material with an irregular microscopic structure as t in finite element software ABAQUS, obtaining a numerical simulation model of the multicellular material microscopic structure with the irregular microscopic structure, and carrying out dynamic numerical simulation under explosion or impact conditions;
obtaining a local strain field of the multicellular material by a local strain field calculation method; the method comprises the following steps:
acquiring a positive strain component at a node and a shear strain component at the node in a non-planar shock wave process of a multi-cellular material microstructure of an irregular microstructure, wherein the two components form an engineering strain tensor and are solved by a Green strain tensor through a continuous medium mechanical method;
characterization of local deformation according to Green strain tensor E in continuous medium mechanics:
(3)
(3) in the formula ,Fis a deformation gradient which is a gradient of the deformation,Iis the tensor of the unit, superscriptTMeans for representing tensors, in the multicellular material, using a discrete local deformation gradient based on a least square method to describe a discrete displacement field resulting from multicellular properties, resulting in a local strain field for the multicellular material;
converting the Green strain tensor into an engineering strain tensor by adopting a continuous medium mechanical calculation method;
obtaining an octahedral equivalent strain field formed by each node of the material by an octahedral equivalent strain calculation method, and determining the position of a non-planar wavefront by software; the method comprises the following steps:
the engineering strain tensor obtained above represents the strain state at the node, but the attribute is the engineering strain tensor, and the mapping to the multi-cell material microstructure of the irregular microstructure is highly dependent on the type and direction of the load; taking positive strain and shear strain effects into consideration, carrying out complete consideration on engineering strain tensors, mapping the engineering strain tensors from 2-order to 0-order, carrying out 0-order characterization on the octahedral equivalent strain by adopting an octahedral equivalent strain calculation method, thereby forming an octahedral equivalent strain field and an octahedral equivalent strainε eff The method comprises the following steps:
(12)
wherein ,ε st is a component of the Green strain tensorSubscript ofsAndttaking 1,2 and 3 respectively.
2. The method of wavefront recognition of a non-planar wave propagating in a multicellular material of claim 1 wherein: for the multicellular material with a regular microscopic structure, a finite element model is periodically constructed by using a finite element method.
3. The method of wavefront recognition of a non-planar wave propagating in a multicellular material of claim 1 wherein: according to the octahedral equivalent strain fields formed by the nodes of the material, interpolation is carried out on the three-dimensional discrete octahedral equivalent strain fields through MATLAB software to obtain three-dimensional continuous strain fields, so that the purpose of identifying non-plane wave fronts is achieved.
CN202010717279.0A 2020-07-23 2020-07-23 Wavefront recognition method for propagation of non-planar wave in multicellular material Active CN111814382B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010717279.0A CN111814382B (en) 2020-07-23 2020-07-23 Wavefront recognition method for propagation of non-planar wave in multicellular material

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010717279.0A CN111814382B (en) 2020-07-23 2020-07-23 Wavefront recognition method for propagation of non-planar wave in multicellular material

Publications (2)

Publication Number Publication Date
CN111814382A CN111814382A (en) 2020-10-23
CN111814382B true CN111814382B (en) 2023-09-22

Family

ID=72862468

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010717279.0A Active CN111814382B (en) 2020-07-23 2020-07-23 Wavefront recognition method for propagation of non-planar wave in multicellular material

Country Status (1)

Country Link
CN (1) CN111814382B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106777813A (en) * 2017-01-16 2017-05-31 中原工学院 A kind of method for improving and controlling Multi-halfwaves Filters energy absorption efficiency
US9915522B1 (en) * 2013-06-03 2018-03-13 Kla-Tencor Corporation Optimized spatial modeling for optical CD metrology
CN108694290A (en) * 2018-06-05 2018-10-23 东北大学 A kind of Soft Tissue Deformation method of the finite element model based on Octree cube
CN110892260A (en) * 2017-06-09 2020-03-17 埃因霍温科技大学 Shear wave viscoelastic imaging using local system identification
CN110948878A (en) * 2019-11-19 2020-04-03 中国科学技术大学 Explosion gradient-imitating polycell bullet and design and manufacture method thereof
CN111125963A (en) * 2020-01-06 2020-05-08 深圳拳石科技发展有限公司 Numerical simulation system and method based on Lagrange integral point finite element

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9915522B1 (en) * 2013-06-03 2018-03-13 Kla-Tencor Corporation Optimized spatial modeling for optical CD metrology
CN106777813A (en) * 2017-01-16 2017-05-31 中原工学院 A kind of method for improving and controlling Multi-halfwaves Filters energy absorption efficiency
CN110892260A (en) * 2017-06-09 2020-03-17 埃因霍温科技大学 Shear wave viscoelastic imaging using local system identification
CN108694290A (en) * 2018-06-05 2018-10-23 东北大学 A kind of Soft Tissue Deformation method of the finite element model based on Octree cube
CN110948878A (en) * 2019-11-19 2020-04-03 中国科学技术大学 Explosion gradient-imitating polycell bullet and design and manufacture method thereof
CN111125963A (en) * 2020-01-06 2020-05-08 深圳拳石科技发展有限公司 Numerical simulation system and method based on Lagrange integral point finite element

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
3D multicellular model of shock wavecell interaction;Dongli Li 等;Acta Biomaterialia;第77卷;282-291 *
基于Livelink数据交互的光纤布拉格光栅-力学结构多物理场感知仿真方法;刘杉 等;光学学报;第40卷(第11期);54-62 *
基于晶体塑性的难变形材料不均匀变形多尺度建模研究进展;詹梅 等;塑性工程学报;第25卷(第01期);1-14 *
基于细观离散元的混凝土端部效应分析;康政 等;哈尔滨工业大学学报;第45卷(12);94-98 *
星形节点周期性蜂窝结构的面内动力学响应特性研究;韩会龙 等;振动与冲击;第36卷(第23期);223-231 *
随机Kelvin胞体开孔泡沫结构模型边界条件及应用;齐敏菊 等;系统仿真学报;第28卷(第01期);190-199 *
颗粒金属材料冲击压缩细观数值模拟;乔良 等;高压物理学报;第27卷(第06期);863-871 *
高速冲击下多胞材料细观结构中波传播规律;李伟 等;装备环境工程;第18卷(第5期);87-93 *

Also Published As

Publication number Publication date
CN111814382A (en) 2020-10-23

Similar Documents

Publication Publication Date Title
Natarajan et al. Numerical integration over arbitrary polygonal domains based on Schwarz–Christoffel conformal mapping
CN108345741B (en) Two-dimensional thermal deformation and thermal stress analysis method based on non-grid RKPM anisotropic material
CN109918712B (en) Calculation method for solving composite stress intensity factor based on p-type finite element method and contour integral method
Daxini et al. A review on recent contribution of meshfree methods to structure and fracture mechanics applications
Kim et al. Analysis of hydroelasticity of floating shiplike structure in time domain using a fully coupled hybrid BEM-FEM
CN110162821B (en) Method for calculating bird-hit high-speed rotating engine blade
He et al. Development of an efficient algorithm to analyze the elastic wave in acoustic metamaterials
Chen et al. Free vibration of moderately thick functionally graded plates by a meshless local natural neighbor interpolation method
CN103699714A (en) Flexible object real-time cutting simulation method based on finite element and meshless coupling
CN106940739A (en) It is a kind of to vibrate the method for quick predicting influenceed on wing conformal phased array antenna electrical property
CN112768006B (en) Method, medium, and apparatus for establishing a star-honeycomb constitutive relation
Zhang et al. MPS-FEM coupled method for fluid–structure interaction in 3d dam-break flows
Lv et al. Element differential method for free and forced vibration analysis for solids
Ma et al. Nonlinear waves in a floating thin elastic plate, predicted by a coupled SPH and FEM simulation and by an analytical solution
Cui et al. Analysis of Mindlin–Reissner plates using cell-based smoothed radial point interpolation method
Sladek et al. Elastic analysis in 3D anisotropic functionally graded solids by the MLPG
CN111814382B (en) Wavefront recognition method for propagation of non-planar wave in multicellular material
Wang et al. Analysis of fatigue crack propagation of an orthotropic bridge deck based on the extended finite element method
Wu et al. The development of hybrid ES-FE-SEA method for mid-frequency vibration analysis of complex built-up structure
Singh et al. Vibroacoustic response of mode localized thin functionally graded plates using physical neutral surface
Xi et al. An efficient hybrid collocation scheme for vibro-acoustic analysis of the underwater functionally graded structures in the shallow ocean
Song et al. An improved material point method using moving least square shape functions
Huang et al. A data exchange algorithm for one way fluid-structure interaction analysis and its application on high-speed train coupling interface
Pu et al. A nesting cell-based smoothed radial point interpolation method with two-level smoothed strains for static, free and forced vibration analysis of solids
Liang et al. Generation of second-mode internal solitary waves during winter in the northern South China Sea

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