Bonjour, cliquez-ici pour vous inscrire et participer au forum.
  • Login:



+ Répondre à la discussion
Affichage des résultats 1 à 2 sur 2

Élément finis sur Matlab

  1. piero57

    Date d'inscription
    janvier 2018
    Messages
    2

    Élément finis sur Matlab

    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

    -----

     


    • Publicité



  2. piero57

    Date d'inscription
    janvier 2018
    Messages
    2

    Programme matlab

    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 ?

    Code:
    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
    Dernière modification par Antoane ; 04/01/2018 à 14h33. Motif: Ajout balises code
     





Sur le même thème :


    301 Moved Permanently

    301 Moved Permanently


    nginx/1.2.1



 

Discussions similaires

  1. Méthode des élement finis 1D
    Par Rawazs dans le forum Physique
    Réponses: 0
    Dernier message: 23/01/2013, 15h17
  2. les élément finis (conditions aux limites de Robin)
    Par narakphysics dans le forum Mathématiques du supérieur
    Réponses: 0
    Dernier message: 22/04/2012, 12h33
  3. calcul d'élément finis sur catia
    Par florian021 dans le forum Technologies
    Réponses: 2
    Dernier message: 21/07/2011, 20h53
  4. modélisation élément finis
    Par yacine_ dans le forum Physique
    Réponses: 1
    Dernier message: 03/04/2011, 13h34
  5. Possible en element finis?
    Par tomtom550 dans le forum Technologies
    Réponses: 1
    Dernier message: 26/03/2008, 00h55