Estudo de pistao, desenho e analise

clear all

clc

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% PROGRAMA PARA ANÁLISE ESTÁTICA POR EF DE VIGAS

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Professor Daniel Colombo

% João Vitor Martinez da Mata 372.103

% Aquilla Colman 372.118

% Luiz Paulo 372.122

% Gabriel Nascimento 372.169

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% ENTRADA DE DADOS

% VALORES DOS COEFICIENTES DE RIGIDEZ DAS MOLAS

%k_x = 1e3*[10 15 10];% VALORES EM N/m

%k_y = 1e3*[10 15 15 10];% VALORES EM N/m

E=200e9; % Módulo de elasticidade para o material em Gpa p/ Inox - 304

h=0.001 % Em metros

b=0.025 % Em metros

L=0.31; % 31 cm para metros

% CALCULO DO MOMENTO DE INERCIA

I=(b*(h^3))/12; % Formula para o calculo do momento de inércia para área rentângular

% MASSAS

m1=0.056 % em Kg

m2=0.018 % em Kg

m3=0.038 % em Kg

m4=0.021 % em Kg

F=m1*9,81; % Em Newton

% NÚMERO DE NÓS

nb_node = 31; % 30 cm do tamanho da régua mais 1 cm do furo onde foi colocado os pesos, para cada cm existe um nó

% MATRIZ DE CONECTIVIDADE

mat_conect=[1 2 3 4; 3 4 5 6; 5 6 7 8; 7 8 9 10; 9 10 11 12; 11 12 13 14; 13 14 15 16; 15 16 17 18; 17 18 19 20; 19 20 21 22; 21 22 23 24; 23 24 25 26; 25 26 27 28; 27 28 29 30; 29 30 31 32; 31 32 33 34; 33 34 35 36; 35 36 37 38; 37 38 39 40; 39 40 41 42; 41 42 43 44; 43 44 45 46; 45 46 47 48; 47 48 49 50; 49 50 51 52; 51 52 53 54; 53 54 55 56; 55 56 57 58; 57 58 59 60; 59 60 61 62; 61 62 63 64]

cond_cont=[1 0 ; 2 0];

% FORÇAS EXTERNAS APLICADAS NOS GDL LIVRES

forcas_aplic= [ 3 0 ; 4 0; 5 0; 6 0; 7 0; 8 0; 9 0; 10 0; 11 0; 12 0; 13 0; 14 0; 15 0; 16 0; 17 0; 18 0; 19 0; 20 0; 21 0; 22 0; 23 0; 24 0; 25 0; 26 0; 27 0; 28 0; 29 0; 30 0; 31 0; 32 0; 33 0; 34 0; 35 0; 36 0; 37 0; 38 0; 39 0; 40 0; 41 0; 42 0; 43 0; 44 0; 45 0; 46 0; 47 0; 48 0; 49 0; 50 0; 51 0; 52 0; 53 0; 54 0; 55 0; 56 0; 57 0; 58 0; 59 0; 60 0; 61 F; 62 0]; % valores em Newtons

[l1 c1]= size(cond_cont)

[l2 c2]= size(forcas_aplic)

nb_ele= 30;

% CONSTRUÇÃO DAS MATRIZES ELEMENTARES E MONTAGEM DA MATRIZ DE RIGIDEZ GLOBAL

K_global=zeros(l1+l2);

for ii=1:nb_ele;

K_elementar= (E*I/L^3)*[12 6*L -12 6*L; 6*L 4*L^2 -6*L 2*L^2; -12 -6*L 12 -6*L; 6*L 2*L^2 -6*L 4*L^2];

mat_ident=eye(l1+l2)

mat_transf=[mat_ident(mat_conect(ii,1),:); mat_ident(mat_conect(ii,2),:);mat_ident(mat_conect(ii,3),:);mat_ident(mat_conect(ii,4),:)];

K_global=K_global+mat_transf'*K_elementar*mat_transf

end

% IMPOSIÇÃO DAS CONDIÇÕES DE CONTORNO PELO MÉTODO DO PARTICIONAMENTO DA MATRIZ DE RIGIDEZ

% IDENTIFICAÇÃO DOS GDL LIVRES E GDL IMPOSTOS

gdl_livres=forcas_aplic(:,1);

gdl_impostos = cond_cont(:,1);

% CONSTRUÇÃO DAS SUBMATRIZEScle DE RIGIDEZ

K_ll=K_global(gdl_livres,gdl_livres);

K_li=K_global(gdl_livres,gdl_impostos);

K_ii=K_global(gdl_impostos,gdl_impostos);

K_il=K_global(gdl_impostos,gdl_livres);

% CONSTRUÇÃO DOS VETORES DE FORÇAS NOS GDL LIVRES E DE DESLOCAMENTOS NOS GDL IMPOSTOS

f_liv=forcas_aplic(:,2);