اصلاح پروژه متلب جابجای گره های المان cstمثلثی سه گرهی خطی دو بعدی
پنج سال پیش منتشر شده
تعداد بازدید: 225
کد پروژه: 230464
شرح پروژه
پروژه متلب جابجای گره های المان cstمثلثی سه گرهی خطی دو بعدی اجزا محدود که بتونه هر شکلی رو محاسبه کند
کد خودم درس کار نمیکنه زمانمم محدوده تا 26 است
clc;
clear all;
nnd = 21 ; % Number of nodes:
nel = 24 ; % Number of elements:
nne = 3 ; % Number of nodes per element:
nodof =2; % Number of degrees of freedom per node
eldof = nne*nodof; % Number of degrees of freedom per element
%
geom = zeros(nnd,2);
%
geom = [ 0, -10; ... % Node 1
0, 0; ... % Node 2
0, 10; ... % Node 3
10, -10; ... % Node 4
10, 0; ... % Node 5
10, 10; ... % Node 6
20, -10; ... % Node 7
20, 0; ... % Node 8
20, 10; ... % Node 9
30, -10; ... % Node 10
30, 0; ... % Node 11
30, 10; ... % Node 12
40, -10; ... % Node 13
40, 0; ... % Node 14
40, 10; ... % Node 15
50, -10; ... % Node 16
50, 0; ... % Node 17
50, 10; ... % Node 18
60, -10; ... % Node 19
60, 0; ... % Node 20
60, 10]; % Node 21
%
disp(geom);
% Element connectivity
%
connec=zeros(nel,3);
connec = [ 1, 4, 2; ...% Element 1
4, 5, 2; ...% Element 2
2, 5, 3; ...% Element 3
5, 6, 3; ...% Element 4
4, 7, 5; ...% Element 5
7, 8, 5; ...% Element 6
5, 8, 6; ...% Element 7
8, 9, 6; ...% Element 8
7, 10, 8; ...% Element 9
10, 11, 8; ...% Element 10
8, 11, 9; ...% Element 11
11, 12, 9; ...% Element 12
10, 13, 11; ...% Element 13
13, 14, 11; ...% Element 14
11, 14, 12; ...% Element 15
14, 15, 12; ...% Element 16
13, 16, 14; ...% Element 17
16, 17, 14; ...% Element 18
14, 17, 15; ...% Element 19
17, 18, 15; ...% Element 20
16, 19, 17; ...% Element 21
19, 20, 17; ...% Element 22
17, 20, 18; ...% Element 23
20, 21, 18]; % Element 24
%
disp(connec);
% Material
%
E = 200000.; % Elastic modulus in MPa
vu = 0.3; % Poisson’s ratio
thick = 5.; % Beam thickness in mm
%
% Form the elastic matrix for plane stress
%
dee =[E/(1-vu^2) E*vu/(1-vu^2) 0;E*vu/(1-vu^2) E/(1-vu^2) 0;0 0 (1-vu)/2];
disp(dee);
%
% Boundary conditions
%
nf = ones(nnd, nodof); % Initialize the matrix nf to 1
nf(19,1) = 0; nf(19,2) = 0; % Prescribed nodal freedom of node 19
nf(20,1) = 0; nf(20,2) = 0; % Prescribed nodal freedom of node 20
nf(21,1) = 0; nf(21,2) = 0; % Prescribed nodal freedom of node 21
%
disp(nf);
% Counting of the free degrees of freedom
%
n=0;
for i=1:nnd
for j=1:nodof
if nf(i,j) ~= 0
n=n+1;
nf(i,j)=n;
end
end
end
%
% loading
%
Nodal_loads= zeros(nnd, 2);
%
Nodal_loads(2,1) = 0.; Nodal_loads(2,2) = -1000.; % Node 2
%
disp(Nodal_loads);
%%%%%%%%%%%% End of input %%%%%%%%%%%%
fg=zeros(n,1);
for i=1: nnd
if nf(i,1) ~= 0
fg(nf(i,1))= Nodal_loads(i,1);
end
if nf(i,2) ~= 0;
fg(nf(i,2))= Nodal_loads(i,2);
end
end
kk = zeros(n, n);
%
x1 = geom(connec(i,1),1); y1 = geom(connec(i,1),2);
x2 = geom(connec(i,2),1); y2 = geom(connec(i,2),2);
x3 = geom(connec(i,3),1); y3 = geom(connec(i,3),2);
%
A = (0.5)*det([1 x1 y1; ...
1 x2 y2; ...
1 x3 y3]);
%
m11 = (x2*y3 - x3*y2)/(2*A);
m21 = (x3*y1 - x1*y3)/(2*A);
m31 = (x1*y2 - y1*x2)/(2*A);
m12 = (y2 - y3)/(2*A);
m22 = (y3 - y1)/(2*A);
m32 = (y1 - y2)/(2*A);
m13 = (x3 - x2)/(2*A);
m23 = (x1 - x3)/(2*A);
m33 = (x2 -x1)/(2*A);
%
bee = [ m12 0 m22 0 m32 0; ...
0 m13 0 m23 0 m33; ...
m13 m12 m23 m22 m33 m32] ;
disp(bee);
l=0;
for k=1:nne
for j=1:nodof
l=l+1;
g(l)=nf(connec(i,k),j);
end
end
for i=1:nel
ke=thick*A*bee'*dee*bee; % Compute stiffness matrix%
for i=1:eldof
if g(i) ~= 0
for j=1: eldof
if g(j) ~= 0
KK(g(i),g(j))= ke(g(i),g(j)) + ke(i,j);
end
end
end
end
end
disp(ke);
disp(kk);
این پروژه شامل 5 فایل مهم است، لطفا قبل از ارسال پیشنهاد حتما نسبت به بررسی این فایل اقدام فرمایید.
مهارت ها و تخصص های مورد نیاز
مهلت برای انجام
1روز
وضعیت مناقصه
بسته
درباره کارفرما
عضویت پنج سال پیش
