How could i extract the mass matrix and the stiffness matrix from equation M *X''+K*X=0 what i am asking for is M(5X5) AND K(5X5) c

36 Ansichten (letzte 30 Tage)
How could i extract the mass matrix and the stiffness matrix from equation M *X''+K*X=0 what i am asking for is M(5X5) AND K(5X5) ci tried but what i get the matrix combiend betweeen M and K to M *X''+K*X=0 so i tried to make 2 equations but i think its worng because the stiffness matrix is not symmetric using the equationsToMatrix order from matlab

Antworten (1)

Vaidyanathan Thiagarajan
Vaidyanathan Thiagarajan am 2 Nov. 2017
Bearbeitet: Vaidyanathan Thiagarajan am 2 Nov. 2017
Hi Mohamed,
I am assuming you want to extract M and K matrices of a free vibration problem in MATLAB's PDE toolbox. You can use 'assembleFEMatrices' function to get the stiffness and mass matrix of a free vibration problem. For more details on 'assembleFEMatrices' please see the following link :
The following sample code extracts the M and K matrices of a free vibration problem using 'assembleFEMatrices' function :
N = 3;
model = createpde(N);
importGeometry(model,'Plate10x10x1.stl');
figure
hc = pdegplot(model,'FaceLabels','on');
hc(1).FaceAlpha = 0.5;
title('Plate with Face Labels');
E = 200e9; % Modulus of elasticity in Pascals
nu = .3; % Poisson's ratio
m = 8000; % Material density in kg/m^3
c = elasticityC3D(E,nu);
a = 0;
specifyCoefficients(model,'m',m,'d',0,'c',c,'a',a,'f',[0;0;0]);
applyBoundaryCondition(model,'mixed','Face',1:4,'u',0,'EquationIndex',3);
hmax = 1.2; % Maximum element length for a moderately fine mesh
generateMesh(model,'Hmax',hmax,'GeometricOrder','quadratic');
figure
pdeplot3D(model);
title('Mesh with Quadratic Tetrahedral Elements');
FEM = assembleFEMatrices(model);
%Mass Matrix
M = FEM.M;
%Stiffness Matrix
K = FEM.K;

Kategorien

Mehr zu Computational Geometry finden Sie in Help 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