Tính toán đáp ứng động lực học của kết cấu

Chương 3 TÍNH TOÁN ĐÁP ỨNG ĐỘNG LỰC HỌC CỦA KẾT CẤU 3.1 .Cọc đơn và kết cấu 1 chiều : Trong phần này đưa ra phương pháp tìm tần số tự nhiên và phân bố khối lượng ảnh hưởng cho cọc đơn có hoặc không có khối lượng ở đầu bị khống chế 3.1.1 . Cọc đơn có hoặc không có khối lượng ở đầu : Bươc 1 : Xác định các đặc trưng hình học vàvật lý của cọc Hình 3.1 Mô hình cọc đơn có khối lượng đầu cọc l : Chiều dài từ đáy mũ đến đất l’ : Chiều dài tính toán l’ = l + la/cosa la : Chiều dài ngàm trong đất d : Chiều sâu nước d’ = d + la a - Góc hợp bỡi phương đứng và trục của cọc . la = 3.5 ¸ 4.5D (Đối với đất sét cứng ) la = 7 ¸ 8.5D ( Đối với phù sa mềm ) Thường lấy la = 6D Tỉ số khối lượng : m = (mw-m)/m Tỉ số chiều sâu : h = d’/l’cosa Bước 2 : Hệ số khối lượng tương đương : b Tra bảng hình 6.2 b= f(h, m) Khối lượng phân bố tương đương Bước 3 : Xác định tần số tự nhiên N Hz (*) Ta có sơ đồ tương đương Hình 3.2 Mô hình cọc tương đương 3.1.2 Cocï đơn bị khống chế bỡi đài : Hình 3.3 - Cọc đơn bị khống chế ở đầu cọc Bước 1 : Xác định các đặc trưng hình học và vật lý của cọc I : Mômen quán tính của cọc l :Chiều dài tính toán E : Modun đàn hồi của cọc m : Khối lượng của cọc trên 1 đơn vị chiều dài trong không khí mw : Khối lượng của cọc trên 1 đơn vị chiều dài trong nước m : Tỉ số khối lượng h : Tỉ số độ sâu Tỉ số độ cứng ngang Tỉ số độ cứng xoay Kx , Kf tùy thuộc vào kết cấu Bước 3 : Khối lượng phân bố tương đương b - được tính từ hình 6.5 phụ thuộc vào h và Cx Bước 4 : Tần số tự nhiên Hz (*) Ac – Tra hình 6.6 phụ thuộc vào Cx ,Cf Bước 5 : Khối lượng phân bố hiệu quả Hình 3.4 – Mô hình cọc ngàm tương đương (*) Những công thức trên được trích từ sách DYNAMIC OF MARINE STRUCTRES –AIT – Asian instite of Technology –M.G.Hallam BScPhD,N.J.Heaf Beng PhD, L.R. Wootton BScPhD MICE MRAeS- 1977 3.1.3.Cọc đơn có khối lượng tập trung và bị khống chế bởi đài : Hình 3.5 – Cọc có khối lượng ở đầu cọc và bị khống chế Hình 3.6 Mô hình cọc lý tưởng Bước 1 : Xác định các đặc trưng hình học và vật lý của cọc I : Mômen quán tính của cọc l :Chiều dài tính toán E : Modun đàn hồi của cọc m : Khối lượng của cọc trên 1 đơn vị chiều dài trong không khí mw : Khối lượng của cọc trên 1 đơn vị chiều dài trong nước m : Tỉ số khối lượng h : Tỉ số độ sâu Bước 2 : Xác định khối lượng M1 , M2 mw M1 M2 m d’ l’ Nếu h>0.5 ® M1 = ml’(1-h)2+mwl’(2h-h2-0.5) M2 = ml’h (1-h)+mwl’(h-0.5)2 + M Nếu h<0.5 ® M1 = ml’(0.5-h2) + mwl’h2 M2 = ml’/4 + M Bước 3 : Xác định độ cứng không thứ nguyên Bước 4 : Xác định các thông số độ mền không thứ nguyên Bước 5 : Xác định w1,2 và N1,2 (Tầng số dao động và tầng số riêng ) Bước 6 : 1) Xác định Mode Shape (dao động) Với n=1,2 Hình 3.7 – Dạng dao động (Mode Shapes) 2 ) Xác định khối lượng hiệu quả n = 1,2 Với 3.2. Kết cấu phẳng và khung : Trong phần này phân tích các loại đáp ứng động để tính toán tần số dao động tự nhiên và khối lượng hiệu quả của kết cấu . Những đại lượng này được sủ dụng để đánh giá độ nhạy của kết cấu khi chịu tác dụng của tải trọng động . 3.2.1 . Khung dao động trong mặt phẳng : Có hai dạng : - Dạng lắc lư như một thể thống nhất . - Các cọc trong kết cấu dao động theo kiểu cung tên Hầu hết hai dạng dao động này đồng thời xảy ra nên không thể xét độc lập 1 . Phương pháp tính tay : Hình 3.8 – Phân tích khung phẳng thành các cọc thành phần Bước 1 : Xác định các đặc trưng hình học và vật lý của cọc và đài (dầm) (l’)j –Chiều dài tính toán của cọc và chiều sâu của nước (d’)j cho mỗi cọc ,Với j là số lượng cọc Xác định đặc trưng mặt cắt ngang Aj của cọc và Ab cho đài cọc hoặc dầm . Tính tỷ số chiều sâu hj cho mỗi cọc, với hj = d’j/ l’j (cosa)j Bước 2 : Tính khối lượng thu gọn đặc trưng của khung . Tính khối lượng phân bố (m)j trong không khí và (mw)j trong nước cho mỗi cọc . Tính khối lượng thu gọn M1,M2 cho mỗi cọc : h j>0.5 (M1)j = (m)j(l’)j(1-(h)j)2 + (mw)j(l’)j(2(h)j – (h)j – 0.5) (M2)j = (m)j(l’)j(hj)j(1-(h)j) + (mw)j(l’)j( (h)j – 0.5)2 Khi h j< 0.5 (M1)j = (m)j(l’)j(0.5-(h)2j) + (mw)j(l’)j( (h)2j (M2)j = (m)j(l’)j /4 Vậy : Trong đó : M là khối lượng của dầm và tải trọng trên dầm Bước 3 : Xác định Kx và Kf ( độ cứng chống ngang và chống xoay) Độ cứng chống ngang Ax : Kf Xác định từ các công thức của hình 6.14 (Phần phụ lục) 3.2.2 Khung dao động ngoài mặt phẳng : Bước 1 : Xác định chiều dài tính toán cho mỗi cọc : (l’)j Bước 2 : Tính mj cho mỗi cọc Khối lượng quy về nút cho mỗi cọc : M1j = mjl’j / 4 M2j = mjl’j / 4 + M Bước 3 : Xác định tổng khối lượng M và momen quán tính quya của khối lượng đối với phương z . IM = ML2b/12 Với Lb – Chiều dài của dầm Bước 4 : Tính độ cứng chống ngang Kyj và độ cứng chống xoay Kf j Trong đó : D1= D-2d d - chiều dày cọc Xj : Toạ độ cọc thứ j Bước 5 : Tính các tần số : Bước 6 : Xác định tần số riêng : (rad/sec) Dạng dao động Bước 7 :Xác định khối lượng hiệu quả với tần số riêng Chuyển vị đầu cọc : yj = Y+ xj q Dùng mode shape (Y=1) Þ yj = 1 + xjq Chuyển vị cọc theo độ sâu : Trong đó aj – góc nghiêng của cọc 3.3 . Trụ cập tàu và kết cấu 3 chiều : 3.3.1 . Phân tích tần số 3 chiều bằng tay : Phương pháp này đưa ra phương pháp tính tần số riêng và khối lượng hiệu dụng theo kiểu lắc lư cho những công trình bao gồm cọc đứng, cọc xiên và đài cứng . Phân tích 3 chiều (x,y,z) : dao động lắc lư x,y và xoắn f Phương pháp đơn giản để tính toán tần số riêng và các dạng dao động (Mode shape ) như sau : Bước 1 : Tính khối lượng tương đương m cho mỗi cọc . Tính tổng khối lượng tương đương cho mỗi cọc theo : Mj = mj l’j Bước 2 : Khối lượng tương đương tại đỉnh mỗi cọc là ¼ Mj Bước 3 : Tính tổng khối lượng, vị trí tâm C và mômen uốn trong mặt phẳng chứa C của đài và của cả khối lượng qui về đầu cọc Mj Khi tính toán, định nghĩa hệ trục OXY trùng với mặt đài và trục Z đi qua tâm C . Bước 4 : Tính độ cứng của cọc Tính độ cứng ngang của cọc Ba dạng cọc được xem xét là : Cọc đứng . Cọc nghiên trong mặt phẳng OXZ 3) Cọc nghiên trong mặt phẳng OYZ Bước 5 : Xác định tâm cứng của kết cấu . X y Kx Kf Ky Dao động khối tâm Tâm độ cứng Y y q Toạ độ khối tâm trước khi dao động Hình 3.9 – Chuyển vị của khối tâm Gọi toạ độ tâm cứng là Lấy mômen theo đường Trong đó : - Độ cứng theo phương x YI : Tọa độ y của cọc thứ i n : Tổng số cọc Tương tự theo đường Bước 6 : Xác định các giá trị độ cứng Kx,Ky,Kf. Trong đó –(kq)I : độ cứng chống xoắn cọc thứ i Bước 7 : Xác định phương trình chuyển động theo phương x,phương y và góc xoay q : Và M - Tổng khối lượng của kết cấu Im – Momen quán tính quay của khối lượng (momen quán tính khối tâm ) Bước 8 : Giải phương trình chuyển động để tìm tần số tự nhiên dao động (Mode shapes) Ma trận viết cho tâm cứng như sau Ma trận chuyển vị a11 a12 a13 a14 a15 a16 x a22 a23 a24 a25 a26 y a33 a34 a35 a36 z a44 a45 a46 qx Đối xứng a55 a56 qy a66 qz Từ ma trận (6x6) trên có thể biến đổi thành ma trận (2x2) Với (a*11-Mw2n) a*12 a*13 a*12 (a*22 – Mw2n) a*23 a*13 a*23 (a*33-Iqw2n) Để tìm tần số dao động w1,w2,w3 ta giả hệ phương trình : Det 3.3.2 Móng cọc đài cao cứng – sơ đồ không gian : yi z 90o Ii ji IIi Oi y x ji IIIi O xi 90o Cọc i Để tính toán móng cọc đài cao cứng có hai phương pháp đó là theo phương pháp của cơ học kết cấu và theo phương pháp ma trận . Trong luận án này chỉ đề cập phương pháp ma trận và tính toán móng cọc đài cao cứng theo sơ đồ không gian theo Spiro .Cho mặt bằng móng cọc như sau : Mặt bằng bố trí móng cọc a b c a g QIIi MIi Ii QIIIi IIi MIIIi Ni MIIi b Hình 3.10 – Các qui ước về chuyển vị và nội lực Trình tự các bước giải : Bước 1 : Xác định tọa độ và góc nghiêng so với phương đứng của từng cọc trong đài Bước 2 : Cách xác định góc l Ví dụ móng có 6 cọc như hình vẽ các góc l được y x 300 6 4 3 2 1 5 xác định như sau : l1 = 0 l2 = 270o l3 = 180o l4 = 90o l5 = 0 l6 = 330o Hình 3.11 – Cách xác định góc l Bước 3 : Xác định ma trận độ cứng Ma trận độ cứng của hệ có thể thành lập như sau : Đặt : ci = cos(aI) cci = cos(lI) si = sin(aI) ssi = sin(aI) cai = -ci. cci cxi = si. cci cbi = - ci . ssi cyi = si . ssi Với a - là góc nghiên của cọc so với phương đứng Trong đó : a1(i) - Ma trận chuyển vị đầu cọc thứ i a2(i) – ma trận biến đổi giữa trục địa phương và toàn cục a1(i)T- ma trận chuyển trí của ma trận a1(i) a2(i)T- ma trận chuyển trí của ma trận a2(i) Đặt Các hệ số của ma trận cho tâm cứng được viết như sau : Ta có : 3.3.2 . Khối lượng hiệu quả liên quan đến tần số riêng (Effective- Mass) : Trình tự xác định và như sau : Bước 1 : Tìm sự chuyển vị của tâm C liên quan đến dạng dao động, bằng cách giải 3 phương trình sau : (a*11 – Mw2n)x + a*12y + a*13q = 0 a*12x + (a*22 – Mw2n)y + a*23q = 0 a*13x + a*23y + (a*33 – Iw2n) q = 0 Từ hệ phương trình trên có thể viết dưới dạng ma trận sau : = Biến đổi ma trận thành ½ tam giác trên : = với q = 1 Þ Bước 2 : Xác định chuyển vị đầu cọc : xj = x - Yjq yj = y + Xjq q = 1 Với : Xj , Yj – là toạ độ đỉnh cọc Bước 3 : Chuyển vị cọc theo độ sâu : Bước 4 : Xác định khối lượng hiệu quả : Với : d- độ sâu nước a- góc nghiêng