Ứng dụng thuật toán di truyền để tối ưu hóa cấu kiện dầm bê tông cốt thép trong xây dựng cầu,luận án thạc sĩ khoa học kỹ thuật chuyên ngành xây dựng cầu hầm

Bộ giáo dục đào tạo Trường đại học Giao thông Vận tải áp dụng thuật toán di truyền để tối ưu hoá cấu kiện dầm BTCT xây dựng cầu Luận văn thạc sỹ kỹ thuật Giáo viên hướng dẫn : GS.TS Nguyễn Viết Trung Học viên: Ưng Thị Việt Lan Hà nội, tháng 01 năm 2006 Bộ giáo dục đào tạo Trường đại học Giao thông Vận tải áp dụng thuật toán di truyền để tối ưu hoá cấu kiện dầm BTCT xây dựng cầu Luận văn thạc sỹ kỹ thuật Chuyên ngành : Xây dựng cầu hầm Mà số : 60.58.25 Giáo viên hướng dẫn : GS.TS Nguyễn Viết Trung Học viên: Ưng Thị Việt Lan Hà nội, tháng 01 năm 2006 Mục lục Chương Bài toán tối ưu Nội dung thuật toán di truyền 1.1 Khái niệm mô hình toán tối ưu hoá 1.2 Phân loại phương pháp giải toán thiết kế tối ưu 1.3 Khái niệm hệ thống thiết kế tối ưu, đặc điểm ràng buộc Hệ 11 thống chương trình Hiệu tht to¸n [1] 1.3.1 Kh¸i niƯm chung 11 1.3.2 ChÕ độ đối thoại kỹ sư máy tính để tìm nghiệm tối ưu 14 1.3.3 Phân tích hiệu thuật toán tối ưu 15 1.4 Khái niệm đặc điểm chđ u cđa tht to¸n di trun (TTDT) 16 1.4.1 Khái niệm TTDT 16 1.4.1.1 Nguyên lý tiến hoá tự nhiên 16 1.4.1.2 Thuật toán di truyền 17 1.4.2 Đặc ®iĨm chđ u cđa TTDT (GA) 18 1.5 Tr×nh tù ¸p dơng TTDT 19 1.6 So s¸nh TTDT víi c¸c thuật toán tìm kiếm khác 21 1.7 Nội dung TTDT 22 1.8 Cơ chế thực Thuật giải di truyền 25 1.9 Nguyên lý hoạt động Thuật giải di truyền 29 1.10 Xử lý ràng buộc 36 Chương :áp dụng thuật toán di truyền (Ga) tính toán tối ưu mặt cắt dầm bê tông cốt thép 2.1 Thiết kế cấu kiện dầm bê tông cốt thép xây dựng cầu 41 42 (22 TCN-272-05) 2.1.1 Giới thiệu 42 2.1.2 Những yêu cầu tính toán 42 2.1.2.1 Trạng thái giới hạn sử dụng 43 2.1.2.2 Trạng thái giới hạn cường độ 43 2.1.2.3 Tải trọng hệ số tải trọng 43 2.1.2.4 Hoạt tải 44 2.1.2.5 Hệ số phân bố ngang cho loại cầu Dầm - Bản 47 2.2 Nội dung kỹ thuật toán 54 2.2.1 Các thông số 55 2.2.2 Biến thiết kế 55 2.2.3 Hàm mục tiêu 56 2.2.4 Các ràng buộc 57 2.2.4.1 Kiểm toán theo trạng thái giới hạn 57 2.2.4.2 Kiểm toán theo trạng thái giới hạn sử dụng 59 2.2.4.3 Ràng buộc kích thước 60 2.3 Mô hình toán thiết kế tối ưu tương ứng 60 2.4 Phương pháp giải 61 2.4.1 Giải theo thuật toán di truyền (GA) 61 2.4.1.1 Phát triển mà nguồn phần mềm 61 2.4.1.2 M· ho¸ biÕn thiÕt kÕ 62 2.4.1.3 Chän lùa c¸c thông số GA 63 2.4.1.4 Kết tối ưu 83 2.4.2 Giải theo phương pháp truyền thống, phổ biến 83 2.4.2.1 Cấu trúc chương trình 84 2.4.2.2 Kết tối ưu 84 2.4.3 So sánh kết tính toán 85 Chương Kết luận kiến nghị 3.1 Kết luận 87 3.2 Kiến nghị 88 Tài liệu tham khảo 90 Chương Bài toán tối ưu Nội dung thuật toán di truyền 1.1 Khái niệm mô hình toán tối ưu hoá Với phát triển kỹ thuật máy tính điện tử, toán tối ưu hoá chế độ làm việc thiết bị; kết cấu công trình; hàm truyền dẫn; hàm điều khiển phi tuyến; thông số điều khiển ngày xem xét hoàn chỉnh, ngày có điều kiện, công cụ giải hữu hiệu, ngày có vai trò, ý nghÜa quan träng thiÕt thùc ®èi víi ®êi sèng x· hội Mục đích toán thiết kế tính toán yếu tố đối tượng, thiết kế cho thoả mÃn yêu cầu quy định Các yếu tố đối tượng thiết kế kích thước hình học, khối lượng, giá thành Một yêu cầu quy định độ bền, độ ổn định Rõ ràng, yếu tố đối tượng thiết kế có giới hạn định mà thông thường người ta gọi giới hạn tiêu chuẩn yếu tố Chính khoảng giới hạn đà thúc nhà thiết kế tìm cách vươn tới giá trị tốt mà yếu tố đối tượng thiết kế đạt Bài toán thiết kế với ý nghĩa, mục đích gọi toán thiết kế tối ưu Yếu tố mà nhà thiết kế định hướng phải đạt đến giá trị tốt gọi Hàm mục tiêu lựa chọn nhà thiết kế Điều này, cách tự nhiên làm xuất khái niệm đơn mục tiêu đa mục tiêu, toán tối ưu hoá đa mục tiêu thường phức tạp khó giải Một yếu tố đối tượng thiết kế phụ thuộc vào nhiều tham số, vÝ dơ u tè h×nh häc phơ thc kÝch th­íc chiều rộng, kích thước chiều dài mặt cắt Các tham số độc lập phụ thuộc lẫn Để xét ảnh hưởng cụ thể tham số (ảnh hưởng độc lập), người ta chọn tham số độc lập để làm biến hàm mục tiêu quan hệ học, toán học Nói chung, quan hệ quan hệ phi tuyến Các tham số nhiều yếu tố đối tượng thiết kế, vai trò biến số nhiều hàm mục tiêu, đồng thời biến số nhiều quan hệ ràng buộc tiêu chuẩn mang ý nghĩa kỹ thuật kinh tế mà người ta thường gọi hàm ràng buộc yếu tố đối tượng thiết kế Nói chung, thông thường quan hệ số độc lập (biến số) với hàm ràng buộc quan hÖ phi tuyÕn Khi ký hiÖu : x j (j= 1, n ) tham số độc lập yếu tố đối tượng thiết kế F hàm mục tiêu yếu tố đối tượng thiết kế R i (i=1, m ) víi m 0) x j phải thoả mÃn quan hệ ràng buộc R i R i gọi điểm cho phép (nghiệm cho phép) Tập hợp x j cho phép gọi miền cho phép G G liên tục rời rạc, miền liên thông miền riêng biệt Điểm X* gọi cực tiểu néi bé (tèi ­u côc bé) nÕu F(X*) ≤ F(x j ) thoả mÃn với điểm nằm G, lân cận X* Điểm X* nằm trªn G (trªn biªn cđa miỊn cho phÐp) Khi ®iĨm X* khiÕn cho F(X*) ≤ F(x j ) ®óng với điểm G điểm X* gọi điểm cực trị toàn miền (tối ưu tuyệt đối) G Tuy nhiên việc tìm điểm X* cực trị toàn miền mục tiêu cao toán thiết kế tối ưu thường không dễ dàng Một công việc coi quan trọng toán thiết kế tối ưu xác lập hàm mục tiêu Trên sở có hàm mục tiêu lần tìm đến tham số liên quan đến yếu tố chọn làm hàm mục tiêu ràng buộc để xây dựng nên mô hình toán học đối tượng thiết kế sau phương hướng giải Tuỳ theo yêu cầu thực tế mà có toán tối ưu đơn mục tiêu (chỉ tối ưu yếu tố đối tượng thiết kế) toán tối ưu đa mục tiêu (tối ưu nhiều yếu tố đối tượng thiết kế) Nói chung việc lập hàm mục tiêu tổng quát toán phức tạp Có nhiều phương pháp khác để giải vấn đề (phương pháp tính cộng mục tiêu; phương pháp xếp hạng mục tiêu theo mức độ quan trọng chúng đưa vào dung sai i cho mục tiêu; phương pháp tự thích nghi để lập hàm mục tiêu ) 1.2 Phân loại phương pháp giải toán thiết kế tối ưu [1] Trên thực tế, toán thiết kế tối ưu dạng tổng quát toán nhiều cực trị, toán thiết kế cụ thể đà cho dù toán cực trị, tính chất phức tạp quan hệ hàm mục tiêu ràng buộc c¸c quan hƯ phi tun víi tham sè thiÕt kÕ, nên, nói chung, coi toán nhiều cực trị, đó, cần tới phương pháp để xác định cực trị toàn miền (cực trị tuyệt đối) Nói chung, đưa phương pháp tổng quát cho phép tìm lời giải xác toán quy hoạch phi tuyến tổng quát sau số hữu hạn bước lặp, đồng thời, việc chøng minh b»ng lý thut vỊ sù héi tơ ®Õn nghiệm tối ưu phương pháp sử dụng tính toán thực tế khó khăn phức tạp Về mặt lý thuyết, toán thiết kế tối ưu tổng quát (bài toán nhiều cực trị), nghiên cứu số trường hợp riêng không nhiều Với phát triển ngày cao kỹ thuật tin học, phương pháp số để giải toán quy hoạch phi tuyến máy tính ngày ứng dụng sâu rộng, đà đưa lại kết đáng kể trình giải toán thiết kế tối ưu tổng quát Có thể biểu thị cách phân loại phương pháp giải toán thiết kế tối ưu mô hình tất định đối tượng thiết kế, tất tham số thay đổi liên tục, theo bảng 1.1 Mỗi mức bảng phân loại biểu thị bảng 1.1 tương ứng với giai đoạn tính toán đó, tương ứng với phương pháp cụ thể sử dụng cho giai đoạn Các đường nét đứt thẳng đứng phân chia giai đoạn tìm kiếm nằm thành phần thủ tục tìm kiếm mức cao Các đường nét liền thẳng đứng phân chia phương án khác Theo bảng 1.1, việc tìm tối ưu toàn miền bao gồm phương pháp tạo điểm ban đầu kết hợp với phương pháp tìm tối ưu cục Mục đích tìm tối ưu toàn miền (tối ưu tuyệt đối) nhìn nhận cách bao quát toàn miền cho phÐp, råi kÕt ln chän ®iĨm tèi ­u toàn miền danh sách điểm tối ưu cục đà tìm Các thuật toán tìm tối ưu toàn miền thường xếp loại tuỳ theo phương pháp tạo điểm ban đầu Mỗi phương pháp mô hình, giai đoạn có ưu nhược điểm khác Việc lựa chọn phương pháp phải bảo đảm phù hợp với mô hình giai đoạn, không phụ thuộc vào khả hiệu mặt lý thuyết phương pháp, mà tuỳ thuộc vào yêu cầu, đối tượng cụ thể toán tối ưu cụ thể //Ki1:=0.075+xmuy(s/2900,0.6)*xmuy(s/l,0.2)*xmuy(Kg/(Lt*ts*ts*ts) ,0.1); Ki2:=0.2+s/3600-sqr(s/10700); //writeln('Kg=',kg); //writeln('Ki1=',ki1:1:4); //writeln('Ki2=',ki2:1:4); Mht:=1635.396*ki1; Qht:=128.017*ki2; Mng:=45.41;Qng:=0; Mlp:=245.79;Qlp:=0; Mtt:=37.845*(5*a*(l/1000)*g_bt+114.66+105.84)/(5*(l/1000)); Qtt:=0; Msd:=Mht+Mng+Mlp+Mtt; Qsd:=Qht+Qng+Qlp+Qtt; //writeln('Msd=',Msd:1:4); //writeln('Qsd=',Qsd:1:4); Mgh1:=Mht*1.75+Mng*1.75+Mlp*1.5+Mtt*1.25; Qgh1:=Qht*1.75+Qng*1.75+Qlp*1.5+Qtt*1.25; //writeln('Mgh1=',Mgh1:1:4); //writeln('Qgh1=',Qgh1:1:4); end; function TGAtktudam.ttgh1:integer; var Mr,mn,aa,c,As1,av,beta1,nx,dx,sx,dv,vn1,vn2,vs,vc,vr:double; teta,alpha,z,dc,Aa1,j,n,k,ro,fs1:double; vp:integer; begin tinhnoiluc; vp:=0; As1:={tongi(So_hang1,thanh_1_hang)*0.25*pi*sqr(d_cot_thep)}16*32*32*pi/4; //writeln('as=',as1:1:3); if fcMgh1) and (Vr>Qgh1) and (c/ds=0.03*fc/fs)and(z/(xmuy(dc*aa1,1/3)) 1e-12 then Result := (1 / ham_muc_tieu)*ham_phat //this minimizes error (find smallest error) else Result := / 1e-12; //prevents divide by zero error} //ShowMessage('v='+IntToStr(ktra_hhoc)+'fit='+FloatToStr(Result)) ; end; function TGAtktudam.Lay_Phuong_an_toi_uu: String; begin Result := FChromosomes[FBestFitnessChromIndex].GetGenesAsStr; end; B chương trình thiết kế tối ưu dầm BTCT theo phương pháp MongteCaclô, ngôn ngữ pascal 103 Chương trình tính const h_min:real=900;h_max:real=2500; hc_min:real=100;hc_max:real=250; bs_min:real=100;bs_max:real=300; hv_min:real=100;hv_max:real=300; bv_min:real=100;bv_max:real=300; hvb_min:real=100;hvb_max:real=200; bvb_min:real=100;bvb_max:real=200; hb_min:real=100;hb_max:real=300; bb_min:real=100;bb_max:real=600; dg_bt:real=70.95; dg_thep:real=5652; dg_vk:real=4.741; type vt=array[0 100] of real; vti=array[0 100] of byte; var h,bc,hdam,hc,bs,hs,hv,bv,bvb,hvb,hb,bb,g_bt,fc,fs,Et,Ec:real; f:text; so_pa:real; bt_baove,d_cot_thep:real; cach_hang:real; So_hang:byte; thanh_1_hang:vti; ds,Ki,s,L,Lt,ts,eg,Kg,ki1,ki2:real; Msd,Qsd,Mgh1,Qgh1,Mht,Qht,Mng,Qng,Mtt,Qtt,Mlp,Qlp:real; A,A_t,A_tmin,A_tmax,ex,ey,Ix,Iy,I_xy,Yb,Yt,Wb,Wt,W,gbt:real; ham_muc_tieu,vk:real; ham_muc_tieu_min:real; h_tu,hc_tu,bs_tu,hs_tu,hv_tu,bv_tu,bvb_tu,hvb_tu,hb_tu,bb_tu:real ; function sonn(dau,cuoi:real):real; var s,n:real;i:integer; begin repeat i:=random(1000); s:=0.015; n:=i*s; until (n>=dau) and (na[j] then min:=a[j]; minds:=min; end; function maxds(n:byte;A:vt):real; var max:real;j:byte; begin max:=a[0]; for j:=1 to n if max=cd)and(x1=bs); end; function tinh_ds:real; var s:vt; i,j,k:byte; begin k:=0; for i:=1 to So_hang 106 begin for j:=1 to thanh_1_hang[i] begin inc(k); if i=1 then s[k]:=bt_baove+d_cot_thep/2 else s[k]:=s[k-1]+cach_hang; end; end; end; procedure dactrung(var A,vk,ex,ey,Ix,Iy,I_xy,Yb,Yt,Wb,Wt,W:real); var x,y:vt;i:byte; ai,sxx,syy,ixx,iyy,ixy:vt; sx,sy,Sixx,Siyy,Sixy:real; begin x[0]:=0; y[0]:=h; x[1]:=0; y[1]:=h-hc; x[2]:=(bc-bs)/2-bv; y[2]:=y[1]; x[3]:=x[2]+bv; y[3]:=y[2]-hv; x[4]:=x[3]; y[4]:=hb+hvb; x[5]:=x[4]-bvb; y[5]:=hb; x[6]:=x[5]; y[6]:=0; x[7]:=x[6]+bb; y[7]:=0; x[8]:=x[7]; y[8]:=y[5]; x[9]:=x[8]-bvb; y[9]:=y[4]; x[10]:=x[9]; y[10]:=y[3]; x[11]:=x[10]+bv; y[11]:=y[2]; x[12]:=bc; y[12]:=y[2]; x[13]:=bc; y[13]:=y[0]; x[14]:=0; y[14]:=y[0]; vk:=0; 107 for i:=0 to 13 vk:=vk+sqrt(sqr(x[i]-x[i+1])+sqr(y[i]-y[i+1])); vk:=vk-bc; vk:=vk/1000; a:=dtdg(14,x,y); a:=a/1000000; for i:=0 to 13 begin ai[i]:=(x[i]*y[i+1]-x[i+1]*y[i])/1000000; sxx[i]:=ai[i]*((y[i]+y[i+1]))/1000; syy[i]:=ai[i]*((x[i]+x[i+1]))/1000; ixx[i]:=ai[i]*(sqr(y[i]+y[i+1])-y[i]*y[i+1])/1000000; iyy[i]:=ai[i]*((sqr(x[i]+x[i+1])-x[i]*x[i+1]))/1000000; ixy[i]:=ai[i]*(((x[i]+x[i+1])*(y[i]+y[i+1])(x[i]*y[i+1]+x[i+1]*y[i])/2))/1000000; end; sx:=tong(14,sxx)/6; sy:=tong(14,syy)/6; Sixx:=tong(14,ixx)/12; Siyy:=tong(14,iyy)/12; Sixy:=tong(14,ixy)/12; ex:=(sy)/a; ey:=(sx)/a; writeln('Sixx=',sixx:1:3); writeln('Siyy=',siyy:1:3); writeln('Sixy=',sixy:1:3); Ix:=Sixx-a*sqr(ey); Iy:=Siyy-a*sqr(ex); I_xy:=Sixy-a*ex*ey; yb:=ey-minds(14,y)/1000; yt:=maxds(14,y)/1000-ey; wb:=ix/yb;wt:=ix/yt; {assign(f,'dam.lsp'); rewrite(f); write(f,'(command "pline" "',x[0],',',y[0],'" "W" "0.0" "'); for i:=1 to 14 write(f,'" "',x[i],',',y[i]); writeln(f,'""")'); close(f);} end; {} function ttgh1:boolean; var Mr,mn,aa,c,As1,av,beta1,nx,dx,sx,dv,vn1,vn2,vs,vc,vr:real; teta,alpha,z,dc,Aa1,j,n,k,ro,fs1:real; begin ttgh1:=false; 108 As1:=a_t*1000000; writeln('as=',as1:1:3); if fcMgh1) and (Vr>Qgh1) and (c/ds=0.03*fc/fs)and(z/(xmuy(dc*aa1,1/3))1000; close(f); 111 {tim bien tk ung voi pa A} assign(f,'ketqua.txt'); reset(f); readln(f,h,hc,bs,hv,bv,bvb,hvb,hb,bb,ham_muc_tieu); ham_muc_tieu_min:=ham_muc_tieu; while not eof(f) begin readln(f,h,hc,bs,hv,bv,bvb,hvb,hb,bb,ham_muc_tieu); if ham_muc_tieu_min>ham_muc_tieu then begin ham_muc_tieu_min:=ham_muc_tieu; h_tu:=h;hc_tu:=hc;bs_tu:=bs;hv_tu:=hv;bv_tu:=bv;bvb_tu:=bvb; hvb_tu:=hvb;hb_tu:=hb;bb_tu:=bb; end; end; close(f); assign(f,'ketqua.txt'); append(f); writeln(f,'Thong so ung voi phuong an toi uu'); writeln(f,h_tu:4:2,' ',hc_tu:4:2,' ',bs_tu:4:2,' ',hv_tu:4:2,' ', bv_tu:4:2,' ',bvb_tu:4:2,' ',hvb_tu:4:2,' ',hb_tu:4:2,' ',bb_tu:4:2,' ',ham_muc_tieu_min:4:3); close(f); {} end KÕt qu¶ tÝnh File kết chạy với số lần lặp 1000, với kết h, hc, bs, hv, bv, hvb, bvb, hb, as, hàm mục tiêu 900.00 169.75 277.00 160.00 163.00 110.50 113.50 280.00 542.50 106.70 1068.00 102.25 283.00 295.00 295.00 151.00 197.50 175.00 355.00 66.16 1092.00 100.00 100.00 265.00 142.00 100.00 134.50 295.00 362.50 68.80 1260.00 210.25 226.00 130.00 265.00 160.00 187.00 253.00 355.00 87.68 972.00 140.50 292.00 280.00 100.00 182.50 118.00 247.00 430.00 63.15 924.00 194.50 139.00 202.00 277.00 151.00 122.50 187.00 580.00 105.13 1140.00 187.75 205.00 274.00 145.00 116.50 154.00 211.00 475.00 119.48 1020.00 205.75 244.00 235.00 172.00 199.00 166.00 202.00 437.50 90.43 1164.00 226.00 295.00 100.00 157.00 152.50 157.00 298.00 497.50 103.19 1068.00 185.50 196.00 121.00 169.00 110.50 145.00 298.00 362.50 96.46 948.00 205.75 175.00 157.00 277.00 161.50 106.00 256.00 490.00 93.51 1092.00 169.75 112.00 262.00 145.00 104.50 148.00 199.00 302.50 83.59 948.00 208.00 241.00 133.00 163.00 161.50 109.00 175.00 280.00 78.21 1284.00 208.00 262.00 262.00 220.00 176.50 173.50 280.00 355.00 81.65 112 948.00 178.75 229.00 100.00 226.00 179.50 112.00 232.00 325.00 67.73 924.00 172.00 100.00 208.00 262.00 164.50 187.00 271.00 557.50 102.47 1140.00 239.50 214.00 262.00 295.00 101.50 190.00 193.00 250.00 91.05 972.00 115.75 184.00 112.00 172.00 187.00 161.50 271.00 535.00 88.54 1092.00 183.25 277.00 229.00 253.00 100.00 103.00 238.00 475.00 92.12 996.00 133.75 184.00 286.00 235.00 181.00 179.50 250.00 542.50 90.39 1188.00 192.25 193.00 229.00 244.00 167.50 160.00 226.00 385.00 85.24 1212.00 178.75 184.00 205.00 268.00 106.00 173.50 229.00 362.50 101.29 1092.00 221.50 262.00 298.00 256.00 166.00 154.00 160.00 595.00 113.51 924.00 232.75 220.00 103.00 190.00 196.00 157.00 250.00 490.00 101.44 1068.00 156.25 181.00 295.00 211.00 134.50 175.00 286.00 527.50 87.46 1068.00 131.50 208.00 229.00 238.00 181.00 176.50 208.00 565.00 106.70 900.00 212.50 250.00 208.00 271.00 178.00 143.50 292.00 475.00 97.78 1092.00 194.50 289.00 100.00 118.00 157.00 199.00 202.00 445.00 100.42 1020.00 185.50 211.00 136.00 256.00 110.50 170.50 196.00 452.50 109.35 900.00 142.75 106.00 172.00 244.00 151.00 158.50 244.00 542.50 96.85 1044.00 140.50 259.00 178.00 187.00 158.50 146.50 259.00 302.50 72.81 1164.00 241.75 112.00 289.00 292.00 187.00 196.00 208.00 512.50 128.12 948.00 136.00 124.00 286.00 151.00 136.00 182.50 160.00 475.00 84.20 1140.00 127.00 133.00 247.00 193.00 152.50 151.00 232.00 325.00 66.41 1236.00 140.50 259.00 193.00 106.00 178.00 109.00 295.00 452.50 74.71 1068.00 138.25 172.00 289.00 184.00 182.50 149.50 241.00 587.50 109.38 1116.00 160.75 163.00 292.00 196.00 152.50 178.00 172.00 400.00 75.74 1044.00 185.50 292.00 262.00 118.00 119.50 151.00 214.00 355.00 87.43 1236.00 158.50 187.00 295.00 292.00 137.50 182.50 268.00 460.00 100.85 1140.00 176.50 118.00 223.00 253.00 161.50 175.00 247.00 512.50 97.07 948.00 104.50 235.00 289.00 235.00 172.00 166.00 292.00 460.00 70.57 1092.00 212.50 103.00 214.00 199.00 167.50 145.00 271.00 542.50 124.84 948.00 136.00 142.00 160.00 106.00 148.00 199.00 211.00 497.50 76.03 1116.00 187.75 241.00 223.00 100.00 125.50 136.00 238.00 392.50 91.06 948.00 214.75 292.00 163.00 250.00 124.00 115.00 223.00 340.00 99.88 996.00 176.50 226.00 250.00 121.00 130.00 124.00 229.00 347.50 72.98 1044.00 230.50 265.00 199.00 277.00 128.50 148.00 229.00 407.50 112.93 1020.00 120.25 187.00 187.00 286.00 191.50 194.50 214.00 317.50 55.23 1116.00 136.00 253.00 223.00 163.00 196.00 145.00 223.00 535.00 91.09 948.00 205.75 175.00 172.00 280.00 160.00 149.50 142.00 415.00 94.70 996.00 104.50 208.00 187.00 271.00 122.50 160.00 265.00 317.50 66.84 972.00 217.00 181.00 223.00 193.00 112.00 103.00 145.00 362.50 95.81 924.00 196.75 271.00 193.00 199.00 136.00 172.00 292.00 325.00 81.70 900.00 147.25 283.00 268.00 205.00 194.50 110.50 217.00 325.00 62.77 900.00 100.00 109.00 193.00 199.00 184.00 161.50 244.00 482.50 76.67 1044.00 239.50 265.00 139.00 115.00 119.50 110.50 298.00 557.50 128.72 1308.00 241.75 118.00 175.00 235.00 107.50 193.00 244.00 490.00 144.27 1044.00 129.25 259.00 169.00 130.00 182.50 160.00 214.00 295.00 59.53 1116.00 208.00 262.00 172.00 100.00 127.00 188.50 286.00 490.00 95.55 1164.00 142.75 163.00 187.00 271.00 140.50 191.50 229.00 182.50 57.25 948.00 127.00 196.00 220.00 202.00 178.00 178.00 292.00 475.00 94.08 972.00 106.75 100.00 289.00 127.00 146.50 185.50 190.00 490.00 75.09 1308.00 178.75 205.00 262.00 115.00 176.50 169.00 268.00 572.50 109.40 1020.00 129.25 214.00 292.00 163.00 116.50 131.50 223.00 257.50 68.03 1140.00 127.00 121.00 169.00 265.00 185.50 140.50 235.00 550.00 86.14 1140.00 239.50 262.00 142.00 124.00 125.50 190.00 223.00 512.50 124.84 113 1140.00 100.00 253.00 199.00 124.00 170.50 179.50 247.00 490.00 89.38 972.00 235.00 217.00 199.00 181.00 145.00 152.50 145.00 265.00 89.19 996.00 167.50 229.00 274.00 181.00 176.50 145.00 157.00 392.50 66.17 900.00 127.00 229.00 130.00 208.00 199.00 104.50 262.00 377.50 57.29 924.00 246.25 169.00 163.00 205.00 104.50 199.00 118.00 572.50 118.10 1308.00 214.75 217.00 160.00 130.00 160.00 158.50 271.00 557.50 122.33 1332.00 196.75 133.00 286.00 244.00 125.50 127.00 274.00 250.00 84.29 996.00 113.50 136.00 295.00 115.00 164.50 116.50 184.00 400.00 72.44 1068.00 113.50 163.00 124.00 259.00 185.50 155.50 244.00 280.00 49.06 948.00 156.25 235.00 262.00 115.00 197.50 149.50 178.00 302.50 63.76 1020.00 109.00 145.00 232.00 250.00 152.50 125.50 235.00 437.50 67.67 1188.00 223.75 226.00 283.00 298.00 155.50 130.00 199.00 512.50 118.00 1164.00 237.25 265.00 160.00 223.00 128.50 158.50 220.00 542.50 106.84 948.00 176.50 166.00 292.00 133.00 172.00 187.00 184.00 557.50 102.23 1308.00 102.25 265.00 205.00 121.00 182.50 193.00 283.00 362.50 72.95 972.00 212.50 295.00 289.00 163.00 175.00 167.50 262.00 385.00 89.26 1092.00 147.25 262.00 175.00 268.00 124.00 158.50 277.00 377.50 95.02 1068.00 174.25 109.00 178.00 268.00 155.50 146.50 229.00 557.50 118.26 1116.00 199.00 283.00 202.00 160.00 193.00 155.50 226.00 325.00 80.99 1092.00 248.50 160.00 262.00 112.00 119.50 134.50 202.00 527.50 128.77 900.00 217.00 196.00 196.00 223.00 197.50 191.50 247.00 580.00 119.60 1332.00 244.00 208.00 157.00 199.00 146.50 199.00 256.00 400.00 111.74 1332.00 230.50 289.00 130.00 172.00 134.50 137.50 298.00 325.00 89.85 1212.00 183.25 289.00 271.00 133.00 164.50 166.00 199.00 542.50 96.62 1260.00 237.25 160.00 223.00 142.00 131.50 116.50 241.00 542.50 127.44 948.00 203.50 163.00 184.00 184.00 143.50 113.50 160.00 497.50 101.05 1068.00 138.25 187.00 142.00 151.00 145.00 143.50 283.00 452.50 97.62 Thong so ung voi phuong an toi uu 1068.00 113.50 163.00 124.00 259.00 49.060 114 185.50 155.50 244.00 280.00