daily pastebin goal
44%
SHARE
TWEET

Untitled

a guest Apr 21st, 2017 42 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. M_g(indice,indice) =...
  2. M_g(indice,indice) + (H' * H * thickness * dw_b * rho_c);
  3.  
  4. M_g(indice+tot_node,indice+tot_node) =...
  5. M_g(indice+tot_node,indice+tot_node) + (H' * H * thickness * I * dw_b * rho_c);
  6.  
  7. M_g(indice+2*tot_node,indice+2*tot_node) =...
  8. M_g(indice+2*tot_node,indice+2*tot_node) + (H' * H * thickness * I * dw_b * rho_c);
  9.  
  10. M_g(indice+3*tot_node,indice+3*tot_node) =...
  11. M_g(indice+3*tot_node,indice+3*tot_node) + (H' * H * thickness * dw_b * rho_c);
  12.  
  13. M_g(indice+4*tot_node,indice+4*tot_node) =...
  14. M_g(indice+4*tot_node,indice+4*tot_node) + (H' * H * thickness * dw_b * rho_c);
  15.  
  16. M_g: global mass matrix
  17. H: shape function matrix
  18. thickness: total thickness of plate
  19. I: thickness^3/12
  20. rho_c: density of the plate
  21. dw_b: detjacob * gaussweights
  22.    
  23. eig(stiffness(activeDof, activeDof), M_g(activeDof, activeDof)
RAW Paste Data
Top