Tính toán dao động riêng của dây văng theo phương pháp nguyên lý cực trị gauss (tt)

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ XÂY DỰNG TRƯỜNG ĐẠI HỌC KIẾN TRÚC HÀ NỘI NGUYỄN THỊ THU HIỀN TÍNH TỐN DAO ĐỘNG RIÊNG CỦA DÂY VĂNG THEO PHƯƠNG PHÁP NGUYÊN LÝ CỰC TRỊ GAUSS LUẬN VĂN THẠC SỸ KỸ THUẬT XÂY DỰNG CƠNG TRÌNH DÂN DỤNG & CƠNG NGHIỆP HÀ NỘI – 04/2018 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ XÂY DỰNG TRƯỜNG ĐẠI HỌC KIẾN TRÚC HÀ NỘI NGUYỄN THỊ THU HIỀN kho¸ 2016-2018; líp ch16x2 TÍNH TỐN DAO ĐỘNG RIÊNG CỦA DÂY VĂNG THEO PHƯƠNG PHÁP NGUYÊN LÝ CỰC TRỊ GAUSS LUẬN VĂN THẠC SỸ CHUYÊN NGÀNH: KỸ THUẬT XÂY DỰNG CƠNG TRÌNH DD&CN MÃ SỐ: 62.58.02.08 NGƯỜI HƯỚNG DẪN KHOA HỌC: TS PHẠM VĂN TRUNG HÀ NỘI – 10/2017 i MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN i DANH MỤC KÝ HIỆU TRONG LUẬN VĂN i DANH MỤC BẢNG BIỂU TRONG LUẬN VĂN i DANH MỤC HÌNH VẼ TRONG LUẬN VĂN i MỞ ĐẦU *Lý chọn đề tài *Mục đích nghiên cứu đề tài *Đối tượng nghiên cứu, phạm vi nghiên cứu *Phương pháp nghiên cứu .2 *Ý nghĩa khoa học thực tiễn đề tài *Kết cấu luận văn: NỘI DUNG CHƢƠNG I TỔNG QUAN VỀ HỆ KẾT CẤU DÂY VĂNG VÀ CÁC LÝ THUYẾT TÍNH TỐN DÂY VĂNG 1.1 Phạm vi áp dụng hệ dây văng [1,2,3] .4 1.1.1 Trong lĩnh vực xây dựng dân dụng & công nghiệp 1.1.2 Trong lĩnh vực xây dựng giao thông vận tải 1.2 Cấu tạo hệ kết cấu dây văng .8 1.2.1 Dây văng 1.2.2 Cột tháp trụ 10 1.2.3 Kết cấu dây văng treo đỡ 11 1.2.4 Kết cấu khối neo giữ 12 1.3 Các lý thuyết tính tốn tĩnh dây văng .14 1.3.1 Phương pháp tính dây theo hai trạng thái 14 1.3.2 Phương pháp tính dây theo trạng thái .18 1.3.3 Phương pháp dây xích .20 1.4 Các lý thuyết tính tốn dao động riêng dây văng 21 ii 1.4.1 Phương pháp lý thuyết 21 1.4.2 Phương pháp thực nghiệm 22 1.5 Phương pháp nguyên lý cực trị Gauss 25 1.5.1 Nguyên lý cực trị Gauss 25 1.5.2 Phương pháp nguyên lý cực trị Gauss 27 CHƢƠNG XÂY DỰNG VÀ GIẢI BÀI TOÁN DAO ĐỘNG RIÊNG DÂY VĂNG 36 2.1 Tóm tắt tốn động lực học cơng trình 36 2.2 Xây dựng giải toán dây văng chịu tác động tĩnh 37 2.2.1 Xây dựng toán 37 2.2.2 Giải toán 41 2.2.3 Thuật tốn tính tốn: 43 2.2.4 Sơ đồ khối tính tốn: 44 2.3 Xây dựng giải toán dao động riêng dây văng 45 2.3.1 Xây dựng toán 45 2.3.2 Giải toán 46 2.3.3 Thuật toán tính tốn: 48 2.3.4 Sơ đồ khối tính tốn: 50 CHƢƠNG CÁC VÍ DỤ TÍNH TỐN 51 3.1 Bài toán dây văng chịu tác động tĩnh: 51 3.1.1 Dây văng chịu tải trọng thân 51 3.1.2 Dây văng chịu tải trọng thân áp lực gió .53 3.1.3 Dây văng chịu tải trọng thân, áp lực gió biến dạng neo .55 3.2 Bài toán dao động riêng dây văng: 57 3.2.1 Đặt toán 57 3.2.2 Kết ví dụ tính tốn 59 KẾT LUẬN VÀ KIẾN NGHỊ 65 Kết luận 65 Kiến nghị .65 PHỤ LỤC i iii i LỜI CAM ĐOAN Tơi xin cam đoan cơng trình nghiên cứu riêng Các số liệu, kết nêu luận văn trung thực chưa cơng bố cơng trình khác Hà nội, tháng năm 2017 Tác giả luận văn Nguyễn Thị Thu Hiền i LỜI CẢM ƠN Tác giả luận văn xin bày tỏ lòng biết ơn chân thành TS Phạm Văn Trung tận tình giúp đỡ cho nhiều dẫn Khoa học có giá trị thường xuyên động viên, tạo điều kiện, giúp đỡ tác giả hoàn thành luận văn nâng cao lực Khoa học tác giả Tác giả xin trân trọng cảm ơn Thầy cô nhà Khoa học trường quan tâm góp ý làm cho luận văn hoàn thiện Tác giả xin trân trọng cảm ơn Thầy, Cô Giáo, văn Bộ môn Sức bền Cơ học kết cấu, Khoa Xây dựng, Khoa Sau Đại học, Trường Đại học Kiến trúc Hà Nội đồng nghiệp tạo điều kiện thuận lợi hợp tác trình nghiên cứu Tác giả luận văn Nguyễn Thị Thu Hiền i DANH MỤC KÝ HIỆU TRONG LUẬN VĂN Ký hiệu Đại lƣợng Đơn vị A Diện tích tiết diện m2 E Mơ dun đàn hồi kháng kéo nén t/m2 F Phiếm hàm F Độ vồng vòm m G Mơ đun đàn hồi kháng cắt t/m2 I Mơ men qn tính tiết diện m4 L Nhịp vòm m M Mơ men vòm t.m N Lực dọc t P Lực tập trung t Q Lực cắt t A Diện tích tiết diện m2 E Mơ dun đàn hồi kháng kéo nén t/m2 F Phiếm hàm F Độ vồng vòm m G Mơ đun đàn hồi kháng cắt t/m2 I Mơ men qn tính tiết diện m4 i DANH MỤC BẢNG BIỂU TRONG LUẬN VĂN Bảng số Nội dung Bảng 3.1 Lực căng cáp chuyển vị điểm nút dây văng chịu tải thân Bảng 3.2 Lực căng cáp chuyển vị điểm nút dây văng chịu tải thân áp lực gió Bảng 3.3 Lực căng cáp chuyển vị điểm nút dây văng chịu tải thân áp lực gió biến dạng neo Bảng 3.4 Lực căng cáp chuyển vị điểm nút dây văng chịu tải thân áp lực gió biến dạng neo biến thiên nhiệt độ i DANH MỤC HÌNH VẼ TRONG LUẬN VĂN Hình số Hình 1.1 Hình 1.2 Hình 1.3 Hình 1.4 Hình 1.5 Hình 1.6 Hình 1.7 Hình 1.8 Hình 1.9 Hình 1.10 Hình 1.11 Hình 1.12 Hình 1.13 Hình 1.14 Hình 1.15 Hình 1.16 Hình 2.1 Hình 2.2 Hình 2.3 Hình 2.4 Hình 3.1 Hình 3.2 Hình 3.3 Hình 3.4 Hình 3.5 Nội dung Sân vận động quốc gia Mỹ đình Sân vận động quốc gia Jeju – Hàn Quốc Cầu Trần Thị Lý, Đà Nẵng Cầu Gia Thiệu, Trung Quốc Các cấu kiện hệ kết cấu dây văng Các loại hình dáng cột tháp trụ Mái sân vận động mỹ đình treo dây văng18 Neo cáp vào đất mố neo Hệ thống neo trọng lực Một số chi tiết cấu tạo khối neo Trạng thái ban đầu dây Trạng thái tính tốn dây Sơ đồ dây tính theo phương pháp trạng thái Sơ đồ phân tích đoạn dây Mơ hình phân tích dao động dây văng Nội lực phân tố Sơ đồ tính toán dây cáp Sơ đồ khối toán tĩnh Sơ đồ tính dao động riêng dây Sơ đồ khối tốn tìm trị riêng Sơ đồ tính dây (vẽ matlab) Dạng ban đầu chuyển vị dây văng chịu tải Dạng ban đầu chuyển vị dây văng chịu tải gió Dạng ban đầu chuyển vị dây văng chịu tải, gió, CV Dạng ban đầu chuyển vị dây văng chịu tải, gió, CV Hình 3.6 Hình 3.7 Hình 3.8 biến thiên nhiệt độ Sơ đồ tính dao động riêng dây Dạng dao động riêng mặt phẳng dây Dạng dao động riêng mặt phẳng dây b04)^2)^(1/2)))/(((a03 - a04)^2 + (b03 - b04)^2)^(1/2)*((x03 - x04)^2 + (y03 - y04)^2 + (z03 - z04)^2)^(1/2))) pz + EA*(((2*z03 - 2*z04)*(((a03 - a04)^2 + (b03 - b04)^2)^(1/2) - ((x03 - x04)^2 + (y03 - y04)^2 + (z03 - z04)^2)^(1/2) + a*t*((a03 - a04)^2 + (b03 b04)^2)^(1/2)))/(((a03 - a04)^2 + (b03 - b04)^2)^(1/2)*((x03 - x04)^2 + (y03 - y04)^2 + (z03 - z04)^2)^(1/2)) - ((2*z04 - 2*z05)*(((a04 - a05)^2 + (b04 - b05)^2)^(1/2) ((x04 - x05)^2 + (y04 - y05)^2 + (z04 - z05)^2)^(1/2) + a*t*((a04 - a05)^2 + (b04 b05)^2)^(1/2)))/(((a04 - a05)^2 + (b04 - b05)^2)^(1/2)*((x04 - x05)^2 + (y04 - y05)^2 + (z04 - z05)^2)^(1/2))) pz + EA*(((2*z04 - 2*z05)*(((a04 - a05)^2 + (b04 - b05)^2)^(1/2) - ((x04 - x05)^2 + (y04 - y05)^2 + (z04 - z05)^2)^(1/2) + a*t*((a04 - a05)^2 + (b04 b05)^2)^(1/2)))/(((a04 - a05)^2 + (b04 - b05)^2)^(1/2)*((x04 - x05)^2 + (y04 - y05)^2 + (z04 - z05)^2)^(1/2)) - ((2*z05 - 2*z06)*(((a05 - a06)^2 + (b05 - b06)^2)^(1/2) ((x05 - x06)^2 + (y05 - y06)^2 + (z05 - z06)^2)^(1/2) + a*t*((a05 - a06)^2 + (b05 b06)^2)^(1/2)))/(((a05 - a06)^2 + (b05 - b06)^2)^(1/2)*((x05 - x06)^2 + (y05 - y06)^2 + (z05 - z06)^2)^(1/2))) pz + EA*(((2*z05 - 2*z06)*(((a05 - a06)^2 + (b05 - b06)^2)^(1/2) - ((x05 - x06)^2 + (y05 - y06)^2 + (z05 - z06)^2)^(1/2) + a*t*((a05 - a06)^2 + (b05 b06)^2)^(1/2)))/(((a05 - a06)^2 + (b05 - b06)^2)^(1/2)*((x05 - x06)^2 + (y05 - y06)^2 + (z05 - z06)^2)^(1/2)) - ((2*z06 - 2*z07)*(((a06 - a07)^2 + (b06 - b07)^2)^(1/2) ((x06 - x07)^2 + (y06 - y07)^2 + (z06 - z07)^2)^(1/2) + a*t*((a06 - a07)^2 + (b06 b07)^2)^(1/2)))/(((a06 - a07)^2 + (b06 - b07)^2)^(1/2)*((x06 - x07)^2 + (y06 - y07)^2 + (z06 - z07)^2)^(1/2))) pz + EA*(((2*z06 - 2*z07)*(((a06 - a07)^2 + (b06 - b07)^2)^(1/2) - ((x06 - x07)^2 + (y06 - y07)^2 + (z06 - z07)^2)^(1/2) + a*t*((a06 - a07)^2 + (b06 b07)^2)^(1/2)))/(((a06 - a07)^2 + (b06 - b07)^2)^(1/2)*((x06 - x07)^2 + (y06 - y07)^2 + (z06 - z07)^2)^(1/2)) - ((2*z07 - 2*z08)*(((a07 - a08)^2 + (b07 - b08)^2)^(1/2) - ((x07 - x08)^2 + (y07 - y08)^2 + (z07 - z08)^2)^(1/2) + a*t*((a07 - a08)^2 + (b07 b08)^2)^(1/2)))/(((a07 - a08)^2 + (b07 - b08)^2)^(1/2)*((x07 - x08)^2 + (y07 - y08)^2 + (z07 - z08)^2)^(1/2))) pz + EA*(((2*z07 - 2*z08)*(((a07 - a08)^2 + (b07 - b08)^2)^(1/2) - ((x07 - x08)^2 + (y07 - y08)^2 + (z07 - z08)^2)^(1/2) + a*t*((a07 - a08)^2 + (b07 b08)^2)^(1/2)))/(((a07 - a08)^2 + (b07 - b08)^2)^(1/2)*((x07 - x08)^2 + (y07 - y08)^2 + (z07 - z08)^2)^(1/2)) - ((2*z08 - 2*z09)*(((a08 - a09)^2 + (b08 - b09)^2)^(1/2) ((x08 - x09)^2 + (y08 - y09)^2 + (z08 - z09)^2)^(1/2) + a*t*((a08 - a09)^2 + (b08 b09)^2)^(1/2)))/(((a08 - a09)^2 + (b08 - b09)^2)^(1/2)*((x08 - x09)^2 + (y08 - y09)^2 + (z08 - z09)^2)^(1/2))) pz + EA*(((2*z08 - 2*z09)*(((a08 - a09)^2 + (b08 - b09)^2)^(1/2) - ((x08 - x09)^2 + (y08 - y09)^2 + (z08 - z09)^2)^(1/2) + a*t*((a08 - a09)^2 + (b08 b09)^2)^(1/2)))/(((a08 - a09)^2 + (b08 - b09)^2)^(1/2)*((x08 - x09)^2 + (y08 - y09)^2 + (z08 - z09)^2)^(1/2)) - ((2*z09 - 2*z10)*(((a09 - a10)^2 + (b09 - b10)^2)^(1/2) ((x09 - x10)^2 + (y09 - y10)^2 + (z09 - z10)^2)^(1/2) + a*t*((a09 - a10)^2 + (b09 b10)^2)^(1/2)))/(((a09 - a10)^2 + (b09 - b10)^2)^(1/2)*((x09 - x10)^2 + (y09 - y10)^2 + (z09 - z10)^2)^(1/2))) pz + EA*(((2*z09 - 2*z10)*(((a09 - a10)^2 + (b09 - b10)^2)^(1/2) - ((x09 - x10)^2 + (y09 - y10)^2 + (z09 - z10)^2)^(1/2) + a*t*((a09 - a10)^2 + (b09 b10)^2)^(1/2)))/(((a09 - a10)^2 + (b09 - b10)^2)^(1/2)*((x09 - x10)^2 + (y09 - y10)^2 + (z09 - z10)^2)^(1/2)) - ((2*z10 - 2*z11)*(((a10 - a11)^2 + (b10 - b11)^2)^(1/2) ((x10 - x11)^2 + (y10 - y11)^2 + (z10 - z11)^2)^(1/2) + a*t*((a10 - a11)^2 + (b10 b11)^2)^(1/2)))/(((a10 - a11)^2 + (b10 - b11)^2)^(1/2)*((x10 - x11)^2 + (y10 - y11)^2 + (z10 - z11)^2)^(1/2))) pz + EA*(((2*z10 - 2*z11)*(((a10 - a11)^2 + (b10 - b11)^2)^(1/2) - ((x10 - x11)^2 + (y10 - y11)^2 + (z10 - z11)^2)^(1/2) + a*t*((a10 - a11)^2 + (b10 b11)^2)^(1/2)))/(((a10 - a11)^2 + (b10 - b11)^2)^(1/2)*((x10 - x11)^2 + (y10 - y11)^2 + (z10 - z11)^2)^(1/2)) - ((2*z11 - 2*z12)*(((a11 - a12)^2 + (b11 - b12)^2)^(1/2) ((x11 - x12)^2 + (y11 - y12)^2 + (z11 - z12)^2)^(1/2) + a*t*((a11 - a12)^2 + (b11 b12)^2)^(1/2)))/(((a11 - a12)^2 + (b11 - b12)^2)^(1/2)*((x11 - x12)^2 + (y11 - y12)^2 + (z11 - z12)^2)^(1/2))) pz + EA*(((2*z11 - 2*z12)*(((a11 - a12)^2 + (b11 - b12)^2)^(1/2) - ((x11 - x12)^2 + (y11 - y12)^2 + (z11 - z12)^2)^(1/2) + a*t*((a11 - a12)^2 + (b11 b12)^2)^(1/2)))/(((a11 - a12)^2 + (b11 - b12)^2)^(1/2)*((x11 - x12)^2 + (y11 - y12)^2 + (z11 - z12)^2)^(1/2)) - ((2*z12 - 2*z13)*(((a12 - a13)^2 + (b12 - b13)^2)^(1/2) ((x12 - x13)^2 + (y12 - y13)^2 + (z12 - z13)^2)^(1/2) + a*t*((a12 - a13)^2 + (b12 b13)^2)^(1/2)))/(((a12 - a13)^2 + (b12 - b13)^2)^(1/2)*((x12 - x13)^2 + (y12 - y13)^2 + (z12 - z13)^2)^(1/2))) pz + EA*(((2*z12 - 2*z13)*(((a12 - a13)^2 + (b12 - b13)^2)^(1/2) - ((x12 - x13)^2 + (y12 - y13)^2 + (z12 - z13)^2)^(1/2) + a*t*((a12 - a13)^2 + (b12 b13)^2)^(1/2)))/(((a12 - a13)^2 + (b12 - b13)^2)^(1/2)*((x12 - x13)^2 + (y12 - y13)^2 + (z12 - z13)^2)^(1/2)) - ((2*z13 - 2*z14)*(((a13 - a14)^2 + (b13 - b14)^2)^(1/2) ((x13 - x14)^2 + (y13 - y14)^2 + (z13 - z14)^2)^(1/2) + a*t*((a13 - a14)^2 + (b13 b14)^2)^(1/2)))/(((a13 - a14)^2 + (b13 - b14)^2)^(1/2)*((x13 - x14)^2 + (y13 - y14)^2 + (z13 - z14)^2)^(1/2))) pz + EA*(((2*z13 - 2*z14)*(((a13 - a14)^2 + (b13 - b14)^2)^(1/2) - ((x13 - x14)^2 + (y13 - y14)^2 + (z13 - z14)^2)^(1/2) + a*t*((a13 - a14)^2 + (b13 b14)^2)^(1/2)))/(((a13 - a14)^2 + (b13 - b14)^2)^(1/2)*((x13 - x14)^2 + (y13 - y14)^2 + (z13 - z14)^2)^(1/2)) - ((2*z14 - 2*z15)*(((a14 - a15)^2 + (b14 - b15)^2)^(1/2) ((x14 - x15)^2 + (y14 - y15)^2 + (z14 - z15)^2)^(1/2) + a*t*((a14 - a15)^2 + (b14 b15)^2)^(1/2)))/(((a14 - a15)^2 + (b14 - b15)^2)^(1/2)*((x14 - x15)^2 + (y14 - y15)^2 + (z14 - z15)^2)^(1/2))) pz + EA*(((2*z14 - 2*z15)*(((a14 - a15)^2 + (b14 - b15)^2)^(1/2) - ((x14 - x15)^2 + (y14 - y15)^2 + (z14 - z15)^2)^(1/2) + a*t*((a14 - a15)^2 + (b14 - b15)^2)^(1/2)))/(((a14 - a15)^2 + (b14 - b15)^2)^(1/2)*((x14 - x15)^2 + (y14 - y15)^2 + (z14 - z15)^2)^(1/2)) - ((2*z15 - 2*z16)*(((a15 - a16)^2 + (b15 - b16)^2)^(1/2) ((x15 - x16)^2 + (y15 - y16)^2 + (z15 - z16)^2)^(1/2) + a*t*((a15 - a16)^2 + (b15 b16)^2)^(1/2)))/(((a15 - a16)^2 + (b15 - b16)^2)^(1/2)*((x15 - x16)^2 + (y15 - y16)^2 + (z15 - z16)^2)^(1/2))) pz + EA*(((2*z15 - 2*z16)*(((a15 - a16)^2 + (b15 - b16)^2)^(1/2) - ((x15 - x16)^2 + (y15 - y16)^2 + (z15 - z16)^2)^(1/2) + a*t*((a15 - a16)^2 + (b15 b16)^2)^(1/2)))/(((a15 - a16)^2 + (b15 - b16)^2)^(1/2)*((x15 - x16)^2 + (y15 - y16)^2 + (z15 - z16)^2)^(1/2)) - ((2*z16 - 2*z17)*(((a16 - a17)^2 + (b16 - b17)^2)^(1/2) ((x16 - x17)^2 + (y16 - y17)^2 + (z16 - z17)^2)^(1/2) + a*t*((a16 - a17)^2 + (b16 b17)^2)^(1/2)))/(((a16 - a17)^2 + (b16 - b17)^2)^(1/2)*((x16 - x17)^2 + (y16 - y17)^2 + (z16 - z17)^2)^(1/2))) pz + EA*(((2*z16 - 2*z17)*(((a16 - a17)^2 + (b16 - b17)^2)^(1/2) - ((x16 - x17)^2 + (y16 - y17)^2 + (z16 - z17)^2)^(1/2) + a*t*((a16 - a17)^2 + (b16 b17)^2)^(1/2)))/(((a16 - a17)^2 + (b16 - b17)^2)^(1/2)*((x16 - x17)^2 + (y16 - y17)^2 + (z16 - z17)^2)^(1/2)) - ((2*z17 - 2*z18)*(((a17 - a18)^2 + (b17 - b18)^2)^(1/2) ((x17 - x18)^2 + (y17 - y18)^2 + (z17 - z18)^2)^(1/2) + a*t*((a17 - a18)^2 + (b17 b18)^2)^(1/2)))/(((a17 - a18)^2 + (b17 - b18)^2)^(1/2)*((x17 - x18)^2 + (y17 - y18)^2 + (z17 - z18)^2)^(1/2))) pz + EA*(((2*z17 - 2*z18)*(((a17 - a18)^2 + (b17 - b18)^2)^(1/2) - ((x17 - x18)^2 + (y17 - y18)^2 + (z17 - z18)^2)^(1/2) + a*t*((a17 - a18)^2 + (b17 b18)^2)^(1/2)))/(((a17 - a18)^2 + (b17 - b18)^2)^(1/2)*((x17 - x18)^2 + (y17 - y18)^2 + (z17 - z18)^2)^(1/2)) - ((2*z18 - 2*z19)*(((a18 - a19)^2 + (b18 - b19)^2)^(1/2) ((x18 - x19)^2 + (y18 - y19)^2 + (z18 - z19)^2)^(1/2) + a*t*((a18 - a19)^2 + (b18 b19)^2)^(1/2)))/(((a18 - a19)^2 + (b18 - b19)^2)^(1/2)*((x18 - x19)^2 + (y18 - y19)^2 + (z18 - z19)^2)^(1/2))) pz + EA*(((2*z18 - 2*z19)*(((a18 - a19)^2 + (b18 - b19)^2)^(1/2) - ((x18 - x19)^2 + (y18 - y19)^2 + (z18 - z19)^2)^(1/2) + a*t*((a18 - a19)^2 + (b18 b19)^2)^(1/2)))/(((a18 - a19)^2 + (b18 - b19)^2)^(1/2)*((x18 - x19)^2 + (y18 - y19)^2 + (z18 - z19)^2)^(1/2)) - ((2*z19 - 2*z20)*(((a19 - a20)^2 + (b19 - b20)^2)^(1/2) ((x19 - x20)^2 + (y19 - y20)^2 + (z19 - z20)^2)^(1/2) + a*t*((a19 - a20)^2 + (b19 b20)^2)^(1/2)))/(((a19 - a20)^2 + (b19 - b20)^2)^(1/2)*((x19 - x20)^2 + (y19 - y20)^2 + (z19 - z20)^2)^(1/2)))]; hien.m x0=[5.000 10.000 15.000 20.000 25.000 30.000 35.000 40.000 45.000 50.000 55.000 60.000 65.000 70.000 75.000 80.000 85.000 90.000 95.000 0.542 1.500 2.875 4.667 6.875 9.500 12.542 16.000 19.875 24.167 28.875 34.000 39.542 45.500 51.875 58.667 65.875 73.500 81.542 0 0 0 0 0 0 0 0 0 0]; % Make a starting guess at the solution options=optimset('Display','iter'); % Option to display output [x,fval] = fsolve('myfun',x0,options); % call optimizer x; LUCCANG=[EA*(s01 - l01 - a*t*l01 )/l01 EA*(s02 - l02 - a*t*l02 )/l02 EA*(s03 - l03 - a*t*l03 )/l03 EA*(s04 - l04 - a*t*l04 )/l04 EA*(s05 - l05 - a*t*l05 )/l05 EA*(s06 - l06 - a*t*l06 )/l06 EA*(s07 - l07 - a*t*l07 )/l07 EA*(s08 - l08 - a*t*l08 )/l08 EA*(s09 - l09 - a*t*l09 )/l09 EA*(s10 - l10 - a*t*l10 )/l10 EA*(s11 - l11 - a*t*l11 )/l11 EA*(s12 - l12 - a*t*l12 )/l12 EA*(s13 - l13 - a*t*l13 )/l13 EA*(s14 - l14 - a*t*l14 )/l14 EA*(s15 - l15 - a*t*l15 )/l15 EA*(s16 - l16 - a*t*l16 )/l16 EA*(s17 - l17 - a*t*l17 )/l17 EA*(s18 - l18 - a*t*l18 )/l18 EA*(s19 - l19 - a*t*l19 )/l19 EA*(s20 - l20 - a*t*l20 )/l20] cvx=[x00 - x01 - a01 x02 - a02 x03 - a03 x04 - a04 x05 - a05 x06 - a06 x07 - a07 x08 - a08 x09 - a09 x10 - a10 x11 - a11 x12 - a12 x13 - a13 x14 - a14 x15 - a15 x16 - a16 x17 - a17 x18 - a18 a00 x19 - a19 x20 - a20] cvy=[y00 - y01 - b01 y02 - b02 y03 - b03 y04 - b04 y05 - b05 y06 - b06 y07 - b07 y08 - b08 y09 - b09 y10 - b10 y11 - b11 y12 - b12 y13 - b13 y14 - b14 y15 - b15 y16 - b16 y17 - b17 y18 - b18 y19 - b19 y20 - b20] cvz=[z00 z01 z02 z03 b00 z04 z05 z06 z07 z08 z09 z10 z11 z12 z13 z14 z15 z16 z17 z18 z19 z20] X=[x00 x01 x02 x03 x04 x05 x06 x07 x08 x09 x10 x11 x12 x13 x14 x15 x16 x17 x18 x19 x20]; Y=[y00 y01 y02 y03 y04 y05 y06 y07 y08 y09 y10 y11 y12 y13 y14 y15 y16 y17 y18 y19 y20]; Z=[z00 z01 z02 z03 z04 z05 z06 z07 z08 z09 z10 z11 z12 z13 z14 z15 z16 z17 z18 z19 z20]; A=[a00 a01 a02 a03 a04 a05 a06 a07 a08 a09 a10 a11 a12 a13 a14 a15 a16 a17 a18 a19 a20]; B=[b00 b01 b02 b03 b04 b05 b06 b07 b08 b09 b10 b11 b12 b13 b14 b15 b16 b17 b18 b19 b20]; C=[0 0 0 0 0 0 0 0 0 0 0]; plot3(A, C, B, X, Z, Y) ... hạn chế Phương pháp nguyên lý cực trị Gauss phương pháp GS TSKH Hà Huy Cương trình bày dựa nguyên lý chuyển vị ảo để nhận biểu thức nguyên lý Gauss gọi phương pháp nguyên lý cực trị Gauss Nhiều... dao động riêng dây văng 21 ii 1.4.1 Phương pháp lý thuyết 21 1.4.2 Phương pháp thực nghiệm 22 1.5 Phương pháp nguyên lý cực trị Gauss 25 1.5.1 Nguyên lý cực trị Gauss ... dây văng .14 1.3.1 Phương pháp tính dây theo hai trạng thái 14 1.3.2 Phương pháp tính dây theo trạng thái .18 1.3.3 Phương pháp dây xích .20 1.4 Các lý thuyết tính tốn dao động