Đang tải... (xem toàn văn)
Neo gia cố là một trong những giải pháp quan trọng được sử dụng gần đây trong xây dựng công trình ngầm, vì nó ảnh hưởng trực tiếp đến ứng suất, biến dạng và độ ổn định của kết cấu cũng như giá thành xây dựng của công trình. Theo nguyên lý làm việc, neo gia cố có tác dụng liên kết vỏ hầm với khối đá xung quanh tạo thành một chỉnh thể làm việc đồng thời, độ cứng khối đá và ổn định của vỏ hầm được tăng lên.
KHOA HỌC CƠNG NGHỆ PHÂN TÍCH ẢNH HƯỞNG NEO GIA CỐ ĐẾN ỨNG SUẤT, BI ẾN DẠNG KẾT CẤU CÔNG TRÌNH NGẦM BẰNG PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN Đào Văn Hưng Trường Đại học Thủy lợi Tóm tắt: Neo gia cố giải pháp quan trọng sử dụng gần xây dựng cơng trình ngầm, ảnh hưởng trực tiếp đến ứng suất, biến dạng độ ổn định kết cấu giá thành xây dựng cơng trình Theo ngun lý làm việc, neo gia cố có tác dụng liên kết vỏ hầm với khối đá xung quanh tạo thành chỉnh thể làm việc đồng thời, độ cứng khối đá ổn định vỏ hầm tăng lên Nghiên cứu sử dụng phương pháp phần tử hữu hạn để phân tích ổn định hầm nhà máy thủy điện với trường hợp phân bố neo gia cố khác Kết tính tốn cho thấy việc nghiên cứu ảnh hưởng neo gia cố đến ứng suất, biến dạng kết cấu cơng trình ngầm cần thiết để đảm bảo điều kiện kỹ thuật kinh tế Từ khóa: cơng trình ngầm, neo gia cố, kết cấu, ứng suất, biến dạng Summary: Anchoring reinforcement (rock bolts) is currently among important reinforcing solution available in the underground engineering, as it directly affects to the stresses, deformations and stability of structures, as well as construction costs The working principle is that rock bolts anchor the lining to the surrounding rock mass forming a simultaneously-working structure, adding strength to the rock mass, and increasing the stability of the tunnel lining This study employs the finite element method (FEM) to analyze the stability of an undergroundhydro power plantin different anchorschemes The calculation results confirmed that it is necessary to investigate the influences of rock bolts to the stresses and deformations of underground structures in order to satisfy technical and economicconditions Keywords: Underground structu res, rock bolts, stress, deformation, FEM ĐẶT VẤN ĐỀ * Neo gia cố trở thành giải pháp quan trọng để tăng ổn định vỏ hầm đá xung quanh Neo xuyên vào khối đá nhằm tăng cường liên kết, tạo thành khối chỉnh thể, đồng thời tăng khả chống uốn, chống cắt tính ổn định khối đá xung [1], [2], [4] Trước thiết kế, thi công quanh neo gia cố chủ yếu theo kinh nghiệm thực tế so sánh loại cơng trình dựa vào địa chất cơng trình đưa tham số thiết kế neo Do kết cấu địa chất phức tạp, trường ứng lực ban đầu không gian kết cấu đá xung quanh thông thường hiệu ứng ba chiều, nên khó chuyển hóa xác loại hiệu ứng này, dẫn đến kết tính khơng thể phản ánh làm việc neo gia cố khối đá xung quanh Đến nay, với phát triển cơng trình ngầm ứng dụng khoa học kỹ thuật mới, có nhiều phương pháp, phần mềm tính tốn neo gia cố đáp ứng yêu cầu thiết kế thi công Ngày nhận bài: 06/01/2017 Ngày thông qua phản biện: 20/2/2017 Ngày duyệt đăng: 28/2/2017 Bài báo dựa theo nguyên lý tác dụng neo gia cố, chuyển hóa phần tử neo hình trụ dạng ẩn mơ hiệu ứng neo gia cố, thơng qua phân tích tính đàn hồi dẻo ba chiều TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ THỦY LỢI SỐ 37 - 2017 KHOA HỌC CÔNG NGHỆ phương pháp phần tử hữu hạn ứng suất, biến dạng cơng trình ngầm CƠ SỞ TÍNH TỐN 2.1 Mơ hình học neo gia cố Dựa vào nhiều tài liệu thí nghiệm quan trắc thự tế, chế hoạt động chủ yếu neo gia cố cơng trình ngầm “tổ hợp cố kết” “hệ thống treo hình cung” Neo xuyên vào đá xung quanh hầm có tác dụng liên kết thành khối thống nhất, tăng thêm lực ma sát khối đá, đồng thời tăng khả chống uốn chống cắt đá xung quanh Khi thi cơng cơng trình ngầm, đá xung quanh hầm bị tác động làm thay đổi trạng thái cân tự nhiên ban đầu, trở nên rời rạc, khả ổn định cơng trình ngầm[1], [4], [5] Do đó, để đảm bảo khả làm việc ổn định cơng trình cần phải gia cố neo xuyên vào khối đá xung quanh M ật độ, số lượng neo đạt đến mức độ định để hệ thống neo gia cố hình thành cung tròn chịu lực, đá xung quanh thành chỉnh thể tăng khả chịu lực bên hạn chế chuyển vị vỏ hầm, làm khối đá tự cân ổn định Theo chế chịu lực neo gia cố đặc điểm tính tốn phần tử hữu hạn ba chiều, báo chọn mơ hình giới hạn neo gia cố khối đá xung quanh hình trụ tròn Khi đó, mơ hình mơ làm việc neo xuyên vào khối đá, cố kết với đá xung quanh hình thành chỉnh thể (Hình 1) M hình phần tử tính tốn có lõi tâm trụ neo, bên khối đá xung quanh, phạm vi ảnh hưởng chiều rộng chiều sâu xác định từ mơ hình thí nghiệm Trong tính tốn hình trụ neo đá xung quanh có xét đến tính dị hướng lớp địa tầng theo thông số lý Từ định luật Hooke tổng quát mối quan hệ ứng suất biến dạng tốn khơng gian hình [3], [5] : tụ (1) Do phần tử tính tốn hình trụ tròn có trục đối xứng , suy quan hệ hướng tâm, nên biến dạng ứng suất biến dạng đơn giản hóa: (2) đó: m, n tham số ma trận đàn hồi (với Er, µr Ez, µz modun đàn hồi, hệ số Poisson theo hướng kính theo trục z) Theo mơ hình học phần tử neo, xây dựng phần tử hình trụ đẳng tham số chiều dọc trục z hướng bán kính r neo (Hình 2), tương đương phần tử tính tốn hình 1, thị điểm nút, phân biệt tọa độ cục nút (z, r) điểm nút tính tốn là: 1(-1, 0); 2(1, 0); 3(-1, 1); 4(1, 1) Hình Mơ hình phần tử neo gia cố khối đá xung quanh TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ THỦY LỢI SỐ 37 - 2017 Hình Mơ hình phần tử hình trụ đẳng tham số KHOA HỌC Để đảm bảo liên tục theo điều kiện biên, hàm dạng Ni phần tử tính tốn phải thỏa mãn Ni = nút i nút khác[3], [6] , đó: CÔNG NGHỆ vật lý hàm dạng, thay vào (3) ta có: C=1/2 Tương tự, xác định biểu thức hàm dạng lại phần tử: (3) Tại nút (N1 = 1, nút lại 2, 3, 4) có r = 0, z = -1, N1(z, r) phải có dạng: (4) Với C số xác định theo ý nghĩa (5) Ta có ma trận độ cứng phần tử trụ neo: (6) đó: |J| định thức ma trận Jacobian; [B] ma trận hình dạng phần tử trụ, với ma trận phần tử [B] là: Dựa vào tích phân hàm Gauss, với t = 2r-1 thay vào (6) ta có: (7) 2.2 Chuyển đổi neo dạng ẩn Neo gia cố cơng trình ngầm thường phụ thuộc vào địa chất nên phân bố không đều, tính tốn phần tử hữu hạn chiều biểu diễn neo thành phần tử kết cấu thể Để tiện lợi cho việc nghiên cứu phương thức phân bố neo khác nhau, mà lại không ảnh hưởng đến chia lưới phần tử hữu hạn, báo ứng dụng phần tử neo dạng ẩn để tiến hành phân tích neo gia cố Gọi phần tử dạng ẩn tức ẩn dấu phần tử neo gia cố phần tử khối đá, chia phần tử khối đá khơng có tác dụng trực tiếp neo phải tăng độ cứng khối đá đến độ cứng chỉnh thể chiều phần tử hữu hạn Từ (7) tính độ cứng phần tử trụ neo, hệ tọa độ hình trụ ma trận độ cứng chiều, biểu thị ma trận phần tử: (8) Trong hệ tọa độ góc vng ma trận độ cứng chiều, chọn hệ trục tọa độ chiều trục neo làm tọa độ cục (hình 3), độ cứng phần tử neo có biểu thị là: (9) Hình Phần tử neo gia cố dạng ẩn Để độ cứng phần tử neo dạng ẩn tăng đến độ cứng tổng thể, trước hết phải xác định độ cứng tọa độ cục theo cơng TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ THỦY LỢI SỐ 37 - 2017 KHOA HỌC CÔNG NGHỆ thức (9) chuyển đổi thành độ cứng hệ tọa độ chỉnh thể M a trận chuyển hóa tọa độ cục tọa độ chỉnh thể chiều phần tử neo là: (10) Trong công thức lij biểu thị cosin phương tọa độ cục i tọa độ chỉnh thể chiều j neo (i, j đại diện phương x, y, z) Theo nguyên lý chuyển đổi ma trận độ cứng phần tử hữu hạn, tọa độ chỉnh thể, ma trận độ cứng phần tử neo là: Giả thiết tải trọng điểm nút phần tử neo dạng ẩn là: (14) Theo nguyên lý tương đương lực tĩnh phần tử hữu hạn, phân bổ tải trọng tương đương điểm nút phần tử neo đến điểm nút i phần tử khối đá (15) Theo phương trình cân phần tử phương pháp phần tử hữu hạn, phân biệt phần tử khối đá neo dạng ẩn phần tử khối đá: (16) (11) Từ công thức (11) xác định ma trận độ cứng phần tử neo dạng ẩn, dựa theo lý thuyết [5] phép nội suy phần tử hữu hạn , tăng độ cứng dạng ẩn lên đến độ cứng phần tử khối đá Từ hình cho thấy, giả thiết giá trị chuyển vị điểm nút phân tử neo dạng ẩn {δ}’ chuyển vị điểm nút phần tử khối đá {δ} tọa độ chỉnh thể tương ứng là: Theo lý thuyết phép nội suy phần tử hữu hạn, suy mối liên hệ chuyển vị điểm nút phần tử neo chuyển vị điểm nút phần tử khối đá: (12) [N] ma trận tham số hình học: (13) với giá trị hàm số hình học điểm nút thứ i phần tử khối đá nằm với nút thứ j phần tử neo gia cố (17) Thay (12) (15) vào (16) ta có (18) Từ (17) (18), sau phân bổ tương đương ma trận độ cứng phần tử neo dạng ẩn đến phần tử khối đá, ma trận độ cứng là: (19) [ke ] [N] phân biệt từ công thức (11) (13) 2.3 Phương pháp phân tích tính đàn hồi dẻo neo gia cố Khi thi cơng cơng trình ngầm làm thay đổi kết cấu tự nhiên ban đầu khối đá, ứng suất khối đá trở thành tải trọng giải phóng tác dụng lên vỏ hầm, làm vỏ hầm xuất biến dạng chuyển vị Do khối đá có tính đàn hồi dẻo, nên sau thi công biến dạng đá xung quanh cần khoảng thời gian định cân trở lại Do đó, sau thi cơng cơng trình ngầm cần phải gia cố neo kịp thời, để có hiệu cao việc ngăn cản biến dạng đá xung quanh hầm Theo đó, báo ứng dụng phương pháp phần tử hữu hạn đàn hồi chiều, gia tăng độ cứng đạt tính dẻo, từ phân tích tác dụng neo gia cố cơng trình ngầm Phương pháp nghĩa phân tải trọng thi công lên kết cầu TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ THỦY LỢI SỐ 37 - 2017 KHOA HỌC cơng trình ngầm thành phần: tải trọng có tính đàn hồi {Re} tải trọng có tính dẻo {Rp} Với tải trọng đàn hồi tác dụng lên kết cấu, tiến hành tính tốn độ cứng kết cấu hầm, sau gia tăng độ cứng neo gia cố để tạo hiệu ứng đến độ cứng chỉnh thể kết cấu hầm; tải trọng tính đàn dẻo tính tốn tương tự Trong phép lặp phân tải, lần lặp lại phải tiến hành tính tốn ứng suất phần tử khối đá xung quanh phần tử neo gia cố Đối với phép lặp ứng suất tăng thêm phần tử khối đá tính theo biến dạng tăng thêm sau: (20) ma trận ứng suất đàn hồi đàn dẻo phần tử đá; S hệ số ứng suất đàn hồi phần tử Với phép lặp ứng suất tăng thêm phần tử neo, chuyển vị nút neo dựa theo chuyển vị điểm nút phần tử khối đá, dùng lý thuyết phép nội suy để tìm chuyển vị điểm cuối neo[4] Do đó, tính biến dạng phần tử neo tính theo cơng thức ứng suất tăng thêm neo: (21) E, modul đàn hồi hệ số poisson vật liệu neo gia cố CƠNG NGHỆ tính tốn phân tích hiệu làm việc neo gia cố Cơng trình nhà máy thủy điện lòng đất có hầm nhà máy chính, hầm biến áp chính, cửa lấy nước, đường hầm chính, đường giao thơng, hầm nước,… hợp thành khơng gian ngầm lòng đất Đá khu nhà máy chủ yếu đá hoa cương vân đen, modul đàn hồi E=35GPa, hệ số poison µ=0.22, lực kết dính c=1,3MPa, góc ma sát =500 Hầm nhà máy nằm độ sâu 220330 m lòng núi, kích thước hầm nhà máy 120mx22mx46m (dài x rộng x cao) Chia lưới phần tử mặt cắt ngang nhà máy phân bố neo hình Để nghiên cứu ảnh hưởng neo gia cố đến ứng suất, biến dạng đá xung quanh, tiến hành tính tốn với trường hợp (TH) sau: - Trường hợp 1: khơng có neo gia cố; - Trường hợp 2: neo đường kính 2.5cm, khoảng cách phân bố neo 3x3m ; chiều dài neo 3m - Trường hợp 3: sở trường hợp 1, gia cố thêm neo 2,5cm trung tâm 3x3m , chiều dài 5,5m - Trường hợp 4: sở trường hợp 1, tăng thêm mật độ phân bố neo 1,5x1,5m , đường kính 2,5cm, chiều dài 3m Dùng phần tử neo dạng ẩn, sử dụng phần tử hữu hạn đàn dẻo ba chiều tiến hành tính tốn Trong q trình lặp, nhập trạng thái dẻo phần tử khối đá phần tử neo, phải dựa theo phương pháp độ cứng đàn dẻo tăng thêm để biến đổi độ cứng, sau thực vòng lặp, tiếp tục lặp phép tính hồn tất NỘI DUNG TÍNH TỐN 3.1 Ví dụ tính tốn Ví dụ tính tốn cho cơng trình thủy điện có hầm nhà máy nằm lòng đất, ứng dụng phần tử hữu hạn đàn dẻo ba chiều tiến hành Hình Bố trí neo chia lưới phần tử TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ THỦY LỢI SỐ 37 - 2017 KHOA HỌC CƠNG NGHỆ 3.2 Phân tích kết tính tốn: Từ kết tính tốn với trường hợp khác nhau, có nhận xét cho thấy tác dụng neo gia cố xây dựng cơng trình ngầm rõ rệt, cụ thể như: - Hạn chế biến dạng đá xung quanh: Từ bảng thấy, sử dụng hệ thống neo gia cố TH2, so với TH1 khơng gia cố chuyển vị lớn đỉnh hầm nhà máy giảm 0,6÷1,2mm, chiếm 15%÷23%; chuyển vị lớn bên tường giảm 2,8÷3,6mm, chiếm 33÷43% So sánh với TH3, mật độ neo tăng thêm TH4, biến dạng đá xung quanh tiếp tục giảm, khơng đáng kể Điều cho thấy, sau thi công xong, neo xuyên vào khối đá, tác dụng liên hợp neo đá xung quanh tạo thành khối thống nhất, đồng thời làm độ cứng chỉnh thể nâng lên ngăn cản biến dạng đá xung quanh Khi mật độ neo tăng đến mức độ định, tiếp tục tăng mật độ neo hiệu thấp Do đó, cần phải nghiên cứu trường hợp gia cố với mật độ hợp lý để cơng trình đạt mục đích vừa an tồn ổn định vừa đảm bảo điều kiện kinh tế Bảng Kết tính chuyển vị xung quanh hầm với trường hợp tính tốn: (Đ/v: mm) Hạng mục Vị trí Đỉnh hầm Nhà Tường mày thượng lưu Tường hạ lưu Đỉnh hầm Phòng Tường biến thượng lưu áp Tường hạ lưu TH1 không gia Có gia cố TH2 Có gia cố TH3 Có gia cố TH4 cố Phương Phương Phương Phương Phương Phương Phương Phương ngang đứng ngang đứng ngang đứng ngang đứng 2.48 4.58 1.83 -3.85 1.65 -3.44 1.59 -3.41 8.55 -0.34 5.73 -0.25 4.88 -0.28 4.69 -3.01 -9.48 1.83 -6.66 1.62 -6.15 1.50 -5.91 1.39 0.21 -4.55 0.35 -3.74 0.25 -3.36 0.29 -3.31 3.26 -2.57 0.273 -1.72 2.09 -1.60 1.99 -1.61 -5.74 0.15 -4.39 0.14 -3.83 0.15 -3.71 0.19 - Giảm trạng thái ứng suất xung quanh hầm - Hạn chế phát triển vùng đàn dẻo đá xung quanh Từ bảng thấy, sử dụng neo gia cố xung quanh hầm, không giảm ứng suất chung mà cường độ ứng s uất tập trung giảm theo Trạng thái ứ ng suất ba chiều có neo so với khơng gia cố có xu hư ớng đồng hóa Nó cho thấy rằng, sau neo đá xung quanh hình thành tổ hợp chỉnh thể, tải trọng tác dụng lên đá xung quanh phần truyền lên neo Do đó, ứng suất đá xung quanh đư ợc cải thiện đồng thời làm tăng khả chịu tải đá quanh hầm So sánh kết biểu đồ phổ màu khu vực đàn dẻo hình hình 6, cho thấy sau gia cố, neo đá xung quanh làm việc đồng thời, lực ma sát đá xung quanh tăng lên tăng khả chống cắt Khu vực đàn dẻo đá xung quanh hầm nhà máy nhà biến áp giảm đáng kể, vùng ảnh hưởng bị thu hẹp So sánh phân bố khu vực đàn dẻo trường hợp hỗ trợ 2, có giảm không đáng kể Điều cho thấy, bố trí neo gia cố cách hợp lý vừa đảm bảo tính ổn định đá xung quanh, vừa đáp ứng yêu cầu thời gian chi phí xây dựng TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ THỦY LỢI SỐ 37 - 2017 KHOA HỌC CÔNG NGHỆ Bảng S ự biến đổi ứng suất chu vi hầm với trường hợp tính tốn (Đơn vị: MPa) Trường hợp Ứng suất ban đầu TH1 Không gia cố Gia cố TH2 Gia cố TH3 Gia cố TH4 x Đỉnh hầm y z Tường thượng lưu x y z Tường hạ lưu x y z Đường hầm x y z -8,2 -6,7 -11,7 -8,6 -7,1 -12,3 -8,4 -6,9 -11,9 -8,6 -7,1 12,5 -12,4 -6,2 -5,1 -2,8 -5,9 -13,2 -3,6 -7,1 -14,7 -5,5 -8,9 -21,3 -11,3 -6,1 -5,9 -3,8 -6,3 -13,4 -4,8 -6,6 -14,4 -5,7 -8,2 -20,4 -11,0 -6,2 -6,3 -4,1 -6,3 -13,2 -5,2 -6,7 -14,2 -6,2 -8,2 -19,9 -11,0 -6,2 -6,5 -4,0 -6,3 -13,2 -5,2 -6,6 -14,1 -6,2 -8,2 -19,7 Ghi chú: trục x hướng thượng lưu hạ lưu; trục y dọc theo trục nhà máy, z hướng thẳng đứng Hình Phân bố khu vực đàn dẻo trường hợp khơng có neo gia cố KẾT LUẬN Dùng phương pháp phẩn tử hữu hạn mơ hình hóa phần tử neo dạng ẩn mơ neo gia cố để tính ứng suất, biến dạng đá xung quanh cơng trình ngầm cho ta kết sát thực tế, phản ánh làm việc đồng thời neo đá xung quanh Khi chuyển hóa phần tử neo dạng ẩn giúp giải toán nhanh chóng, tiện lợi kể trường hợp phân bố neo không đều, địa chất phức tạp, kết cho thấy hiệu neo gia cố ổn định đá xung quanh cơng trình ngầm Ngun lý làm việc neo gia cố thông qua tác dụng liên hợp neo đá xung Hình Phân bố khu vực đàn dẻo trường hợp có neo gia cố TH1 quanh, có hiệu ngăn cản đá xung quanh biến dạng, cải thiện trạng thái ứng suất Khi neo xuyên vào khối đá, khối đá neo làm việc đồng thời nên hạn chế phát triển khu đàn hồi dẻo đá xung quanh hầm, từ tăng khả chịu lực bên nâng cao tính ổn định vỏ hầm Phân bố neo gia cố khác tính ổn định đá xung quanh khác Khi gia tăng mật độ neo độ cứng tính chỉnh thể khối đá nâng lên Nhưng sau mật độ neo gia tăng đến mức độ định, tiếp tục gia tăng mật độ neo tác dụng làm giảm ứng suất, biến dạng, ổn định kết cấu TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ THỦY LỢI SỐ 37 - 2017 KHOA HỌC CƠNG NGHỆ cơng trình Do đó, việc nghiên cứu biện pháp, mật độ neo gia cố việc làm cần thiết để cơng trình ngầm đạt u cầu kỹ thuật kinh tế./ TÀI LIỆU THAM KHẢO [1] Trần Thanh Giám, Tạ Tiến Đạt Tính tốn thiết kế cơng trình ngầm; NXB Xây dựng; 2002 [2] Võ Trọng Hùng, Phùng M ạnh Đắc, (2005), Cơ học đá ứng dụng xây dựng cơng trình ngầm khai thác mỏ NXB Khoa học Kỹ thuật [3] Trần Ích Thịnh, Ngô Như Khoa, (2007), Phương pháp phần tử hữu hạn, NXB Khoa học Kỹ thuật [4] Guo Lingyun, Xiao M ing, Ren Yi, (2007), Numerical Simulation and M echanical Analysis of End-Anchored Bolts, Chinese Journal of Rock Machanics and Engineering, Vol.26, No.2, p4221-4226 (Trung văn) [5] Ru Zhongliang, nnt, (2005), Parallel Finite Element M ethod Anlysis of Bolted Rock M ass in Underground Engineering, Chinese Journal of Rock Machanics and Engineering, Vol.24, No.1, p13-17 (Trung văn) [6] A A Elsayed, (2011), Study of Rock-Lining Interaction for Circular Tunnels Using Finite Element Analysis, Jordan Journal of Civ il Engineering, Volume 5, No 1, p5063 (Anh văn) TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ THỦY LỢI SỐ 37 - 2017 ... hướng thẳng ứng Hình Phân bố khu vực đàn dẻo trường hợp khơng có neo gia cố KẾT LUẬN Dùng phương pháp phẩn tử hữu hạn mơ hình hóa phần tử neo dạng ẩn mơ neo gia cố để tính ứng suất, biến dạng đá... phương trình cân phần tử phương pháp phần tử hữu hạn, phân biệt phần tử khối đá neo dạng ẩn phần tử khối đá: (16) (11) Từ công thức (11) xác định ma trận độ cứng phần tử neo dạng ẩn, dựa theo... Tường biến thượng lưu áp Tường hạ lưu TH1 khơng gia Có gia cố TH2 Có gia cố TH3 Có gia cố TH4 cố Phương Phương Phương Phương Phương Phương Phương Phương ngang ứng ngang ứng ngang ứng ngang đứng