Scribd
Upload
169 views4 pages

Metode Numerice - Tema 1 Matlab

This document defines the parameters, geometry, and material properties for a structural analysis problem involving three beam elements. Stiffness matrices are formulated for each element and combined into a global structural stiffness matrix. The matrix is rearranged and partitioned to solve for displacements and reactive forces due to an applied load. Axial forces are then calculated for each element. Key results include a displacement of -0.332 for element 2, a reactive force of 6.919e3 for element 1, and axial stresses of 5.58, 11.415, and -10.147 MPa for elements 1, 2, and 3 respectively.

Uploaded by

Bogdan Bîngă
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
169 views4 pages

Metode Numerice - Tema 1 Matlab

This document defines the parameters, geometry, and material properties for a structural analysis problem involving three beam elements. Stiffness matrices are formulated for each element and combined into a global structural stiffness matrix. The matrix is rearranged and partitioned to solve for displacements and reactive forces due to an applied load. Axial forces are then calculated for each element. Key results include a displacement of -0.332 for element 2, a reactive force of 6.919e3 for element 1, and axial stresses of 5.58, 11.415, and -10.147 MPa for elements 1, 2, and 3 respectively.

Uploaded by

Bogdan Bîngă
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 4

ORIGIN 1

n 21
3

P1 ( 5 0 .9 n ) 1 0

P2 1.5P1
l1 ( 1.4 0.5n
) 10

l2 1.2l1

l3 0.8l1

A1 ( 4 0.4n
) 1 0

A2 1.2A1

E 2 1 0

mm

A3 1.5A1

l1 1 .1 9 1 0

A1 1 .2 4 1 0

l2 1 .4 2 8 1 0
3

l3 9 .5 2 1 0

1.Force displacement relation

1 1 A1
E
1 1 l1

1 1 A2
E
1 1 l2

k1

k2

2 .0 84 1 04 2 .0 84 1 04

4
4
2 .0 84 1 0 2 .0 84 1 0

1 1 A3
E
1 1 l3

k3

2 .0 84 1 04 2 .0 84 1 04

k2

4
4
2 .0 84 1 0 2 .0 84 1 0

k1

3 .9 08 1 04 3 .9 08 1 04

k3

4
4
3 .9 08 1 0 3 .9 08 1 0

2.Expanded matrices

2 .0 84 1 04 2 .0 84 1 04

4
4
k 1exp 2 .0 84 1 0 2 .0 84 1 0

0
0

0
0

0 0

0 0
0 0

0 0

0
0
0
0 2 .0 8 4 1 04 2 .0 8 4 1 04
k 2exp
0 2 .0 8 4 1 04 2 .0 8 4 1 04
0
0
0

0
0
k 3exp
0
0

4
4
0 3 .9 08 1 0 3 .9 08 1 0

4
4
0 3 .9 08 1 0
3 .9 08 1 0
0

3.Structural stiffness matrix

2 .0 84 1 04 2 .0 84 1 04

0
0

4
4
4
2 .0 84 1 0 4 .1 68 1 0 2 .0 84 1 0
0

K k 1exp k 2exp k 3exp

4
4
4
0
2 .0 84 1 0 5 .9 92 1 0 3 .9 08 1 0

4
4

0
0
3 .9 08 1 0 3 .9 08 1 0

4.Rearranging the matrix

2 .0 84 1 04 2 .0 84 1 04

0
0

4
4
0
0
3 .9 08 1 0 3 .9 08 1 0

K1 stack ( sub matrix( K 1114) sub matrix( K 4414) sub matrix( K 2314) )

4
4
4
0
2 .0 84 1 0 4 .1 68 1 0 2 .0 84 1 0

4
4
4

0
2 .0 84 1 0 5 .9 92 1 0 3 .9 08 1 0

4
2 .0 84 1 04

0
2 .0 84 1 0
0

4
4
0
3 .9 08 1 0
0
3 .9 08 1 0

K2 aug ment( sub matrix( K11411) sub matrix( K11444) sub matrix( K11423) )

4
4
4
0
4 .1 68 1 0 2 .0 84 1 0
2 .0 84 1 0

4
4
4

0
3 .9 08 1 0 2 .0 84 1 0 5 .9 92 1 0

4
2 .0 84 1 04

0
2 .0 84 1 0
0

4
4
0
3 .9 08 1 0
0
3 .9 08 1 0

K2

4
4
4
0
4 .1 68 1 0 2 .0 84 1 0
P1
2 .0 84 1 0

Pn

4
4
4

P2 4
0
3 .9 08 1 0 2 .0 84 1 0 5 .9 92 1 0

2 .3 9 1 0

Pn

4
3 .5 8 5 1 0

5.Partitioning the rearranged structural stiffness matrix, K2

2 .0 84 1 04

Krn sub matrix( K21 234)

4
3 .9 08 1 0
0

4 .1 68 1 04 2 .0 84 1 04

Kn n sub matrix( K23434)

4
4
2 .0 84 1 0 5 .9 92 1 0
6.Determination of the free displacements and reactive forces

0.332

0.483

Pn

Dn Knn

6 .9 1 9 1 03

Pr Krn Dn

4
1 .8 8 7 1 0
D1 0

D2 0.332
D3 0.483
D4 0

*Checking up:
Pr

1 1

Pr

2 1

Pn

1 1

Pn

2 1

12

7.276 10

7.The axial force in the finite elements:

D1 0

D2 0.332

dn1
3

N1 6 .9 191
0

6 .9 19 1 03

3
6 .9 19 1 0

Fn1 k 1 d n1

x1

N1
A1

5.58

D2 0.332

D3 0.483

dn2

1 .6 98 1 04

4
1 .6 98 1 0

Fn2 k 2 d n2

N2 Fn2

2 1

x2

N2
A2

1.698 10

11.415

D3 0.483

D4 0

dn3

1 .8 87 1 04

Fn3 k 3 d n3

4
1 .8 87 1 0
N3 Fn3

2 1

x3

N3
A3

1.887 10

10.147

You might also like