CN112784460A - Method for analyzing stability of mechanical metamaterial compression bar - Google Patents
Method for analyzing stability of mechanical metamaterial compression bar Download PDFInfo
- Publication number
- CN112784460A CN112784460A CN202110118958.0A CN202110118958A CN112784460A CN 112784460 A CN112784460 A CN 112784460A CN 202110118958 A CN202110118958 A CN 202110118958A CN 112784460 A CN112784460 A CN 112784460A
- Authority
- CN
- China
- Prior art keywords
- metamaterial
- compression bar
- target
- finite element
- element analysis
- 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
Images
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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Computing Systems (AREA)
- General Physics & Mathematics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention provides a method for analyzing the stability of a mechanical metamaterial compression bar, which comprises the following steps: obtaining a sample period unit of a sample metamaterial compression bar, and constructing a finite element analysis model according to the sample period unit; acquiring a target period unit of a target metamaterial compression bar, and performing finite element analysis on the target period unit according to the finite element analysis model to acquire a plurality of characteristic stiffness parameters; obtaining an integral rigidity matrix of the target metamaterial compression bar by combining the characteristic rigidity parameters according to the deformation coordination condition of the target metamaterial compression bar; establishing a buckling control equation, substituting the set boundary condition of the compression bar and the integral rigidity matrix into the buckling control equation, and solving to obtain a critical buckling load; and comparing the design load with the critical instability load, and judging the stability of the target metamaterial compression bar. The method shortens the time required by judging the structural stability of the metamaterial compression bar, and greatly improves the structural design efficiency.
Description
Technical Field
The invention relates to the technical field of metamaterial structure performance analysis, in particular to a method for analyzing the stability of a mechanical metamaterial compression bar.
Background
The metamaterial refers to a composite material which has an artificially designed structure and shows extraordinary physical properties which natural materials do not have, and the metamaterial breaks through the limitation of certain apparent natural laws so as to obtain extraordinary material functions. The metamaterial has wide application prospect and comprises the fields of electronic engineering, condensed state physics, microwave, optoelectronics, classical optics, material science, semiconductor science, nanotechnology and the like. The compression bar structure made of the metamaterial can be applied to the manufacturing industries of aircraft automobile manufacturing, engineering machinery manufacturing and the like.
In the prior art, for mechanical stability analysis of a compression bar structure constructed by using metamaterials, a 3D printing prototype experimental verification or a finite element analysis method based on a full-scale model is generally adopted. However, these two methods require a long time and have low working efficiency.
Disclosure of Invention
In view of the above, it is necessary to provide a method for analyzing the stability of a mechanical metamaterial compression bar.
A method for analyzing the stability of a mechanical metamaterial compression bar comprises the following steps: obtaining a sample period unit of a sample metamaterial compression bar, and constructing a finite element analysis model according to the sample period unit; acquiring a target period unit of a target metamaterial compression bar, and performing finite element analysis on the target period unit according to the finite element analysis model to acquire a plurality of characteristic stiffness parameters; obtaining an integral rigidity matrix of the target metamaterial compression bar by combining the characteristic rigidity parameters according to the deformation coordination condition of the target metamaterial compression bar; establishing a buckling control equation, substituting the set boundary condition of the compression bar and the integral rigidity matrix into the buckling control equation, and solving to obtain a critical buckling load; and comparing the design load with the critical instability load, and judging the stability of the target metamaterial compression bar.
In one embodiment, the target cycle unit for acquiring the target metamaterial compression bar specifically includes: the target metamaterial compression bar is formed by periodically arranging basic units, and the basic units forming the target metamaterial compression bar are used as the target periodic units.
In one embodiment, the performing finite element analysis on the target cycle unit according to the finite element analysis model to obtain a plurality of characteristic stiffness parameters specifically includes: and performing four finite element analyses on the target periodic unit, loading displacement loads in a specific form in the four finite element analyses, wherein the displacement loads comprise strain and curvature, acquiring strain energy and reaction force in a finite element analysis model according to the displacement loads, and acquiring the characteristic stiffness parameters according to the strain energy and the reaction force.
In one embodiment, the quartic finite element analysis comprises;
in the first finite element analysis, the first strain and the first curvature are set to ∈ {0,0, ∈ respectively3Acquiring first strain energy U, wherein k is {0,0,0}1And a first reaction torque according to the first strain energy and the formula U1=C3ε3 2Acquiring a first characteristic stiffness parameter C3sObtaining a second characteristic stiffness parameter H according to the ratio of the first reaction torque to the first strainsIf the metamaterial is achiral, then HsIs 0;
in the second finite element analysis, the second strain and the second curvature are set to ∈ {0,0,0} and κ {0,0, τ } respectively, and the second strain energy U is acquired2According to the second strain energy and formula U2=D3τ2Acquiring a third characteristic stiffness parameter D3s;
In the third finite element analysis, the third strain and the third curvature are respectively set to be ∈ ═ epsilon { (epsilon)10,0 and k ═ 0,0,0, and a third strain energy U is obtained3According to the third strain energy and formula U3=C1s(1-Hs 2/C3sD3s)ε1 2Acquiring a fourth characteristic stiffness parameter C1sAnd acquiring a fifth characteristic stiffness parameter C according to the symmetry2s;
In the fourth finite element analysis, the fourth strain and the fourth curvature are set to e ═ 0,0,0, and k ═ k { κ { (0, 0, 0) } respectively10,0}, and obtaining fourth strain energy U4According to the fourth strain energy and equation U4=D1s(1-Hs 2/C3sD3s)κ1 2Acquiring a sixth characteristic stiffness parameter D1sAnd obtaining a seventh characteristic stiffness parameter D according to the symmetry2s。
In one embodiment, the obtaining an overall stiffness matrix of the metamaterial compression bar according to the deformation coordination condition of the target metamaterial compression bar by combining the plurality of characteristic stiffness parameters specifically includes: if the target metamaterial compression bar is formed by repeatedly arranging m × n periodic units, and the size of each target periodic unit is a, constructing a mechanical model by combining the characteristic stiffness parameters, and obtaining integral stiffness matrixes C, B and D of the target metamaterial compression bar according to the mechanical model, wherein the integral stiffness matrixes C, B and D are as follows:
wherein, H is mnHs,C1=C2=mnC1s,C3=mnC3s,
In one embodiment, the buckling control equation is:
simultaneous equations (1) and (2), the general solution is obtained as:
in one embodiment, the boundary condition includes:
for the cantilever beam of the target metamaterial compression bar, the boundary conditions are as follows:
θ1(0)=θ2(0)=θ′1(L)-ηPθ2(L)=θ′2(L)+ηPθ1(L)=0;
for the simply supported beam of the target metamaterial compression bar, the boundary conditions are as follows:
θ1′(0)-ηPθ2(0)=θ′2(0)+ηPθ1(0)=θ′1(L)-ηPθ2(L)=θ′2(L)+ηPθ1(L)=0。
in one embodiment, the comparing the design load with the critical buckling load to determine the stability of the target metamaterial compression bar specifically includes: if the design load is larger than the critical buckling load, the target metamaterial compression bar is not stable, and buckling can occur to cause buckling; and if the design load is less than or equal to the critical buckling load, the target metamaterial is stable, and buckling cannot occur.
Compared with the prior art, the invention has the advantages and beneficial effects that:
1. according to the invention, the critical instability load of the metamaterial compression bar can be rapidly calculated or verified through a repeatable periodic unit, the design time required by judging the structural stability of the metamaterial compression bar is shortened, and the structural design efficiency is greatly improved.
2. The stiffness matrix of the metamaterial is obtained through characterization of a plurality of characteristic stiffness parameters, and the metamaterial can be suitable for chiral metamaterials or non-chiral metamaterials, so that the application range is expanded.
Drawings
FIG. 1 is a schematic flow chart of a method for analyzing the stability of a mechanical metamaterial compression bar according to an embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings by way of specific embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
In one embodiment, as shown in fig. 1, a method for analyzing the stability of a mechanical metamaterial compression bar is provided, which comprises the following steps:
and S101, obtaining a sample period unit of the sample metamaterial compression bar, and constructing a finite element analysis model according to the sample period unit.
Specifically, since the mechanical metamaterial compression bar structure is generally periodically arranged by the basic units, the basic units constituting the metamaterial compression bar can be used as the periodic units, and therefore, the periodic units can be obtained by analyzing the basic units of the metamaterial of the compression bar structure. Sample period units of a plurality of sample metamaterial compression bars can be obtained, and a finite element analysis model is constructed according to the sample period units.
And S102, acquiring a target period unit of the target metamaterial compression bar, and performing finite element analysis on the target period unit according to a finite element analysis model to acquire a plurality of characteristic stiffness parameters.
Specifically, a target period unit of the target metamaterial compression bar is obtained, and finite element analysis is carried out on the target period unit for multiple times according to the constructed finite element analysis model, so that a plurality of characteristic rigidity parameters are obtained.
In this embodiment, four times of finite element analysis needs to be performed on the target periodic unit according to the finite element analysis model to obtain seven characteristic stiffness parameters, so as to determine the stability of the target metamaterial compression bar structure according to the seven characteristic stiffness parameters. In the process of four times of finite element analysis, a specific form of displacement load is required to be loaded respectively, wherein the displacement load comprises strain and curvature, so that strain energy and reaction force are obtained, and characteristic rigidity parameters are obtained according to the strain energy and the reaction force.
And S103, obtaining an overall stiffness matrix of the target metamaterial compression bar by combining a plurality of characteristic stiffness parameters according to the deformation coordination condition of the target metamaterial compression bar.
The deformation coordination condition is a condition which can still keep the integrity and continuity of the target metamaterial compression bar after the continuous solid is deformed. The deformation coordination condition of the target metamaterial compression bar can be obtained according to the parameters when the metamaterial compression bar is designed.
Specifically, a mechanical model can be constructed according to deformation coordination conditions of a single target period unit in the overall deformation of the target metamaterial compression bar, and the mechanical model is obtained under the condition that the strain energy is equal. For example, when the strut is axially compressed to a certain strain, each periodic unit also needs to compress the corresponding strain, so the strain energy of the strut as a whole can be expressed as a multiple of the strain energy of a single target periodic unit. Wherein the bending and twisting principle of the pressure rod is similar to the compression principle of the pressure rod. Therefore, the integral rigidity matrix of the target metamaterial compression bar is obtained by combining a plurality of characteristic rigidity parameters according to the deformation coordination condition of a single target period unit.
And step S104, constructing a buckling control equation, substituting the set boundary condition of the compression bar and the integral rigidity matrix into the buckling control equation, and solving to obtain the critical instability load.
Specifically, a constructed buckling control equation is substituted according to boundary conditions of the target metamaterial compression bar and an integral rigidity matrix, the equation is solved, and the critical instability load of the target metamaterial compression bar is obtained.
And step S105, comparing the design load with the critical instability load, and judging the stability of the target metamaterial compression bar.
Specifically, since the design load is set during the design of the compression bar, the design load is compared with the critical instability load, so as to judge the stability of the target metamaterial compression bar. If the design load is larger than the critical instability load, the target metamaterial compression bar is not stable, and buckling can occur to cause instability; otherwise, the target metamaterial compression bar is stable, and buckling cannot be sent to cause instability.
In the embodiment, a sample period unit of the sample metamaterial compression bar is obtained, and a finite element analysis model is constructed according to the sample period unit; carrying out finite element analysis on a target periodic unit of the target metamaterial compression bar according to the finite element analysis model to obtain a plurality of characteristic rigidity parameters, obtaining an integral rigidity matrix of the target metamaterial compression bar according to a plurality of characteristic rigidity parameters and the deformation coordination condition of the compression bar structure, substituting the integral rigidity matrix and the boundary condition of the compression bar into a buckling control equation, solving and obtaining a critical destabilizing load, comparing the temporary destabilizing load with a design load, thereby judging the stability of the target metamaterial, representing the whole compression bar structure through repeatable periodic units, the critical destabilization load of the target metamaterial compression bar structure is calculated, so that the stability of the target metamaterial compression bar structure is judged, the time for calculating the critical destabilization load is shortened, the structural design efficiency is improved, the rigidity matrix is suitable for chiral or achiral materials, and the application range is expanded.
Wherein, step S102 specifically includes: the target metamaterial compression bar is formed by periodically arranging basic units, and the basic units forming the target metamaterial compression bar are used as the target periodic units.
Wherein, step S102 further includes: and carrying out quartic finite element analysis on the target periodic unit, loading a displacement load in a specific form in the quartic finite element analysis, wherein the displacement load comprises strain and curvature, acquiring strain energy and reaction force in a finite element analysis model according to the displacement load, and acquiring the characteristic stiffness parameter according to the strain energy and the reaction force.
Wherein the fourth finite element analysis comprises:
in the first finite element analysis, the first strain and the first curvature are set to ∈ {0,0, ∈ respectively3Acquiring first strain energy U, wherein k is {0,0,0}1And a first reaction torque according to the first strain energy and the formula U1=C3ε3 2Acquiring a first characteristic stiffness parameter C3sObtaining a second characteristic stiffness parameter H according to the ratio of the first reaction torque to the first strainsIf the metamaterial is achiral, then HsIs 0;
in the second finite element analysis, the second strain and the second curvature are set to ∈ {0,0,0} and κ {0,0, τ } respectively, and the second strain energy U is acquired2According to the second strain energy and formula U2=D3τ2Acquiring a third characteristic stiffness parameter D3s;
In the third finite element analysis, the third strain and the third curvature are respectively set to be ∈ ═ epsilon { (epsilon)10,0 and k ═ 0,0,0, and a third strain energy U is obtained3According to the third strain energy and formula U3=C1s(1-Hs 2/C3sD3s)ε1 2Acquiring a fourth characteristic stiffness parameter C1sAnd acquiring a fifth characteristic stiffness parameter C according to the symmetry2s;
In the fourth finite element analysis, the fourth strain and the fourth curvature are set to e ═ 0,0,0, and k ═ k { κ { (0, 0, 0) } respectively10,0}, and obtaining fourth strain energy U4According to the fourth strain energy and formula U4=D1s(1-Hs 2/C3sD3s)κ1 2Acquiring a sixth characteristic stiffness parameter D1sAnd obtaining a seventh characteristic stiffness parameter D according to the symmetry2s。
In particular, the fourth characteristic stiffness parameter C1sAnd a fifth characteristic stiffness parameter C2sHas symmetry therebetween, then has C1s=C2s. Similarly, the sixth characteristic stiffness parameter D1sAnd a seventh characteristic stiffness parameter D2sBetween, there is D1s=D2s。
Wherein, step S103 specifically includes: if the target metamaterial compression bar is formed by repeatedly arranging m × n periodic units, the size of each periodic unit is a, a mechanical model is constructed by combining a plurality of characteristic stiffness parameters, and the overall stiffness matrixes C, B and D of the metamaterial compression bar obtained according to the mechanical model are as follows:
Wherein, the buckling control equation in step S104 is:
simultaneous equations (1) and (2), the general solution is obtained as:
the boundary conditions in step S104 include:
for the cantilever beam of the target metamaterial compression bar, the boundary conditions are as follows:
θ1(0)=θ2(0)=θ′1(L)-ηPθ2(L)=θ′2(L)+ηPθ1(L)=0;
for the simply supported beam of the target metamaterial compression bar, the boundary conditions are as follows:
θ′1(0)-ηPθ2(0)=θ′2(0)+ηPθ1(0)=θ′1(L)-ηPθ2(L)=θ′2(L)+ηPθ1(L)=0。
wherein, step S105 specifically includes: if the design load is larger than the critical instability load, the target metamaterial compression bar is unstable, and buckling can occur to cause instability; if the design load is less than or equal to the critical buckling load, the target metamaterial compression bar is stable, and buckling cannot occur.
The foregoing is a more detailed description of the present invention that is presented in conjunction with specific embodiments, and the practice of the invention is not to be considered limited to those descriptions. For those skilled in the art to which the invention pertains, several simple deductions or substitutions can be made without departing from the spirit of the invention, and all shall be considered as belonging to the protection scope of the invention.
Claims (8)
1. A method for analyzing the stability of a mechanical metamaterial compression bar is characterized by comprising the following steps:
obtaining a sample period unit of a sample metamaterial compression bar, and constructing a finite element analysis model according to the sample period unit;
acquiring a target period unit of a target metamaterial compression bar, and performing finite element analysis on the target period unit according to the finite element analysis model to acquire a plurality of characteristic stiffness parameters;
obtaining an integral rigidity matrix of the target metamaterial compression bar by combining the characteristic rigidity parameters according to the deformation coordination condition of the target metamaterial compression bar;
establishing a buckling control equation, substituting the set boundary condition of the compression bar and the integral rigidity matrix into the buckling control equation, and solving to obtain a critical buckling load;
and comparing the design load with the critical instability load, and judging the stability of the target metamaterial compression bar.
2. The method for analyzing the stability of the mechanical metamaterial compression bar according to claim 1, wherein the obtaining of the target period unit of the target metamaterial compression bar specifically comprises:
the target metamaterial compression bar is formed by periodically arranging basic units, and the basic units forming the target metamaterial compression bar are used as the target periodic units.
3. The method for analyzing the stability of the mechanical metamaterial compression bar according to claim 1, wherein the step of performing finite element analysis on the target periodic unit according to the finite element analysis model to obtain a plurality of characteristic stiffness parameters specifically comprises:
and performing four finite element analyses on the target periodic unit, loading displacement loads in a specific form in the four finite element analyses, wherein the displacement loads comprise strain and curvature, acquiring strain energy and reaction force in a finite element analysis model according to the displacement loads, and acquiring the characteristic stiffness parameters according to the strain energy and the reaction force.
4. The method for analyzing the stability of the compression bar of the mechanical metamaterial according to claim 3, wherein the quartic finite element analysis includes;
in the first finite element analysis, the first strain and the first curvature are set to ∈ {0,0, ∈ respectively3Acquiring first strain energy U, wherein k is {0,0,0}1And a first reaction torque according to the first strain energy and the formula U1=C3ε3 2Acquiring a first characteristic stiffness parameter C3sObtaining a second characteristic stiffness parameter H according to the ratio of the first reaction torque to the first strainsIf the metamaterial is achiral, then HsIs 0;
in the second finite element analysis, the second strain and the second curvature are set to ∈ {0,0,0} and κ {0,0, τ } respectively, and the second strain energy U is acquired2According to the second strain energy and formula U2=D3τ2Acquiring a third characteristic stiffness parameter D3s;
In the third finite element analysis, the third strain and the third curvature are respectively set to be ∈ ═ epsilon { (epsilon)10,0 and k ═ 0,0,0, and a third strain energy U is obtained3According to the third strain energy and formula U3=C1s(1-Hs 2/C3sD3s)ε1 2Acquiring a fourth characteristic stiffness parameter C1sAnd acquiring a fifth characteristic stiffness parameter C according to the symmetry2s;
In a fourth finite element analysis, a fourth strain and a fourth curvature are dividedRespectively set as epsilon ═ {0,0,0} and kappa ═ { kappa { (kappa) }10,0}, and obtaining fourth strain energy U4According to the fourth strain energy and equation U4=D1s(1-Hs 2/C3sD3s)κ1 2Acquiring a sixth characteristic stiffness parameter D1sAnd obtaining a seventh characteristic stiffness parameter D according to the symmetry2s。
5. The method for analyzing the stability of the mechanical metamaterial compression bar according to claim 4, wherein the obtaining of the overall stiffness matrix of the metamaterial compression bar according to the deformation coordination condition of the target metamaterial compression bar by combining the plurality of characteristic stiffness parameters specifically comprises:
if the target metamaterial compression bar is formed by repeatedly arranging m × n periodic units, and the size of each target periodic unit is a, constructing a mechanical model by combining the characteristic stiffness parameters, and obtaining integral stiffness matrixes C, B and D of the target metamaterial compression bar according to the mechanical model, wherein the integral stiffness matrixes C, B and D are as follows:
wherein, H is mnHs,C1=C2=mnC1s,C3=mnC3s,
7. the method for analyzing the stability of the mechanical metamaterial compression bar according to claim 6, wherein the boundary conditions include:
for the cantilever beam of the target metamaterial compression bar, the boundary conditions are as follows:
θ1(0)=θ2(0)=θ′1(L)-ηPθ2(L)=θ′2(L)+ηPθ1(L)=0;
for the simply supported beam of the target metamaterial compression bar, the boundary conditions are as follows:
θ′1(0)-ηPθ2(0)=θ′2(0)+ηPθ1(0)=θ′1(L)-ηPθ2(L)=θ′2(L)+ηPθ1(L)=0。
8. the method for analyzing the stability of the mechanical metamaterial compression bar according to claim 1, wherein the comparing the design load with the critical destabilizing load to judge the stability of the target metamaterial compression bar specifically comprises:
if the design load is larger than the critical buckling load, the target metamaterial compression bar is not stable, and buckling can occur to cause buckling;
and if the design load is less than or equal to the critical buckling load, the target metamaterial is stable, and buckling cannot occur.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110118958.0A CN112784460B (en) | 2021-01-28 | 2021-01-28 | Stability analysis method for mechanical metamaterial compression bar |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110118958.0A CN112784460B (en) | 2021-01-28 | 2021-01-28 | Stability analysis method for mechanical metamaterial compression bar |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112784460A true CN112784460A (en) | 2021-05-11 |
CN112784460B CN112784460B (en) | 2023-05-16 |
Family
ID=75759386
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110118958.0A Active CN112784460B (en) | 2021-01-28 | 2021-01-28 | Stability analysis method for mechanical metamaterial compression bar |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112784460B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115252246A (en) * | 2022-09-28 | 2022-11-01 | 微创神通医疗科技(上海)有限公司 | Intracranial vascular interventional device and preparation method thereof |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105845121A (en) * | 2016-04-19 | 2016-08-10 | 黄礼范 | Sound-insulating, ventilating and heat transfer enhancing acoustic metamaterial unit, composite structure and preparation method |
CN108153981A (en) * | 2017-12-26 | 2018-06-12 | 中航沈飞民用飞机有限责任公司 | A kind of composite material fuselage Material Stiffened Panel Post-Buckling Analysis of Structures method based on finite element analysis |
CN109657284A (en) * | 2018-11-27 | 2019-04-19 | 华中科技大学 | A kind of equal geometry Topology Optimization Method towards Meta Materials |
CN110096808A (en) * | 2019-04-30 | 2019-08-06 | 中国石油工程建设有限公司北京设计分公司 | A kind of spherical shell vault Finite Element Simulation Analysis method with ribbing under multiple spot load |
CN110287637A (en) * | 2019-07-03 | 2019-09-27 | 西南交通大学 | A kind of elastic-plastic buckling bearing capacity computation method |
CN111259592A (en) * | 2020-01-20 | 2020-06-09 | 湖南工业大学 | Vibration energy collection piezoelectric metamaterial sheet material topology optimization method |
CN111597661A (en) * | 2020-06-18 | 2020-08-28 | 南昌航空大学 | Method for controlling stability of coupling processing of aluminum alloy thin-wall component |
-
2021
- 2021-01-28 CN CN202110118958.0A patent/CN112784460B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105845121A (en) * | 2016-04-19 | 2016-08-10 | 黄礼范 | Sound-insulating, ventilating and heat transfer enhancing acoustic metamaterial unit, composite structure and preparation method |
CN108153981A (en) * | 2017-12-26 | 2018-06-12 | 中航沈飞民用飞机有限责任公司 | A kind of composite material fuselage Material Stiffened Panel Post-Buckling Analysis of Structures method based on finite element analysis |
CN109657284A (en) * | 2018-11-27 | 2019-04-19 | 华中科技大学 | A kind of equal geometry Topology Optimization Method towards Meta Materials |
CN110096808A (en) * | 2019-04-30 | 2019-08-06 | 中国石油工程建设有限公司北京设计分公司 | A kind of spherical shell vault Finite Element Simulation Analysis method with ribbing under multiple spot load |
CN110287637A (en) * | 2019-07-03 | 2019-09-27 | 西南交通大学 | A kind of elastic-plastic buckling bearing capacity computation method |
CN111259592A (en) * | 2020-01-20 | 2020-06-09 | 湖南工业大学 | Vibration energy collection piezoelectric metamaterial sheet material topology optimization method |
CN111597661A (en) * | 2020-06-18 | 2020-08-28 | 南昌航空大学 | Method for controlling stability of coupling processing of aluminum alloy thin-wall component |
Non-Patent Citations (3)
Title |
---|
GAOJIAN LIN等: ""Buckling of lattice columns made from three-dimensional chiral mechanical metamaterials"", 《INTERNATIONAL JOURNAL OF MECHANICAL SCIENCES》 * |
JIAN HUA等: ""Parameters analysis and optimization of a typical multistable mechanical metamaterial"", 《EXTREME MECHANICS LETTERS》 * |
苏继龙等: ""蜂窝结构力学超材料弹性及抗冲击性能的研究进展"", 《材料工程》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115252246A (en) * | 2022-09-28 | 2022-11-01 | 微创神通医疗科技(上海)有限公司 | Intracranial vascular interventional device and preparation method thereof |
CN115252246B (en) * | 2022-09-28 | 2022-12-16 | 微创神通医疗科技(上海)有限公司 | Intracranial vascular interventional device and preparation method thereof |
Also Published As
Publication number | Publication date |
---|---|
CN112784460B (en) | 2023-05-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10909282B2 (en) | Method for rigidity enhancement and weight reduction using laser peening | |
CN106777626B (en) | A kind of trusses with discrete variables Multidisciplinary systems optimum design method | |
CN111191378A (en) | Nonlinear constitutive relation analysis method, system and device of negative Poisson ratio structure | |
CN111881531B (en) | Elastic parameter calculation and dimensionless design method for four-side concave pyramid lattice structure | |
CN107563005B (en) | A kind of tension and compression different-stiffness Instantaneous method for optimally controlling | |
CN112784460A (en) | Method for analyzing stability of mechanical metamaterial compression bar | |
Fazli et al. | Design of tensegrity structures for supporting deployable mesh antennas | |
Wang et al. | Customizable plateau in face-centered cubic hierarchical lattices achieved by self-similar embedded design | |
He et al. | Buckling analysis of thin-walled members with open-branched cross section via semi-analytical finite strip transfer matrix method | |
CN113722879A (en) | Method and system for determining mechanical property of microstructure based on homogenization theory | |
Siad et al. | Comparison of explicit and implicit finite element simulations of void growth and coalescence in porous ductile materials | |
Hart et al. | Applications of the projective-iteration versions of finite-element method in problems of damage for engineering structures | |
CN110688784A (en) | Finite element analysis method for stress performance of high-strength double-angle-steel combined section component | |
CN105350646A (en) | Two-dimensional tensegrity structure unit based on hexagon geometry | |
CN113190980B (en) | Structural topology optimization design method based on response surface method | |
CN112084683B (en) | Plastic deformation analysis method of aluminum alloy material hardware fitting for power transmission line | |
Liu et al. | Theoretical and numerical study on the in-plane mechanics of an anti-tetrachiral structure | |
Zheng et al. | Maximizing the load carrying capacity of a variable stiffness composite cylinder based on the multi-objective optimization method | |
CN115596793B (en) | Negative poisson ratio honeycomb metamaterial with multi-step deformation effect | |
El-Shaer | Topology optimization of a pressure chamber design using computer aided design software | |
CN112560164B (en) | Automobile B-pillar section design method and computer readable storage medium | |
CN117033870A (en) | Pile foundation bias reinforcement batch calculation method and computer storage medium thereof | |
Chiorean | A fast incremental-iterative procedure for ultimate strength analysis of composite cross-sections of arbitrary shape | |
Modin et al. | Deformation properties of metal woven mesh packets. Experiment and simulation (modeling) | |
Ohta | Approximating Bending Stiffness for Structural Optimization of Double-skin Hollowed Car Body Panels |
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 |