matlab

[invitef0337726]

Bonsoir à tous !

Alors j'ai un problème concernant cet exercice, je ne vois pas du tout comment faire, si quelqu'un pourrait m'aider sa serait sympa !
En faite, c est mur soumis à une charge p3 et à une force F à son extrémité.

J'ai les codes sur matlab ci-joint :

function [x,U] = beam_approx(N, param, bc)

    % vector of equidistant x coordinates
    x = % COMPLETE HERE...

    % system matrix and rhs vector
    K = zeros(% COMPLETE HERE,% COMPLETE HERE); % the coefficient matrix
    p = zeros(% COMPLETE HERE,1);   % the right hand side

    % assemble finite element contributions
    for ii = 1:N-1
        idx = 3*(ii-1)+[1:6];
        [Ke, pe] = beam_matrix_vector(x(ii),x(ii+1),param);
        K(idx,idx) = K(idx,idx) + Ke;
        p(idx)     = p(idx)     + pe;
    end

    % add static boundary conditions
    % horizontal force
    p(bc.F1_idx) = p(bc.F1_idx) + bc.F1_val;
    % vertical force
    p(bc.F3_idx) = p(bc.F3_idx) + bc.F3_val;
    % moment
    p(bc.M_idx)  = p(bc.M_idx)  + bc.M_val ;

    % impose kinematic boundary conditions
    % horizontal displacement
    K(bc.u_idx,:) = 0;
    K(bc.u_idx,bc.u_idx) = eye(length(bc.u_idx));
    p(bc.u_idx) = bc.u_val;
    % vertical displacement
    K(bc.w_idx,:) = 0;
    K(bc.w_idx,bc.w_idx) = eye(length(bc.w_idx));
    p(bc.w_idx) = bc.w_val;
    % rotation
    K(bc.phi_idx,:) = 0;
    K(bc.phi_idx,bc.phi_idx) = eye(length(bc.phi_idx));
    p(bc.phi_idx) = bc.phi_val;

    % solve the resulting system of equations
    u = K\p;

    % return solved degrees of freedom row-wise per type
    U = [ u(1:3:3*N)'; u(2:3:3*N)'; u(3:3:3*N)' ];
end

et le 2eme pour la matrice raideur :

% provide finite element matrix and load vector

%% function: [x,U] = beam_matrix_vector(...)
%
% input:
%    xA    = x coordinate of A
%    xB    = x coordinate of B
%    param = parameter set of the problem
%
% output:
%    Ke    = element stiffness matrix
%    pe    = element load vector
%
function [Ke,pe] = beam_matrix_vector(xA, xB, param)

    % COMPLETE HERE...

end

Le sujet est joint.

Merci pour tout type d'aide
CDT
-----