The finite element discretization of a plane stress problem gives the stiffness

Question

The finite element discretization of a plane stress problem gives the stiffness matrix [K] and the load vector {F} as presented. If the displacement vector {x} is related to the stiffness matrix and the load vector by the equation: [K].{x}={F}

- Solve the displacement vector {x} by using the naive Gauss elimination method with scaling. ( MATLAB CODE )

The finite element discretization of a plane stress problem gives the stiffness matrix [K] and the load vector {F} as presented. If the displacement vector {x} is related to the stiffness matrix and the load vector by the equation: [K].{x}={F} - Solve the displacement vector {x} by using the naive Gauss elimination method with scaling. ( MATLAB CODE ) - Submit your program code together with your results for the numerical values of {x}.

Explanation / Answer

function [B pr]=spp(A,b)

B=[A b];

[n np1]=size(B);

pr=1:1:n;

s=max(abs(A'))';

for i=1:n [v ind]=max(abs(B(i:n,i))./s(i:n));

p=i+ind-1; [pr(i) pr(p)]=swap(pr(i),pr(p));

[B(i,:) B(p,:)]=swap(B(i,:),B(p,:));

[s(i) s(p)]=swap(s(i),s(p));

for j=i+1:n

d=B(j,i)/B(i,i);

B(j,:)=B(j,:)-d*B(i,:);

end

end

return

function [a,b]=swap(a,b)

temp=a;

a=b;

b=temp;

return

Related Questions

Navigate

• Browse (All)
• Browse T
• Subjects

---

---

Integrity-first tutoring: explanations and feedback only — we do not complete graded work. Learn more.