Élément finis sur Matlab

[invitede1c8495]

Bonjour à tous,

Je souhaite programmer un fichier qui permet de retourner les matrices raideurs de chaque éléments d'une structure et qui les assemblent afin de donner la matrice raideur globale de celle-ci.
Pour finir, il faut également retourner la matrice force pour résoudre le système Ku=f

J'ai déjà programmé la structure sur matlab(graphiquement), qui s'agira d'un portique, or je ne sais pas comment programmer les parties que je viens d'évoquer précédemment.

Je vous remercie d'avance pour votre aide car je bloque vraiment.

En vous souhaitant une très bonne journée.

Pierre
-----

[invitede1c8495]

Bonjour à tous,

Je souhaiterai savoir si ce programme peut retourner la matrice raideur d'une structure ? Et comment pourrais-je faire pour avoir la matrice raideur d'un seul élément par contre ?

classdef FeAnalysisStatic < mafe.FeAnalysis
    % static analysis of finite element problem

    properties
        %
    end

    methods
        %% constructor
        function obj = FeAnalysisStatic(fep)
            if nargin > 0
                obj.fep = fep;
            end
        end
        %% calculate static response
        function analyse(self)
            % make sure the problem is initialised
            self.fep.init();
            % get system stiffness matrix and force vector
            K = self.fep.assembleSystemStiffnessMatrix();
            p = self.fep.assembleSystemForceVector();
            % apply constraints
            [f          ] = self.fep.applyConstraintsNeumann(p);
            [f, dbc, act] = self.fep.applyConstraintsDirichlet(K, f);
            % allocate solution vector
            u = zeros(self.fep.ndofs,1);
            % solve system of equations (active indices only)
            u(act) = K(act,act)\f(act);
            u(dbc) = f(dbc);
            % compute reaction forces
            r = K * u - p;
            % feedback of solution values (u and r) to nodal dofs
            self.fep.updateNodalDofs(u, r);
            % postprocess elements
            self.fep.postprocessElements(u);
        end
    end
end

Je vous remercie amplement pour l'aide apportée à mon problème