Xác định hệ số phản ứng nhiệt của tường xây dựng ở điều kiện Việt Nam theo phương pháp hồi qui miền tần số

Mục đích của nghiên cứu này là xác định các hệ số phản ứng nhiệt của vách tường xây dựng trong điều kiện Việt Nam để sử dụng cho tính toán tải nhiệt. Trong nghiên cứu này, một chương trình Matlab được phát triển dựa trên phương pháp hồi quy miền tần số để thiết lập hàm truyền đa thức s dẫn nhiệt tức thời qua tường xây dựng.

Tạp Chí Khoa Học Giáo Dục Kỹ Thuật Số 39 (12/2016) 56 Trường Đại Học Sư Phạm Kỹ Thuật TP Hồ Chí Minh XÁC ĐỊNH HỆ SỐ PHẢN ỨNG NHIỆT CỦA TƯỜNG XÂY DỰNG Ở ĐIỀU KIỆN VIỆT NAM THEO PHƯƠNG PHÁP HỒI QUI MIỀN TẦN SỐ DETERMINATION OF THERMAL RESPONSE FACTORS OF BUILDING WALL IN VIET NAMCONDITION BASED ON FREQUENCY DOMAIN REGRESSION METHOD Lê Minh Nhựt Trường Đại học Sư phạm Kỹ thuật TPHCM Ngày tòa soạn nhận 05/9/2016, ngày phản biện đánh giá 6/10/2016, ngày chấp nhận đăng 01/11/2016 TĨM TẮT Mục đích nghiên cứu xác định hệ số phản ứng nhiệt vách tường xây dựng điều kiện Việt Nam để sử dụng cho tính tốn tải nhiệt Trong nghiên cứu này, chương trình Matlab phát triển dựa phương pháp hồi qui miền tần số để thiết lập hàm truyền đa thức s dẫn nhiệt tức thời qua tường xây dựng Các hệ số phản ứng nhiệt đạt dựa vào biến đổi Laplace hàm truyền đa thức s Kết tính tốn cho thấy hàm truyền đa thức s bậc truyền nhiệt phía ngồi, xun qua vào phía tường xây dựng tìm 55 điểm N=11 (7-2)+1 điểm tần số dải tần số 𝜔𝑘 từ 10-7 đến 10-2 radian/s phần trăm sai số hệ số truyền nhiệt U tổng hệ số phản ứng nhiệt truyền nhiệt phía ngồi, xun qua phía vách tường theo cơng thức 7.4410-9%, 5.8310-9%, 5.8610-9% Từ khóa: Hệ số phản ứng nhiệt; tường xây dựng; phương pháp hồi qui miền tần số; dẫn nhiệt tức thời; tải nhiệt ABSTRACT The objective of this research is to determine the thermal response factors of building wall to use for thermal load calculation of building in Viet Nam condition In this work, a MATLAB program was developed based on the frequency domain regression (FDR) method to estimate the polynomial s-transfer function of transient heat conduction of building wall The values of thermal response factors were obtained by applying Laplace transforms on the polynomial s-transfer function The calculated results showed that the fifth-order polynomial s-transfer functions of the outside heat conduction G X(s), across heat conduction G Y(s) and inside heat conduction G Z(s) of building wall were found at the 55 first points of N=11(7-2)+1 frequency points within the frequency range of 10-7 to 10-2 radian.s-1, and the error percentage between the thermal transmittance coefficient U and the sum of thermal response factors of the outside heat conduction, across heat conduction and inside heat conduction of the building wall were 7.4410-9%, 5.8310-9%, 5.8610-9%, respectively Keywords: Thermal response factors; building wall; frequency domain regression method; transien heat conduction; thermal load Tạp Chí Khoa Học Giáo Dục Kỹ Thuật Số 39 (12/2016) Trường Đại Học Sư Phạm Kỹ Thuật TP Hồ Chí Minh GIỚI THIỆU Trong kinh tế phát triển, tòa nhà cao tầng, khách sạn, khu nghỉ dưỡng nhà máy xí nghiệp ngày xây dựng nhiều địi hỏi bổ sung phương pháp tính tốn nhiệt vật liệu để đáp ứng yêu cầu xây dựng cơng trình xanh Việc tính toán tải nhiệt nhiệm vụ cốt yếu thiết kế mơ q trình nhiệt qua kết cấu bao che cơng trình Đối với tịa nhà khí hậu nhiệt đới, dịng nhiệt qua kết cấu xây dựng chiếm phần lớn phụ tải nhiệt, việc tính tốn xác dịng nhiệt điều quan trọng Thơng thường, tính tốn dịng nhiệt qua kết cấu xây dựng thường tính theo trạng thái ổn định đơn giản có nhược điểm độ xác khơng cao Thực tế dịng nhiệt qua cấu trúc xây dựng dịng nhiệt khơng ổn định(phụ thuộc thời gian) tính tốn dịng nhiệt khơng ổn định có ưu điểm dự đốn xác tải nhiệt tải đỉnh có nhược điểm tính tốn phức tạp Do nhà khoa học q trình xây dựng phương pháp tính tốn đơn giản hóa tốn dịng nhiệt tức thời qua kết cấu xây dựng dựa vào phương pháp chuyển đổi Laplace Laplace ngược đề xuất sớm Churchill Brisken[1,2] Dựa nghiên cứu Churchill Brisken, Stephenson Mitalas[3-6] cải tiến đề xuất phương pháp tính dịng nhiệt tức thời qua cấu trúc xây dựng dựa vào hệ số phản ứng nhiệt Trong phương pháp này, dòng nhiệt tức thời liên quan đến giá trị nhiệt độ giá trị trước Tuy nhiên, chuỗi giá trị hệ số phản ứng nhiệt cấu trúc xây dựng thường dài, trình tìm nghiệm hàm hyperbolic phức tạp dẫn đến tính tốn sai nghiệm, đặc biệt nghiệm gần Một phương pháp cải tiến tìm nghiệm mẫu số hàm đa thức dựa vào đổi dấu nghiệm trục 57 tọa độ đề xuất, phương pháp có nhược điểm tăng độ phức tạp tính tốn[7] Phương pháp khơng gian trạng thái đề xuất để tránh nghiệm, ưu điểm phương pháp hệ số phản ứng nhiệt cấu trúc xây dựng tìm mà khơng cần tìm nghiệm hàm truyền đa thức[8] Gần đây, Chen cộng [9] trình bày phương pháp điểm lưới dựa kỹ thuật nhận dạng hệ thống để xác định hàm truyền cấu trúc xây dựng kín từ liệu thí nghiệm Tác giả kết luận kỹ thuật tìm hệ số phản ứng nhiệt cấu trúc xây dựng không đồng Phương pháp hồi qui miền tần số phát triển để thiết lập hàm truyền dẫn nhiệt cấu trúc xây dựng nhiều lớp theo đặc trưng phản ứng tần số lý thuyết tính tốn dựa phép nhân nhiều ma trận bên phạm vi tần số xem xét[10] Wang Chen [11-13] đề xuất mơ hình tính tải nhiệt cho cấu trúc xây dựng dựa phương pháp hồi qui miền tần số để xác định hệ số hàm truyền dẫn nhiệt Theo kết tính tốn, tác giả khẳng định hàm truyền đa thức s đạt từ lý thuyết phản ứng tần số cấu trúc xây dựng dựa vào phương pháp hồi qui miền tần số tương đương hàm truyền hyperbolic s theo số hạng tần số đặc trưng Để giảm bớt bước tính tốn cho cấu trúc xây dựng có nhiều hệ số phản ứng nhiệt, phương pháp tính hệ số phản ứng nhiệt theo chu kỳ thời gian đề xuất, nghiên cứu này, hệ số phản ứng nhiệt theo chu kỳ xác định theo cực phần dư hàm truyền đa thức s[14] Để rút ngắn khoảng cách thời gian hai giá trị hệ số phản ứng nhiệt nhằm tăng độ xác phù hợp với giá trị thời gian thực tế thí nghiệm, Wang cộng [15] nghiên cứu kết luận cấu trúc xây dựng nhẹ trung bình hệ số hàm truyền dẫn nhiệt tính với bước nhảy thời gian 60 giây cấu trúc xây Tạp Chí Khoa Học Giáo Dục Kỹ Thuật Số 39 (12/2016) 58 Trường Đại Học Sư Phạm Kỹ Thuật TP Hồ Chí Minh dựng nặng 300 giây Tóm lại, để tính dịng nhiệt tức thời qua vách tường xây dựng hệ số phản ứng nhiệt đóng vai trị quan trọng Hiện phương pháp hồi qui miền tần số sử dụng rộng rãi để tính tốn hệ số phản ứng nhiệt có độ xác cao Tuy nhiên, loại vật liệu xây dựng khác hệ số phản ứng nhiệt hội tụ khác đặc trưng tần số khác khơng đổi, giá trị nhiệt độ T(x,0)=0, phương trình (1) (2) biểu diễn qua biến Laplace s liên quan đến dòng nhiệt nhiệt độ hai bề mặt cơng thức sau[1-7]: Mục đích nghiên cứu sử dụng phương pháp hồi qui miền tần số để tính tốn hệ số phản ứng nhiệt vách tường xây dựng theo tiêu chuẩn Việt Nam nhằm xây dựng sở liệu để sử dụng cho việc tính tốn dịng nhiệt tức thời tính tốn tải nhiệt, tải lạnh mơ dòng nhiệt tức thời để xác nhận lý thuyết thực nghiệm  A( s) B( s)  M (s)     M n ( s) M ( s) C ( s) D( s) Ti (s)  To (s)   A(s) B(s)  To (s)  q (s)  M (s)q (s)  C (s) D(s) q (s)  o   i   o   (3) M(s) ma trận tích tồn vách tường tính công thức (4): (4) Ở Mi(s) ma trân lớp thứ i tính sau:  Ai M i ( s)   Ci Bi  D i  (i=1,2…n) (5) CƠ SỞ LÝ THUYẾT phân tố biểu diễn qua hàm hyperbolic với biến Laplace s sau: 2.1 Phương pháp hồi qui miền tần số Ai  Di  cosh(Li s / ) (6) Bi   Ri sinh(Li s / ) /( Li s / ) (7) Ci   Li s / sinh(Li s / ) / Ri (8) Thông thường, cấu trúc xây dựng tịa nhà thường có nhiều hai lớp Trong nghiên cứu này, vật liệu lớp tường xây dựng xem đồng tính chất, hệ số dẫn nhiệt nhiệt dung riêng v.v…Phương trình vi phân dịng nhiệt truyền qua lớp tường xây dựng xem chiều có dạng tính sau:  2T ( x, t )  T ( x, t )  x  c p t (1) Ở đây, T trường nhiệt độ phụ thuộc tọa độ x, thời gian t Còn  ,  c p hệ số dẫn nhiệt, khối lượng riêng nhiệt dung riêng vật liệu Trong Li, Ri  (i / i c pi ) chiều dày, nhiệt trở hệ số dẫn nhiệt độ lớp thứ i tường xây dựng Công thức (3) viết lại dựa liên quan nhiệt độ dòng nhiệt hai bề mặt sau: qo ( s)  GX ( s) GY ( s)  To ( s)  q ( s )     G ( s ) G ( s )  T ( s )  Y Z  i   i   (9) Trong GX (s) , GY (s) GZ (s) Dòng nhiệt vị trí tọa độ x thời gian t tính sau: hàm truyền nhiệt phía ngồi, xun qua phía vách tường Đây hàm siêu việt Khi A(s) D(s)  C(s) B(s)  T ( x, t ) q( x, t )   x GX (s)  A(s) / B(s) (10) GY (s)  / B(s) (11) GZ (s)  D(s) / B(s) (12) (2) Khi thừa nhận tính chất vật liệu tường hệ số dẫn nhiệt, khối lượng riêng nhiệt dung riêng lớp tường GX (s) , GY (s) GZ (s) tính sau: Tạp Chí Khoa Học Giáo Dục Kỹ Thuật Số 39 (12/2016) Trường Đại Học Sư Phạm Kỹ Thuật TP Hồ Chí Minh Khi thay s  j vào công thức (10-12) 59 Các giá trị N điểm tầng số tìm [10 n1 ,10 n2 ] hàm trở thành hàm phức GX ( j ) , GY ( j ) GZ ( j ) , đặc N  10(n1  n2 )  tần số xác định theo trưng tần số truyền nhiệt phía ngồi, xun qua phía vách tường  k  10  n  ( k 1)( n  n 1 2) /( N 1) , (k=1,2…N).Hàm ~ G( jk )  Pk  jQk pháp hồi qui miền tần số sau[10-12]: ~    s   s    r s ~ B (s) G ( s)   ~  1s   s    m s m  A(s) vi phức công thức (14) biểu diễn thành hai thành phần phần thực phần ảo điểm thứ k cơng thức (15): Khi thừa nhận hàm truyền vách tường biểu diễn qua hàm truyền ~ hyperbolic G ( s) dạng đa thức dựa phương phạm (15) Các hệ số hàm đa thức (13,14) tìm dựa hàm chuẩn[11] sau: r (13) J ( )  ( H *  g * )T ( H  g ) Ở đây: i and  i hệ số thực, r (16) Ở dấu * lượng liên hợp phức,  T , m bậc tử số mẫu số Khi thay s  jk (k=1,2…N) vào cơng thức (13) ta gT có cơng thức tính sau:  T  [ 1     1    ] (17) ~    j   ( j )    r ( jk ) r B ( jk ) ~ (14) G( jk )  k k2  ~  1 jk   ( j k )    m ( jk ) m  A( j k ) g T  [G1 G2 GN ] 1  1 H  1   1  j1 ( j1 ) ( j1 ) ( j1 )   j1G1 H tính sau[11]:  ( j1 ) G1  ( j1 ) G1 j2 ( j2 ) ( j2 ) ( j2 )  j2G2  ( j2 ) G2  ( j2 ) G2 j3  j3G3  ( j3 ) G3  ( j3 ) G3 ( j3 ) ( j3 ) ( j3 ) j N ( j N ) ( j N ) ( j N )  (18)  ( j1 ) G1    ( j2 ) G2   ( j3 ) G3     ( j N ) GN     j N GN  ( j N ) GN  ( j N ) GN (19)  Khi thừa nhận tồn hàm   cực tiểu hàm mục tiêu   minj( ), Đặt   real ( H *T H ) 2  real ( H *T g  H T g * ) cơng thức (20) có dạng: hàm mục tiêu biểu diễn:        1  2H *T H   H *T g  H T g*  V0    V2  V     W   S2 W  S   W5    V V S W S W S V V V W S W S W (20) V V S W S W S V  V  V  W S W S W       (21) Giá trị   tính sau: W S W S W U  U    U S W S W S U U U W S W S W U U S W S W S U W S W S W  U 6 0 U U 8 U 10                    (22) Tạp Chí Khoa Học Giáo Dục Kỹ Thuật Số 39 (12/2016) 60 Trường Đại Học Sư Phạm Kỹ Thuật TP Hồ Chí Minh T  S0 W1  S2 W3 S4  U2 U4  (23) Giá trị Vi, Wi Ui Si tính sau: Vi  Ở đây: X j ( W / m2 K )là hệ số phản ứng nhiệt bên (j=0,1,2 n), Y j ( W / m2 K )là hệ số phản ứng nhiệt xuyên N N qua(j=0,1,2 n) Z j ( W / m2 K )là hệ số phản k 1 k 1 ứng nhiệt bên trong(j=0,1,2 n), Ti , To N N  ki , Si   ki Pk , Wi  ki Qk , Ui  ki (Qk2  Pk2 ) k 1 (24) k 1 2.2 Sử dụng kết đạt từ (22) để tính đặc trưng lý thuyết tầng số (10-12), giá trị hệ số phản ứng nhiệt xuyên qua vách tường xây dựng đạt từ cơng thức (11) tính sau: Y0  U  m i   (1  e nhiệt độ bề mặt bên bên si  Cấu trúc tường xây dựng Trong nghiên cứu này, tường đơn (bề dày qui ước: 110mm) gạch đất sét nung (QCVN 09:2013/BXD) hình chọn để tính tốn hệ số phản ứng nhiệt loại tường xây dựng phổ biến Việt Nam (25) ) i 1 Yj   m i   (1  e si  ( j 1) si  (26) ) e i 1 Ở U nhiệt trở vách tường xây dựng, thời gian 𝑡 = 𝑗∆𝜏(j=1, 2, 3…),  i phần dư GY (s) / s tính sau: ~  B (s ) i  i (27) si A( si ) ~ Ở A( si ) đạo hàm A( s) nghiệm thứ i si ( i=1,2,3 m) nghiệm ~ thứ m mẫu số  A( s) Giá trị hệ số phản ứng nhiệt bên bên X j Z j (j=1, 2, 3…) ~ tính tốn giống (23) (24) cho GX ( s) ~ GZ ( s) Dòng nhiệt tức thời qua vách tường bề mặt vách liên quan thời gian, nhiệt độ khứ tính sau[1-11]: n n qi (t )  Z oTi,t   Z jTi,t  j  YoTo,t   Y jTo,t  j j 1 (28) j 1 n n j 1 j 1 qo (t )   X oTo,t   X jTo,t  j  YoTi ,t   Y jTi ,t  j (29) Hình Cấu trúc tường đơn Chi tiết tính chất tường xây dựng thể bảng Bảng Chiều dày tính chất tường Mơ tả Chiều Hệ số Khối Nhiệt Nhiệt dày dẫn lượng dung trở δ(m) nhiệt riêng riêng R ρ cp (m2K/W) λ (W/mK) (kg/m3) (J/kgK) Lớp vữa 15 0.93 1800 840 0.0161 Lớp gạch 105 0.52 1300 880 0.2019 Lớp vữa 15 0.93 1800 840 0.0161 Trong nghiên cứu này, chương trình Matlab phát triển dựa vào mơ hình hóa tốn học bên để tính tốn giá trị hệ số phản ứng nhiệt vách tường xây dựng, sơ đồ thuật tốn mơ tả hình Tạp Chí Khoa Học Giáo Dục Kỹ Thuật Số 39 (12/2016) Trường Đại Học Sư Phạm Kỹ Thuật TP Hồ Chí Minh 61 N lựa chọn phụ thuộc vào tính chất vật liệu tường vật liệu có cấu trúc nặng, trung bình nhẹ Để tăng độ xác q trình tính tốn, đặc tính tầng số có giá trị từ N1 đến N2 bên N điểm tần số lựa chọn với điều kiện N1 1 N2  N để xây dựng hàm truyền đa thức s Trong suốt q trình tính tốn, hệ số hàm truyền đa thức s n, m, n1 , n2 , N1 N thay Bắt đầu Đọc: kích thước tính chất vách tường Tính: tầng số n1,n2, N, ωk Tính ma trận tích lớp vách Tính đặc tính tầng số GX,GY,GZ đổi thỏa mãn điều kiện hội tụ (công thức 30) Điều kiện để kiểm tra giá trị hệ số phản ứng nhiệt hay sai tổng hệ số phản ứng nhiệt phải hệ số truyền nhiệt U(W/m2K) vách tường cơng thức sau[7-15]: Tính: Hàm truyền đa thức s, cực phần dư hàm truyền Tính: Hệ số phản ứng nhiệt Xj,Yj,Zj Kiểm tra hội tụ M Yes  Stop Xj  j 1 No Hình Sơ đồ thuật tốn tính hệ số phản ứng nhiệt Trong tính tốn này, việc lựa chọn xác số điểm tần số quan trọng Các giá trị dãy tần số n1,n2 số điểm tần số M  Yj  j 1 M Z j U (30) j 1 Ở U hệ số truyền nhiệt vách tường tính sau: U (31) Rvua  Rgach  Rvua Trong đó: Rvữa(m2K/W), Rgạch(m2K/W) nhiệt trở lớp vữa lớp gạch tường 5.76596 E  2s  20.52861 s  1.57038 E  1s  3.27745 E  4s  1.65566 E  7s  9.33195 E  12 ~ G X ( s)  (32)  3.73796 E  3s  3.74384 E  5s  5.00301 E  7s  2.91567 E  9s  1.30533 E  11s  2.53427 E  14 ~ GY (s)  (33) s  1.16256 E  1s  1.76683 E  3s  7.48867 E  6s  9.05539 E  9s  2.18158 E  12 s  6.39119 E  3s  2.57115 E  5s  4.20341 E  8s  3.09404 E  11s  5.93646 E  15 5.76596 E  2s  20.52862 s  1.57038 E  1s  3.27745 E  4s  1.65566 E  7s  9.33195 E  12 ~ GZ ( s )  s  1.16256 E  1s  1.76683 E  3s  7.48867 E  6s  9.05539 E  9s  2.18158 E  12 (34) Bảng Hệ số phản ứng nhiệt vách tường xây dựng j X(j) (W/m2K) 17.692758884088555 -11.557014807303762 -1.277732832644018 -0.401232480732582 -0.127340528173747 -0.040421149410789 -0.012830741967018 -0.004072817033803 -0.001292819904215 Y(j) (W/m2K) 0.358414062908375 2.065745782193964 1.185332497637180 0.422490663005649 0.152298587250957 0.054940194440758 0.019816814201323 0.007147935853234 0.002578264071633 Z(j) (W/m2K) 17.692758892782006 -11.557014818820837 -1.277732833811779 -0.401232478023299 -0.127340527298692 -0.040421149132950 -0.012830741878825 -0.004072817005808 -0.001292819895328 Tạp Chí Khoa Học Giáo Dục Kỹ Thuật Số 39 (12/2016) 62 Trường Đại Học Sư Phạm Kỹ Thuật TP Hồ Chí Minh 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 -0.000410375249091 -0.000130263963695 -0.000041349229212 -0.000013125339564 -0.000004166330109 -0.000001322503429 -0.000000419797585 -0.000000133254863 -0.000000042298620 -0.000000013426701 -0.000000004261990 -0.000000001352868 -0.000000000429436 -0.000000000136314 -0.000000000043270 ∑ 4.270219489599711 4.270219489032414 4.270219489531991 4.270219489281749 4.270219489281749 4.270219489281749 U 0.000929981133487 0.000335444658397 0.000120995055470 0.000043642976812 0.000015742043488 0.000005678162932 0.000002048116200 0.000000738756534 0.000000266469850 0.000000096115808 0.000000034669020 0.000000012505133 0.000000004510608 0.000000001626978 0.000000000586852 -0.000410375246270 -0.000130263962799 -0.000041349228927 -0.000013125339473 -0.000004166330081 -0.000001322503420 -0.000000419797582 -0.000000133254862 -0.000000042298619 -0.000000013426701 -0.000000004261990 -0.000000001352868 -0.000000000429436 -0.000000000136314 -0.000000000043270 miền tần số phát triển để tính tốn Kết tính tốn theo phương pháp hồi hệ số phản ứng nhiệt vách tường xây qui miền tần số chương trình Matlab dựng điều kiện Việt Nam Hàm truyền đa phát triển cho thấy hàm truyền đa thức s truyền nhiệt phía ngồi, xun qua phía xác định phạm vi ~ thức s bậc truyền nhiệt phía ngồi GX ( s) , dải tần số n1, n2 số điểm tần số N 9,2 ~ ~ xuyên qua GY ( s) phía GZ ( s) vách 55 Sai số hệ số truyền nhiệt U tường thể công thức(32-34) tổng hệ số phản ứng nhiệt truyền Trong tính tốn này, giá trị phạm vi dãi nhiệt phía ngồi, xuyên qua phía tần số n1, n2 số điểm tần số N xác định vách tường 7.44x109% , 5.83x109% , 9,2 55( N  11(n1  n2 )  ) Hệ số phản 5.86x109% Sai số nhỏ chứng tỏ kết KẾT QUẢ VÀ THẢO LUẬN ứng nhiệt tìm dựa vào cơng thức(25-26) tương ứng cho truyền nhiệt phía ngồi, xun qua phía thể bảng 2, hội tụ tổng giá trị hệ số phản ứng nhiệt giá trị j=24 Phần trăm sai số hệ số truyền nhiệt U tổng hệ số phản ứng nhiệt dẫn nhiệt phía ngồi, xun qua phía vách tường theo công thức (28) 7.44x109 % , 5.83x109% , 5.86x109% KẾT LUẬN Trong nghiên cứu này, chương trình Matlab dựa phương pháp hồi qui tính tốn có độ xác cao Kết nghiên cứu sử dụng làm liệu để tính tốn mơ dịng nhiệt tức thời xuyên qua loại vách tường xây dựng phổ biến Việt Nam LỜI CẢM ƠN Tác giả chân thành cảm ơn tài trợ kinh phí nghiên cứu khn khổ đề tài nghiên cứu khoa học cấp trường trọng điểm T2016-65TĐ trường Đại học Sư phạm Kỹ thuật Tp.HCM Tạp Chí Khoa Học Giáo Dục Kỹ Thuật Số 39 (12/2016) Trường Đại Học Sư Phạm Kỹ Thuật TP Hồ Chí Minh 63 TÀI LIỆU THAM KHẢO [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] R.V Churchill, Operational mathematics, 2nded New York: McGraw-Hill, W N Brisken and S.G Reque, Heat load calculations by thermal response, ASHRAE Transaction,62, pp.391-424, 1956 G P Mitalas, Calculations of transient heat flow through walls and roofs, ASHRAE Transaction, 74, pp.182-188, 1968 G P Mitalas and D.G Stephenson, Room thermal response factors, ASHRAE Transaction, 73(2), III.1.1-1.7, 1967 G P Mitalas and D.G Stephenson, Cooling load calculation by thermal response method, ASHRAE Transaction, 73(1), pp.125-134, 1967 G P Mitalas and D.G Stephenson, Calculations of heat conduction transfer functions for multi-layer slabs, ASHRAE Transaction, 77, pp.117-126, 1971 D C Hittle and R Bishop, An improved root-finding procedure for use in calculating transient heat flow through multilayered slab, Int.J Heat and Mass Transfer, 26, pp.1685-1693, 1983 K Ouyang and F Haghighat, A procedure for calculating thermal response factors of multi-layer walls-state space method, Building and Environment, 26, pp.173-177, 1991 Y M Chen and Z K Chen, A neural – network based experimental technique for determining z-transfer function coefficients of a building envelope, Building and Environment, 35, pp.181-189, 2000 Y M Chen and S W Wang, Frequency-domain regression method for estimating CTF models of building multilayer constructions, Applied Mathematical Modelling, 25, pp.579-592, 2001 S W Wang and Y M Chen, A novel and simple building load calculation model for building and system dynamic, Applied Thermal Engineering, 21, pp.683-702, 2001 S W Wang and Y M Chen, A simple procedure for calculating thermal response factors and conduction transfer functions of multilayer, Applied Thermal Engineering, 22, pp.333-338, 2002 S W Wang and Y M Chen, Transient heat flow calculation for multilayer construction using a frequency domain regression, Building and Environment, 38, pp.45-61, 2003 Y M Chen and S W Wang, A new procedure for calculating periodic response factors based on frequency domain regression method, In.l Journal of Thermal Sience, 44, pp.382-392, 2005 J Wang, S Wang and et al, Short time step heat flow calculation of building constructions based on frequency domain regression method, Int Journal of Thermal Sience, 48, pp.2355-2364, 2009 Tác giả chịu trách nhiệm viết: Lê Minh Nhựt Trường Đại học Sư phạm Kỹ thuật TP.HCM Email: nhutlm@hcmute.edu.vn ... cho cấu trúc xây dựng có nhiều hệ số phản ứng nhiệt, phương pháp tính hệ số phản ứng nhiệt theo chu kỳ thời gian đề xuất, nghiên cứu này, hệ số phản ứng nhiệt theo chu kỳ xác định theo cực phần... -0.000000000043270 miền tần số phát triển để tính tốn Kết tính tốn theo phương pháp hồi hệ số phản ứng nhiệt vách tường xây qui miền tần số chương trình Matlab dựng điều kiện Việt Nam Hàm truyền... cứu sử dụng phương pháp hồi qui miền tần số để tính toán hệ số phản ứng nhiệt vách tường xây dựng theo tiêu chuẩn Việt Nam nhằm xây dựng sở liệu để sử dụng cho việc tính tốn dịng nhiệt tức thời