RAMA Z PRETEM

exec('/pub/programy/calfem/bar1e.sci');

exec('/pub/programy/calfem/bar2e.sci');

exec('/pub/programy/calfem/bar1s.sci');

exec('/pub/programy/calfem/bar2s.sci');

exec('/pub/programy/calfem/plani4s.sci'); 

exec('/pub/programy/calfem/plani4e.sci'); 

exec('/pub/programy/calfem/statcon.sci');

exec('/pub/programy/calfem/planqs.sci');

exec('/pub/programy/calfem/planqe.sci'); 

exec('/pub/programy/calfem/eldraw2.sci'); 

exec('/pub/programy/calfem/eldisp2.sci'); 

exec('/pub/programy/calfem/coordxtr.sci');

exec('/pub/programy/calfem/solveq.sci'); 

exec('/pub/programy/calfem/hooke.sci'); 

exec('/pub/programy/calfem/plante.sci'); 

exec('/pub/programy/calfem/assem.sci'); 

exec('/pub/programy/calfem/extract.sci');

exec('/pub/programy/calfem/plants.sci');

  Edof=[1 1 2 3 4 9 10 7 8;

2 3 4 5 6 11 12 9 10;

3 7 8 9 10 19 20 17 18;

4 9 10 11 12 21 22 19 20;

5 11 12 13 14 23 24 21 22;

6 13 14 15 16 25 26 23 24;]

 Edof1=[7 15 16 27 28]

ex=[0 2 2 0;2 4 4 2;0 2 2 0;2 4 4 2;4 6 6 4;6 8 8 6;]

ex1=[8,8]

ey=[0 0 4 4;0 0 4 4;4 4 6 6; 4 4 6 6;4 4 6 6;4 4 6 6;]

ey1=[4,0]

eldraw2(ex,ey,Edof(:,1))

eldraw2(ex1,ey1,Edof1(:,1))

ep=[2,0.2]

D=hooke(2,30000,0.16)

ep1=[30000,0.2]

K=zeros(28,28)

Ke1=planqe(ex(1,:),ey(1,:),ep,D);

Ke2=planqe(ex(2,:),ey(2,:),ep,D);

Ke3=planqe(ex(3,:),ey(3,:),ep,D);

Ke4=planqe(ex(4,:),ey(4,:),ep,D);

Ke5=planqe(ex(5,:),ey(5,:),ep,D);

Ke6=planqe(ex(6,:),ey(6,:),ep,D);

Ke7=bar2e(ex1(1,:),ey1(1,:),ep1)

K=assem(Edof(1,:),K,Ke1)

K=assem(Edof(2,:),K,Ke2)

K=assem(Edof(3,:),K,Ke3)

K=assem(Edof(4,:),K,Ke4)

K=assem(Edof(5,:),K,Ke5)

K=assem(Edof(6,:),K,Ke6)

K=assem(Edof1(1,:),K,Ke7)

f=zeros(28,1)

bc=[1 0;2 0;3 0;4 0;5 0;6 0;27 0;28 0]

f(25)=50

a=solveq(K,f,bc)

ed=extract(Edof,a)

ed1=extract(Edof1,a)

eldisp2(ex,ey,ed)

eldisp2(ex1,ey1,ed1)

N=bar2s(ex1,ey1,ep1,ed1) // LICZY SIŁĘ W PRĘCIEQ-bar

es=zeros(6,3)

for i=1:6

  [es,et]=planqs(ex(i,:),ey(i,:),ep,D,ed(i,:))

end

[es,et]=planqs(ex(1,:),ey(1,:),ep,D,ed)