Bài báo trình bày thuật toán phần tử hữu hạn (PTHH) và phương pháp phân tích động lực học của hệ liên hợp giàn thép không gian – bể chứa trên nền san hô chịu tác dụng đồng thời của tải trọng sóng và gió với việc sử dụng mô hình kết cấu và nền san hô làm việc đồng thời. Các kết quả nghiên cứu số là cơ sở khoa học cho phép lựa chọn các giải pháp hợp lý đối với công trình bố trí tại các bãi cạn ven đảo, góp phần giảm mật độ công trình trên đảo, tăng khả năng tác chiến trên đảo.
Tuyển tập công trình Hội nghị Khoa học toàn quốc Cơ học Vật rắn biến dạng lần thứ XII, Đà Nẵng, 6-7/8/2015 Tập ISBN 978-604-913-458-6 Đơn vị tổ chức ĐẠI HỌC DUY TÂN HỘI CƠ HỌC VẬT RẮN BIẾN DẠNG Các nhà khoa học tham gia phản biện Phạm Hoàng Anh, Vũ Quốc Anh, Vũ Mai Ba, Đào Huy Bích, Nguyễn Đăng Bích, Lương Xuân Bính, Phạm Đức Chính, Nguyễn Thái Chung, Nguyễn Tiến Chương, Nguyễn Mạnh Cường, Đào Văn Dũng, Nguyễn Thế Dương, Phạm Tiến Đạt, Nguyễn Đình Đức, Đặng Vũ Hiệp, Hoàng Phương Hoa, Vũ Công Hòa, Nguyễn Thái Tất Hoàn, Đặng Xuân Hùng, Lê Thái Hùng, Nguyễn Xuân Huy, Vũ Lê Huy, Lê Xuân Huỳnh, Ngô Như Khoa, Nguyễn Việt Khoa, Phan Bùi Khôi, Nguyễn Đình Kiên, Nguyễn Trung Kiên (ĐHSPKTTPHCM), Bùi Hải Lê, Trần Văn Liên, Vũ Đỗ Long, Nguyễn Thị Hiền Lương, Hoàng Xuân Lượng, Đào Như Mai, Dương Phạm Tường Minh, Nguyễn Văn Phó, Trần Hữu Quốc, Mai Phú Sơn, Nguyễn Xuân Thành, Vũ Duy Thắng, Trần Ích Thinh, Nguyễn Võ Thông, Vũ Thanh Thủy, Nguyễn Đắc Trung, Trương Tích Thiện, Trần Thế Truyền, Đỗ Văn Trường, Trần Minh Tú, Hoàng Văn Tùng, Lã Đức Việt, Phạm Chí Vĩnh Mục lục Tập Mục lục i Đào Huy Bích Về hoạt động khoa học ngành Cơ học Vật rắn biến dạng hai năm 2014–2015 xviii Lê Bá Anh Thermal behavior analysis of concrete box girder under the influence of temperature in Vietnam Lê Bá Anh, Trần Thế Truyền and Hồ Xuân Nam Analysis of compressive behavior of concrete with different aggregates by multiscale approaches Lê Dương Hùng Anh, Trương Quang Tri Ngô Kiều Nhi Thử nghiệm nâng cao độ xác gia công máy phay CNC lý thuyết mạng thần kinh nhân tạo 13 Le Thi Ngoc Anh, Nguyen Quang Huan, Nguyen Dinh Kien Free vibration of a tapered functionally graded sandwich beam based on the third–order shear deformation theory 21 Lê Thị Việt Anh, Nguyễn Việt Khoa, Cao Văn Mai Đào Như Mai Phân tích mỏi cho kết cấu mảnh chịu tác động tải trọng gió có kể đến ảnh hưởng dòng rối 29 Lê Thị Việt Anh, Nguyễn Việt Khoa, Cao Văn Mai Đào Như Mai Phân tích mỏi cho kết cấu khơi chịu tác động tải trọng sóng kể đến yếu tố ngẫu nhiên 37 Pham Hoang Anh and Nguyen Xuan Thanh Solving engineering optimization problems by constrained differential evolution with nearest neighbor comparison 45 Vũ Quốc Anh Nguyễn Hải Quang Tính toán dầm hình khớp dẻo chịu tải trọng động có xét đến ảnh hưởng tĩnh tải 53 Nguyễn Quốc Bảo Nghiên cứu biện pháp điều chỉnh nội lực cầu dây văng thi công đúc hẫng cân sử dụng phương pháp căng chỉnh lần căng chỉnh hai lần 60 Nguyễn Quốc Bảo, Nguyễn Văn Mỹ, Vũ Quốc Hán Lê Hoàng Bảo Đánh giá độ tin cậy cầu dây văng chịu tải trọng gió xét đến dao động flutter 68 Mục lục iii Nguyễn Thái Chung Tương tác kết cấu đường hầm san hô đảo chịu tác dụng sóng xung kích 177 Nguyễn Thái Chung, Lê Hoàng Anh Nguyễn Thanh Hưng Nghiên cứu phản ứng động công trình DKI chịu tác dụng tải trọng sóng sử dụng mô hình hệ không gian san hô làm việc đồng thời 185 Nguyễn Thái Chung Lê Hải Châu Phân tích ổn định tĩnh vỏ composite áp điện có gân gia cường 193 Nguyễn Thái Chung, Phạm Tiến Đạt Trần Văn Bình Nghiên cứu xác định vận tốc lan truyền sóng môi trường đá san hô đảo Song Tử Tây – Quần đảo Trường Sa thực nghiệm 201 Nguyễn Thái Chung Nguyễn Văn Đăng Nghiên cứu phản ứng động đường ray cong chịu tác dụng tải trọng đoàn tàu gây lý thuyết thực nghiệm 207 Nguyễn Thái Chung Hoàng Hải Nghiên cứu điều khiển dao động cho nhà cao tầng chịu tác dụng động đất thiết bị tiêu tán lượng TMD sử dụng mô hình liên hợp không gian - 215 Nguyễn Thái Chung, Nguyễn Trang Minh Lê Phạm Bình Phân tích động lực học dầm có vết nứt chịu tác dụng tải trọng di động phương pháp phần tử hữu hạn 223 Nguyễn Thái Chung Trần Trung Thành Nghiên cứu ổn định tĩnh composite áp điện có gân gia cường 231 Nguyễn Thái Chung, Lê Xuân Thùy Lê Hoàng Anh Phân tích động lực học hệ liên hợp giàn không gian – bể chứa san hô chịu tác dụng đồng thời tải trọng sóng gió 239 Nguyễn Tiến Chương Nguyễn Hải Quang Ảnh hưởng độ cứng liên kết đến phản ứng động lực khung thép có liên kết nửa cứng đàn – dẻo chịu tải trọng động 247 Nguyễn Tiến Chương Phạm Thu Hiền Ảnh hưởng độ lớn lực chảy dẻo giằng chống oằn đến làm việc khung thép chịu động đất 253 Phạm Đức Cương Nguyễn Xuân Chính Thiết kế kết cấu Bê tông cốt thép dạng khung theo độ tin cậy phần tử cho trước điều kiện Việt Nam 260 Han Van Cuong and Vu Cong Hoa Thermal performance in Electronics Control Unit 267 239 Hội nghị Khoa học toàn quốc Cơ học Vật rắn biến dạng lần thứ XII Đại học Duy Tân, TP Đà Nẵng, 7/8/2015 Phân tích động lực học hệ liên hợp giàn không gian - bể chứa san hô chịu tác dụng đồng thời tải trọng sóng gió Nguyễn Thái Chung1, Lê Xuân Thùy2, Lê Hoàng Anh3 1,2 Học viện Kỹ thuật quân sự, 236 Hoàng Quốc Việt, Bắc Từ Liêm - Hà Nội Đại học Công nghệ GTVT, 54 Triều Khúc, Thanh Xuân - Hà Nội Email: Thaichung1273@gmail.com, thuylxmta@gmail.com Tóm tắt Bài báo trình bày thuật toán phần tử hữu hạn (PTHH) phương pháp phân tích động lực học hệ liên hợp giàn thép không gian – bể chứa san hô chịu tác dụng đồng thời tải trọng sóng gió với việc sử dụng mô hình kết cấu san hô làm việc đồng thời Các kết nghiên cứu số sở khoa học cho phép lựa chọn giải pháp hợp lý công trình bố trí bãi cạn ven đảo, góp phần giảm mật độ công trình đảo, tăng khả tác chiến đảo Từ khóa: giàn không gian, bể chứa, san hô, sóng, gió, tương tác Đặt vấn đề Các đảo san hô xa bờ giữ vị trí địa trị quan trọng, có ý nghĩa to lớn việc khẳng định chủ quyền toàn vẹn lãnh thổ Quốc gia, vậy, việc nghiên cứu, đầu tư, xây dựng để sử dụng hiệu đảo điều quan trọng Trong việc nghiên cứu đó, vấn đề tận dụng bãi cạn nhằm xây dựng công trình phục vụ cho đời sống, sinh hoạt nhiều nhà khoa học quan tâm Do nhu cầu sử dụng điều kiện chật hẹp đảo nổi, ngày sử dụng bãi cạn ven đảo để giảm mật độ phân bố công trình đảo cấp thiết, theo đó, yêu cầu đưa kho tàng, bể chứa từ đảo bãi cạn Giải pháp kết cấu móng cọc giải pháp khả thi, kết cấu dạng liên hợp hệ móng cọc bể chứa, kho chứa dạng điển hình Các công trình xây dựng bãi cạn thường phải chịu tác dụng hai dạng tải trọng chủ yếu sóng gió Đây hai dạng tải trọng biến đổi ngẫu nhiên theo thời gian không gian, cộng thêm biến đổi tính đa dạng san hô ứng với vùng thềm khác nên tăng thêm mức độ phức tạp cho toán tương tác công trình – san hô chịu tác dụng tải trọng sóng gió Do vậy, với khu vực thềm san hô khác cần có nghiên cứu cụ thể toán với thông số cụ thể Sử dụng kết nghiên cứu đề tài cấp nhà nước, mã số KC.09.26/11-15 tính chất lý san hô san hô đảo Song Tử Tây – quần đảo Trường Sa, nhóm tác giả hướng đến giải toán tính tương tác công trình giàn thép không gian – bể chứa san hô chịu tác dụng tải trọng sóng gió, với mong muốn cung cấp cho nhà thiết kế thông tin tham khảo bổ ích thiết kế công trình xây dựng khu vực quần đảo Phương pháp phần tử hữu hạn sử dụng với loại phần tử phần tử không gian (mô tả kết cấu giàn) phần tử khối (mô tả san hô) Bài toán giải sở giả thiết: Vật liệu kết cấu đồng nhất, đẳng hướng đàn hồi tuyến tính, biến dạng bé; trình hệ làm việc, lớp không tách trượt tương nhau; bỏ qua ảnh hưởng rối, dòng chảy tác động môi trường Mô hình toán hình 240 Nguyễn Thái Chung, Lê Xuân Thùy, Lê Hoàng Anh Hình Mô hình toán Mô phần tử hữu hạn phương trình 2.1 Quan hệ biến dạng – chuyển vị Trường hợp tổng quát, chuyển vị theo phương trực giao x, y, z điểm vật rắn biến dạng thời điểm t tương ứng u = u(x,y,z,t), v = v(x,y,z,t), w = w(x,y,z,t) Biến dạng tỷ đối điểm xác định theo biểu thức sau (Bathe K J., 1976): x = u x xy = ,y = u y v x v y , z = , yz = w z v z w y , zx = Hay dạng ma trận: w x u z (1) = Bq , (2) [B] ma trận quan hệ biến dạng – chuyển vị, {q} véc tơ chuyển vị 2.2 Quan hệ ứng suất - biến dạng Quan hệ ứng suất – biến dạng thể công thức (Bathe K J., 1976): = D 0 , (3) đó: véc tơ biến dạng, véc tơ ứng suất, 0 véc tơ biến dạng ban đầu, [D]6×6 = [D(E,)] (4) ma trận hệ số đàn hồi 2.3 Mô hình phần tử hữu hạn 2.3.1 Phần tử chiều (hình 2) a) Các quan hệ ứng xử: Ma trận độ cứng Ke phần tử: Ke b b 1212 e T B b D Ve e e , b B b dV (5) với [Be]b, [De]b tương ứng ma trận quan hệ biến dạng – chuyển vị ma trận số đàn hồi phần tử Phân tích động lực học hệ liên hợp giàn không gian – bể chứa san hô chịu tác dụng đồng thời tải trọng sóng gió Ma trận khối lượng Me b phần tử: Me b Ne N dV , b b b 241 (6) Ve 1212 với b khối lượng riêng phần tử b) Véc tơ tải trọng sóng: Áp dụng phương trình Morison, tải trọng tác dụng lên phẩn tử theo phương x (Jamaloddin, N., Samsul, I.B., Mohammad, S.J., 2005, Khosro, B., Hosseini, S R, Mohammad, H.T., Hesam S., 2011, Ktrine, V.R., 2005): f D2 , CD D u u C1 ax (7) khối lượng riêng nước, CD C1 hệ số cản hệ số quán tính, u a x chuyển vị gia tốc y theo phương x Cường độ vận tốc hạt nước lên phương pháp tuyến: Wn u2 v2 cx u cy v , (8) thành phần vận tốc theo phương x, y, z: Unx u cx cx u cy v , Uny v cy cx u cy v , Unz cz cx u cy v , (9) cx = sincos, cy = cos, cz = sinsin, 5 cosh nky sinh nky u Gn cosn kx- t , v G n sin n kx- t , k i 1 sinh nky k i 1 sinh nky (10) với d – khoảng cách từ đáy biển đến mặt nước tĩnh, k – số sóng, - tần số góc sóng a nx a x cx cx a x cy a y , Các thành phần gia tốc theo phương x, y, z: a ny a y cy cx a x cy a y , (11) a nz cz cx a x cy a y , Áp dụng công thức Morison, ta có thành phần tải trọng sóng phân bố chiều dài sau (Jamaloddin, N., Samsul, I.B., Mohammad, S.J., 2005): f e w f x f y 0,5 CD D Wn f z Véc tơ tải trọng quy nút: Unx a nx , 2 Uny 0, 25 C1D a ny U a nz nz e Le N f P t w T e = e w dl , (12) (13) c) Véc tơ tải trọng gió: Với cấu tạo công trình biển dạng hệ nói chung, tải trọng gió tác dụng vào hai phần chủ yếu kết cấu: phần diện tích chắn gió tổng cộng phần công tác phần diện tích chắn gió thành phần Trong trường hợp tổng quát, áp lực gió tác dụng lên đơn vị diện tích chắn gió kết cấu xác định theo biểu thức sau (Kim Byoung-Wan, Kim Woon-Hak and Lee In-Won, 2002): (14) pwin t Cpair Uwin t cos , đó: pwin(t) hàm áp lực gió phân bố, Cp hệ số áp lực gió, air khối lượng riêng không khí, Uwin(t) hàm vận tốc gió theo thời gian, góc hợp Uwin t pháp tuyến mặt chắn gió 242 Nguyễn Thái Chung, Lê Xuân Thùy, Lê Hoàng Anh Tương tự, lực gió phân bố theo chiều dài thanh: qwin t BCpair Uwin t 2 coswin , (15) B bề rộng mặt cắt ngang thanh, trường hợp hình trụ B = D (đường kính mặt cắt ngang thanh), win góc lệch trục với trục y hệ toạ độ tổng thể Trường hợp vận tốc gió không đổi Uwin(t) = U0, gió thổi theo phương y, biểu thức áp lực gió (14) (15) trở thành: pwin t p0win Cpair U0 (16) qwin t q0win BCpair U0 (17) Và lúc gió tác dụng lên kết cấu tương tự tải trọng xung tức thời, theo thống kê mức độ nguy hiểm kết cấu trường hợp không trường hợp vận tốc gió hàm thời gian Véc tơ tải trọng nút phần tử áp lực gió gây nên xác định theo biểu thức sau (Nguyễn Thái Chung, Nguyễn Văn Chình, 2013, Bathe K J., 1976): P t win e Le T = N b pwin t dl 0 (18) 2.3.2 Phần tử khối nút Phần tử dùng mô tả san hô, xét phần tử khối lục diện nút đẳng tham số, nút có bậc tự chuyển vị dài theo phương x, y z hệ trục toạ độ tổng thể Hình mô tả phần tử hệ tọa độ tổng thể hệ tọa độ cục sau chuẩn hóa đơn vị v2 v1 θx2 θz2 y x u2 θx1 θy2 w2 θz1 w1 u1 z θy1 a, Trong hệ toạ độ chung b, Trong hệ toạ độ cục Hình Phần tử lục diện điểm nút Hình Phần tử chiều Hình dạng hình học phần tử cho (Bathe K J., 1976): 8 x i 1 Ni x i , y i 1 Ni yi , z Ni zi , (19) i 1 x1 y 1 z1 x Viết lại (19) dạng ma trận: y N s x , z z8 (20) Phân tích động lực học hệ liên hợp giàn không gian – bể chứa san hô chịu tác dụng đồng thời tải trọng sóng gió 243 [N]s ma trận hàm dạng phần tử khối nút, xi, yi zi toạ độ nút i hệ toạ độ vuông góc tổng thể Các hàm chuyển vị theo phương x, y z hệ trục tổng thể biểu diễn: u N u , 8 v i i i 1 N v , w i i N w , i (21) i i 1 i 1 với ui, vi wi bậc tự nút i Ma trận độ cứng phần tử: Ke s 2424 B D B dV , T s s (22) s Ve [B]s, [D]s tương ứng ma trận quan hệ biến dạng – chuyển vị ma trận số đàn hồi phần tử khối điểm nút Ma trận khối lượng [M]e phần tử xác định sau: Me s N N dV , (23) s T s s Ve 2424 với s khối lượng riêng phần tử 2.4 Phương trình động lực học Theo phương pháp PTHH, tác dụng tải trọng động, phương trình vi phân mô tả dao động phần tử hệ – phần tử khối sau: M e s Me q C C q K K q P (t) P (t) P (t) , e b e e s e e b e s e e w e win e (24) b Me , Me , Ce , Ce , K e , K e ma trận khối lượng, ma trận cản ma s b s b s b trận độ cứng phần tử khối (mô tả san hô) phần tử 2.5 Thuật toán phần tử hữu hạn giải toán Sau ghép nối ma trận, véc tơ tải trọng phần tử thành ma trận, véc tơ tải trọng tổng thể khử biên, thu phương trình vi phân mô tả dao động hệ sau: Mq Cq Kq P(t) , (25) với M Me Me , K K e K e , C M K , s b s b e e hệ số cản Rayleigh (Bathe K J., 1976) Pháp độ cứng trực tiếp, ma trận số sơ đồ Skyline sử dụng để tập hợp ma trận tổng thể Giải phương trình vi phân tuyến tính (25), tác giả sử dụng phương pháp tích phân trực tiếp Newmark lập trình tính môi trường Ansys Chương trình tính có tên 3D_Offshore_Structures, kết cấu mô tả dạng phần tử: PIPE16, BEAM188, SOLID45 Khảo sát số 3.1 Bài toán Thông số kết cấu: Kích thước hình bao khối chân đế giàn thép không gian: a x a x h = 1,5mx1,5mx4,0m, cọc thép ống: Ф50x3,0mm, giằng có kích thước Ф42x3,0mm Kích thước sàn công tác: b x b = 2,5x2,5m, sàn làm từ thép hộp, tiết diện ngang: 244 Nguyễn Thái Chung, Lê Xuân Thùy, Lê Hoàng Anh 50x50x3mm Vật liệu giàn thép, có mô đun đàn hồi E = 2,11011N/m2, hệ số Poisson = 0,3, khối lượng riêng = 7800kg/m3 Chọn sai số lặp d = 0,5%, xác định kích thước miền nghiên cứu H×B1×L1 = 5,5m×8,0m×8,0m Nền gồm lớp, đánh số 1, 2, 3, 4, với đặc trưng lý bảng Lớp Bảng Đặc trưng vật liệu san hô Ef (N/cm2) f f (kg/m3) 0,17 0,20105 1,28103 0,19 0,2510 1,29103 0,24 0,31105 2,16103 0,29 0,4310 2,7103 Bề dày lớp (m) 0,8 1,2 1,5 2,0 Tỷ số cản 0,05 Thông số tải trọng: sóng biển có chiều cao Hw = 3,56m, độ sâu nước dw = 1,5m, khối lượng riêng nước w = 1000kg/m3, chu kỳ sóng Tw = 7,83s, hệ số lực cản CD = 0,5, hệ số quán tính C1 = 2,0, hệ số áp lực gió Cp = 1, khối lượng riêng không khí air = 1,225kg/m3 Khối lượng tổng cộng bể chứa quy đổi P = Giản đồ vận tốc gió U win t hình (Kim Byoung-Wan, Kim Woon-Hak and Lee In-Won, 2002) P b Uwin(t) sóng biển GIAN DO VAN TOC GIO THEO THOI GIAN 50 a H1 45 Van toc Uw in(t) [m/s] 40 H2 H 35 H3 30 H4 25 B1 20 15 10 20 30 40 50 60 70 80 90 100 L1 Thoi gian t [s] Hình Giản đồ vận tốc gió Uwin(t) với Umax = 46,35m/s Hình Mô hình toán Với chương trình lập, tính toán với hai mô hình: xét đồng thời kết cấu – nền, xét riêng kết cấu Kết so sánh cho trường hợp chuyển vị ngang đỉnh giàn mô men uốn chân cọc, thể bảng hình 6, Bảng Giá trị lớn chuyển vị đỉnh giàn mô men uốn chân cọc Mô hình tính max Ux [cm] max M y [Nm] Xét riêng kết cấu 7,04118 4861,01 Xét kết cấu – 7,52944 3982,20 Sai số [%] 6,9 - 18,1 Phân tích động lực học hệ liên hợp giàn không gian – bể chứa san hô chịu tác dụng đồng thời tải trọng sóng gió 3000 MH tuong tac MH khong tuong tac MH tuong tac MH khong tuong tac Momen uon mat cat chan cot [Nm] 2000 Chuyen vi ngang U [cm] 245 -2 -4 -6 1000 -1000 -2000 -3000 -4000 -8 10 Thoi gian t[s] 15 -5000 20 Hình Đáp ứng chuyển vị ngang Ux theo thời gian t đỉnh giàn 10 Thoi gian t[s] 15 20 Hình Đáp ứng mô men uốn My theo thời gian t mặt cắt chân cọc 3.2 Khảo sát số yếu tố đến đáp ứng động giàn thép 3.2.1 Ảnh hưởng tải trọng bể chứa: Đánh giá mức độ ảnh hưởng tải trọng bể chứa đến dao động hệ giàn không gian – bể chứa, tác giả khảo sát đáp ứng động hệ với giá trị P thay đổi: P1 = 500kg, P2 = 1200kg, P3 = 2000kg, P4 = 4000kg, P5 = 6000kg Giá trị cực trị chuyển vị ngang đỉnh mômen uốn chân cọc thể bảng Bảng Giá trị lớn chuyển vị đỉnh giàn mô men uốn chân cọc ứng với Pi P [kg] 500 1200 2000 4000 6000 max [cm] 7,52941 7,52942 7,52944 7,52946 7,52949 max M y [Nm] 3982,24 3982,23 3982,20 3982,11 3982,02 Ux Đáp ứng cụ thể hình 8, 3982.3 Momen uon chan cot lon nhat Mymax [Nm] Chuyen vi ngang lon nhat U max [cm] 7.5295 7.5295 7.5294 7.5294 7.5293 1000 2000 3000 4000 Tai be chua P [kg] 5000 6000 3982.25 3982.2 3982.15 3982.1 3982.05 3982 1000 2000 3000 4000 Tai be chua P [kg] 5000 6000 Hình Đáp ứng chuyển vị ngang U max theo x Hình Đáp ứng mô men uốn M max theo y tải trọng P bể chứa tải trọng P bể chứa 3.2.2 Ảnh hưởng vận tốc gió: Đánh giá ảnh hưởng tải trọng gió đến đáp ứng động hệ, tác giả khảo sát thay đổi vận tốc gió so với giản đồ vận tốc gió ban đầu, biến đổi từ 25% đến 125% Các giá trị thay đổi chuyển vị ngang lớn đỉnh giàn mômen uốn lớn chân cột thể bảng hình 8, 246 Nguyễn Thái Chung, Lê Xuân Thùy, Lê Hoàng Anh max Bảng Giá trị U x Uwin(t) [%] (so với BTCB) max Ux [cm] max M y [Nm] max đỉnh giàn mô men M y 25 50 75 100 125 0,463 1,877 4,232 7,529 11,769 241,966 986,149 2234,49 3982,20 6229,19 8000 Momen uon chan cot lon nhat Mymax [Nm] Chuyen vi ngang lon nhat U max [cm] 16 14 12 10 20 chân cọc thay đổi giản đồ vận tốc 40 60 80 100 Ty le [%] gio so voi gian van toc BTCB 120 Hình Chuyển vị ngang U max theo Uwin(t) x 7000 6000 5000 4000 3000 2000 1000 20 40 60 80 100 Ty le [%] gio so voi gian van toc BTCB 120 Hình Mô men uốn M max theo Uwin(t) y Kết luận Bài báo đạt nội dung sau: Xây dựng thuật toán chương trình phân tích động lực kết cấu giàn không gian – bể chứa san hô chịu tải trọng sóng gió Phân tích, đánh giá ảnh hưởng mô hình tính số yếu tố khác đến đáp ứng hệ, đưa nhận xét, cụ thể: Đánh giá theo mô hình tương tác kết cấu – nền, hệ mềm so với xét riêng kết cấu (chuyển vị điểm tính kết cấu tăng mô men uốn chân cọc giảm đi) Với thông số kết cấu mô hình tính cho, dải tải trọng bể chứa khảo sát ảnh hưởng đến chuyển vị ngang mô men uốn điểm tính kết cấu Khi tăng tải trọng bể chứa, chuyển vị ngang lớn tăng, mô men uốn chân cọc lớn giảm, lượng tăng, giảm không đáng kể Tăng vận tốc gió, chuyển vị nội lực kết cấu tăng phi tuyến Lời cảm ơn: Các tác giả trân trọng cảm ơn Ban chủ nhiệm đề tài NCKH cấp Nhà nước, mã số KC.09.26/11-15 tạo điều kiện thuận lợi hỗ trợ để tác giả thực nội dung nghiên cứu Tài liệu tham khảo Nguyễn Thái Chung, Nguyễn Văn Chình, 2013, Ảnh hưởng số yếu tố đến phản ứng động công trình biển tác dụng tải trọng sóng gió, Tạp chí Khoa học Công nghệ biển, Tập 13, Số 2; 2013, tr.33-40 Bathe K J (1976), Numerical Methods in Finite Element Analysis, Prentice - Hall, Inc., New Jersey Jamaloddin, N., Samsul, I.B., Mohammad, S.J., 2005, Waleed Abdul Malik Thanoon1 and Shahrin Mohammad, Simulation of wave and current forces on template offshore structures, Suranaree J Sci Technol 12(3), pp 193-210 Khosro, B., Hosseini, S R, Mohammad, H.T., Hesam S., 2011, Seismic Response of a Typical Fixed Jacket-Type Offshore Platform (SPD1) under Sea Waves, Open Journal of Marine Science, 2011, 1, pp 36-42, doi:10.4236/ojms.2011.12004 Published Online July 2011 (http://www.Sci-RP.org/journal/ojms) Kim Byoung-Wan, Kim Woon-Hak and Lee In-Won, 2002, Three-dimensional Plate Analyses of Wind loaded Structures, Department of Civil Engineering, Korea Advanced Institute of Science and Technology, 3731 Guseong-dong, Yuseong -gu, Daejeon, 305-701, Korea Ktrine, V.R., 2005, Dynamic behaviour of jackets esposed to wave-in-deck forces, Doctor of Philosophy, University of Stavanger, Norway