LÝ Thuyết đàn hồi Báo cáo Matlap

TRƯỜNG ĐẠI HỌC BÁCH KHOA TP.HCM KHOA KHOA HỌC ỨNG DỤNG BÁO CÁO BÀI TẬP LỚN LÝ THUYẾT ĐÀN HỒI ĐỀ TÀI: GIẢI BÀI TOÁN LÝ THUYẾT BIẾN DẠNG BẰNG MATLAP Giáo viên hướng dẫn: Lớp: Thực hiện: Sinh viên thực hiện: Vũ Cơng Hòa CKT16 Nhóm Lài Anh Tuấn Lê Minh Thuận Nguyễn Bá Hưng Chiếng Quang Phước 1613882 1613442 1611431 1612724 Viết code Năm học: 2017-2018 Ngày 25 tháng năm 2018 BÀI TẬP LỚN LÝ THUYẾT ĐÀN HỒI – NHÓM Nội dung Đề tài: Giải toán liên quan đến Lý Thuyết Biến Dạng Ý tưởng nhóm: Dùng matlap để tính tốn biến dạng chính,Tenso biến dạng lệch,vẽ đường tròn Morh,… Thiết kế: Matlap Mục lục Mục lục …………………………………………………………………… I.CODE ………………………………………………………………… II.VÍ DỤ MINH HỌA …………………………………………………… BÀI TẬP LỚN LÝ THUYẾT ĐÀN HỒI – NHÓM 2 I.CODE function BTLltdhnhom2 T=input('nhap ma tran tenso bien dang T= '); % kiem tra ma tran T co la ma tran doi xung hay khong !! for i=1:size(T,1) for j = 1:size(T,2) while T(i,j) ~= T(j,i) T= input(' xin loi day khong phai la ma tran doi xung!! \n vui long nhap lai ma tra T= '); break end end end j1=T(1,1)+T(2,2)+T(3,3); j2=T(1,1)*T(2,2)+T(2,2)*T(3,3)+T(3,3)*T(1,1)-T(1,2)*T(2,1)T(2,3)*T(3,2)-T(1,3)*T(3,1); j3=det(T); %giai phuong trinh bien dang p=[1 -j1 j2 -j3]; a=roots(p); % dung vong lap de tao lai a=[e1;e2;e3] voi e1xt a4=[xt:0.0000001:k2]; b4=(yt/(xt-k2))*a4-k2*yt/(xt-k2); plot(a4,b4,'blue') else a4=[k2:0.0000001:xt]; b4=(yt/(xt-k2))*a4-k2*yt/(xt-k2); plot(a4,b4,'blue') end a5=[e1:0.000001:xt]; b5=(yt/(xt-e1))*a5-e1*yt/(xt-e1); plot(a5,b5,'LineWidth',2) %%%%%%%%% r5=sqrt((e4-k3)^2+(y/2)^2); syms X mm=solve(r1^2-(X-k1)^2==r5^2-(X-k3)^2); xt1=double(mm); a3=[xt1:0.00001:e4]; b3=sqrt(r5^2-(a3-k3).^2); plot(a3,b3,'LineWidth',2) yt1=double(sqrt(r5^2-(xt1-k3)^2)); if k3