1 / 15

Finite Element Method

Finite Element Method. Equations assembling Displacement boundary conditions Gaussian elimination. Monday, 11/11/2002. Corrections to Homework (node coordinates). We have 6 nodes, their coordinates are coor=(0.1,0.1; 0,0.6; 0.6,0; 0.5,0.4; 1.2 ,0.2; 1.2 ,0.5);.

fagan
Download Presentation

Finite Element Method

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Finite Element Method Equations assembling Displacement boundary conditions Gaussian elimination Monday, 11/11/2002

  2. Corrections to Homework(node coordinates) We have 6 nodes, their coordinates are coor=(0.1,0.1; 0,0.6; 0.6,0; 0.5,0.4; 1.2,0.2; 1.2,0.5);

  3. Corrections to Homework(element area) • In MATLAB with the above definite of elmt and coor, the area can be calculated with • A=0.5*det([1, coor(elmt(e,1),1), coor(elmt(e,1),2); … • 1, coor(elmt(e,2),1), coor(elmt(e,2),2); … • 1, coor(elmt(e,3),1), coor(elmt(e,3),2)])

  4. System Equations

  5. System Equations (matrix form)

  6. Assembled System Equations N: number of nodes

  7. Symmetric Stiffness Matrix

  8. Proof of Symmetric Stiffness Matrix

  9. MATLAB form If T_x, F_x, T_y, and F_y are column vector F = [T_x+F_x; T_y+F_y] K = [Kxx,Kxy; Kyx,Kyy]

  10. MATLAB Matrix Equation

  11. Displacement Boundary Conditions To solve: F = K * u u = gaussElim(K,F) u(5) u(6) u(11) u(12)

  12. Introducing Boundary Conditions

  13. Keep Stiffness Matrix Symmetric % apply displacement boundary condition for node 5 % in y direction ibound = 5 + nNode; for i = 1:2*nNode if(i ~= ibound) K(i,ibound)=0; K(ibound,i)=0; end end F(ibound)=0;

  14. Gaussian Elimination MATLAB Tutorial, Tutorial on subroutines. gaussElim.m:

  15. Tractions on Displacement Boundaries

More Related