Phân tích ảnh hưởng của từ biến đến cầu treo dây văng

ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA CAO BÁ DỖN PHÂN TÍCH ẢNH HƯỞNG CỦA TỪ BIẾN ĐẾN CẦU TREO DÂY VĂNG CHUYÊN NGÀNH: XÂY DỰNG CẦU, HẦM Mã số: 60 58 25 LUẬN VĂN THẠC SĨ TP HỒ CHÍ MINH, THÁNG 07 NĂM 2010 TĨM TẮT Trong cơng việc thiết kế, thơng số khó ước lượng thay đổi tình trang ban đầu cơng trình với thời gian Trong suốt thời gian tồn mình, cơng trình chịu tác động tải trọng khác làm thay đổi hình dạng hình học, chí thay đổi đặc tính học Trong yếu tố ảnh hưởng đến thay đổi cơng trình từ biến yếu tố quan trọng Sự tác động tránh khỏi thời gian lâu Do ảnh hưởng từ biến mà biến dạng cơng trình có khuynh hướng tăng lên tác dụng tải trọng ban đầu qua thời gian, biến dạng cuối đạt trị số gấp hay lần trị số ban đầu Trong suốt trình thiết kế, phát triển biến dạng theo thời gian thật khó để xác định đặc biệt kết cấu siêu tĩnh bậc cao cầu dây văng nơi ứng suất có liên hệ mât thiết tới hình dạng hình học cơng trình Trong cầu dây văng dầm bê tơng, việc chuyển vị dầm thay đổi theo thời gian làm thay đổi lực căng dây ban đầu Điều làm thay đổi đáng kể nội lực dầm Mà nói trên, nguyên nhân chủ yếu tạo nên thay đổi tượng từ biến bê tông Chương luận văn giới thiệu tổng quan cầu dây văng, từ biến việc nghiên cứu từ biến cầu dây văng Chương minh họa giải thích phương pháp xấp xỉ sử dụng để ước lượng ảnh hưởng từ biến đến kết cấu bê tông Phương pháp giới thiệu luận văn “Phương Pháp Mơđun Đàn Hồi Hữu Hiệu” Nó chọn phương pháp khác sử dụng phổ biến xác Trong chương kết hợp phương pháp môđun đàn hồi hữu hiệu phương pháp lực để nghiên cứu kết cấu siêu tĩnh ảnh hưởng từ biến Trong chương lý thuyết trước sử dụng kết hợp với ngơn ngữ lập trình matlab để nghiên cứu thay đổi lực căng dây nội lực cầu treo dây văng Chương cuối đưa kết luận ảnh hưởng từ biến đến cầu treo dây văng phương pháp dung để nghiên cứu ABSTRACT In design of structure, one of the most difficult parameter to assess is how the initial condition of structure will change with time During its life span, a structure can be subjected to different load, changes in geometric configuration and even changes in its mechanical properties All among factors that affect the geometrical reconfiguration(i e settlement, overloads) creep is one of the most important given its inevitability and because of its persistent effects( normally estimated as ten thousand days) Due to effect of creep, a structure tends to amplify its deformation under a given load condition over time, and final deformation can even reach values five or six time grater than the initial ones During the design, the fact that deformations grow with time can be difficult condition especially for highly indeterminate structure like cable stayed bridges where stress relate to geometric configuration of the structure itself In concrete cable stayed bridges, in fact, the increase deformation of the deck and the pylon overtime leads to changes in the initial tension in stay with a obvious difficulty in design phase of the structure The first chapter of this thesis will present generally about cable stayed bridges, creep as well as studying creep effects on cable stayed bridges The second chapter of this thesis will illustrate and explain one approximate method used to estimate the effect of creep on a concrete structure The method proposed in this thesis is the “Age – Adjusted Effective Modulus Method” It was chosen among others because it is one of the most commonly used, and because it is highly accessible In the third chapter, the Age – Adjusted Effective Modulus Method will be used in conjunction with the force method to study indeterminate structure under the effect of creep In the fourth chapter, the theory developed in the previous section will be used to study the changes of tension in the stay in a cable stayed bridges The last chapter of this thesis will give conclusions about the effects of creep on cable stayed bridges and method used to estimate the effect of creep on cable stayed bridges MỤC LỤC Mở đầu…………………………………………………………………………3 Ý nghĩa thực tiễn đề tài……………………………………………… Phạm vi phương pháp nghiên cứu …………………………………… Chương : Tổng quan cầu dây văng từ biến………………………….4 1.1 Tổng quan cầu dây văng………………………………………………4 1.2 Tổng quan từ biến …………………………………………… ……6 1.2.1 Khái niệm……………………………………………………………6 1.2.2 Các mơ hình tính tốn từ biến phổ biến…………………………… Sơ lược nghiên cứu từ biến cầu dây văng …………………8 Chương Nghiên cứu lý thuyết tính tốn từ biến kết cấu bê tông Từ biến kết cấu bê tơng có ứng suất khơng đổi ………………….8 2 Từ biến kết cấu bê tơng có ứng suất thay đổi theo thời gian … 15 Phương pháp bước để tìm hàm chùng kết cấu bê tông…… 20 Phương pháp môđun đàn hồi hữu hiệu theo độ tuổi……………………22 Chương Phân tích ảnh hưởng từ biến đến cầu treo dây văng………23 Thiết lập biểu thức toán……………………………………………… 23 1 Biểu thức toán xem cáp đàn hồi tuyến tính…………………….23 Biểu thức tốn xét đến tượng chùng cáp………………….26 Ví dụ ………………………………………………………………… 29 Chương Xây dựng chương trình tính tốn từ biến cầu dây văng 39 Sơ lược phương pháp phần tử hữu hạn(PTHH)…………………… 39 1 Khái niệm phương pháp PTHH……………………………………39 Trình tự thực toán theo phương pháp PTHH…………… 40 Hàm xấp xỉ phép nội suy phương pháp PTHH……………41 Hàm xấp xỉ ……………………………………………………41 Phép nội suy………………………………………………… 41 3 Dạng đa thức xấp xỉ………………………………………… 42 4 Chọn bậc đa thức xấp xỉ…………………………………… 43 Ma trận hàm dạng …………………………………………….43 Mơ hình tốn ngơn ngữ lập trình matlab…………………… 45 Chương Kết luận……………………………………………………………49 Tài liệu tham khảo .51 MỞ ĐẦU Ý nghĩa thưc tiễn đề tài Trong suốt trình tồn mình, cầu dây văng chịu nhiều loại tải trọng khác thay đổi hình dạng hình học Trong tất yếu tố ảnh hưởng đến thay đổi hình dạng hình học cầu dây văng từ biến yếu tố quan trọng Do ảnh hưởng từ biến, chuyển vị dầm cầu treo dây văng thay đổi làm thay đổi lực căng dây văng nội lực dầm, gây khó khăn cho q trình thiết kế Điều cho thấy việc nghiên cứu ảnh hưởng từ biến đến cầu dây văng bê tông cần thiết, ảnh hưởng trực tiếp đến phương án thiết kế hiệu sử dụng công trình Phạm vi nghiên cứu phương pháp nghiên cứu - Phạm vi nghiên cứu cầu dây văng dầm bê tông giai đoạn khai thác - Phương pháp nghiên cứu + Nghiên cứu lý thuyết từ biến kết cấu bê tơng + Dựa vào ta xây dựng cơng thức tính tốn từ biến cầu dây văng + Mơ hình tốn Matlab CHƯƠNG TỔNG QUAN VỀ CẦU DÂY VĂNG VÀ TỪ BIẾN 1Tổng quan cầu dây văng Cùng với cầu treo, cầu dây văng có lịch sử phát triển lâu đời Năm 1817 ý tưởng Poet thực Anh cho cầu có nhịp 33 5m Hệ dầm mặt cầu đỡ dây văng xuất phát từ đỉnh tháp cầu, phía đối diện bố trí dây neo Thành phần lực ngang dây văng truyền vào tay vịn lan can Đó tiền đề cho hệ cầu thành phần lực ngang truyền cho hệ dầm cứng [9] Với ý tưởng với việc phát triển hoàn thiện hệ dàn dây Giselard theo hướng tạo hệ bất biến hình với dây xiên chịu kéo dầm cứng chịu uốn cầu dây văng áp dụng phát triển Từ cầu đầu tiên, cầu Stromsund xây dựng Thụy Điển năm 1955, đến cầy dây văng đả ứng dụng rông rãi giới đạt thành tựu rực rỡ Cuối kỹ 20 kỹ lục chiều dài nhịp cầu dây văng liên tục bị phá, nhiều cầu trở thành di sản văn hóa biểu tượng kiến trúc, đánh dấu phát triển khoa học kỹ thuật thời đại Cầu Sromsund Thụy Điển Vào kỷ 21 nhiều cầu dây văng nhịp 1000m nghiên cứu áp dụng cần qua vịnh Messine(Italia);cầu qua vịnh Storebely(Đan Mạch) Trong cầu dây văng dầm có độ cứng lón hệ làm việc giống dầm cứng tựa gối đàn hồi dây văng, tăng số lượng gối đàn hồi khối lương dây văng lực nén dầm tăng không đáng kể lại làm giảm đáng kể momen uốn dầm Do cầu dây văng vượt nhịp lớn mà lượng vật liệu tăng khơng đáng kể Do cầu dây văng có ưu điểm nhiều mặt nên nhanh chóng áp dụng tồn giới Tính đa dạng đặc điểm quan trọng cầu dây văng mà hệ thống cầu khác khơng thể có Tính đa dạng thể hiệ chiều dài nhịp, số lượng nhịp khả thay đổi dáng vẻ, sơ đồ phân bố dây số lượng mặt phẳng dây ……mà đảm bảo tiêu kinh tế kỹ thuật tốt Ở nước ta cầu dây văng xây dựng vào năm 1976 Quảng Trị sau vào năm 1999 cầu bi sập gỉ, đứt neo Hiện xây dựng cầu dây văng khác với hỗ trợ nước cầu Mỹ Thuận, cầu qua sông Hàn Cầu Mỹ Thuận Đặc biệt ta khánh thành đưa vào sử dụng cầu dây văng thành phố cầu Phú Mỹ Cầu Phú Mỹ Tổng quan từ biến Khái niệm Từ biến tượng biến dạng cấu kiện tăng dần theo thời gian ảnh hưởng ứng suất không đổi tác dụng lâu dài Nguyên nhân tượng xếp lại cấu trúc vật liệu Giá trị biến dạng từ biến đạt đến lần biến dạng đàn hồi [5] Từ biến gây hiệu ứng có lợi làm giảm tác động cưỡng biến dạng không phận kết cấu Tuy nhiên từ biến gây nhiều hiệu ứng bất lợi tăng độ võng không mong muốn kết cấu, mát dự ứng lực, xuất thêm nội lực thứ cấp biến dạng bổ xung, vv… 2 Các mô hình tính tốn từ biến phổ biến Tổng biến dạng cấu kiện bê tông εc tác dụng σc tổng biến dạng đàn hồi εc, el (t, t0)và biến dạng từ biến εcc(t, t0) ε c (t , t0 ) = ε c ,el (t , t0 ) + ε cc (t , t0 ) Biến dạng đàn hồi xác định theo công thức quen thuộc σ ε c ,el (t , t0 ) = c (Ec mô-đun đàn hồi) Ec Khi σc ≤ 45fck (với fck cường độ chịu nén hình lăng trụ ) lúc biến dạng từ biến tính theo biến dạng đàn hồi theo công thức sau: ε cc (t , t0 ) = ϕ (t , t0 ).ε c ,el (t , t0 ) Với : t tuổi bê tông thời điểm xét t0 tuổi bê tông có ứng suất σc (tuổi chất tải) φ(t, t0) hàm từ biến Khi t hàm từ biến cấu kiện bê tông tăng : - Độ ẩm khơng khí giảm; - Độ mịn xi măng tăng; - Tỷ lệ nước /xi măng tăng; - Độ đặc bê tơng giảm; - Kích thước mặt cắt ngang giảm; - Tuổi chất tải bê tông giảm Tùy theo tiêu chuẩn mà hàm φ(t, t0) tính theo cơng thức khác Sau mơ hình tính tốn phổ biến : *Theo tiêu chuẩn AASTO tiêu chuẩn 22TCN-272-05 [9] Tiêu chuẩn AASHTO áp dụng công thức thực nghiệm Collins Mitchell kiến nghị : H ⎞ −0.118 ⎡ (t − to )0.6 ⎤ ⎛ ϕ (t , t ) = 3.5 × kc × k f × ⎜1.58 − ⎢ ⎟ × t0 0.6 ⎥ 120 ⎠ ⎝ ⎣⎢10 + ( t − t0 ) ⎦⎥ Với : H độ ẩm tương đối môi trường(%) kc hệ số điều chỉnh từ biến phản ánh ảnh hưởng quan hệ thể tích /diện tích Giá trị kc thường xác định theo biểu đồ thực nghiệm hay theo công thức sau : kc = Ac / u (Ac diện tích mặt cắt, u chu vi mặt cắt phần tiếp xúc với mơi trường) Diện tích bề mặt dung để xác định tỷ lệ thể tích bề mặt tính diện tích tiếp xúc với khí Đối với mật hộp kín mà khả thơng gió tính 50% diện tích bề mặt bên hộp kf = 62 42 + f ck f ck cường độ chịu nén hình lăng trụ tuổi 28 ngày bê tông) *Theo tiêu chuẩn CEP- FIB CODE [8] K1=smoframelem(A, L, I, alpha); K=K1*Ect0(ie); KOS=smosoframelem(KOS, K, ix); end [KOS, f]=proores(KOS, f, ixres, vodof); disp=KOS\f; for j=1:nodofos B(j, i)=disp(j) end f=zeros(nodofos, 1); KOS=zeros(nodofos, nodofos); ix=zeros(nonpe*nodofpn, 1); K=zeros(nonpe*nodofpn, nonpe*nodofpn); end B(:, 1)=[]; B(:, noe/2)=[]; B(:, non-2)=[]; B(1, :)=[]; B(1, :)=[]; B(1, :)=[]; B(3*noe/2-2, :)=[]; B(3*noe/2-2, :)=[]; B(3*noe/2-2, :)=[]; B((non-2)*3-2, :)=[]; B((non-2)*3-2, :)=[]; B((non-2)*3-2, :)=[]; deltac=zeros(non-3, non-3); for i=1:noe/2-1 alpha(i)=beta(i+1); end for i=noe/2:non-3 alpha(i)=beta(i+2); end for j=1:non-3 for i=1:non-3 deltac(i, j)=B(3*i-1, j)*sin(alpha(i))+B(3*i-2, j)*cos(alpha(i)) 77 end end for j=1:non-3 for i=1:non-3 deltac0(i, j)=B(3*i-1, j) end end for i=1:noe/2-1 beta(i+1)=alpha(i) end for i=noe/2:non-3 beta(i+2)=alpha(i) end f=zeros(nodofos, 1); KOS=zeros(nodofos, nodofos); ix=zeros(nonpe*nodofpn, 1); K=zeros(nonpe*nodofpn, nonpe*nodofpn); for i=1:non if i==1 f(3*i-2)=0; f(3*i-1)=0; elseif i==noe/2+1 f(3*i-2)=0; f(3*i-1)=0; elseif i==non f(3*i-2)=0; f(3*i-1)=0; else f(3*i-2)=cos(beta(i)); f(3*i-1)=sin(beta(i)); end for ie=1:noe endoe(1)=elem(ie, 1); endoe(2)=elem(ie, 2); x1=coord(endoe(1), 1); y1=coord(endoe(1), 2); 78 x2=coord(endoe(2), 1); y2=coord(endoe(2), 2); L=sqrt((x2-x1)^2+(y2-y1)^2); if (x2-x1)==0 alpha=pi/2; else alpha=atan((y2-y1)/(x2-x1)); end ix=indexos(endoe, nonpe, nodofpn); K1=smoframelem(A, L, I, alpha); K=K1*Ect1(ie); KOS=smosoframelem(KOS, K, ix); end [KOS, f]=proores(KOS, f, ixres, vodof); disp1=KOS\f for j=1:nodofos C(j, i)=disp1(j) end f=zeros(nodofos, 1); KOS=zeros(nodofos, nodofos); ix=zeros(nonpe*nodofpn, 1); K=zeros(nonpe*nodofpn, nonpe*nodofpn); end C(:, 1)=[]; C(:, noe/2)=[]; C(:, non-2)=[]; C(1, :)=[]; C(1, :)=[]; C(1, :)=[]; C(3*noe/2-2, :)=[]; C(3*noe/2-2, :)=[]; C(3*noe/2-2, :)=[]; C((non-2)*3-2, :)=[]; C((non-2)*3-2, :)=[]; C((non-2)*3-2, :)=[]; deltac1=zeros(non-3, non-3); 79 for i=1:noe/2-1 alpha(i)=beta(i+1); end for i=noe/2:non-3 alpha(i)=beta(i+2); end for j=1:non-3 for i=1:non-3 deltac1(i, j)=C(3*i-1, j)*sin(alpha(i))+C(3*i-2, j)*cos(alpha(i)) end end for i=1:noe/2-1 beta(i+1)=alpha(i) end for i=noe/2:non-3 beta(i+2)=alpha(i) end f=zeros(nodofos, 1); KOS=zeros(nodofos, nodofos); ix=zeros(nonpe*nodofpn, 1); K=zeros(nonpe*nodofpn, nonpe*nodofpn); for i=1:non if i==1 f(3*i-2)=0; f(3*i-1)=0; elseif i==noe/2+1 f(3*i-2)=0; f(3*i-1)=0; elseif i==non f(3*i-2)=0; f(3*i-1)=0; else f(3*i-2)=cos(beta(i)); f(3*i-1)=sin(beta(i)); end for ie=1:noe 80 endoe(1)=elem(ie, 1); endoe(2)=elem(ie, 2); x1=coord(endoe(1), 1); y1=coord(endoe(1), 2); x2=coord(endoe(2), 1); y2=coord(endoe(2), 2); L=sqrt((x2-x1)^2+(y2-y1)^2); if (x2-x1)==0 alpha=pi/2; else alpha=atan((y2-y1)/(x2-x1)); end ix=indexos(endoe, nonpe, nodofpn); K1=smoframelem(A, L, I, alpha); K=K1*Ect(ie); KOS=smosoframelem(KOS, K, ix); end [KOS, f]=proores(KOS, f, ixres, vodof); disp2=KOS\f for j=1:nodofos D(j, i)=disp2(j) end f=zeros(nodofos, 1); KOS=zeros(nodofos, nodofos); ix=zeros(nonpe*nodofpn, 1); K=zeros(nonpe*nodofpn, nonpe*nodofpn); end D(:, 1)=[]; D(:, noe/2)=[]; D(:, non-2)=[]; D(1, :)=[]; D(1, :)=[]; D(1, :)=[]; D(3*noe/2-2, :)=[]; D(3*noe/2-2, :)=[]; D(3*noe/2-2, :)=[]; 81 D((non-2)*3-2, :)=[]; D((non-2)*3-2, :)=[]; D((non-2)*3-2, :)=[]; deltac2=zeros(non-3, non-3); for i=1:noe/2-1 alpha(i)=beta(i+1); end for i=noe/2:non-3 alpha(i)=beta(i+2); end for j=1:non-3 for i=1:non-3 deltac2(i, j)=D(3*i-1, j)*sin(alpha(i))+D(3*i-2, j)*cos(alpha(i)) end end nores1=3; ixres1(1)=1; vodof(1)=0; ixres1(2)=2; vodof(2)=0; ixres1(3)=(noe+1)*3-1; vodof((noe+1)*3-1)=0; f=zeros(nodofos, 1); KOS=zeros(nodofos, nodofos); ix=zeros(nonpe*nodofpn, 1); K=zeros(nonpe*nodofpn, nonpe*nodofpn); gamabt=input('nhap luong rieng'); q=gamabt*A for i=1:non if i==1 f(3*i-1)=-q*l(i)/2 f(3*i)=-q*l(i)^2/12 elseif i==non f(3*i-1)=-q*l(i-1)/2 f(3*i)=q*l(i-1)^2/12 else 82 f(3*i-1)=-(q*l(i-1)/2+q*l(i)/2) f(3*i)=q*l(i-1)^2/12-q*l(i)^2/12 end end for ie=1:noe endoe(1)=elem(ie, 1); endoe(2)=elem(ie, 2); x1=coord(endoe(1), 1); y1=coord(endoe(1), 2); x2=coord(endoe(2), 1); y2=coord(endoe(2), 2); L=sqrt((x2-x1)^2+(y2-y1)^2); if (x2-x1)==0 alpha=pi/2; else alpha=atan((y2-y1)/(x2-x1)); end ix=indexos(endoe, nonpe, nodofpn); K1=smoframelem(A, L, I, alpha); K=K1*Ect0(ie); KOS=smosoframelem(KOS, K, ix); end [KOS, f]=proores1(KOS, f, ixres1, vodof); disp3=KOS\f f=zeros(nodofos, 1); KOS=zeros(nodofos, nodofos); ix=zeros(nonpe*nodofpn, 1); K=zeros(nonpe*nodofpn, nonpe*nodofpn); for i=1:non if i==1 f(3*i-1)=-q*l(i)/2; f(3*i)=-q*l(i)^2/12; elseif i==non f(3*i-1)=-q*l(i-1)/2; f(3*i)=q*l(i-1)^2/12; else 83 f(3*i-1)=-(q*l(i-1)/2+q*l(i)/2); f(3*i)=q*l(i-1)^2/12-q*l(i)^2/12; end end for ie=1:noe endoe(1)=elem(ie, 1); endoe(2)=elem(ie, 2); x1=coord(endoe(1), 1); y1=coord(endoe(1), 2); x2=coord(endoe(2), 1); y2=coord(endoe(2), 2); L=sqrt((x2-x1)^2+(y2-y1)^2); if (x2-x1)==0 alpha=pi/2; else alpha=atan((y2-y1)/(x2-x1)); end ix=indexos(endoe, nonpe, nodofpn); K1=smoframelem(A, L, I, alpha); K=K1*Ect1(ie); KOS=smosoframelem(KOS, K, ix); end [KOS, f]=proores1(KOS, f, ixres1, vodof); disp4=KOS\f f=zeros(nodofos, 1); KOS=zeros(nodofos, nodofos); ix=zeros(nonpe*nodofpn, 1); K=zeros(nonpe*nodofpn, nonpe*nodofpn); f((noe/2+1)*3-1)=1; for ie=1:noe endoe(1)=elem(ie, 1); endoe(2)=elem(ie, 2); x1=coord(endoe(1), 1); y1=coord(endoe(1), 2); x2=coord(endoe(2), 1); y2=coord(endoe(2), 2); 84 L=sqrt((x2-x1)^2+(y2-y1)^2); if (x2-x1)==0 alpha=pi/2; else alpha=atan((y2-y1)/(x2-x1)); end ix=indexos(endoe, nonpe, nodofpn); K1=smoframelem(A, L, I, alpha); K=K1*Ect0(ie); KOS=smosoframelem(KOS, K, ix); end [KOS, f]=proores1(KOS, f, ixres1, vodof); disp5=KOS\f Ht0=-disp3((noe/2+1)*3-1)/disp5((noe/2+1)*3-1); KOS=zeros(nodofos, nodofos); ix=zeros(nonpe*nodofpn, 1); K=zeros(nonpe*nodofpn, nonpe*nodofpn); for ie=1:noe endoe(1)=elem(ie, 1); endoe(2)=elem(ie, 2); x1=coord(endoe(1), 1); y1=coord(endoe(1), 2); x2=coord(endoe(2), 1); y2=coord(endoe(2), 2); L=sqrt((x2-x1)^2+(y2-y1)^2); if (x2-x1)==0 alpha=pi/2; else alpha=atan((y2-y1)/(x2-x1)); end ix=indexos(endoe, nonpe, nodofpn); K1=smoframelem(A, L, I, alpha); K=K1*Ect1(ie); KOS=smosoframelem(KOS, K, ix); end [KOS, f]=proores1(KOS, f, ixres1, vodof); 85 disp6=KOS\f KOS=zeros(nodofos, nodofos); ix=zeros(nonpe*nodofpn, 1); K=zeros(nonpe*nodofpn, nonpe*nodofpn); for ie=1:noe endoe(1)=elem(ie, 1); endoe(2)=elem(ie, 2); x1=coord(endoe(1), 1); y1=coord(endoe(1), 2); x2=coord(endoe(2), 1); y2=coord(endoe(2), 2); L=sqrt((x2-x1)^2+(y2-y1)^2); if (x2-x1)==0 alpha=pi/2; else alpha=atan((y2-y1)/(x2-x1)); end ix=indexos(endoe, nonpe, nodofpn); K1=smoframelem(A, L, I, alpha); K=K1*Ect(ie); KOS=smosoframelem(KOS, K, ix); end [KOS, f]=proores1(KOS, f, ixres1, vodof); disp7=KOS\f Ht=(-disp3((noe/2+1)*3-1)-disp4((noe/2+1)*3-1)+(disp7((noe/2+1)*3-1)disp6((noe/2+1)*3-1))*Ht0)/(disp5((noe/2+1)*3-1)+disp7((noe/2+1)*3-1)) for i=1:noe/2-1 alpha(i)=beta(i+1); end for i=noe/2:non-3 alpha(i)=beta(i+2); end disp3(1, :)=[]; disp3(1, :)=[]; disp3(1, :)=[]; disp3(3*noe/2-2, :)=[]; 86 disp3(3*noe/2-2, :)=[]; disp3(3*noe/2-2, :)=[]; disp3((non-2)*3-2, :)=[]; disp3((non-2)*3-2, :)=[]; disp3((non-2)*3-2, :)=[] deltab1=zeros(non-3, 1); for i=1:non-3 deltab1(i, 1)=disp3(3*i-1, 1)*sin(alpha(i)) end disp4(1, :)=[]; disp4(1, :)=[]; disp4(1, :)=[]; disp4(3*noe/2-2, :)=[]; disp4(3*noe/2-2, :)=[]; disp4(3*noe/2-2, :)=[]; disp4((non-2)*3-2, :)=[]; disp4((non-2)*3-2, :)=[]; disp4((non-2)*3-2, :)=[] deltab2=zeros(non-3, 1); for i=1:non-3 deltab2(i, 1)=disp4(3*i-1, 1)*sin(alpha(i)) end disp5(1, :)=[]; disp5(1, :)=[]; disp5(1, :)=[]; disp5(3*noe/2-2, :)=[]; disp5(3*noe/2-2, :)=[]; disp5(3*noe/2-2, :)=[]; disp5((non-2)*3-2, :)=[]; disp5((non-2)*3-2, :)=[]; disp5((non-2)*3-2, :)=[] deltab3=zeros(non-3, 1); for i=1:non-3 deltab3(i, 1)=disp5(3*i-1, 1)*sin(alpha(i)) end disp6(1, :)=[]; 87 disp6(1, :)=[]; disp6(1, :)=[]; disp6(3*noe/2-2, :)=[]; disp6(3*noe/2-2, :)=[]; disp6(3*noe/2-2, :)=[]; disp6((non-2)*3-2, :)=[]; disp6((non-2)*3-2, :)=[]; disp6((non-2)*3-2, :)=[] deltab4=zeros(non-3, 1); for i=1:non-3 deltab4(i, 1)=disp6(3*i-1, 1)*sin(alpha(i)) end disp7(1, :)=[]; disp7(1, :)=[]; disp7(1, :)=[]; disp7(3*noe/2-2, :)=[]; disp7(3*noe/2-2, :)=[]; disp7(3*noe/2-2, :)=[]; disp7((non-2)*3-2, :)=[]; disp7((non-2)*3-2, :)=[]; disp7((non-2)*3-2, :)=[] deltab5=zeros(non-3, 1); for i=1:non-3 deltab5(i, 1)=disp7(3*i-1, 1)*sin(alpha(i)) end disp0=-disp3-Ht0*disp5; disp00=zeros(non-3, 1) for i=1:non-3 disp00(i, 1)=disp0(3*i-1, 1) end Es=input('nhap modun dan hoi cua cap=') for i=1:non-3 As(i)=input('nhap dien tich cap=') end for i=2:noe/2 Lb(1)=l(1); 88 Lb(i)=Lb(i-1)+l(i) end for i=1:noe/2-1 Lc(i)=(Lb(noe/2)-Lb(i))/cos(alpha(i)) end for i=noe/2:non-3 Lc(i)=Lc(noe-i-1) end deltas=zeros(non-3, non-3) for i=1:non-3 deltas(i, i)=Lc(i)/(Es*As(i)) end Xt0=zeros(non-3, 1) k=input('nhap su lua chon') if k==1 for i=1:non-3 Xt0(i, 1)=input('nhap luc cang ban dau') end elseif k==2 Xt0=deltac0\disp00 end cvbd=(deltac+deltas)*Xt0+deltab1+deltab3*Ht0 Xtt=(deltac+deltac2+deltas)\(cvbd-deltab1-deltab2-deltab4*Ht0+deltab5*Ht0deltab3*Ht-deltab5*Ht-deltac1*Xt0+deltac2*Xt0) gamaca=input('nhap luong rieng cua cap') M=deltac+deltac2+deltas N=zeros(non-3, 1) for i=1:non-3 N(i, 1)=gamaca^2*Lc(i)^3*cos(alpha(i))^2*As(i)^2/12; end N deltasp=zeros(non-3, 1) for i=1:non-3 deltasp(i, 1)=(gamaca^2*Lc(i)^3*cos(alpha(i))^2*As(i)^2)/(12*Xt0(i)^2); end 89 Q=deltab2+deltab4*Ht0-deltab5*Ht0+deltab3*Ht+deltab5*Ht+deltac1*Xt0deltac2*Xt0-deltac*Xt0-deltas*Xt0-deltab3*Ht0-deltasp end e=ones(3, 1) n=1 X0=zeros(non-3, 1) for i=1:non-3 X0(i, 1)=Xt0(i, 1); end while e>0 000001 for i=1:non-3 sum=0 for j=1:non-3 sum=sum+M(i, j)*X0(j, 1) K(i, j)=M(i, j) if i==j K(i, j)=M(i, j)-2*N(i)*X0(i, 1)^(-3) end end f(i)=sum+N(i)*X0(i, 1)^(-2)+Q(i) U(i, 1)=f(i) end delta=-K\U; Xpt=X0+delta; X0=Xpt for i=1:3 e(i)=abs(f(i)); end n=n+1 if n>100 break end end 90 TÓM TẮT LÝ LỊCH KHOA HỌC Lý lịch Họ tên : Cao Bá Doãn Ngày sinh : 27/4/1985 Giới tính : Nam Nơi sinh : Quảng Ngãi Địa liên lạc - 133/66/10 Ni Sư Huỳnh Liên P 10 Q Tân Bình Tp Hồ Chí Minh -ĐT : 01686863817 Q trình học tập: Từ năm 2003 – 2008 :Học Đại học trường Đại học Bách Khoa Tp Hồ Chí Minh chuyên ngành Xây Dựng Cầu Đường Từ năm 2008- : Học Cao học trường Đại học Bách Khoa Tp Hồ Chí Minh chuyên ngành Xây Dựng Cầu Hầm 91 ... chuyển vị nội lực kết cấu từ biến Việc thực chương 24 CHƯƠNG PHÂN TÍCH ẢNH HƯỞNG CỦA TỪ BIẾN ĐẾN CẦU TREO DÂY VĂNG Như ta biết, cầu dây văng hệ siêu tĩnh bậc cao ảnh hưởng từ biến thông số quan trọng... cứu từ biến cầu dây văng Như ta biết ảnh hưởng từ biến đến nội lực cầu dây văng yếu tố cần phải xem xét q trình thiết kế Do nghiên cứu từ biến cầu dây văng tiến hành nhằm ước lượng cách xác phân. .. dây văng) ảnh hưởng từ biến làm tăng giảm lực căng dây +Làm tăng lực căng dây ảnh hưởng từ biến trọng lượng thân lớn lực căng dây +Ngược lại làm giảm lực căng dây ảnh hưởng từ biến lực căng dây