Tối ưu tiến độ về thời gian và chi phí dựa trên cơ sở phân phối tài nguyên phù hợp với năng lực thi công của nhà thầu

ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA  TRẦN KIẾN TƯỜNG Đề tài: TỐI ƯU TIẾN ĐỘ VỀ THỜI GIAN VÀ CHI PHÍ DỰA TRÊN CƠ SỞ PHÂN PHỐI TÀI NGUYÊN PHÙ HP VỚI NĂNG LỰC THI CÔNG CỦA NHÀ THẦU CHUYÊN NGÀNH: CÔNG NGHỆ & QUẢN LÝ XÂY DỰNG MÃ SỐ NGÀNH: 60.58.90 LUẬN VĂN THẠC SỸ PHẦN THUYẾT MINH Tp Hồ Chí Minh,03/2007 TRƯỜNG ĐẠI HỌC BÁCH KHOA CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM PHÒNG ĐÀO TẠO SĐH ĐỘC LẬP – TỰ DO – HẠNH PHÚC Tp HCM, ngày tháng năm 2007 NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên : TRẦN KIẾN TƯỜNG Phái: Nam Ngày, tháng, năm sinh : 24/07/1975 .Nôi sinh: Tp.Hồ Chí Minh Chuyên ngành : Công nghệ & quản lý xây dựng Mã số ngành: 60.58.90 I TÊN ĐỀ TÀI : TỐI ƯU TIẾN ĐỘ VỀ THỜI GIAN VÀ CHI PHÍ DỰA TRÊN CƠ SỞ PHÂN PHỐI TÀI NGUYÊN PHÙ HP VỚI NĂNG LỰC THI CÔNG CỦA NHÀ THẦU II NHIỆM VỤ VÀ NOÄI DUNG: Nhiệm vụ: Tối ưu tiến độ thi công dựa sở phân phối tài nguyên phù hợp với lực thi công (nguồn lực) nhà thầu Noäi dung: Chương 1: Mở Ñaàu Chương 2: Tổng Quan Chương 3: Phương Pháp Nghiên Cứu Chương 4: Tối Ưu Hóa Tiến Độ Mạng Theo Tài Nguyên Chương 5: Tối Ưu Hóa Tiến Độ Mạng Sao Cho Thời Gian Là Cực Tiểu Với Điều Kiện Tài Nguyên Giới Hạn Chương 6: Tối Ưu Hóa Tiến Độ Mạng Sao Cho Chi Phí Là Cực Tiểu Với Thời Gian Quy Ñònh Chương 7: Ví Dụ Minh Hoïa Chương 8: Kết Luận Và Kiến Nghò III NGÀY GIAO NHIỆM VỤ: IV NGÀY HOÀN THÀNH NHIỆM VỤ: V HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN: TS NGÔ QUANG TƯỜNG CÁN BỘ HƯỚNG DẪN CHỦ NHIỆM NGÀNH BỘ MÔN QUẢN LÝ CHUYÊN NGÀNH TS NGÔ QUANG TƯỜNG TS NGÔ QUANG TƯỜNG TS NGÔ QUANG TƯỜNG Nội dung đề cương luận văn thạc só Hội đồng chuyên ngành thông qua Ngày .tháng năm 2007 PHÒNG ĐÀO TẠO SĐH KHOA QUẢN LÝ NGÀNH CÔNG TRÌNH ĐƯC HOÀN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐẠI HỌC QUỐC GIA TP.HỒ CHÍ MINH Cán hướng dẫn khoa học: TS NGÔ QUANG TƯỜNG Cán chấm nhận xeùt 1: Cán chấm nhận xét 2: Luận văn thạc sỹ bảo vệ 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 LỜI CẢM ƠN Luận văn thạc sỹ hoàn thành nhờ quan tâm giúp đỡ từ gia đình, thầy cô, bạn bè đồng nghiệp Xin bày tỏ lòng biết ơn chân thành đến thầy TS.Ngô Quang Tường tận tình giúp đỡ, động viên tạo điều kiện thuận lợi suốt trình thực luận văn Xin gởi lời cảm ơn đến Thầy Cô ngành Công nghệ & Quản lý Xây Dựng, Phòng Đào tạo Sau Đại học trường Đại học Bách Khoa Tp.Hồ Chí Minh hết lòng giảng dạy, hướng dẫn hỗ trợ tận tình suốt khóa học Xin chân thành cảm ơn bạn KS.Nguyễn Khánh Chí, ThS.KS.Trịnh Tuấn, tất bạn bè đồng nghiệp nhiệt tình giúp đỡ, động viên để hoàn thành luận văn Tp.Hồ Chí Minh, ngày 05 tháng 03 năm 2007 Người thực Trần Kiến Tường TÓM TẮT NỘI DUNG Mục tiêu thực dự án đạt thành cao với điều kiện dự án phải hoàn thành khoảng thời gian ngân sách giới hạn Đối với đơn vị chủ đầu tư, thành cao chất lượng công trình bao gồm công năng, thẫm mỹ, độ bền Đối với đơn vị tư vấn, thành nghóa chất lượng công trình hài lòng chủ đầu tư Đối với nhà thầu, thành cao sử dụng hiệu nguồn lực (năng lực) đảm bảo chất lượng công trình Nguồn lực nhà thầu bao gồm nhiều thành phần, từ công nhân có tay nghề khác đến loại máy thi công, lực tài Do có nhiều thành phần nên việc sử dụng hiệu nguồn lực vấn đề khó khăn cho nhà thầu Trong luận văn này, tác giả nghiên cứu biện pháp phân phối tài nguyên bao gồm nhân công, máy móc thi công cho tổ đội, công tác cho suất lao động cao nhất, đồng thời tổng tài nguyên cần thiết đáp ứng tiến độ không vượt qua nguồn lực giới hạn nhà thầu tiến độ đảm bảo tiêu kinh tế kỹ thuật ban đầu (tài nguyên, thời gian, chi phí) Luận văn giải vấn đề sau: (1)Với nguồn lực (năng lực) có giá trị giới hạn đó, phân phối tài nguyên cho tổ đội để thời gian hoàn thành công trình ngắn nhất; (2)Với thời gian hoàn thành giới hạn mà chủ đầu tư quy định, phân phối tài nguyên cho tổ đội để chi phí thực công trình nhỏ nhất; (3)Sau phân phối tài nguyên cho tổ đội, xác định thời điểm khởi công tác cho tổng tài nguyên sử dụng biểu đồ tài nguyên điều hòa đảm bảo thời gian thi công công trình ABSTRACT The goal of executing a project is obtaining the best achievement so that the project was finished in restricted duration and limited budget Following the Owner, the best achievement is the good quality of project Following the Consultants, it is the good quality of project and satisfying the Owner Following the Contractors, it is using their capacities effectively Their capacities consist of many kinds of workers, many kind of machine and their finance Therefor, using capacity effectively is a difficult problem of Contractor In this essay, I reseach the method to distribute resource that consist of workers, machine to each team and each task so that productivity is highest, at the same time general necessary resource is less than the limited capacity of the Contractor, but the schedule will reach the target of resource, duration or expense There are some problems solved in the essay as below: (1) With in the limtited capacity, how to distribute resource to each team so that complection duration of the project is minimum; (2) With in the restricted duration, how to distribute resource to each team so that the expense of the project is minimum; (3) After each team is distributed resource, selecting starting time of each task so that general necessary resource is the least or chart of resouce becomes more regular Muïc Luïc Luận án cao học- GVHD: TS.Ngô Quang Tường MỤC LỤC _ Trang Trang bìa Nhiệm vụ luận án Lời cảm ơn Tóm tắt nội dung luận văn Muïc luïc Danh sách ký hiệu Chương 1: MỞ ÑAÀU 1.1 Đặt vấn đề 1.2 Mục tiêu nghiên cứu 1.3 Phạm vi nghiên cứu 1.4 Đối tượng nghiên cứu 1.5 Phương pháp công cụ nghiên cứu Chương 2: TỔNG QUAN 2.1 Kỹ thuật lập tiến độ thi công 2.2 Những nghiên cứu liên quan đến tiến độ thi công 2.3 Những phương pháp toán kinh tế áp dụng việc quản lý xây dựng Chương 3: PHƯƠNG PHÁP NGHIÊN CỨU 3.1 Lý thuyết sơ đồ mạng song lập (PDM) 3.2 Toán học ứng dụng máy tính việc tính toán sơ đồ mạng 10 3.3 Ngôn ngữ lập trình MatLab 17 3.4 Trình tự tối ưu tiến độ thi coâng 19 3.5 Dữ liệu toán 23 Chương 4: TỐI ƯU HÓA TIẾN ĐỘ MẠNG THEO TÀI NGUYÊN 27 4.1 Đặt vấn ñeà 27 4.2 Giải vấn ñeà 27 Chương 5: TỐI ƯU HÓA TIẾN ĐỘ MẠNG SAO CHO THỜI GIAN LÀ CỰC TIỂU VỚI ĐIỀU KIỆN TÀI NGUYÊN GIỚI HẠN 54 5.1 Đặt vấn đề 54 5.2 Giải vấn đề 54 Chương 6: TỐI ƯU HÓA TIẾN ĐỘ MẠNG SAO CHO CHI PHÍ LÀ CỰC TIỂU VỚI THỜI GIAN QUY ĐỊNH 65 6.1 Đặt vấn đề 65 6.2 Giải vấn đề 65 HVTH: KS.Trần Kiến Tường Mục Lục Luận án cao học- GVHD: TS.Ngô Quang Tường Chương 7: VÍ DỤ MINH HOÏA 78 7.1 Dữ liệu chung 78 7.2 Dữ liệu riêng 90 7.3 Kết 98 Chương 8: KẾT LUẬN VÀ KIẾN NGHỊ 111 8.1 Kết luận 111 8.2 Kiến nghị 112 TÀI LIỆU THAM KHẢO 113 TOÙM TẮT LÝ LỊCH CÁ NHÂN HVTH: KS.Trần Kiến Tường Luận án cao học- GVHD: TS.Ngô Quang Tường CÁC KÝ HIỆU, MA TRẬN VÀ HÀM SỐ CÁC KÝ HIỆU: : số ca thực công tác i ngày; CBQL : số lượng trung bình cán quản lý công trình (người); : số lượng công nhân trung bình làm việc công trường (người); : tổng chi phí làm lán trại cho công nhân cho công trình(VND); : tổng chi phí quản lý công trình(VND); : tổng chi phí tăng ca công nhân làm ca/ngày công trình(VND); : tổng chi phí tuyển dụng công nhân cho công trình(VND); : tổng chi phí thuê mướn máy móc cho công trình(VND); : tổng chi phí tài nguyên (nhân công, máy thi công) trường hợp không kể chi phí tăng ca, thuê mướn công trình(VND); : tổng chi phí vận chuyển thiết bị vào công trường cho công trình(VND); : chi phí lán trại cho công nhân có tay nghề I/7(VND/người); CNTB CPLT CPQL CPTC CPTD CPTM CPTN CPVC DGLTNI DGQL DGTDNI : đơn giá quản lý, tức mức lương trung bình cán quản lý công trường(VND/tháng); : đơn giá tuyển dụng, tính theo bậc thợ (VND/người); DGTM MI : chi phí thuê mướn máy thi công có số thứ tự I bảng chủng loại máy DGTN MI thi công (VND/ca); : chi phí khấu hao (đơn giá) máy thi công có số thứ tự I bảng chủng DGTN NI loại máy thi công, cập nhật theo thời gian (VND/ca); : đơn giá công nhân bậc I/7, cập nhật theo thời gian (VND/công); DGTNVtuI : đơn giá vật tư có số thứ tự MI bảng chủng loại vật tư, cập nhật theo thời gian (VND/đơn vị tính); DGToSX i : chi phí cho tổ sản xất thực công tác i ngaøy (VND/ngaøy); DGToSX N & M ,i : chi phí nhân công máy thi công cho tổ sản xất thực DGVCMI công tác i ngày (VND/ngày); : chi phí vận chuyển máy thi công có số thứ tự I bảng chủng loại máy dmVtuI ,i thi công từ nơi cho thuê đến công trường(VND/ca); : định mức tiêu thụ vật tư thứ I (trong bảng chủng loại vật tư) công tác i; DM NI ,i : số lượng công nhân bậc I/7 tổ đội sản xuất sử dụng để DM MI ,i thực công tác i ca; : số lượng ca máy loại I tổ đội sản xuất sử dụng để thực công tác i ca; HVTH: KS.Trần Kiến Tường Luận án cao học- GVHD: TS.Ngô Quang Tường DM VtuI ,i : số lượng vật tư thứ I tổ đội sản xuất thực công tác i k1,i ca; : hệ số bất điều hòa tài nguyên thứ i, không kể vật tư; k2,i : hệ số phân bố tài nguyên tài nguyên thứ i, không kể vật tư; Ki Ki : thời điểm kết thúc công tác i; Khi : thời điểm khởi công công tác i; KhM i : thời điểm khởi muộn công tác i, thời điểm công tác i bắt đầu KhSi muộn mà kết thúc không ảnh hưởng đến thời điểm kết thúc sơ đồ mạng; : thời điểm khởi sớm công tác i, thời điểm mà công tác i bắt KM i đầu sớm nhất; : thời điểm kết muộn công tác i, thời điểm công tác i kết thúc KSi muộn mà không ảnh hưởng đến kết thúc công tác cuối cùng; : thời điểm kết sớm công tác i, thời điểm mà công tác i kết thúc sớm nhất; MaxBDTN MI : giá trị cực đại số lượng máy móc thi công loại I biểu đồ tài nguyên; MaxBDTN NI : giá trị cực đại số công nhân bậc I biểu đồ tài nguyên; MaxTN MI : tổng số máy móc thi công loại I lớn sử dụng công MaxTN NI trường; : tổng số công nhân bậc I lớn làm việc công trường; MaxTongCN : tổng số công nhân lớn làm việc công trường; q0i : khối lượng công tác i tổ sản xuất thực ca làm Qi việc (m, m2, m3, tấn, /ca); : khối lượng công tác i (m, m2, m3, tấn, ); Si : tổng số công (ca) tài nguyên thứ i, không kể vật tư (công/ca); S du ,i : tổng số công (ca) vượt trội so với đường trung bình tài nguyên thứ i, ti không kể vật tư (công/ca); : khoảng thời gian thi công công tác i (ngày); : thời gian hoàn thành công trình (ngày); : thời gian quy định hoàn thành công trình (ngày); T ( Ln ) T ( Ln ) : thời gian đường xuyên mạng cần tìm hệ số găng; T [T] Tg ( Ln ) : tổng thời gian thi công công tác găng nằm đường này; Tg : thời gian đường găng; TN NI ,i : số lượng công nhân bậc I/7 thực công tác i; TN MI ,i : số lượng thiết bị máy móc có số thứ tự I bảng chủng loại máy thi TNVtuI ,i công sử dụng để thực công tác i hàng ngày; : lượng vật tư có số thứ tự I bảng chủng loại vật tư sử dụng cho HVTH: KS.Trần Kiến Tường Chương Phụ lục Luận án cao học- GVHD: TS.Ngô Quang Tường function [Kh, BDTNi, k1, k2] = baitoan2() %Bai toan 2: Xac dinh thoi diem khoi cua cac cong tac cho bieu tai nguyen dieu %hoa nhat voi dieu kien phan cac to san xuat co ban cua cac to doi da duoc xac dinh %(tuc da cho phuong an phan phoi tai nguyen) va moi quan he giua cac cong tac khong %thay doi global SoTN; global N; global M; global SoCtac; %Du lieu chung: N = 3; M = 8; %Du lieu rieng: Q = data_Q(); q0 = data_q0(); ToSX = data_ToSX(); a = data_a(); SoCtac = size(Q,1); DM = data_dinhmuc(); [SS, FS, FF] = quanhe(SoCtac); SoTN = N + M; %Xac dinh cac thong so cua ham patternsearch: TSCT = thongsoctac(Q, q0, ToSX, a, FS, SS, FF); KhS = TSCT(:,2); KhM = TSCT(:,3); [A, B] = dkkhoi(Q, q0, ToSX, a, FS, SS, FF); Kh = []; %Xac dinh ma tran Kh bieu tai nguyen thu i, TNi, dieu hoa: err = 0; while err == st = sprintf('Yeu cau dieu hoa so cong nhan tay nghe nao? (1 %d)', N); answer = inputdlg(st, '', 1); if isempty(answer) err = 1; return; end answer = strtrim(answer); TNi = str2double(answer); if 0 plot_bdtn(Q, q0, ToSX, a, DM, Kh, TNi); else continue; end end HVTH: KS.Traàn Kiến Tường Trang 32 Chương Phụ lục Luận án cao học- GVHD: TS.Ngô Quang Tường 2.3 CÁC CHƯƠNG TRÌNH HÀM TRONG CHƯƠNG 6: function [TCP] = tongchiphi(Q, q0, a, DM, FS, SS, FF, TenToSX, ToDoi, DayChuyen, DGTN, DGLT, DGVC, DGQL, DGTD, DGTM, DGTNTC,TNCPhep) %Ham xac dinh tong chi phi thi cong cua cong trinh da biet so to san xuat co ban cua %cac to doi, tuc da chon phung an phan phoi tai nguyen cho cac cong tac global N; global M; DayChuyen(:, 2) = ToDoi; SoCtac = size(Q,1); %Xac dinh gia tri ma tran [ToSX]: for i=1:1:SoCtac [row, col] = find(DayChuyen(:, 1) == TenToSX(i, 1)); ToSX(i, 1) = DayChuyen(row, 2); end %Xac dinh ma tran cuc tieu dai dien ma tran [MinMaxBDTN]: [MinMaxBDTN, Khoi] = minmaxbdtn(Q, q0, ToSX, a, DM, FS, SS, FF); GTMaxBDTN = DGTN*MinMaxBDTN; [i, j] = min(GTMaxBDTN); Kh = Khoi(:, j); BDTN = bdtn(Q, q0, ToSX, a, DM, Kh); %Xac dinh ma tran bieu tai nguyen nhan cong va may thi cong tu [BDTN]: SoTN = N + M; T = size(BDTN, 2); BDNC = BDTN(1:N, :); BDMM = BDTN(N + 1:SoTN, :); MaxBDNC = max(BDNC')'; MaxBDMM = max(BDMM')'; %Xac dinh chi phi: CPTN = DGTN*BDTN*repmat(1,T,1); CPLT = DGLT*MaxBDNC; CPVC = DGVC*MaxBDMM; CPQL = 10*DGQL*(repmat(1,1,N)*BDNC*repmat(1,T,1)/30)/100; CPTD = 0; for NI=1:1:N if MaxBDNC(NI,1) > TNCPhep(NI, 1) HVTH: KS.Trần Kiến Tường Trang 33 Chương Phụ lục Luận án cao học- GVHD: TS.Ngô Quang Tường CPTD = CPTD + DGTD(1, NI)*(MaxBDNC(NI, 1) - TNCPhep(NI, 1)); end end CPTM = 0; for MI=1:1:M if MaxBDMM(MI,1) > TNCPhep(N + MI, 1) CPTM = CPTM + DGTM(1,MI)*(MaxBDMM(MI, 1) - TNCPhep(N + MI, 1)); end end CPTC = 0; TN = []; for i=1:1:SoCtac if a(i, 1) >= TN(i, :) = DM(i, :)*ToSX(i, 1); CPTC = CPTC + (a(i, 1)-1)*DGTNTC*TN(i,:)'; end end TCP = CPTN + CPLT + CPVC + CPQL + CPTC + CPTD + CPTM; HVTH: KS.Trần Kiến Tường Trang 34 Chương Phụ lục Luận án cao học- GVHD: TS.Ngô Quang Tường function [T, DayChuyen, BDTN, TSCT, MinTongCP, Kh, a] = baitoan4() %Bai toan 3: Xac dinh phuong an phan phoi tai nguyen cho cac cong tac cho chi phi %thi cong cong trinh la cuc tieu voi dieu kien moi quan he giua cac cong tac khong thay %doi va thoi gian hoan cong trinh khong lon hon thoi gian cho phep global N; global M; global SoTN; %Du lieu chung: N = 3; M = 8; DGTN = data_DGTN(); DGTD = data_DGTD(); DGTM = data_DGTM(); DGVC = data_DGVC(); DGLT = data_DGLT(); DGQL = data_DGQL(); DGTNTC = [1.5*DGTN(:, 1:N) DGTN(:, (N+1):(N+M))]; %Du lieu rieng: Q = data_Q(); q0 = data_q0(); ToSX = data_ToSX(); TenToSX = data_TenToSX(); ToSXMin = data_ToSXMin(); ToSXMax = data_ToSXMax(); a = data_a(); SoCtac = size(Q,1); DM = data_dinhmuc(); [SS, FS, FF] = quanhe(SoCtac); SoTN = N + M; %Gioi han cho phep: TNCPhep = data_tnguyenchophep(); TCPhep = data_thoigianchophep(); DayChuyen = daychuyen(TenToSX, ToSXMin, ToSXMax); err = 0; while err == [MinT, DayChuyen, BDTN, TSCT, Kh] = cuctieuthoigian(Q, q0, FS, SS, FF, ToSXMin, ToSXMax, TenToSX, a, DM, DGTN, TNCPhep); ToDoiMin = DayChuyen(:, 3); ToDoiMax = DayChuyen(:, 4); HVTH: KS.Trần Kiến Tường Trang 35 Chương Phụ lục Luận án cao học- GVHD: TS.Ngô Quang Tường ToDoiMinT = DayChuyen(:,2); options = psoptimset('TolCon',0.5,'TolMesh',0.5,'TolFun',1000,'TolX',1,'Cache','On', 'CompletePoll','On','CompleteSearch','On','SearchMethod','GPSPositiveBasis2N'); %Kiem tra dieu kien ve thoi gian yeu cau voi kha nang thuc hien: if MinT = DayChuyen(:, 2) = (ToDoi); for i=1:1:SoCtac [row, col] = find(DayChuyen(:, 1) == TenToSX(i, 1)); ToSX(i, 1) = DayChuyen(row, 2); end [MinMaxBDTN, Khoi] = minmaxbdtn(Q, q0, ToSX, a, DM, FS, SS, FF); GTMaxBDTN = DGTN*MinMaxBDTN; [i, j] = min(GTMaxBDTN); Kh = Khoi(:, j); BDTN = bdtn(Q, q0, ToSX, a, DM, Kh); TSCT = thongsoctac(Q, q0, ToSX, a, FS, SS, FF); T = max(TSCT(:,5)); %Ve tien & bieu tai nguyen tuong ung vo ma tran [Kh] plot_td(Q, q0, ToSX, a, FS, SS, FF, Kh); for TNi=1:1:SoTN DMTNi = DM(:,TNi); if max(DMTNi')>0 plot_bdtn(Q, q0, ToSX, a, DM, Kh, TNi); else continue; end end return; else err = 2; msgbox('Bai toan khong hoi tu!'); break; end else HVTH: KS.Trần Kiến Tường Trang 36 Phụ lục Luận án cao học- GVHD: TS.Ngô Quang Tường Chương button = questdlg('Thoi gian quy dinh nho hon thoi gian cho phep, co dong y huy dong them cong nhan & may moc khong? (Y/N)', '', 'Yes', 'No', 'Yes'); if strcmpi(button, 'Yes') objfun = @(ToDoi) tongchiphi(Q, q0, a, DM, FS, SS, FF, TenToSX, ToDoi, DayChuyen, DGTN, DGLT, DGVC, DGQL, DGTD, DGTM, DGTNTC,TNCPhep); nonlcon = @(ToDoi) dkthoigian(Q, q0, a, FS, SS, FF, TenToSX, ToDoi, DayChuyen, TCPhep); [ToDoi, MinTongCP, exitflag] = patternsearch(objfun, ToDoiMax, [], [], [], [], ToDoiMin, ToDoiMax, nonlcon, options); if exitflag >= DayChuyen(:, 2) = round(ToDoi); for i=1:1:SoCtac [row, col] = find(DayChuyen(:, 1) == TenToSX(i, 1)); ToSX(i, 1) = DayChuyen(row, 2); end [MinMaxBDTN, Khoi] = minmaxbdtn(Q, q0, ToSX, a, DM, FS, SS, FF); GTMaxBDTN = DGTN*MinMaxBDTN; [i, j] = min(GTMaxBDTN); Kh = Khoi(:, j); BDTN = bdtn(Q, q0, ToSX, a, DM, Kh); TSCT = thongsoctac(Q, q0, ToSX, a, FS, SS, FF); T = max(TSCT(:,5)); %Ve tien & bieu tai nguyen tuong ung voi ma tran [Kh]: plot_td(Q, q0, ToSX, a, FS, SS, FF, Kh); for TNi=1:1:SoTN DMTNi = DM(:,TNi); if max(DMTNi')>0 plot_bdtn(Q, q0, ToSX, a, DM, Kh, TNi); else continue; end end return; else err = 2; msgbox('Bai toan khong hoi tu!'); end else button = questdlg('Thoi gian quy dinh qua nho, co dong y giai phap tang ca hay khong? (Y/N)', '', 'Yes', 'No', 'Yes'); if strcmpi(button, 'Yes') HVTH: KS.Trần Kiến Tường Trang 37 Phụ lục Luận án cao học- GVHD: TS.Ngô Quang Tường Chương %Tao day doi ISO chua STT & DGToSXTC cua cac cong tac co so ca lam viec