Scribd
Cargar
160 vistas4 páginas

Programación de Elementos Finitos 3D

Este documento presenta el código para resolver un problema de elementos finitos de una estructura tridimensional de vigas (truss) utilizando el método de los elementos finitos. Se define la geometría, condiciones de contorno y cargas de la estructura, y se ensambla la matriz de rigidez global considerando las matrices de rigidez locales de cada elemento viga. Luego, se resuelve el sistema para obtener los desplazamientos nodales, y con esto los esfuerzos y fuerzas internas en cada elemento. Finalmente, se muestran los resultados obten

Cargado por

Karla Villegas
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
160 vistas4 páginas

Programación de Elementos Finitos 3D

Este documento presenta el código para resolver un problema de elementos finitos de una estructura tridimensional de vigas (truss) utilizando el método de los elementos finitos. Se define la geometría, condiciones de contorno y cargas de la estructura, y se ensambla la matriz de rigidez global considerando las matrices de rigidez locales de cada elemento viga. Luego, se resuelve el sistema para obtener los desplazamientos nodales, y con esto los esfuerzos y fuerzas internas en cada elemento. Finalmente, se muestran los resultados obten

Cargado por

Karla Villegas
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

ESCUELA POLITÉCNICA NACIONAL

FACULTAD DE INGENIERÍA MECÁNICA


ELEMENTOS FINITOS

Villegas Ponce Karla Elizabeth


Velasco Carrera Xavier Efren
Gr 2
10/01/2018

DEBER

TEMA: Programación por elementos finitos para el elemento barra 3D (Truss 3D).

PROBLEMA:

CODIGO:

clear
clc
%% Datos de armadura 3D
coord=[0 0 2000;4000 0 0;0 0 -2000;0 3000 0];
conect=[1 2;2 3;2 4];
f=[0 0 0 0 -10000 -12000 0 0 0 0 0 0]';
nfijos=[1 2 3 7 8 9 10 11 12]';
U=[0 0 0 0 0 0 0 0 0 0 0 0]';
area=[100;100;100];
mod_elas=[200e3;200e3;200e3];
k1=[1 -1;-1 1]; % Matriz de rigidez elemental en ejes locales

%% Rigidez axial
r_axial=area.*mod_elas;

%% Ensamble Matriz Global


nnodos=size(coord,1);
num_elem=size(conect,1);
GDL=3*nnodos;
K=zeros(GDL);
for ii=1:num_elem
% Nodos del elementos
i=conect(ii,1);
j=conect(ii,2);
ESCUELA POLITÉCNICA NACIONAL
FACULTAD DE INGENIERÍA MECÁNICA
ELEMENTOS FINITOS

% coordenadas nodales del elemento


xn=[coord(i,1) coord(j,1)]';
yn=[coord(i,2) coord(j,2)]';
zn=[coord(i,3) coord(j,3)]';

% Coordenadas nodales del elemento


xi=xn(1); yi=yn(1); zi=zn(1);
xj=xn(2); yj=yn(2); zj=zn(2);
L=sqrt((xj-xi)^2+(yj-yi)^2+(zj-zi)^2);

cx=(xj-xi)/L; %Coseno director X


cy=(yj-yi)/L; %Coseno director Y
cz=(zj-zi)/L; %Coseno director Z

% matrices de transformacion
R=[cx cx cz 0 0 0; 0 0 0 cx cy cz];
X=mod_elas(1)*[-1/L(1 ) 1/L(1 )]*R;

EA=r_axial(ii);
k(ii)=(EA/L);
%matriz de rigidez del elemento en coordenadas globales
ke=(R')*((EA/L)*k1*R);

in=[3*i-2 3*i-1 3*i 3*j-2 3*j-1 3*j];


K(in,in)=K(in,in)+ke;
end
disp(K)

%% SOLUCION DEL SISTEMA


%Aplicacion de condiciones de contorno
ic=nfijos;
id=[1:GDL];%Define los indices de los GDL desconocidos
id(ic)=[];
Kmod=K([id, ic'],[id, ic']);
dim=size(id, 2);
K11=Kmod([1:dim],[1:dim]);
UD=K11\f(id);

U(4)=UD(1);
U(5)=UD(2);
U(6)=UD(3);
U
Fuerzas=Kmod*U
L=[447.21 500 447.21];

%% Esfuerzos
X1=[-40.0000 -40.0000 20.0000 40.0000 0 -20.0000];
X2=[40.0000 40.0000 20.0000 -40.0000 0 -20.0000];
X3=[32 32 0 -32 24 0];
Um=U(1:6);
Esf1=X1*Um
Esf2=X2*Um
Esf3=X3*Um

%% R obtenidos de cada elemento


R1=[0.8944 0.8944 -0.4472 0 0 0; 0 0 0 0.8944 0 -0.4472];
R2=[-0.8944 -0.8944 -0.4472 0 0 0; 0 0 0 -0.8944 0 -0.4472];
ESCUELA POLITÉCNICA NACIONAL
FACULTAD DE INGENIERÍA MECÁNICA
ELEMENTOS FINITOS

R3=[-0.8000 -0.8000 0 0 0 0 ; 0 0 0 -0.8000 0.6000 0];

%% Fuerzas elementales
% Para:
disp('FUERZAS ELEMENTALES')
Elemento1=k(1)*R1*Um
Elemento2=k(2)*R2*Um
Elemento3=k(3)*R3*Um

RESULTADOS

Matriz de rigidez

Desplazamientos, fuerzas y esfuerzos: [mm] [N] [Pa] respectivamente


ESCUELA POLITÉCNICA NACIONAL
FACULTAD DE INGENIERÍA MECÁNICA
ELEMENTOS FINITOS

También podría gustarte