CN106338431B - Method and its application design method for determining plate-fin heat exchanger mechanics parameter - Google Patents
Method and its application design method for determining plate-fin heat exchanger mechanics parameter Download PDFInfo
- Publication number
- CN106338431B CN106338431B CN201610985659.6A CN201610985659A CN106338431B CN 106338431 B CN106338431 B CN 106338431B CN 201610985659 A CN201610985659 A CN 201610985659A CN 106338431 B CN106338431 B CN 106338431B
- Authority
- CN
- China
- Prior art keywords
- heat exchanger
- plate
- fin heat
- fin
- stress
- 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.)
- Expired - Fee Related
Links
- 238000013461 design Methods 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000002474 experimental method Methods 0.000 claims abstract description 19
- 239000007787 solid Substances 0.000 claims abstract description 9
- 230000035882 stress Effects 0.000 claims description 71
- 239000002184 metal Substances 0.000 claims description 37
- 229910052751 metal Inorganic materials 0.000 claims description 37
- 238000004458 analytical method Methods 0.000 claims description 13
- 230000032683 aging Effects 0.000 claims description 9
- 230000000704 physical effect Effects 0.000 claims description 9
- 230000003321 amplification Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 6
- 239000012528 membrane Substances 0.000 claims description 6
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 6
- 208000025599 Heat Stress disease Diseases 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 4
- 230000006735 deficit Effects 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 claims description 4
- 230000014759 maintenance of location Effects 0.000 claims description 4
- 238000005452 bending Methods 0.000 claims description 3
- 239000012141 concentrate Substances 0.000 claims description 3
- 238000010008 shearing Methods 0.000 claims description 3
- 238000000205 computational method Methods 0.000 abstract description 8
- 230000008569 process Effects 0.000 abstract description 5
- 238000012938 design process Methods 0.000 abstract 1
- 239000000463 material Substances 0.000 description 16
- 238000007789 sealing Methods 0.000 description 8
- 238000000265 homogenisation Methods 0.000 description 7
- 239000002131 composite material Substances 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000005476 soldering Methods 0.000 description 4
- 230000000737 periodic effect Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005219 brazing Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 239000001307 helium Substances 0.000 description 2
- 229910052734 helium Inorganic materials 0.000 description 2
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 description 2
- 229910000679 solder Inorganic materials 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 150000002739 metals Chemical class 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010181 polygamy Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000007634 remodeling Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000008646 thermal stress Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
- G01N2203/0071—Creep
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
- G01N2203/0073—Fatigue
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Heat-Exchange Devices With Radiators And Conduit Assemblies (AREA)
Abstract
Method and its application design method for determining plate-fin heat exchanger mechanics parameter, belongs to technical field of heat exchangers.It is characterized in that：Comprise the following steps：Step a, multiple shape identical plate wing cell elements are divided into by plate fin structure；Step b, plate wing cell element is equivalent into uniform solid plate；Step c, obtains the equivalent mechanical parameter of plate wing cell element, so as to obtain the equivalent mechanical parameter of whole fin heat exchanger core.This is used for the method for determining plate-fin heat exchanger mechanics parameter, and equivalent mechanical parameter, the problem of process is complicated cumbersome can only be obtained by experiment or finite element modelling by solving prior art.The result of the computational methods can be used for the high-temperature strength design of plate-fin heat exchanger, simplify design process；The application design method of this plate-fin heat exchanger can carry out effective life prediction to the plate-fin heat exchanger being on active service under high temperature, alternate load, and effective method is provided for the high-temperature strength design of plate-fin heat exchanger.
Description
Technical field
Method and its application design method for determining plate-fin heat exchanger mechanics parameter, belongs to heat exchanger technology neck
Domain.
Background technology
With the development of science and technology, energy resource consumption is also more and more, and energy-saving and emission-reduction, improving the utilization rate of the energy just turns into people
Focus of attention.Heat transmission equipment, as core a member of high-temperature systems, is not only required nothing more than with efficient heat exchange property, and
It is required that there is compact structure.But, existing heat transmission equipment, mostly shell-and-tube heat exchanger take volume greatly, heat exchange efficiency is low,
It is difficult to meet the requirement in fields such as Aero-Space, HTGR, gas turbines.Plate-fin heat exchanger has compact conformation,
The characteristics of heat exchange efficiency is high, research plate-fin heat exchanger has prospect very much.But it is due to plate-fin heat exchanger periodicity labyrinth
And result in and be difficult to carry out plate-fin heat exchanger finite element modelling problem, and existing method complexity is cumbersome, need to spend a large amount of
Manpower, material resources and financial resources, it is difficult to use, limit with finite element analysis software to plate-fin heat exchanger carry out elevated temperature strength set
Meter, has had a strong impact on the development of plate-fin heat exchanger.
The content of the invention
The technical problem to be solved in the present invention is：The deficiencies in the prior art are overcome to be changed there is provided the convenient plate-fin that calculates of one kind
The method and its application design method for being used to determine plate-fin heat exchanger mechanics parameter of hot device core body equivalent mechanical parameter.
The technical solution adopted for the present invention to solve the technical problems is：This is used to determine plate-fin heat exchanger mechanics parameter
Method, it is characterised in that：Comprise the following steps：
Step a, multiple shape identical plate wing cell elements are divided into by fin heat exchanger core；
Step b, each plate wing cell element is equivalent into uniform solid plate；
Step c, obtains the equivalent mechanical parameter of any one plate wing cell element, so as to obtain whole plate-fin heat exchanger core
The equivalent mechanical parameter of body.
It is preferred that, described fin heat exchanger core includes being provided between flat board and fin, every two pieces of adjacent flat boards
One piece of fin, so as to form multiple runners between per two pieces of adjacent flat boards.
It is preferred that, described equivalent mechanical parameter includes anisotropy equivalent elastic modulus, equivalent shear modulus, Poisson
Than.
It is preferred that, the calculation formula of described anisotropy equivalent elastic modulus is as follows：
Origin is in the midpoint of bottom on front side of using plate-fin heat exchanger, and x is in the direction using in horizontal plane with runner diameter parallel
Axle, is y-axis perpendicular to the direction of runner axis, and vertical direction is that z-axis sets up coordinate system,
,
,
,
Wherein,、、Respectively x-axis, y-axis, the equivalent elastic modulus in z-axis direction,
For the modulus of elasticity of mother metal,
D is the width on the upside of the runner of fin heat exchanger core,
For runner side and vertical plane angle,
For the length of the side of runner,
T is the thickness of fin heat exchanger core flat board,
δ is the thickness of the fin of fin heat exchanger core.
It is preferred that, the calculation formula of Poisson's ratio is as follows：
,
,
,
Wherein,The ratio strained for x-axis direction under being acted in y-axis direction load with y-axis direction,
The ratio strained for x-axis direction under being acted in z-axis direction load with z-axis direction,
The ratio strained for y-axis direction under being acted in z-axis direction load with z-axis direction,
For the Poisson's ratio of mother metal.
It is preferred that, the calculation formula of described equivalent shear modulus is as follows：
Origin is in the midpoint of bottom on front side of using plate-fin heat exchanger, and x is in the direction using in horizontal plane with runner diameter parallel
Axle, is y-axis perpendicular to the direction of runner axis, and vertical direction is that z-axis sets up coordinate system,
,
,
,
Wherein,、The respectively Poisson's ratio and modulus of elasticity of mother metal,
D is the width on the upside of the runner of fin heat exchanger core,
For runner side and vertical plane angle,
For the length of the side of runner,
T is the thickness of fin heat exchanger core flat board,
δ is the thickness of the fin of fin heat exchanger core.
The application design method of plate-fin heat exchanger, it is characterised in that：Comprise the following steps：
Step 1, Preliminary design is carried out to plate-fin heat exchanger structure according to design temperature, design pressure requirement, and clearly
Operating temperature, operation cycle-index and the service life of plate-fin heat exchanger；
Step 2, plate fin structure primary stress analysis is carried out by finite element software, identified sign concentrates position, and determines
Allowable stress；
Step 3, judge whether the stress level at stress concentration position meets following condition：
；；
Wherein,For a membrane stress,For local membrane stress,For primary bending stress,To have with the time
The allowable stress of pass,Span be 1.05 ~ 1.16；
If meeting condition, step 4 is performed；If primary stress evaluation is unsatisfactory for condition, change plate-fin heat exchanger core
The structure of body, sheet metal thickness, return to step 2；
Step 4, under arms under environment, creep rupture experiment and fatigue experiment are carried out to plate fin structure, to aging mother metal
Creep rupture experiment and fatigue experiment are carried out, stress amplification coefficient is calculatedWith strain amplification coefficient, and tied according to experiment
Fruit is modified to the design fatigue curve and creep rupture design curve of mother metal；
,,
Wherein,、The creep rupture strength of mother metal and plate fin structure under respectively identical creep fracture time,
、The macro-strain scope of mother metal and plate fin structure under respectively identical fatigue life；
Step 5, the equivalent of plate fin structure is calculated using the above-mentioned method for being used to determine plate-fin heat exchanger mechanics parameter
Mechanics parameter, while the equivalent thermal physical property parameter of plate fin structure is obtained, so as to carry out heat fatigue finite element to plate-fin heat exchanger
Analysis, draws macro-stress of the fin heat exchanger core along short transverseTime history, calculate fillet at overall strain,
,
Wherein,To analyze the gained range of stress by primary stressThe range of strain drawn；
The macro-stress obtained by Thermal FatigueMaxima and minima difference and plate-fin heat exchanger core
The ratio of body short transverse modulus of elasticity；
Step 6, the fatigue damage of fin heat exchanger core is calculatedAnd creep impairment,
,
Wherein,For fatigue period number,
For on revised design fatigue curve, corresponding fatigue life when range of strain is ε；
,
Wherein,For fatigue period number,
To strain the retention time,
For the macro-stress of t,
Corresponding life-span of creep rupture when for revised creep rupture design curve upper stress being σ；
Step 7, ifLess than 1, then step 8 is performed；IfMore than or equal to 1, then step 1 is performed；
Step 8, plate-fin heat exchanger design is completed.
Compared with prior art, the present invention is had an advantageous effect in that：
1st, the equivalent mechanical parameter computational methods of this plate-fin heat exchanger are by the uneven fin heat exchanger core of structure
Multiple structure same plate wing cell elements are divided into, plate wing cell element is equivalent into uniform solid plate, obtain the Equivalent Mechanical of plate wing cell element
Parameter, so as to obtain the equivalent mechanical parameter of fin heat exchanger core, is solved because plate-fin heat exchanger is periodically multiple
The problem of being difficult to carry out finite element modelling to plate-fin heat exchanger caused by miscellaneous structure, and then facilitate subsequently through finite element
Analysis software carries out Thermal Fatigue to plate-fin heat exchanger.
2nd, anisotropy equivalent elastic modulus, equivalent shear modulus and the Poisson's ratio of fin heat exchanger core pass through
Calculated plate wing cell element is equivalent into uniform solid plate, calculating process is convenient, it is not easy to mistake occur.
3rd, the equivalent mechanical parameter of fin heat exchanger core is provided in the way of analytic expression, and solving can only pass through in the past
The problem of complicated computer simulation or experimental method obtain equivalent parameters, facilitates the calculating of equivalent parameters, greatly improves
The efficiency of plate-fin heat exchanger design.
4th, the application design method of this plate-fin heat exchanger has considered the brazing process of plate-fin heat exchanger, military service ring
Border, invalid position etc. influence, and have carried out equivalent homogenization to fin heat exchanger core, calculate fin heat exchanger core
Equivalent mechanical parameter and equivalent thermal physical property parameter, solving plate-fin heat exchanger can not be direct due to complicated periodic structure
The problem of effectively carrying out high-temperature strength design with finite element software, reason is provided for the high-temperature strength design of plate-fin heat exchanger
By basis, so as to carry out effective life prediction to the plate-fin heat exchanger being on active service under high temperature, alternate load, so as to be
The design of the plate-fin heat exchanger of high-temperature service provides effective method.
Brief description of the drawings
Fig. 1 is the schematic front view of fin heat exchanger core.
Fig. 2 is the schematic front view of plate wing cell element.
Fig. 3 is that creep fatigue life evaluates line chart.
In figure：1st, flat board 2, fin 3, runner.
Embodiment
Fig. 1 ~ 3 are highly preferred embodiment of the present invention, and 1 ~ 3 the present invention will be further described below in conjunction with the accompanying drawings.
As shown in Figure 1 and 2：Plate-fin heat exchanger core body includes being provided between flat board 1 and fin 2, every two pieces of adjacent flat boards 1
Fin 2, is spaced superposition soldering by polylith flat board 1 and fin 2 and forms, so as to form multiple between per two pieces of adjacent flat boards 1
Runner 3, the section of runner 3 is isosceles trapezoid.
The flat board 1 and fin 2 of fin heat exchanger core are formed by soldering, compact conformation, different from conventional material
Material, fin heat exchanger core has a feature of periodicity loose structure, the periodicity of fin heat exchanger core structure and multiple
Polygamy, it is difficult to directly carry out finite element analysis, need to carry out finite element analysis using the method for equivalent homogenization.In order to plate
Fin type heat exchanger core body carries out finite element analysis, introduces homogenization method.
Homogenization method refers to that composite has rule or the structure of approximate regulation, the heterogeneous material of this quite rule
Material assume that as with periodic structure, it should be emphasised that, compared to the size scale of complex, these non-homogeneous materials
Matter is very little.In consideration of it, the material of these types is sometimes referred to as the composite with periodicity microstructure.But
Even if it is also extremely difficult to analyze these boundary value problems comprising a large amount of heterogeneous materials with Modern High-Speed computer.Gram
Take this difficulty to be accomplished by finding a kind of method, replace composite using a kind of equivalent material model, this process is claimed
For homogenization.The essence of homogenization is the composite that equivalent material is replaced to periodic structure, obtains the performance of equivalent material
Parameter, this is the committed step homogenized.
This is used to determine that the method for plate-fin heat exchanger mechanics parameter comprises the following steps：
Step a, multiple shape identical plate wing cell elements are divided into by fin heat exchanger core；
In the present embodiment, structure as shown in Figure 2 is plate wing cell element, so as to by fin heat exchanger core
Regard multiple plate wing cell elements as to combine.
Step b, each plate wing cell element is equivalent into uniform solid plate；
Because the structure of plate wing cell element is not uniform, plate wing cell element is regarded as homogeneous material, i.e., by plate wing cell element etc.
Uniform solid plate is imitated into, uneven plate fin structure is replaced with equivalent solid state plate.
Step c, obtains the equivalent mechanical parameter of any one plate wing cell element, so as to obtain whole plate-fin heat exchanger core
The equivalent mechanical parameter of body.Here equivalent mechanical parameter can also be obtained by finite element analysis software or the method for experiment.
Equivalent mechanical parameter includes anisotropy equivalent elastic modulus, equivalent shear modulus, Poisson's ratio.With plate fin heat-exchanging
The midpoint of bottom is origin on front side of device, and x-axis is in the direction using in horizontal plane with the diameter parallel of runner 3, perpendicular to the axis of runner 3
Direction is y-axis, and vertical direction is that z-axis sets up coordinate system, so as to be counted to the equivalent mechanical parameter of fin heat exchanger core
Calculate.
The computational methods of the equivalent elastic modulus of fin heat exchanger core are as follows：
The equivalent elastic modulus in z-axis direction are calculated, power suffered on flat board 1 and the power suffered by the vertical component of fin 2 is utilized
Balance,
,
The equivalent elastic modulus in x-axis direction are calculated, using equivalent stress and the concept of actual strain,
,
The equivalent elastic modulus in y-axis direction are calculated, using equivalent stress and the concept of actual strain,
,
Wherein,、、Respectively x-axis, y-axis, the equivalent elastic modulus in z-axis direction,
For the modulus of elasticity of mother metal,
D is the width on the upside of the runner 3 of fin heat exchanger core,
For the side of runner 3 and the angle of vertical plane,
For the length of the side of runner 3,
T is the thickness of fin heat exchanger core flat board 1,
δ is the thickness of the fin 2 of fin heat exchanger core.
The computational methods of the Poisson's ratio of fin heat exchanger core are as follows：
Due to flat board booster action, calculate, first calculate, i.e., first calculate the y-axis direction under the effect of x-axis direction load
The ratio strained with x-axis direction, then can draw according to the relation of modulus of elasticity and Poisson's ratio,
I.e.：By,
Draw,
It can similarly obtain：
,
,
Wherein,The ratio strained for x-axis direction under being acted in y-axis direction load with y-axis direction,
The ratio strained for x-axis direction under being acted in z-axis direction load with z-axis direction,
The ratio strained for y-axis direction under being acted in z-axis direction load with z-axis direction,
For the Poisson's ratio of mother metal.
The computational methods of the equivalent shear modulus of fin heat exchanger core are as follows：
Calculate, to isotropism homogeneous material, have,
So as to draw：,
CalculateWith,、It is the ratio of equivalent shearing stress and actual shear strain,
So as to draw：
,
,
Wherein,、The respectively Poisson's ratio and modulus of elasticity of mother metal.
A kind of application design method of plate-fin heat exchanger comprises the following steps：
Step 1, Preliminary design is carried out to plate-fin heat exchanger structure according to design temperature, design pressure requirement, and clearly
Operating temperature, operation cycle-index and the service life of plate-fin heat exchanger；
Run the product that cycle-index is design life and annual machine stop times；Service life is projected life.
Step 2, plate fin structure primary stress analysis is carried out by finite element software, identified sign concentrates position, and determines
The allowable stress at stress concentration position；
Do not consider influence of the soldered seam to structural stress during analysis, and assume that soldering bonding rate is 100%.
Then thermal ageing, Service Environment of the material in brazing process are considered（Helium environment is to the strength of materials in such as HTGR
Influence）, seal structure, be determined via laboratory analysis allowable stress.
Allowable stressFor the allowable stress relevant with the time, allowable stressIncluding fin region allowable stressAnd
Strip of paper used for sealing area allowable stress.Obtaining fin region allowable stressWhen, for soldering high temperature and military service ambient influnence, to aging
Mother metal is under arms in environment（Such as helium environment）Carry out uniaxial tension and creep rupture experiment；Obtaining strip of paper used for sealing area allowable stressWhen, it is necessary to carry out solder tensile strength experiment.Finally according to ASME design criterias and correction result, fin is determined respectively
Area's allowable stressAnd strip of paper used for sealing area allowable stress。
Fin region allowable stress takes the minimum value of the following four factor：
1. design temperature lower yield stress *（1/1.1）* p,
2. the 67%*q of the minimum stress of creep rupture is caused,
3. the 80%*q for the minimum stress that tertiary creep starts,
4. overall strain is reached（Elasticity, plasticity, creep）1% minimum stress,
It is defined in seal structure strain, permissible takes 1/3rd of solder fracture elongation, and strip of paper used for sealing area is allowable
Stress takes the minimum value of the three below factor：
1. mother metal yield stress * under design temperature（1/1.1）* p,
2. permissible * mother metals modulus of elasticity,
3. the 67%*q of the minimum stress of creep rupture is caused,
Wherein, p is thermal ageing mother metal and the ratio of non-thermal ageing mother metal yield stress；
Q is the ratio of aging mother metal creep rupture strength and the creep rupture strength of unaged mother metal.
Step 3, judge whether the stress level at stress concentration position meets following condition：
；；
Wherein,For a membrane stress,For local membrane stress,For primary bending stress,To have with the time
The allowable stress of pass,Span be 1.05 ~ 1.16；
If meeting condition, step 4 is performed；If primary stress evaluation is unsatisfactory for condition, change plate-fin heat exchanger core
The structure of body, sheet metal thickness, return to step 2.
The condition that above-mentioned Rule of judgment is set up is：When fin region stress concentration position stress reaches the allowable stress of fin region
When, plate fin structure failure；When seal structure stress concentration position reaches strip of paper used for sealing area allowable stress, seal structure failure.Plate wing
Any one failure of structure and seal structure thinks that fin heat exchanger core fails, and need to redesign reduction structural stress water
It is flat.
When carrying out stress appraisal, fin region allowable stress is utilized respectivelyAnd the allowable stress in strip of paper used for sealing areaEnter
Row is evaluated, when fin region stress meets allowable stressAnd strip of paper used for sealing area stress meets allowable stressWhen, that is, once should it make
Power evaluates qualified judgement, when fin region stress is unsatisfactory for allowable stressOr the stress in strip of paper used for sealing area is unsatisfactory for allowable stress
When, then make the judgement for being unsatisfactory for primary stress evaluation.
Step 4, under arms under environment, creep rupture experiment and fatigue experiment are carried out to plate fin structure, to aging mother metal
Creep rupture experiment and fatigue experiment are carried out, stress amplification coefficient is calculatedWith strain amplification coefficient, and tied according to experiment
Fruit is modified to the design fatigue curve and creep rupture design curve of mother metal；
,,
Wherein,、The creep rupture strength of mother metal and plate fin structure under respectively identical creep fracture time,
、The macro-strain scope of mother metal and plate fin structure under respectively identical fatigue life；
Step 5, the equivalent of plate fin structure is calculated using the above-mentioned method for being used to determine plate-fin heat exchanger mechanics parameter
Mechanics parameter, while the equivalent thermal physical property parameter of plate fin structure is obtained, so as to carry out heat fatigue finite element fraction to plate-fin heat exchanger
Analysis.Wherein equivalent mechanical parameter includes anisotropy equivalent elastic modulus, equivalent shear modulus, Poisson's ratio.
The acquisition methods of equivalent thermal physical property parameter are identical with the acquisition methods of equivalent mechanical parameter, by the uneven plate of structure
Fin type heat exchanger core body is divided into multiple structure same plate wing cell elements, and plate wing cell element is equivalent into uniform solid plate, obtains plate wing
The equivalent thermal physical property parameter of cell element, so as to obtain the equivalent thermal physical property parameter of fin heat exchanger core.Equivalent heat physical property is joined
Number includes Equivalent Thermal Conductivities, fiber yarn, equivalent density and equivalent specific heat.
The computational methods of the Equivalent Thermal Conductivities of fin heat exchanger core are as follows：
Minimum thermal resistance force method then, is also referred to as rule in parallel, when heat is transmitted in object, and hot-fluid can lead to along resistance is minimum
Road is transmitted, or passage is in minimum thermal resistive state when flowing through orientation heat flow, and the entire thermal resistance of respective channel is minimum thermal resistance,
Also referred to as equivalent thermal resistance.Equivalent Thermal Conductivities rule, as long as the cell cube of composite is having equal ratio equivalent thermal resistance with overall,
No matter cell cube size, when only considering heat transfer, this cell cube is equal with overall Equivalent Thermal Conductivities.
As known from the above, it is desirable to the overall thermal conductivity factor of fin heat exchanger core, it is only necessary to seek plate wing cell element
Equivalent Thermal Conductivities.
,
,
,
Wherein,、、Respectively x-axis direction, y-axis direction, the Equivalent Thermal Conductivities in z-axis direction,
、The respectively thermal conductivity factor of mother metal and air.
The computational methods of the fiber yarn of fin heat exchanger core are as follows：
Because the flat board 1 of the upper side and lower side of a plate wing cell element can free wxpansion, the part of horizontal part of fin 2 and inclination
Partially due to the swell increment in z-axis direction is different, there is interaction, so as to draw：
,
Because flat board 1 is identical with the thermal expansion amount in x-axis direction in y-axis direction with fin 2, so having
,
,
Wherein,For equivalent mother metal thermal coefficient of expansion.
The equivalent density of fin heat exchanger core and the computational methods of equivalent specific heat are as follows：
,
,
,
,
Wherein,、Respectively mother metal volume fraction and volume of air fraction,
、For mother metal density and atmospheric density,
、The respectively specific heat of mother metal and air,
、Respectively equivalent specific heat and equivalent density.
Using the equivalent mechanical parameter and equivalent thermal physical property parameter of calculating, by finite element analysis software, with homogenization
Method carries out the anisotropic elasticity analysis of heat fatigue, from the result of thermal-stress analysis, draws plate fin structure Direction of superposition, i.e. z
The macro-stress of direction of principal axisTime history, the difference of macro-stress maxima and minima and z-axis direction modulus of elasticity
Ratio, as range of strain；By primary stress scopeThe range of strain drawn, so as to calculate total at fillet
Strain,
。
Step 6, the fatigue damage of fin heat exchanger core is calculatedAnd creep impairment,
,
Wherein,For fatigue period number,
For on revised design fatigue curve, corresponding fatigue life when range of strain is ε；
,
Wherein,For fatigue period number,
To strain the retention time,
For the macro-stress of t,
Corresponding life-span of creep rupture when for revised creep rupture design curve upper stress being σ.
Fatigue damageFor fin heat exchanger core each point in design fatigue curve of the aging mother metal in maximum temperature
Period is run with allowing the ratio of period.Calculate creep impairmentWhen, during due to the retention time, it may occur that stress relaxation,
So can preferably draw the stress relaxation curve of structure, i.e.,Change curve.
Step 7, ifLess than 1, then step 8 is performed；IfMore than or equal to 1, then step 1 is performed；
According to creep-fatigue assessment of impairments criterion in ASME, as shown in figure 3, being damaged according to the total creep calculatedWith it is total
Fatigue damage, withFor abscissa,For ordinate,For the envelope of crack initiation, if
Less than 1, i.e.,Under envelope, then illustrate that plate-fin heat exchanger under design temperature and pressure, entirely designs the longevity
It will not be failed in life, meet design requirement, perform step 8.IfMore than or equal to 1, illustrate that plate-fin heat exchanger is discontented with
Sufficient design requirement, now needs improved structure, more conversion materials, and military service pressure, temperature are reduced under enabled condition, step is re-executed
1, untilUnder envelope, the high-temperature strength design of plate-fin heat exchanger is completed.
Step 8, plate-fin heat exchanger design is completed.
The above described is only a preferred embodiment of the present invention, being not the limitation for making other forms to the present invention, appoint
What those skilled in the art changed or be modified as possibly also with the technology contents of the disclosure above equivalent variations etc.
Imitate embodiment.But it is every without departing from technical solution of the present invention content, the technical spirit according to the present invention is to above example institute
Any simple modification, equivalent variations and the remodeling made, still fall within the protection domain of technical solution of the present invention.
Claims (4)
1. the application design method of plate-fin heat exchanger, it is characterised in that：Comprise the following steps：
Step 1, Preliminary design, and clear and definite plate wing are carried out to plate-fin heat exchanger structure according to design temperature, design pressure requirement
Operating temperature, operation cycle-index and the service life of formula heat exchanger；Fin heat exchanger core includes flat board（1）And fin
（2）, per two pieces of adjacent flat boards（1）Between provided with one piece of fin（2）, so that per two pieces of adjacent flat boards（1）Between formed
Multiple runners（3）；
Step 2, plate fin structure primary stress analysis is carried out by finite element software, identified sign concentrates position, and determines allowable
Stress；
Step 3, judge whether the stress level at stress concentration position meets following condition：
；；
Wherein,For a membrane stress,For local membrane stress,For primary bending stress,To be relevant with the time
Allowable stress,Span be 1.05 ~ 1.16；
If meeting condition, step 4 is performed；If primary stress evaluation is unsatisfactory for condition, change fin heat exchanger core
Structure, sheet metal thickness, return to step 2；
Step 4, under arms under environment, creep rupture experiment and fatigue experiment is carried out to plate fin structure, aging mother metal is also carried out
Creep rupture experiment and fatigue experiment, calculate stress amplification coefficientWith strain amplification coefficient, and according to experimental result pair
The design fatigue curve and creep rupture design curve of mother metal are modified；
,,
Wherein,、The creep rupture strength of mother metal and plate fin structure under respectively identical creep fracture time,
、The macro-strain scope of mother metal and plate fin structure under respectively identical fatigue life；
Step 5, calculate plate fin structure equivalent mechanical parameter, equivalent mechanical parameter include anisotropy equivalent elastic modulus, etc.
Imitate modulus of shearing, Poisson's ratio；The calculation procedure of equivalent mechanical parameter is as follows：
Step a, multiple shape identical plate wing cell elements are divided into by fin heat exchanger core；
Step b, each plate wing cell element is equivalent into uniform solid plate；
Step c, obtains the equivalent mechanical parameter of any one plate wing cell element, so as to obtain whole fin heat exchanger core
Equivalent mechanical parameter；
The equivalent thermal physical property parameter of plate fin structure is obtained simultaneously, so that heat fatigue finite element analysis is carried out to plate-fin heat exchanger,
Draw macro-stress of the fin heat exchanger core along short transverseTime history, calculate fillet at overall strain,
,
Wherein,To analyze the gained range of stress by primary stressThe range of strain drawn；
The macro-stress obtained by Thermal FatigueMaxima and minima difference and fin heat exchanger core it is high
Spend the ratio of direction modulus of elasticity；
Step 6, the fatigue damage of fin heat exchanger core is calculatedAnd creep impairment,
,
Wherein,For fatigue period number,
For on revised design fatigue curve, corresponding fatigue life when range of strain is ε；
,
Wherein,For fatigue period number,
To strain the retention time,
For the macro-stress of t,
For on revised creep rupture design curve, corresponding life-span of creep rupture when stress is σ；
Step 7, ifLess than 1, then step 8 is performed；IfMore than or equal to 1, then step 1 is performed；
Step 8, plate-fin heat exchanger design is completed.
2. the application design method of plate-fin heat exchanger according to claim 1, it is characterised in that：Described anisotropy
The calculation formula of equivalent elastic modulus is as follows：
Origin is in midpoint using bottom on front side of plate-fin heat exchanger, with horizontal plane and runner（3）The direction of diameter parallel is x
Axle, perpendicular to runner（3）The direction of axis is y-axis, and vertical direction is that z-axis sets up coordinate system,
,
,
,
Wherein,、、Respectively x-axis, y-axis, the equivalent elastic modulus in z-axis direction,
For the modulus of elasticity of mother metal,
D is the runner of fin heat exchanger core（3）The width of upside,
For runner（3）Side and the angle of vertical plane,
For runner（3）Side length,
T is fin heat exchanger core flat board（1）Thickness,
δ is the fin of fin heat exchanger core（2）Thickness.
3. the application design method of plate-fin heat exchanger according to claim 2, it is characterised in that：The calculating of Poisson's ratio is public
Formula is as follows：
,
,
,
Wherein,The ratio strained for x-axis direction under being acted in y-axis direction load with y-axis direction,
The ratio strained for x-axis direction under being acted in z-axis direction load with z-axis direction,
The ratio strained for y-axis direction under being acted in z-axis direction load with z-axis direction,
For the Poisson's ratio of mother metal.
4. the application design method of plate-fin heat exchanger according to claim 1, it is characterised in that：Described equivalent shearing
The calculation formula of modulus is as follows：
Origin is in midpoint using bottom on front side of plate-fin heat exchanger, with horizontal plane and runner（3）The direction of diameter parallel is x
Axle, perpendicular to runner（3）The direction of axis is y-axis, and vertical direction is that z-axis sets up coordinate system,
,
,
,
Wherein,、The respectively Poisson's ratio and modulus of elasticity of mother metal,
D is the runner of fin heat exchanger core（3）The width of upside,
For runner（3）Side and the angle of vertical plane,
For runner（3）Side length,
T is fin heat exchanger core flat board（1）Thickness,
δ is the fin of fin heat exchanger core（2）Thickness.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610985659.6A CN106338431B (en) | 2016-11-09 | 2016-11-09 | Method and its application design method for determining plate-fin heat exchanger mechanics parameter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610985659.6A CN106338431B (en) | 2016-11-09 | 2016-11-09 | Method and its application design method for determining plate-fin heat exchanger mechanics parameter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106338431A CN106338431A (en) | 2017-01-18 |
| CN106338431B true CN106338431B (en) | 2017-09-29 |
Family
ID=57841204
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610985659.6A Expired - Fee Related CN106338431B (en) | 2016-11-09 | 2016-11-09 | Method and its application design method for determining plate-fin heat exchanger mechanics parameter |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106338431B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106503390B (en) * | 2016-11-09 | 2017-08-25 | 中国石油大学(华东) | A kind of creep fatigue strength design of plate-fin heat exchanger |
| CN107842979B (en) * | 2017-11-30 | 2023-07-14 | 苏州惠林节能材料有限公司 | Energy-saving plate-fin type heat recovery core and manufacturing method thereof |
| CN108319798B (en) * | 2018-03-13 | 2021-11-23 | 合肥通用机械研究院有限公司 | Design method of large-scale heat exchanger tube plate under complex working conditions |
| CN116090307A (en) * | 2023-02-15 | 2023-05-09 | 哈尔滨锅炉厂有限责任公司 | Finite element modeling method of PCHE core based on orthotropic constitutive |
| CN116757027B (en) * | 2023-06-08 | 2024-03-26 | 西安交通大学 | Construction method and system of general database of equivalent mechanical properties of heat exchanger channel |
| CN116822167B (en) * | 2023-06-08 | 2024-04-16 | 西安交通大学 | Heat exchanger thermal coupling performance multi-scale analysis method, system, medium and equipment |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005227065A (en) * | 2004-02-12 | 2005-08-25 | Hitachi Ltd | Display method for estimated life duration of structure, method of creating life duration distribution, and image display |
| CN1943973A (en) * | 2006-11-07 | 2007-04-11 | 管保清 | Method for producing aluminum plate-fin heat exchanger |
| CN103761365A (en) * | 2013-12-28 | 2014-04-30 | 合肥通用机械研究院 | High-temperature pressure vessel creep fatigue strength design method based on service life |
| CN205607214U (en) * | 2016-03-16 | 2016-09-28 | 中国石油大学(华东) | Chevron shape trompil corrugated fin type finned plate heat exchanger |
-
2016
- 2016-11-09 CN CN201610985659.6A patent/CN106338431B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005227065A (en) * | 2004-02-12 | 2005-08-25 | Hitachi Ltd | Display method for estimated life duration of structure, method of creating life duration distribution, and image display |
| CN1943973A (en) * | 2006-11-07 | 2007-04-11 | 管保清 | Method for producing aluminum plate-fin heat exchanger |
| CN103761365A (en) * | 2013-12-28 | 2014-04-30 | 合肥通用机械研究院 | High-temperature pressure vessel creep fatigue strength design method based on service life |
| CN205607214U (en) * | 2016-03-16 | 2016-09-28 | 中国石油大学(华东) | Chevron shape trompil corrugated fin type finned plate heat exchanger |
Non-Patent Citations (1)
| Title |
|---|
| 《钎料厚度对不锈钢板翅结构抗拉强度影响的有限元分析与试验研究》;巩建鸣等;《真空科学与技术学报》;20101031;第30卷(第5期);第561-566页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106338431A (en) | 2017-01-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106503390B (en) | A kind of creep fatigue strength design of plate-fin heat exchanger | |
| CN106338431B (en) | Method and its application design method for determining plate-fin heat exchanger mechanics parameter | |
| CN106355042B (en) | The homogenization design method of plate-fin heat exchanger | |
| Wang et al. | State of the art in the aeroelasticity of wind turbine blades: Aeroelastic modelling | |
| Hu et al. | Integrating variable wind load, aerodynamic, and structural analyses towards accurate fatigue life prediction in composite wind turbine blades | |
| CN107451308B (en) | Multi-scale calculation method for equivalent heat conduction coefficient of complex composite material structure | |
| CN107451309B (en) | Method for multi-scale calculation of equivalent thermal expansion coefficient of complex composite material structure | |
| Liao et al. | Residual fatigue life analysis and comparison of an aluminum lithium alloy structural repair for aviation applications | |
| CN104711923A (en) | Method for simulating rigidity of suspension bridge aeroelastic model main beam through trapezoidal springs | |
| CN108009311A (en) | A kind of creep Parameters of constitutive model recognition methods for creep test | |
| Ahn et al. | Analysis of cracked aluminum plates with one-sided patch repair using p-convergent layered model | |
| Li et al. | Life prediction model of creep-rupture and creep-buckling of a pyramidal lattice truss panel structure by analytical and finite element study | |
| Fossum et al. | Design and fatigue performance of large utility-scale wind turbine blades | |
| Haghdoust et al. | Investigation of shape memory alloy embedded wind turbine blades for the passive control of vibrations | |
| Herath et al. | Development of a shape-adaptive composite propeller using bend-twist coupling characteristics of composites | |
| Lotfi et al. | 3D fluid-structure interaction (FSI) simulation of new type vortex generators in smooth wavy fin-and-elliptical tube heat exchanger | |
| Sun et al. | Study on impact resistance and parameter optimization of patch‐repaired plain woven composite based on multi‐scale analysis | |
| Heinrich et al. | A numerical approach to simulate printed electronics on fiber reinforced composites | |
| Liu et al. | Impact of saddle membrane structure by hail with combined particle sizes: Numerical simulation and experimental investigation | |
| CN202381164U (en) | Intercooler with structurally novel core | |
| Tu et al. | Creep of brazed plate-fin structures in high temperature compact heat exchangers | |
| Jang et al. | Computation of energy release rates for composite beam through cross-sectional analysis and virtual crack closure technique | |
| CN103593528A (en) | Equivalent standard solid model of viscoelastic damper | |
| Zhang et al. | STRESS AND FATIGUE ANALYSIS OF AIR-TO-REFRIGERANT HEAT EXCHANGERS WITH NON-ROUND TUBE SHAPES | |
| Jin et al. | Life prediction of wind turbine blades |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170929 |