Phân tích và chẩn đoán dầm đàn hồi có nhiều vết nứt

ĐẠ I HỌ C QUỐ C GIA THÀNH PHỐHỒCHÍMINH TRƯỜNG ĐẠI HỌC BÁCH KHOA  LÊ XUÂN HÀNG PHÂN TÍCH VÀ CHẨN ĐOÁN DẦM ĐÀN HỒI CÓ NHIỀU VẾT NỨT CHUYÊ N NGÀNH: XÂ Y DƯ Ï NG DÂ N DỤ NG VÀCÔ NG NGHIỆP MÃSỐNGÀNH : 60.58.20 LUẬN VĂN THẠC SĨ TP HỒCHÍMINH, THÁNG NĂ M 2007 CÔ NG TRÌNH ĐƯ Ơ Ï C HOÀN THÀNH TẠ I TRƯ Ơ ØNG ĐẠ I HỌ C BÁCH KHOA TP HỒCHÍ MINH ĐẠ I HỌ C QUỐ C GIA TP HỒCHÍ MINH Cán bộhư ớng dẫn khoa học: PGS TS Nguyễn Thị Hiền Lương Cán bộchấm nhận xét 1: … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … Cán bộchấm nhận xét 2: … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … Luận văn thạc sóđư ợc bảo vệtại: HỘI ĐỒ NG CHẤ M BẢO VỆLUẬN VĂ N THẠ C SĨ TRƯ Ơ ØNG ĐẠ I HỌ C BÁCH KHOA, ngà y … … … tháng … … … năm 2007 TRƯ Ơ ØNG ĐẠ I HỌ C BÁCH KHOA PHÒNG ĐÀO TẠO SAU ĐẠI HỌC CỘNG HÒA XÃHỘI CHỦNGHĨA VIỆT NAM Độc Lập – Tư ïDo - Hạnh Phúc -o0o - NHIỆM VỤ LUẬN VĂN THẠC SĨ Họvàtên học viên: LÊ XUÂN HÀNG Phái: Ngà y, tháng, năm sinh: 27/09/1980 Nam Nơi sinh: Bình Định Chuyên ngà nh: XÂY DỰNG DD & CN Mãsố: 23.04.10 Khoá: K15 MSHV: 02104544 I- TÊN ĐỀ TÀI: PHÂN TÍCH VÀ CHẨN ĐOÁN DẦM ĐÀN HỒI CÓ NHIỀU VẾT NỨT II- NHIỆM VỤ VÀ NỘI DUNG: Phân tích vàchẩn đoán dầm đà n hồi có nhiều vết nư ùt: - Bà i toán thuận: khảo sát ảnh hư ởng vết nư ùt đến tần sốdao động tư ï nhiên dầm - Bà i toán ngư ợc: xác định vị trí vàchiều sâu vết nư ùt dầm sư û dụng thuật toán di truyền III- NGÀY GIAO NHIỆM VỤ: 05/02/2007 IV- NGÀY HOÀN THÀNH NHIỆM VỤ: 05/07/2007 V- HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN: PGS TS NGUYỄN THỊHIỀ N LƯ Ơ NG CÁN BỘ HƯỚNG DẪN CHỦ NHIỆM NGÀNH BỘ MÔN QUẢN LÝ NGÀNH PGS TS NGUYỄN THỊHIỀ N LƯ ƠNG Nội dung vàđềcư ơng luận văn thạc sóđãđư ợc Hội Đồng Chuyên Ngà nh thông qua Ngày PHÒNG ĐÀO TẠ O SAU ĐẠ I HỌ C tháng năm 2007 KHOA QUẢN LÝNGÀNH LỜI CẢM ƠN Trước tiên, xin có lời cảm ơn chân thành đến giáo viên hướng dẫn, PGS TS Nguyễn Thị Hiền Lương, người tận tình hướng dẫn suốt trình thực luận văn Tôi xin chân thành cảm ơn Thầy, Cô truyền đạt kiến thức quý giá năm qua, kiến thức tảng để hoàn thành luận văn Cảm ơn Trường Đại Học Bách Khoa Tp Hồ Chí Minh, trường nuôi dưỡng đào tạo khóa học vừa qua mà suốt thời gian sinh viên Tôi xin gởi lời cảm ơn đến tác giả có tài liệu mà sử dụng thực luận văn Cuối xin bày tỏ lòng biết ơn đến gia đình, bạn bè đồng nghiệp quan tâm động viên giúp đỡ cho trình thực luận văn Chân thành cảm ơn! MỤC LỤC LỜI CẢM ƠN Chương 1: TỔNG QUAN 1.1 Ýnghóa thờ i sư ïcủa đềtà i 1.2 Lý thuyết dao động dầm có vết nư ùt 1.3 Thuật toán di truyeàn 1.4 Mục tiêu luận văn 1.5 Bốcục luận văn Chương 2: MÔ HÌNH PHẦN TỬ HỮU HẠN CỦA DẦM CÓ VẾT NỨT 2.1 Phư ơng pháp phần tư û hư õu hạn 2.2 Phần tư û dầm phẳng vết nư ùt 2.3 Phần tư û dầm phẳng có vết nư ùt 11 2.3.1 Môhình phần tư û dầm cóvết nư ùt 11 2.3.2 Ma trận độcư ùng phần tư û dầm phẳng có vết nư ùt 14 2.4 Xác định tần sốdao động tư ïnhiên 15 Chương 3: THUẬT TOÁN DI TRUYỀN 18 3.1 Giới thiệu 18 3.2 Khái quát vềthuật toán di truyền 19 3.2.1 Mãhoá 19 3.2.2 Các toán tư û di truyền 19 3.2.3 Đặc điểm thuật toán di truyền 22 3.3 Cô chếthư ïc thuật toán di truyền 23 3.3.1 Mãhoá cách biểu diễn 24 3.3.2 Khởi tạo quần thể 25 3.3.3 Toaùn tư û chọn lọc 25 3.3.4 Toán tư û lai gheùp 26 3.3.5 Toán tư û đột biến 27 3.4 Nguyên lý hoạt động toán di truyền 27 3.5 Thuaät toán di truyền Matlab 32 3.5.1 Dạng quần thể 32 3.5.2 Kích thư ớc quần thể 32 3.5.3 Khởi tạo hà m 32 3.5.4 Quần thểban đầu 32 3.5.5 Mư ùc độđánh giá ban đầu 32 3.5.6 Khoảng chênh lệch ban đầu 33 3.5.7 Thay đổi độthích nghi – Hà m đánh giá 33 3.5.8 Lư ïa chọn cá thể 34 3.5.9 Choïn loïc tinh hoa 34 3.5.10 Xác suất lai 35 3.5.11 Hà m đột biến 35 3.5.12 Haø m lai gheùp 36 3.5.13 Sư ïdi trú 37 Chương 4: ẢNH HƯỞNG CỦA VẾT NỨT ĐẾN TẦN SỐ DAO ĐỘNG TỰ NHIÊN CỦA DẦM 38 4.1 Dầm công xôn 39 4.1.1 Ảnh hư ởng vết nư ùt dầm đến tần sốthư ù 39 4.1.2 Ảnh hư ởng vết nư ùt dầm đến tần sốthư ù hai 45 4.1.3 Ảnh hư ởng vết nư ùt dầm đến tần sốthư ù ba 52 4.1.4 Kieåm tra độtin cậy môhình tính 57 4.2 Dầm tư ïa đơn hai đầu 59 4.2.1 Ảnh hư ởng vết nư ùt dầm đến tần sốthư ù 59 4.2.2 Ảnh hư ởng vết nư ùt dầm đến tần sốthư ù hai 63 4.2.3 Ảnh hư ởng vết nư ùt dầm đến tần sốthư ù ba 67 4.2.4 Kiểm tra độtin cậy môhình tính 70 4.3 Keát luaän 71 Chương 5: XÁC ĐỊNH VẾT NỨT TRONG DẦM SỬ DỤNG THUẬT TOÁN DI TRUYỀN 73 5.1 Thiết lập bà i toán tối ö u 75 5.1.1 Các biến số 75 5.1.2 Hà m mục tiêu 75 5.1.3 Các điều kiện rà ng buộc 76 5.2 Ảnh hư ởng thông sốdi truyền đến kết bà i toán 78 5.2.1 Phân tích ảnh hư ởng sốcá thểtrong quần thể 78 5.2.2 Phân tích ảnh hư ởng sốthếhệ 80 5.2.3 Phân tích ảnh hư ởng khoảng chênh lệch ban đầu 81 5.2.4 Phân tích ảnh hư ởng sốcá thểchọn lọc tinh hoa 84 5.2.5 Phân tích ảnh hư ởng sốthếhệsau lần di trú 85 5.2.6 Phân tích ảnh hư ởng xác suất lai 87 5.2.7 Phân tích ảnh hư ởng hệsốScale 89 5.2.8 Phân tích ảnh hư ởng hệsốShrink 91 5.2.9 Phân tích ảnh hư ởng thông sốphạt mu 93 5.2.10 Phân tích ảnh hư ởng giá trịmumax 95 5.3 Độ xác nghiệm trình thiết lập bộthông sốthích hợp thuật toán di truyền 97 5.4 Bộthông sốthích hợp thuật toán di truyền 97 5.5 Xác định vết nư ùt dầm sư û dụng thuật toán di truyền 99 5.5.1 Dầm công xoân 100 5.5.2 Dầm tư ïa đơn hai đầu 101 5.6 Kết luận 103 Chương 6: KẾT LUẬN VÀ KIẾN NGHỊ 104 6.1 Kết luận 104 6.2 Kiến nghịvàhư ớng phát triển 105 TÀI LIỆU THAM KHẢO 107 PHUÏ LUÏC 111 TOÅNG QUAN CHƯƠNG TỔNG QUAN 1.1 Ý NGHĨA THỜI SỰ CỦA ĐỀ TÀI Khuyết tật công trình làmột vấn đềmàchúng ta phải đối mặt suốt trình tồn công trình Đây làmột vấn đềrất đáng lư u ý cù ng với sư ïphát triển nhanh chóng khoa học kỹthuật Khuyết tật công trình thư ng nguyên nhân sau [31]: - Sai sót trình thiết kế, thi công vàsư û dụng công trình - Ư Ùng dụng thiết kếmới, vật liệu ng kết thu đư ợc không mong muốn Khuyết tật nguyên nhân thư ù có khả khắc phục đư ợc Đó õng lỗi thông thư ng thợkhông nh nghề, vật liệu không thích hợp hay không chuẩn, phân tích ùng suất sai, quản lý … có thểkhắc phục dư ïa o kinh nghiệm hay kỹthuật thích hợp Tuy vậy, việc khắc phục nà y không hềđơn giản Khuyết tật nguyên nhân thư ù hai khó khắc phục nhiều Công nghệ có nhiều u điểm ng có nhiều vấn đề tiềm ẩn màngư i thiết kế không thểdư ïbáo trư ớc đư ợc Nhìn chung, nên áp dụng công nghệ hay vật liệu sau đãphân tích vàthư û nghiệm kỹcà ng Như hạn chếđư ợc khuyết tật ng không thểhạn chếtriệt đểđư ợc số yếu tốquan trọng có thểđãbịbỏ qua trình phân tích vàthư û nghiệm Do đó, công việc phân tích vàchẩn đoán khuyết tật có vai tròrất quan trọng việc phát khuyết tật công trình đểtư øđó có biện pháp xư û lý thích hợp Chúng ta sẽtiết kiệm nhiều nhân lư ïc, vật lư ïc phát ngăn ngư ø a kịp thờ i hư hỏng có thểgây nguy hại cho công trình Điều nà y Trang TỔNG QUAN có ý nghóa có sở lý thuyết tin cậy kết hợp việc thư û nghiệm thư ïc tế Với đề tà i nà y, tác giả mong muốn đóng góp o việc xây dư ïng sở lý thuyết phân tích vàchẩn đoán khuyết tật kết cấu đểtư øđó có thểmở rộng phạm vi ùng dụng thư ïc tế, ngăn ngư ø a hư hỏng không phát kịp thờ i có thểdẫn đến hư hỏng, sụp đổcông trình 1.2 LÝ THUYẾT DAO ĐỘNG CỦA DẦM CÓ VẾT NỨT Việc sư û dụng dư õliệu dao động đểphát vết nư ùt đãđư ợc nghiên cư ùu rộng rãi nhiều năm vàthu đư ợc õng bư ớc phát triển quan trọng lónh vư ïc nà y Vết nư ùt cấu kiện kết cấu đư ợc môhình hoá với giả thuyết độmềm cục m ảnh hư ởng đến đáp ùng dao động kết cấu Đặc tính nà y thư ng đư ợc dù ng đểxác định sư ïtồn vết nư ùt kết cấu bao gồm vị trí lẫn chiều sâu Chondros vàDimarogonas [29] xem ảnh hư ởng vết nư ùt lên biến dạng dầm làmột khớp dẻo Gouraris vàDimarogonas [11] đư a mô hình phần tư û hư õu hạn màtrong hai hà m dạng khác tiêu biểu cho hai đoạn dầm, tư øđó khảo sát sư ïbất liên tục vềbiến dạng vết nư ùt Cawley Adams [4] đãchỉ sư ïphân bốư ùng suất kết cấu dao động làkhông đồng vàkhác dạng dao động Vì vậy, vết vết nư ùt nà o ảnh hư ởng đến dạng dao động làkhác phụthuộc o vị trí vết nư ùt Chondros vàDimarogonas [5] đãsư û dụng phư ơng pháp lư ợng vàlý thuyết dầm liên tục bị nư ùt đểphân tích dao động ngang dầm bị nư ùt Trong nghiên cư ùu giải tích vềvấn đềnà y, có hai phư ơng pháp đư ợc sư û dụng đểđịnh lư ợng độ mềm cục vết nư ùt Trong phư ơng pháp thư ù nhất, Rizos đồng sư ï[27] xây dư ïng ma trận độ cư ùng tiết diện nư ùt cách xem vết Trang PHUÏ LUÏC for iel=(nel1+nel2+1):(nel1+nel2+nel3) for i=1:(nnel*ndof) index(i)=(iel-1)*(nnel-1)*ndof+i; end [k,m]=uncrack3(le3,E,J,ro,b,h); kk=ghep(kk,k,index);%lap ghep ma tran cung tong the mm=ghep(mm,m,index);%lap ghep ma tran khoi luong tong the end iel=(nel1+nel2+nel3+nel4); for i=1:(nnel*ndof) index(i)=(iel-1)*(nnel-1)*ndof+i; end [k,m]=crack4(d2,b,h,nuy,E,J,le4,ro,F); kk=ghep(kk,k,index);%lap ghep ma tran cung tong the mm=ghep(mm,m,index);%lap ghep ma tran khoi luong tong the for iel=(nel1+nel2+nel3+nel4+1):(nel1+nel2+nel3+nel4+nel5) for i=1:(nnel*ndof) index(i)=(iel-1)*(nnel-1)*ndof+i; end [k,m]=uncrack5(le5,E,J,ro,b,h); kk=ghep(kk,k,index);%lap ghep ma tran cung tong the mm=ghep(mm,m,index);%lap ghep ma tran khoi luong tong the end kk=kk([3:sdof],[3:sdof]);%dieu kien bien mm=mm([3:sdof],[3:sdof]);%dieu kien bien %TAN SO DAO DONG TU NHIEN w2=eig(kk,mm); w=sqrt(w2); %f=w/2/pi; tanso(o,p) = w(1)/w0(1) end end %VE DO THI plot(x2c/L,tanso(1,:),'^-','LineWidth',1,'Color',[0 0.5 0]) hold on; plot(x2c/L,tanso(2,:),'o ','LineWidth',1,'Color',[1 0]) Trang 114 PHUÏ LUÏC hold on; plot(x2c/L,tanso(3,:),'+ ','LineWidth',1,'Color',[0 1]) hold on; legend('a2=0.05h','a2=0.10h','a2=0.15h',-1) grid on; xlabel('Vi tri vet nut (x2/L)','Fontname','Times New Roman','Color','b','Fontsize',13) ylabel('Tan so thu ba (w3/w03)','Fontname','Times New Roman','Color','b','Fontsize',13) hold off; function [k,m]=uncrack1(le1,E,J,w,b,h) k=E*J/le1^3*[12 6*le1 -12 6*le1; 6*le1 4*le1^2 -6*le1 2*le1^2; -12 -6*le1 12 -6*le1; 6*le1 2*le1^2 -6*le1 4*le1^2]; m=w*b*h*le1/420*[156 22*le1 54 -13*le1; 22*le1 4*le1^2 13*le1 -3*le1^2; 54 13*le1 156 -22*le1; -13*le1 -3*le1^2 -22*le1 4*le1^2]; function [k,m] = crack2(d1,b,h,nuy,E,J,le2,ro,F) s = d1/h; S = sin(pi*s); C = cos(pi*s); T = tan(pi*s); se = sec(pi*s); B1 = pi*(1-nuy^2)/(E*b*h^2); B2 = (0.1215854204e-7/(18*b))*(-139124880*S^3 - 6259043520*S*C^2 + 6259043520*log(se+T)*C^2 + 398772120*T^2*C^2 + 6259043520*log(C)*C^2709379280*S^5 - 2039972880*S^3*C^2 + 842641620*S^6 + 365374820*S^4*C^2 +2730749640*S^2*C^2 - 798356160*S^5*C^2 + 33264800*S^8 + 348488800*S^6*C^2 - 95042400*S^9 - 95042400*S^7*C^2 + 11880300*S^10 + 11880300*S^8*C^2 + 25456707 - 25456707*C^2)/C^2; B3 = (1/b)*(-0.01178181818*s^11 + 0.01368*s^10 + 0.1101488889*s^9 - 0.356885*s^8 + 0.082612*s^7+1.483734833*s^6 - 3.1577768*s^5 + 2.66192*s^4 - 0.2274253333*s^3 - 0.341138*s^2-0.682276*s - 0.682276*log(abs(s-1))); C11 = le2^3/(3*E*J) + 2*B1*(9*le2^2*B2 + h^2*B3); C22 = le2/(E*J) + 72*B1*B2; Trang 115 PHUÏ LUÏC C12 = le2^2/(2*E*J) + 36*le2*B1*B2; C21 = C12; k = (1/(C11*C22-C12*C21))*[ C22 C22*le2-C21 -C22 C21; C22*le2-C21 C22*le2^2-C21*le2-C12*le2+C11 -C22*le2+C12 C21*le2-C11; -C22 -C22*le2+C21 C22 -C21; C12 C12*le2-C11 -C12 C11]; m = (ro*F*le2/420)*[156 22*le2 54 -13*le2; 22*le2 4*le2^2 13*le2 -3*le2^2; 54 13*le2 156 -22*le2; -13*le2 -3*le2^2 -22*le2 4*le2^2]; function [kk]=ghep(kk,k,index) edof=4;%edof=nnel*ndof; %edof=length(index); for i=1:edof ii=index(i); for j=1:edof jj=index(j); kk(ii,jj)=kk(ii,jj)+k(i,j); end end Chương trình tính tần số riêng dầm tựa đơn hai đầu clc clear format long; L=10;%m(chieu dai dam) E=2.8e10;%N/m2(modun dan hoi cua vat lieu) ro=2350;%kg/m3(trong luong rieng cua vat lieu thep) nuy=0.3;%he so Poisson b=0.2;%m(chieu rong tiet dien) h=0.6;%m(chieu cao tiet dien) F=b*h; J=b*h^3/12;%m4(momen quan tinh tiet dien dam) nnel=2;%so nut mot phan tu ndof=2;%so bac tu tren nut %BAI TOAN DAM KHONG NUT sel0=10; le0=L/sel0; Trang 116 PHUÏ LUÏC snode0=(nnel-1)*sel0+1;%tong so nut dam sdof0=snode0*ndof;%tong so bac tu cua dam kk0=zeros(sdof0,sdof0);%ma tran cung ban dau mm0=zeros(sdof0,sdof0);%ma tran khoi luong ban dau index0=zeros(sel0*ndof,1);%vecto chi so dau vao for iel=1:sel0 for i=1:(nnel*ndof) index0(i)=(iel-1)*(nnel-1)*ndof+i; end [k0,m0]=uncrack0(le0,E,J,ro,b,h); kk0=ghep(kk0,k0,index0);%lap ghep ma tran cung tong the mm0=ghep(mm0,m0,index0);%lap ghep ma tran khoi luong tong the end kk0([1 (sdof0-1)],:)=[];%dieu kien bien mm0([1 (sdof0-1)],:)=[];%dieu kien bien kk0(:,[1 (sdof0-1)])=[];%dieu kien bien mm0(:,[1 (sdof0-1)])=[];%dieu kien bien %TAN SO DAO DONG TU NHIEN w02=eig(kk0,mm0); w0=sqrt(w02); %f0=w0/2/pi; %BAI TOAN DAM CO NHIEU VET NUT d2c = h*[0.05 0.1 0.15]; x2c = L*[0.02 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.98]; tanso = zeros(length(d2c),length(x2c)); for p =1:length(x2c) for o = 1:length(d2c) x1t=0.35*L; d1t=0.1*h; d2t = d2c(o); x2t = x2c(p); %CHIA PHAN TU if x2t>=x1t x1=x1t; d1=d1t; Trang 117 PHUÏ LUÏC x2=x2t; d2=d2t; else x1=x2t; d1=d2t; x2=x1t; d2=d1t; end %CHIA PHAN TU l2=0.01*L;%KICH THUOT PHAN TU l4=l2;% l1=x1-l2/2; l2; l3=(x2-x1)-(l2+l4)/2; l4; l5=L-(l1+l2+l3+l4); nel1=max(round(l1/l2),1); le1=l1/nel1; le2=l2; nel2=1; nel3=max(round(l3/l2),1); le3=l3/nel3; le4=l4; nel4=1; nel5=max(round(l5/l2),1); le5=l5/nel5; %MA TRAN DO CUNG & KHOI LUONG sel=nel1+nel2+nel3+nel4+nel5; snode=(nnel-1)*sel+1;%tong so nut dam sdof=snode*ndof;%tong so bac tu cua dam kk=zeros(sdof,sdof);%ma tran cung ban dau mm=zeros(sdof,sdof);%ma tran khoi luong ban dau index=zeros(sel*ndof,1);%vecto chi so dau vao Trang 118 PHUÏ LUÏC for iel=1:nel1 for i=1:(nnel*ndof) index(i)=(iel-1)*(nnel-1)*ndof+i; end [k,m]=uncrack1(le1,E,J,ro,b,h); kk=ghep(kk,k,index);%lap ghep ma tran cung tong the mm=ghep(mm,m,index);%lap ghep ma tran khoi luong tong the end iel=(nel1+nel2); for i=1:(nnel*ndof) index(i)=(iel-1)*(nnel-1)*ndof+i; end [k,m]=crack2(d1,b,h,nuy,E,J,le2,ro,F); kk=ghep(kk,k,index);%lap ghep ma tran cung tong the mm=ghep(mm,m,index);%lap ghep ma tran khoi luong tong the for iel=(nel1+nel2+1):(nel1+nel2+nel3) for i=1:(nnel*ndof) index(i)=(iel-1)*(nnel-1)*ndof+i; end [k,m]=uncrack3(le3,E,J,ro,b,h); kk=ghep(kk,k,index);%lap ghep ma tran cung tong the mm=ghep(mm,m,index);%lap ghep ma tran khoi luong tong the end iel=(nel1+nel2+nel3+nel4); for i=1:(nnel*ndof) index(i)=(iel-1)*(nnel-1)*ndof+i; end [k,m]=crack4(d2,b,h,nuy,E,J,le4,ro,F); kk=ghep(kk,k,index);%lap ghep ma tran cung tong the mm=ghep(mm,m,index);%lap ghep ma tran khoi luong tong the for iel=(nel1+nel2+nel3+nel4+1):(nel1+nel2+nel3+nel4+nel5) for i=1:(nnel*ndof) index(i)=(iel-1)*(nnel-1)*ndof+i; end [k,m]=uncrack5(le5,E,J,ro,b,h); kk=ghep(kk,k,index);%lap ghep ma tran cung tong the mm=ghep(mm,m,index);%lap ghep ma tran khoi luong tong the end kk([1 (sdof-1)],:)=[];%dieu kien bien Trang 119 PHUÏ LUÏC mm([1 (sdof-1)],:)=[];%dieu kien bien kk(:,[1 (sdof-1)])=[];%dieu kien bien mm(:,[1 (sdof-1)])=[];%dieu kien bien %TAN SO DAO DONG TU NHIEN w2=eig(kk,mm); w=sqrt(w2); %f=w/2/pi; tanso(o,p) = w(3)/w0(3) end end %VE DO THI plot(x2c/L,tanso(1,:),'^-','LineWidth',1,'Color',[0 0.5 0]) hold on; plot(x2c/L,tanso(2,:),'o ','LineWidth',1,'Color',[1 0]) hold on; plot(x2c/L,tanso(3,:),'+ ','LineWidth',1,'Color',[0 1]) hold on; legend('a2=0.05h','a2=0.10h','a2=0.15h',-1) grid on; xlabel('Vi tri vet nut (x2/L)','Fontname','Times New Roman','Color','b','Fontsize',13) ylabel('Tan so thu ba (w3/w03)','Fontname','Times New Roman','Color','b','Fontsize',13) hold off; BÀI TOÁN NGƯC Chương trình xác định vết nứt dầm công xôn clear; clc; format long btthuandamcongxon for i=1:6 fre(i)=f(i); end %DU LIEU DAU VAO Nvar=4;%so bien so Nind=30;%so ca the quan the Trang 120 PHUÏ LUÏC Maxgen=100;%so the he tinh toan %Matrix=randn(Nind,Nvar);%quan the dau tien ngau nhien Matrix=repmat([vtri1 0.0075 vtri2 0.0075],Nind,1);%quan the dau gan vung kha thi mientri=[0 0 0;0.1 0.01 0.1 0.01]; %mientri=[0;0.1]; Scale=1;%xac xuat dot bien ban dau Shrink=1;%xac xuat dot bien pc=0.8;%xac xuat lai nElite=3;%so ca the chon loc tinh hoa nInterval=30;%sau interval thi cac ca the se di tru %CAC THAM SO PHAT (xu ly cac rang buoc) mu=10;%gan truoc gia tri cho tham so phat mu mumax=1e9;%gia tri mu lon nhat cho phep pStep=100;%hang so nang cao muc phat sau moi buoc lap pTol=1e-10;%gia tri danh gia su hoi tu cua phuong phap phat %CAC THONG SO TRONG QUA TRINH GIAI options=gaoptimset('Display','iter', 'Generation',Maxgen, %so the he tinh toan 'InitialPopulation',Matrix, %quan the khoi tao ban dau 'PlotFcns',{@gaplotrange,@gaplotbestf}, %bieu dien ket qua duoi dang thi 'EliteCount',nElite, %so ca the chon loc tinh hoa 'PopInitRange',mientri, %khoang chenh lech ban dau 'MigrationInterval',nInterval, %sau bao nhieu the he thi di tru 'PopulationSize',Nind, %kich thuoc quan the 'SelectionFcn',@selectionroulette, %ham lua chon quan the moi 'CrossoverFcn',@crossoverscattered, %ham lai 'MutationFcn',{@mutationgaussian,Scale,Shrink}, %ham dot bien 'CrossoverFraction',pc, %xac xuat lai 'TimeLimit',Inf, %khong che thoi gian tinh toan 'StallGenLimit',Inf);%khong che so the he tinh toan %QUA TRINH LAP ctr=0;%bien dem so buoc lap flag=true;%dieu kien hoi tu while flag ctr=ctr+1; [x,cFval,e,o,p,s]=ga({@fitness_congxon,mu,fre},Nvar,options);%danh gia the he hien tai bang thuat toan di truyen cons=constrain(x);%danh gia cac rang buoc tai gia tri x vua tim Trang 121 PHUÏ LUÏC ncon=0; for i=1:length(cons) if cons(i)>0 ncon=ncon+1;%dem so rang buoc bi vi pham end end ncon=ncon; fval(ctr)=cFval;%gia tri cua ham muc tieu if cons(i)>0 fval(ctr)=fval(ctr)-mu*sum(cons(i).^2);%gia tri ham muc tieu nhung cong them mot luong phat else fval(ctr)=fval(ctr); end if (ctr>1 && abs(fval(ctr)-fval(ctr-1))