CN109858068A - Multilyer armor heat transfer model research method based on ant group algorithm - Google Patents
Multilyer armor heat transfer model research method based on ant group algorithm Download PDFInfo
- Publication number
- CN109858068A CN109858068A CN201811452775.7A CN201811452775A CN109858068A CN 109858068 A CN109858068 A CN 109858068A CN 201811452775 A CN201811452775 A CN 201811452775A CN 109858068 A CN109858068 A CN 109858068A
- Authority
- CN
- China
- Prior art keywords
- model
- thickness
- layer
- temperature
- algorithm
- 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.)
- Pending
Links
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000012546 transfer Methods 0.000 title claims abstract description 22
- 238000011160 research Methods 0.000 title claims abstract description 12
- 238000005457 optimization Methods 0.000 claims abstract description 24
- 230000001681 protective effect Effects 0.000 claims abstract description 15
- 238000009826 distribution Methods 0.000 claims abstract description 14
- 230000008569 process Effects 0.000 claims description 6
- 238000009413 insulation Methods 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 4
- 230000002159 abnormal effect Effects 0.000 claims description 3
- 230000003321 amplification Effects 0.000 claims description 3
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 3
- 230000036760 body temperature Effects 0.000 claims description 2
- 238000002834 transmittance Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims 1
- 238000013178 mathematical model Methods 0.000 abstract description 13
- 230000008859 change Effects 0.000 abstract description 10
- 210000003491 skin Anatomy 0.000 description 12
- 230000005855 radiation Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 239000004744 fabric Substances 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 239000003016 pheromone Substances 0.000 description 3
- 238000012827 research and development Methods 0.000 description 3
- 241000257303 Hymenoptera Species 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 210000002615 epidermis Anatomy 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000007773 growth pattern Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000008929 regeneration Effects 0.000 description 1
- 238000011069 regeneration method Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
The present invention relates to protective garment studying technological domains, more particularly to the multilyer armor heat transfer model research method based on ant group algorithm, include the following steps: optimization (3) thickness optimal solution of foundation (2) model of (1) model；The present invention model uses multilayer heat transfer model, and each layer is refined, keeps the situation of change of temperature spatially more specific；Model is optimized using traditional linear optimization, dichotomy exhaustion and ant group algorithm, so that thickness estimated value is more accurate, close to actual conditions, traditional algorithm and modern optimization algorithm are compared simultaneously, more conducively solve corresponding practical problem, calculate the optimal thickness of Temperature Distribution and each layer of protective garment, entire invention clear thinking, the temperature variations on the outside of dummy's skin are determined with mathematical model, it reduces costs, shortens the R&D cycle, reduce certain error range, the time solved the problems, such as is reduced, practicability is high.
Description
Technical field
The present invention relates to protective garment studying technological domains, more particularly to the transmitting of the multilyer armor heat based on ant group algorithm
Model research method.
Background technique
Ant group algorithm is also known as ant algorithm, is a kind of probability type algorithm for finding path optimizing in figure.Its inspiration
The behavior in path is found during search of food from ant.Ant group algorithm is a kind of simulated evolutionary algorithm, and preliminary grinds
Study carefully and shows that the algorithm has many excellent properties.For PID controller parameter optimization design problem, by ant colony algorithm for optimization design
Result and the result of genetic Algorithm Design compare, Numerical Simulation Results show that ant group algorithm has a kind of new mould
The validity and application value of quasi- evolution optimization method.
Human body thermal protection is the important research direction of public safety field.Rationally human body-clothes under assessment hot environment
Dress-environment heat exchange, the protective garment for designing optimal thickness is the important means for ensureing their life securities.
Existing protective garment is generally treble cloths material composition, and every layer of thickness is not designed particularly, so that
The heat insulation of protective garment is undesirable, incomplete to the protection of human body；In order to preferably design this spe-cial-purpose uniform, we
Dummy is placed in the hot environment in laboratory, the temperature on the outside of dummy's skin is measured.To simulate true man in actual environment
In situation variation, the temperature variations on the outside of dummy's skin are determined using mathematical model, calculate Temperature Distribution and anti-
Shield takes the optimal thickness of each layer.
Summary of the invention
The purpose of the present invention is to provide the multilyer armor heat transfer model research methods based on ant group algorithm, with solution
Certainly above-mentioned technical problem.
The present invention using following technical scheme in order to solve the above technical problems, realized:
Include the following steps:
(1) foundation of model
A. imported data to first by step-length group forming criterion, according to related physical calorifics principle and protective garment layering into
Row derives, it is established that each layer of model I；
B. according to the relationship of each layer contact area, model I is improved, the model II after being optimized；
C. protective garment is divided into I, II, III layers by ecto-entad, there is also certain between layer III and human skin
Gap, using this gap as Section IV layer；
D. it obtains temperature and specifically divides obtained Temperature Distribution according to the number of plies at any time；
E. partial differential equation are adjusted using calculus of finite differences；
F. Matlab software is utilized, to 54 groups of data have been extracted in attachment, every group of 10 collecting samples carry out data
Screening removes abnormal data and is supplemented to obtain relevant table answer；
(2) optimization of model
A. estimated value is obtained by exhaustion by the way of traditional linear optimization and dichotomy；
B. using the ant group algorithm in self start type algorithm, the concrete condition of the model is iterated, is obtained more accurate
Thickness estimated value；
C. most suitable II thickness degree is determined；
(3) Optimum Solution
A. from single thickness, amplification to two thickness, respectively II layer with IV layers；
B. estimated value is obtained by exhaustion by the way of traditional linear optimization and dichotomy；
C. the iteration that the concrete condition of the model carries out is obtained more smart using the ant group algorithm in self start type algorithm
True thickness estimated value；
D. the optimal thickness of II layers He IV layers is obtained；
Preferably, the multilyer armor mainly includes shell, waterproof layer and thermal insulation layer.
Preferably, in the step (3) for environment temperature be 80 DEG C when, it is to be ensured that work 30 minutes when, dummy's skin
Temperature outside is no more than 47 DEG C, and the time more than 44 DEG C is no more than 5 minutes.
Preferably, during the step (1), body temperature control is placed on to the height in laboratory in 37 DEG C of dummy
In warm environment, the temperature on the outside of dummy's skin is measured, variation the case where to simulate true man in the actual environment.
Preferably, the material of the multilyer armor is uniform, and structure is constant in transmittance process, and buckling factor is normal
Number.
The beneficial effects of the present invention are:
The present invention model uses multilayer heat transfer model, and each layer is refined, makes the situation of change of temperature in sky
Between it is upper more specific；Model is optimized using traditional linear optimization, dichotomy exhaustion and ant group algorithm, so that
Thickness estimated value is more accurate, compares close to actual conditions, while to traditional algorithm and modern optimization algorithm, more conducively solves
Certainly corresponding practical problem, calculates the optimal thickness of Temperature Distribution and each layer of protective garment, and entire invention clear thinking uses number
Model is learned to determine the temperature variations on the outside of dummy's skin, reduces costs, shortens the R&D cycle, reduce certain
Error range, practicability are high.
Detailed description of the invention
Fig. 1 is simulation system schematic diagram of the invention；
Fig. 2 is time and dummy's shell temperature schematic diagram of the invention；
Fig. 3 is that Multi-layer thermal of the invention transmits mathematical model analogue data schematic diagram；
Fig. 4 is ant group algorithm analog case schematic diagram of the invention；
Specific embodiment
In order to be easy to understand the technical means, the creative features, the aims and the efficiencies achieved by the present invention, below
In conjunction with specific embodiments, the present invention is further explained, and however, the following embodiments are merely preferred embodiments of the present invention, not entirely
Portion.Based on the implementation example in the implementation mode, those skilled in the art obtain it without making creative work
Its embodiment, belongs to protection scope of the present invention.Experimental method in following embodiments is unless otherwise specified routine
Method, the materials, reagents and the like used in the following examples are commercially available unless otherwise specified.
By the information given, show that professional heat-protective clothing-air layer-skin simulation system as shown in Figure 1 is shown
It is intended to, it is thus understood that the profession heat-protective clothing mainly includes shell, waterproof layer and thermal insulation layer, and the air between fabric and skin
Layer has certain thickness.
Case study
In order to calculate the distribution of temperature, and generate corresponding table.We are understood to the temperature according to each layer
Distribution situation is spent, each layer is again careful to divide the number of plies according to certain standard.We are from table it can be found that in preceding 15s
It is interior, dummy's skin temperature outside be always maintained at it is constant, then with clothes absorb heat increase, temperature is begun to ramp up, but is worked as
Temperature reaches 48.08 DEG C, have passed through at this time for a long time, ambient temperature remains at 75 DEG C, and the shell temperature of dummy is simultaneously
There is no variations, and according to these data, by above analysis, we have tried progress to take table every the interval of 10s
In corresponding time and skin temperature, to obtain corresponding 540 groups of data, following table provides the parameter situation of some fixations.
| The number of plies | I | II | III | IV |
| Thickness (mm) | 0.6 | 0.6-25 | 3.6 | 0.6-6.4 |
The specific thickness and ambit table of 1 number of plies of table
The related given parameters of 2 number of plies of table
Multilyer armor heat transfer model research method based on ant group algorithm, characterized by the following steps:
(1) foundation of model:
A. imported data to first by step-length group forming criterion, according to related physical calorifics principle and protective garment layering into
Row derives, it is established that each layer of model I；
B. according to the relationship of each layer contact area, model I is improved, the model II after being optimized；
C. protective garment is divided into I, II, III layers by ecto-entad, there is also certain between layer III and human skin
Gap, using this gap as Section IV layer；
D. it obtains temperature and specifically divides obtained Temperature Distribution according to the number of plies at any time；
E. partial differential equation are adjusted using calculus of finite differences；
F. Matlab software is utilized, to 54 groups of data have been extracted in attachment, every group of 10 collecting samples carry out data
Screening removes abnormal data and is supplemented to obtain relevant table answer；
The understanding of logical heat transfer form, it is known that we are primary concern is that heat transfer and the correlation of heat radiation
Knowledge, by heat transfer law in calorifics, we carry out the model foundation of the problem.
Fourier heat equation is mainly explained as follows, in Heat Conduction Phenomenon, by uniform section in the unit time
Heat, direct proportion in perpendicular to the rate of temperature change and area of section on the interface direction, and the direction of heat transfer then with temperature
It spends raised contrary.According to correlation it is assumed that it is contemplated that at any time perpendicular to the temperature of dummy surface and every layer of fabric
Situation of change.Using Fourier law expression formula, the heat transfer formula of our available stationary states, formula is as follows:
Wherein, Q indicates the heat of conduction；K is the pyroconductivity of substance, the i.e. coefficient of heat conduction, with substance nature,
The temperature of substance is related.Due to we assuming that in have been described above influence of the temperature to fabric itself under experimental situation little,
So certain situation of change of k has been ignored, S refers to the cross-sectional area of object；It is temperature gradient, indicates that temperature becomes
The intensity of change, it is equal to temperature T dividing partially to coordinate x.
We also indicate the intensity of heat transfer using the intensity of heat flow simultaneously:
Spread speed due to giving heat in Fourier law is infinitely great, if that is, Q, S, k are it is known that distribution temperature
Degree T (x) is function about x, unrelated with time t, does not meet situation described in the century of the topic, so we need pair
The model carries out certain amendment, also obtains relevant model.
The heat that clothes absorbs, the heat that dummy absorbs, the heat that dummy is distributed, the heat that clothes is distributed, between
Energy consumption and energy acquisition, meet conservation of energy formula, thus we are it can be concluded that conservation of energy equation, from
And providing relevant equation for model foundation is proved.
QIt absorbs=QIt distributes
X-Y scheme as shown in Figure 2 transmits mathematical model by Multi-layer thermal and carries out the specific Temperature Distribution feelings of hierarchical solving
Condition.
The basic theory of two above is simultaneously subject to corresponding chart, we obtain an approximate mathematical model, finally
Certain optimization is being carried out to the model, to obtain following course of solving questions.
In view of sample data 5400, data volume is huge, and from table come carry out it is certain it has been observed that
After the variation of certain time interval, the amplitude of variation very little of data, so we are using 10s as interval, by 5400 groups
Data take out 540 groups every time, the case where by this more than 500 a data, carry out relevant solution.
In order to calculate the distribution of temperature, and generate corresponding table.The problem can be interpreted as each layer by we
Profiling temperatures.Then we are again more careful to each layer divides the number of plies according to certain standard.
Using Fourier heat equation, we separately consider different layers of situations corresponding to different-thickness, obtain
To its every layer of PDE model.
We are expressed as γ specific each layer of dividing condition, the case where by from known dummy's skin surface temperature,
Start boundary's environment derivation from inside to outside with Section IV layer, according to table one it is found that I layers of range γ1For (L4+L3+L2, L4+L3+L2+
L1) i.e. [0.146m, 0.152m], II layers of range be γ2(L4+L3, L4+L3+L2) i.e. [0.086m, 0.146], III layers of model
Enclose γ3For (L4, L4+L3) i.e. [0.05m, 0.086m], the range of Section IV layer is γ4(0, L4) i.e. [0,0.05].
Firstly, tier I is due to being only merely to contact with external environment, we are according to certain equation of heat conduction and suction
It receives heat and releases heat and meet law of conservation of energy relationship, it is contemplated that the external world is in the environment and heat radiation of a constant temperature
The case where, list the mathematical model of tier I.
(x, t)=γ1×(0,texp) (1)
Wherein T indicates the temperature (DEG C) solved required for us, and t is the time, and x indicates corresponding number of plies area where it
Domain is embodied as pyroconductivity corresponding to first layer, indicates that the sensible heat of corresponding tier I holds, sensible heat appearance passes through related data
It can be seen that its calculating process is as follows, it is desirable that correspond to the density of the number of plies, and the specific heat capacity of the corresponding number of plies is
Can, specific equation is as follows:
And the F on the right of (1) formula equationlAnd FRIt is expressed as the difference of External Heat Flux, i.e., left radiation-right radiation calculating side
Formula, specific expansion are as follows:
Wherein β is radiation absorption constant, can be learnt substantially by pertinent literature 40, K-1, σ be Stefan-bohr hereby
Graceful Changshu, 5.670*10-8w/(m2*k4)
It can be seen that I layers of primary condition meets:
T (x, 0)=T1(x),x∈(0,Lfab)
Wherein LfabMeet L1+L2+L3, and correspond to the tier I right boundary condition we it can be concluded that the following formula into
Row boundary value solves
(1-ξ1)FL(0,t)+ξ1σT4(0, t)=FR(0,t)
In the above-mentioned formula, ξ1It is expressed as a radiance of I layers of left margin, unit is w/ (s*m2),qairIt is expressed as
The heat radiation heat flow density of external environment, unit w/ (m2),kairIt is expressed as the pyroconductivity of air in external environment.
Tier ii, it is contemplated that the case where fabric on its both sides is all solid-state, so we are according to the equation of heat conduction and energy
Law of conservation is measured, the mathematical model of tier ii is listed.
WhereinThe sensible heat for being still expressed as corresponding tier ii holds, and x is still the thickness area of the corresponding number of plies
γ2, the t expression corresponding time, k2It is expressed as the pyroconductivity of corresponding tier ii.
Same settling mode still meets following equation for the contact surface of tier ii and layer III
(1-ξ2)FL(L1,t)+ξ2σT4(L1, t) and=FR(L1,t)
Layer III, it is contemplated that the case where third layer is solid-state on one side, and another side is air layer, we draw both sides
Boundary condition is different, and boundary condition is separated and lists different Soil-structure Interaction Models, while centre is still according to heat transfer and energy
Conservation lists the model of layer III.
WhereinThe sensible heat for being expressed as corresponding layer III holds, and x is still the thickness area γ of the corresponding number of plies3, t
Indicate corresponding time, k3Be expressed as the pyroconductivity of corresponding layer III, required method be still also temperature for when
Between partial differential equation.
According to III layers and IV layers of contact surface, still meet the following formula:
(qconv+qrad)|X=0=hc,fl(Tg-T|X=0)
In formula: qconvFor the thermal convection density of external environment to I, and qradFor external environment to I layers of heat radiation density；
he,flThermal convection density between external environment and the outer surface I, TgFor ambient temperature.
Section IV layer, it is contemplated that Section IV layer, both bounded sides are fabric and dummy's epidermis respectively, while inside is also
Gas, we carry out certain processing to the model of solid before, obtain the mathematical model for belonging to gas Section IV layer, while right
Both sides contact surface lists corresponding boundary condition model respectively.
Wherein, γ4=(LfabrLfab+Lair), and Lfab=L1+L2+L3, while LairIV layers of corresponding thickness are expressed as, by
Smaller in the thickness of air layer, we can be by being regarded as a regular different shape objects, thus more
The relative theory of limited heat transfer and thermal convection using in this space, before and is solved, comprehensively considers air
Pass to and to professional protective garment tool have a certain impact, so that it should be one that we, which can release this IV layers of air,
It is the phenomenon that stable state, as follows to the primary condition of the air layer:
T (x, 0)=T1(x),x∈γ4
We carry out integration arrangement to formula (1,2,3,4), while the correlative encountered in the above problem being brought into,
Finally obtain our Multi-layer thermal transmitting mathematical model.
It is therein mainly to need condition as follows:
T (x, 0)=T1(x),x∈(0,Lfab)
T (x, 0)=T1(x),x∈γ4
Using the separation of variable, we can solve the specific discrete formula and phase of above four layers of different distributions model
The boundary condition and primary condition of pass, so that generating corresponding program optimizes calculating.But analogue data at the beginning is anti-
And the temperature value amplitude of variation that 0.1mm is found out before ambient temperature is very big, and former topic is by error, therefore to above situation into
The certain model refinement of row.
Current model refinement, we are more than and divide according to thickness, we can be with from data given in table
It was found that do not influenced on dummy's skin in 15s seconds first, temperature be always maintained at it is constant, then as clothes absorbs the increasing of temperature
Add, temperature is begun to ramp up, but when 48.08 DEG C of temperature arrival, have passed through at this time for a long time, ambient temperature remains at 75
DEG C, and the shell temperature of dummy does not change.
According to these data, by above analysis, we utilize Fourier heat equation, are reacted according to fig. 2
Function growth pattern, it has been found that the situation of change of temperature value is relatively slow, so, we are by thickness Li(i=
1,2,3,4) and thickness range, γi(i=1,2,3,4) divides value according to the range of 0.1mm, i.e., by overall thickness Lsum=
152mm is divided into 152 parts, and 153 number values are discussed according to the situation of every portion since 0 part, therefore in total, simultaneously
The case where I layers and environment contact surface is added column model, successively II layers and III layers of consideration of contact surface establishes model later, then examines
The contact surface considered between IV layers and III layers accounts for, and by the model of the above contact surface, is added to our well-established moulds
Type carries out related optimization, and last model result is as follows:
The completely new Multi-layer thermal transmitting mathematical model successively established using us, by calculus of finite differences and the separation of variable,
Partial differential equation are carried out to certain adjustment, obtain corresponding discrete equation, then according to thickness Li(i=1,2,3,4) locating
Different range γi(i=1,2,3,4) brings mathematical model corresponding to the above-mentioned different number of plies into, finally by the time according to taking
10s is a boundary, obtains corresponding every layer of temperature value, it is contemplated that the data volume of temperature value is excessive, we have chosen every layer
Situation of change of the middle dividing value in corresponding time 50min, while obtaining corresponding mathematical image, such as Fig. 3.
The substantially numerical value that we will be calculated, generates corresponding form attachment, and table 3 only does part displaying, as asks
Solution corresponding to topic one, Temperature Distribution therein are refined into the calculating under related different levels again.
Table 3 changes over time local table based on the portion temperature that model solution comes out
(2) optimization of model
By regular hour exhaustion by the way of traditional linear optimization and dichotomy, certain estimated value is obtained,
We are carried out the iteration of certain number to the concrete condition of the model, are obtained using the ant group algorithm in self start type algorithm later
More accurate thickness estimated value, to reach the estimated value for being more nearly true actual conditions.
It needs to be determined that optimal tier ii depth information, we are by reference to relevant parameter provided by table 2, Ke Yichu
Step recognizes the basic demand of tier ii thickness itself, so that we obtain first relevant constraint condition (1), constrains item
Part is as follows:
L2∈[0.6mm,25mm]
According to system of unit scaling results, it is known that thickness L itself2Meet in 0.6*10-3To 25*10-3Between m, with this
As first constraint condition.
According to the relevant information that second and third constraint condition and topic provide, our available following formulas
Son:
TDummy surface≤47℃
The constraint condition that can determine that now according to us, and start given certain initial value, we release three kinds
The different specific optimisation strategies for method for solving.
The method that we think deeply first is directly according to exhaustive method, by L2All situations enumerated one by one, it
First successively is brought into as known conditions afterwards and asks that released Multi-layer thermal transmits mathematical model, to obtain corresponding temperature
And the time, according to known two or three constraint condition, to substantially obtain an estimated value of Section II thickness degree.
The process of the method is comparatively fairly simple, and still, specific each value requires to solve experiment, i.e., using steady
The mode that fixed step-length is 0.1mm empirical value in given section, the process of answer is still considerably complicated, Er Qieji
It spends the time.
So far we release second set of calculation, for first constraint condition, it is proposed that dichotomy is asked
Solution, i.e., intercept median, i.e. median in given II thickness section every time, is conducted into first and asks to solve
Model, carry out verifying specific TDummy surfaceActual temp distribution situation, thus by the time complexity of first set calculation method
The specific time of half, but the data volume being to solve for are optimized, is still sizable, it is possible thereby to illustrate, most of tradition
Linear programming equation optimisation strategy, as can be seen that all data are enumerated from a large amount of data plane, also or
The half situation of original data is solved, is all to spend time and very big period.
Although dichotomy it is specific solution have certain complexity, after all still obtained constraints above condition it
Under obtained best L2, optimum thickness that we obtain in 9.2mm or so but subtle optimization we have no idea more into one
Step reduces error, therefore we carry out related summary to both the above method.
Both the above resolution policy, all requiring to reduce research and development cost to problem itself and shortening the R&D cycle is not phase
Symbol.To this, we are proposed the calculation method of the third mainstream, iteration efficiency repeatedly can be improved, while simplifying and asking
The specific consideration situation of topic.
It is proposed that the constraint condition of the linear programming is optimized and arranged using ant group algorithm, main thought
It is as follows:
We utilize the simple rule of conduct of ant itself, guarantee that endless loop will not be entered when Food Recruiment In Ants, thus
The L that will be obtained desired by us2More rapid and convenient, so as to significantly avoid many duplicate calculating from either leading
The case where some numerical value different from common sense occurred during entering model.It is both the exploration thinking for simulating ant, Yi Xieyu
The bias for actually wanting to solution value is excessive, i.e., the excessive value of error we can use the algorithm and effectively avoid, thus by former
There is the calculation amount of data half to further reduce, improve whole performance, as more convenient optimization algorithm, we are detected
Some samples generated at random, the sample are generated by correlation random number, so that it is original to find that the speed of its operation is substantially better than
Dichotomy or exhaustive thinking, and then using the algorithm to this problem solving.
We obtain more accurate value, specific mathematical thought is as follows by the algorithm:
By probing into, above-mentioned experiment has reacted ant in one of group behavior information positive feedback phenomenon.Ant
Food is searched for by this Information exchange mechanism between body.And the chemical factor for being used to AC regeneration is referred to as by us now
For --- " pheromones ".
Then relevant mathematical model is established, firstly, the sum and in the past one of the pheromones at the both ends of professional thermal protective clothing
It is directly proportional (i.e. every ant pheromone release ability all having the same) by the ant number of the bridge in the section time；Secondly,
Assuming that certain moment ant selected according to the number of residual risk amount on every layer model wherein most suitable II thickness degree formed
Path, more by the ant number in path corresponding to this thickness, then the residual risk element total amount on the II thickness is got over
Greatly.Then after all M ants pass through obstacle, the probability of (M+1) ant selection path A (i.e. optimum thickness) are as follows:
PB(m)=1-PA(m)
Wherein parameter h (expecting factor) and k (heuristic factor) simulates matched data, and (M+1) ant is according to above-mentioned
Formula calculates probability, generates equally distributed random number w on [0,1] later, selects the path A if w≤P (A), otherwise select B
Path, finally draws associated picture according to above AB value, and iteration obtains an approximation.The result that we calculate should II
Optimal thickness is 9.28mm.
Iteration diagram therein is as shown in figure 4, and relevant error analysis figure
Using we at the beginning dichotomy find out come proximity values, can simulate that we want repeatedly is most suitable for
II layer thickness, obtained result pass through successive ignition, due to be basically stable at a definite value do not occur it is subsequent
Iteration variation, suspend current the number of iterations, the content obtained is as shown in Figure 4, it can be seen that carry out Current Function
Value be exactly that can finally stablize in 9.28mm or so in fluctuation up and down repeatedly.
(3) Optimum Solution
A. from single thickness, amplification to two thickness, respectively II layer with IV layers；
B. estimated value is obtained by exhaustion by the way of traditional linear optimization and dichotomy；
C. the iteration that the concrete condition of the model carries out is obtained more smart using the ant group algorithm in self start type algorithm
True thickness estimated value；
D. the optimal thickness of II layers He IV layers is obtained；
According to the available the following formula of general contents of step (2), asks and still set up at this, and can be used as about
Beam condition considers that main constraints are as follows:
L2∈[0.6mm,25mm]
TDummy surface≤47℃
We can be clearly understood that the optimal of the more than tier ii that needs solve, therefore according in identical table 2
Master data, our available following new restraint conditions:
L4∈[0.6mm,6.4mm]
After the completion of above four kinds of situation simultaneous, then the multi-level heat transfer model established with step at the beginning (1) into
The certain optimization of row, while we can be obtained again according to the basic ideas of step (2) with fundamental inference: exhaustion etc. is traditional
The case where calculating tuning of Optimized model and equation in this topic is more difficult, and conventional method is more dependent on initial solution is evolved
Algorithm has more opposite global search ability, and practicability is wide, it would be desirable to which more convenient and function is more powerful
Alternative manner solves this problem.
In view of current problem is more than simple dimension word variation, mainly by L4And L2The two dimension of the two composition
Matrix is changed, and simple inference simulation can not be made.So far it is proposed that a kind of thought algorithm of traversal and optimization, considers
Whether thinking processes optimization can be continued in the case where conventional method can not play certain effect.
This method mainly utilizes ant group algorithm to simulate the Main change under binary restraint condition, main algorithm thought and step
Suddenly the case where (2), is consistent, is only merely the complexity that the complexity of unitary variant is extended to two variables, by above
Related thinking, we have to carry out according to the multilayer heat transfer model come given by step (1) the temperature value in certain
Analogue data is carried out, by our obtained related data deposit correlation graphs, is optimized followed by Matlab software
Analog image, image are deduced according to current intelligence, finally obtain the data of some relatively correlations.
Thus the result that we deduce out substantially, II layers of value, that is, L2=18.3mm, and IV layers of value, that is, L4=
3.9mm。
The basic principles, main features and advantages of the present invention have been shown and described above.The technology of the industry
For personnel it should be appreciated that the present invention is not limited to the above embodiments, described in the above embodiment and specification is only the present invention
Preference, be not intended to limit the invention, without departing from the spirit and scope of the present invention, the present invention also has respectively
Kind changes and improvements, these changes and improvements all fall within the protetion scope of the claimed invention.The claimed scope of the invention by
Appended claims and its equivalent thereof.
Claims (5)
1. the multilyer armor heat transfer model research method based on ant group algorithm, characterized by the following steps:
(1) foundation of model
A. it is imported data to first by step-length group forming criterion, the layering according to related physical calorifics principle and protective garment is pushed away
It leads, it is established that each layer of model I；
B. according to the relationship of each layer contact area, model I is improved, the model II after being optimized；
C. protective garment is divided into I, II, III layers by ecto-entad, there is also gaps between layer III and human skin, by this sky
Gap is as Section IV layer；
D. it obtains temperature and specifically divides obtained Temperature Distribution according to the number of plies at any time；
E. partial differential equation are adjusted using calculus of finite differences；
F. to 54 groups of data are extracted in attachment, every group of 10 collecting samples screen data, remove abnormal data and are mended
It fills to obtain relevant table answer；
(2) optimization of model
A. estimated value is obtained by exhaustion by the way of traditional linear optimization and dichotomy；
B. using the ant group algorithm in self start type algorithm, the concrete condition of the model is iterated, more accurate thickness is obtained
Spend estimated value；
C. most suitable II thickness degree is determined；
(3) thickness optimal solution
A. from single thickness, amplification to two thickness, respectively II layer with IV layers；
B. estimated value is obtained by exhaustion by the way of traditional linear optimization and dichotomy；
C. the iteration that the concrete condition of the model carries out is obtained more accurate using the ant group algorithm in self start type algorithm
Thickness estimated value；
D. the optimal thickness of II layers He IV layers is obtained.
2. the multilyer armor heat transfer model research method according to claim 1 based on ant group algorithm, feature exist
In: the multilyer armor mainly includes shell, waterproof layer and thermal insulation layer.
3. the multilyer armor heat transfer model research method according to claim 1 based on ant group algorithm, feature exist
In: in the step (3) for environment temperature be 80 DEG C when, it is to be ensured that work 30 minutes when, dummy's skin temperature outside does not surpass
47 DEG C are crossed, and the time more than 44 DEG C is no more than 5 minutes.
4. the multilyer armor heat transfer model research method according to claim 1 based on ant group algorithm, feature exist
In: during the step (1), by body temperature control in the hot environment that 37 DEG C of dummy is placed on laboratory, measurement
Temperature on the outside of dummy's skin, variation the case where to simulate true man in the actual environment.
5. the multilyer armor heat transfer model research method according to claim 1 based on ant group algorithm, feature exist
In: the material of the multilyer armor is uniformly that structure is constant in transmittance process, and buckling factor is constant.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811452775.7A CN109858068A (en) | 2018-11-30 | 2018-11-30 | Multilyer armor heat transfer model research method based on ant group algorithm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811452775.7A CN109858068A (en) | 2018-11-30 | 2018-11-30 | Multilyer armor heat transfer model research method based on ant group algorithm |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109858068A true CN109858068A (en) | 2019-06-07 |
Family
ID=66890518
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811452775.7A Pending CN109858068A (en) | 2018-11-30 | 2018-11-30 | Multilyer armor heat transfer model research method based on ant group algorithm |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109858068A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111783281A (en) * | 2020-06-12 | 2020-10-16 | 永基众合(山东)智能科技有限公司 | Method for optimizing thickness of multilayer heat-insulating medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN2833585Y (en) * | 2005-07-19 | 2006-11-01 | 东华大学 | Tester for heat conduction performance of protection fabric |
| JP2014091890A (en) * | 2012-11-06 | 2014-05-19 | Teijin Ltd | Cloth for thermoconductive activity protection garment and thermoconductive activity protection garment using the same |
-
2018
- 2018-11-30 CN CN201811452775.7A patent/CN109858068A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN2833585Y (en) * | 2005-07-19 | 2006-11-01 | 东华大学 | Tester for heat conduction performance of protection fabric |
| JP2014091890A (en) * | 2012-11-06 | 2014-05-19 | Teijin Ltd | Cloth for thermoconductive activity protection garment and thermoconductive activity protection garment using the same |
Non-Patent Citations (7)
| Title |
|---|
| 上海交通大学: "《机电词典》", 31 October 1991 * |
| 刘国熠等: "消防防护服热湿传递模型理论分析", 《成都纺织高等专科学校学报》 * |
| 吴养会等: "《农村金融计量研究方法及应用》", 30 November 2015 * |
| 周远等: "结合粒子群算法的一类双层纺织材料厚度设计反问题", 《纺织学报》 * |
| 李宁舟等: "高温作业服装材料厚度的设计", 《新型工业化》 * |
| 田苗等: "数值模拟在热防护服装性能测评中的应用", 《纺织学报》 * |
| 瞿伟希等: "《克·火灾与爆燃事故》", 31 May 2014 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111783281A (en) * | 2020-06-12 | 2020-10-16 | 永基众合(山东)智能科技有限公司 | Method for optimizing thickness of multilayer heat-insulating medium |
| CN111783281B (en) * | 2020-06-12 | 2024-05-07 | 永基众合(山东)智能科技有限公司 | Optimization method for thickness of multilayer heat insulation medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | Numerical analysis and optimization study on shell-side performances of a shell and tube heat exchanger with staggered baffles | |
| Liu et al. | Inferring gene regulatory networks using the improved Markov blanket discovery algorithm | |
| Lazzús | Neural network-particle swarm modeling to predict thermal properties | |
| CN110442974B (en) | Method and device for optimizing performance of regenerative chamber of horseshoe flame glass kiln | |
| Zhou et al. | Zone modeling of radiative heat transfer in industrial furnaces using adjusted Monte-Carlo integral method for direct exchange area calculation | |
| Sun et al. | An inverse method for the reconstruction of thermal boundary conditions of semitransparent medium | |
| CN103544356B (en) | Heat treatment furnace heating model calculation method based on graphic processing unit (GPU) | |
| CN116186978B (en) | Data and physical dual-drive temperature field prediction method for complex geometric area of aircraft | |
| Bryden et al. | Optimization of heat transfer utilizing graph based evolutionary algorithms | |
| Malek et al. | A partition of unity finite element method for three-dimensional transient diffusion problems with sharp gradients | |
| Goethals et al. | Sensitivity of night cooling performance to room/system design: surrogate models based on CFD | |
| CN109858068A (en) | Multilyer armor heat transfer model research method based on ant group algorithm | |
| CN114036861A (en) | Three-dimensional temperature field prediction method based on infoGAN | |
| Berger et al. | Intelligent co-simulation: neural network vs. proper orthogonal decomposition applied to a 2D diffusive problem | |
| CN109799007A (en) | Each layer temperature profile analysis model of high-temperature operation clothes based on Fourier's therorem | |
| Zhang et al. | Building energy performance optimization: a new multi-objective particle swarm method | |
| Xu et al. | The improvement of fruit fly optimization algorithm-using bivariable function as example | |
| CN111444619A (en) | Online analysis method and equipment for injection mold cooling system | |
| Xu et al. | Multiple parameter determination in textile material design: A Bayesian inference approach based on simulation | |
| Adomaitis | A reduced-basis discretization method for chemical vapor deposition reactor simulation | |
| Bachir Cherif et al. | Metaheuristics as a solving approach for the infrared heating in the thermoforming process | |
| CN113362911A (en) | IPSO-HKELM-based blast furnace molten iron silicon content prediction method | |
| Li et al. | Research of steel plate temperature prediction based on the improved PSO-ANN algorithm for roller hearth normalizing furnace | |
| Fujarewicz et al. | Spatiotemporal sensitivity of systems modeled by cellular automata | |
| Hadjiski et al. | Application of neural networks to the optimization of the thermal treatment process of wood materials |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190607 |
|
| RJ01 | Rejection of invention patent application after publication |