belka MES matlab

Ćwiczenie projektowe 3

Wprowadzenie danych

E=210e3

E =

210000

I=4.0000e-005

I =

4.0000e-005

l=0.6

l =

0.6000

P1=10

P1 =

P2=6

P2 =

q=2

q =

le=4

le =

N=10

N =

lw=le+1

lw =

Ee=[E,E,E,E]

Ee =

210000 210000 210000 210000

Ie=[I,I,I,I]

Ie =

1.0e-004 *

0.4000 0.4000 0.4000 0.4000

lE=[2*l,2*l,2*l,l]

lE =

1.2000 1.2000 1.2000 0.6000

ien=[1 2;2 3;3 4;4 5]'

ien =

1 2 3 4

2 3 4 5

lw=le+1

lw =

k1=sztywelbelki(1,Ee,Ie,lE)

k1 =

58.3333 35.0000 -58.3333 35.0000

35.0000 28.0000 -35.0000 14.0000

-58.3333 -35.0000 58.3333 -35.0000

35.0000 14.0000 -35.0000 28.0000

k2=sztywelbelki(2,Ee,Ie,lE)

k2 =

58.3333 35.0000 -58.3333 35.0000

35.0000 28.0000 -35.0000 14.0000

-58.3333 -35.0000 58.3333 -35.0000

35.0000 14.0000 -35.0000 28.0000

k3=sztywelbelki(3,Ee,Ie,lE)

k3 =

58.3333 35.0000 -58.3333 35.0000

35.0000 28.0000 -35.0000 14.0000

-58.3333 -35.0000 58.3333 -35.0000

35.0000 14.0000 -35.0000 28.0000

k4=sztywelbelki(4,Ee,Ie,lE)

k4 =

466.6667 140.0000 -466.6667 140.0000

140.0000 56.0000 -140.0000 28.0000

-466.6667 -140.0000 466.6667 -140.0000

140.0000 28.0000 -140.0000 56.0000

K=zeros (N,N)

K =

0 0 0 0 0 0 0 0 0 0

lm=pomocagregacji (2,2,lw,ien,le)

lm =

1 3 5 7

2 4 6 8

3 5 7 9

4 6 8 10

K=agregacjamacglob (K,1,k1,2,2,le,lw,lm);

K=agregacjamacglob (K,2,k2,2,2,le,lw,lm);

K=agregacjamacglob (K,3,k3,2,2,le,lw,lm);

K=agregacjamacglob (K,4,k4,2,2,le,lw,lm)

K =

Columns 1 through 7

58.3333 35.0000 -58.3333 35.0000 0 0 0

35.0000 28.0000 -35.0000 14.0000 0 0 0

-58.3333 -35.0000 116.6667 0 -58.3333 35.0000 0

35.0000 14.0000 0 56.0000 -35.0000 14.0000 0

0 0 -58.3333 -35.0000 116.6667 0 -58.3333

0 0 35.0000 14.0000 0 56.0000 -35.0000

0 0 0 0 -58.3333 -35.0000 525.0000

0 0 0 0 35.0000 14.0000 105.0000

0 0 0 0 0 0 -466.6667

0 0 0 0 0 0 140.0000

Columns 8 through 10

0 0 0

35.0000 0 0

14.0000 0 0

105.0000 -466.6667 140.0000

84.0000 -140.0000 28.0000

-140.0000 466.6667 -140.0000

28.0000 -140.0000 56.0000

fr1=[-P1/2; -P1*lE(1)/8;-P1/2;P1*lE(1)/8]

fr1 =

-5.0000

-1.5000

-5.0000

1.5000

fr2=[-P2/2; -P2*lE(2)/8;-P2/2;P1*lE(2)/8]

fr2 =

-3.0000

-0.9000

-3.0000

1.5000

fr3=[-q*lE(3)/2; -q*lE(3)*lE(3)/12; -q*lE(3)/2; q*lE(3)*lE(3)/12]

fr3 =

-1.2000

-0.2400

-1.2000

0.2400

fr4=[-q*lE(4)/2; -q*lE(4)*lE(4)/12; -q*lE(4)/2; q*lE(4)*lE(4)/12]

fr4 =

-0.6000

-0.0600

-0.6000

0.0600

V=zeros (N,1);

V=agregacjawekobc (V,1,2,2,lm,fr1)

V =

-5.0000

-1.5000

-5.0000

1.5000

V=agregacjawekobc (V,2,2,2,lm,fr2)

V=agregacjawekobc (V,3,2,2,lm,fr3)

V =

-5.0000

-1.5000

-5.0000

1.5000

-1.2000

-0.2400

-1.2000

0.2400

V=agregacjawekobc (V,4,2,2,lm,fr4)

V =

-5.0000

-1.5000

-5.0000

1.5000

-1.2000

-0.2400

-1.8000

0.1800

-0.6000

0.0600

P=[0 0 0 0 0 0 0 0 0 0 ]

P =

0 0 0 0 0 0 0 0 0 0

K (1,:)=[]

K =

Columns 1 through 7

35.0000 28.0000 -35.0000 14.0000 0 0 0

-58.3333 -35.0000 116.6667 0 -58.3333 35.0000 0

35.0000 14.0000 0 56.0000 -35.0000 14.0000 0

0 0 -58.3333 -35.0000 116.6667 0 -58.3333

0 0 35.0000 14.0000 0 56.0000 -35.0000

0 0 0 0 -58.3333 -35.0000 525.0000

0 0 0 0 35.0000 14.0000 105.0000

0 0 0 0 0 0 -466.6667

0 0 0 0 0 0 140.0000

Columns 8 through 10

0 0 0

35.0000 0 0

14.0000 0 0

105.0000 -466.6667 140.0000

84.0000 -140.0000 28.0000

-140.0000 466.6667 -140.0000

28.0000 -140.0000 56.0000

K(:,1)=[]

K =

Columns 1 through 7

28.0000 -35.0000 14.0000 0 0 0 0

-35.0000 116.6667 0 -58.3333 35.0000 0 0

14.0000 0 56.0000 -35.0000 14.0000 0 0

0 -58.3333 -35.0000 116.6667 0 -58.3333 35.0000

0 35.0000 14.0000 0 56.0000 -35.0000 14.0000

0 0 0 -58.3333 -35.0000 525.0000 105.0000

0 0 0 35.0000 14.0000 105.0000 84.0000

0 0 0 0 0 -466.6667 -140.0000

0 0 0 0 0 140.0000 28.0000

Columns 8 through 9

0 0

-466.6667 140.0000

-140.0000 28.0000

466.6667 -140.0000

-140.0000 56.0000

K (2,:)=[]

K =

Columns 1 through 7

28.0000 -35.0000 14.0000 0 0 0 0

14.0000 0 56.0000 -35.0000 14.0000 0 0

0 -58.3333 -35.0000 116.6667 0 -58.3333 35.0000

0 35.0000 14.0000 0 56.0000 -35.0000 14.0000

0 0 0 -58.3333 -35.0000 525.0000 105.0000

0 0 0 35.0000 14.0000 105.0000 84.0000

0 0 0 0 0 -466.6667 -140.0000

0 0 0 0 0 140.0000 28.0000

Columns 8 through 9

0 0

-466.6667 140.0000

-140.0000 28.0000

466.6667 -140.0000

-140.0000 56.0000

K (:,2)=[]

K =

Columns 1 through 7

28.0000 14.0000 0 0 0 0 0

14.0000 56.0000 -35.0000 14.0000 0 0 0

0 -35.0000 116.6667 0 -58.3333 35.0000 0

0 14.0000 0 56.0000 -35.0000 14.0000 0

0 0 -58.3333 -35.0000 525.0000 105.0000 -466.6667

0 0 35.0000 14.0000 105.0000 84.0000 -140.0000

0 0 0 0 -466.6667 -140.0000 466.6667

0 0 0 0 140.0000 28.0000 -140.0000

Column 8

140.0000

28.0000

-140.0000

56.0000

K (3,:)=[]

K =

Columns 1 through 7

28.0000 14.0000 0 0 0 0 0

14.0000 56.0000 -35.0000 14.0000 0 0 0

0 14.0000 0 56.0000 -35.0000 14.0000 0

0 0 -58.3333 -35.0000 525.0000 105.0000 -466.6667

0 0 35.0000 14.0000 105.0000 84.0000 -140.0000

0 0 0 0 -466.6667 -140.0000 466.6667

0 0 0 0 140.0000 28.0000 -140.0000

Column 8

140.0000

28.0000

-140.0000

56.0000

K (:,3)=[]

K =

28.0000 14.0000 0 0 0 0 0

14.0000 56.0000 14.0000 0 0 0 0

0 14.0000 56.0000 -35.0000 14.0000 0 0

0 0 -35.0000 525.0000 105.0000 -466.6667 140.0000

0 0 14.0000 105.0000 84.0000 -140.0000 28.0000

0 0 0 -466.6667 -140.0000 466.6667 -140.0000

0 0 0 140.0000 28.0000 -140.0000 56.0000

K (4,:)=[]

K =

28.0000 14.0000 0 0 0 0 0

14.0000 56.0000 14.0000 0 0 0 0

0 14.0000 56.0000 -35.0000 14.0000 0 0

0 0 14.0000 105.0000 84.0000 -140.0000 28.0000

0 0 0 -466.6667 -140.0000 466.6667 -140.0000

0 0 0 140.0000 28.0000 -140.0000 56.0000

K (:,4)=[]

K =

28.0000 14.0000 0 0 0 0

14.0000 56.0000 14.0000 0 0 0

0 14.0000 56.0000 14.0000 0 0

0 0 14.0000 84.0000 -140.0000 28.0000

0 0 0 -140.0000 466.6667 -140.0000

0 0 0 28.0000 -140.0000 56.0000

K (5,:)=[]

K =

28.0000 14.0000 0 0 0 0

14.0000 56.0000 14.0000 0 0 0

0 14.0000 56.0000 14.0000 0 0

0 0 14.0000 84.0000 -140.0000 28.0000

0 0 0 28.0000 -140.0000 56.0000

K (:,5)=[]

K =

28.0000 14.0000 0 0 0

14.0000 56.0000 14.0000 0 0

0 14.0000 56.0000 14.0000 0

0 0 14.0000 84.0000 28.0000

0 0 0 28.0000 56.0000

V(1)=[]

V =

-1.5000

-5.0000

1.5000

-1.2000

-0.2400

-1.8000

0.1800

-0.6000

0.0600

V(2)=[]

V =

-1.5000

1.5000

-1.2000

-0.2400

-1.8000

0.1800

-0.6000

0.0600

V(3)=[]

V =

-1.5000

1.5000

-0.2400

-1.8000

0.1800

-0.6000

0.0600

V(4)=[]

V =

-1.5000

1.5000

-0.2400

0.1800

-0.6000

0.0600

V(5)=[]

V =

-1.5000

1.5000

-0.2400

0.1800

0.0600

U=K\V

U =

-0.0792

0.0512

-0.0186

0.0059

-0.0019

u1=[0,U(1),0, U(2)]'

u1 =

-0.0792

0.0512

F1=-fr1+k1*u1

F1 =

4.0213

0.0000

5.9787

-1.1744

u2=[0,U(2),0,U(3)]'

u2 =

0.0512

-0.0186

F2=-fr2+k2*u2

F2 =

4.1433

2.0744

1.8567

-1.3024

u3=[0,U(3),0,U(4)]'

u3 =

-0.0186

0.0059

F3=-fr3+k3*u3

F3 =

0.7555

-0.1976

1.6445

-0.3359

u4=[0,U(4),0,U(5)]'

u4 =

0.0059

-0.0019

F4=-fr4+k4*u4

F4 =

1.1598

0.3359

0.0402

0.0000

Wyszukiwarka

Podobne podstrony:
3 BELKA MES
Belka MES id 82481 Nieznany
Belka MES MRS
LAB1 MES MATLAB statyka 1 id 25 Nieznany
belka mes
Belka MES
BELKA DRGANIA MATLAB
3 BELKA MES
5 3 Zał 1 MES Belka na gruncie Zadania pomocnicze
belka matlab
Matlab cw1 2 zaoczni
cz 1, Matlab moj
Image Processing with Matlab 33
Hala CECHOWANIE BELKA SPRĘŻONA ok
MATLAB graf(1)
kod matlab
Cw08 Matlab2

więcej podobnych podstron