brick_stress.m

结构力学中的有限元例子,包含了7个分类文件夹

function [SIGsys] = BRICK_stress (in_data, resp)

% stress field for 'brick' 3D FEM model
% ---------------------------------------------------------------------------
% INPUT:
%        ND  -  nodes matrix
%        D   -  computed nodal displacements vector
%        EL  -  elements matrix
% ---------------------------------------------------------------------------
fj = size(in_data.EL);
SIGsys = zeros(fj(1)*6,1);
for i=1:fj(1) 
   
   if in_data.EL(i,2)==6  % "6" - BRICK element (8-nodes) 
      node1 = find(in_data.ND(:,1)==in_data.EL(i,3));
      node2 = find(in_data.ND(:,1)==in_data.EL(i,4));
      node3 = find(in_data.ND(:,1)==in_data.EL(i,5));
      node4 = find(in_data.ND(:,1)==in_data.EL(i,6));
      node5 = find(in_data.ND(:,1)==in_data.EL(i,7));
      node6 = find(in_data.ND(:,1)==in_data.EL(i,8));
      node7 = find(in_data.ND(:,1)==in_data.EL(i,9));
      node8 = find(in_data.ND(:,1)==in_data.EL(i,10));
      Em    = in_data.EL(i,11);
      miu_1 = in_data.EL(i,12);
      rho_  = in_data.EL(i,13);
      
      [Klc,Bsys,Esys,V] = D3_LB (in_data.ND(node1,2),in_data.ND(node1,3),in_data.ND(node1,4),...
         in_data.ND(node2,2),in_data.ND(node2,3),in_data.ND(node2,4),in_data.ND(node3,2),in_data.ND(node3,3),...
         in_data.ND(node3,4),in_data.ND(node4,2),in_data.ND(node4,3),in_data.ND(node4,4),in_data.ND(node5,2),...
         in_data.ND(node5,3),in_data.ND(node5,4),in_data.ND(node6,2),in_data.ND(node6,3),in_data.ND(node6,4),...
         in_data.ND(node7,2),in_data.ND(node7,3),in_data.ND(node7,4),in_data.ND(node8,2),in_data.ND(node8,3),...
         in_data.ND(node8,4), Em, miu_1);
      Dlocal = [resp.static.D(node1*3-2) resp.static.D(node1*3-1) resp.static.D(node1*3) ...
              resp.static.D(node2*3-2) resp.static.D(node2*3-1) resp.static.D(node2*3)  ...
              resp.static.D(node3*3-2) resp.static.D(node3*3-1) resp.static.D(node3*3) ...
              resp.static.D(node4*3-2) resp.static.D(node4*3-1) resp.static.D(node4*3) ...
              resp.static.D(node5*3-2) resp.static.D(node5*3-1) resp.static.D(node5*3) ...
              resp.static.D(node6*3-2) resp.static.D(node6*3-1) resp.static.D(node6*3) ...
              resp.static.D(node7*3-2) resp.static.D(node7*3-1) resp.static.D(node7*3) ...
              resp.static.D(node8*3-2) resp.static.D(node8*3-1) resp.static.D(node8*3)];
      
      SIGlocal = Esys*Bsys*Dlocal';
      SIGsys((i*6-5):(i*6)) = SIGlocal;
   end;
   
end;