BELKA

////belka

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

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

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

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

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

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

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

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

Edof=[1,1,2;2,2,3;3,3,4]

ex1=[0,30]

ex2=[30,60]

ex3=[60,90;]

ep1=[200000000,0.0004]

ep2=[200000000,0.0004]

ep3=[200000000,0.0004]

K=zeros(4,4)

Ke1=bar1e(ex1,ep1)

Ke2=bar1e(ex2,ep2)

Ke3=bar1e(ex2,ep2)

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

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

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

f=zeros(4,1)

bc=[1,0;4,0]

f(3)=-25

[a,R]=solveq(K,f,bc)

ex=[ex1,ex2,ex3]

ep=[ep1,ep2,ep3]

ed=extract(Edof,a)

es1=bar1s(ex1,ep,ed(1,:))

es2=bar1s(ex2,ep,ed(2,:))

es3=bar1s(ex3,ep,ed(3,:))