76 end 77 end
78 % DEFINE LOADSAND SUPPORTS(HALF MBB-BEAM) 79 F(2,1) = -1;
80 fixeddofs = union([1:2:2*(nely+1)], [2*(nelx+1)*(nely+1)]);
81 alldofs = [1:2*(nely+1)*(nelx+1)]; 82 freedofs = setdiff(alldofs,fixeddofs); 83 % SOLVING
84 U(freedofs,:) = K(freedofs,freedofs) \\ F(freedofs,:);
85 U(fixeddofs,:)= 0;
86 %%%%%%%%%% ELEMENT STIFFNESS MATRIX %%%%%%% 87 function [KE]=lk 88 E = 1.; 89 nu = 0.3;
90 k=[ 1/2-nu/6 1/8+nu/8 -1/4-nu/12 -1/8+3*nu/8 ... 91 -1/4+nu/12 -1/8-nu/8 nu/6 1/8-3*nu/8]; 92 KE = E/(1-nu^2)*
[ k(1) k(2) k(3) k(4) k(5) k(6) k(7) k(8) 93 k(2) k(1) k(8) k(7) k(6) k(5) k(4) k(3) 94 k(3) k(8) k(1) k(6) k(7) k(4) k(5) k(2) 95 k(4) k(7) k(6) k(1) k(8) k(3) k(2) k(5) 96 k(5) k(6) k(7) k(8) k(1) k(2) k(3) k(4) 97 k(6) k(5) k(4) k(3) k(2) k(1) k(8) k(7) 98 k(7) k(4) k(5) k(2) k(3) k(8) k(1) k(6) 99 k(8) k(3) k(2) k(5) k(4) k(7) k(6) k(1)];
