Bối cảnh và vấn đề nghiên cứu
Công trình cầu là bộ phận huyết mạch trong mạng lưới hạ tầng giao thông vận tải. Khi xảy ra động đất, các dao động nền đất truyền qua kết cấu mố trụ và móng có thể gây biến dạng nghiêm trọng, nứt vỡ, xô lệch kết cấu nhịp hoặc dẫn đến sụp đổ cục bộ đến toàn phần. Tại Việt Nam, nguy cơ phát sinh động đất đã được chứng minh qua nhiều ghi nhận lịch sử và khảo sát địa vật lý. Tuy nhiên, trong thực tiễn tư vấn thiết kế và quản lý xây dựng công trình giao thông tại thời điểm nghiên cứu, hệ thống quy chuẩn kỹ thuật kháng chấn cho ngành cầu đường còn thiếu tính đồng bộ.
Bối cảnh thực tiễn cho thấy các cơ quan trong nước chủ yếu áp dụng Tiêu chuẩn ngành 22 TCN 221-95 do Bộ Giao thông Vận tải ban hành năm 1995. Đồng thời, các dự án hạ tầng sử dụng vốn vay quốc tế lại áp dụng các tiêu chuẩn nước ngoài như Tiêu chuẩn thiết kế cầu của Hoa Kỳ (AASHTO-1996, AASHTO-1998) hoặc các quy trình của Pháp, Úc, Nhật Bản. Thực tế này nảy sinh nhiều bất cập:
- Việc áp dụng đa dạng các tiêu chuẩn ngoại lai nhưng chưa nắm vững bản chất động lực học và mức độ an toàn chênh lệch giữa các quy trình dẫn đến khó khăn trong công tác thẩm tra, thiết kế và quản lý kỹ thuật.
- Các phần mềm chuyên dụng tính toán kháng chấn đi kèm các dự án quốc tế có chi phí bản quyền rất cao, rào cản ngôn ngữ và hệ điều hành phức tạp (chẳng hạn các phần mềm của Nhật Bản), hầu như không được chuyển giao lại cho kỹ sư Việt Nam.
- Việc tính toán dao động động lực của kết cấu cầu theo Tiêu chuẩn ngành 22 TCN 221-95 đòi hỏi khối lượng xử lý toán học lớn đối với chuyển vị của các chất điểm, dễ xảy ra sai số nếu không có công cụ tính toán tự động hóa hỗ trợ.
Xuất phát từ các vấn đề trên, nghiên cứu xác lập các mục tiêu và nhiệm vụ cụ thể sau:
- Giới thiệu một số đặc trưng cơ bản của hiện tượng động đất, đánh giá nguy cơ xảy ra động đất trên thế giới và tại lãnh thổ Việt Nam.
- Trình bày hệ thống phương pháp lý thuyết tính toán động đất áp dụng cho công trình cầu, bao gồm các phương pháp phân tích động đất trong các quy trình hiện hành tại Việt Nam.
- Dựa trên các chương trình máy tính sẵn có, xây dựng các phần mềm phụ trợ nhằm đơn giản hóa quy trình tính toán kháng chấn và thực hành tính toán trên các trường hợp cụ thể.
- Nghiên cứu so sánh các nguyên lý cơ bản, phương pháp thực hành phân tích, sơ bộ đánh giá mức độ an toàn và nội lực kháng chấn của kết cấu cầu theo Tiêu chuẩn ngành 22 TCN 221-95, Tiêu chuẩn thiết kế cầu 22 TCN 272-01 của Bộ Giao thông Vận tải và Tiêu chuẩn AASHTO-1996 của Hoa Kỳ.
- Đưa ra các khuyến nghị kỹ thuật về việc áp dụng các phần mềm máy tính sẵn có trong tính toán động đất cho kết cấu cầu.
- Đưa ra các ý kiến đánh giá, định hướng áp dụng thống nhất các tiêu chuẩn thiết kế kháng chấn cho công trình cầu tại Việt Nam.
Về đối tượng và phạm vi nghiên cứu, đề tài tập trung vào công trình cầu chịu tác động của tải trọng động đất. Về mặt vật lý, nghiên cứu giới hạn tính chất vật liệu tuân theo các định luật đàn hồi cổ điển, khảo sát các chuyển vị nhỏ quanh vị trí cân bằng ổn định. Do những hạn chế trong bản thân tiêu chuẩn 22 TCN 221-95, phạm vi so sánh chỉ dừng lại ở việc đánh giá các nguyên lý cơ bản, phương pháp thực hành và sơ bộ đánh giá kết quả nội lực kháng chấn cho một số dạng kết cấu cầu cơ bản khi sử dụng AASHTO-1996, 22 TCN 272-01 so với 22 TCN 221-95.
Cơ sở lý thuyết và phương pháp
Khung lý thuyết và mô hình tính toán
Nghiên cứu được xây dựng trên nền tảng cơ học kết cấu, động lực học công trình và lý thuyết địa chấn công trình:
- Thuyết kiến tạo mảng (Plate Tectonics): Giải thích nguồn gốc sinh chấn từ sự chuyển dịch của 11 vĩ mảng thạch quyển (dày 70 km ở đại dương, 140 km ở lục địa) trượt trên lớp dung nham dẻo, gồm 5 dạng chuyển động cơ bản: phân ly, dũi ngầm, trườn, va chạm đàn hồi và rúc đồng quy.
- Lý thuyết dao động đàn hồi kết cấu: Khảo sát dao động của hệ cơ học đàn hồi dưới kích động nền đất với các phương trình vi phân chuyển động có cản. Phân tích kết cấu theo mô hình hệ một bậc tự do (SDOF) và hệ nhiều bậc tự do (MDOF).
- Lý thuyết phổ phản ứng (Response Spectrum Method): Khái niệm do M. Bio khởi xướng năm 1933 và được Cotrinski hoàn thiện. Phương pháp dựa trên phản ứng cực đại (chuyển vị, vận tốc, gia tốc) của hệ một bậc tự do chịu kích động động đất để phân tích cho hệ kết cấu phức tạp.
- Nguyên lý tương đương năng lượng: Chuyển đổi hệ kết cấu nhiều bậc tự do ($n$ bậc) về $n$ hệ tương đương một bậc tự do dựa trên điều kiện bảo toàn chu kỳ dao động riêng ($T_i$), lực cắt đáy ($F_i$), mômen uốn đáy ($M_i$) và động năng ($E_i$).
- Hệ thống tiêu chuẩn và thang đo quy ước: Tiêu chuẩn 22 TCN 221-95 (Bộ GTVT ban hành theo Quyết định số 3008 QĐ/KH-KT ngày 30/5/1995), Tiêu chuẩn thiết kế cầu 22 TCN 272-01, Tiêu chuẩn AASHTO-1996; thang chấn cấp MSK-64 (Medvedev-Sponheuer-Karnik), thang EMS-1992, thang Mercalli cải biên (MM), thang JMA (Nhật Bản) và thang Rossi-Forel.
Phương pháp nghiên cứu thực tế
Tác giả triển khai kết hợp các phương pháp giải tích kỹ thuật, phương pháp phân tích số và thống kê địa chấn:
- Phương pháp tính toán tĩnh tương đương (Static Analysis): Dựa trên lý thuyết của Ohmori và Sano, giả định công trình là vật rắn tuyệt đối chịu lực quán tính phân bố tỉ lệ với gia tốc nền thông qua hệ số động đất $K_s$.
- Phương pháp phổ ứng xử (Response Spectrum Analysis): Sử dụng phương pháp ma trận độ cứng, giải bài toán trị riêng để xác định các tần số dao động riêng ($\omega_i$) và các véc-tơ dạng dao động riêng (${\phi}_i$). Tách các phương trình vi phân liên kết thành các phương trình độc lập theo tọa độ chính nhờ tính trực giao của các dạng dao động, xác định chuyển vị và lực quán tính cực đại của từng mode dao động, sau đó tổ hợp phản ứng toàn phần (theo công thức Merchant hoặc R. Clough).
- Phương pháp phân tích theo lịch sử thời gian (Time-History Analysis): Áp dụng phương pháp sai phân từng bước để tích phân trực tiếp phương trình vi phân chuyển động theo thời gian đối với các kết cấu có khối lượng tập trung ở đỉnh (như trụ cao, tháp cầu).
- Phương pháp thống kê và giải tích xác suất: Xử lý chuỗi số liệu động đất lịch sử nhằm thiết lập quy luật lặp lại động đất $\lg N^* = a - b M_S$ và lập bản đồ phân vùng chấn động lãnh thổ.
Nguồn dữ liệu sử dụng
- Số liệu địa chấn lịch sử, điều tra thực địa và quan sát máy móc tại Việt Nam từ năm 1900 đến 1995 (với hơn 700 trận động đất có $M_S < 7.0$).
- Đề tài nghiên cứu khoa học cấp Nhà nước: "Cơ sở dữ liệu cho các giải pháp giảm nhẹ hậu quả động đất ở Việt Nam" (bản đồ nguồn phát sinh động đất, bản đồ tần suất chấn động $I \ge 6, 7, 8$ theo thang MSK-64, xác suất xuất hiện $10%$ trong các chu kỳ 20, 50 và 100 năm).
- Băng ghi gia tốc thực tế của các trận động đất mạnh trên thế giới, tiêu biểu là số liệu trận động đất El Centro (1940).
- Các công thức thực nghiệm suy giảm gia tốc nền toàn cầu và khu vực: công thức Donovan (1973), Cornell, Esteva (1974), Shebalin (1968), Blake và công thức dự báo $M_{\max}$ của Nguyễn Đình Xuyên.
Nội dung chính theo từng chương
Chương 1: Khái niệm về động đất
Chương này tập trung trình bày bản chất vật lý, cơ chế phát sinh, các đại lượng đặc trưng của động đất, đánh giá độ nguy hiểm động đất trên thế giới và phân vùng địa chấn tại Việt Nam.
1. Cơ chế sinh chấn và phân loại chuyển động:
Động đất phát sinh do 4 nguyên nhân chính: va chạm thiên thạch, thử nổ hạt nhân ngầm, sập hang ngầm và vận động kiến tạo. Trong đó, vận động kiến tạo chiếm $95%$ các trận động đất trên thế giới. Dưới tác động trượt tương đối giữa các mảng thạch quyển (đặc biệt là chuyển động dũi ngầm và trườn), biến dạng thế năng tích lũy dọc theo các đới đứt gãy tới hạn sẽ giải phóng đột ngột thành động năng dưới dạng sóng địa chấn. Ba loại sóng chính gồm: sóng dọc ($P$, vận tốc lan truyền $V_P = 7 \div 8\text{ km/s}$), sóng ngang ($S$, sóng cắt, $V_S = 4 \div 5\text{ km/s}$) và sóng mặt (vận tốc $1{,}5 \div 15\text{ km/s}$ trên nền cứng, $0{,}5 \div 1{,}5\text{ km/s}$ trên nền yếu). Thời gian chênh lệch tiếp nhận tín hiệu giữa sóng $P$ và sóng $S$ ($T_{PS}$) là cơ sở xác định tọa độ chấn tiêu và chấn tâm.
2. Phân biệt Magnitude và cấp động đất:
Văn bản phân định ranh giới khoa học giữa hai khái niệm:
- Magnitude ($M$): Đại lượng biểu thị mức năng lượng giải phóng dưới dạng sóng đàn hồi (độ Richter), gồm: Magnitude địa phương ($M_L$, thiết lập bởi Richter năm 1935 trên máy Wood-Anderson cho khoảng cách chấn tâm $\Delta \le 800\text{ km}$), Magnitude theo sóng mặt ($M_S$, áp dụng cho $\Delta = 200 \div 1600\text{ km}$, độ sâu chấn tiêu $h < 50\text{ km}$ theo chuẩn Zurich 1967) và Magnitude theo sóng khối ($M_B$). Năng lượng động đất liên hệ với Magnitude qua phương trình $\lg E = 11{,}5 + 1{,}5 M_S$ (trong đó $E$ tính bằng erg). Một trận động đất $M_L = 7{,}3$ tương đương năng lượng giải phóng từ vụ nổ 50 triệu tấn thuốc nổ TNT (50 megaton).
- Cấp động đất ($I$): Đại lượng quy ước chỉ cường độ chấn động và mức độ phá hoại trên mặt đất. Thang MSK-64 phân chia thành 12 cấp, đánh giá mức độ hư hại từ cấp I (không cảm thấy) đến cấp XII (thay đổi địa hình hoàn toàn) trên 3 nhóm công trình (Kiểu A: nhà đất/vật liệu địa phương; Kiểu B: nhà gạch, kết cấu khối lớn; Kiểu C: nhà bê tông cốt thép, khung gỗ tốt).
Bảng 1: Đặc trưng quy ước của dao động nền đất theo cấp độ mạnh (thang MSK-64)
| Cấp độ mạnh MSK-64 |
Gia tốc nền $\alpha$ ($\text{cm/s}^2$) (chu kỳ $0{,}1 \div 0{,}5\text{ s}$) |
Vận tốc nền $v$ ($\text{cm/s}$) (chu kỳ $0{,}5 \div 2{,}0\text{ s}$) |
Biên độ khối lượng con lắc $x_0$ ($\text{mm}$) ($T_0=0{,}25\text{ s}$, cản $0{,}5$) |
| Cấp V |
$12 \div 25$ |
$1{,}0 \div 2{,}0$ |
$0{,}5 \div 1{,}0$ |
| Cấp VI |
$25 \div 50$ |
$2{,}1 \div 4{,}0$ |
$1{,}1 \div 2{,}0$ |
| Cấp VII |
$50 \div 100$ |
$4{,}1 \div 8{,}0$ |
$2{,}1 \div 4{,}0$ |
| Cấp VIII |
$100 \div 200$ |
$8{,}1 \div 16{,}0$ |
$4{,}1 \div 8{,}0$ |
| Cấp IX |
$200 \div 400$ |
$16{,}1 \div 32{,}0$ |
$8{,}1 \div 16{,}0$ |
| Cấp X |
$400 \div 800$ |
$32{,}1 \div 64{,}0$ |
$16{,}1 \div 32{,}0$ |
3. Đặc trưng dao động nền và các công thức gia tốc cực đại:
Gia tốc cực đại nền đất ($a_{\max}$) là thông số trực tiếp trong tính toán kháng chấn. Do thiếu dữ liệu ghi chấn trực tiếp, các công thức thực nghiệm thống kê được áp dụng rộng rãi:
- Công thức Donovan (1973): $a_{\max} = \frac{1080 \cdot e^{0{,}5 M}}{(R + 25)^{1{,}32}}\text{ }(\text{cm/s}^2)$
- Công thức Cornell: $a_{\max} = \frac{1080 \cdot e^{0{,}8 M}}{(R + 25)^{1{,}32}}\text{ }(\text{cm/s}^2)$
- Công thức Esteva (1974): $a_{\max} = \frac{1230 \cdot e^{0{,}8 M}}{(R + 25)^2}\text{ }(\text{cm/s}^2)$
Tác giả nhận xét rằng ở các mức cấp chấn động $I = 5 \div 7$, công thức của Cornell và Donovan cho kết quả phù hợp với điều kiện nền đất tại Việt Nam. Quá trình dao động nền của động đất mạnh thường kéo dài $10 \div 40\text{ s}$, gồm 3 pha: pha đầu (biên độ nhỏ, tần số cao), pha chính (biên độ lớn nhất, chu kỳ $0{,}1 \div 0{,}5\text{ s}$) và pha cuối (biên độ giảm dần).
4. Phân vùng động đất tại Việt Nam:
Số liệu điều tra cho thấy đến năm 1995, Việt Nam ghi nhận hơn 700 trận động đất có $M_S < 7{,}0$. Các trận động đất mạnh tập trung dọc theo các đới đứt gãy kiến tạo sâu hoạt động:
- Động đất Điện Biên (tháng 01/1935) trên đứt gãy Sông Mã ($M_S = 6{,}8$, cấp $I_0 = 8 \div 9$).
- Động đất Tuần Giáo (ngày 24/6/1983) có $M_S = 6{,}7$, độ sâu chấn tiêu $h = 23 \div 25\text{ km}$, gây chấn động cấp $I_0 = 8 \div 9$ trên vùng tâm chấn rộng $2500\text{ km}^2$, chấn động cấp 7 lan rộng $13000\text{ km}^2$, nứt đất kéo dài hàng chục km.
- Động đất Lục Yên (tháng 11/1954) trên đứt gãy Sông Chảy.
- Động đất Bắc Giang (tháng 12/1961) trên đứt gãy Đông Triều.
Giới hạn trên của Magnitude cực đại dự báo trên các đứt gãy được xác định theo công thức của Nguyễn Đình Xuyên:
$$M_{\max} = 2\lg L(\text{km}) + 1{,}48$$
(trong đó $L$ là chiều dài đoạn đứt gãy nguyên vẹn).
Bảng 2: Tần suất động đất trong các vùng lãnh thổ ở Việt Nam (số lần xuất hiện trong một năm)
| Vùng lãnh thổ |
$M_S \ge 4{,}5$ ($I_0 \ge 6$) |
$M_S \ge 5{,}0$ ($I_0 \ge 7$) |
$M_S \ge 5{,}5$ ($I_0 \ge 7 \div 8$) |
$M_S \ge 6{,}0$ ($I_0 \ge 8$) |
$M_S \ge 6{,}5$ ($I_0 \ge 8 \div 9$) |
| Miền Bắc Việt Nam |
$0{,}6017$ |
$0{,}2254$ |
$0{,}0873$ |
$0{,}0338$ |
$0{,}0136$ |
| Miền Nam Việt Nam (địa khối Indosini) |
$0{,}1585$ |
$0{,}0525$ |
$0{,}0174$ |
- |
- |
| Vùng Đông Bắc Việt Nam (nền hoạt động Hoa Nam) |
$0{,}0803$ |
$0{,}0283$ |
$0{,}0103$ |
$0{,}0037$ |
- |
| Vùng Tây Bắc Việt Nam (miền uốn nếp Bắc Việt Nam) |
$0{,}5029$ |
$0{,}1988$ |
$0{,}0734$ |
$0{,}0300$ |
$0{,}0111$ |
Khu vực Tây Bắc (dọc đới Lai Châu - Điện Biên, Sông Mã, Sơn La) có tần suất chấn động cao nhất, với tần suất chấn động cấp 8 là $B_{I \ge 8} < 0{,}005$ và cấp 7 là $B_{I \ge 7} < 0{,}01$.
Chương 2: Một số lý thuyết kháng chấn trong thiết kế công trình
Chương 2 trình bày các nền tảng toán học và động lực học kết cấu phục vụ xác định nội lực và chuyển vị của công trình khi chịu động đất.
1. Tiến trình phát triển các phương pháp tính:
- Phương pháp tĩnh (Ohmori và Sano): Xuất phát sau các trận động đất Nobi (1891) và San Francisco (1906). Lực động đất được quy về tải trọng ngang tác dụng lên bộ phận thứ $k$:
$$F_k = K_s Q_k$$
(với $Q_k$ là trọng lượng bộ phận, $K_s$ là hệ số động đất thực nghiệm). Hạn chế cốt lõi của phương pháp tĩnh là coi công trình như vật rắn tuyệt đối, bỏ qua biến dạng đàn hồi và không phản ánh đúng trạng thái làm việc động lực học của kết cấu.
- Phương pháp động lực học: Hình thành từ những năm 1940, mô hình hóa công trình thành hệ đàn hồi có độ cứng $K$, hệ số cản $C$, chịu kích động từ hàm gia tốc chuyển động nền $\ddot{x}_0(t)$.
2. Phổ ứng xử của hệ một bậc tự do (SDOF):
Mô hình hóa kết cấu cầu (như một trụ cầu đỡ kết cấu nhịp) thành thanh đàn hồi một đầu ngàm, đỉnh mang khối lượng tập trung $m$. Phương trình vi phân dao động dưới kích động gia tốc nền $\ddot{x}_0(t)$:
$$m\ddot{x}(t) + c\dot{x}(t) + kx(t) = -m\ddot{x}_0(t)$$
Đặt $\omega^2 = k/m$ và $2\beta = c/m$, phương trình chuẩn hóa thành:
$$\ddot{x}(t) + 2\beta \dot{x}(t) + \omega^2 x(t) = -\ddot{x}_0(t)$$
Tỉ số cản tới hạn $\nu = c/c_0 = \beta/\omega$. Với kết cấu công trình thông thường, $\nu < 0{,}2$. Nghiệm chuyển vị tương đối bình ổn được xác định thông qua tích phân Duhamel:
$$x(t) = -\frac{1}{\omega} \int_0^t \ddot{x}_0(\tau) e^{-\nu\omega(t-\tau)} \sin\omega(t-\tau) d\tau$$
Từ đó định nghĩa các phổ ứng xử cực đại:
- Phổ chuyển vị tương đối: $S_d = |x(t)|_{\max}$
- Phổ vận tốc tương đối: $S_v = |\dot{x}(t)|_{\max}$
- Phổ gia tốc tuyệt đối: $S_a = |\ddot{x}0(t) + \ddot{x}(t)|{\max}$
Giữa các á phổ tồn tại mối quan hệ giải tích:
$$S_d = \frac{S_v}{\omega} = \frac{S_a}{\omega^2}$$
Lực động đất cực đại tác dụng lên hệ SDOF được tính trực tiếp từ phổ gia tốc:
$$F_{\max} = m S_a = m \omega S_v = m S_v \frac{2\pi}{T}$$
Đường bao các phổ ứng xử qua nhiều trận động đất trong điều kiện địa chất tương tự tạo thành phổ chuẩn (phổ thiết kế), tiêu biểu là phổ trận El Centro (1940). Phổ vận tốc $S_v$ biến đổi mạnh khi chu kỳ dao động riêng $T < 0{,}8\text{ s}$ và gần như không đổi khi $0{,}8\text{ s} \le T \le 2{,}8\text{ s}$.
3. Phổ ứng xử của hệ nhiều bậc tự do (MDOF):
Mô hình hóa kết cấu có $n$ bậc tự do chịu kích động nền được biểu diễn bằng hệ phương trình vi phân ma trận:
$$[M]{\ddot{x}} + [C]{\dot{x}} + [K]{x} = -[M]{1}\ddot{x}_0(t)$$
Phương trình tần số dao động riêng nhận được từ điều kiện triệt tiêu định thức:
$$|[K] - \omega^2 [M]| = 0$$
Giải phương trình thu được $n$ tần số dao động riêng $\omega_1 < \omega_2 < \dots < \omega_n$ và $n$ chu kỳ dao động riêng tương ứng $T_i = 2\pi/\omega_i$. Nhờ tính chất trực giao của các véc-tơ dạng dao động riêng ${\phi}_i$, hệ phương trình ma trận được tách thành $n$ phương trình vi phân độc lập theo từng dạng dao động chính thứ $i$.
Chuyển vị tương đối và lực động đất cực đại tại mức khối lượng $m_k$ trong dạng dao động thứ $i$:
$$x_{ki} = \eta_{ki} S_{di}$$
$$F_{ki} = \eta_{ki} m_k S_{ai}$$
Trong đó $\eta_{ki}$ là hệ số hình dáng:
$$\eta_{ki} = \phi_{ki} \frac{\sum_{j=1}^n m_j \phi_{ji}}{\sum_{j=1}^n m_j \phi_{ji}^2}$$
Tổng lực cắt ở đáy công trình trong dạng dao động thứ $i$:
$$F_i = S_{ai} \sum_{k=1}^n \eta_{ki} m_k$$
Hệ số phân phối lực dọc theo chiều cao:
$$P_{ki} = \frac{F_{ki}}{F_i} = \frac{Q_k \phi_{ki}}{\sum_{j=1}^n Q_j \phi_{ji}}$$
Do đó: $F_{ki} = P_{ki} F_i$ (với $Q_k = m_k g$).
4. Chuyển đổi hệ MDOF về các hệ SDOF tương đương:
Dựa trên nguyên lý tương đương năng lượng, hệ $n$ bậc tự do được quy đổi thành $n$ hệ một bậc tự do độc lập. Khối lượng tương đương của hệ một bậc tự do cho mode thứ $i$ xác định theo công thức:
$$m_e(i) = \frac{\left(\sum_{k=1}^n m_k \phi_{ki}\right)^2}{\sum_{k=1}^n m_k \phi_{ki}^2} = m \cdot \varepsilon(i)$$
Trong đó $\varepsilon(i)$ là hệ số tương đương khối lượng:
$$\varepsilon(i) = \frac{\left(\sum_{k=1}^n m_k \phi_{ki}\right)^2}{m \sum_{k=1}^n m_k \phi_{ki}^2}$$
Lực cắt đáy cực đại được tính đơn giản qua mô hình SDOF tương đương:
$$F_i = m_e(i) S_{ai} = \varepsilon(i) \frac{Q}{g} S_{ai}$$
5. Quy tắc tổ hợp phản ứng cực đại toàn phần:
Để xác định lực động đất tổng hợp từ các dạng dao động riêng, các tác giả điểm lại công thức của H. Merchant:
$$F_{\max} \approx \sqrt{\sum |F_i|^2 + \left|\sum F_i\right|^2}$$
và công thức của R. Clough:
$$F_{\max} \approx F_1 + \gamma |F_2|$$
(với $\gamma = 0$ khi chu kỳ dao động cơ bản $T_1 < 0{,}8\text{ s}$ và $\gamma = 0{,}25$ khi $T_1 > 0{,}8\text{ s}$). Đối với công trình thông thường, ứng xử chịu ảnh hưởng chủ yếu của 3 dạng dao động đầu tiên; riêng công trình có kết cấu mảnh (tháp cầu, trụ cầu rất cao), cần xét đến 5 dạng dao động đầu.
6. Phương pháp phân tích lịch sử thời gian (Time-History):
Áp dụng cho các kết cấu có khối lượng tập trung ở đỉnh (trụ cao, tháp cầu). Phương trình vi phân dao động được giải bằng thuật toán sai phân hữu hạn từng bước thời gian: chuyển vị tương đối $x^{(t+\Delta t)}$ tại thời điểm tiếp theo được tính toán truy hoàn từ các giá trị chuyển vị tại thời điểm $t$ và $t-\Delta t$ cùng giá trị gia tốc nền $\ddot{x}_0(t)$.
Kết quả và đóng góp
Kết quả chính kèm số liệu thực chứng
- Tổng hợp số liệu địa chấn phục vụ ngành cầu: Hệ thống hóa chuỗi số liệu hơn 700 trận động đất tại Việt Nam (giai đoạn 1900–1995), xác lập các thông số chấn động cực đại cho các trận động đất điển hình (Điện Biên 1935 đạt $M_S = 6{,}8$, Tuần Giáo 1983 đạt $M_S = 6{,}7$ với độ sâu chấn tiêu $23 \div 25\text{ km}$, vùng phá hoại cấp 8–9 rộng $2500\text{ km}^2$).
- Xác lập bảng tần suất động đất theo vùng: Cung cấp bảng dữ liệu định lượng về tần suất xuất hiện động đất phân theo các cấp Magnitude ($M_S \ge 4{,}5$ đến $M_S \ge 6{,}5$) cho 4 vùng kiến tạo: Miền Bắc, Miền Nam, Đông Bắc và Tây Bắc.
- Đánh giá tính thích ứng của công thức suy giảm gia tốc: Xác định công thức thực nghiệm của Donovan (1973) và Cornell phù hợp để tính toán gia tốc nền dao động ($a_{\max}$) ở các vùng động đất cấp 5 đến cấp 7 tại Việt Nam.
- Khép kín hệ thống công thức giải tích kháng chấn cho kết cấu cầu: Cung cấp hệ thống giải tích hoàn chỉnh từ phương trình vi phân chuyển động, hệ số hình dáng $\eta_{ki}$, hệ số phân phối tải trọng $P_{ki}$ đến hệ số tương đương khối lượng $\varepsilon(i)$ phục vụ đưa bài toán dao động không gian phức tạp của cầu về hệ một bậc tự do tương đương.
Giải pháp và kiến nghị của tác giả
- Kiến nghị Bộ Giao thông Vận tải và các cơ quan chuyên môn sớm nghiên cứu, xây dựng và ban hành tiêu chuẩn thống nhất về thiết kế kháng chấn cho công trình cầu tại Việt Nam, khắc phục tình trạng phụ thuộc tản mát vào nhiều quy trình nước ngoài trái ngược nhau.
- Đề xuất khai thác các chương trình máy tính sẵn có để phát triển các phần mềm phụ trợ chuyên dụng bằng tiếng Việt, chạy trên các hệ điều hành thông dụng, giúp kỹ sư trong nước tự chủ tính toán chuyển vị và nội lực động đất mà không phải phụ thuộc vào việc thuê chuyên gia nước ngoài hay mua phần mềm ngoại đắt đỏ.
- Khuyến nghị trong tính toán thiết kế cầu thông thường cần xét tối thiểu 3 dạng dao động đầu tiên; đối với các kết cấu thanh mảnh (trụ cầu quá cao, tháp cầu dây văng/treo), bắt buộc phải tính toán đến 5 dạng dao động đầu tiên.
Đóng góp kỹ thuật
Tài liệu làm rõ bản chất động lực học công trình và cơ chế làm việc của kết cấu cầu chịu tải trọng động đất. Bằng việc phân tích so sánh đối chiếu giữa Tiêu chuẩn ngành 22 TCN 221-95 với Tiêu chuẩn AASHTO-1996 và 22 TCN 272-01, nghiên cứu cung cấp cơ sở đánh giá độ an toàn và độ lệch nội lực kháng chấn, tạo tiền đề khoa học cho quá trình hội nhập và thống nhất quy phạm thiết kế cầu đường tại Việt Nam.
Hạn chế và hướng nghiên cứu tiếp
Hạn chế được ghi nhận trong văn bản
- Nghiên cứu giới hạn tính chất vật liệu trong phạm vi đàn hồi tuyến tính cổ điển, chưa khảo sát ứng xử phi tuyến hình học và phi đàn hồi (biến dạng dẻo) của vật liệu khi xảy ra chấn động vượt ngưỡng.
- Giả thiết chuyển vị của công trình là chuyển vị nhỏ quanh vị trí cân bằng bền.
- Mô hình tính toán chưa đi sâu vào hệ kết cấu vô hạn bậc tự do (kết cấu bản hoặc tường chịu lực) do đặc thù kết cấu cầu chủ yếu làm việc dạng thanh và khung chịu lực tập trung.
- Mối tương quan gia tốc cực đại $a_{\max}$ theo khoảng cách chấn tiêu và Magnitude của các công thức thực nghiệm (Donovan, Cornell) chưa phản ánh rõ quy luật ở vùng quá gần tâm chấn ($\Delta < 15 \div 20\text{ km}$).
- Do những hạn chế nội tại của tiêu chuẩn 22 TCN 221-95, việc so sánh nội lực kháng chấn giữa các tiêu chuẩn mới chỉ dừng lại ở mức độ sơ bộ trên một số dạng kết cấu cơ bản.
Hướng nghiên cứu tiếp
- Mở rộng phân tích ứng xử phi đàn hồi của mố, trụ và gối cầu dưới tác dụng của tải trọng động đất mạnh.
- Hoàn thiện thuật toán và giao diện cho các phần mềm phụ trợ tính toán động đất bằng ngôn ngữ tiếng Việt để chuyển giao rộng rãi cho các đơn vị tư vấn thiết kế cầu đường trong nước.
- Tiếp tục nghiên cứu thực nghiệm dao động nền cục bộ tại các địa bàn xây dựng cầu trọng điểm nhằm hiệu chỉnh đường cong phổ phản ứng thiết kế phù hợp điều kiện địa chất cụ thể ở Việt Nam.
Giá trị tham khảo
Đối tượng thụ hưởng
- Kỹ sư thiết kế, chuyên gia thẩm tra kết cấu công trình cầu đường tại các công ty tư vấn giao thông.
- Cán bộ quản lý kỹ thuật, quản lý dự án xây dựng cơ bản ngành giao thông vận tải.
- Giảng viên, học viên cao học và sinh viên chuyên ngành Xây dựng Cầu đường, Công trình Cầu - Hầm, Xây dựng Công trình Thủy và Công trình ngầm.
Các nội dung có giá trị tham khảo cao
- Phương pháp luận phân tích phổ ứng xử và phương pháp tương đương năng lượng chuyển đổi hệ MDOF về SDOF trong Chương 2: Cung cấp đầy đủ các biểu thức toán học tường minh, thuận tiện cho việc lập bảng tính Excel hoặc viết code giải thuật tự động.
- Dữ liệu địa chấn học và bản đồ phân vùng tần suất động đất tại Chương 1: Cung cấp số liệu đầu vào về chấn cấp cực đại, tần suất lặp lại chấn động theo chu kỳ 20, 50, 100 năm tại các vùng lãnh thổ Việt Nam để phục vụ lập báo cáo đầu tư và thiết kế cơ sở.
- Các công thức thực nghiệm xác định gia tốc cực đại nền ($a_{\max}$) và mối tương quan giữa thang MSK-64 với các thông số dao động vật lý ($a, v, x_0$).
Câu hỏi thường gặp
1. Đại lượng Magnitude và cấp động đất (Intensity) khác nhau như thế nào?
Magnitude là đại lượng vật lý đặc trưng cho toàn bộ mức năng lượng được giải phóng từ chấn tiêu dưới dạng sóng đàn hồi (đo bằng độ Richter qua các thang $M_L, M_S, M_B$), mỗi trận động đất chỉ có duy nhất một giá trị Magnitude. Ngược lại, cấp động đất (như thang MSK-64 gồm 12 cấp) biểu thị cường độ chấn động và mức độ phá hoại thực tế trên mặt đất đối với nhà cửa, công trình và địa hình tại một vị trí cụ thể; cường độ này thay đổi phụ thuộc vào khoảng cách từ điểm quan sát đến tâm chấn và điều kiện địa chất nền đất tại chỗ.
2. Dựa trên cơ sở cơ học nào để chuyển hệ kết cấu cầu nhiều bậc tự do (MDOF) về hệ một bậc tự do (SDOF) tương đương?
Việc chuyển đổi dựa trên nguyên lý tương đương về năng lượng giữa hệ thực và hệ tương đương trong từng dạng dao động chính thứ $i$. Bốn điều kiện bảo toàn bắt buộc gồm:
- Chu kỳ dao động riêng bằng nhau ($T_i = T_{(i)}$);
- Lực cắt ở chân công trình bằng nhau ($F_i = F_{(i)}$);
- Mômen uốn ở chân công trình bằng nhau ($M_i = M_{(i)}$);
- Động năng của hai hệ trong dạng dao động thứ $i$ bằng nhau ($E_i = E_{(i)}$).
Thông qua các điều kiện này, khối lượng của hệ được quy đổi về khối lượng tương đương $m_e(i) = m \cdot \varepsilon(i)$.
3. Quy tắc tổ hợp dao động theo R. Clough xác định lực động đất cực đại toàn phần như thế nào?
Theo nghiên cứu của R. Clough trên các công trình chịu động đất ở vùng California, ứng xử cực đại toàn phần được xác định bằng cách lấy ứng xử cực đại ở dạng dao động cơ bản cộng với một phần ứng xử cực đại ở dạng dao động thứ hai:
$$F_{\max} \approx F_1 + \gamma |F_2|$$
Hệ số thực nghiệm $\gamma$ phụ thuộc vào độ dài chu kỳ dao động cơ bản $T_1$:
- Với công trình có chu kỳ dao động $T_1 < 0{,}8\text{ s}$: lấy $\gamma = 0$ (chỉ xét dạng dao động thứ nhất);
- Với công trình có chu kỳ dao động $T_1 \ge 0{,}8\text{ s}$: lấy $\gamma = 0{,}25$.
4. Tại sao các kỹ sư Việt Nam gặp khó khăn khi ứng dụng các phần mềm tính toán kháng chấn nhập khẩu từ nước ngoài?
Nguyên nhân chủ yếu do các phần mềm chuyên dụng phục vụ các quy trình như AASHTO hay tiêu chuẩn Nhật Bản có giá thành bản quyền rất cao nên hầu như không được chuyển giao sau dự án. Ngoài ra, một số phần mềm của Nhật Bản sử dụng ngôn ngữ tiếng Nhật và vận hành trên các hệ điều hành chuyên biệt khác với các hệ điều hành thông dụng tại Việt Nam, đòi hỏi chuyên gia nước ngoài trực tiếp vận hành.
5. Những trận động đất mạnh có Magnitude lớn nhất từng ghi nhận trên lãnh thổ Việt Nam xảy ra ở khu vực nào?
Văn bản ghi nhận hai trận động đất mạnh nhất gồm: trận động đất Điện Biên (tháng 01/1935) xảy ra trên đới đứt gãy Sông Mã với $M_S = 6{,}8$ (chấn tâm cấp 8–9 theo thang MSK-64) và trận động đất Tuần Giáo (ngày 24/6/1983) có $M_S = 6{,}7$, độ sâu chấn tiêu $23 \div 25\text{ km}$, gây chấn động cấp 8–9 trên diện tích $2500\text{ km}^2$ và cấp 7 trên vùng rộng $13000\text{ km}^2$.
Kết luận
Tài liệu cung cấp hệ thống cơ sở lý thuyết động lực học công trình và các phương pháp tính toán kết cấu cầu chịu tải trọng động đất một cách mạch lạc, chặt chẽ. Nội dung chuyên khảo đã liên kết thành công các dữ liệu thực chứng địa chấn tại Việt Nam với các phương pháp giải tích hiện đại như phổ ứng xử và tích phân lịch sử thời gian. Công trình đặt nền móng kỹ thuật quan trọng cho việc so sánh, đánh giá sự tương quan giữa tiêu chuẩn ngành 22 TCN 221-95 với các tiêu chuẩn quốc tế (AASHTO-1996, 22 TCN 272-01), đồng thời định hướng việc tự chủ công cụ tính toán kháng chấn cho đội ngũ kỹ sư cầu đường trong nước.
