CN101476753A - Construction thermal environment and construction energy-saving control method in heating construction - Google Patents

Construction thermal environment and construction energy-saving control method in heating construction Download PDF

Info

Publication number
CN101476753A
CN101476753A CN 200810246324 CN200810246324A CN101476753A CN 101476753 A CN101476753 A CN 101476753A CN 200810246324 CN200810246324 CN 200810246324 CN 200810246324 A CN200810246324 A CN 200810246324A CN 101476753 A CN101476753 A CN 101476753A
Authority
CN
China
Prior art keywords
control index
wall
saving
building
heat preservation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN 200810246324
Other languages
Chinese (zh)
Other versions
CN101476753B (en
Inventor
黄衡栋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
XIAOGAN COLLEGE
Original Assignee
XIAOGAN COLLEGE
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by XIAOGAN COLLEGE filed Critical XIAOGAN COLLEGE
Priority to CN 200810246324 priority Critical patent/CN101476753B/en
Publication of CN101476753A publication Critical patent/CN101476753A/en
Application granted granted Critical
Publication of CN101476753B publication Critical patent/CN101476753B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Building Environments (AREA)

Abstract

The invention relates to a method for controlling building thermal environment and building energy conservation in a heating building, which comprises the following steps: calculating out the comprehensive temperature average value and the corresponding wave amplitude value thereof of a horizontal external surface and the west wall external surface of a local building roof according to data information that the temperature average value t0 of the outdoor air in a local heating building in continuous 5 days in the lately 10 years is less than or equal to 5 DEG C; determining the thermal transmission resistance control index of a roof and a wall in the process of heat preservation in the heating building; determining the heat preservation and energy conservation three-element control index of the wall, the roof and a window in the heating building; and controlling the heat preservation and energy conservation performances of the wall, the roof and the window according to the heat preservation and energy conservation three-element control index. The method gives full play to the use value of the existing meteorological data, adopts a simple method to make the computation of the control indexes of each energy conservation aspect simple and fast, and has small investment cost. The method applied to the designs of the building thermal environment and the building energy conservation can greatly improve the energy conservation efficiency of the building.

Description

Building thermal environments and construction energy-conserving control method in the heating building
Technical field
The present invention relates to building thermal environments and construction energy-conserving control method in the heating building.
Background technology
" 21 century course " passed through in the United Nations's " environmental development " in 1992 conference; China makes an immediate response and has proposed China's Agenda 21; " energy savings; reduce coal consumption; the exploitation renewable sources of energy " fundamental state policy and vital task as China's " the protection environment is adhered to sustainable development ".
In various countries' total energy consumption, architectural energy consumption accounts for 30%~40%, and China is no exception.Simultaneously, China's production of energy total amount annual average rate of increase only is 1/2 of a national total output value annual average rate of increase.Production of energy far lags behind development and national economy.China's architectural exterior-protecting construction (roof, wall and window) thermal property is not high in addition, and heating system efficient is low, causes China's residential architecture heating energy consumption to be three times more than of developed country.Therefore, improve indoor thermal environment, improve quality of residents'life, save energy for building, alleviate national energy and produce short pressure and have urgency.
In China Ministry of Construction " building energy conservation the ninth five-year plan and rule in 2010 ", claim: newly-built heating residential architecture, till the end of the year 1996, on 1980~1981 years local universal design energy consumption level bases, generally reduce by 30%, as the phase I; 1997~2004 years energy-conservation again 30% on energy-conservation basis of phase I, as second stage; Rose in 2005, and required on the energy-conservation basis of second stage, energy-conservation again 30%, as the phase III.Along with China's building energy conservation is carried out and is implemented, country issued " code for thermal design of civil buildings " (GB-50176-93) and " hot summer and warm winter region energy Saving Design of Residential Buildings standard " (JGJ-25-2003).
" code for thermal design of civil buildings " is with China architecture science research institute Hu phosphorus, the Jiang Jian of Zhejiang University is bright, the Chen Qi of Chongqing Univ. of Architecture professor (Architecural Physics doctor tutor) is older generation's building technology person of representative, the study Soviet Union, in conjunction with China's national situation, solved China Xia Redong cold-peace hot summer and warm winter region, under architecture indoor gravity-flow ventilation in the summer condition, the heat insulation control design of indoor thermal environment and building enclosure and China's severe cold, cold and Xia Redong cryogenic region, under the winter architecture indoor heating condition, indoor thermal environment and building enclosure preservation and controlling PRACTICE OF DESIGN problem.Wherein " the heat insulation control of building enclosure under the gravity-flow ventilation condition " is based on multilayer planomural periodic instability heat transfer theory, adopt outdoor climate characteristic parameter statistical value averaging method, temperature profile parameter inequality is as indoor thermal environment controlled condition below having proposed under gravity-flow ventilation (open the door, window) condition:
θ i·max≤t o·max [1.1-1]
In the formula, θ ImaxUnder-indoor natural ventilating (open the door, the window) condition, the maximum of building enclosure internal surface temperature (characteristic parameter) [℃]:
t Omax-local outside air temperature (characteristic parameter) t oMaximum [℃].
By multilayer planomural unsteady heat transfer theory, following formula becomes:
θ i · max ≤ θ ‾ i + A tsa v 0 + A t 0 v i ≤ t o · max - - - [ 1.1 - 2 ]
In the formula, θ iThe mean temperature of-building enclosure inner surface [℃];
A Tsa-building enclosure outer surface (consideration solar radiation) integrated temperature wave amplitude [℃];
v o-outdoor integrated temperature heat wave is reached the attenuation multiple of building enclosure inner surface by outer surface;
A ToThe wave amplitude of-outside air temperature harmonic wave [℃];
v i-outside air temperature harmonic wave is by the indoor attenuation multiple that reaches the building enclosure inner surface.
Obvious formula [1.1-1], the Thermal Environment Control condition under promptly indoor no air-conditioning (gravity-flow ventilation) condition.It can only make within the scope of " can restrain oneself thermal balance " of indoor " people's body heat balance " as close as possible comfortable balance.Also can only make this expedient regulation, with the overheating conditions of control indoor thermal environment.
By formula [1.1-2], at indoor natural ventilating, under the condition of windowing of opening the door, building enclosure internal surface temperature maximum θ Imax, be at internal surface temperature mean value θ iOn the basis, get two-way heat wave effect, promptly outdoor integrated temperature heat wave acts on the building enclosure outer surface and the outside air temperature heat wave acts on the building enclosure inner surface, and decays to the synthetic result of two ripples on the inner surface respectively.According to (Soviet Union) AM ш к л о в e р multilayer planomural unsteady heat transfer formula, two wavelength-division supplementary biographies are to the attenuation multiple v of inner surface oAnd v iDetermine by following formula:
v 0 = 0.9 e ΣD / 2 s 1 + α i s 1 + y 1,0 s 2 + y 1,0 s 2 + y 2,0 · · · s n + y n - 1,0 s n + y n , 0 α 0 + y n , 0 α 0 - - - [ 1.1 - 3 ]
v i = 0.95 α i + y i · f α i - - - [ 1.1 - 4 ]
In the formula: the thermal characteristics parameter-heat inertia index of D-multilayer planomural
s 1, s 2..., s n-from inside to outside, the heat storage coefficient [w/m of layers of material 2K] (air space heat storage coefficient S=0).
y 1,0, y 2,0, y N, 0-layers of material outer surface heat storage coefficient [w/m from inside to outside 2K].
α i, α o-be respectively the coefficient of heat transfer [w/m on the inside and outside surface of planomural 2K].
y If-planomural inner surface heat storage coefficient [w/m 2K].
Under the multilayer planomural condition of limited thickness, material layer is subjected to the temperature harmonic wave of cyclic fluctuation and does the time spent, and the temperature fluctuation size that it is surperficial is not only relevant with the physical property of material itself, and relevant with boundary condition.Promptly along the temperature wave direction of advance, its surface temperature degree of fluctuation had influence with the thermal property and the radiating condition of the contacted medium of this material layer (another material layer or air).
Know that by heat transfer theory when material surface was subjected to the heat effect of temperature harmonic wave, the size of its surface temperature fluctuation amplitude depended primarily on the material thermophysical property in " big ups and downs layer " scope." the physical computing thickness " that this " big ups and downs layer " scope is exactly this layer-be thermal inertia index D=1.0.So, to the limited thickness material layer, when " big ups and downs layer " also do not exceed this material layer scope, i.e. the heat inertia index D of material layer 〉=1.0 o'clock, then its surface temperature fluctuation is main relevant with the thermophysical property of material layer, at this moment, can think y=s approx; When D<1.0 of material layer, i.e. " big ups and downs layer " when having exceeded this material ranges, the heat storage coefficient y on material layer back of the body corrugated then, fluctuation has influence to surface temperature as boundary condition, this moment y ≠ s.Thus, for two kinds of situations of material layer " physical computing thickness " D 〉=1.0 and D<1.0, determine the material surface heat storage coefficient by following recurrence formula respectively:
[1], material layer outer surface heat storage coefficient calculates: from inboard ground floor, promptly inversion degree wavefront advances direction, outwards successively calculates successively.
The 1st layer of outer surface heat storage coefficient:
Work as D 1〉=1.0 o'clock, y 1,0=s 1[1.1-5a]
Work as D 1<1.0 o'clock, y 1,0 = R 1 S 1 2 + α i 1 + R 1 α i - - - [ 1.1 - 5 b ]
In the formula, α i-be the influence of the border factor (with the degree of room air heat exchange) of inboard the 1st layer of planomural, R 1, S 1All are influences of this layer material thermal characteristic parameter of reflection.
From the second layer, later arbitrary layer outer surface heat storage coefficient:
Work as D m〉=1.0 o'clock, y 1,0=s m[1.1-6a]
Work as D m<1.0 o'clock, y m , 0 = R m S m 2 + y m - 1,0 1 + R m y m - 1,0 , ( m = 2,3,4 . . . n ) - - - [ 1.1 - 6 b ]
Y in the formula M-1,0Represent of the influence of the thermal storage effect of m-1 layer to the m layer.
The heat storage coefficient of outermost external surface promptly is the heat storage coefficient of planomural outer surface, that is:
y n,0=y 0,f [1.1-6c]
[2], planomural inner surface heat storage coefficient calculates: determine planomural inner surface heat storage coefficient y IfThe time, divide following several situation to carry out.
Work as D 1〉=1.0 o'clock, y I, f=s 1[11-7]
Work as D 1<1.0, but D 1+ D 2〉=1.0 o'clock, y i , f = R 1 S 1 2 + S 2 1 + R 1 S 1 - - - [ 1.1 - 8 ]
Work as D 1+ D 2+ ... D M-1<1.0, and D 1+ D 2+ ... D m, then ask it at 〉=1.0 o'clock by following recurrence formula:
y m - 1 , i = R m - 1 S m - 1 2 + S m 1 + R m - 1 S m y m - 2 , i = R m - 2 S m - 2 2 + y m - 1 , i 1 + R m - 2 y m - 1 , i y 1 , i = y i , f = R 1 S 1 2 + y 2 , i 1 + R 1 y 2 , i - - - [ 1.1 - 9 ]
When the ∑ D of whole planomural<1.0, then at first obtain the inner surface heat storage coefficient of last one deck:
y n · i = R n S n 2 + α o 1 + R n α o - - - [ 1.1 - 10 ]
In the formula, R nResistance of heat transfer [the m of-planomural outermost one deck 2K/w]
S nHeat storage coefficient [the m of-outermost layer of material 2K/w]
α oThe coefficient of heat transfer [the w/m of-planomural outer surface 2K] (by " standard " α 0=23.0)
Try to achieve y NiAfter, obtain the heat storage coefficient y of planomural inner surface again by formula [1.1-9] recursion If
From the above mentioned, want to obtain θ Imax, then must be by the top various complicated calculation of carrying out.While at first will be known the material characteristic parameter value that " structure " of building enclosure multilayer planomural had: the heat storage coefficient S of each material layer, resistance of heat transfer R and heat inertia index D etc.Therefore, first supposition multilayer planomural " structure ", i.e. the S of each material layer, R, D is known, the various θ that obtains at last above the substitution ImaxAfter, substitution formula [1.1-1] is judged.If formula [1.1-1] is false, need suppose " structure " scheme again, double counting again is till formula [1.1-1] is set up fully.This computational process is called " tentative calculation proof method repeatedly " (supposition of " structure " scheme need be carried out by rule of thumb).
Hence one can see that, the heat insulation control method of this building enclosure, in theory, science, complete, rigorous.But not enough below existing:
[1], computational process is numerous and diverse, and must preestablish tentative calculation checking repeatedly behind the characterisitic parameter of building enclosure (roof, wall), can not directly obtain its characteristic parameter value and control the indoor thermal environment situation;
[2], indoor thermal environment control is inapplicable at indoor thermal environment and the Energy Saving Control under room conditioning (window of the closing the door) condition set up under gravity-flow ventilation (open the door and the window) condition.Therefore must inherit above-mentioned advantage, the more new development of its deficiency.
In sum, " code for thermal design of civil buildings " is (GB50176-93) only suitable, and China Xia Redong was cold at that time, the preservation and controlling of building enclosure under the heat insulation control of the building enclosure (roof, wall) under hot summer and warm winter region gravity-flow ventilation in the summer condition and three Norths cold, severe cold and Xia Redong cryogenic region indoor heating in the winter condition does not relate to for the Energy Saving Control milli; Both fail to unite on heat transfer theory simultaneously; After needing to suppose the building enclosure characteristic parameter value earlier on the computational methods, just can carry out experimental checking computations.
" hot summer and warm winter region energy Saving Design of Residential Buildings standard " be " code for thermal design of civil buildings " (GB50176-93) in, on the heat insulation control of the building enclosure basis, China south of foundation is than under low latitudes coastal " hot summer and warm winter region " room conditioning in the summer condition " energy conservation standard " under indoor natural ventilating in summer (open the door and the window) condition.It emphasizes that aggregate measures such as indoor natural ventilating and sunshade solve the Energy Saving Control of the building enclosure in summer.Therefore, it has under the indoor natural ventilating condition in above-mentioned " standard " middle summer, the deficiency in the heat insulation control of building enclosure.
According to top brief introduction, for making each department " design standard for energy efficiency of buildings " have comparativity, indoor thermal environment control under summer room conditioning and the winter indoor heating condition and be converted to the heat insulation and preservation and controlling of building enclosure and the correlation between the Energy Saving Control, promptly from then on the controlled condition of three aspects and control parameter and local climate feature are controlled correlation between parameter, building enclosure (window, wall) correlation between the Characteristics Control parameter is set out, find out the inner link between their correlations, unifying, science, under the complete theoretical direction, create science, the shortcut Calculation method, unified formulation each department with Chinese characteristics " design standard for energy efficiency of buildings " are current urgency major issues to be solved.
By above-mentioned prior art brief introduction as can be known: thermal environment under " southern energy conservation standard " under " indoor natural ventilating " condition and the thermal environment under " indoor natural ventilating " condition and heat insulation " standard " and " indoor heating condition " and insulation " standard " etc., have the disunity on the heat transfer theory: for example, " the southern energy conservation standard " under the thermal environment under " indoor natural ventilating " condition and heat insulation " standard " and " the indoor natural ventilating condition " adopts building enclosure multilayer planomural periodic instability heat transfer theory; Secondly, also there is the disunity on the computational methods in they: for example, the thermal environment under " indoor natural ventilating " condition and heat insulation " standard " and " southern energy conservation standard ", both have all adopted " tentative calculation proof method repeatedly "; These computational methods all belong to empirical assignment method.Therefore, cause building enclosure (roof, wall) the energy conservation characteristic control index in above-mentioned three " energy conservation standards " to lose " corresponding comparativity ".
In addition, above-mentioned several " energy conservation standard " all with corresponding conditions under " standard " lost the correlation contact, promptly fail corresponding building enclosure Energy Saving Control in " energy conservation standard ", get in touch with the heat insulation and preservation and controlling three of thermal environment in " standard " and building enclosure (roof and wall) and unite.Simultaneously, " energy conservation standard " fails the correlation contact between the Characteristics Control index of the Characteristics Control index of building enclosure (roof and wall) and building enclosure (window), Characteristics Control index by wall, directly calculate the Characteristics Control index of locking window, thereby caused the artificial assignment arbitrarily of Characteristics Control index of window.Make three " energy conservation standards " middle corresponding Characteristics Control index of building enclosure window also lose " corresponding comparativity ".
Summary of the invention
Problem to be solved by this invention is, building thermal environments and construction energy-conserving control method in a kind of heating building that can solve heat preservation energy-saving problem in the heating building preferably is provided.
Technical scheme provided by the invention is:
Building thermal environments and construction energy-conserving control method in the heating building may further comprise the steps:
One, according to the local nearest 10 years continuous outside air temperature mean value t more than 5 days of heating building o≤ 5 ℃ outside air temperature maximum t OmaxWith minimum of a value t OminAnd changing value t 0-l, the local building horizontal outer surface in roof and north wall outer surface solar irradiance daily mean I H, I N, maximum I H max, I N maxCalculate local building roof water
The integrated temperature mean value t of flat outer surface and north wall outer surface H-saAnd t N-saWave amplitude A corresponding with it H-tsaAnd A N-tsa
(1), calculates continuous outside air temperature mean value t more than 5 days in local nearest 10 years by the measured value of local weather bureau o, amplitude A 0And initial phase
Figure A200810246324D00091
t o≤ 5 ℃,
[1], outside air temperature mean value: t 0
[2], outside air temperature wave amplitude:
A 0 = 1 2 ( t 0 . max - t 0 . min )
[3], t 0.maxPairing time τ 0By changing value t O-τDetermine by interpolation method:
τ o = t 014 + t 0 . max - t 014 t 014 - t 020 × 6
In the formula: t 014, t 020Be respectively continuous outside air temperature mean value t more than 5 days in nearest 10 years oThe mean value of the outside air temperature when≤5 ℃ 14 and 20;
With initial time τ 0Corresponding initial phase
Figure A200810246324D00094
For:
(2) calculate t H-saAnd t N-saAnd A H-tsaAnd A N-tsa
[1], calculates t H-saAnd t N-sa:
t ‾ H - sa = t ‾ 0 + ρ I ‾ H α 0 , t ‾ N - sa = t ‾ 0 + ρ I ‾ N α 0
ρ is the absorption coefficient of horizontal outer surface in roof or north wall outer surface material, ρ=0.7;
α 0Be the coefficient of heat transfer of horizontal outer surface in roof or north wall outer surface, α 0=23.0[w/m 2K];
[2] local building horizontal outer surface in roof and north wall outer surface solar radiation equivalent temperature wave amplitude:
A IH · max = ρ I H max - I ‾ H α 0 , A IN · max = ρ I N max - I ‾ N α 0
[3], calculate A H-tsaAnd A N-tsa: by outside air temperature wave amplitude A 0And phase angle With solar radiation equivalent temperature wave amplitude and phase angle
Figure A200810246324D000910
With
Figure A200810246324D000911
Calculate by following vector superposition formula:
Figure A200810246324D00101
The resistance of heat transfer control index of roof and wall when two, determining to be incubated in the heating building
The resistance of heat transfer of wall control index [R when [1], being incubated o] N min: determine by locality building the eighties universal design wall heat transfer resistance;
The resistance of heat transfer on roof control index [R when [2], being incubated o] H min:
[ R o ] H min = 18 - t &OverBar; H - sa 18 - &theta; &OverBar; i R i ; (t d-1)<θ i<18 R i=0.11[m 2k/w];
t dFor indoor when adopting 18 ℃ of dim temperature, the dew point under the heating building average relative humidity in local winter;
Three, determine in the heating building not heat preservation energy-saving ternary control index: the resistance of heat transfer of wall control index [R during heat preservation energy-saving with the wall of window O-E] N min, thermal inertia control index [D 0-E] N min, heat flow density control index [q 0-E] N max
[ R o - E ] N min = [ R o ] N min 1 - &epsiv; ; [ D 0 - E ] N min = 2.13 ln ( 1.46 A N - tsa R i [ R 0 - E ] N min ) , R i = 0.11 [ m 2 k / w ] ; [ q 0 - E ] N max = 18 - t &OverBar; N - sa [ R o ] N min ( 1 - &epsiv; )
Be the heat preservation energy-saving efficient of wall for ε in the formula,
Four, the heat preservation energy-saving ternary on roof control index: the resistance of heat transfer on roof control index [R during heat preservation energy-saving O-E] H min, thermal inertia control index [D 0-E] H min, heat flow density control index [q 0-E] H max
[ R o - E ] H min = [ R o ] H min 1 - &epsiv; ; [ D o ] H min = 2.13 ln ( 1.46 A H - tsa R i [ R o - E ] H min ) , R i = 0.11 [ m 2 k / w ] ; [ q o - E ] H max = 18 - t &OverBar; H - sa [ R o ] H min ( 1 - &epsiv; ) ;
Five, determine that whole face is the heat preservation energy-saving ternary control index of window in the heating building
The heat preservation energy-saving resistance of heat transfer control index [R of window S-E] Min:
[ R s - E ] N min = 1 2 [ ( &xi; - s ) + ( &xi; - s ) 2 + 4 s ] [ R o - E ] N - min
In the formula, s=F s/ F o, F sWindow ara, F oWall area, ξ=1/ (1+s);
The heat preservation energy-saving thermal inertia control index [D of window S-E] N min:
[ D s - E ] N min = 2.13 ln ( 1.46 A N - tsa R i [ R s - E ] N min )
The heat preservation energy-saving heat flow density control index [q of window S-E] N max:
[ q s - E ] N max = 18 - t &OverBar; N - sa [ R s - E ] N min
Six, determine in the heating building heat preservation energy-saving ternary control index of wall and window in the wall of band window
The heat preservation energy-saving resistance of heat transfer of wall control index according to
Figure A200810246324D00113
Determine, the heat preservation energy-saving thermal inertia control index of wall according to [ D 0 - E ] N min = 2.13 ln ( 1.46 A N - tsa R i [ R 0 - E ] N min / &alpha; ) , The heat preservation energy-saving heat flow density of wall control index according to [ q 0 - E ] N max = 18 - t &OverBar; N - sa [ R o ] N min ( 1 - &epsiv; ) &times; &alpha; Determine;
The heat preservation energy-saving resistance of heat transfer of window control index according to
Figure A200810246324D00116
Determine; The heat preservation energy-saving thermal inertia of window control index according to [ D s - E ] N min = 2.13 ln ( 1.46 A N - tsa R i [ R s - E ] N min / &alpha; ) , The heat preservation energy-saving heat flow density of window control index according to [ q s - E ] N max = 18 - t &OverBar; N - sa [ R s - E ] N min / &alpha; Determine;
Wherein &alpha; = [ R s - E ] min &CenterDot; [ R o - E ] W min ( 1 + s ) [ R s - E ] min + s [ R o - E ] W min [ R o - E ] W min
Control the heat preservation energy-saving performance of wall, roof and window according to the heat preservation energy-saving ternary control index (resistance of heat transfer control index, heat inertia index control index and heat flow density control index) of above-mentioned wall, roof and window.
ε in the above-mentioned steps three and four=heating building energy-saving efficiency * 70%.The heating building energy-saving efficiency can be by the pertinent regulations value in the national standard.Be not less than 50% as CNS regulation heating building energy-saving efficiency, then ε=0.350 in the step 3 and four; For the heating building energy-saving efficiency is 65%, then ε=0.455 in the step 3 and four
The present invention's's " relevant control method of building thermal environments " advantage with energy-saving design in construction:
According to the front building thermal environments and building energy conservation are controlled the argumentation and the demonstration of each aspect control problem, obviously the present invention has following advantage:
[1] given full play to the use value of China's meteorological data that accumulate decades, especially adopt " the assembly average method of outdoor climate feature control parameter " (seeing step 1), make the calculating of the control index of each aspect of building energy conservation, simple and fast, input cost is little.
[2] exempted the relevant area ratio of window to wall of known " typical building " computation model, bodily form coefficient, parameter limit such as building enclosure heat inertia index have been eliminated all artificial randomness, simultaneously, make the architect bring into play initiative in the architectural design creation.
[3] unified indoor thermal environment controlled condition under the indoor heating condition in winter at Gai Jin And on " code for thermal design of civil buildings " basis, made Thermal Environment Control condition unification under the indoor heating condition in winter in summer to the multilayer planomural periodic instability heat transfer theory basis.
[4] " multilayer planomural periodic instability heat transfer formula of reduction " (thermal inertia control index computing formula) especially innovation proposed, make indoor thermal environment controlled condition, change over " directly calculating lock method " by original " tentative calculation proof method repeatedly ", exempted complicated calculation, simplified computational process widely, be convenient to the architect and grasp.
[5] on indoor heating Thermal Environment Control basis, to the limit value (promptly controlling index) of the insulation thermal technology Characteristics Control parameter on the preservation and controlling aspect on building enclosure (wall, roof), " directly calculating locking ", clear concept calculates simple and direct.
[6] because the difference of the area ratio of window to wall of fenestrate building enclosure, cause its Chuan Qiang And to join equivalent resistance of heat transfer changes thereupon, make the overall Energy Saving Control of fenestrate building enclosure can not reach the control requirement of the energy-saving efficiency that building enclosure should have fully, can only reach energy-saving efficiency and control within desired 88%~100% scope and change.Wait true, accurate, the reliability of statistics building energy conservation efficient for national Chuan Qiang And connection, " imitating resistance of heat transfer form invariance principle " (step 6) then proposed, by energy-saving efficiency percent value in the determined above-mentioned scope of the various area ratio of window to wall of building, with the resistance of heat transfer that removes window and wall, to adjust the resistance of heat transfer control index of window and wall, make fenestrate building enclosure can reach the requirement of the energy-saving efficiency 100% of national regulation fully.
The specific embodiment
2-2 specifically calculates
1.1, local winter meteorological data calculating
1. press the nearest 10 years the coldest month air themperature t in Hankow by Hubei Province weather bureau 0≤ 5[℃] measured value (as following table 01) counting chamber outer air temperature-averaging value
Figure A200810246324D0012162137QIETU
, wave amplitude A 0With the first phase parallactic angle
Figure A200810246324D00121
(initial time τ 0)
Table 01 Hankow t 0≤ 5[℃] (winter) outside air temperature measured value
Figure A200810246324D00122
[1] outdoor air air mean value by table 01 value is:
Figure A200810246324D00123
[2] outside air temperature wave amplitude:
[3] the pairing initial time of outside air temperature maximum is defined as with the Nei Chafa formula by table 01:
&tau; o = 14 + 6.2802 - 5.0083 5.0082 - 3.3220 &times; 6 = 18.5255 [ h ]
Corresponding first phase parallactic angle:
Figure A200810246324D00131
When the time origin of coordinates was taken at high noon 12, the outside air temperature harmonic wave was accurate to the instantaneous value of first-harmonic:
t 0-τ=2.78+3.07cos[15(τ-12.53)]
Or t 0-τ=2.78+3.07cos (15 τ-187.88)
2. press Hankow t in the coldest month of nearest 10 years of winter by Hubei Province weather bureau 0≤ 5[℃] solar irradiance observation (as table 02).
The table 02. Hankow per day total amount of solar radiation in the coldest month of winter, daily mean, day maximum and time of occurrence thereof
Figure A200810246324D00132
[1] solar radiation equivalent temperature mean value I[w/m 2] calculate (result is shown in table 03) by table 02 value;
[2] outdoor integrated temperature mean value is by formula Calculate by (result is shown in table 03);
[3] solar radiation equivalent temperature wave amplitude (by " people have the building thermal technique design specification ") is by formula A I . max = &rho; I max - I &OverBar; &alpha; 0 Calculate by (result is shown in table 03);
[4] outdoor integrated temperature wave amplitude calculates: by outside air temperature wave amplitude A O, the phase angle
Figure A200810246324D00135
With solar radiation equivalent temperature wave amplitude A ImaxAnd phase angle
Figure A200810246324D00136
Calculate (following all parameter calculation procedures and result of calculation are shown in table 03) by following vector superposition formula:
1. make progress integrated temperature wave amplitude and phase angle of level:
Figure A200810246324D00137
Figure A200810246324D00138
2. Dong Qiang is to integrated temperature wave amplitude and phase angle:
Figure A200810246324D00139
Figure A200810246324D001310
3. comprehensive wave amplitude of western metope and phase angle:
Figure A200810246324D001311
Figure A200810246324D00141
4. southern metope integrated temperature wave amplitude and phase angle:
Figure A200810246324D00142
Figure A200810246324D00143
5. north wall surfaces integrated temperature wave amplitude and phase angle %
Figure A200810246324D00145
Mean value, wave amplitude and the phase angle of table 03 solar radiation equivalent temperature, outside air temperature and integrated temperature
Figure A200810246324D00146
1.2, the heat-holding dampproof controlled condition and the insulation ternary control index of building enclosure-roof and wall
[1] indoor heat insulating anti humid heat environment controlled condition
It is indoor that to adopt dim temperature be t i=18[℃] time, indoor heat insulating anti humid heat environment controlled condition:
In 78 cities, Hubei Province, 59 city relative humidity in winter are arranged
Figure A200810246324D00151
Get this value as city (winter) the building enclosure inner surface of the fraction 79% moistureproof Thermal Environment Control of indoor heat insulating of dewfall not.Promptly ought indoorly adopt dim temperature t i=18[℃] time, indoor saturated vapor partial pressure is: E=15.47[mmHg], work as relative humidity
Figure A200810246324D00152
The time, indoor actual partial pressure of water vapor:
Figure A200810246324D00153
As this actual partial pressure of water vapor e, when reaching capacity water vapor pressure, corresponding dew point is: t d=13.5[℃] (or get indoor when adopting 18 ℃ of dim temperature, the dew point t under the heating building average relative humidity in local winter d).Therefore, require internal surface temperature mean value by dewfall controlled condition not:
i]≥t d-1=13.5-1=12.5[℃]
[2] under the indoor heating condition in winter, the determining of the insulation ternary of roof and wall control index: with north wall not dewfall be benchmark, the known universal design eighties: the resistance of heat transfer [R of wall 0] Min=0.5[m 2K/w], check whether dewfall of north wall inner surface:
Figure A200810246324D00154
(inner surface is dewfall not)
Requiring roof and north wall inner surface mean temperature θ IN=15.6[℃] condition under, the roof heat insulation resistance of heat transfer:
[ R 0 &prime; ] min = 18 - 15.5 18 - 15.6 &times; 0.11 = 0.11 [ m 2 k / w ]
The insulation ternary control index of wall:
[ R o &prime; ] min = 0 . 5 [ m 2 k / w ]
[ D o &prime; ] min = 2.13 ln ( 1.46 A N - tsa R i [ R o ] N min ) = 2.13 ln ( 1.46 &times; 0.72 &times; 0.11 0.5 ) = 0
[ q o &prime; ] max = 18 - 7.2 0.5 = 21.6 [ w / m 2 ]
The insulation ternary control index on roof:
[ R 0 &prime; ] min = 0 . 11 [ m 2 k / w ] [ D 0 &prime; ] min = 2.13 ln ( 1.46 A N - tsa &prime; R i [ R 0 &prime; ] min ) = 2.13 ln ( 1.46 &times; 0.75 &times; 0.11 0.11 ) = 5.1 ; [ q 0 &prime; ] max = 18 - 15.5 18 - 15.6 &times; 0.11 = 0.1 [ w / m 2 ]
1.3, the heat preservation energy-saving ternary control index of building enclosure-roof, wall determines:
The heat preservation energy-saving ternary control index condition on wall and roof:
[1] heat preservation energy-saving of wall control index, in the heating building by national regulation energy-conservation 50%, heat preservation energy-saving efficient ε=0.5 * 0.7=0.35 then has:
[ R o - E &prime; ] min = [ R o &prime; ] min 1 - &epsiv; = 0.5 1 - 0.35 = 0 . 77 [ m 2 k / w ] [ D o - E &prime; ] min = 2.13 ln ( 1.46 &times; 0.72 &times; 0.11 0 . 77 ) = 0 [ q o - E &prime; ] max = 18 - 7 . 2 0.77 = 14.0 [ w / m 2 ]
[2] thermal insulation of roof Energy Saving Control index:
[ R o - E &prime; ] min = 0.11 0.65 = 0.17 [ m 2 k / w ]
[ D o - E &prime; ] min = 2.13 ln ( 1.46 &times; 0.75 &times; 0.11 0 . 17 ) = 4.2
[ q o - E &prime; ] max = 18 - 15.5 0.17 = 1.47 [ w / m 2 ]
1.4, the heat preservation energy-saving ternary of window control index:
By " window wall correlation principle ", under known area ratio of window to wall s condition, then have:
Figure A200810246324D00164
In the formula, The resistance of heat transfer of-wall;
Figure A200810246324D00166
-under known area ratio of window to wall s condition, the resistance of heat transfer of window;
Figure A200810246324D00167
Window wall relative coefficient during-corresponding window-wall ratio s.Now getting s=1.0 is example:
When s=1.0, Then the heat preservation energy-saving of window is controlled index:
[ R s - E &prime; ] min = 0.78 [ R o - E &prime; ] min = 0.78 &times; 0.77 = 0.60 [ m 2 k / w ]
[ D s - E &prime; ] min = 2.13 ln ( 1.46 &times; 0.72 &times; 0.11 0 . 60 ) = 0
[ q s - E &prime; ] max = &Delta;t [ R s - E &prime; ] min 18 - 7 . 2 0 . 60 = 18 . 0 [ w / m 2 ]
Thus, the heat preservation energy-saving ternary of above-mentioned building enclosure-wall, roof, window control index is shown in following table 04
The heat preservation energy-saving ternary control index of 04 roof, wall and window.
Figure A200810246324D001612
Annotate: after considering solar irradiance winter, energy-conservation more favourable thermal insulation of roof.
1.5, the adjustment of building enclosure heat preservation energy-saving ternary control index:
By " window wall equivalent resistance of heat transfer form invariance principle in parallel ", then have:
R E = [ R &prime; s - E ] min [ R o - E &prime; ] min ( 1 + s ) [ R s - E &prime; ] min + s [ R o - E &prime; ] min = 0.60 &times; 0.77 ( 1 + 1 ) 0.60 + 1 &times; 0.77 = 0.924 1.37 = 0.67 < 0.77
Get &alpha; = 0.67 0.77 = 0.87 , [ R s - E &prime; ] min = 0.60 0.87 = 0.69 ; [ R o - E &prime; ] min = 0 . 77 0.87 = 0 . 88 ; Then adjust the equivalent resistance of heat transfer [R in back E]/α, promptly
R E &prime; = 0.69 &times; 0.88 &times; ( 1 + 1 ) 0.69 + 1 &times; 0.88 = 0.77 [ m 2 k / w ]
At last, winter indoor adopting under the dim condition, the heat preservation energy-saving ternary of building enclosure control index is shown in following table 05
The insulation of 05 adjusted roof, wall, window-energy-conservation ternary control index
Figure A200810246324D00172

Claims (4)

1. building thermal environments and construction energy-conserving control method in the heating building may further comprise the steps:
One, according to continuous outside air temperature mean value more than 5 days in local nearest 10 years of the heating building
Figure A200810246324C00021
Outside air temperature maximum t O-maxWith minimum of a value t O-mimAnd changing value t 0-τ, the local building horizontal outer surface in roof and north wall outer surface solar irradiance daily mean I H, I N, maximum I Hmax, I NmaxCalculate the integrated temperature mean value t of local building horizontal outer surface in roof and north wall outer surface H-saAnd t N-saWave amplitude A corresponding with it H-tsaAnd A N-tsa
(1), calculates continuous outside air temperature mean value t more than 5 days in local nearest 10 years by the measured value of local weather bureau o, amplitude A 0And initial phase
Figure A200810246324C00022
Figure A200810246324C00023
[1], outside air temperature mean value: t 0
[2], outside air temperature wave amplitude:
A 0 = 1 2 ( t 0 . max - t 0 .min )
[3], t 0.maxPairing time τ 0By changing value t O-τDetermine by interpolation method:
&tau; o = t 014 + t 0 . max - t 014 t 014 - t 020 &times; 6
In the formula: t 014, t 020Be respectively continuous outside air temperature mean value t more than 5 days in nearest 10 years o≤ 5 ℃
The mean value of the outside air temperature when 14 and 20;
With initial time τ 0Corresponding initial phase
Figure A200810246324C00026
For:
Figure A200810246324C00027
(2) calculate t H-saAnd t N-saAnd A H-tsaAnd A N-tsa
[1], calculates t H-saAnd t N-sa:
t &OverBar; H - sa = t &OverBar; 0 + &rho; I &OverBar; H &alpha; 0 , t &OverBar; N - sa = t &OverBar; 0 + &rho; I &OverBar; N &alpha; 0
ρ is the absorption coefficient of horizontal outer surface in roof or north wall outer surface material, ρ=0.7;
α 0Be the coefficient of heat transfer of horizontal outer surface in roof or north wall outer surface, α 0=23.0;
[2] local building horizontal outer surface in roof and north wall outer surface solar radiation equivalent temperature wave amplitude:
A IH . max = &rho; I H max - I &OverBar; H &alpha; 0 , A IN . max = &rho; I N max - I &OverBar; N &alpha; 0
[3], calculate A H-tsaAnd A N-tsa: by outside air temperature wave amplitude A 0And phase angle
Figure A200810246324C0003121928QIETU
, with solar radiation equivalent temperature wave amplitude and phase angle
Figure A200810246324C00033
With
Figure A200810246324C00034
Calculate by following vector superposition formula:
Figure A200810246324C00035
The resistance of heat transfer control index of roof and wall when two, determining to be incubated in the heating building
The resistance of heat transfer of wall control index [R when [1], being incubated o] Nmin: determine by locality building the eighties universal design wall heat transfer resistance;
The resistance of heat transfer on roof control index [R when [2], being incubated o] Hmin:
[ R o ] H min = 18 - t &OverBar; H - sa 18 - &theta; &OverBar; i R i ; (t d-1)<θ i<18 R i=0.11;
t dFor indoor when adopting 18 ℃ of dim temperature, the dew point under the heating building average relative humidity in local winter;
Three, determine in the heating building not heat preservation energy-saving ternary control index: the resistance of heat transfer of wall control index [R during heat preservation energy-saving with the wall of window O-E] Nmin, thermal inertia control index [D 0-E] Nmin, heat flow density control index [q 0-E] Nmax
[ R o - E ] N min = [ R o ] N min 1 - &epsiv; ; [ D 0 - E ] N min = 2.13 ln ( 1.46 A N - tsa R i [ R 0 - E ] N min ) , R i = 0.11 ; [ q 0 - E ] N max = 18 - t &OverBar; N - sa [ R o ] N min ( 1 - &epsiv; )
Be the heat preservation energy-saving efficient of wall for ε in the formula,
Four, the heat preservation energy-saving ternary on roof control index: the resistance of heat transfer on roof control index [R during heat preservation energy-saving O-E] Hmin, thermal inertia control index [D 0-E] Hmin, heat flow density control index [q 0-E] Hmax
[ R o - E ] H min = [ R o ] H min 1 - &epsiv; ; [ D o ] H min = 2.13 ln ( 1.46 A H - tsa R i [ R o - E ] H min ) , R i = 0.11 ; [ q o - E ] H max = 18 - t &OverBar; H - sa [ R o ] H min ( 1 - &epsiv; ) ;
Five, determine that whole face is the heat preservation energy-saving ternary control index of window in the heating building
The heat preservation energy-saving resistance of heat transfer control index [R of window S-E] Min:
[ R s - E ] N min = 1 2 [ ( &xi; - s ) + ( &xi; - s ) 2 + 4 s ] [ R o - E ] N - min
In the formula, s=F s/ F o, F sWindow ara, F oWall area, ξ=1/ (1+s);
The heat preservation energy-saving thermal inertia control index [D of window S-E] Nmin:
[ D s - E ] N min = 2.13 ln ( 1.46 A N - tsa R i [ R s - E ] N min )
The heat preservation energy-saving heat flow density control index [q of window S-E] N Max:
[ q s - E ] N max = 18 - t &OverBar; N - sa [ R s - E ] N min
Six, determine in the heating building heat preservation energy-saving ternary control index of wall and window in the wall of band window
The heat preservation energy-saving resistance of heat transfer of wall control index according to
Figure A200810246324C00044
Determine, the heat preservation energy-saving thermal inertia control index of wall according to [ D 0 - E ] N min = 2.13 ln ( 1.46 A N - tsa R i [ R 0 - E ] N min / &alpha; ) , The heat preservation energy-saving heat flow density of wall control index according to [ q 0 - E ] N max = 18 - t &OverBar; N - sa [ R o ] N min ( 1 - &epsiv; ) &times; &alpha; Determine;
The heat preservation energy-saving resistance of heat transfer of window control index according to
Figure A200810246324C00047
Determine; The heat preservation energy-saving thermal inertia of window control index according to [ D s - E ] N min = 2.13 ln ( 1.46 A N - tsa R i [ R s - E ] N min / &alpha; ) , The heat preservation energy-saving heat flow density of window control index according to [ q s - E ] N max = 18 - t &OverBar; N - sa [ R s - E ] N min / &alpha; Determine;
Wherein &alpha; = [ R s - E ] min &CenterDot; [ R o - E ] W min ( 1 + s ) [ R s - E ] min + s [ R o - E ] W min [ R o - E ] W min
Control the heat preservation energy-saving performance of wall, roof and window according to the heat preservation energy-saving ternary control index of above-mentioned wall, roof and window.
2. method according to claim 1 is characterized in that: ε in the step 3 and four=heating building energy-saving efficiency * 70%.
3. method according to claim 2 is characterized in that: ε in the step 3 and four=0.350.
4. method according to claim 2 is characterized in that: ε in the step 3 and four=0.455.
CN 200810246324 2008-12-31 2008-12-31 Construction thermal environment and construction energy-saving control method in heating building Expired - Fee Related CN101476753B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 200810246324 CN101476753B (en) 2008-12-31 2008-12-31 Construction thermal environment and construction energy-saving control method in heating building

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 200810246324 CN101476753B (en) 2008-12-31 2008-12-31 Construction thermal environment and construction energy-saving control method in heating building

Publications (2)

Publication Number Publication Date
CN101476753A true CN101476753A (en) 2009-07-08
CN101476753B CN101476753B (en) 2011-05-11

Family

ID=40837501

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 200810246324 Expired - Fee Related CN101476753B (en) 2008-12-31 2008-12-31 Construction thermal environment and construction energy-saving control method in heating building

Country Status (1)

Country Link
CN (1) CN101476753B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101894183A (en) * 2010-07-02 2010-11-24 重庆星能建筑节能技术发展有限公司 Method for analyzing and designing green building
CN102095591A (en) * 2010-12-01 2011-06-15 西安建筑科技大学 Energy-saving performance estimating method of driven solar heating building
CN103033534A (en) * 2012-12-20 2013-04-10 深圳市建筑科学研究院有限公司 Detection method for heat insulation performance of building enclosure and related device
CN103088937A (en) * 2013-01-28 2013-05-08 南京航空航天大学 Design method for inner side and outer side air layer thickness of condensation-preventing external thermal insulation wall body
CN104456698A (en) * 2014-11-05 2015-03-25 长沙理工大学 Building thermal environment and building energy saving control method for heating building
CN104634484A (en) * 2015-02-28 2015-05-20 中南大学 Method for obtaining winter heating running peak load based on changing thermal resistance of enclosure structure
CN111859514A (en) * 2020-07-30 2020-10-30 西安建筑科技大学 Method and system for optimizing thermal performance of building envelope under multi-working-condition operation
CN114719318A (en) * 2022-04-18 2022-07-08 亢天之 Low-temperature hot water radiation heating floor

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100385356C (en) * 2006-03-21 2008-04-30 纪义盛 Intelligent control method for public building heating and energy saving
DE102006062659A1 (en) * 2006-12-22 2008-06-26 Dr. Riedel Automatisierungstechnik Gmbh Ambient temperature and/or air quality regulating method for use by court, involves using reference variable, and using lateral characteristics of area brightness of area that is heated, as reference values

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101894183A (en) * 2010-07-02 2010-11-24 重庆星能建筑节能技术发展有限公司 Method for analyzing and designing green building
CN102095591A (en) * 2010-12-01 2011-06-15 西安建筑科技大学 Energy-saving performance estimating method of driven solar heating building
CN102095591B (en) * 2010-12-01 2013-10-23 西安建筑科技大学 Energy-saving performance estimating method of driven solar heating building
CN103033534A (en) * 2012-12-20 2013-04-10 深圳市建筑科学研究院有限公司 Detection method for heat insulation performance of building enclosure and related device
CN103088937A (en) * 2013-01-28 2013-05-08 南京航空航天大学 Design method for inner side and outer side air layer thickness of condensation-preventing external thermal insulation wall body
CN104456698A (en) * 2014-11-05 2015-03-25 长沙理工大学 Building thermal environment and building energy saving control method for heating building
CN104634484A (en) * 2015-02-28 2015-05-20 中南大学 Method for obtaining winter heating running peak load based on changing thermal resistance of enclosure structure
CN104634484B (en) * 2015-02-28 2017-07-07 中南大学 Heat supply in winter based on building enclosure change thermal resistance runs peak load preparation method
CN111859514A (en) * 2020-07-30 2020-10-30 西安建筑科技大学 Method and system for optimizing thermal performance of building envelope under multi-working-condition operation
CN111859514B (en) * 2020-07-30 2024-01-30 西安建筑科技大学 Method and system for optimizing thermal performance of building envelope under multi-station operation
CN114719318A (en) * 2022-04-18 2022-07-08 亢天之 Low-temperature hot water radiation heating floor
CN114719318B (en) * 2022-04-18 2024-02-09 亢天之 Low-temperature hot water radiation heating floor

Also Published As

Publication number Publication date
CN101476753B (en) 2011-05-11

Similar Documents

Publication Publication Date Title
CN101476754B (en) Construction thermal environment and construction energy-saving control method in air-conditioning construction
CN101476753B (en) Construction thermal environment and construction energy-saving control method in heating building
Andersson et al. The impact of building orientation on residential heating and cooling
CN104680004A (en) Building energy-saving rate calculation method
Wang et al. Energy performance and applicability of naturally ventilated double skin façade with Venetian blinds in Yangtze River Area
Hong et al. Annual energy performance simulation of solar chimney in a cold winter and hot summer climate
CN104680001A (en) Method for calculating architectural energy saving rate based on human body thermal adaptive model
CN102828565A (en) Building envelope structure for realizing energy saving by using wall heat storage and design method of building envelope structure
Liu et al. Heating load reduction characteristics of passive solar buildings in Tibet, China
Feng et al. Research on energy efficiency design key parameters of envelope for nearly zero energy buildings in cold area
CN109057002B (en) Building envelope structure containing heat accumulating layer and heat accumulating layer calculating method
CN110737937A (en) Existing building outer wall optimization application method based on reflective coating and vertical greening
Kowalski et al. Computational and the real energy performance of a single-family residential building in Poland–an attempt to compare: a case study
CN112035924B (en) Method for optimizing thickness of insulation layer of orientation-differentiated building enclosure
Siddique et al. Optimum insulation thickness for walls and roofs for reducing peak cooling loads in residential buildings in Lahore
Ibrahim et al. Optimizing insulation-thermal mass wall layer distribution from maximum time lag and minimum decrement factor point of view
Duan et al. Climate adaptive thermal characteristics of envelope of residential passive house in China
von Tils Effect of trees and greening of buildings on the indoor heating and cooling load–microscale numerical experiment
Liu et al. The research on the influence of building air tightness to energy consumption of residential building in a hot summer and cold winter zone in China
Takase et al. Practice of net positive energy house in the suburb of Tokyo
Li Application of energy-saving structural design under numerical simulation in solar heating buildings
Csáky et al. Analysis of degree day and cooling energy demand in educational buildings
Liu et al. Technical and Economic Analysis of External Wall Thermal Insulation for Nearly Zero Energy Buildings Based on Life Cycle Cost
Wang et al. Analysis of thermal energy storage system for energy saving reconstruction of building in region with heating provision and high sunshine
Djana et al. A Case Study on PV-Aided Net Zero-Energy Building: the Daycare in IKCU

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C53 Correction of patent for invention or patent application
CB03 Change of inventor or designer information

Inventor after: Huang Hengdong

Inventor before: Huang Hengdong

CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110511

Termination date: 20181231

CF01 Termination of patent right due to non-payment of annual fee