CN105572161B - A kind of method and used test device of non-steady state Determination of conductive coefficients - Google Patents
A kind of method and used test device of non-steady state Determination of conductive coefficients Download PDFInfo
- Publication number
- CN105572161B CN105572161B CN201610012630.XA CN201610012630A CN105572161B CN 105572161 B CN105572161 B CN 105572161B CN 201610012630 A CN201610012630 A CN 201610012630A CN 105572161 B CN105572161 B CN 105572161B
- Authority
- CN
- China
- Prior art keywords
- temperature
- temperature sensor
- steady state
- thermal coefficient
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/18—Investigating or analyzing materials by the use of thermal means by investigating thermal conductivity
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)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
A kind of method of non-steady state Determination of conductive coefficients, this approach includes the following steps：1）Thermal perturbation is applied to sample, obtains the non-steady state heat transfer temperature data of material internal；2）By temperature data, function of many variables fitting is carried out to non-steady state thermal coefficient using least square method；3）Based on least square finite element method according to the non-steady state process of thermal coefficient of the temperature value Inversion Calculation material of each node, non-steady state thermal coefficient is obtained.The method and used test device of a kind of non-steady state Determination of conductive coefficients provided by the invention, can solve the problems, such as traditional steady state test method can not Accurate Determining go out the non-steady state thermal coefficient of material, under the premise of not destroying sample, obtain the temperature data of the material internal specified point in non-steady state heat transfer process, and then the thermal coefficient of material is obtained using least square finite element reverse calculation algorithms by these temperature datas put, test result is accurate.
Description
Technical field
The present invention relates to the measurement of the non-steady state thermal coefficient of multiphase porous media in Engineering Thermodynamics field, especially one
The method and used test device of the non-steady state Determination of conductive coefficients of kind.
Background technology
It is universal thermodynamics Heat Conduction Phenomenon that heat is transmitted to lower from temperature higher position, and the size of substance thermal conductivity is usual
It can be characterized with thermal coefficient.Thermal coefficient refers to the material of 1m thickness under the conditions of steady heat transfer, and the temperature difference of both side surface is 1
Degree, in 1 second, by 1 square metre of area transmit heat, unit be watt/ meter Du (W/ (mK)), to indicate unit
The power of thermal conductivity under gradient.It is exploited in Temperature Controlling of Mass Concrete, temperature stress calculation, underground heat, and High-temperature cooling etc. is many
Field has a wide range of applications.The Measured Results of Thermal Conductivity method of traditional common is mostly to measure steady state thermal coefficient, based in Fu
The test philosophy of leaf heat transfer law,The thermal coefficient measured is a certain range of average value.
And hydraulic large-volume concrete temperature control or underground heat exploitation in material medium majority be comprising hole, moisture and
The mixing porous media of other impurity, the ingredients of these materials real-time change at any time, causes thermal coefficient nor steady state to change,
Using traditional steady state test method can not Accurate Determining go out the non-steady state thermal coefficient of material.
Invention content
Technical problem to be solved by the invention is to provide a kind of methods and used test of non-steady state Determination of conductive coefficients
Device, can solve the problems, such as traditional steady state test method can not Accurate Determining go out the non-steady state thermal coefficient of material, not
Under the premise of destroying sample, the temperature data of the material internal specified point in non-steady state heat transfer process is obtained, and then passes through this
The temperature data put a bit obtains the thermal coefficient of material using least square finite element reverse calculation algorithms, and test result is accurate.
In order to solve the above technical problems, the technical solution adopted in the present invention is：A kind of non-steady state Determination of conductive coefficients
Method, this approach includes the following steps：
1) thermal perturbation is applied to sample, obtains the non-steady state heat transfer temperature data of material internal；
2) by temperature data, function of many variables fitting is carried out to non-steady state thermal coefficient using least square method；
3) based on least square finite element method according to the non-perseverance of thermal coefficient of the temperature value Inversion Calculation material of each node
State process obtains non-steady state thermal coefficient.
Sample and multiple sensors are arranged symmetrically, for sample application thermal perturbation, the circulating-heating in certain temperature range,
It allows temperature to be uniformly diffused into sample, then acquires the temperature data under different moments different radii, obtain the non-steady state heat of material internal
Transition temperature data.
The non-steady state function of many variables fitting for the least square method that step 2) uses.
Least square method by minimize error quadratic sum come search with give Data Matching degree best function.It is thought
Think to be exactly to find out to be fitted best function with real data, i.e., the square distance of point and observation point on function and reach minimum,
The function balances the matching degree between each equation, is an optimal solution of the overall situation so that geometrical curve described in function
It is more close with actual conditions.Specific method is, if variable y and n variable x1,x2,…,xn(n >=1) inner link is linear
, that is, there is following relational expressionIf xjIth measured value be xij, corresponding functional value is yi(i=1,
2 ..., m), then sum of square of deviations isThe inverse of non-steady state thermal coefficient is just
It is a kind of function of many variables fitting, if thermal coefficient K=α1+α2t+α3t2+...+αntn-1, with α1, α2... ..., αnA undetermined parameter
It is related.In order to make s be minimized to obtain equation group：
I.e.
By test data (xij, yi) (2) formula is substituted into, up to a0,a1,…,am。Xij,yij
The step of step 3) is：
3-1) in solid thermal conduction problem, it is assumed that density p and specific heat capacity c are not changed over time, arbitrary in computational domain
At point, meet the equation of heat conduction
In formula, T is temperature (DEG C)；K is thermal coefficient (kJ/mh DEG C)；For endogenous pyrogen (kJ)；T times (h), formula (3)
The finite element differential form without endogenous pyrogen that is steady state heat transfer；
3-2) to step 3-1) formula (3) that obtains is discrete using the golden method of gal the Liao Dynasty in spatial domain, in time-domain using poor
Divide method, it is as follows to obtain finite element scheme：
{ T }=(1- ξ) { Tn}+ξ{Tn-1} (5)
Wherein, [C] is heat capacity matrix；[H] is heat transfer matrix；{ P } temperature load array；{ T } node temperature array；Node temperature is to time-derivative array；
3-3) assume that thermal coefficient and the relational expression of time are：
K=α1+α2t+α3t2+...+αntn-1 (7)
Fundamental unknown variables are transformed to each undetermined parameter (α of description thermal coefficient by node temperature value1,α2,…,αn), therefore formula
(4) it is transformed to：
[L] { α }={ f } (8)
Wherein, matrix [L] and { P } are integrated by the corresponding matrix element of unit, i.e.,
In formula,It is that unit treats the contribution for seeking unknown quantity coefficient matrix after converting fundamental unknown variables, is provided by formula (11)：
3-4) formula (8) is with parameter alpha1,α2,…,αnFor the super equation group of fundamental unknown variables, sought entirely with least square method
Office's optimal solution：
ψ=Wk(Lkjαj-fk)2 (12)
In formula, WkThe weight of k-th of equation is obtained since formula (8) is to include the super equation group of multiple unknown numbers when being cumulative
The solution arrived is not unique, and formula (12) is turned to because this introduces least square method to formula (8)；
Formula (13) is the full rank equation group obtained to each coefficient derivation to be asked, solve this equation group can ask to obtain it is undetermined
Parameter alphaiValue, then non-steady state thermal coefficient K can be calculated by formula (7).
This method further includes step 4)：It writes programed temperature acquisition to calculate in real time with thermal coefficient, obtains non-steady state
The graph of relation of thermal coefficient and time.
A kind of test device tested using the above method in non-steady state thermal coefficient for acquiring specimen temperature data,
It includes Die tester, temperature control system and data collecting system；
In Die tester, circular shaped foam plate is arranged in die trial bottom, and die trial center is inserted with steel pipe；
In temperature control system, the single-ended heating rod being connect with relay is arranged in steel pipe center, and relay connects with power supply
It connects；
Data collecting system includes multiple temperature sensor groups, in multiple temperature sensor groups, the first temperature sensor cloth
It is placed in steel duct, second temperature sensor and third temperature sensor are arranged symmetrically the temperature sensor group of composition, the 4th temperature
Temperature sensor group, the 6th temperature sensor and the 7th temperature that degree sensor and the 5th temperature sensor are arranged symmetrically composition pass
The temperature sensor group that sensor is arranged symmetrically composition is arranged in inside sample respectively from the inside to the outside along different radii, multiple temperature
Each temperature sensor of sensor group is connect with temperature collecting module respectively, and temperature collecting module connects with power supply and computer respectively
It connects；
Each temperature sensor and the first temperature sensor of multiple temperature sensor groups are passed through by each lead and are opened in circle
The mode in the hole of shape cystosepiment is fixed on circular shaped foam plate；
The bottom seal of steel pipe, steel pipe is interior to inject cold water, and the first temperature sensor is for monitoring the interior water temperature of pipe.
Second temperature sensor and third temperature sensor are symmetrically arranged in 30 centimeters of radius of sample；
4th temperature sensor and the 5th temperature sensor are arranged in 43.5 centimeters of radius of sample；
6th temperature sensor and the 7th temperature sensor are symmetrically arranged in 72.5 centimeters of radius of sample.
A method of applied automatically to what specimen temperature disturbed using the realization of above-mentioned test device, this method is：
The original state of relay is to disconnect when on-test is arranged, and temperature collecting module acquires test coupon centre distance
The temperature of farthest the 6th temperature sensor or the 7th temperature sensor is T, works as T<Relay is closed when the minimum temperature of setting,
Single-ended heating rod powers on, and begins to warm up the water in steel pipe；Then the temperature of each temperature point is begun to ramp up immediately, when the 6th
The temperature T of temperature sensor or the 7th temperature sensor>When the maximum temperature of setting, relay disconnects, single-ended heating rod and electricity
The channel in source is cut off, and stops heating, as a heat cycles, realizes the automatic application disturbed to specimen temperature.
In order to control the time of one cycle, and improve the discrimination of each temperature point temperature value, setting it is minimum
Temperature is that room temperature adds 5 DEG C, and maximum temperature is that room temperature adds 12 DEG C.
The method and used test device of a kind of non-steady state Determination of conductive coefficients provided by the invention, advantage is such as
Under：
1, can solve the problems, such as traditional steady state test method can not Accurate Determining go out the non-steady state thermal coefficient of material,
Industry blank is compensated for, under the premise of not destroying sample, obtains the material internal specified point in non-steady state heat transfer process
Temperature data, and then the heat conduction system of material is obtained using least square finite element reverse calculation algorithms by the temperature data of these points
Number, test result are accurate.
2, each temperature sensor and the first temperature sensor of multiple temperature sensor groups are passed through by each lead and are opened in
The mode in the hole of circular shaped foam plate is fixed on the mode on circular shaped foam plate, plays fixed sensor, and reduces heat and dissipate
The effect of mistake, preferably heat-insulation and heat-preservation make whole system form One-dimensional heat transfer so that the test data measured is more accurate.
The present invention obtains SanXia University postgraduate scientific research and innovation fund assistance (project number 3) in 2015.
Description of the drawings
The invention will be further described with reference to the accompanying drawings and examples：
Fig. 1 is the structural schematic diagram of test device of the present invention；
Fig. 2 is the arrangement signal of each temperature sensor of multiple temperature sensor groups used by test device of the present invention
Figure；
Fig. 3 is that the present invention realizes the control flow chart applied automatically disturbed to specimen temperature using test device；
Fig. 4 writes program Inversion Calculation material thermal coefficient in the method for the non-steady state Determination of conductive coefficients of the present invention is non-
The block diagram of steady state process；
Passes of the Fig. 5 for step 4) non-steady state thermal coefficient and time in the method for the non-steady state Determination of conductive coefficients of the present invention
It is curve graph.
Specific implementation mode
Embodiment one
A kind of method of non-steady state Determination of conductive coefficients, this approach includes the following steps：
1) thermal perturbation is applied to sample, obtains the non-steady state heat transfer temperature data of material internal；
2) by temperature data, function of many variables fitting is carried out to non-steady state thermal coefficient using least square method；
3) based on least square finite element method according to the non-perseverance of thermal coefficient of the temperature value Inversion Calculation material of each node
State process obtains non-steady state thermal coefficient.
Sample and multiple sensors are arranged symmetrically, for sample application thermal perturbation, the circulating-heating in certain temperature range,
It allows temperature to be uniformly diffused into sample, then acquires the temperature data under different moments different radii, obtain the non-steady state heat of material internal
Transition temperature data.
The step of step 3) is：
3-1) in solid thermal conduction problem, it is assumed that density p and specific heat capacity c are not changed over time, arbitrary in computational domain
At point, meet the equation of heat conduction
In formula, T is temperature (DEG C)；K is thermal coefficient (kJ/mh DEG C)；For endogenous pyrogen (kJ)；T times (h), formula (3)
The finite element differential form without endogenous pyrogen that is steady state heat transfer；
3-2) to step 3-1) formula (3) that obtains is discrete using the golden method of gal the Liao Dynasty in spatial domain, in time-domain using poor
Divide method, it is as follows to obtain finite element scheme：
{ T }=(1- ξ) { Tn}+ξ{Tn-1} (5)
Wherein, [C] is heat capacity matrix；[H] is heat transfer matrix；{ P } temperature load array；{ T } node temperature array；Node temperature is to time-derivative array；
3-3) assume that thermal coefficient and the relational expression of time are：
K=α1+α2t+α3t2+...+αntn-1 (7)
Fundamental unknown variables are transformed to each undetermined parameter (α of description thermal coefficient by node temperature value1,α2,…,αn), therefore formula
(4) it is transformed to：
[L] { α }={ f } (8)
Wherein, matrix [L] and { P } are integrated by the corresponding matrix element of unit, i.e.,
In formula,It is that unit treats the contribution for seeking unknown quantity coefficient matrix after converting fundamental unknown variables, is provided by formula (11)：
3-4) formula (8) is with parameter alpha1,α2,…,αnFor the super equation group of fundamental unknown variables, sought entirely with least square method
Office's optimal solution：
ψ=Wk(Lkjαj-fk)2 (12)
In formula, WkThe weight of k-th of equation is obtained since formula (8) is to include the super equation group of multiple unknown numbers when being cumulative
The solution arrived is not unique, and formula (12) is turned to because this introduces least square method to formula (8)；
Formula (13) is the full rank equation group obtained to each coefficient derivation to be asked, solve this equation group can ask to obtain it is undetermined
Parameter alphaiValue, then non-steady state thermal coefficient K can be calculated by formula (7).
This method further includes step 4)：It writes programed temperature acquisition to calculate in real time with thermal coefficient, obtains non-steady state
The graph of relation of thermal coefficient and time, as shown in figure 5, it can be seen that concrete thermal coefficient with the age of concrete and
Non- steady state variation, with the increase of age number of days, concrete thermal coefficient constantly reduces.
Embodiment two
As depicted in figs. 1 and 2, a kind of to be tested in non-steady state thermal coefficient using the above method for acquiring specimen temperature number
According to test device,
It includes Die tester, temperature control system and data collecting system；
In Die tester, circular shaped foam plate 5 is arranged in 6 bottom of die trial, and 6 center of die trial is inserted with steel pipe 9；
In temperature control system, the single-ended heating rod 10 being connect with relay 4 is arranged in 9 center of steel pipe, relay 4 and electricity
Source 3 connects；
Data collecting system includes multiple temperature sensor groups 8, in multiple temperature sensor groups 8, the first temperature sensor
11 are arranged in inside steel pipe 9, and second temperature sensor 12 and third temperature sensor 13 are arranged symmetrically the temperature sensor of composition
Group, the 4th temperature sensor 14 and the 5th temperature sensor 15 are arranged symmetrically the temperature sensor group of composition, the 6th temperature sensing
The temperature sensor group that device 16 and the 7th temperature sensor 17 are arranged symmetrically composition is arranged in inside sample 7 edge from the inside to the outside respectively
At different radii, each temperature sensor of multiple temperature sensor groups 8 is connect with temperature collecting module 2 respectively, temperature acquisition mould
Block 2 is connect with power supply 3 and computer 1 respectively；
Each temperature sensor of multiple temperature sensor groups 8 and the first temperature sensor 11 are passed through by each lead to be opened in
The mode in the hole of circular shaped foam plate 5 is fixed on circular shaped foam plate 5；
The bottom seal of steel pipe 9, steel pipe 9 is interior to inject cold water, and the first temperature sensor 11 is for monitoring the interior water temperature of pipe.
Second temperature sensor 12 and third temperature sensor 13 are symmetrically arranged in 30 centimeters of radius of sample 7；
4th temperature sensor 14 and the 5th temperature sensor 15 are arranged in 43.5 centimeters of radius of sample 7；
6th temperature sensor 16 and the 7th temperature sensor 17 are symmetrically arranged in 72.5 centimeters of radius of sample 7.
Cystosepiment 5 is EPS circular shaped foams plate 5, and steel pipe 9 is stainless-steel pipe, and single-ended heating rod 10 heats for single-end DC
Stick, power supply 3 are Double-way dc voltage-stabilizing source, and temperature collecting module 2 is temperature acquisition card.
The temperature of measuring point where temperature acquisition card acquires each temperature sensor adopts control time interval and is set as 60s, i.e., every
60s acquires a temperature data, and is shown in the row of real time temperature one.Automated collection systems use C# window writing routines, temperature acquisition
Module 2, relay 4 are communicated by RS-232 serial ports with computer 1, therefore respectively there are one serial ports, and are controlled by baud rate
The frequency of communication.
While filling sample to die trial, embedding each temperature sensing at temperature measuring point is preset in steel pipe 9 and in sample 7
Device.To ensure temp probe burial place accurately and not shifting, first filling material is tamped, temp probe is fixed on finger
Positioning is set, then fills upper layer of material, and checks the position of temp probe at any time, if any displacement, in time correction or in detail note case.Most
Sample 7 is placed in vibration compacting on the platform that vibrates afterwards.
Temperature control system is auto temperature controlled system, during primary heating recycles, when the temperature of temperature sensor
Value is within certain period when kept stable, it is believed that system reaches Heat Conduction process.Multiple heat cycles constitute non-perseverance
State heat transfer process.
Embodiment three
As shown in figure 3, a kind of realizing the method applied automatically disturbed to specimen temperature, the party using above-mentioned test device
Method is：
The original state of relay 4 is to disconnect when on-test is arranged, and temperature collecting module 2 acquires 7 center of test coupon
The temperature of the 6th farthest temperature sensor 16 of distance or the 7th temperature sensor 17 is T, works as T<When the minimum temperature of setting after
Electric appliance 4 is closed, and single-ended heating rod 10 powers on 3, begins to warm up the water in steel pipe 9；Then the temperature of each temperature point is immediately
It begins to ramp up, as the 6th temperature sensor 16 or the temperature T of the 7th temperature sensor 17>When the maximum temperature of setting, relay 4
It disconnects, single-ended heating rod 10 and the channel of power supply 3 are cut off, and stop heating, as a heat cycles, are realized to sample temperature
The automatic application for spending disturbance carries out multiple heat cycles and constitutes non-steady state heat transfer process.
In order to control the time of one cycle, and improve the discrimination of each temperature point temperature value, setting it is minimum
Temperature is that room temperature adds 5 DEG C, and maximum temperature is that room temperature adds 12 DEG C.
Example IV
The non-steady state process of thermal coefficient of program Inversion Calculation material can also be write, it is specific as shown in figure 4, using
Fortran finite element program calculation of thermal conductivity.Its major sub programs is as follows：
1) mesh subprograms.Generate One-dimensional heat transfer grid wg.dxf files, temperature wd.dxt files and gridding information
Wg.dxt files.
2) input subprograms.Wg.dxt files are read, node, unit information are read in；Wd.dxt files are read, temperature is read in
Spend information.
3) ask subprograms.Heat transfer matrix is assembled, thermal coefficient always rigid matrix is formed.
4) LSM subprograms.Principle of least square method is introduced, least square normal equation is generated.
5) LSLRG subprograms.The function for calling IMSL numerical analysis programs library, solves above-mentioned least square full rank equation
Group.
6) output subprograms.Output thermal coefficient inverse result test.txt, output check1.txt,
Check2.txt, checks whether equation residual value meets required precision.
Claims (6)
1. a kind of method of non-steady state Determination of conductive coefficients, it is characterised in that this approach includes the following steps：
1) thermal perturbation is applied to sample, obtains the non-steady state heat transfer temperature data of material internal；
2) by temperature data, function of many variables fitting is carried out to non-steady state thermal coefficient using least square finite element method；
3) based on least square finite element method according to the non-steady state mistake of thermal coefficient of the temperature value Inversion Calculation material of each node
Journey obtains non-steady state thermal coefficient；
The step of step 3) is：
3-1) in solid thermal conduction problem, it is assumed that density p and specific heat capacity c are not changed over time, in computational domain at arbitrary point,
Meet the equation of heat conduction
In formula, T is temperature (DEG C)；K is thermal coefficient (kJ/mh DEG C)；T be the time (h), formula (3) be steady state heat transfer without
The finite element differential form of endogenous pyrogen；
3-2) to step 3-1) formula (3) that obtains is discrete using the golden method of gal the Liao Dynasty in spatial domain, use difference side in time-domain
It is as follows to obtain finite element scheme for method：
{ T }=(1- ξ) { Tn}+ξ{Tn-1} (5)
Wherein, [C] is heat capacity matrix；[H] is heat transfer matrix；{ P } temperature load array；{ T } node temperature array；Section
Point temperature is to time-derivative array；N is node number, and ζ is the parameter that value range is [0,1]；
{TnBe n-th of node temperature array, { Tn-1Be (n-1)th node temperature array；Δ t is time variation amount；
3-3) assume that thermal coefficient and the relational expression of time are：
K=α1+α2t+α3t2+...+αntn-1 (7)
Fundamental unknown variables are transformed to each undetermined parameter α of description thermal coefficient by node temperature value1, α2..., αn, therefore formula (4) becomes
It is changed to：
[L] { α }={ f } (8)
Matrix [L] is the temperature matrices by face domain Ω, and { f } is that the set of all temperature variations is write a Chinese character in simplified form, and { α } is description heat conduction
Each undetermined parameter α of coefficient1, α2..., αnSet write a Chinese character in simplified form；
Wherein, matrix [L] and { f } are integrated by the corresponding matrix element of unit, i.e.,
LijFor the i-th row jth column element of matrix [L],For [L] matrix the i-th row jth column element of unit e；
In formula,For temperature term, T is temperature；It is temperature to the derivative of time,For heat source item,To pass in and out heat flux, qi is heat flux,For convective exchange item；fiFor i-th of element of { f } matrix；
In formula, N is weight, and t is time, NiFor i row weights, NkFor k row weights, TkFor node temperature, tj-1For step-length j-1 when
Between；
3-4) formula (8) is with parameter alpha1, α2..., αnFor the super equation group of fundamental unknown variables, sought with least square finite element method
Seek globally optimal solution：
ψ=Wk(Lkjαj-fk)2 (12)
In formula, WkThe weight of k-th of equation, L when being cumulativekjFor the row k jth column element of matrix [L], fkIt is the of { f } matrix
K element, αjFor j-th of element of parameter { α }, αiFor i-th of element of parameter { α }；
Since formula (8) is to include the super equation group of multiple unknown numbers, obtained solution is not unique, because this introduces minimum to formula (8)
Least square method turns to formula (12), and ψ is least square method optimal solution；
Formula (13) is the full rank equation group obtained to each coefficient derivation to be asked, and solving this equation group can ask to obtain undetermined parameter
αiValue, then non-steady state thermal coefficient K can be calculated by formula (7).
2. a kind of method of non-steady state Determination of conductive coefficients according to claim 1, it is characterised in that specific in step 1)
For：Sample and multiple sensors are arranged symmetrically, and apply thermal perturbation for sample, the circulating-heating in certain temperature range allows temperature
Degree is uniformly diffused into sample, then acquires the temperature data under different moments different radii, obtains the non-steady state heat transfer of material internal
Temperature data.
3. a kind of method of non-steady state Determination of conductive coefficients according to claim 1, it is characterised in that this method further includes
Step 4)：It writes programed temperature acquisition to calculate in real time with thermal coefficient, obtains the relationship of non-steady state thermal coefficient and time
Curve graph.
4. a kind of tested in non-steady state thermal coefficient using any one of the claims 1-3 methods for acquiring specimen temperature
The test device of data, it is characterised in that：
It includes Die tester, temperature control system and data collecting system；
In Die tester, circular shaped foam plate (5) is arranged in die trial (6) bottom, and die trial (6) center is inserted with steel pipe (9)；
In temperature control system, the single-ended heating rod (10) being connect with relay (4) is arranged in steel pipe (9) center, relay (4)
It is connect with power supply (3)；
Data collecting system includes multiple temperature sensor groups (8), in multiple temperature sensor groups (8), the first temperature sensor
(11) it is arranged in steel pipe (9) inside, second temperature sensor (12) and third temperature sensor (13) are arranged symmetrically the temperature of composition
Degree sensor group, the 4th temperature sensor (14) and the 5th temperature sensor (15) be arranged symmetrically composition temperature sensor group,
The temperature sensor group that 6th temperature sensor (16) and the 7th temperature sensor (17) are arranged symmetrically composition is arranged in examination respectively
Sample (7) is internal from the inside to the outside along different radii, each temperature sensors of multiple temperature sensor groups (8) respectively with temperature acquisition
Module (2) connects, and temperature collecting module (2) is connect with power supply (3) and computer (1) respectively；
Each temperature sensor and the first temperature sensor (11) of multiple temperature sensor groups (8) are passed through by each lead and are opened in
The mode in the hole of circular shaped foam plate (5) is fixed on circular shaped foam plate (5)；
The bottom seal of steel pipe (9), steel pipe (9) is interior to inject cold water, and the first temperature sensor (11) is for monitoring the interior water temperature of pipe；
Second temperature sensor (12) and third temperature sensor (13) are symmetrically arranged in 30 centimeters of radius of sample (7)；
4th temperature sensor (14) and the 5th temperature sensor (15) are arranged in 43.5 centimeters of radius of sample (7)；
6th temperature sensor (16) and the 7th temperature sensor (17) are symmetrically arranged in 72.5 centimeters of radius of sample (7).
5. a kind of test device using described in claim 4 realizes the method applied automatically disturbed to specimen temperature, special
Sign is that this method is：
The original state of relay (4) is to disconnect when on-test is arranged, and temperature collecting module (2) acquires in test coupon (7)
The temperature of heart distance farthest the 6th temperature sensor (16) or the 7th temperature sensor (17) is T, works as T<The lowest temperature of setting
Relay (4) is closed when spending, and single-ended heating rod (10) powers on (3), begins to warm up the water in steel pipe (9)；Then each temperature
The temperature of measuring point is begun to ramp up immediately, as the 6th temperature sensor (16) or the temperature T of the 7th temperature sensor (17)>Setting
When maximum temperature, relay (4) disconnects, and single-ended heating rod (10) and the channel of power supply (3) are cut off, and stop heating, as
One heat cycles realizes the automatic application disturbed to specimen temperature.
6. the test device realization according to claim 5 using described in claim 4 disturbs specimen temperature automatic
The method of application, it is characterised in that：In order to control the time of one cycle, and improve the differentiation of each temperature point temperature value
Degree, the minimum temperature set add 5 DEG C as room temperature, and maximum temperature is that room temperature adds 12 DEG C.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610012630.XA CN105572161B (en) | 2016-01-08 | 2016-01-08 | A kind of method and used test device of non-steady state Determination of conductive coefficients |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610012630.XA CN105572161B (en) | 2016-01-08 | 2016-01-08 | A kind of method and used test device of non-steady state Determination of conductive coefficients |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105572161A CN105572161A (en) | 2016-05-11 |
| CN105572161B true CN105572161B (en) | 2018-09-11 |
Family
ID=55882534
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610012630.XA Active CN105572161B (en) | 2016-01-08 | 2016-01-08 | A kind of method and used test device of non-steady state Determination of conductive coefficients |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105572161B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109408926B (en) * | 2018-10-12 | 2023-04-07 | 大连理工大学 | Method for solving complex structure multi-dimensional transient nonlinear heat conduction inverse problem |
| CN114544700A (en) * | 2022-02-22 | 2022-05-27 | 三峡大学 | Device and method for testing heat conductivity coefficient of mass concrete |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101556256B (en) * | 2009-05-11 | 2012-06-27 | 天津科技大学 | Dual-flat plate thermal conductivity coefficient measuring instrument of thermal insulation materials |
| CN104040327A (en) * | 2011-12-23 | 2014-09-10 | 西格里碳素欧洲公司 | Method for measuring thermal conductivity |
| CN103940845B (en) * | 2014-03-11 | 2017-01-18 | 工业和信息化部电子第五研究所 | Measuring method for thermal conductivity of metal at high temperature |
| CN104237301B (en) * | 2014-09-18 | 2017-05-24 | 毕文明 | In-situ thermal response testing method for layered rock and soil thermophysical properties |
| CN104597078A (en) * | 2015-01-14 | 2015-05-06 | 北京科技大学 | Method for measuring anisotropic material heat conductivity based on small-plane heat source |
-
2016
- 2016-01-08 CN CN201610012630.XA patent/CN105572161B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105572161A (en) | 2016-05-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101126729B (en) | Double heat flux gauge steady state method for measuring material heat conductivity | |
| Heitman et al. | A test of coupled soil heat and water transfer prediction under transient boundary temperatures | |
| CN103954648B (en) | A kind of hemisphere face thermal source steady state method heat conductivity measuring device and method | |
| CN114068051B (en) | Method for calculating temperature flow of nuclear reactor main pipeline coolant based on ultrasonic array | |
| CN104864984B (en) | Microresponse device thermometry based on neutral net | |
| CN104215660B (en) | A kind of method and system simultaneously can surveying solid material heat conductivity and thermal diffusivity | |
| Pătrăşcioiu et al. | Prediction of the outlet temperatures in triple concentric—tube heat exchangers in laminar flow regime: case study | |
| CN107421980B (en) | Heating impedance compensation type thermo-resistance measurement method | |
| CN105572161B (en) | A kind of method and used test device of non-steady state Determination of conductive coefficients | |
| CN106248725A (en) | A kind of porous media Equivalent Thermal Conductivities measuring method | |
| CN113990539B (en) | Ultrasonic measurement reconstruction method for temperature distribution of coolant in hot section of main pipeline of nuclear reactor | |
| CN111795993A (en) | Rock thermophysical property transient test system and method under high temperature and high pressure | |
| Yuan et al. | Multiscale and multilayer structural modeling and simulation on drying of grain packing porous media | |
| CN100585391C (en) | A kind of soil thermal parameter automatic measuring equipment | |
| CN114813828B (en) | Micro-thermal test method for determining thermophysical property parameters of aquifer | |
| Hedayati-Dezfooli et al. | A design of experimental apparatus for studying coupled heat and moisture transfer in soils at high-temperature conditions | |
| CN205280627U (en) | Test testing arrangement who is used for gathering sample temperature data among non - permanent attitude coefficient of heat conductivity | |
| Bhowmik et al. | Analyses of natural convection heat transfer from a heated cylinder mounted in vertical duct | |
| Cortellessa et al. | Experimental and numerical analysis in heat flow sensors calibration | |
| Arvanitidis et al. | A model based study of the drying and shrinkage behavior of a ceramic green body | |
| CN204086184U (en) | A kind of system simultaneously can surveying solid material heat conductivity and thermal diffusivity | |
| CN114544700A (en) | Device and method for testing heat conductivity coefficient of mass concrete | |
| Cheng et al. | A double cell triaxial apparatus for testing unsaturated soil under heating and cooling | |
| Behrem et al. | Estimate of heat transfer coefficient during quenching steel in water | |
| CN104792817A (en) | Measuring system for heat storage coefficient of material |
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 |