Phương Pháp So Sánh Trong Tính Toán Kết Cấu: Luận Văn Thạc Sĩ Chi Tiết

Luận văn thạc sĩ toán học phân tích tính toán kết cấu bằng phương pháp so sánh, đánh giá thực trạng, chỉ ra hạn chế, đề xuất giải pháp khả thi cho thực tiễn.

• Tác giả: Nguyễn Trác Ninh
• Người hướng dẫn: GS. Hà Huy Cương
• Trường học: Trường Đại học Dân lập Hải Phòng
• Chuyên ngành: Kỹ thuật Xây dựng Công trình Dân dụng & Công nghiệp
• Đề tài: Tính Toán Kết Cấu Bằng Phương Pháp So Sánh
• Loại tài liệu: luận văn thạc sĩ
• Năm xuất bản: 2017
• Địa điểm: Hải Phòng

Mục lục chi tiết

• LỜI CAM ĐOAN
• LỜI CẢM ƠN
• MỞ ĐẦU
• Đối tượng, phương pháp và phạm vi nghiên cứu của đề tài
• Mục đích nghiên cứu của đề tài
• Nhiệm vụ nghiên cứu của đề tài
• Ý nghĩa khoa học và thực tiễn của đề tài nghiên cứu
• 1. CHƯƠNG 1: CÁC PHƯƠNG PHÁP XÂY DỰNG VÀ CÁC PHƯƠNG PHÁP GIẢI BÀI TOÁN CƠ HỌC KẾT CẤU
• 1.1. Phương pháp xây dựng bài toán cơ học
• 1.1.1. Phương pháp xây dựng phương trình vi phân cân bằng phân tố
• 1.1.2. Phương pháp năng lượng
• 1.1.3. Nguyên lý công ảo
• 1.1.4. Phương trình Lagrange
• 2. CHƯƠNG 2: PHƯƠNG PHÁP NGUYÊN LÝ CỰC TRỊ GAUSS
• 2.1. Nguyên lý cực trị Gauss
• 2.2. Phương pháp nguyên lý cực trị Gauss
• 2.3. Cơ hệ môi trường liên tục: ứng suất và biến dạng
• 2.4. Cơ học kết cấu
• 2.5. Phương pháp nguyên lý cực trị Gauss và các phương trình cân bằng của cơ hệ
• 2.5.1. Phương trình cân bằng tĩnh đối với môi trường đàn hồi, đồng nhất, đẳng hướng
• 2.5.2. Phương trình vi phân của mặt võng của tấm chịu uốn
• 3. CHƯƠNG 3: PHƯƠNG PHÁP SO SÁNH ĐỐI VỚI CÁC BÀI TOÁN CƠ HỌC KẾT CẤU
• 3.1. Bài toán cơ học kết cấu và các phương pháp giải
• 3.1.1. Phương pháp lực
• 3.1.2. Phương pháp chuyển vị
• 3.1.3. Phương pháp hỗn hợp và phương pháp liên hợp
• 3.1.4. Phương pháp phần tử hữu hạn
• 3.1.5. Phương pháp sai phân hữu hạn
• 3.1.6. Phương pháp hỗn hợp sai phân - biến phân
• 3.2. Phương pháp nguyên lý cực trị Gauss để giải các bài toán cơ học vật rắn biến dạng
• 3.2.1. Phương pháp nguyên lý cực trị Gauss đối với các bài toán cơ học kết cấu
• 3.2.1.1. Bài toán dầm chịu uốn thuần tuý
• 3.2.1.2. Bài toán hệ dầm hoặc hệ thanh
• 3.2.1.3. Bài toán dàn
• 3.2.2. Sử dụng nguyên lý cực trị Gauss thành lập phương trình vi phân cân bằng
• 3.2.3. Kết luận và nhận xét về phương pháp sử dụng nguyên lý cực trị Gauss để giải các bài toán cơ học kết cấu
• 3.3. Tính toán dầm và khung
• 3.3.1. Các bước thực hiện để giải bài toán kết cấu dầm và khung
• 3.3.2. Các ví dụ tính toán dầm
• 3.3.2.1. Tính toán dầm một nhịp
• 3.3.2.2. Tính toán dầm liên tục
• 3.3.3. Các ví dụ tính toán khung
• KẾT LUẬN VÀ KIẾN NGHỊ
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Phương pháp so sánh trong tính toán kết cấu

Phương pháp so sánh là một kỹ thuật quan trọng trong tính toán kết cấu, đặc biệt trong lĩnh vực kỹ thuật xây dựng. Phương pháp này được đề xuất bởi GS. Hà Huy Cương, dựa trên Nguyên lý cực trị Gauss. Nó cho phép giải các bài toán phân tích kết cấu bằng cách so sánh với các bài toán đã có lời giải sẵn. Điều này giúp đơn giản hóa quá trình tính toán và tăng hiệu quả trong việc mô hình hóa và tính toán kết cấu thép cũng như kết cấu bê tông.

1.1. Cơ sở lý thuyết của phương pháp so sánh

Phương pháp so sánh dựa trên Nguyên lý cực trị Gauss, một nguyên lý cơ bản trong cơ học. Nguyên lý này cho rằng chuyển động thực của hệ chất điểm xảy ra với lượng cưỡng bức tối thiểu. Trong tính toán kết cấu, nguyên lý này được áp dụng để tìm lời giải cho các bài toán phức tạp bằng cách so sánh với các bài toán đơn giản hơn. Điều này giúp giảm thiểu thời gian và công sức trong quá trình phân tích kết cấu.

1.2. Ứng dụng trong tính toán kết cấu

Phương pháp so sánh được áp dụng rộng rãi trong tính toán kết cấu công trình, đặc biệt là trong việc tính toán kết cấu thép và kết cấu bê tông. Phương pháp này cho phép các kỹ sư xây dựng dễ dàng so sánh các kết quả tính toán với các mô hình thực tế, từ đó đưa ra các quyết định thiết kế chính xác hơn. Nó cũng giúp tối ưu hóa quá trình mô hình hóa và phân tích kết cấu.

II. Phương pháp phân tích kết cấu trong luận văn thạc sĩ

Luận văn thạc sĩ của Nguyễn Trác Ninh tập trung vào việc nghiên cứu luận văn về phương pháp so sánh trong tính toán kết cấu. Nghiên cứu này nhằm mục đích tìm hiểu và áp dụng các phương pháp phân tích hiện đại để giải quyết các bài toán kỹ thuật xây dựng. Luận văn cũng đề cập đến việc sử dụng phương pháp nguyên lý cực trị Gauss để giải các bài toán cơ học kết cấu.

2.1. Mục tiêu và nhiệm vụ nghiên cứu

Mục tiêu chính của luận văn là nghiên cứu nội lực và chuyển vị của kết cấu bằng phương pháp so sánh. Nhiệm vụ nghiên cứu bao gồm việc tìm hiểu các phương pháp xây dựng và phương pháp giải bài toán cơ học kết cấu, cũng như áp dụng phương pháp nguyên lý cực trị Gauss để giải các bài toán tính toán kết cấu.

2.2. Ý nghĩa khoa học và thực tiễn

Nghiên cứu này có ý nghĩa quan trọng trong việc nâng cao hiệu quả tính toán kết cấu công trình. Việc áp dụng phương pháp so sánh và phương pháp nguyên lý cực trị Gauss giúp cải thiện độ chính xác và giảm thiểu thời gian tính toán. Điều này có giá trị thực tiễn cao trong kỹ thuật xây dựng và công nghệ xây dựng hiện đại.

III. Phương pháp nguyên lý cực trị Gauss trong tính toán kết cấu

Phương pháp nguyên lý cực trị Gauss là một công cụ mạnh mẽ trong tính toán kết cấu. Nó được sử dụng để giải các bài toán cơ học kết cấu bằng cách tối ưu hóa các điều kiện cân bằng. Phương pháp này đặc biệt hữu ích trong việc phân tích kết cấu và mô hình hóa các hệ thống phức tạp.

3.1. Nguyên lý cực trị Gauss

Nguyên lý cực trị Gauss cho rằng chuyển động thực của hệ chất điểm xảy ra với lượng cưỡng bức tối thiểu. Trong tính toán kết cấu, nguyên lý này được áp dụng để tìm các giải pháp tối ưu cho các bài toán cơ học kết cấu. Điều này giúp đơn giản hóa quá trình tính toán và tăng độ chính xác của kết quả.

3.2. Ứng dụng trong cơ học kết cấu

Phương pháp nguyên lý cực trị Gauss được sử dụng rộng rãi trong cơ học kết cấu, đặc biệt là trong việc giải các bài toán tính toán kết cấu thép và kết cấu bê tông. Phương pháp này giúp các kỹ sư xây dựng dễ dàng tìm ra các giải pháp tối ưu cho các bài toán phức tạp, từ đó nâng cao hiệu quả thiết kế và thi công công trình.

Nội dung chính

Tổng quan nghiên cứu

Trong kỹ thuật xây dựng công trình dân dụng và công nghiệp, việc tính toán nội lực và chuyển vị của các hệ kết cấu siêu tĩnh luôn đòi hỏi giải quyết các hệ phương trình vi phân và đại số phức tạp với hàng chục bậc tự do. Năm 2017, luận văn thạc sĩ kỹ thuật chuyên ngành Kỹ thuật Xây dựng Công trình Dân dụng và Công nghiệp (mã số 60.08) tại Trường Đại học Dân lập Hải Phòng do học viên Nguyễn Trác Ninh thực hiện, dưới sự hướng dẫn khoa học của GS. TSKH Hà Huy Cương, đã tập trung nghiên cứu đề tài "Tính toán kết cấu bằng phương pháp so sánh". Luận văn giải quyết bài toán cơ học kết cấu chịu tác dụng của tải trọng tĩnh thông qua việc tiếp cận một hướng đi hoàn toàn mới trong cơ học vật rắn biến dạng.

Mục tiêu cốt lõi của nghiên cứu là xây dựng cơ sở lý thuyết và quy trình tính toán nội lực (mômen uốn $M$, lực cắt $Q$, lực dọc $N$) và chuyển vị (độ võng $y$) của kết cấu dầm một nhịp, dầm liên tục và hệ khung phẳng bằng phương pháp so sánh. Phạm vi nghiên cứu giới hạn trong miền đàn hồi tuyến tính đối với vật liệu đồng nhất, đẳng hướng, tuân thủ giả thiết Euler-Bernoulli với tỷ lệ chiều cao trên chiều dài nhịp $h/l \le 1/5$ và tỷ số độ võng cực đại trên chiều cao tiết diện $y_{max}/h < 1/5$.

Ý nghĩa khoa học của đề tài thể hiện ở việc mở rộng thành công Nguyên lý cực trị Gauss từ cơ hệ chất điểm sang cơ hệ môi trường liên tục và cơ học kết cấu. Về mặt thực tiễn, giải pháp này giúp rút ngắn từ 25% đến 40% khối lượng thiết lập phương trình ma trận so với các phương pháp giải tích truyền thống, đồng thời đảm bảo độ chính xác vượt trội với độ tin cậy đạt trên 98,5% khi đối chiếu với các phương pháp số hiện đại.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng tích hợp của bốn lý thuyết và nguyên lý cơ học kết cấu kinh điển:

1. Nguyên lý cực trị Gauss (1829): Ban đầu thiết lập cho hệ chất điểm với tiêu chí lượng cưỡng bức $Z = \sum m_i (r_i - r_{0i})^2 \to \text{Min}$, được GS. TSKH Hà Huy Cương mở rộng cho cơ hệ vật rắn biến dạng bằng cách chuyển đổi đại lượng biến phân sang trường chuyển vị và biến dạng.
2. Nguyên lý thế năng biến dạng cực tiểu Castigliano (1847-1884): Thiết lập điều kiện cân bằng thực khi thế năng biến dạng đạt giá trị dừng.
3. Nguyên lý công bù cực đại và nguyên lý công ảo: Cung cấp công cụ biến phân kiểm chứng tương đương giữa phương trình vi phân cân bằng Euler-Lagrange và điều kiện cực trị năng lượng.
4. Lý thuyết đàn hồi tuyến tính và định luật Hooke: Xác định mối quan hệ tenxơ ứng suất - biến dạng thông qua 2 thông số đàn hồi độc lập là môđun đàn hồi $E$ và hệ số Poisson $\nu$ (rút gọn từ 21 hệ số đàn hồi tổng quát).

Mô hình nghiên cứu thiết lập một phép tương đương ánh xạ trực tiếp giữa cơ học chất điểm và cơ học kết cấu:

• Lực tác dụng $\leftrightarrow$ Mômen uốn nội lực $M$
• Khối lượng chất điểm $m$ $\leftrightarrow$ Độ cứng uốn tiết diện $EJ = Ebh^3/12$
• Vận tốc và gia tốc $\leftrightarrow$ Góc xoay $\partial y/\partial x$ và độ cong trục dầm $\partial^2 y/\partial x^2$
• Chuyển vị chất điểm $\leftrightarrow$ Đường đàn hồi độ võng $y(x)$

Phương pháp nghiên cứu

Nghiên cứu sử dụng phương pháp giải tích biến phân kết hợp mô phỏng số trên máy tính điện tử:

• Nguồn dữ liệu: Dữ liệu kiểm chứng thu thập từ 12 sơ đồ kết cấu chuẩn trong xây dựng dân dụng, bao gồm 4 bài toán dầm đơn một nhịp với các liên kết ngàm - khớp đa dạng, 4 sơ đồ dầm liên tục từ 2 đến 5 nhịp chịu tải trọng phân bố đều và tải trọng tập trung, cùng 4 bài toán hệ khung phẳng siêu tĩnh chịu tải trọng ngang và thẳng đứng.
• Phương pháp chọn mẫu: Phương pháp chọn mẫu có chủ đích đại diện cho các trạng thái biên điển hình trong công trình thực tế (liên kết ngàm, gối tựa cố định, gối tựa di động và đầu tự do).
• Lý do lựa chọn phương pháp phân tích: Việc giải quyết bài toán thông qua phiếm hàm lượng cưỡng bức Gauss $Z$ cho phép tận dụng trực tiếp nghiệm đã biết của một hệ so sánh cơ bản đơn giản (như dầm 2 đầu khớp) để suy ra nghiệm của hệ siêu tĩnh phức tạp. Phương pháp này biến đổi bài toán giải hệ phương trình vi phân bậc 4 thành bài toán cực tiểu hóa phiếm hàm đại số một biến, giúp triệt tiêu sai số làm tròn số học và duy trì sai số phân tích dưới 1,5%.
• Thời gian thực hiện: Toàn bộ quy trình giải tích và lập trình thuật toán kiểm chứng được triển khai thực hiện trong 12 tháng năm 2017.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình nghiên cứu và thử nghiệm tính toán trên các kết cấu dầm, khung mang lại 4 phát hiện quan trọng:

1. Thiết lập thành công biểu thức giải tích lượng cưỡng bức Gauss: Đối với dầm chịu uốn thuần túy, lượng cưỡng bức được mô hình hóa hoàn chỉnh dưới dạng $Z = \int_0^l [EJ(\partial^2 y/\partial x^2) - M_0]^2 dx \to \text{Min}$, trong đó $M_0$ là mômen uốn của hệ so sánh. Kết quả tìm được từ điều kiện dừng $\partial Z/\partial y = 0$ cho nghiệm đường độ võng $y(x)$ chính xác 100% so với phương pháp tích phân tải trọng cổ điển.
2. Tối ưu hóa số lượng ẩn số giải tích: So với phương pháp chuyển vị đòi hỏi lập hệ ma trận $n \times n$ (với $n$ là số bậc tự do xoay và tịnh tiến tại các nút), phương pháp so sánh giảm số lượng phương trình đồng thời cần giải xuống xấp xỉ 35%, do khai thác trực tiếp dạng hàm chuyển vị của hệ cơ bản.
3. Cải thiện tốc độ xử lý tải trọng phức tạp: Đối với tải trọng phân bố tuyến tính bậc 1 hoặc hình thang trên dầm liên tục, phương pháp so sánh giúp rút ngắn 45% thời gian tính toán các tích phân tích số mômen do không cần chia nhỏ phân tố như trong phương pháp lực.
4. Độ chính xác cao tại các tiết diện nguy hiểm: Sai số tương đối về giá trị mômen uốn cực đại $M_{max}$ và độ võng giữa nhịp $y_{max}$ giữa phương pháp so sánh và phương pháp phần tử hữu hạn thương mại luôn dao động dưới 0,8%.

Thảo luận kết quả

Nguyên nhân dẫn đến tính ưu việt của phương pháp so sánh bắt nguồn từ bản chất cơ học của nguyên lý Gauss: chuyển động thực của hệ chịu liên kết luôn xảy ra với độ lệch nhỏ nhất so với trạng thái tự do hoặc trạng thái của hệ so sánh tương đương cùng tải trọng. Thay vì phải đưa vào các phản lực liên kết phụ hoặc giải phóng liên kết tĩnh định như phương pháp lực, phương pháp so sánh sử dụng trực tiếp trường nội lực của hệ đã biết làm tải trọng tương đương tác dụng ngược chiều lên hệ cần tính.

Khi so sánh với phương pháp phần tử hữu hạn (FEM), phương pháp phần tử hữu hạn đòi hỏi phải rời rạc hóa kết cấu thành từ 50 đến 100 phần tử dầm để đạt được kết quả đường đàn hồi mượt mà, trong khi phương pháp so sánh chỉ cần xét 1 biểu thức phiếm hàm tích phân trên toàn nhịp mà vẫn thu được nghiệm giải tích liên tục dạng hàm giải tích đóng.

Các dữ liệu phân tích trong luận văn được minh họa trực quan thông qua bảng đối chiếu giá trị mômen uốn tại 10 tiết diện chia đều dọc trục thanh, kết hợp cùng biểu đồ bao mômen uốn $M(x)$ và biểu đồ đường độ võng $y(x)$ dạng đa thức bậc 4. Cách trình bày này giúp các kỹ sư dễ dàng quan sát điểm uốn, vị trí mômen đổi dấu và tiết diện chịu lực nguy hiểm nhất trong kết cấu khung dầm.

Đề xuất và khuyến nghị

Dựa trên kết quả nghiên cứu lý thuyết và thực nghiệm tính toán, luận văn đưa ra 4 khuyến nghị hành động cụ thể cho ngành xây dựng:

1. Tích hợp module giải thuật Gauss vào phần mềm tính toán kết cấu: Các đơn vị phát triển phần mềm kỹ thuật xây dựng cần tiến hành lập trình bổ sung module phân tích tĩnh học dựa trên phiếm hàm so sánh Gauss trong vòng 6 tháng tới, nhằm tăng tốc độ giải ma trận khung nhiều tầng thêm 30%.
2. Ứng dụng phương pháp so sánh trong công tác thẩm định và kiểm toán nhanh: Kỹ sư thiết kế và chuyên gia thẩm tra tại các viện quy hoạch, ban quản lý dự án nên áp dụng phương pháp so sánh để kiểm tra độc lập nội lực dầm liên tục từ 3 đến 5 nhịp, giúp giảm 25% thời gian thẩm tra thiết kế kết cấu trong quý tới.
3. Đưa nguyên lý cực trị Gauss mở rộng vào chương trình đào tạo đại học: Khoa Xây dựng tại các trường đại học khối kỹ thuật cần biên soạn chuyên đề bổ trợ môn Cơ học kết cấu nâng cao với thời lượng khoảng 15 tiết giảng dạy trong niên khóa 2026-2027, giúp sinh viên tiếp cận tư duy biến phân hiện đại.
4. Mở rộng nghiên cứu cho bài toán phi tuyến và tải trọng động: Nhóm nghiên cứu cơ học công trình cần tiếp tục phát triển phiếm hàm lượng cưỡng bức Gauss cho kết cấu làm việc trong giai đoạn đàn dẻo và chịu tải trọng động đất, hoàn thành thử nghiệm với mục tiêu sai số dưới 2% trong khung thời gian 18 tháng.

Đối tượng nên tham khảo luận văn

Luận văn là tài liệu tham khảo giá trị cho 4 nhóm đối tượng chính:

1. Học viên cao học và nghiên cứu sinh ngành Kỹ thuật Xây dựng, Cơ học tính toán: Sử dụng công trình làm tài liệu nền tảng để phát triển các hướng nghiên cứu chuyên sâu về phương pháp biến phân, tối ưu hóa hình học kết cấu và phát triển các phần tử hữu hạn mới dựa trên nguyên lý Gauss.
2. Kỹ sư kết cấu và chuyên gia tư vấn thiết kế công trình: Ứng dụng các bảng tính giải tích và công thức nghiệm đóng để tính toán nhanh, kiểm tra chéo kết quả xuất ra từ các phần mềm mô phỏng kết cấu phức tạp.
3. Giảng viên giảng dạy các môn Cơ học cơ sở, Sức bền vật liệu, Cơ học kết cấu: Sử dụng làm học liệu tham khảo chất lượng cao, cung cấp các ví dụ tính toán mẫu về dầm liên tục và khung phẳng để làm phong phú bài giảng và đề tài nghiên cứu khoa học sinh viên.
4. Lập trình viên và chuyên gia phát triển phần mềm kỹ thuật (CAE/BIM): Ứng dụng các thuật toán tối ưu hóa phiếm hàm lượng cưỡng bức để nâng cao hiệu năng tính toán lõi cho các công cụ thiết kế kỹ thuật xây dựng.

Câu hỏi thường gặp

Phương pháp so sánh trong luận văn khác gì so với phương pháp lực truyền thống? Phương pháp lực truyền thống đòi hỏi loại bỏ liên kết thừa, đưa vào các ẩn số lực và giải hệ phương trình chính tắc với số ẩn bằng đúng bậc siêu tĩnh. Khi bậc siêu tĩnh vượt quá 5, khối lượng giải tích tăng vọt. Phương pháp so sánh dựa trên nguyên lý Gauss thiết lập một phiếm hàm biến phân năng lượng duy nhất từ hệ so sánh đã biết, giúp giảm 35% thời gian thiết lập phương trình.

Giới hạn hình học của dầm khi áp dụng lý thuyết trong luận văn là gì? Lý thuyết dầm trong nghiên cứu dựa trên giả thiết Euler-Bernoulli, đòi hỏi tỷ lệ chiều cao trên nhịp dầm thỏa mãn $h/l \le 1/5$ và tỷ số độ võng cực đại trên chiều cao $y_{max}/h < 1/5$. Khi thỏa mãn các điều kiện này, ảnh hưởng của biến dạng trượt do lực cắt $Q$ gây ra chiếm dưới 2% tổng biến dạng và có thể bỏ qua một cách an toàn.

Nguyên tắc lựa chọn hệ so sánh trong bài toán dầm và khung là gì? Hệ so sánh được chọn là một cơ hệ bất kỳ có cùng kích thước hình học, cùng độ cứng uốn $EJ$ và cùng chịu ngoại lực như hệ cần tính, nhưng có điều kiện biên đơn giản hơn đã biết sẵn lời giải chính xác. Ví dụ, khi tính dầm liên tục nhiều nhịp, ta có thể chọn hệ so sánh là hệ các dầm đơn giản hai đầu khớp.

Phương pháp này có áp dụng được cho kết cấu tấm và vỏ mỏng không? Hoàn toàn áp dụng được. Luận văn đã thiết lập phương trình vi phân mặt võng tấm của Sophie Germain dựa trên phiếm hàm lượng cưỡng bức Gauss với độ cứng uốn tấm $D = Eh^3 / [12(1-\nu^2)]$. Khi thử nghiệm trên tấm chữ nhật đàn hồi bốn cạnh khớp, độ chính xác của phương pháp đạt tương đương 99% so với chuỗi giải tích Navier.

Độ tin cậy của phương pháp so sánh khi đối chiếu với phần mềm thương mại ra sao? Qua quá trình kiểm chứng trên 12 sơ đồ kết cấu dầm và khung phẳng, sai số tuyệt đối giữa kết quả tính theo phương pháp so sánh và phần mềm phần tử hữu hạn chuyên dụng luôn duy trì dưới mức 1,2%. Tại các tiết diện có mômen uốn cực đại, độ chênh lệch chỉ dao động trong khoảng từ 0,3% đến 0,8%.

Kết luận

• Luận văn hệ thống hóa toàn diện 4 phương pháp xây dựng và 5 phương pháp giải kinh điển trong cơ học kết cấu công trình.
• Mở rộng thành công Nguyên lý cực trị Gauss từ cơ hệ chất điểm sang cơ hệ vật rắn biến dạng chịu uốn.
• Thiết lập quy trình tính toán nội lực và độ võng cho dầm đơn, dầm liên tục và khung phẳng bằng phương pháp so sánh.
• Kiểm chứng độ tin cậy vượt trội với sai số tính toán dưới 1,2% so với phần mềm phần tử hữu hạn hiện đại.
• Mở ra hướng đi mới trong việc tối ưu hóa giải thuật phân tích kết cấu tự động cho ngành xây dựng.

Kế hoạch phát triển tiếp theo là xây dựng thư viện mã nguồn mở cho thuật toán so sánh trong vòng 6 tháng tới và tích hợp vào các nền tảng thiết kế tự động hóa trong giai đoạn 2026-2027.

Quý độc giả, kỹ sư kết cấu và các nhà nghiên cứu quan tâm có thể tra cứu toàn văn luận văn tại thư viện Trường Đại học Dân lập Hải Phòng hoặc liên hệ tác giả để trao đổi chuyên môn và ứng dụng giải thuật vào các dự án công trình thực tế ngay hôm nay.

Trích đoạn nội dung tài liệu

chương I đã trình bay bốn đường lối xây dựng bài toán cơ học và các phương pháp giải hiện nay thường dùng trong các giáo trình, tài liệu trong và ngoài nước. Khác với chương l, chương này trình bảy nguyên ly Gauss, sau đó trình bày phương pháp mới dựa trên nguyên lý cực trị Gauss để xây dựng và giải các bài toán cơ học dưới dạng tổng quát, chủ yếu là của cơ hệ vật rắn biến dạng. Để đạt mục tiêu trên, trong chương còn giới thiệu các khái niệm ứng suất và biến dạng của cơ hệ môi trường liên tục và của cơ học kết cấu. Cuối cùng, để làm ví dụ, trình bày việc áp dụng phương pháp mới để nhận được các phương trình vi phân cân bằng của cơ hệ.

Nguyên lý cực trị Gauss Năm 1829 nhà toán học người Đức K. Gauss đã đưa ra nguyên lý sau đây đối với cơ hệ chat diém [1,tr. 171]: “Chuyển động thực của hệ chất điểm có liên kết tùy ý chịu tác động bat ki ở mỗi thời điểm xảy ra một cách phù hợp nhất có thể với chuyển động của hệ đó khi hoàn toàn tự do, nghĩa là chuyển động thực xảy ra với lượng cưỡng bức tối thiểu nếu như số đo lượng cưỡng bức lấy bằng tổng các tích khối lượng chất điểm với bình phương độ lệch vị trí chất điểm so với vị trí khi chúng hoàn toàn tự do”. Goi m, là khối lượng chất điểm, 4; là vị trí của nó, Ö, là vị trí sau thời đoạn vô cùng bé do tác động lực ngoài và do vận tốc ở đầu thời đoạn gây ra, Cc, là vị trí có thể ( bị ràng buộc bởi liên kết) thì lượng cưỡng bức được viết như sau: Z=>m(BC) > Min (2.1) Dấu tổng trong (2.1) lấy theo số chất điểm.

15 Sử dụng nguyên lý vận tốc ảo và nguyên lý D “Alembert, xét hệ ở trạng thái cân bằng và cho rằng có lực với độ lớn tỉ lệ với độ dài Z,C, tác dụng theo chiều từ C, đến B,, Gauss đã chứng minh nguyên lý cua minh [1,tr. Để có thê sử dụng nguyên lý Gauss cần biết đại lượng biến phân của nó. 889], Gibbs (nam 1879) va Appell (năm 1899) đi từ các lập luận khác nhau đều nhận được nguyên lý Gauss và chỉ ra rằng đại lượng biến phân của nguyên lý này là gia tốc. Điều này có nghĩa là: ôn =0; ðz¡=0; ðï;z 0 (2.2) ở đây ö là kí hiệu biến phân ( lấy vi phân khi có định thời gian ), r,, #¡ và ?; lần lượt là vectơ toạ độ, vectơ vận tốc và vectơ gia tốc của điểm i.

Chuyén dich của chất điểm của hệ có liên kết đưới tác dụng của luc F; sau thời đoạn dt tính theo công thức sau đây: r na (2.3) Vì õr;=0 và z¡=0 nên chuyền dịch của chất điểm hoàn toàn tự do (có thể hình dung ở đầu thời đoạn dt liên kết được giải phóng nhưng vẫn giữ lực tác dụng) sau thời đoạn dt là : r, sedt+t iar (2.3) cho ta độ lệch vị trí của chất điểm có liên kết so với vị trí của nó khi hoàn toàn tự do. Có thể xem dt là hằng thì lượng cưỡng bức Z theo (2.1) được viết dưới dạng lực như sau (với độ chính xác bằng thừa số đt/ 4) : m fi z-Ym( 4 -z] —> Min (2.5a) ™ mM, Khi tính lượng cưỡng bức theo (2.5) cần xem gia tốc là đại lượng biến phân (biến phân kiểu Gauss theo cách nói của Boltzmann ). Như vậy, phương pháp tìm cực tiểu của các bài toán cơ học được xây đựng theo nguyên lý (2.5) không thể là bất kỳ mà phải là (khi không có ràng buôc nào khác): OZ oF =0 (2.6) sẽ cho ta phương trình cân bằng. Thật vậy, áp dụng (2.5) ta nhận được phương trình cân bằng của hệ ( ở đây lực tác dụng bằng lực quán tính).

Appell và Boltzmann (năm 1897) còn cho biết nguyên lý Gauss đúng cho hệ liên kết holonom và cả hệ liên kết không holonom [1,tr.5) có dạng của phương pháp bình phương tối thiểu là phương pháp cũng do Gauss đưa ra và được dùng rộng rãi trong toán học hiện đại, trong giải tích cũng như trong lời giải số. Có lẽ vì vậy nguyên lý Gauss thu hut su chú ý của nhiều nhà khoa học, thí đụ, Hertz (năm 1894) dựa trên ý tưởng lượng cưỡng bức đưa ra nguyên lý đường thắng nhất (đường có độ cong nhỏ nhất) hoặc Prigogine (năm 1954) và Gyarmati (năm 1965) đã xây dựng được lượng cưỡng bức của các quá trình không hồi phục trong nhiệt động lực học [2]. Các tài liệu giáo khoa về cơ học thường giới thiệu nguyên lý Gauss đưới dạng (2.5) là dạng dùng được dé tinh toán. Nhưng nguyên lý (2.5) với đại lượng biến phân là gia tốc chỉ là một biểu thị của nguyên lý Gauss (2.1) bởi vì đại lượng biến phân trong cơ học còn có thể là chuyyén vị và vận tốc như trình bảy sau đây.

Phương pháp nguyên lý cực trị Gauss Trong bài viết của mình Gauss nêu nhận xét rằng nguyên lý vận tốc ảo biến vấn đề tĩnh học thành vấn đề toán học thuần tuý, còn nguyên lý D’Alembert dua bai toán động lực học về bài toán tĩnh học và mọi nguyên lý của cơ học hoặc nhiều hoặc ít đều có thể trực tiếp rút ra từ hai nguyên lý trên. Dưới đây trình bày phương pháp dựa trên nguyên lý chuyền vị ảo để nhận được biểu thức (2.1) của nguyên lý Gauss. Xét hệ chất điểm có liên kết tuỳ ý ở một thời điểm bất kì nào đó có nghĩa là phải đưa lực quán tính f; của hệ tại thời điểm đó tác dụng lên hệ. Đối với hệ hoan toan tu do luc quan tinh fo; cua nd bang với ngoại lực (chỉ số “0° ở chân kí tự chỉ rằng kí tự đó thuộc hệ so sánh, trường hợp này là hệ hoàn toàn tự do có cùng khối lượng và cùng chỊu tác dụng lực ngoài giống như hệ có liên kết).

Như vậy, các lực tác dụng lên hệ có liên kết gồm các lực ñ= m¡#; và các lực fo = mị Fø (thay cho ngoại lực). Theo nguyên lý chuyên vị ảo đối với liên kết giữ (liên kết đưới dạng đẳng thức) và không giữ (liên kết dưới dạng bat đăng thức) điều kiện cần và đủ dé hệ ở trạng thái cân bằng la [1,tr.7) cũng được Fourier (năm 1798 ) và Ostrogradsky ( năm 1838) độc lập đưa ra. Có thể nhận xét ngay rằng phần trong ngoặc đơn của (2.7) biểu thị lực tác dụng lên hệ nên phải bằng không để hệ ở trạng thái cân bằng. Trong biểu thức (2.7) cần xem các chuyển vị r¡ độc lập đối với lực tác dụng.

Cho nên từ (2.7) có thể viết: Z=>,- fun > Min (2.8) r¡ là các biến độc lập cần tìm để bảo đảm cho Z cực tiểu. Vì chuyển vị ro của hệ hoàn toàn tự do đã biết nên biểu thức (2.8) tương đương với các 18 biểu thức dưới đây: Z =) (F, — foi) — nụ) > Min (2.8a) i hoặc Z =) maa (4 -%;) > Min (2.8b) m 1 i Dễ dàng nhận thấy (2.8b) là tích của khối lượng m; với bình phương độ lệch vị trí chất điểm và do đó Z.8) là lượng cưỡng bức của nguyên lý Gauss (với độ chính xác bằng thừa số đ/ 2).5), lượng cưỡng bức Z, xác định theo (2.8) biểu thị đầy đủ và rõ ràng tư tưởng của nguyên lý Gauss thể hiện ở chỗ, thứ nhất, nó cho phép so sánh hệ có liên kết với hệ hoàn toàn tự do, thứ hai, đại lượng không biết (đại lượng biến phân) trong (2.8) là chuyền vị giống như trong (2. Cực tiêu của (2.8) cần và phải được tìm từ điều kiện (khi không có các ràng buộc nào khác): OZ or, i =0 (2.9) ap dung vào (2.8) cho ta phương trình cân bằng của cơ hệ. Vi du 1 Vi dụ này lấy từ [3,tr.

Viết phương trình chuyên động của khối lượng m chạy trên đường cong y= bx? trong mặt phẳng (xy), không có lực ma sát, dưới tác dụng của trường gia tốc g (Hình I.1 Các lực tác dụng lên khối lượng m bao gồm: lực quán tính theo chiều y, lực trọng trường theo chiều âm của y, lực quán tính theo x. Chọn hệ so sánh là hệ có 19 cùng khối lượng m nằm trong trường gia tốc g nhưng hoàn toàn tự do. Lượng cưỡng bức được viết theo (2.8) như sau: Z =(mÿ+mg)y+(mĐx -—> Min (a) Thế y=øx? vào (a) ta có Z =(n+mg)bx°+(m3x -> Min (b) Xem chuyên vị x 1a bién độc lập và từ điều kiện = =0 nhận được: X 2bxÿ + 2bgx + š =0 (c) Thay ÿ = 2øxïr+25+” vào (c) nhận được phương trình chuyển động của khối lượng m (4b?x? +1)#+4b?xš? + 2bgx =0 (d) Phương trình (đ) là kết quả cần tìm. Như nhận xét của Gauss nêu trên, có thể nói biéu thức (2.7) đã biến van dé tĩnh học (cân bằng lực) thành vấn đề toán học thuần tuý.

Thật vậy, nếu ta dùng gia tốc là đại lượng biến phân thì tương tự như (2.10) với điều kiện gia tốc ?¡là đại lượng độc lập đối với lực tác dụng.10) có thể viết Z=}È, -f,)?i¡ > Min (2.11) cần xem gia tốc z;¡ là đại lượng biến phân đề bảo đảm cho Z cực tiểu. Vì gia tốc Zo¡ của hệ hoàn toàn tự đo đã biết nên biểu thức (2.11) tương đương với các biểu thức dưới đây: Z=}È, (ƒ-f„)(r Fo) = > Min (2.1 Ib) trùng với (2. Các gia tốc # phải thỏa mãn các liên kết nếu có và điều kiện cực tiểu của (2.11) là biểu thức (2. Làm lại ví dụ 1 (Hình 1) theo nguyên lí (2.5) hoặc biểu thức (2.11) Khối lượng m vừa chuyển động theo x, vừa chuyển động theo y, nhưng do có liên kết y= bx? nên chỉ có một bậc tự do, thí dụ là x.

Các lực tác đụng lên m bao gồm: Lực quán tính theo chiều y, lực trọng trường theo chiều âm của y, lực quán tính theo x. Lượng cưỡng bức Z viét theo (2.5) là: Z= mS +3) +m? > Min (a) Lấy đạo hàm ràng buộc y=bx2 theo thời gian hai lần ta có : 3= 2bx# + 2bx? (b) Thay ÿ trong (a) bằng (b), nhận được Z= (g + 2bx¥ + 2b%?)? +¥2 > Min (c) Xem gia tốc š là biến độc lập và từ điều kiện 2Z/Øï=0 ta có phương trình chuyền động của khối lượng m như sau : (4b?x? +1)#+4b”xš” + 2bgx =0 (d) Phương trình (đ) là kết quả cần tìm. Tương tự, cũng có thể dùng vận tốc ¿ là đại lượng biến phân, khi đó lượng cưỡng bức Z. được viết : Z =),(ƒ-f,) => Min (2.12) với điều kiện vận tốc ” là biến độc lập và thoả mãn các liên kết nếu có.

Trong trường hợp này điều kiện cực tiểu của nguyên lý(2.