Código elementos finitos matlab treliças

clear all

clc

E=207e9; A=0.01;

L1=PlaneTrussElementLength(1.82,2.86,0,0);

theta1=231.47;

L2=PlaneTrussElementLength(0,0,182,0);

theta2= 0;

L3=PlaneTrussElementLength(1.82,0,1.82,2.286);

theta3=90;

L4=PlaneTrussElementLength(1.82,2.286,3.64,2.286);

theta4=0;

L5=PlaneTrussElementLength(3.64,2.286,1.82,0);

theta5=231.47;

L6=PlaneTrussElementLength(1.82,0,3.64,0.76);

theta6=22.7;

L7=PlaneTrussElementLength(3.64,0.76,3.64,2.286);

theta7=90;

L8=PlaneTrussElementLength(3.64,2.286,5.46,2.286);

theta8=0;

L9=PlaneTrussElementLength(5.46,2.286,3.64,0.76);

theta9=219.86;

L10=PlaneTrussElementLength(3.64,0.76,5.46,1.52);

theta10=22.7;

L11=PlaneTrussElementLength(5.46,1.52,5.46,2.286);

theta11=90;

L12=PlaneTrussElementLength(5.46,2.286,7.28,2.286);

theta12=0;

L13=PlaneTrussElementLength(7.28,2.286,5.46,1.52);

theta13=202.82;

k1=PlaneTrussElementStiffness(E,A,L1,theta1)

k2=PlaneTrussElementStiffness(E,A,L2,theta2)

k3=PlaneTrussElementStiffness(E,A,L3,theta3)

k4=PlaneTrussElementStiffness(E,A,L4,theta4)

k5=PlaneTrussElementStiffness(E,A,L5,theta5)

k6=PlaneTrussElementStiffness(E,A,L6,theta6)

k7=PlaneTrussElementStiffness(E,A,L7,theta7)

k8=PlaneTrussElementStiffness(E,A,L8,theta8)

k9=PlaneTrussElementStiffness(E,A,L9,theta9)

k10=PlaneTrussElementStiffness(E,A,L10,theta10)

k11=PlaneTrussElementStiffness(E,A,L11,theta11)

k12=PlaneTrussElementStiffness(E,A,L12,theta12)

k13=PlaneTrussElementStiffness(E,A,L13,theta13)

K=zeros(16,16)

K=PlaneTrussAssemble(K,k1,1,5)

K=PlaneTrussAssemble(K,k2,5,6)

K=PlaneTrussAssemble(K,k3,6,1)

K=PlaneTrussAssemble(K,k4,1,2)

K=PlaneTrussAssemble(K,k5,2,6)

K=PlaneTrussAssemble(K,k6,6,7)

K=PlaneTrussAssemble(K,k7,7,2)

K=PlaneTrussAssemble(K,k8,2,3)

K=PlaneTrussAssemble(K,k9,3,7)

K=PlaneTrussAssemble(K,k10,7,8)

K=PlaneTrussAssemble(K,k11,8,3)

K=PlaneTrussAssemble(K,k12,3,4)

K=PlaneTrussAssemble(K,k13,4,8)

k=[K(1:8) K(11) K(13:16);

K(2,1:8) K(2,11) K(2,13:16);

K(3,1:8) K(3,11) K(3,13:16);

K(4,1:8) K(4,11) K(4,13:16);

K(5,1:8) K(5,11) K(5,13:16);

K(6,1:8) K(6,11) K(6,13:16);

K(7,1:8) K(7,11) K(7,13:16);

K(8,1:8) K(8,11) K(8,13:16);

K(11,1:8) K(11,11) K(11,13:16);

K(13,1:8) K(13,11) K(13,13:16);

K(14,1:8) K(14,11) K(14,13:16);

K(15,1:8) K(15,11) K(15,13:16);