PHÂN TÍCH HỆ THANH BÊ TÔNG CỐT THÉP CÓ XÉT ĐẾN ĐẶC TRƯNG BIẾN DẠNG CỦA BÊ TÔNG

Chia sẻ: Nguyễn Văn Thắng | Ngày: | Loại File: DOC | Số trang:6

Thêm vào BST

Báo xấu

280
lượt xem 131
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Phân tích hệ thanh bê tông cốt thép có xét đến đặc trưng biến dạng của bê tông: biến dạng dẻo, khe nứt... nhằm đánh giá sự phân bố lại nội lực, để kết quả phân tích chính xác hơn và sử dụng vật liệu hiệu quả hơn.

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: PHÂN TÍCH HỆ THANH BÊ TÔNG CỐT THÉP CÓ XÉT ĐẾN ĐẶC TRƯNG BIẾN DẠNG CỦA BÊ TÔNG

PHÂN TÍCH HỆ THANH BÊ TÔNG CỐT THÉP CÓ XÉT ĐẾN ĐẶC TRƯNG BIẾN DẠNG CỦA BÊ TÔNG ANALYSIS OF REINFORCED CONCRETE FRAME STRUCTURE CONSIDERING THE DEFORMATION CHARACTERISTIC OF CONCRETE TRỊNH QUANG THỊNH Trường Đại học Bách khoa, Đại học Đà Nẵng TÓM TẮT Phân tích hệ thanh bê tông cốt thép có xét đến đặc trưng biến dạng của bê tông: biến dạng dẻo, khe nứt... nhằm đánh giá sự phân bố lại nội lực, để kết quả phân tích chính xác hơn và sử dụng vật liệu hiệu quả hơn. Lập thuật toán xây dựng ma trận độ cứng và vectơ ứng lực nút của phần tử thanh bê tông cốt thép. Kết hợp với phương pháp phần tử hữu hạn để phân tích kết cấu bê tông cốt thép trên máy tính điện tử. ABSTRACT Analysis of reinforced – concrete frame structure with consideration to deformation characteristic of concrete: plastic strain, cracks, and others is used to determine the redistribution of stress, to achieve more accurate analysis results and to use materials more effectively. The paper also presents the establishment of algorithm to build matrix of stiffness and node stress vector of reinforced – concrete frame elements, and in combination with the Finite Element Method to analyze reinforced-concrete structures on computer. 1. Mở đầu Biến dạng dẻo và sự có mặt của các khe nứt trong bê tông làm cho sự phân bố nội lực trong các kết cấu bê tông cốt thép khá phức tạp. Các phần mềm để phân tích kết cấu công trình phần nhiều là sử dụng phương pháp phần tử hữu hạn. Tuy nhiên việc xem xét một cách đầy đủ và chính xác ảnh hưởng của biến dạng phức tạp của bê tông trong việc xác định độ cứng và ứng lực nút của các phần tử mẫu cho phương pháp phần tử hữu hạn còn hạn chế. Bài viết này muốn đề xuất một bc’ biện pháp kết hợp các qui tắc tính toán hiện có để giải quyết khó khăn này. Fa’ hc’ 2. Cơ sở lý thuyết x h0 Xét dầm chịu uốn với tải trọng h tăng dần: lúc đầu chưa bị nứt, toàn bộ tiết diện bê tông chịu ứng suất, dầm có độ a cứng lớn. Tải trọng tăng, bê tông có biến hc dạng dẻo, vết nứt xuất hiện, tại tiết diện nứt mô men quán tính giảm làm giảm rõ rệt độ cứng của dầm. b Fa bc Hình 1: Tiết diện dầm tại vị trí khe nứt
2.1. Các đặc trưng hình học và biến dạng của cấu kiện bằng bê tông cốt thép thường sau khi xuất hiện khe nứt * Diện tích miền BT chịu nén qui đổi: n Fbqđ=(bc’ - b).hc’ + .Fa’ + b.x. (1 - 1) ν h0 Z1 * Độ cứng chống uốn của TD: B =  ψa  ψb   (1 - 2) +  E a Fa ν .Eb Fbqd    Trong đó: ν là hệ số đàn hồi của bê tông. Ζ 1 khoảng cách từ hợp lực vùng nén đến trọng tâm cốt thép chịu kéo. ψb: Hệ số xét đến biến dạng không đều của thớ bê tông chịu nén ngoài cùng, (lấy = 0.9 trong các trường hợp bình thường, khi chịu tải trọng rung động lấy = 1). ψa: Hệ số xét sự làm việc chịu kéo của bê tông; (ψa ≤ 1) * Độ cứng của thanh khi chịu kéo (nén): EF Cấu kiện chịu kéo lệch tâm có lực dọc Nc nằm trong phạm vi giữa các trọng tâm ' Ea .( Fa + Fa ) cốt thép Fa và Fa’: EF = ; (1 - ψa 3) ( E a .Fa + ν .Eb .Fbqd ) Cấu kiện chịu nén lệch tâm: EF = ; (1 - 4) ψa 2.2. Ma trận độ cứng và vectơ ứng lực nút tương đương của phần tử trong phương pháp phàn tử hữu hạn 3 Tổng quát với phần tử thanh ϕ1 - ϕ2 - lăng trụ trong bài toán phẳng (hệ thanh u1 - N1 M1 M2 phẳng): 1 v -Q u2 - N2 2 v2 - Q2 Gọi q1, q2 là tải trọng phân bố 1 1 theo phương trục 1, trục 2, m là mô men phân bố. Hình 2: Phần tử thanh phẳng Gọi t1, t2 là nhiệt độ tăng lên ở mặt dưới và mặt trên của phần tử. Các thành phần chuyển vị {V} = {u, v, ϕ} (u: chuyển vị thắng dọc trục 1, v: chuyển vị thắng theo trục 2, ϕ: chuyển vị Q xoay quanh trục 3) N Và nội lực {F} = {N, Q, M.} (qui ước dấu theo hình vẽ bên) M - Vectơ chuyển vị nút {δ } ≡ [δ 1 δ 2.. δ 6]T ≡ [u1 v1 ϕ 1 u2 v2 ϕ 2]T; - Vectơ ứng lực nút {F} ≡ [F1 F2.. F6]T ≡ [N1 Q1 M1 N2 Q2 M2]T; - Các quan hệ vi phân giữa chuyển vị và nội lực:
 du   dx   dv   1   t1 + t 2    0 0 0 E.F 0 0  u  α . 2   dx  0 v    dϕ  0 1 0 0 0    0  1   t −t   dx  0 0 0 0 0 ϕ    α . 2 1   dN  =  E.J  x   +  2.l  ; (2 - 1)   0 0 0 0 0 0  N   − q1   dx    Q     dQ  0 0 0 0 0 0     − q2   dx  0 M   0 0 0 −1 0      −m     dM     dx  dZ Hay = [A]. {Z} + {B}; (2 - 2) dx Theo phương pháp phần tử hữu hạn - Mô hình chuyển vị, có được phương trình cơ bản của phần tử hữu hạn: [K] {δ } = {F}; [K] -Ma trận độ cứng của phần tử. {δ } -Vectơ các chuyển vị nút. {F} -Vectơ các lực nút tương đương của phần tử. Với phần tử thanh có các đặc trưng tiết diện không thay đổi, có thể với các phép khai triển tích phân thông thường, hoặc với các tiết diện quen thuộc có thể sử dụng các bảng tra lập sẵn để xác định các phần tử của ma trận độ cứng và vectơ ứng lực nút tương đương của phần tử. 3. Phương pháp giải với cấu kiện bê tông cốt thép Với cấu kiện bê tông cốt thép, do các đặc trưng của tiết diện: độ cứng dọc trục EF, độ cứng chống uốn B (tương đương với EJ) xác định theo (1-3), (1-12) là thay đổi nên chưa thể giải được nghiệm đúng của hệ phương trình trên. Trong trường hợp này, có thể giải theo phương pháp gần đúng Runge-Kutta (bậc 4): Phương trình (2 - 2) được viết lại:  Z1 ( x )   a11 ( x) a12 ( x) .. a1n ( x)   Z1 ( x )   b1 ( x)   Z ( x)  a ( x) a ( x) .. a ( x)  Z ( x)  b ( x)  d  2   21 22 2n   2   2   ..  =  .. .. .. ..  x  ..  +  ..  ; (2 - 3) dx          ..   .. .. .. ..   ..   ..   Z n ( x)    an1 ( x) a n 2 ( x) .. a nn ( x)   Z n ( x ) bn ( x)       Nghiệm của hệ phương trình tại x: Z(x) = α(x).Z(0) + β(x); (2 - 4) Với khoảng chia h, tại x+h có: Z(x+h) = α(x+h).Z(0) + β(x+h); (2 - 5) Với các công thức truy hồi, có được sự liên hệ của các giá trị α và β ở 2 bước tích phân x và x+h như sau: β i ( x + h) = α ij ( x) + 6 ( 1 ' ) K i + 2( K i' ' + K i' ' ' ) + K iIV ;
1 ' α ij ( x + h) = α ij ( x ) + 6 ( ' ' IV ) K ij + 2( K ij' + K ij' ' ) + K ij ;  n  Với:  (  ) K i' = h. bi ( x) + ∑ β p ( x).aip ( x )  ;  p =1   h n   1  h  K i' ' = h. bi ( x + ) + ∑   β p ( x) + K 'p .aip ( x + )   ;   2 p =1   2  2     h n   1 ' h  K i' ' ' = h. bi ( x + ) + ∑   β p ( x) + K 'p .aip ( x + )   ;   2 p =1   2  2      n (( ) K iIV = h. bi ( x + h) + ∑ β p ( x ) + K 'p' .aip ( x + h) ' )  ;  p =1  n Và: ' ( K ij = h. ∑ aip ( x).α pj ( x) ; ) p =1 n  h  1  K ij' = h. ∑  aip ( x + ). α pj ( x ) + K 'pj   ; '   p =1  2  2  n  h  1 '  K ij' ' = h. ∑  aip ( x + ). α pj ( x ) + K 'pj   ; '   p =1  2  2  IV n ( ( K ij = h. ∑ aip ( x + h). α pj ( x) + K 'pj ; '' )) p =1 Vậy hàm Z(x) được xác định với mọi điểm nếu biết Z(0) và với các giá trị [ α(0)] = [I] và {β(0)} = {0}; ϕ 1 - M1 ϕ 2 - M2 Với phần tử (ik), (2 – 4) sẽ là: u1 - N1 i k Z(k) = α(k).Z(i) + β(k); u2 - N2 v1 - Q1 v2 - Q2 {Vk }   [α1 ] [α 2 ]  {Vi }  { β1}   = . + Hay { Fk }  [α 3 ] [α 4 ]  K { Fi }  { β 2 }  K      ⇒ {Vk} = [α1]k.{Vi} + [α2]k.{Fi} + {β 1}k; (2 - 6) {Fk} = [α3]k.{Vi} + [α4]k.{Fi} + {β 2}k; (2 - 7) (2 - 6) ⇒ {Fi} = [α 2 ] −1 .{Vk} - [α 2 ] −1 .[α1]k.{Vi} - [α 2 ] −1 .{β 1}k; k k k (2 - 7) ⇒{Fk}=([α3]k-[α4]k. [α 2 ] −1 .[α1]k).{Vi}+[α4]k. [α 2 ] −1 .{Vk}-[α4]k. [α 2 ] −1 .{β 1}k)+{β 2}k; k k k Mặt khác: {F} = [K] {δ } + {F0};
[ − { Fi }   K ik , i ] [ K ik , k ]. {Vi }  + {F i0 };   Hay:  = [  { Fk }   K ki, i ] [ Kk ]  {Vk }  {F }    0  k  Vậy, ma trận độ cứng và vectơ ứng lực nút tương đương của phần tử đang xét: [ K ik , i ] = [α 2 ] −1 .[α ] ; k 1 k [ K ik , k ] = - [α 2 ] −1 ; k {F } = [α 0 i ]−1 2 k .{β 1}k; [ K ki,i ] = [α ] - [α ] . [α 2 ] −1 .[α ] ; [ K ki, k ] = [α ] .( [α 2 ] −1 ; 3 k 4 k k 1 k k4 k {F } = {β } - [α ] . [α 0 k 2 k 4 k ]−1 2 k .{β 1}k; 4. Kết quả và thảo luận 2φ 18 Xét dầm có sơ đồ và tiết diện như hình vẽ, q=1.5T/m 500 1 2 3 2φ 2 6m 200 0 Với tiết diện ban đầu xem bê tông là vật liệu đàn hồi, độ cứng của dầm là EJ= 5.103 T.m2; Nội lực xác định được như sơ đồ. Xác định độ cứng của dầm bê tông cốt thép theo (1-3), có được sự phân bố lại nội lực như sơ đồ sau: - Theo kết quả phân tích, nếu có xét biến dạng dẻo và nứt trong bê tông thì độ cứng của dầm thay đổi đáng kể. Tại tiết diện có nội lực xấp xỉ giới hạn chịu lực, độ cứng có thể giảm 2-3 lần độ cứng ban đầu. - Giá trị của hệ số đàn hồi của bê tông giảm nhiều khi ứng suất tăng, do đó độ cứng của dầm không những thay đổi nhiều hơn, mà còn khác nhau nhiều giữa các tiết diện có nội lực chênh lệch lớn. Nếu xét đến yếu tố này sự phân bố lại nội lực càng rõ rệt. - Phân tích kết cấu bê tông cốt thép nếu không xét sự phân bố lại nội lực do độ cứng của các cấu kiện thay đổi sẽ làm mất ý nghĩa của những cố gắng chính xác hóa khi tính toán tiết diện bê tông cốt thép.
5. Kết luận - Với thuật toán trên, các đặc trưng biến dạng của bê tông là các số liệu rời rạc có được từ thực nghiệm có thể sử dụng để xác định độ cứng của từng tiết diện của cấu kiện tùy theo giá trị nội lực. - Khi phân tích kết cấu bê tông cốt thép nếu xét quan hệ giữa đặc trưng biến dạng và ứng lực trong bê tông, phải xác định đặc trưng biến dạng của bê tông khi chưa biết nội lực. Điều này sẽ làm cho khối lượng tính toán trở nên rất lớn nếu tính theo phương pháp lặp, vì vậy cần tìm một phương pháp tốt hơn để xác định quan hệ này. TÀI LIỆU THAM KHẢO [1] Ngô Thế Phong (Chủ biên), Kết cấu Bê tông cốt thép (Phần cấu kiện cơ bản), Tính toán cấu kiện Bê tông cốt théptheo trạng thái giới hạn thứ 2 1995: 83-90. [2] Nguyễn Quốc Tuấn, Bài giảng Áp dụng phương pháp phần tử hữu hạn phân tích kết cấu công trình. [3] NguyễnViết Trung, Thiết kế Kết cấu Bê tông cốt thép hiện đại theo tiêu chuẩn ACI, Chương 9: Khả năng khai thác. 1999: 384-405. [4] Arthur H. Nilson, Design of Concrete Structures-Twelfth Edition. Chapter 6 Serviceability 1997: 211-230.