[Code] Matlab Tính các giá trị thủy tĩnh của Tàu

Xem chủ đề cũ hơn Xem chủ đề mới hơn Go down

default [Code] Matlab Tính các giá trị thủy tĩnh của Tàu

Bài gửi by tnvule_ship on Mon Jan 02, 2012 11:16 am

clc
clear all
p= input('nhap gamma : ') ;

r= input('nhap chieu rong tau : ');

g = input(' nhap vao so duong nuoc: ');

for i=1:g
fprintf(' nhap vao chieu cao duong nuoc DN %d: ',i);
h(i) = input('');
end;
for i=1:g
fprintf('nhap vector x%d cua duong nuoc: ',i);
x = input('');
fprintf('nhap vector y%d cua duong nuoc: ',i);
y = input('');
giongnhau = 1;
for (k=1:length(y)-1)
if (y(k) ~= y(k+1))
giongnhau=0;
end;
end;
l(i) = length(x);
dau = x(1);
cuoi = x(length(x));
mx(i,1:l(i)) = x;
my(i,1:l(i)) = y;
if (giongnhau == 0)
xx = linspace(dau,cuoi,101);
yy = spline(x,y,xx);
else
xx = linspace(dau,cuoi,101);
yy = linspace(y(1),y(1),101);
end;

f(i) = simp1(100,xx,yy);
Aw(i) = f(i)*2;
Moy(i) = Momenttinh(100,xx,yy);
a(i) = Moy(i)/f(i);
IL(i) = Momentquantinh(100,xx,yy);
IT(i) = Mqtn(100,xx,yy);
L(i)=cuoi-dau;
Trim(i) = (p*IL(i))/L(i);
ILa(i)=IL(i)-a(i)*a(i)*f(i);
Cw(i)=f(i)/(L(i)*r);
end;
AwT = Aw
LCF = a
ILT = IL
ITT = IT
LT = L
TrimT = Trim
ILaT = ILa
giongnhau =1;
for (k=1:length(Aw)-1)
if (y(k) ~= y(k+1))
giongnhau=0;
end;
end;

for i=1:g
if (giongnhau ==0) %khong giong nhau
x=AwT;
y= h;
a=0;
b=h(i);
yy=linspace(a,b,101);
xx=spline(y,x,yy);
n = 100;
V(i)=simp1(100,yy,xx);
else %giong nhau
V(i) = Aw(i) * h(i);
end;

end;

VC = V

for i=1:g
Cm(i) = Aw(i)/(L(i)*r);
CmT = Cm

Cb(i)=V(i)/(L(i)*r*h(i));
CbT=Cb

CP(i) = Cb(i)/Cm(i);
CPT=CP

BM(i) = IT(i)/V(i);
BMT=BM

BML(i) = ILa(i)/V(i);
BMLT = BML

Disp(i)=V(i)*p;
DispT = Disp
end;

tnvule_ship
New Member
New Member

Tổng số bài gửi : 1
Age : 26
Điểm : 3
Uy tín : 0
Registration date : 02/01/2012

Xem lý lịch thành viên

Về Đầu Trang Go down

default Re: [Code] Matlab Tính các giá trị thủy tĩnh của Tàu

Bài gửi by shipdesign on Thu Jan 05, 2012 11:37 am

Rất tốt, cần phát huy hơn nữa những module thế này để hoàn thiện một chương trình hoàn thiện hay chí ít phục vụ cho một tác vụ nào đó.

Mình cũng có một bài toán rất hay phát sinh từ thực tế công việc.

Bài toán thực tế là: Tính toán cẩu cho một sàn công tác có khối lượng 100 Tấn, dự kiến cẩu bằng 4 sợi cáp có chiều dài khác nhau, cho rằng trọng tâm sàn trùng với hình chiếu của móc cẩu trên mặt phẳng nằm ngang. Vị trí 4 tai cẩu trên sàn hợp thành 1 tứ giác bất kỳ. Chiều dài các cạnh(tính từ tai cẩu này đến tai cẩu kia) cho trước(đo được trên thực tế). Chiều dài 4 sợi cáp cẩu cho trước. Tính lực căng trên mỗi sợi dây.

Để giải bài toán này ta cẩn phải giải 2 bài toán hóc búa:
1. Tính lực thưo phương thẳng đứng tác dụng lên mỗi tai cẩu. Do đây là bài toán lực không gian nên khá phức tạp. Qua tìm hiểu, chỉ có thể giải bằng phương pháp gần đúng là quy lực phân bố sàn về dạng lực phân bố theo đường, xem sàn là 1 dầm chịu tác dụng bởi 3 lực F1, , F2(lực tổng cộng của 2 dây cáp cùng phía) và G(trọng lượng sàn). Tính được F1, F2 ta lại giải tiếp bài toán lực trên 2 gối tựa để tính lực thực tế trên các móc(ví dụ lực F1 sẽ phân thành lực P1 và P2, coi là bài toán dầm phẳng chịu tác dụng của 3 lực F1, P1, P2 với tỉ lệ khoảng cách theo như khoảng cách từ tai cẩu đến tâm chiếu xuống phương của 2 móc cẩu, tương tự tính được F3 = P3+P4)
2. Tính đường cao của hình chóp SABCD, có AB=a, BC=b, CD=c, DA=d, SA=l1, SB=l2, SC=l3,SD=l4. khi có chiều dài đường cao, ta có thể tính được góc nghiêng của cáp nâng qua quan hệ tam giác lượng, từ đó tính được lực căng dây sau khi đã biết lực tác dụng theo phương thẳng đứng trên mỗi tai cẩu (P1, P2, P3, P4).

Ở bài toán thứ 2 nảy sinh ra một hệ phương trình 28 phương trình 28 ẩn chứa các hàm lượng giác. Để giải được bài này phải dùng đến MATLAB. Nhưng vấn đề là nếu lập ma trận để giải theo các phương pháp khử(Gauss...)thì không giải quyết được hàm lượng giác. Trong MATLAB mình cũng tìm thấy hàm Slove nhưng hình như hàm này chạy không ổn định. Đã thử với hệ 5 pt 5 ẩn kết quả chương trình báo lỗi.

Nếu bạn nào sành về MATLAB thì nên thử bài toán thực tế này. Xin được hồi âm.

shipdesign
Admin
Admin

Tổng số bài gửi : 174
Age : 35
Đến từ: : shipdesign company
Nghề nghiệp/sở thióch : research
Điểm : 340
Uy tín : 3
Registration date : 24/05/2008

Xem lý lịch thành viên http://www.esnips.com/user/tuyen

Về Đầu Trang Go down

default Re: [Code] Matlab Tính các giá trị thủy tĩnh của Tàu

Bài gửi by hoangtukimcuong_tvn on Wed Oct 24, 2012 10:29 am

Chào các anh em. Các bậc tiền bối trong ngành đóng tàu!
có bài ví dụ nào cho đoạn code trên không các anh em. tôi mới học về matlab nên không có hiểu nhiều về phần này mong các anh em giúp đỡ chỉ thêm cho. cảm ở anh em rất nhiều.
Chúc anh em và các bậc tiền bối có một ngày vui.


hoangtukimcuong_tvn
New Member
New Member

Tổng số bài gửi : 1
Age : 27
Điểm : 1
Uy tín : 0
Registration date : 24/10/2012

Xem lý lịch thành viên

Về Đầu Trang Go down

default Re: [Code] Matlab Tính các giá trị thủy tĩnh của Tàu

Bài gửi by Sponsored content


Sponsored content


Về Đầu Trang Go down

Xem chủ đề cũ hơn Xem chủ đề mới hơn Về Đầu Trang

- Similar topics

 
Permissions in this forum:
Bạn không có quyền trả lời bài viết