(Luận văn thạc sĩ) đo lường ứng suất bề mặt trụ của thép cán dùng nhiễu xạ tia x

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ BÙI LÊ QUỐC CƯỜNG ĐO LƯỜNG ỨNG SUẤT BỀ MẶT TRỤ CỦA THÉP CÁN DÙNG NHIỄU XẠ TIA X NGÀNH: CÔNG NGHỆ CHẾ TẠO MÁY - 605204 S KC 0 Tp Hồ Chí Minh, tháng 11 năm 2010 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH  LUẬN VĂN THẠC SĨ BÙI LÊ CƯỜNG QUỐC ĐO LƯỜNG ỨNG SUẤT BỀ MẶT TRỤ CỦA THÉP CÁN DÙNG NHIỄU XẠ TIA X NGÀNH: CÔNG NGHỆ CHẾ TẠO MÁY - 605204 TP Hồ Chí Minh, tháng 11 năm 2010 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH  LUẬN VĂN THẠC SĨ BÙI LÊ CƯỜNG QUỐC ĐO LƯỜNG ỨNG SUẤT BỀ MẶT TRỤ CỦA THÉP CÁN DÙNG NHIỄU XẠ TIA X NGÀNH: CÔNG NGHỆ CHẾ TẠO MÁY - 605204 Hướng dẫn khoa học: TS LÊ CHÍ CƯƠNG TP Hồ Chí Minh, tháng 11 năm 2010 LÝ LỊCH CÁ NHÂN -I LÝ LỊCH SƠ LƢỢC: Họ & tên: Bùi Lê Cƣờng Quốc Giới tính: Nam Ngày, tháng, năm sinh: 15-11-1981 Nơi sinh: Bình Thuận Quê quán: Bình Thuận Dân tộc: Kinh Chỗ riêng địa liên lạc: Hàm Mỹ, Hàm Thuận Nam, Bình Thuận Điện thoại quan: Điện thoại nhà riêng: 0908.472.984 Fax: E-mail:Quoc_cuong1986@ymail.com II QUÁ TRÌNH ĐÀO TẠO: Đại học: Hệ đào tạo: Chính Quy Thời gian đào tạo từ 10/2000 đến7/2005 Nơi học (trƣờng, thành phố): ĐẠI HỌC SƢ PHẠM KỸ THUẬT TP.HCM Ngành học: CƠ KHÍ MÁY Tên đồ án, luận án môn thi tốt nghiệp: THIẾT KẾ VÀ TÍNH TỐN MÁY DẬP BÁNH Nơi bảo vệ đồ án tốt nghiệp: KHOA CƠ KHÍ MÁY Ngƣời hƣớng dẫn: Thạc Sĩ TRẦN QUỐC HÙNG III Q TRÌNH CƠNG TÁC CHUYÊN MÔN KỂ TỪ KHI TỐT NGHIỆP ĐẠI HỌC: Thời gian 2006 2008 Nơi công tác TRƢỜNG TRUNG HỌC PHỔ THƠNG NGUYỄN VĂN LINH TRƢỜNG TRUNG CẤP NGHỀ BÌNH THUẬN Công việc đảm nhiệm GIÁO VIÊN GIÁO VIÊN LỜI CAM ĐOAN Tơi cam đoan cơng trình nghiên cứu Các số liệu, kết nêu luận văn trung thực chƣa đƣợc công bố cơng trình khác Tp Hồ Chí Minh, ngày … tháng … năm 2010 Ký tên ghi rõ họ tên LỜI CẢM ƠN  Sau hai năm theo học chƣơng trình đào tạo sau đại học trƣờng Đại học Sƣ phạm Kỹ thuật Thành phố Hồ Chí Minh, em đúc kết đƣợc kiến thức bổ ích cho chun mơn Với đề tài nghiên cứu dƣới hình thức luận văn thạc sĩ, em vận dụng kiến thức mà đƣợc trang bị để tiến hành giải tốn thực tiễn Vì đề tài luận văn nghiên cứu giải vấn đề mẻ dựa sở tính tốn lý thuyết chun sâu lĩnh vực vật liệu dùng kỹ thuật nhiễu xạ X–quang, nên lúc đầu tiếp cận em gặp nhiều bỡ ngỡ khó khăn Nhƣng với tận tình thầy hƣớng dẫn TS Lê Chí Cƣơng, với hỗ trợ từ phía gia đình, bạn bè đồng nghiệp, luận văn em đạt đƣợc kết nhƣ mong muốn Đến đây, cho phép em gửi lời tri ân sâu sắc đến: - Ban Giám Hiệu trƣờng Đại học Sƣ phạm Kỹ thuật Thành phố Hồ Chí Minh - Thầy TS Lê Chí Cƣơng – Khoa Cơ khí máy - trƣờng Đại học Sƣ phạm Kỹ thuật Thành phố Hồ Chí Minh - Tiến sĩ NGUYỄN ĐỨC THÀNH - Giám Đốc – trung tâm hạt nhân thành phố Hồ Chí Minh - Q thầy khoa Cơ khí máy - trƣờng Đại học Sƣ phạm Kỹ thuật Thành phố Hồ Chí Minh - Phòng Đào tạo - Sau Đại học phòng khoa trƣờng Đại học Sƣ phạm Kỹ thuật Thành phố Hồ Chí Minh - Gia đình, bạn bè, đồng nghiệp anh, chị lớp cao học Cơng Nghệ Chế Tạo Máy, khóa 2008- 2010 Một lần em xin chân thành cảm ơn giúp đỡ, hỗ trợ, động viên quý báu tất ngƣời Xin trân trọng cảm ơn! Thành phố Hồ Chí Minh, tháng 11 năm 2010 Học viên thực luận văn MỤC LỤC Trang tựa Trang Quyết định giao đề tài Lý lịch cá nhân i Lời cam đoan ii Lời cám ơn iii Tóm tắt iv Mục lục v Danh sách ký hiệu vi Danh sách bảng vii Danh sách hình viii Chƣơng I DẪN NHẬP 1.1 Đặt vấn đề 1.2 Nội dung nghiên cứu 1.3 Nhiệm vụ nghiên cứu giới hạn đề tài 1.4 Phƣơng pháp nghiên cứu 1.5 Điểm luận văn 1.6 Giá trị thực tiễn luận văn Chƣơng II NGUỒN GỐC VÀ SỰ PHÁT TRIỂN TIA X 2.1 Giới thiệu tia X 2.1.1 Lịch sử tia X 2.1.2 Lịch sử phát triển tia X 2.1.3 Ứng dụng tia X 2.1.4 Tạo tia X [2] 2.1.5 Đặc điểm đƣờng xạ [3] 2.1.6 Nhiễu xạ tia X 10 a Hiện tƣợng nhiễu xạ Tia X [4] 10 b Định luật Bragg 11 2.1.7 Giới hạn bƣớc sóng tƣợng quang phổ 13 2.2 Nguyên lý cấu tạo thiết bị [1] 15 2.2.1 Phƣơng pháp chụp ảnh 15 a Cấu tạo 15 b Nguyên lý phƣơng pháp 15 2.2.2 Phƣơng pháp nhiễu xạ kế (diffractometer) 16 a Cấu tạo 16 b Nguyên lý phƣơng pháp 16 2.2.3 Chiều sâu thấm tia X 19 2.2.4 Tính ứng suất 19 2.2.5 Xác định biến dạng 19 2.2.6 Xác định mối quan hệ ứng suất - biến dạng 22 2.2.7 Ƣu điểm phƣơng pháp đo nhiễu xạ kế so với phƣơng pháp chụp ảnh 25 Chƣơng III HỆ SỐ HẤP THU TRÊN CÁC BỀ MẶT 26 3.1 Các phƣơng pháp đo máy nhiễu xạ đơn tinh thể 26 3.1.1 Phƣơng pháp đo kiểu  27 a Phƣơng pháp đo kiểu  cố định  28 b Phƣơng pháp đo kiểu  cố định o 29 3.1.2 Phƣơng pháp đo kiểu  30 a Phƣơng pháp đo kiểu  cố định  30 b Phƣơng pháp đo kiểu  cố định O 31 3.2 Hệ số hấp thụ ảnh hƣởng tới cƣờng độ nhiễu xạ mặt phẳng 31 3.3 Hệ số hấp thụ không giới hạn vùng nhiễu xạ giới hạn vùng nhiễu xạ [7] 33 3.4 Hàm số hấp thụ bề mặt trụ phƣơng pháp đo kiểu  cố định o [11] 35 3.5 Hàm số hấp thu bề mặt trụ phƣơng pháp đo kiểu  cố định  o[4]: (Lê Minh Tấn 2005) 37 3.6 Khảo sát hàm hấp thu bề mặt ellipsoid phƣơng pháp đo kiểu , cố định góc  0[5]: (Nguyễn Thị Hồng 2009) 40 Chƣơng IV THIẾT BỊ NHIỄU XẠ TIA X, CHUẨN BỊ MẪU, ĐO ĐẠC MẪU VÀ TÍNH ỨNG SUẤT 42 4.1 Thiết bị nhiễu xạ tia X X’Pert Pro 42 4.2 Chuẩn bị mẫu, đo đạc mẫu 44 4.3 Đo mẫu hệ máy nhiễu xạ X’Pert Pro 45 4.4 Tính tốn ứng suất 46 4.4.1 Tính tốn ứng suất hàm hấp thu A phẳng 46 4.4.2 Tính tốn ứng suất cho mẫu cịn lại 50 4.5 Ứng suất dùng hàm hấp thụ bề măt trụ 54 4.5.1 Tính tốn ứng suất hàm hấp thu A trụ A phẳng 54 4.6 Nguyên nhân sai số 57 Chƣơng V KẾT LUẬN VÀ HƢỚNG PHÁT TRIỂN CỦA ĐỀ TÀI 59 5.1 Tóm tắt kết đề tài 59 5.2 Đánh giá kết đề tài 61 5.3 Hƣớng phát triển đề tài 61 LUẬN VĂN THẠC SĨ Chương I  DẪN NHẬP 1.1 Đặt vấn đề: Ứng suất dƣ nhân tố quan trọng ảnh hƣởng đến độ bền, tuổi thọ chi tiết máy Ứng suất dƣ đƣợc tạo trình gia công cơ, gia công áp lực, xử lý nhiệt…là nguyên nhân gây hƣ hỏng, biến dạng vật liệu Chính việc xác định ứng suất dƣ có vai trị quan trọng q trình xử lý cải thiện điều kiện làm việc chi tiết máy Hiện nay, giới có nhiều phƣơng pháp xác định ứng suất dƣ bề mặt chi tiết nhƣ phƣơng pháp: đục lỗ, cắt tiết diện, siêu âm, nhiệt đàn hồi, nhiễu xạ Neutron, nhiễu xạ X–quang… Trong đó, phƣơng pháp nhiễu xạ X–quang có nhiều ƣu điểm: xác định xác ứng suất, dễ dàng tự động hóa… mà khơng phá hủy chi tiết mẫu Các nghiên cứu hàm hấp thu tia X Cullity, Koistinen, Ch.Genzel… áp dụng cho phƣơng pháp đo đơn giản vật liệu đẳng hƣớng Luận văn Thạc sĩ KS Lê Minh Tấn – trƣờng ĐH Sƣ Phạm Kỹ Thuật Tp HCM – năm 2008, phân tích hình dạng bề mặt trụ ảnh hƣởng đến hàm hấp thu tổng quát tia X vật liệu đẳng hƣớng Luận văn thạc sĩ KS Nguyễn Thị Hồng Trƣờng ĐH SPKT Tp HCM năm 2009 – xác định hàm hấp thu tổng quát cho bề mặt ellipsoid đo ứng suất dùng nhiễu xạ X quang Tuy nhiên, tính tốn dựa sở lý thuyết chƣa đƣa đƣợc kiểm định cách cụ thể hay giá trị đo định máy công cụ nào, dƣới giúp đỡ thầy TS Lê Chí Cƣơng, tác giả chọn lĩnh vực để làm sở nghiên cứu thực đề tài: “Đo lường ứng suất bề mặt trụ thép cán dùng nhiễu xạ tia X.” Trang Bảng 2: Dữ liệu thực nghiệm để tính tốn mẫu 12mm 2Ɵ 115.7 115.74 115.78 115.82 115.86 115.9 115.94 115.98 116.02 116.06 116.1 116.14 116.18 116.22 116.26 116.3 116.34 116.38 116.42 116.46 116.5 116.54 116.58 116.62 116.66 116.7 Ψ=00 113 112 137 145 157 158 151 164 180 187 185 195 200 190 190 200 189 180 193 180 174 180 180 170 155 159 Ψ=50 117 135 137 162 181 162 195 194 198 193 196 214 205 211 227 219 218 226 216 215 203 184 188 190 175 160 Ψ=100 137 139 133 150 154 161 190 195 187 202 205 195 207 212 208 211 209 214 199 210 197 183 186 183 183 159 Ψ=150 125 137 139 145 168 158 165 171 178 189 199 190 189 208 198 194 207 201 201 194 193 180 180 172 164 154 Ψ=200 128 134 135 159 156 154 169 170 180 176 183 180 189 182 176 175 182 178 171 164 162 163 162 159 132 141 Ψ=250 141 140 141 148 153 152 157 160 155 155 160 161 155 158 153 155 149 150 148 151 139 143 135 136 126 123 Ứng với góc ψ=00 Phương trình nội suy bậc hai y = -230.65x2 + 53700x – 308 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9202 khoảng cách d = 0.906299 Ứng với góc ψ=50 Phương trình nội suy bậc hai y = -295.23x2 + 68734x – 407 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9247 khoảng cách d = 0.906312 Ứng với góc ψ=100 Phương trình nội suy bậc hai y = -252.19x2 + 58707x – 306 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9169 khoảng cách d = 0.906378 Ứng với góc ψ=150 Phương trình nội suy bậc hai y = -232.87x2 + 54210x – 300 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9385 khoảng cách d = 0.906373 Ứng với góc ψ=200 Phương trình nội suy bậc hai y = -205.81x2 + 47903x – 301 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9306 khoảng cách d = 0.906465 Ứng với góc ψ=250 Phương trình nội suy bậc hai y = -104.72x2 + 24372x – 106 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9269 khoảng cách d = 0.90651 Bảng 3: Dữ liệu thực nghiệm để tính tốn mẫu 16mm 2Ɵ 115.7 115.74 115.78 115.82 115.86 115.9 115.94 115.98 116.02 116.06 116.1 116.14 116.18 116.22 116.26 116.3 116.34 116.38 116.42 116.46 116.5 116.54 116.58 116.62 116.66 116.7 Ψ=00 163 163 159 174 180 202 205 211 236 240 242 245 233 244 240 252 229 225 233 228 215 217 193 185 170 181 Ψ=50 148 166 164 184 211 214 218 210 235 232 235 238 222 233 239 239 228 215 222 196 198 192 186 175 151 162 Ψ=100 158 163 170 185 183 209 217 219 224 222 239 235 220 250 239 225 226 229 214 214 212 189 177 174 166 163 Ψ=150 170 173 190 190 192 204 214 208 210 212 209 216 222 223 218 204 201 203 200 198 197 184 178 156 133 127 Ψ=200 156 156 155 175 187 197 199 190 207 208 211 204 208 211 204 215 214 200 194 191 200 184 175 177 166 162 Ψ=250 145 149 161 166 170 165 171 173 178 183 180 193 186 191 188 189 191 188 192 185 189 184 175 170 170 167 Ứng với góc ψ=00 Phương trình nội suy bậc hai y = -336.98x2 + 78409x – 506 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9103 khoảng cách d = 0.9066 Ứng với góc ψ=50 Phương trình nội suy bậc hai y = -343.22x2 + 79855x – 508 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9247 khoảng cách d = 0.90668 Ứng với góc ψ=100 Phương trình nội suy bậc hai y = -328.98x2 + 76537x – 406 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9332 khoảng cách d = 0.906716 Ứng với góc ψ=150 Phương trình nội suy bậc hai y = -192.57x2 + 44797x – 306 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9194 khoảng cách d = 0.906772 Ứng với góc ψ=200 Phương trình nội suy bậc hai y = -223.44x2 + 51975x – 300 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9098 khoảng cách d = 0.906808 Ứng với góc ψ=250 Phương trình nội suy bậc hai y = -130.36x2 + 30319x – 206 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9295 khoảng cách d = 0.906892 Bảng 4: Dữ liệu thực nghiệm để tính tốn mẫu 18mm 2Ɵ 115.7 115.74 115.78 115.82 115.86 115.9 115.94 115.98 116.02 116.06 116.1 116.14 116.18 116.22 116.26 116.3 116.34 116.38 116.42 116.46 116.5 116.54 116.58 116.62 116.66 116.7 Ψ=00 150 155 162 171 186 197 200 212 241 224 241 225 256 248 250 240 238 241 257 230 225 232 241 217 196 195 Ψ=50 115 149 144 155 167 191 179 221 220 200 240 245 250 260 247 240 265 254 230 242 243 212 210 217 196 190 Ψ=100 124 159 147 177 200 185 217 234 213 251 267 260 255 267 281 266 251 252 237 254 235 233 216 190 194 168 Ψ=150 131 152 147 177 170 174 185 212 223 238 230 240 230 250 260 262 250 240 250 230 228 245 240 215 205 193 Ψ=200 147 165 190 174 187 199 200 220 241 240 233 240 243 233 240 237 240 223 239 222 230 230 210 194 186 170 Ψ=250 130 142 145 168 173 182 183 189 208 199 206 208 209 211 219 230 224 227 229 225 218 220 206 205 209 190 Ứng với góc ψ=00 Phương trình nội suy bậc hai y = -300.29x2 + 69834x – 406 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9181 khoảng cách d = 0.906952 Ứng với góc ψ=50 Phương trình nội suy bậc hai y = -377.59x2 + 87813x – 506 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9218 khoảng cách d = 0.906936 Ứng với góc ψ=100 Phương trình nội suy bậc hai y = -439.92x2 + 102295x – 606 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9282 khoảng cách d = 0.907013 Ứng với góc ψ=150 Phương trình nội suy bậc hai y = -340.54x2 + 79182x – 506 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9192 khoảng cách d = 0.907043 Ứng với góc ψ=200 Phương trình nội suy bậc hai y = -304.18x2 + 70722x – 406 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9273 khoảng cách d = 0.907086 Ứng với góc ψ=250 Phương trình nội suy bậc hai y = -215.16x2 + 50012x – 306 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.915 khoảng cách d = 0.907228 Bảng 5: Dữ liệu thực nghiệm để tính tốn mẫu 20mm 2Ɵ 115.7 115.74 115.78 115.82 115.86 115.9 115.94 115.98 116.02 116.06 116.1 116.14 116.18 116.22 116.26 116.3 116.34 116.38 116.42 116.46 116.5 116.54 116.58 116.62 116.66 116.7 Ψ=00 144 152 178 180 203 211 233 237 251 220 240 230 260 243 259 261 258 237 239 239 218 210 209 188 184 172 Ψ=50 153 194 179 201 242 248 230 255 260 270 274 254 245 260 259 227 251 218 229 234 194 197 175 181 152 154 Ψ=100 182 176 183 191 207 225 242 258 255 240 240 275 252 268 270 273 245 264 239 239 237 237 190 180 164 162 Ψ=150 153 154 165 177 168 192 218 200 232 220 240 230 250 221 251 241 258 251 255 228 240 227 219 205 186 175 Ψ=200 164 158 195 180 184 193 215 205 215 220 228 229 228 226 206 208 216 209 200 180 180 184 166 153 152 152 Ψ=250 146 139 154 161 159 167 169 169 192 191 184 189 201 195 215 193 206 195 207 201 193 190 181 190 185 172 Ứng với góc ψ=00 Phương trình nội suy bậc hai y = -386.08x2 + 89760x – 506 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9245 khoảng cách d = 0.907229 Ứng với góc ψ=50 Phương trình nội suy bậc hai y = -387.26x2 + 90008x – 500 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9147 khoảng cách d = 0.90727 Ứng với góc ψ=100 Phương trình nội suy bậc hai y = -372.8x2 + 86643x – 480 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9197 khoảng cách d = 0.907301 Ứng với góc ψ=150 Phương trình nội suy bậc hai y = -307.89x2 + 71549x – 460 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9124 khoảng cách d = 0.907368 Ứng với góc ψ=200 Phương trình nội suy bậc hai y = -250.92x2 + 58306x – 306 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9204 khoảng cách d = 0.907408 Ứng với góc ψ=250 Phương trình nội suy bậc hai y = -171.94x2 + 39944x – 206 Tỷ số tương quan hàm nội suy bậc hai: χ = 0.9029 khoảng cách d = 0.907546  Chương trình tính ứng suất tìm hàm hấp thụ, ứng suất trụ góc 2Ɵ =1100÷1200 viết phần mềm Matlab: clc clear all format long; %nhap gia tri sìnsy sin2sy(1)=0; sin2sy(2)=0.0075; sin2sy(3)=0.0301; sin2sy(4)=0.0669; sin2sy(5)=0.1169; sin2sy(6)=0.1786; %nhap gia tri cua dinh D d(1)=0.90625; d(2)=0.90622; d(3)=0.906263; d(4)=0.906316; d(5)=0.906326; d(6)=0.90645; Sxx=0; for i=1:length(sin2sy) Sxx=sin2sy(i)*sin2sy(i)+ Sxx end Sx=0; for i=1:length(sin2sy) Sx=sin2sy(i)+ Sx end Sxy=0; for i=1:length(sin2sy) Sxy=sin2sy(i)*d(i)+ Sxy end Sy=0; for i=1:length(sin2sy) Sy=d(i)+ Sy end D=Sx*Sx - length(sin2sy)*Sxx a=(1/D)*(Sx*Sy-length(sin2sy)*Sxy) b=(1/D)*(Sx*Sxy - Sxx*Sy) x=0:0.01:0.2; for i=1:length(x) y(i)=a*x(i)+b; end plot(x,y,'>-','LineWidth',2); legend('y=0.0011x+0.9062'); ylabel('d'); grid on hold on plot(sin2sy,d,'hr','LineWidth',2) clear a b Sx Sxx Sy Sxy x y d %================================== d(1)=0.906299; d(2)=0.906312; d(3)=0.906378; d(4)=0.906373; d(5)=0.906465; d(6)=0.90651; Sxx=0; for i=1:length(sin2sy) Sxx=sin2sy(i)*sin2sy(i)+ Sxx end Sx=0; for i=1:length(sin2sy) Sx=sin2sy(i)+ Sx end Sxy=0; for i=1:length(sin2sy) Sxy=sin2sy(i)*d(i)+ Sxy end Sy=0; for i=1:length(sin2sy) Sy=d(i)+ Sy end D=Sx*Sx - length(sin2sy)*Sxx a=(1/D)*(Sx*Sy-length(sin2sy)*Sxy) b=(1/D)*(Sx*Sxy - Sxx*Sy) x=0:0.01:0.2; for i=1:length(x) y(i)=a*x(i)+b; end plot(x,y,'d-g','LineWidth',2); ylabel('d'); legend('y=0.0012x+0.9063'); grid on hold on plot(sin2sy,d,'sm','LineWidth',2) % tru 16 clear a b Sx Sxx Sy Sxy x y d %================================== d(1)=0.9066; d(2)=0.90668; d(3)=0.906716; d(4)=0.906772; d(5)=0.906808; d(6)=0.906892; Sxx=0; for i=1:length(sin2sy) Sxx=sin2sy(i)*sin2sy(i)+ Sxx end Sx=0; for i=1:length(sin2sy) Sx=sin2sy(i)+ Sx end Sxy=0; for i=1:length(sin2sy) Sxy=sin2sy(i)*d(i)+ Sxy end Sy=0; for i=1:length(sin2sy) Sy=d(i)+ Sy end D=Sx*Sx - length(sin2sy)*Sxx a=(1/D)*(Sx*Sy-length(sin2sy)*Sxy) b=(1/D)*(Sx*Sxy - Sxx*Sy) x=0:0.01:0.2; for i=1:length(x) y(i)=a*x(i)+b; end plot(x,y,' rs','LineWidth',2); ylabel('d'); legend('y=0.0014x+0.9067'); grid on hold on plot(sin2sy,d,'pm','LineWidth',2) % axis([0 0.2 0.9066 0.9070]) % tru 18 clear a b Sx Sxx Sy Sxy x y d %================================== d(1)=0.906952; d(2)=0.906936; d(3)=0.907013; d(4)=0.907043; d(5)=0.907086; d(6)=0.907228; Sxx=0; for i=1:length(sin2sy) Sxx=sin2sy(i)*sin2sy(i)+ Sxx end Sx=0; for i=1:length(sin2sy) Sx=sin2sy(i)+ Sx end Sxy=0; for i=1:length(sin2sy) Sxy=sin2sy(i)*d(i)+ Sxy end Sy=0; for i=1:length(sin2sy) Sy=d(i)+ Sy end D=Sx*Sx - length(sin2sy)*Sxx a=(1/D)*(Sx*Sy-length(sin2sy)*Sxy) b=(1/D)*(Sx*Sxy - Sxx*Sy) x=0:0.01:0.2; for i=1:length(x) y(i)=a*x(i)+b; end plot(x,y,'p-m','LineWidth',2); ylabel('d'); % -Chương trình tính ứng suất trụ clc clear all format short legend('y=0.0015x+0.9069'); grid on hold on plot(sin2sy,d,'bx','LineWidth',3) % axis([0 0.2 0.9069 0.9073]) % tru 20 clear a b Sx Sxx Sy Sxy x y d %================================== d(1)=0.907229; d(2)=0.90727; d(3)=0.907301; d(4)=0.907368; d(5)=0.907408; d(6)=0.907546; Sxx=0; for i=1:length(sin2sy) Sxx=sin2sy(i)*sin2sy(i)+ Sxx end Sx=0; for i=1:length(sin2sy) Sx=sin2sy(i)+ Sx end Sxy=0; for i=1:length(sin2sy) Sxy=sin2sy(i)*d(i)+ Sxy end Sy=0; for i=1:length(sin2sy) Sy=d(i)+ Sy end D=Sx*Sx - length(sin2sy)*Sxx a=(1/D)*(Sx*Sy-length(sin2sy)*Sxy) b=(1/D)*(Sx*Sxy - Sxx*Sy) x=0:0.01:0.2; for i=1:length(x) y(i)=a*x(i)+b; end plot(x,y,'-ok','LineWidth',2); ylabel('d'); legend('y=0.0016x+0.9072'); legend('mau 10','mau 12','mau 16' ,'mau 18','mau 20'); grid on hold on plot(sin2sy,d,'hm','LineWidth',2) close all sin2sy(1)=0; sin2sy(2)=0.0075; sin2sy(3)=0.0301; sin2sy(4)=0.0669; sin2sy(5)=0.1169; sin2sy(6)=0.1786; % %tru 10 d(1)=0.90625; d(2)=0.90622; d(3)=0.906263; d(4)=0.906316; d(5)=0.906326; d(6)=0.90645; %tru 12 d(1)=0.906299; d(2)=0.906312; d(3)=0.906378; d(4)=0.906373; d(5)=0.906465; d(6)=0.90651; %tru 16 d(1)=0.9066; d(2)=0.90668; d(3)=0.906716; d(4)=0.906772; d(5)=0.906808; d(6)=0.906892; %tru 18 d(1)=0.906952; d(2)=0.906936; d(3)=0.907013; d(4)=0.907043; d(5)=0.907086; d(6)=0.907228; % tru 20 d(1)=0.907229; d(2)=0.90727; d(3)=0.907301; d(4)=0.907368; d(5)=0.907408; d(6)=0.907546; Sxx=0; for i=1:length(sin2sy) Sxx=sin2sy(i)*sin2sy(i)+ Sxx; end Sx=0; for i=1:length(sin2sy); Sx=sin2sy(i)+ Sx; end Sxy=0; for i=1:length(sin2sy); Sxy=sin2sy(i)*d(i)+ Sxy; end Sy=0; for i=1:length(sin2sy); Sy=d(i)+ Sy; end D=Sx*Sx - length(sin2sy)*Sxx; a=(1/D)*(Sx*Sy-length(sin2sy)*Sxy); b=(1/D)*(Sx*Sxy - Sxx*Sy); x=0:0.01:0.2; for i=1:length(x) y(i)=a*x(i)+b; end S_phang=1.625*a/d(1) %======================= % omega=1; teta_do=110:2:120; %112:1:120; beta_do=60:5:85; teta=pi*teta_do/180; beta=pi*beta_do/180; sy_do=0:5:25; gama=teta+beta; singama=sin(gama) R=1:0.2:2; muy=2490.6; for i=1:length(R) for j=1:length(R) A(i,j)=2*muy*R(i)-singama(j) end end % %======================== % teta_do=51:1:60; %112:1:120; sy_do=0:5:25; sy=pi*sy_do/180; teta=2*pi*teta_do/180; for i=1:length(sy) Ap(i)=1-(cos(teta)/sin(teta))*(sin(sy(i))/cos(sy(i))) end for i=1:length(Ap) for j=1:length(Ap) S_tru(i,j)= S_phang*A(i,j)/Ap(j) end end hold on grid on plot(teta_do,S_tru(1,:),'o-g','LineWidth',2.5) plot(teta_do,S_tru(2,:),'*-k','LineWidth',2.5) plot(teta_do,S_tru(4,:),'d-r','LineWidth',2.5) plot(teta_do,S_tru(5,:),'p-m','LineWidth',2.5) plot(teta_do,S_tru(6,:),'s-c','LineWidth',2.5) legend('R=10','R=12','R=16','R=18','R=20') clear all clc close all %%%%%%%%%%%%%%%%%%%%%%% % UNG SUAT % %%%%%%%%%%%%%%%%%%%%%% x=xlsread('a Quoc.xls','trucu-0', 'B96:B121'); y=xlsread('a Quoc.xls','trucu-0', 'C96:C121'); [d0,z]=Dspacing(x,y); figure (1) %duongnhieuxa(x,y) plot(x,y,'ok','markersize',5,'linewidth',2) axis on hold on grid on [d,z]=Dspacing(x,y); y=ptduongnhieuxa(z); x=[115.7:0.04:116.7]; y=subs(y,x); xlabel('theta') ylabel('count') plot(x,y,'linewidth',2) hold on %%%%%%%%%%%%%%%%%%%% x=xlsread('a Quoc.xls','trucu-5', 'B92:B121'); y=xlsread('a Quoc.xls','trucu-5', 'C92:C121'); [d5,z]=Dspacing(x,y); figure (2) plot(x,y,'ok','markersize',5,'linewidth',2) axis on hold on grid on [d,z]=Dspacing(x,y); y=ptduongnhieuxa(z); x=[115.54:0.04:116.7]; y=subs(y,x); xlabel('theta') ylabel('count') plot(x,y,'linewidth',2) hold on %duongnhieuxa(x,y) %%%%%%%%%%%%%%%%%%%%% x=xlsread('a Quoc.xls','trucu-10', 'B94:B124'); y=xlsread('a Quoc.xls','trucu-10', 'C94:C124'); [d10,z]=Dspacing(x,y); figure (3) plot(x,y,'ok','markersize',5,'linewidth',2) axis on hold on grid on [d,z]=Dspacing(x,y); y=ptduongnhieuxa(z); x=[115.62:0.04:116.82]; y=subs(y,x); xlabel('theta') ylabel('count') plot(x,y,'linewidth',2) hold on %duongnhieuxa(x,y) %%%%%%%%%%%%%%%%%%%%% x=xlsread('a Quoc.xls','trucu-15', 'B94:B124'); y=xlsread('a Quoc.xls','trucu-15', 'C94:C124'); d15=Dspacing(x,y); [d15,z]=Dspacing(x,y); figure (4) plot(x,y,'ok','markersize',5,'linewidth',2) axis on hold on grid on [d,z]=Dspacing(x,y); y=ptduongnhieuxa(z); x=[115.62:0.04:116.82]; y=subs(y,x); xlabel('theta') ylabel('count') plot(x,y,'linewidth',2) hold on %duongnhieuxa(x,y) %%%%%%%%%%%%%%%%%%%%% x=xlsread('a Quoc.xls','trucu-20', 'B93:B123'); y=xlsread('a Quoc.xls','trucu-20', 'C93:C123'); d20=Dspacing(x,y); [d20,z]=Dspacing(x,y); figure (5) plot(x,y,'ok','markersize',5,'linewidth',2) axis on hold on grid on [d,z]=Dspacing(x,y); y=ptduongnhieuxa(z); x=[115.58:0.04:116.78]; y=subs(y,x); xlabel('theta') ylabel('count') plot(x,y,'linewidth',2) hold on % %duongnhieuxa(x,y) % % % %%%%%%%%%%%%%%%%%%%% x=xlsread('a Quoc.xls','trucu-25', 'B95:B125'); y=xlsread('a Quoc.xls','trucu-25', 'C95:C125'); d20=Dspacing(x,y); [d25,z]=Dspacing(x,y); figure (6) % duongnhieuxa(x,y) plot(x,y,'ok','markersize',5,'linewidth',2) %axis on grid on hold on [d,z]=Dspacing(x,y); y=ptduongnhieuxa(z); x=[115.66:0.04:116.86]; y=subs(y,x); plot(x,y,'linewidth',2) xlabel('theta') ylabel('count') hold on ... ? ?Đo lường ứng suất bề mặt trụ thép cán dùng nhiễu x? ?? tia X. ” Trang LUẬN VĂN THẠC SĨ Nội dung nghiên cứu: Nghiên cứu sở đo lƣờng ứng suất dùng tia X Phân tích tính tốn ứng suất bề mặt trụ thép. .. mẫu lên bàn gá máy đo nhiễu x? ?? đơn tinh thể, điều chỉnh mặt phẳng đo ứng suất trùng với mặt phẳng chứa tia X tới tia nhiễu x? ??, đƣa máy vị trí chuẩn bị đo (tia X tới tia X nhiễu x? ?? gần sát tức góc... Trang 41 z2 Rb2 LUẬN VĂN THẠC SĨ Chương IV  THIẾT BỊ NHIỄU X? ?? TIA X, CHUẨN BỊ MẪU, ĐO ĐẠC MẪU VÀ TÍNH ỨNG SUẤT 4.1 Thiết bị nhiễu x? ?? tia X X’Pert Pro: Hệ máy nhiễu x? ?? tia X X’Pert Pro hãng