Đang tải... (xem toàn văn)
BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA TPHCM BÁO CÁO BÀI TẬP LỚN VẬT LÝ 1 Tên đề tài XÁC ĐỊNH CÔNG CỦA HỆ TRONG CÁC QUÁ TRÌNH CÂN BẰNG TỪ GIẢN ĐỒ (p,V) Giảng viên hướng dẫn ThS Dương Thị Như[.]
BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA TPHCM BÁO CÁO BÀI TẬP LỚN VẬT LÝ Tên đề tài: XÁC ĐỊNH CÔNG CỦA HỆ TRONG CÁC QUÁ TRÌNH CÂN BẰNG TỪ GIẢN ĐỒ (p,V) Giảng viên hướng dẫn: ThS Dương Thị Như Tranh Sinh viên thực hiện: Hồng Minh Qn Võ Trí Nhân Lê Dương Qua 4.Võ Thị Thao Vũ Minh Phương Nguyễn Ngọc Huỳnh Như Vũ Nguyễn Tố Như Đỗ Tuấn Thanh Ngô Phú Lộc TpHCM, 12 - 2013 Luan van LỜI CẢM ƠN Trong suốt trình thực tiểu luận nói trên, nhóm chúng tơi nhận nhiều quan tâm ủng hộ, giúp đỡ tận tình thầy cơ, anh chị em bè bạn Ngồi ra, nhóm xin gửi lời tri ân chân thành đến cô Dương Thị Như Tranh, giảng viên hướng dẫn cho đề tài matlab Nhờ có hết lịng bảo mà nhóm hoàn thành tiểu luận tiến độ giải tốt vướng mắc gặp phải Sự hướng dẫn cô kim nam cho hành động nhóm phát huy tối đa mối quan hệ hỗ trợ thầy trị mơi trường giáo dục Lời cuối, xin lần gửi lời biết ơn sâu sắc đến cá nhân, thầy dành thời gian dẫn cho nhóm Đây niềm tin, nguồn động lực to lớn để nhóm đạt kết Nhóm thực đề tài matlab Hoàng Minh Quân Võ Trí Nhân Lê Dương Qua 4.Võ Thị Thao Vũ Minh Phương Nguyễn Ngọc Huỳnh Như Vũ Nguyễn Tố Như Ngô Phú Lộc Đỗ Tuấn Thanh Luan van MỤC LỤC TRANG LỜI CẢM ƠN MỤC LỤC DANH MỤC CÁC BẢNG BIỂU VÀ HÌNH VẼ Chương 1: NGUYÊN LÍ THỨ I CỦA NHIỆT ĐỘNG LỰC HỌC - CÁC KHÁI NIỆM LIÊN QUAN, CÔNG CỦA HỆ VÀ CƠNG TRONG CHU TRÌNH 1.1 KHÁI NIỆM VỀ HỆ NHIỆT ĐỘNG, TRẠNG THÁI CÂN BẮNG - QUÁ TRÌNH CÂN BẰNG 1.1.1 Hệ nhiệt động 1.1.2 Trạng thái cân - trình cân 1.1.2.1 Khái niệm trạng thái cân 1.1.2.2 Khái niện trình cân 1.1.2.3 Khái niệm trình thuận nghịch 1.2 NGUYÊN LÍ CỦA NHIỆT ĐỘNG LỰC HỌC, ỨNG DỤNG CỦA NGUN LÍ ĐỂ TÍNH CƠNG 1.2.1 Ngun lí Nhiệt động lực học 1.2.2 Khái niệm công 1.2.3 Ứng dụng định luật để tính cơng thực q trình 1.2.4 Cơng q trình đẳng áp, đẳng tích, đẳng nhiệt Chương 2: GIẢN ĐỒ PV VÀ CÁC HÀM CƠ BẢN TRONG MATLAB 2.1 GIẢN ĐỒ TRẠNG THÁI Luan van 2.1.1 Khái niệm giản đồ trạng thái 2.1.2 Phương trình trạng thái khí lí tưởng - giản đồ pV 2.2 CÁC HÀM CƠ BẢN MATLAB ĐƯỢC SỬ DỤNG 2.2.1 Tổng quan MATLAB 2.2.2 Các hàm matlab sử dụng toán Chương 3: THỨ TỰ GIẢI BÀI TỐN TÌM CƠNG TỪ GIẢN ĐỔ PV TRÊN MATLAB VÀ CODE HOÀN CHỈNH 3.1 THỨ TỰ GIẢI BÀI TỐN TÌM CƠNG TỪ GIẢN ĐỒ PV 11 3.2 CODE HOÀN CHỈNH (KÈM LỜI DỊCH CODE) 12 DANH MỤC CÁC TÀI LIỆU THAM KHẢO Luan van DANH MỤC CÁC BẢNG BIỂU VÀ HÌNH VẼ Tên bảng biểu hình vẽ Trang Hình 1.1: Hình 2.1: Hình 2.2: Hình 2.3: Luan van LỜI MỞ ĐẦU Vật lý đại cương môn học đại cương có tầm quan trọng sinh viên ĐH Bách Khoa TPHCM nói riêng sinh viên ngành khối khoa học kỹ thuật – cơng nghệ nói chung Do đó, việc dành cho mơn học khối lượng thời gian định thực hành điều tất yếu để giúp cho sinh viên có sở vững môn KHTN làm tiền đề để học tốt mơn khác chương trình đào tạo Sự phát triển toán tin đời hỗ trợ lớn trình phát triển môn học vật lý Việc ứng dụng tin học q trình giải thích sở liệu vật lý, giải toán vật lý làm cho thời gian bỏ rút ngắn lại mang hiệu cao Như ta biết, phần mềm ứng dụng Matlab giải vấn đề Vì việc tìm hiểu matlab ứng dụng matlab việc thực hành môn học vật lý đại cương quan trọng có tính cấp thiết cao Ở tập lớn này, nhóm thực nội dung “Tính cơng q trình cân dựa giản đồ pV” thông qua phần mềm Matlab Đây dạng toán quan trọng phần Nhiệt học, ứng dụng nhiều môn nghiên cứu khác nhứ: Hóa Lý, Vậy Lý học, nghiên cứu giản đồ chuyển pha, nhiệt học… Sau nội dung tìm hiểu tập lớn nhóm! Luan van Chương 1: NGUYÊN LÍ THỨ I CỦA NHIỆT ĐỘNG LỰC HỌC - CÁC KHÁI NIỆM LIÊN QUAN, CÔNG CỦA HỆ VÀ CƠNG TRONG CHU TRÌNH 1.1 KHÁI NIỆM VỀ HỆ NHIỆT ĐỘNG, TRẠNG THÁI CÂN BẮNG - QUÁ TRÌNH CÂN BẰNG 1.1.1 Hệ nhiệt động Một tập hợp vật xác định hồn tồn thơng số vĩ mô, độc lập với nhau, gọi hệ vĩ mô hay hệ nhiệt động (gọi tắt hệ) Các vật hệ ngoại vật hệ hay môi trường xung quanh hệ Nếu hệ môi trường khơng trao đổi nhiệt hệ lập ngoại vật phương diện nhiệt: ta nói hệ ngoại vật có vỏ cách nhiệt Nếu hệ ngoại vật trao đổi nhiệt không sinh cơng nén dãn nở hệ cô lập ngoại vật phương diện học Hệ gọi lập hồn tồn khơng tương tác trao đổi lượng với mơi trường ngồi 1.1.2 Trạng thái cân - q trình cân 1.1.2.1 Khái niệm trạng thái cân Trong nhiệt động lực học khái niệm trạng thái cân hệ trạng thái đại lượng vĩ mô (p, V, T) xác định trạng thái hệ không thay đổi Những đại lượng xác định trạng thái vật gọi thông số trạng thái Ở trạng thái cân nhiệt động lực học xảy tượng truyền nhiệt, phản ứng hóa học, biến đổi trạng thái khí, lỏng, rắn Trạng thái cân nhiệt động lực học khác với trạng thái cân học chỗ đại lượng vĩ mô đặc trưng cho hệ không đổi phần tử cấu tạo nên hệ không ngừng chuyển động hỗn loạn Chẳng hạn hệ gồm chất lỏng, đựng Luan- 1van bình kín, mặt chất lỏng có bão hồ Hệ trạng thái cân nên đại lượng p, V, T khơng đổi Vậy rút định nghĩa đầy đủ trạng thái cân nhiệt động lực học Ðó trạng thái hệ mà thông số trạng thái hệ không thay đổi trạng thái hệ không thay đổi, hệ khơng xảy q trình dẫn nhiệt, khuếch tán, phản ứng hóa học, chuyển pha.v.v 1.1.2.2 Khái niện trình cân Khi hệ biến đổi từ trạng thái sang trạng thái khác, chuổi trạng thái nối tiếp xảy ra, tạo nên trình Những trạng thái nối tiếp trạng thái cân biến thiên thông số trạng thái theo thời gian đủ chậm so với khoảng thời gian hai trạng thái chọn tuỳ ý Một q trình diễn biến vơ chậm gọi trình chuẩn cân (chuẩn tĩnh) coi dãy nối tiếp trạng thái cân Những trình xảy thực tế khơng phải q trình chuẩn cân chúng xảy chậm gần trình chuẩn cân nhiêu 1.1.2.3 Khái niệm trình thuận nghịch Trong nhiệt động lực học, xét q trình nói chung mà ta cần ý đến chiều diễn biến trình Vì ta xét khái niệm trình thuận nghịch Quá trình thuận nghịch trình diễn biến theo hai chiều, lúc đầu q trình diễn theo chiều (chiều thuận) sau lại diễn theo chiều ngược lại để trở trạng thái ban đầu hệ qua trạng thái giống lúc hệ diễn biến theo chiều thuận hệ trở trạng thái ban đầu khơng gây biến đổi cho ngoại vi Luan- 2van Mọi trình thuận nghịch q trình chuẩn cân Ta biểu diễn trình thuận nghịch đồ thị đường cong liền nét trình chuẩn cân 1.2 NGUN LÍ CỦA NHIỆT ĐỘNG LỰC HỌC, ỨNG DỤNG CỦA NGUYÊN LÍ ĐỂ TÍNH CƠNG 1.2.1 Ngun lí Nhiệt động lực học Nguyên lý thứ nhiệt động học trường hợp riêng định luật bảo toàn chuyển hóa lượng vận dụng vào q trình vĩ mơ Phát biểu: “Độ biến thiên lượng tồn phần ΔW hệ trình biến đổi vĩ mơ có giá trị tổng cơng A nhiệt Q mà hệ nhận trình đó” ΔW = A + Q Ở ta giả thuyết hệ không đổi (Wđ + Wt = const) ΔW = ΔU nên viết lại: ΔU = A + Q “Trong trình biến đổi, độ biến thiên nội hệ có giá trị tổng cơng nhiệt mà hệ nhận q trình đó” 1.2.2 Khái niệm công Công dạng truyền lượng làm tăng mức độ chuyển động có trật tự vật Điều xảy tương tác vật vĩ mơ, nghĩa vật có kích thước lớn kích thước phân tử nhiều Thí dụ: Khí giãn nở xylanh làm piston dịch chuyển sinh công (làm quay bánh xe, quay máy phát điện ), lượng hệ giảm Ta nói khối khí sinh cơng (A V2 > V1 δA A có giá trị dương Ngược lại hệ giảm thể tích tức dV < V2 < V1 δA A có giá trị âm Trường hợp biểu thị đồ thị δA > A > (chú ý lập luận nói ta coi áp suất p đại lượng vô hướng) Cơng q trình thuận nghịch lớn cơng q trình khơng thuận nghịch Để hiểu rõ điều ta xét nén giãn khí xi lanh: Gọi p áp suất ngoại vật tác dụng lên chất khí, px áp suất chất khí tác dụng lên ngoại vật Luan- 4van Nếu q trình nén giãn chất khí thuận nghịch p = pk Vậy cơng ngoại vật tác dụng lên chất khí q trình nén thuận nghịch (từ thể tích V2 giảm xuống thể tích V1) ta kí hiệu ANTN 𝑉1 𝑉1 𝐴NTN = ∫𝑉2 𝑝𝑑𝑉 = ∫𝑉2 pk𝑑𝑉 Cịn cơng ngoại vật nhận chất khí q trình giãn thuận nghịch (từ thể tích V1 tăng đến thể tích V2) ta ký hiệu AGTN, 𝑉2 𝑉2 AGTN = ∫𝑉2 𝑝𝑑𝑉 = ∫𝑉1 pk𝑑𝑉 So sánh hai công ta thấy chúng giá trị Nếu q trình nén giãn chất khí không thuận nghịch, hiểu nén giãn nhanh) nén chất khí, phải nén nhanh nên p > px cịn giãn chất khí phải giãn nhanh, nên p < px Vậy ta thấy q trình khơng thuận nghịch, để nén khí ta phải dùng cơng có giá trị lớn so với trình thuận nghịch, để khí dãn cơng mà ta thu lại nhỏ so với trình thuận nghịch Vậy so với trình thuận nghịch trình khơng thuận nghịch ta khơng lợi cơng Nói cách khác trình thực thực gần thuận nghịch ta lợi công nhiêu Dưới ta tính cơng số q trình Ta cần ý có hai cách tính cơng: dựa vào cơng thức, dựa vào biểu thức giải tích nguyên lý thứ Tùy trường hợp cụ thể để thuận tiện cho việc tính tốn, ta chọn hai cách nói 1.2.4 Cơng q trình đẳng áp, đẳng tích, đẳng nhiệt Ðể tính cơng q trình này, ta dựa vào cơng thức * Ðối với q trình đẳng tích (dV = 0): 𝐴 = ∫ 𝑝𝑑𝑉 = Luan- 5van * Ðối với trình đẳng áp(p = const): 𝑉2 𝑉2 𝐴 = ∫𝑉1 𝑝𝑑𝑉 = 𝑝 ∫𝑉2 𝑑𝑉 = 𝑝(𝑉2 − 𝑉1) * Ðối với trình đẳng nhiệt (T = const): 𝑉2 𝑉2 𝑇 𝑉2 𝐴 = ∫𝑉1 𝑝𝑑𝑉 = ∫𝑉1 𝑛𝑅 𝑑𝑉 = 𝑛𝑅𝑇𝑙𝑛 𝑉 𝑉1 Vậy truyền cho khí phần lượng hình thức cơng sau khí phải truyền phần lượng cho ngoại vật hình thức nhiệt Hai phần lượng mà khí nhận vào truyền có giá trị Chương 2: GIẢN ĐỒ PV VÀ CÁC HÀM CƠ BẢN TRONG MATLAB 2.1 GIẢN ĐỒ TRẠNG THÁI 2.1.1 Khái niệm giản đồ trạng thái Đồ thị biểu diễn mối quan hệ thông số trạng thái hệ cân nhiệt động Giản đồ trạng thái thường xây dựng sở số liệu thực nghiệm Dựa vào giản đồ người ta xác định điều kiện cân pha nhiều thông tin khác nên gọi giản đồ pha, sử dụng rộng rãi nhiệt động học phân tích hố lí Giản đồ trạng thái đơn giản hệ cấu tử biểu diễn toạ độ phẳng: áp suất - nhiệt độ (đồ thị p - T); áp suất thể tích (đồ thị p - V), thể tích - nhiệt độ (đồ thị V - T) 2.1.2 Phương trình trạng thái khí lí tưởng giản đồ pV (giản đồ Clapeyron) Phương trình xác định mối liên hệ ba thơng số trạng thái chất khí gọi phương trình trạng thái khí lí tưởng Để lập phương trình xét, ta lượng khí từ trạng thái (p1, V1, T1) sang trạng thái (p2, V2, T2) qua trạng thái trung gian 1' (p', V2, T1) đẳng trình học trước Ta dễ dàng chứng minh được: p1V1 T1 = p2V2 T2 hay Luan- 6van pV T = số Phương trình nhà vật lí người Pháp Cla-pê-rôn (Clapeyron, 1799 1864) đưa năm 1834 gọi phương trình trạng thái khí lí tưởng hay phương trình Cla-pê-rơn phát biểu: “Đối với khối lượng khí cho, tích áp suất thể tích chia cho nhiệt độ tuyệt đối đại lượng khơng đổi” Hình 2.1: Giản đồ pV 2.2 CÁC HÀM CƠ BẢN MATLAB ĐƯỢC SỬ DỤNG 2.2.1 Tổng quan MATLAB MATLAB (Matrix Laboratory) phần mềm khoa học thiết kế để cung cấp việc tính tốn số hiển thị đồ họa ngơn ngữ lập trình cấp cao MATLAB cung cấp tính tương tác tuyệt vời cho phép người sử dụng thao tác liệu linh hoạt dạng mảng ma trận để tính tốn quan sát Các liệu vào MATLAB nhập từ "Command line" từ "mfiles", tập lệnh cho trước MATLAB MATLAB cung cấp cho người dùng toolbox tiêu chuẩn tùy chọn Người dùng tạo hộp cơng cụ riêng gồm "mfiles" viết cho ứng dụng cụ thể Luan- 7van 2.2.2 Các hàm matlab sử dụng tốn 2.2.2.1 Lệnh FPRINTF Cơng dụng lệnh FPrint đưa liệu có định dạng file hình Cú pháp: fprintf(FILE *f, const char *format, varexpr) Giải thích: + Format: chuỗi định dạng (giống với định dạng hàm printf()) + Varexpr: danh sách biểu thức, biểu thức cách dấu phẩy (,) + Ví dụ : fprintf('Giai doan #%d \n',iPoint-1) ( Có %d hiển thị kết số ngun Có \n in xong xuống hàng, đưa dấu nhắc lệnh đầu dòng kế tiếp) Kết hiển thị : Giai doan #1 2.2.2.2 Lệnh SPRINTF( Tương tự lệnh FPRINTF) 2.2.2.3 Lệnh MENU Công dụng nhằm tạo menu để chọn chức Cú pháp: tên biến = menu (‘Tên menu’,‘chức 1’,‘chức 2’, … , ‘chức n’) Giải thích: - Tên menu: tiêu đề menu - Tên biến: nơi cất giá trị nhận sau chọn chức menu Chức 1, 2, …., n: chọn chức tên biến có giá trị số thứ tự chức Vd: PathType = menu(sprintf('Giai doan %d: Chon qua trinh ke tiep',iPoint-1), 'Dang ap (Constant P)', 'Dang tich (Constant V)', 'Dang nhiet (Select new V)','QUIT'); Kết hiển thị : 2.2.2.3 Lệnh WHILE Hình 2.2 Cú pháp: while Biểu thức điều kiện Các lệnh thực thi End Luan- 8van Các lệnh thực biểu thức điều kiện 2.2.2.5 Lệnh IF Cú pháp: if Biểu thức điều kiện Lệnh thực thi Đk elseif Biểu thức điều kiện Lệnh thực thi Đk end 2.2.2.6 Lệnh FOR Cú pháp: for var = start : step : end Lệnh thực thi End 2.2.2.7 Lệnh PLOT Công dụng lệnh Plot: Phần lớn câu lệnh để vẽ đồ thị mặt phẳng lệnh plot Lệnh plot vẽ đồ thị mảng liệu hệ trục thích hợp nối điểm đường thẳng Cú pháp: plot(x,y) plot(x,y,’linetype’) Giải thích: x,y: vẽ giá trị x theo giá trị y Ví dụ: plot(V,P,'o',VPlot,PPlot,'-'); Giải thích: dấu'o': dạng dấu chấm trịn Nét vẽ'-': nét liền thẳng Kết hiển thị: Hình 2.3 Luan- 9van 2.2.2.8 Lệnh AXIS Công dụng lệnh Axis chia lại trục tọa độ với cú pháp: axis([xmin xmax ymin ymax zmin zmax]) Giải thích: - xmin, ymin, zmin: giá trị nhỏ trục x, y, z - xmax, ymax, zmax: giá trị lớn trục x, y, z - Vd: axis([0 1.5*max(V) 1.5*max(P)]); Ví dụ giải thích sau: + 1.5*max(V): Gía trị nhỏ trục V 0, lớn 1.5 lần giá trị lớn giá trị V nhâp vào + 1.5*max(P): Gía trị nhỏ trục P 0, lớn 1.5 lần giá trị lớn giá trị P nhâp vào 2.2.2.9 Lệnh XLABEL, YLABEL Công dụng đặt tên cho trục toạ độ Cú pháp: xlabel(' viết tiêu đề trục’) Vd: xlabel('Volume (m^3)') ylabel('Pressure (Pa)') Giải thích: Đặt tên cho trục x Volume(m3) Đặt tên cho trục y Pressure(Pa) Kết hiển thị: Hình 2.3 2.2.2.10 Lệnh TEXT Cơng dụng: văn vị trí trục đồ thị Cú pháp: text(x,y,’chuỗi’) Giải thích: x,y: toạ độ ‘ ‘ văn cần ghi toạ độ Vd: text(V(i),P(i),sprintf(' %d',i)) Giải thích: V(i),P(i): toạ độ V,P sprintf(' %d',i): ghi số thứ tự trình 2.2.2.11 Lệnh TITLE Công dụng hàm Title đặt tiêu đề cho đồ thị Cú pháp: Title(‘text 1’,’text 2’, ) Giải thích: text 1,2 tiêu đề Luan- 10van Vd: title(sprintf('Work; Last leg = %.3f J, Total = %.3f J',W,WTotal) ); Với - sprintf('Work; Last leg = %.3f J, Total = %.3f J',W,WTotal):in tên đồ thị - %.3: chuỗi định dạng hàm sprintf, hiển thị làm tròn chữ số thập phân 2.2.2.12 Lệnh DRAWNOW Qua lượt lặp, ta phải gọi drawnow để MATLAB thực hiển thị hình vẽ Cịn khơng, đợi đến tận bạn vẽ xong hết hình cập nhật trình Chương 3: THỨ TỰ GIẢI BÀI TỐN TÌM CƠNG TỪ GIẢN ĐỔ PV TRÊN MATLAB VÀ CODE HOÀN CHỈNH 3.1 THỨ TỰ GIẢI BÀI TỐN TÌM CƠNG TỪ GIẢN ĐỒ PV a/ Khởi tạo biến: nhập số mol khí, nhập áp suất ban đầu, nhập V ban đầu b/ Đặt số khí R (J/mol), đặt nhiệt độ ban đầu T1 theo phương trình trạng thái khí lí tưởng c/ Ghi áp suất thể tích ban đầu cho biểu đồ Nhập tổng công thực được, nhập số điểm ban đầu để vẽ đường cong đẳng nhiệt d/ Lập lệnh QUIT chọn vào menu với QUIT lựa chọn thứ menu, gán giá trị Pathtype vào vòng lặp (PathType ~= QuitType) e/ Chọn trình đẳng áp, đẳng tích hay đẳng nhiệt QUIT Tạo điểm Lập trình chi tiết cho menu g/ Nếu trình đẳng áp, xác định thông số P, T Nhập biểu thức tính cơng q trình đẳng áp Thêm thể tích, áp suất vào liệu đồ thị h/ Nếu trình đẳng tích, xác định thơng số P, V, T Nhập tổng cơng thực Thêm thể tích áp suất vào liệu đồ thị Luan- 11van i/ Nếu q trình đẳng tích, xác định thông số P, V, T Nhập biểu thức tính cơng thực Tính tốn giá trị P V trình đẳng nhiệt thêm vào liệu đồ thị k/ Tạo điểm để vẽ đồ thị (tạo ma trận ngang tất điểm) l/ Vẽ đồ thị thích (chia lại trục, đặt tên trục x, y…) 3.2 CODE HOÀN CHỈNH (KÈM LỜI DỊCH CODE) Dựa vào yêu cầu mục tiêu toán đặt ra, dựa kết làm được, sau đoạn code hoàn chỉnh để giúp giải tốn tìm cơng q trình cân từ giản đồ pV: % xac dinh cong cua he qua trinh can bang tu gian (P,V) clear all; %@ khoi tao cac bien nMoles = input('Nhap so mol '); P(1) = input('Nhap ap suat ban dau (Pa): '); V(1) = input('Nhap the tich ban dau (m^3): '); R = 8.314; % Hang so (J/mole) T(1) = P(1)*V(1)/(nMoles*R); % Nhiet ban dau PPlot = P(1); % Ghi ap suat ban dau cho bieu VPlot = V(1); % Ghi the tich ban dau cho bieu WTotal = 0; iPoint = 1; % Tong cong thuc hien duoc (J) % diem ban dau NCurve = 100; % So diem duoc su dung de ve duong cong dang nhiet %@Lap cho den lenh QUIT duoc chon tu menu QuitType = 4; % QUIT la lua chon thu menu PathType = 0; % Gia tri gia cua PathType vao vong lap while(PathType ~= QuitType) Luan- 12van %@ Chon cac qua trinh (Dang ap, Dang tich, dang nhiet) or QUIT iPoint = iPoint + 1; % Diem ke tiep fprintf('Giai doan #%d \n',iPoint-1); %Co \n thi in xong xuong hang, dau dau nhac lenh ve dau dong ke tiep ; %d: Ghi so nguyen PathType = menu(sprintf('Giai doan %d: Chon qua trinh ke tiep',iPoint-1), 'Dang ap (Constant P)', 'Dang tich (Constant V)', 'Dang nhiet (Select new V)','QUIT'); %@ Neu qua trinh ke tiep la dang ap (Constant P) if( PathType == ) %@ Xac dinh thong so moi:Ap suat, nhiet V(iPoint) = input('Nhap the tich moi: '); P(iPoint) = P(iPoint-1); % Ap suat moi giong ap suat cu T(iPoint) = P(iPoint)*V(iPoint)/(nMoles*R); % Nhiet moi %@ Tinh cong qua trinh dang ap W = P(iPoint)*( V(iPoint) - V(iPoint-1) ); %@ Them the tich va ap suat vao du lieu cua thi VPlot = [VPlot V(iPoint)]; % Them diem vao du lieu the tich de ve PPlot = [PPlot P(iPoint)]; % Them diem vao du lieu ap suat de ve %@ neu qua trinh ke tiep la dang tich (Constant V) elseif( PathType == ) %@ Xac dinh thong so moi:P,V,T P(iPoint) = input('Nhap ap suat moi: '); V(iPoint) = V(iPoint-1); T(iPoint) = P(iPoint)*V(iPoint)/(nMoles*R); %@ Tong cong thuc hien duoc =0 W = 0; Luan- 13van %@ Them the tich va ap suat vao du lieu cua thi VPlot = [VPlot V(iPoint)]; % Them diem vao du lieu the tich de ve PPlot = [PPlot P(iPoint)]; % Them diem vao du lieu ap suat de ve %@ Neu qua trinh ke tiep la dang nhiet (Select new V) elseif( PathType == ) %@ Xac dinh cac thong so moi:P,V,T V(iPoint) = input('Nhap the tich moi: '); T(iPoint) = T(iPoint-1); P(iPoint) = nMoles*R*T(iPoint)/V(iPoint); %@ Cong thuc hien duoc W = nMoles*R*T(iPoint)*log(V(iPoint)/V(iPoint-1)); %@ Tinh toan cac gia tri cua P V qua trinh dang nhiet va them vao %du lieu thi for i=1:NCurve (i) = V(iPoint-1) + (i-1)/(NCurve-1)*(V(iPoint)-V(iPoint-1)); PNew(i) = nMoles*R*T(iPoint)/VNew(i); end VPlot = [VPlot VNew]; PPlot = [PPlot PNew]; % Them diem vao du lieu the tich de ve % Them diem vao du lieu ap suat de ve end %@ Draw the total path on the PV diagram, adding this leg if( PathType ~= QuitType ) WTotal = WTotal + W; % Tong cong plot(V,P,'o',VPlot,PPlot,'-'); Luan- 14van ... NGUYÊN LÍ THỨ I CỦA NHIỆT ĐỘNG LỰC HỌC - CÁC KHÁI NIỆM LIÊN QUAN, CÔNG CỦA HỆ VÀ CÔNG TRONG CHU TRÌNH 1.1 KHÁI NIỆM VỀ HỆ NHIỆT ĐỘNG, TRẠNG THÁI CÂN BẮNG - QUÁ TRÌNH CÂN BẰNG 1.1.1 Hệ nhiệt động... Chương 2: GIẢN ĐỒ PV VÀ CÁC HÀM CƠ BẢN TRONG MATLAB 2.1 GIẢN ĐỒ TRẠNG THÁI 2.1.1 Khái niệm giản đồ trạng thái Đồ thị biểu diễn mối quan hệ thông số trạng thái hệ cân nhiệt động Giản đồ trạng thái... (đồ thị p - T); áp suất thể tích (đồ thị p - V), thể tích - nhiệt độ (đồ thị V - T) 2.1.2 Phương trình trạng thái khí lí tưởng giản đồ pV (giản đồ Clapeyron) Phương trình xác định mối liên hệ