CN107145626A - A kind of negative poisson's ratio structure energy-absorption box and its multidisciplinary collaboration optimization method - Google Patents
A kind of negative poisson's ratio structure energy-absorption box and its multidisciplinary collaboration optimization method Download PDFInfo
- Publication number
- CN107145626A CN107145626A CN201710201160.6A CN201710201160A CN107145626A CN 107145626 A CN107145626 A CN 107145626A CN 201710201160 A CN201710201160 A CN 201710201160A CN 107145626 A CN107145626 A CN 107145626A
- Authority
- CN
- China
- Prior art keywords
- negative poisson
- ratio
- energy
- single cell
- hypotenuse
- 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
- 238000010521 absorption reaction Methods 0.000 title claims abstract description 113
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000005457 optimization Methods 0.000 title claims abstract description 23
- 230000006835 compression Effects 0.000 claims abstract description 24
- 238000007906 compression Methods 0.000 claims abstract description 24
- 238000006073 displacement reaction Methods 0.000 claims abstract description 23
- 230000002068 genetic effect Effects 0.000 claims abstract description 7
- 230000006698 induction Effects 0.000 claims description 56
- 230000004044 response Effects 0.000 claims description 54
- 238000013461 design Methods 0.000 claims description 22
- 239000000463 material Substances 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 7
- 238000002474 experimental method Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 5
- 230000004075 alteration Effects 0.000 claims description 3
- 230000004323 axial length Effects 0.000 claims description 3
- 239000012141 concentrate Substances 0.000 claims description 3
- 230000008676 import Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 description 5
- 235000008331 Pinus X rigitaeda Nutrition 0.000 description 4
- 235000011613 Pinus brutia Nutrition 0.000 description 4
- 241000018646 Pinus brutia Species 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 238000004140 cleaning Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000001939 inductive effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Vibration Dampers (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of negative poisson's ratio structure energy-absorption box and its multidisciplinary collaboration optimization method, the negative poisson's ratio structure energy-absorption box is made up of energy-absorption box box body, front mounting plate, rear installing plate and two-dimentional negative poisson's ratio structure inner core.Because two-dimentional negative poisson's ratio structure inner core is made up of a large amount of negative poisson's ratio single cell structures, there is very big influence to energy-absorption box energy absorption in the parameter of negative poisson's ratio single cell structure, the present invention proposes multidisciplinary collaboration optimization method based on the negative poisson's ratio structure energy-absorption box, build using quality as main system, peak value impact force, compression displacement, total energy-absorbing optimize computation model for the multidisciplinary collaboration of the negative poisson's ratio structure energy-absorption box of subsystem, and archipelago genetic algorithm is used to main system, subsystem carries out multidisciplinary collaboration optimization using sequential quadratic programming algorithm.
Description
Technical field
The invention belongs to vehicle passive safety protection field, and in particular to a kind of negative poisson's ratio structure energy-absorption box and its many
Section's cooperative optimization method.
Background technology
Common energy-absorption box is hollow square structure, when a vehicle is in a collision, and it is unstable that energy-absorption box often has deformation
Fixed and insufficient the problem of, cause energy-absorption box poor to the assimilation effect of collision energy, less amount of collision energy can only be absorbed.
The energy-absorbing buffering that energy-absorption box serves is not obvious, other vehicle parts beyond remaining most of energy demand energy-absorption box
To absorb, therefore the maintenance cost of costliness is caused, very big injury more seriously can be also caused to occupant.
Negative poisson's ratio structural material shows the idio-morphosis performance different with common material when being acted on by load, makes
Obtain it and show more excellent performance in terms of energy absorption than common material, therefore negative poisson's ratio structural material can be filled
Negative poisson's ratio structure energy-absorption box is formed in traditional energy-absorption box, so as to solve common energy-absorption box well when vehicle collides
The shortcomings of deformation of presence is unstable, energy-absorbing effect is poor.
The parameter of the energy absorption and negative poisson's ratio single cell structure of negative poisson's ratio structure energy-absorption box is closely bound up, different
The energy absorption for the energy-absorption box that the single cell structure of geometric parameter is constituted also is differed, it is therefore desirable to negative poisson's ratio unit cell knot
Structure optimizes design, so as to further improve the energy-absorbing effect of energy-absorption box.
The content of the invention
It is an object of the invention to the deficiency for above-mentioned background technology, it is proposed that a kind of negative poisson's ratio structure energy-absorption box and
Its multidisciplinary collaboration optimization method.
The present invention solves its technical problem and is achieved through the following technical solutions:
A kind of negative poisson's ratio single cell structure, including two symmetrical parallel bases, the same side on two bases passes through phase
The first hypotenuse and the second hypotenuse even is attached, the equal length of the first hypotenuse and the second hypotenuse, and the first of the same side is oblique
Side and the second hypotenuse slope inwardly;The thickness of first hypotenuse, the second hypotenuse and base is t, and 0.6mm≤t≤1.2mm;
The width of first hypotenuse, the second hypotenuse and base is b, and 2.2mm≤b≤3mm;Every hypotenuse and the angle on adjacent base are
D, and 55 °≤d≤75 °;The length on two bases is a, and 12mm≤a≤16mm;Vertical range between two bases is h,
And 8mm≤h≤13mm.
Further, the thickness t of the first hypotenuse, the second hypotenuse and base is 1.1575mm;First hypotenuse, the second hypotenuse and
The width b on base is 2.2009mm;Every hypotenuse and the angle d on adjacent base are 74.684 °;The length a on two bases is
13.764mm;Vertical range h between two bases is 12.877mm.
A kind of two-dimentional negative poisson's ratio structure inner core based on above-mentioned negative poisson's ratio single cell structure, including more than one negative pool
Pine is than single cell structure, and the negative poisson's ratio single cell structure carries out array arrangement, the negative poisson's ratio unit cell along width b bearing of trends
Two ends bearing of trend of the structure along the vertical range h between two bases carries out array arrangement, and the negative poisson's ratio single cell structure
Array arrangement is carried out along the length a on two bases two ends bearing of trend;On width b bearing of trend, negative poisson's ratio unit cell
The stacking of structure symmetrically;On length a bearing of trend, the first hypotenuse of negative poisson's ratio single cell structure and adjacent negative Poisson
The second hypotenuse than single cell structure is overlapped;On the bearing of trend apart from h, the base phase mutual respect of adjacent negative poisson's ratio single cell structure
Close.
A kind of negative poisson's ratio structure energy-absorption box, it is internally provided with two-dimentional negative poisson's ratio structure inner core, this kind of negative poisson's ratio
Structure energy-absorption box include energy-absorption box box body (1), front mounting plate (2) and rear installing plate (5), described energy-absorption box box body (1) one end and
Front mounting plate (2) is connected, and energy-absorption box box body (1) other end is connected with rear installing plate (5);The front mounting plate (2) is used for
It is connected by bolt with car insurance bar crossbeam, the rear installing plate (5) is used to connect by the longeron of bolt and body of a motor car.
Further, the energy-absorption box box body (1) is that section is trapezoidal hollow prismatic structures, the energy-absorption box box
The whole surface of body (1) includes upper surface, lower surface, left surface and right flank, the trapezoidal energy-absorption box box body of the upper surface correspondence
(1) upper bottom, the bottom of the trapezoidal energy-absorption box box body (1) of the lower surface correspondence, the left surface trapezoidal suction corresponding with right flank
Two sides of energy box box body (1);
Two induction grooves one (41) are symmetrically arranged with the left surface and right flank, the induction groove one (41) is located at
At the trisection point of energy-absorption box box body (1) axial length, the induction groove one away from front mounting plate (2) is the first induction groove
One, the induction groove one away from rear installing plate (5) is the second induction groove one, and the induction groove one (41) each extends to upper surface;
The upper and lower surface is symmetrically arranged with three induction grooves two (42), before the induction groove two (42) of upper surface first is located at
The middle that installing plate (2) and the second induction groove one are projected in upper surface, upper surface Article 2 induction groove two (42) is located at first
The middle that the induction induction groove one of groove one and second is projected in upper surface, upper surface Article 3 induction groove two (42) is lured positioned at first
The middle that guide groove one and rear installing plate (5) are projected in upper surface;The induction groove two (42) and induction groove one (41) are interior
Concavity, and the depth of the induction groove two (42) is more than the depth for inducing groove one (41);From front mounting plate (2) to rear installing plate
(5) on direction, the depth of three induction grooves two (42) is gradually reduced, and the depth of two induction grooves one (41) is also gradually reduced.
A kind of multidisciplinary collaboration optimization method of negative poisson's ratio structure energy-absorption box, built in the negative poisson's ratio structure energy-absorption box
There is a kind of two-dimentional negative poisson's ratio structure inner core based on negative poisson's ratio single cell structure, the negative poisson's ratio single cell structure includes symmetrical
Two parallel bases, the same side on two bases is attached by connected the first hypotenuse and the second hypotenuse, and same
The first hypotenuse and the second hypotenuse of side slope inwardly;The thickness of first hypotenuse, the second hypotenuse and base is t, the first hypotenuse,
Second hypotenuse and the width on base are b, every hypotenuse and the angle d on adjacent base, and the length on two bases is a, two
Vertical range between base is h;The two-dimentional negative poisson's ratio structure inner core includes more than one negative poisson's ratio single cell structure, institute
State negative poisson's ratio single cell structure and carry out array arrangement along width b bearing of trends, the negative poisson's ratio single cell structure is along two bases
Between vertical range h two ends bearing of trend carry out array arrangement, and length of the negative poisson's ratio single cell structure along two bases
The two ends bearing of trend for spending a carries out array arrangement;On width b bearing of trend, negative poisson's ratio single cell structure is symmetrical
Stack;On length a bearing of trend, the of the first hypotenuse of negative poisson's ratio single cell structure and adjacent negative poisson's ratio single cell structure
Two hypotenuses are overlapped;On the bearing of trend apart from h, the base of adjacent negative poisson's ratio single cell structure overlaps, and this method includes
Following steps:
Step 1), first in ISIGHT optimization softwares, orthogonal experiment design method is chosen, in the change of each design variable
In the range of it is uniform choose N groups design sample point, the design variable is respectively the base length a, often of negative poisson's ratio single cell structure
The thickness t of vertical range h, negative poisson's ratio single cell structure between bar hypotenuse and the angle d on adjacent base, two bases and negative pool
Width b of the pine than single cell structure;
Step 2), according to the N groups design sample point of selection, set up in PROE softwares in the two-dimentional negative poisson's ratio structure of N groups
The CAD model of core;
The specific forming process of the CAD model of two-dimentional negative poisson's ratio structure inner core is:First in negative poisson's ratio single cell structure
On the basis of carry out X-direction array change to form negative poisson's ratio Multi cell structure;Then Y side is carried out to negative poisson's ratio Multi cell structure
Change to form described two-dimentional negative poisson's ratio structure inner core to the array of, Z-direction;
Step 3), the CAD model of two-dimentional negative poisson's ratio structure inner core is imported in HYPERMESH softwares, it carried out several
What cleaning and mesh generation, and the material and thickness of two-dimentional negative poisson's ratio structure inner core are set;
Step 4), by traditional energy-absorption box shell model of the inside without two-dimentional negative poisson's ratio structure inner core and for test collisions
Rigid wall model import in HYPERMESH, and two-dimentional negative poisson's ratio structure inner core is filled in traditional energy-absorption box shell, if
Put the impact velocity between rigid wall and negative poisson's ratio structure energy-absorption box, during the energy-absorption box collision of constraint negative poisson's ratio structure not with just
Property wall contact 6 frees degree of an end node, while defining contact between rigid wall and negative poisson's ratio structure energy-absorption box and defeated
Go out;
Step 5), the limit element artificial module of the N group negative poisson's ratio structure energy-absorption boxes of foundation is imported into dynamic analysis soft
Part LS-DYNA carry out solution calculating, and in HYPERGRAPH and HYPERVIEW gather energy-absorption box collision process in quality m,
The data such as peak value impact force P, compression displacement S, total energy-absorbing W;
Step 6), a kind of high-order response surface model is chosen in ISIGHT softwares, with N group negative poisson's ratios single cell structure correspondence
Base length a, the vertical range h between hypotenuse and the angle d on base, two bases, negative poisson's ratio single cell structure thickness t with
And the width b of negative poisson's ratio single cell structure is used as input, the corresponding quality m of N group negative poisson's ratio structure energy-absorption boxes, peak value impact force
P, compression displacement S, total energy-absorbing W build quality m response surface model, peak value impact force P response face mould as output parameter
Type, compression displacement S response surface model, total energy-absorbing W response surface model;
Step 7), error-tested is carried out to four response surface models fitted, i.e., calculates fitted four respectively
The coefficient R of individual response surface model2With root-mean-square error σRMSE, and with the correlation coefficient value and root-mean-square error that pre-set
Value is contrasted;The correlation coefficient value pre-set is 0.93, and the root-mean-square error value pre-set is 0.1.
Step 8), if the coefficient R of four response surface models2It is all higher than equal to the correlation coefficient value pre-set,
Square error σRMSERespectively less than it is equal to the root-mean-square error value pre-set, then performs step 9);Otherwise step 1 is re-executed)
To step 7), until the coefficient R of four response surface models2It is all higher than being equal to the correlation coefficient value pre-set, root mean square
Error σRMSERespectively less than it is equal to the root-mean-square error value pre-set;
Step 9), built in ISIGHT using quality m as main system, peak value impact force P, compression displacement S, total energy-absorbing W are
The multidisciplinary collaboration optimization computation model of the negative poisson's ratio structure energy-absorption box of subsystem;
Step 10), computation model is optimized according to the multidisciplinary collaboration of foundation, then main system calculated using archipelago heredity
Method, subsystem uses sequential quadratic programming algorithm, to the base length a, hypotenuse and the angle on base of negative poisson's ratio single cell structure
The thickness t of height h, negative poisson's ratio single cell structure between d, two bases, the width b of negative poisson's ratio single cell structure carry out multidisciplinary association
With optimization, and the Pareto disaggregation after optimize, finally concentrate one group of optimal solution of selection from Pareto solutions;Described archipelago is lost
The parameter of propagation algorithm is set to:Sub- Population Size is 10, and island quantity is 10, and genetic algebra is 10, and crossing-over rate is 1, aberration rate
With mobility equal 0.01.
Further, step 1) described in N groups design sample point in N values be 128, and step 1) described in it is orthogonal
The orthogonal arrage of experiment is L128 (45), i.e. the experiment of 5 factor, 4 level 128 times, and the excursion of each design variable is respectively:a
∈ [12,16], b ∈ [2.2,3], h ∈ [8,13], d ∈ [55 °, 75 °], t ∈ [0.6,1.2].
Further, step 4) described in rigid wall and negative poisson's ratio structure energy-absorption box between impact velocity be
16km/h。
Further, step 6) described in high-order response surface model exponent number be quadravalence, the response face mould of the quality m
Type, peak value impact force P response surface model, compression displacement S response surface model, total energy-absorbing W response surface model are respectively:
1) quality m response surface model:
M=211.3606+49.6045a+24.1501b+5.3261h+0.6283d+21.2653t-5. 3558a2-
12.2847b2-0.7646h2-0.0202d2-32.3845t2-0.0408ab-0.005471ah+0.004664ad-0.1212at-
0.004321bh-0.005127bd+0.2367bt-0.001471hd+0.01352ht-0.01615dt+0.25515a3+
2.8967b3+0.05h3+0.0002531d3+25.1231t3-0.004538a4-0.2558b4-0.001219h4-1.1354d4-
7.2777t4
2) peak value impact force P response surface model:
P=4668216.0363-524643.9041a-7708169.0354b+180530.2357h+1 16513.3631d+
676777.0333t+53084.1734a2+4440679.005b2-27221.9961h2-2741.6867d2-
1112772.4445t2-764.2867ab+92.2161ah+2.71103ad-3377.2085at-887.0286bh+0.6144bd
+8016.8117bt-13.102hd-2325.5156ht-166.4317dt-2366.5567a3-1130853.1258b3+
1820.4986h3+28.5245d3+888866.9604t3+39.4309a4+107640.8348b4-44.8933h4-0.1104d4-
257318.5951t4
3) compression displacement S response surface model:
S=3164.2166-2619.5605a-3295.4876b-1502.5493h+753.841d+10 92.6702t+
269.7967a2+1835.7762b2+219.7779h2-17.4459d2-1698.359t2+1.6865ab+0.7515ah+
0.0458ad-5.3776at-0.5554bh-0.3566bd+14.0272bt-0.06261hd-1.1016ht-0.1667dt-
12.3464a3-452.8848b3-14.1708h3+0.1794d3+1199.2214t3+0.211a4+41.6366b4+0.3391h4-
0.0006909d4-308.9682t4
4) total energy-absorbing W response surface model:
W=34302933.709+644211.1649a-50669721.7989b+2764796.5129h-
344895.9615d+7488596.9982t-77930.7763a2+30603398.0454b2-417962.829h2+
7361.5903d2-12384142.8218t2-5224.2092ab+1267.0536ah-161.2317ad+24530.8658at-
1716.4935bh+1638.7714bd-117348.4422bt+30.2318hd+2382.9861ht-4541.5054dt+
4175.07a3-8154173.7925b3+27621.0081h3-68.1648d3+9091335.913t3-83.9117a4+
809729.5918b4-677.5896h4+0.2309d4-2449116.2989t4。
Further, step 9) described in structure using quality m as main system, peak value impact force P, compression displacement S, total inhale
Energy W optimizes computation model for the multidisciplinary collaboration of the negative poisson's ratio structure energy-absorption box of subsystem:
1) quality m main systems computation model:
In formula:12≤a≤16,2.2≤b≤3,8≤h≤13,55 °≤d≤75 °, 0.6≤t≤1.2, relaxation factor ε=
0.001
2) peak value impact force P subsystem computation models:
3) compression displacement S subsystems computation model:
4) total energy-absorbing W subsystem computation models:
Beneficial effects of the present invention are:
1st, the present invention on the basis of common energy-absorption box by the two-dimentional negative poisson's ratio structure with preferable energy absorbing efficiency
Inner core is filled in wherein formation negative poisson's ratio structure energy-absorption box, so as to solve common energy-absorption box well when vehicle collides
The shortcomings of deformation of presence is unstable, energy-absorbing effect is poor.
2nd, the parameter of the energy absorption of negative poisson's ratio structure energy-absorption box and negative poisson's ratio single cell structure is closely bound up, right
Negative poisson's ratio single cell structure parameter is optimized after design, can further improve the energy-absorbing effect of energy-absorption box.
Brief description of the drawings
Fig. 1 is the structural representation of negative poisson's ratio structure energy-absorption box of the present invention;
Fig. 2 is the CAD model schematic diagram of two-dimentional negative poisson's ratio structure inner core of the invention;
Fig. 3 is the schematic diagram of negative poisson's ratio single cell structure;
Fig. 4 is the multidisciplinary collaboration optimization method schematic flow sheet of negative poisson's ratio structure energy-absorption box of the present invention.
Description of reference numerals:
The two-dimentional negative poisson's ratio structure inner core of 1- energy-absorption boxes box body, 2- front mounting plates, 3-, 41- inductions groove one, 42- induction grooves
2nd, installing plate after 5-.
Embodiment
Below by specific embodiment, the invention will be further described, and following examples are descriptive, is not limit
Qualitatively, it is impossible to which protection scope of the present invention is limited with this.
It will be understood to those skilled in the art that unless otherwise defined, all terms used herein have and the present invention
The general understanding identical meaning of those of ordinary skill in art.It should also be understood that fixed in such as general dictionary
Those terms of justice should be understood that with the meaning consistent with the meaning in the context of prior art, and unless as this
In equally define, will not be explained with idealization or excessively formal implication.
A kind of negative poisson's ratio single cell structure, as shown in figure 3, including two symmetrical parallel bases, two bases it is same
Side is attached by connected the first hypotenuse and the second hypotenuse, the equal length of the first hypotenuse and the second hypotenuse, and same
The first hypotenuse and the second hypotenuse of side slope inwardly;The thickness of first hypotenuse, the second hypotenuse and base is t, and 0.6mm≤
t≤1.2mm;The width of first hypotenuse, the second hypotenuse and base is b, and 2.2mm≤b≤3mm;Every hypotenuse and adjacent bottom
The angle on side is d, and 55 °≤d≤75 °;The length on two bases is a, and 12mm≤a≤16mm;Hanging down between two bases
Straight distance is h, and 8mm≤h≤13mm.This 5 parameters determine the complete characteristic of negative poisson's ratio single cell structure, while also determining
Its change in size.Negative poisson's ratio single cell structure profile be indent hexagonal cell structure, when bearing uniaxial compression, structure it is oblique
While occuring bending and deformation to produce Negative poisson's ratio.
Optimal solution of the present invention is:The thickness t of first hypotenuse, the second hypotenuse and base is 1.1575mm;First hypotenuse, second
Hypotenuse and the width b on base are 2.2009mm;Every hypotenuse and the angle d on adjacent base are 74.684 °;The length on two bases
A is 13.764mm;Vertical range h between two bases is 12.877mm.
A kind of two-dimentional negative poisson's ratio structure inner core based on above-mentioned negative poisson's ratio single cell structure, including more than one negative pool
Pine is than single cell structure, and negative poisson's ratio single cell structure carries out array arrangement along width b bearing of trends, and negative poisson's ratio single cell structure is along two
The two ends bearing of trend of vertical range h between individual base carries out array arrangement, and negative poisson's ratio single cell structure is along two bases
Length a two ends bearing of trend carries out array arrangement;On width b bearing of trend, negative poisson's ratio single cell structure is symmetrical
Stacking;On length a bearing of trend, the first hypotenuse of negative poisson's ratio single cell structure and adjacent negative poisson's ratio single cell structure
Second hypotenuse is overlapped;On the bearing of trend apart from h, the base of adjacent negative poisson's ratio single cell structure overlaps.Whole two dimension
Negative poisson's ratio structure inner core 3 is designed as the square endergonic structure being made up of 22*9*32=6336 negative poisson's ratio single cell structure.
A kind of negative poisson's ratio structure energy-absorption box, as shown in figure 1, it is internally provided with two-dimentional negative poisson's ratio structure inner core, should
Planting negative poisson's ratio structure energy-absorption box includes energy-absorption box box body 1, front mounting plate 2 and rear installing plate 5, and the one end of energy-absorption box box body 1 is with before
Installing plate 2 is connected, and the other end of energy-absorption box box body 1 is connected with rear installing plate 5.Front mounting plate 2 is used to pass through four bolts and automobile
Bumper bar crossbeam is connected, and rear installing plate 5 is used to connect by the longeron of four bolts and body of a motor car.Due to two-dimentional negative poisson's ratio
Structure inner core shows the idio-morphosis performance different with common material when being acted on by load so that it exists than common material
More significant performance is shown in terms of energy absorption, so that the energy absorption of energy-absorption box is lifted well.
Energy-absorption box box body 1 is that section is trapezoidal hollow prismatic structures, and the whole surface of energy-absorption box box body 1 is including upper
Surface, lower surface, left surface and right flank, the upper bottom of the trapezoidal energy-absorption box box body 1 of upper surface correspondence, the trapezoidal energy-absorbing of lower surface correspondence
The bottom of box box body 1, two sides of left surface trapezoidal energy-absorption box box body 1 corresponding with right flank.
Two induction grooves 1 are symmetrically arranged with left surface and right flank, induction groove 1 is located at energy-absorption box box body 1
At the trisection point of axial length.Induction groove one away from front mounting plate 2 is the first induction groove one, luring away from rear installing plate 5
Guide groove one is the second induction groove one, and induction groove 1 each extends to upper surface.Upper and lower surface is symmetrically arranged with three
Bar induces groove 2 42, and the induction groove 2 42 of upper surface first is located at front mounting plate 2 and second and induces what groove one was projected in upper surface
Middle, upper surface Article 2 induction groove 2 42 is located at the center that the first induction induction groove one of groove one and second is projected in upper surface
Between, Article 3 induction groove 2 42 in upper surface is located at the middle that the first induction groove one and rear installing plate 5 are projected in upper surface.Induction
Groove 2 42 is inner concavity with induction groove 1, and the depth of induction groove 2 42 is more than the depth for inducing groove 1.In the past install
Plate 2 is on the direction of rear installing plate 5, and the depth of three induction grooves 2 42 is gradually reduced, the depth of two induction grooves 1 also by
It is decrescence small, more stablize and abundant so that energy-absorption box is deformed when colliding, so as to improve the energy absorption capacity of energy-absorption box
Energy.
As shown in figure 4, a kind of multidisciplinary collaboration optimization method of negative poisson's ratio structure energy-absorption box, negative poisson's ratio structure energy-absorbing
Box is built-in with a kind of two-dimentional negative poisson's ratio structure inner core based on negative poisson's ratio single cell structure, and negative poisson's ratio single cell structure includes pair
Claim two parallel bases, the same side on two bases is attached by connected the first hypotenuse and the second hypotenuse, and together
The first hypotenuse and the second hypotenuse of side slope inwardly;The thickness of first hypotenuse, the second hypotenuse and base is t, and first is oblique
The width on side, the second hypotenuse and base is b, every hypotenuse and the angle d on adjacent base, and the length on two bases is a, two
Vertical range between individual base is h;Two-dimentional negative poisson's ratio structure inner core includes more than one negative poisson's ratio single cell structure, bears pool
Pine than single cell structure along width b bearing of trends carry out array arrangement, negative poisson's ratio single cell structure along between two bases it is vertical away from
Two ends bearing of trend from h carries out array arrangement, and length a of the negative poisson's ratio single cell structure along two bases two ends extension side
To progress array arrangement;On width b bearing of trend, the stacking of negative poisson's ratio single cell structure symmetrically;In prolonging for length a
Stretch on direction, the first hypotenuse of negative poisson's ratio single cell structure is overlapped with the second hypotenuse of adjacent negative poisson's ratio single cell structure;Away from
On bearing of trend from h, the base of adjacent negative poisson's ratio single cell structure overlaps, and this method comprises the following steps:
Step 1), first in ISIGHT optimization softwares, orthogonal experiment design method is chosen, in the change of each design variable
In the range of it is uniform choose N groups design sample point, design variable be respectively the base length a of negative poisson's ratio single cell structure, every it is oblique
The thickness t and negative poisson's ratio of vertical range h, negative poisson's ratio single cell structure between the angle d on side and adjacent base, two bases
The width b of single cell structure.
Wherein the orthogonal arrage of orthogonal test is L128 (45), i.e. 5 factor, 4 level, 128 experiments, so N values are 128, respectively
The excursion of design variable is respectively:A ∈ [12,16], b ∈ [2.2,3], h ∈ [8,13], d ∈ [55 °, 75 °], t ∈
[0.6,1.2]。
Step 2), according to the 128 of selection groups of design sample points, 128 groups of two-dimentional negative poisson's ratio knots are set up in PROE softwares
The CAD model of structure inner core.
The CAD model of two-dimentional negative poisson's ratio structure inner core is as shown in Figure 2:Initially set up negative poisson's ratio single cell structure such as Fig. 2
In Step1 shown in;The array that X-direction is carried out on the basis of negative poisson's ratio single cell structure changes to form many born of the same parents' knots of negative poisson's ratio
Structure is as shown in the Step2 in Fig. 2;Then the array change such as the Step3 in Fig. 2 of Y-direction are carried out to negative poisson's ratio Multi cell structure
It is shown;The array for finally carrying out Z-direction changes to form two-dimentional negative poisson's ratio structure inner core as shown in the Step4 in Fig. 2.
Step 3), the CAD model of two-dimentional negative poisson's ratio structure inner core is imported in HYPERMESH softwares, it carried out several
What cleaning and mesh generation, and the material and thickness of two-dimentional negative poisson's ratio structure inner core are set.
Step 4), by traditional energy-absorption box shell model of the inside without two-dimentional negative poisson's ratio structure inner core and for test collisions
Rigid wall model import in HYPERMESH, and two-dimentional negative poisson's ratio structure inner core is filled in traditional energy-absorption box shell, if
Put the impact velocity between rigid wall and negative poisson's ratio structure energy-absorption box, during the energy-absorption box collision of constraint negative poisson's ratio structure not with just
Property wall contact 6 frees degree of an end node, while defining contact between rigid wall and negative poisson's ratio structure energy-absorption box and defeated
Go out;
Impact velocity between rigid wall and negative poisson's ratio structure energy-absorption box is 16km/h, defined negative poisson's ratio structure
The output data of energy-absorption box includes:Quality m, peak value impact force P, compression displacement S, total energy-absorbing W.
Step 5), the limit element artificial module of the N group negative poisson's ratio structure energy-absorption boxes of foundation is imported into dynamic analysis soft
Part LS-DYNA carry out solution calculating, and in HYPERGRAPH and HYPERVIEW gather energy-absorption box collision process in quality m,
The data such as peak value impact force P, compression displacement S, total energy-absorbing W;
Step 6), a kind of high-order response surface model is chosen in ISIGHT softwares, with N group negative poisson's ratios single cell structure correspondence
Base length a, the vertical range h between hypotenuse and the angle d on base, two bases, negative poisson's ratio single cell structure thickness t with
And the width b of negative poisson's ratio single cell structure is used as input, the corresponding quality m of N group negative poisson's ratio structure energy-absorption boxes, peak value impact force
P, compression displacement S, total energy-absorbing W build quality m response surface model, peak value impact force P response face mould as output parameter
Type, compression displacement S response surface model, total energy-absorbing W response surface model;
The exponent number of high-order response surface model is quadravalence, quality m response surface model, peak value impact force P response surface model,
Compression displacement S response surface model, total energy-absorbing W response surface model are respectively:
1) quality m response surface model:
M=211.3606+49.6045a+24.1501b+5.3261h+0.6283d+21.2653t-5. 3558a2-
12.2847b2-0.7646h2-0.0202d2-32.3845t2-0.0408ab-0.005471ah+0.004664ad-0.1212at-
0.004321bh-0.005127bd+0.2367bt-0.001471hd+0.01352ht-0.01615dt+0.25515a3+
2.8967b3+0.05h3+0.0002531d3+25.1231t3-0.004538a4-0.2558b4-0.001219h4-1.1354d4-
7.2777t4
2) peak value impact force P response surface model:
P=4668216.0363-524643.9041a-7708169.0354b+180530.2357h+1 16513.3631d+
676777.0333t+53084.1734a2+4440679.005b2-27221.9961h2-2741.6867d2-
1112772.4445t2-764.2867ab+92.2161ah+2.71103ad-3377.2085at-887.0286bh+0.6144bd
+8016.8117bt-13.102hd-2325.5156ht-166.4317dt-2366.5567a3-1130853.1258b3+
1820.4986h3+28.5245d3+888866.9604t3+39.4309a4+107640.8348b4-44.8933h4-0.1104d4-
257318.5951t4
3) compression displacement S response surface model:
S=3164.2166-2619.5605a-3295.4876b-1502.5493h+753.841d+10 92.6702t+
269.7967a2+1835.7762b2+219.7779h2-17.4459d2-1698.359t2+1.6865ab+0.7515ah+
0.0458ad-5.3776at-0.5554bh-0.3566bd+14.0272bt-0.06261hd-1.1016ht-0.1667dt-
12.3464a3-452.8848b3-14.1708h3+0.1794d3+1199.2214t3+0.211a4+41.6366b4+0.3391h4-
0.0006909d4-308.9682t4
4) total energy-absorbing W response surface model:
W=34302933.709+644211.1649a-50669721.7989b+2764796.5129h-
344895.9615d+7488596.9982t-77930.7763a2+30603398.0454b2-417962.829h2+
7361.5903d2-12384142.8218t2-5224.2092ab+1267.0536ah-161.2317ad+24530.8658at-
1716.4935bh+1638.7714bd-117348.4422bt+30.2318hd+2382.9861ht-4541.5054dt+
4175.07a3-8154173.7925b3+27621.0081h3-68.1648d3+9091335.913t3-83.9117a4+
809729.5918b4-677.5896h4+0.2309d4-2449116.2989t4
Step 7), error-tested is carried out to four response surface models fitted, i.e., calculates fitted four respectively
The coefficient R of individual response surface model2With root-mean-square error σRMSE, and with the correlation coefficient value and root-mean-square error that pre-set
Value is contrasted;
Coefficient R2It is as follows with root-mean-square error calculation formula:
Wherein N is sample points, and P is multinomial item number, and i is the i-th sample point, fiFor the finite element fraction of i-th of sample point
Analysis value,For the response surface model calculated value of i-th of sample point,For the finite element analysis average of all sample points;
The correlation coefficient value pre-set is 0.93, and the root-mean-square error value pre-set is 0.1.
Step 8), if the coefficient R of four response surface models2It is all higher than equal to the correlation coefficient value pre-set,
Square error σRMSERespectively less than it is equal to the root-mean-square error value pre-set, then performs step 9);Otherwise step 1 is re-executed)
To step 7), until the coefficient R of four response surface models2It is all higher than being equal to the correlation coefficient value pre-set, root mean square
Error σRMSERespectively less than it is equal to the root-mean-square error value pre-set;
Step 9), built in ISIGHT using quality m as main system, peak value impact force P, compression displacement S, total energy-absorbing W are
The multidisciplinary collaboration optimization computation model of the negative poisson's ratio structure energy-absorption box of subsystem;
1) quality m main systems computation model:
In formula:12≤a≤16,2.2≤b≤3,8≤h≤13,55 °≤d≤75 °, 0.6≤t≤1.2, relaxation factor ε=
0.001
2) peak value impact force P subsystem computation models:
3) compression displacement S subsystems computation model:
4) total energy-absorbing W subsystem computation models:
Step 10), computation model is optimized according to the multidisciplinary collaboration of foundation, then main system calculated using archipelago heredity
Method, subsystem uses sequential quadratic programming algorithm, to the base length a, hypotenuse and the angle on base of negative poisson's ratio single cell structure
The thickness t of height h, negative poisson's ratio single cell structure between d, two bases, the width b of negative poisson's ratio single cell structure carry out multidisciplinary association
With optimization, and the Pareto disaggregation after optimize, finally concentrate one group of optimal solution of selection from Pareto solutions.Archipelago genetic algorithm
Parameter be set to:Sub- Population Size is 10, and island quantity is 10, and genetic algebra is 10, and crossing-over rate is 1, aberration rate and migration
Rate equal 0.01.
Described above is only the preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, under the premise without departing from the principles of the invention, some improvements and modifications can also be made, these improvements and modifications also should
It is considered as protection scope of the present invention.
Claims (10)
1. a kind of negative poisson's ratio single cell structure, it is characterised in that:The negative poisson's ratio single cell structure includes symmetrical parallel two
Base, the same side on two bases is attached by connected the first hypotenuse and the second hypotenuse, the first hypotenuse and second oblique
The equal length on side, and first hypotenuse and the second hypotenuse of the same side slope inwardly;First hypotenuse, the second hypotenuse and base
Thickness is t, and 0.6mm≤t≤1.2mm;The width of first hypotenuse, the second hypotenuse and base is b, and 2.2mm≤b≤
3mm;Every hypotenuse and the angle on adjacent base are d, and 55 °≤d≤75 °;The length on two bases is a, and 12mm≤a≤
16mm;Vertical range between two bases is h, and 8mm≤h≤13mm.
2. a kind of two-dimentional negative poisson's ratio structure inner core of the negative poisson's ratio single cell structure based on described in claim 1, its feature exists
In:Including more than one negative poisson's ratio single cell structure, the negative poisson's ratio single cell structure carries out array along width b bearing of trends
Arrangement, two ends bearing of trend of the negative poisson's ratio single cell structure along the vertical range h between two bases carries out array arrangement, and
Length a of the negative poisson's ratio single cell structure along two bases two ends bearing of trend carries out array arrangement;In width b extension
On direction, the stacking of negative poisson's ratio single cell structure symmetrically;On length a bearing of trend, negative poisson's ratio single cell structure
First hypotenuse is overlapped with the second hypotenuse of adjacent negative poisson's ratio single cell structure;On the bearing of trend apart from h, adjacent negative poisson's ratio
The base of single cell structure overlaps.
3. a kind of negative poisson's ratio structure energy-absorption box, it is internally provided with two-dimentional negative poisson's ratio structure as claimed in claim 2
Core, it is characterised in that:Including energy-absorption box box body (1), front mounting plate (2) and rear installing plate (5), the energy-absorption box box body (1) one
End is connected with front mounting plate (2), and energy-absorption box box body (1) other end is connected with rear installing plate (5);The front mounting plate (2)
For being connected by bolt with car insurance bar crossbeam, the rear installing plate (5) is used for the longeron by bolt and body of a motor car
Connection.
4. negative poisson's ratio structure energy-absorption box as claimed in claim 3, it is characterised in that:The energy-absorption box box body (1) is section
For trapezoidal hollow prismatic structures, the whole surface of the energy-absorption box box body (1) includes upper surface, lower surface, left surface
And right flank, the upper bottom of the trapezoidal energy-absorption box box body (1) of the upper surface correspondence, the trapezoidal energy-absorption box box body of the lower surface correspondence
(1) bottom, two sides of the left surface trapezoidal energy-absorption box box body (1) corresponding with right flank;
Two induction grooves one (41) are symmetrically arranged with the left surface and right flank, the induction groove one (41) is located at described
At the trisection point of energy-absorption box box body (1) axial length, the induction groove one away from front mounting plate (2) is the first induction groove one, far
Induction groove one from rear installing plate (5) is the second induction groove one, and the induction groove one (41) each extends to upper surface;On described
Surface and lower surface are symmetrically arranged with three induction grooves two (42), and the induction groove two (42) of upper surface first is located at front mounting plate
(2) and the second middle for project in upper surface of induction groove one, upper surface Article 2 induction groove two (42) is located at first and induces groove
One and second induces the middle that groove one is projected in upper surface, and upper surface Article 3 induction groove two (42) is located at the first induction groove one
The middle projected with rear installing plate (5) in upper surface;It is described to induce groove two (42) and induction groove one (41) they are inner concavity, and
The depth of the induction groove two (42) is more than the depth of induction groove one (41);From front mounting plate (2) to the direction of rear installing plate (5)
On, the depth of three induction grooves two (42) is gradually reduced, and the depth of two induction grooves one (41) is also gradually reduced.
5. a kind of negative poisson's ratio single cell structure as claimed in claim 1, it is characterised in that:First hypotenuse, the second hypotenuse and bottom
The thickness t on side is 1.1575mm;The width b of first hypotenuse, the second hypotenuse and base is 2.2009mm;Every hypotenuse and adjacent bottom
The angle d on side is 74.684 °;The length a on two bases is 13.764mm;Vertical range h between two bases is 12.877mm.
6. a kind of multidisciplinary collaboration optimization method of negative poisson's ratio structure energy-absorption box, the negative poisson's ratio structure energy-absorption box is built-in with
A kind of two-dimentional negative poisson's ratio structure inner core based on negative poisson's ratio single cell structure, the negative poisson's ratio single cell structure includes symmetrical flat
Two capable bases, the same side on two bases is attached by connected the first hypotenuse and the second hypotenuse, and the same side
The first hypotenuse and the second hypotenuse slope inwardly;The thickness of first hypotenuse, the second hypotenuse and base is t, the first hypotenuse,
Two hypotenuses and the width on base are b, every hypotenuse and the angle d on adjacent base, and the length on two bases is a, two bottoms
Vertical range between side is h;The two-dimentional negative poisson's ratio structure inner core includes more than one negative poisson's ratio single cell structure, described
Negative poisson's ratio single cell structure carries out array arrangement along width b bearing of trends, and the negative poisson's ratio single cell structure is along between two bases
Vertical range h two ends bearing of trend carry out array arrangement, and length a of the negative poisson's ratio single cell structure along two bases
Two ends bearing of trends carry out array arrangement;On width b bearing of trend, the symmetrical heap of negative poisson's ratio single cell structure
It is folded;On length a bearing of trend, the first hypotenuse of negative poisson's ratio single cell structure and the second of adjacent negative poisson's ratio single cell structure
Hypotenuse is overlapped;On the bearing of trend apart from h, the base of adjacent negative poisson's ratio single cell structure overlaps, it is characterised in that:Bag
Include following steps:
Step 1), first in ISIGHT optimization softwares, orthogonal experiment design method is chosen, in the excursion of each design variable
Interior uniform selection N groups design sample point, the design variable be respectively the base length a of negative poisson's ratio single cell structure, every tiltedly
The thickness t and negative poisson's ratio of vertical range h, negative poisson's ratio single cell structure between the angle d on side and adjacent base, two bases
The width b of single cell structure;
Step 2), according to the N groups design sample point of selection, the two-dimentional negative poisson's ratio structure inner core of N groups is set up in PROE softwares
CAD model;
The specific forming process of the CAD model of two-dimentional negative poisson's ratio structure inner core is:First in the base of negative poisson's ratio single cell structure
The array that X-direction is carried out on plinth changes to form negative poisson's ratio Multi cell structure;Then Y-direction, Z are carried out to negative poisson's ratio Multi cell structure
The array in direction changes to form described two-dimentional negative poisson's ratio structure inner core;
Step 3), the CAD model of two-dimentional negative poisson's ratio structure inner core is imported in HYPERMESH softwares, geometry is carried out to it clear
Reason and mesh generation, and the material and thickness of two-dimentional negative poisson's ratio structure inner core are set;
Step 4), by traditional energy-absorption box shell model of the inside without two-dimentional negative poisson's ratio structure inner core and for the firm of test collisions
Property wall model import in HYPERMESH, and two-dimentional negative poisson's ratio structure inner core is filled in traditional energy-absorption box shell, sets firm
Property wall and negative poisson's ratio structure energy-absorption box between impact velocity, during constraint negative poisson's ratio structure energy-absorption box collision not with rigid wall
6 frees degree of an end node are contacted, while defining the contact and output between rigid wall and negative poisson's ratio structure energy-absorption box;
Step 5), the limit element artificial module of the N group negative poisson's ratio structure energy-absorption boxes of foundation is imported into dynamics analysis software LS-
DYNA carries out solution calculating, and quality m, peak value in HYPERGRAPH and HYPERVIEW in collection energy-absorption box collision process
The data such as impact force P, compression displacement S, total energy-absorbing W;
Step 6), a kind of high-order response surface model is chosen in ISIGHT softwares, with the corresponding bottom of N group negative poisson's ratio single cell structures
The thickness t of vertical range h, negative poisson's ratio single cell structure between edge lengths a, hypotenuse and the angle on base d, two bases and negative
The width b of Poisson's ratio single cell structure is used as input, the corresponding quality m of N group negative poisson's ratio structure energy-absorption boxes, peak value impact force P, pressure
Shifting S, total energy-absorbing W condense as output parameter, quality m response surface model, peak value impact force P response surface model, pressure is built
Condense and move S response surface model, total energy-absorbing W response surface model;
Step 7), error-tested is carried out to four response surface models fitted, i.e., calculates four fitted sound respectively
Answer the coefficient R of surface model2With root-mean-square error σRMSE, and enter with the correlation coefficient value and root-mean-square error value that pre-set
Row contrast;The correlation coefficient value pre-set is 0.93, and the root-mean-square error value pre-set is 0.1;
Step 8), if the coefficient R of four response surface models2It is all higher than being equal to the correlation coefficient value pre-set, root mean square mistake
Poor σRMSERespectively less than it is equal to the root-mean-square error value pre-set, then performs step 9);Otherwise step 1 is re-executed) to step
7), until the coefficient R of four response surface models2It is all higher than being equal to the correlation coefficient value pre-set, root-mean-square error σRMSE
Respectively less than it is equal to the root-mean-square error value pre-set;
Step 9), built in ISIGHT using quality m as main system, peak value impact force P, compression displacement S, total energy-absorbing W are subsystem
The multidisciplinary collaboration optimization computation model of the negative poisson's ratio structure energy-absorption box of system;
Step 10), computation model is optimized according to the multidisciplinary collaboration of foundation, archipelago genetic algorithm, son are then used to main system
System uses sequential quadratic programming algorithm, to the base length a, hypotenuse and the angle d on base, two bottoms of negative poisson's ratio single cell structure
The thickness t of height h, negative poisson's ratio single cell structure between side, the width b progress multidisciplinary collaborations of negative poisson's ratio single cell structure are excellent
Change, and the Pareto disaggregation after optimize, finally concentrate one group of optimal solution of selection from Pareto solutions;
The parameter of described archipelago genetic algorithm is set to:Sub- Population Size is 10, and island quantity is 10, and genetic algebra is 10,
Crossing-over rate is 1, aberration rate and mobility equal 0.01.
7. a kind of multidisciplinary collaboration optimization method of negative poisson's ratio structure energy-absorption box as claimed in claim 6, it is characterised in that:
Step 1) described in N groups design sample point in N values be 128, and step 1) described in the orthogonal arrage of orthogonal test be
L128 (45), i.e. the experiment of 5 factor, 4 level 128 times, the excursion of each design variable is respectively:A ∈ [12,16], b ∈
[2.2,3], h ∈ [8,13], d ∈ [55 °, 75 °], t ∈ [0.6,1.2].
8. a kind of multidisciplinary collaboration optimization method of negative poisson's ratio structure energy-absorption box as claimed in claim 6, it is characterised in that:
Step 4) described in rigid wall and negative poisson's ratio structure energy-absorption box between impact velocity be 16km/h.
9. a kind of multidisciplinary collaboration optimization method of negative poisson's ratio structure energy-absorption box as claimed in claim 6, it is characterised in that:
Step 6) described in the exponent number of high-order response surface model be quadravalence, the response surface model of the quality m, peak value impact force P
Response surface model, compression displacement S response surface model, total energy-absorbing W response surface model are respectively:
1) quality m response surface model:
M=211.3606+49.6045a+24.1501b+5.3261h+0.6283d+21.2653t-5. 3558a2
-12.2847b2-0.7646h2-0.0202d2-32.3845t2-0.0408ab-0.005471ah
+0.004664ad-0.1212at-0.004321bh-0.005127bd+0.2367bt-0.001471hd
+0.01352ht-0.01615dt+0.25515a3+2.8967b3+0.05h3+0.0002531d3+25.1231t3
-0.004538a4-0.2558b4-0.001219h4-1.1354d4-7.2777t4
2) peak value impact force P response surface model:
P=4668216.0363-524643.9041a-7708169.0354b+180530.2357h+1 16513.3631d
+676777.0333t+53084.1734a2+4440679.005b2-27221.9961h2-2741.6867d2
-1112772.4445t2-764.2867ab+92.2161ah+2.71103ad-3377.2085at-887.0286bh
+0.6144bd+8016.8117bt-13.102hd-2325.5156ht-166.4317dt-2366.5567a3
-1130853.1258b3+1820.4986h3+28.5245d3+888866.9604t3+39.4309a4
+107640.8348b4-44.8933h4-0.1104d4-257318.5951t4
3) compression displacement S response surface model:
S=3164.2166-2619.5605a-3295.4876b-1502.5493h+753.841d+10 92.6702t
+269.7967a2+1835.7762b2+219.7779h2-17.4459d2-1698.359t2+1.6865ab
+0.7515ah+0.0458ad-5.3776at-0.5554bh-0.3566bd+14.0272bt-0.06261hd
-1.1016ht-0.1667dt-12.3464a3-452.8848b3-14.1708h3+0.1794d3+1199.2214t3
+0.211a4+41.6366b4+0.3391h4-0.0006909d4-308.9682t4
4) total energy-absorbing W response surface model:
W=34302933.709+644211.1649a-50669721.7989b+2764796.5129h -344895.9615d
+7488596.9982t-77930.7763a2+30603398.0454b2-417962.829h2+7361.5903d2
-12384142.8218t2-5224.2092ab+1267.0536ah-161.2317ad+24530.8658at
-1716.4935bh+1638.7714bd-117348.4422bt+30.2318hd+2382.9861ht-4541.5054dt
+4175.07a3-8154173.7925b3+27621.0081h3-68.1648d3+9091335.913t3-83.9117a4
+809729.5918b4-677.5896h4+0.2309d4-2449116.2989t4。
10. a kind of multidisciplinary collaboration optimization method of negative poisson's ratio structure energy-absorption box as claimed in claim 6, its feature exists
In:Step 9) described in structure using quality m as main system, peak value impact force P, compression displacement S, total energy-absorbing W is subsystems
The multidisciplinary collaboration of negative poisson's ratio structure energy-absorption box optimizes computation model:
1) quality m main systems computation model:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>M</mi>
<mi>i</mi>
<mi>n</mi>
<mi>i</mi>
<mi>m</mi>
<mi>i</mi>
<mi>z</mi>
<mi>e</mi>
<mo>:</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>Z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>,</mo>
<mi>h</mi>
<mo>,</mo>
<mi>d</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mi>P</mi>
</msub>
<mo>&le;</mo>
<mi>&epsiv;</mi>
<mo>,</mo>
<msub>
<mi>R</mi>
<mi>S</mi>
</msub>
<mo>&le;</mo>
<mi>&epsiv;</mi>
<mo>,</mo>
<msub>
<mi>R</mi>
<mi>W</mi>
</msub>
<mo>&le;</mo>
<mi>&epsiv;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
</mtable>
</mfenced>
In formula:12≤a≤16,2.2≤b≤3,8≤h≤13,55 °≤d≤75 °, 0.6≤t≤1.2, relaxation factor ε=
0.001
2) peak value impact force P subsystem computation models:
3) compression displacement S subsystems computation model:
4) total energy-absorbing W subsystem computation models:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710201160.6A CN107145626B (en) | 2017-03-30 | 2017-03-30 | Negative poisson ratio structure energy-absorbing box and multidisciplinary collaborative optimization method thereof |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710201160.6A CN107145626B (en) | 2017-03-30 | 2017-03-30 | Negative poisson ratio structure energy-absorbing box and multidisciplinary collaborative optimization method thereof |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107145626A true CN107145626A (en) | 2017-09-08 |
CN107145626B CN107145626B (en) | 2023-11-07 |
Family
ID=59783626
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710201160.6A Active CN107145626B (en) | 2017-03-30 | 2017-03-30 | Negative poisson ratio structure energy-absorbing box and multidisciplinary collaborative optimization method thereof |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107145626B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109829247A (en) * | 2019-03-04 | 2019-05-31 | 南京航空航天大学 | A kind of novel bionic bumper assembly and its multidisciplinary design optimization method |
CN111516627A (en) * | 2020-04-15 | 2020-08-11 | 吉利汽车研究院(宁波)有限公司 | Energy absorption box, energy absorption performance optimization method and automobile |
CN111707402A (en) * | 2020-06-28 | 2020-09-25 | 中国人民解放军国防科技大学 | Explosion shock wave energy passive measurement sensor based on negative Poisson ratio structure |
DE102020201772A1 (en) * | 2019-04-05 | 2020-10-08 | Ford Global Technologies, Llc | Crash box for a motor vehicle |
CN112406756A (en) * | 2020-12-01 | 2021-02-26 | 大连理工大学 | Anticollision roof beam assembly based on jump single cell structure of bullet |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103573891A (en) * | 2013-11-14 | 2014-02-12 | 马正东 | Negative Poisson ratio structural component |
CN106184085A (en) * | 2016-06-30 | 2016-12-07 | 南京航空航天大学 | A kind of closed loop negative poisson's ratio construction car security system and method for designing |
CN106202758A (en) * | 2016-07-15 | 2016-12-07 | 南京航空航天大学 | A kind of bionical negative poisson's ratio construction car guard system and method for designing |
CN207529378U (en) * | 2017-03-30 | 2018-06-22 | 南京航空航天大学 | Negative poisson's ratio single cell structure, two-dimentional negative poisson's ratio structure inner core and energy-absorption box |
-
2017
- 2017-03-30 CN CN201710201160.6A patent/CN107145626B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103573891A (en) * | 2013-11-14 | 2014-02-12 | 马正东 | Negative Poisson ratio structural component |
CN106184085A (en) * | 2016-06-30 | 2016-12-07 | 南京航空航天大学 | A kind of closed loop negative poisson's ratio construction car security system and method for designing |
CN106202758A (en) * | 2016-07-15 | 2016-12-07 | 南京航空航天大学 | A kind of bionical negative poisson's ratio construction car guard system and method for designing |
CN207529378U (en) * | 2017-03-30 | 2018-06-22 | 南京航空航天大学 | Negative poisson's ratio single cell structure, two-dimentional negative poisson's ratio structure inner core and energy-absorption box |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109829247A (en) * | 2019-03-04 | 2019-05-31 | 南京航空航天大学 | A kind of novel bionic bumper assembly and its multidisciplinary design optimization method |
DE102020201772A1 (en) * | 2019-04-05 | 2020-10-08 | Ford Global Technologies, Llc | Crash box for a motor vehicle |
CN111516627A (en) * | 2020-04-15 | 2020-08-11 | 吉利汽车研究院(宁波)有限公司 | Energy absorption box, energy absorption performance optimization method and automobile |
CN111707402A (en) * | 2020-06-28 | 2020-09-25 | 中国人民解放军国防科技大学 | Explosion shock wave energy passive measurement sensor based on negative Poisson ratio structure |
CN112406756A (en) * | 2020-12-01 | 2021-02-26 | 大连理工大学 | Anticollision roof beam assembly based on jump single cell structure of bullet |
Also Published As
Publication number | Publication date |
---|---|
CN107145626B (en) | 2023-11-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107145626A (en) | A kind of negative poisson's ratio structure energy-absorption box and its multidisciplinary collaboration optimization method | |
CN106934175B (en) | Negative poisson ratio structure energy-absorbing box and multi-objective optimization method thereof | |
Wang et al. | Structure design and multi-objective optimization of a novel NPR bumper system | |
CN208149239U (en) | Negative poisson's ratio single cell structure and three-dimensional negative poisson's ratio structure inner core and energy-absorption box | |
CN207529378U (en) | Negative poisson's ratio single cell structure, two-dimentional negative poisson's ratio structure inner core and energy-absorption box | |
CN106709151B (en) | Design method for front end structure of automobile body with 25% overlap ratio collision on front side of automobile | |
Nian et al. | Crashworthiness design of self-similar graded honeycomb-filled composite circular structures | |
US20180265023A1 (en) | Additively manufactured lattice core for energy absorbers adaptable to different impact load cases | |
CN106202758B (en) | Bionic negative Poisson ratio structure automobile protection system and design method | |
CN109670252B (en) | Force and rigidity equivalence-based head turning shrinkage model construction method and head turning shrinkage model | |
CN106184085A (en) | A kind of closed loop negative poisson's ratio construction car security system and method for designing | |
CN113642211A (en) | Composite material energy absorption box with negative Poisson's ratio structure and design method thereof | |
CN206090407U (en) | Anticollision traffic safety guardrail | |
CN103425848A (en) | Front longitudinal beam anti-collision performance design method for passenger vehicle | |
CN103303329B (en) | A kind of rail truck bottom bracket supporting endergonic structure and collision performance analogy method thereof | |
CN103569236A (en) | Method for designing front structure of body of automobile and method for designing same | |
CN115027397A (en) | Negative poisson ratio filling inner core energy absorption box based on animal horn bionic structure | |
Li et al. | Energy absorption characteristics of crash box of new honeycomb core structure with foam-filled | |
CN103770840B (en) | Method for designing front longitudinal girder | |
CN110941903A (en) | Automobile front bumper beam anti-collision performance optimization method based on DOE | |
CN107415654A (en) | Imitative ox horn vehicle door anti-collision joist and its optimization method based on negative poisson's ratio structure | |
Yang et al. | Compressive mechanical properties and dynamic behaviour of origami-inspired tri-directional auxetic metastructure | |
CN106055849A (en) | Suction and energy-storage anti-collision negative poisson's ratio structure automobile frame and design method | |
CN103279598B (en) | A kind of Variable Selection method of body of a motor car multivariate minibus optimal design | |
CN214939641U (en) | Corrugated energy absorption structure |
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 |