Multidimensional Arrays Global stiff matrix

Do anyone know how to make a Multidimensional Arrays GSM, cause I have already make a 2x2 Local Stiff matrix. But I am confuse on make a GSM with for loop. And here is my code for LSM and vaeiables. Thank you
%Step 1 Discretize
%Assign Variables
E = young;
A = area;
L = leng;
F1 = F;
Num_elements = no_e;
l = L/Num_elements ;
K = E*A/l ;
%step 2 Local Behavior
% Form local stiffness Matrixes
LSM = zeros(2, 2, Num_elements);
for x=1: 1: Num_elements
LSM(1,1,x)= K(x,1);
LSM(1,2,x)= -K(x,1);
LSM(2,1,x)= -K(x,1);
LSM(2,2,x)= K(x,1);
end
%Step 3
%Assemble Global Stiffness Matrix (GSM)

 Akzeptierte Antwort

Sulaymon Eshkabilov
Sulaymon Eshkabilov am 8 Sep. 2020

0 Stimmen

Substitute this portion of your code:
LSM = zeros(2, 2, Num_elements);
for x=1: Num_elements
LSM(1,1)= K;
LSM(1,2)= -K;
LSM(2,1)= -K;
LSM(2,2)= K;
end
With this one:
LSM(1,1)= K;
LSM(1,2)= -K;
LSM(2,1)= -K;
LSM(2,2)= K;

Weitere Antworten (1)

Sulaymon Eshkabilov
Sulaymon Eshkabilov am 8 Sep. 2020

0 Stimmen

Here is the GSM created in loops:
LSM(1,1)= K;
LSM(1,2)= -K;
LSM(2,1)= -K;
LSM(2,2)= K;
GSM = zeros(Num_elements);
for ii=1:Num_elements-1
ROW = ii + [0 1];
COL = ROW;
GSM(ROW, COL) = GSM(ROW, COL)+LSM;
end

1 Kommentar

do the COL and ROW need to put the number i want. cause when I run the code, it give me an error say that
lab_2_part_3(50, 60, 55, 60,5)
Unable to perform assignment because the size of the left side is 2-by-2 and the size of the right side is 2-by-2-by-5.
Error in lab_2_part_3 (line 38)
GSM(ROW, COL) = GSM(ROW, COL)+LSM;
my whole code is like this ;
function [ displacements, external_Forces] = lab_2_part_3(E ,A ,L ,F1 ,Num_elements)
%Step 1 Discretize
%Assign Variables
E = young;
A = area;
L = len;
F1 = force;
Num_elements = no_ele;
l = L/Num_elements ;
K = E*A/l ;
%step 2 Local Behavior
% Form local stiffness Matrixes
LSM = zeros(2, 2, Num_elements);
for x=1: Num_elements
LSM(1,1)= K;
LSM(1,2)= -K;
LSM(2,1)= -K;
LSM(2,2)= K;
end
%Step 3
%Assemble Global Stiffness Matrix (GSM)
GSM = zeros(Num_elements);
for y=1:Num_elements-1
ROW = y + [0 1];
COL = ROW;
GSM(ROW, COL) = GSM(ROW, COL)+LSM;
end
%Step 4
%Apply Boundary Conditions
F1 = [60 ; 0 ; 0 ];
LSM = GSM (1:end-1 ,1:end-1);
% Step 5
%solve for unknown displacements
mini_U = LSM \ external_Forces;
%step 6
% Analyze and find Reaction Forces
displacements = mini_U;
external_Forces = [60 ; 0 ;0 ; -60];
end

Melden Sie sich an, um zu kommentieren.

Kategorien

Mehr zu Numerical Integration and Differential Equations finden Sie in Hilfe-Center und File Exchange

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by