پروژه شبیه سازی سازه بین ۳ تا ۶ درجه آزادی

پروژه شبیه سازی سازه بین ۳ تا ۶ درجه آزادی

توضیحات پروژه: گروه متلب پروژه در ادامه پروژه های آماده پروژه شبیه سازی سازه بین ۳ تا ۶ درجه آزادی که توسط اساتید متلب پروژه تهیه شده است را برای فروش در اختیار کاربران قرار داده است این پروژه از دروس رشته مهندسی هوافضا می باشد.

clc

clear

close all

 ماتریس جرم سختی وارد می شود.

ماتریس ها باید مربعی باشد.

 

;M=input(‘Enter the mass matrix: ‘)

;[n,o]=size(M)

if n~=o

;error(‘M matrix must be square!’)

end

;K=input(‘Enter the sti§ness matrix: ‘)

;[n,o]=size(K)

if n~=o

;;error(‘K matrix must be square!’)

end

% Giving the option to input damping via

% modal damping ratios or through the multiplying

% factors Cm and Ck.

%

بر اساس میزان دمپ مودال یا ضرایب معادله حل شود.

1 مودال دمپ و گرنه استفاه از ضرایب

fprintf(‘Press 1 to enter modal damping ratios, ‘)

;qz=input(‘or anything else to enter Cm and Ck. ‘)

برای هر مود میزان دمپ وارد می شود

if qz==1

for iz=1:n

;zeta(iz)=input([‘Enter the damping for mode ‘,num2str(iz),’: ‘])

end

else

;fprintf(‘Given that [c]=Cm[M]+Ck[K], ‘)

;Cm=input(‘Enter the factor Cm: ‘)

;Ck=input(‘Enter the factor Ck: ‘)

end

;qu=0

ماتریس مقادیر و بردار ویژه

;[u,l]=eig(K,M)

% Using ìeigîin this way allows us to subtract M*w^2

% from K, instead of I*w^2 (where I is the n by n identity

% The output from ìeigîgives unit-length eigenvectors

% We need to scale them with respect to M

%

مقدا ر تغیر مکان و سرعت اولیه هر استوری وارد می شود.

for s=1:n

;alfa=sqrt(u(:,s)’*M*u(:,s))

;u(:,s)=u(:,s)/alfa

end

;x0=input(‘Enter the initial displacement column vector: ‘);

xd0=input(‘Enter the initial velocity column vector: ‘)

;tf=input(‘Enter the final time: ‘)

;p0=input(‘Enter the p0’)

;t=0:0.1:tf; q=tf/0.1;x=zeros(size(n,q))

خروجی برنامه در محیط متلب