1-Unknowns:: A - The System Stiffness Matrix

Download as pdf or txt
Download as pdf or txt
You are on page 1of 6

Benha University Department: Civil Engineering

Benha Faculty of Engineering Program: Structural

Semester: 1 (2020-2021) Subject: Structural Analysis


Examiner:Prof. Ashraf Abourayan Code: CE1311
Total Marks: 75

Question 4:
𝑳𝟏 = 𝑳𝟐 = 𝟒𝒎
𝑬𝑰 = 𝟏𝟎𝟔 𝑲𝑵/𝒎𝟐 , 𝒍𝒆𝒕 𝑬𝑨 = 𝟐𝟎𝟎𝟎𝟎𝑲𝑵
𝑬𝑨 𝟒
→ 𝜶&𝜷 𝒊𝒔 𝒄𝒐𝒏𝒔𝒕, 𝜷 = 𝑳 = 𝟎. 𝟑𝟐
𝑬𝑰
a- The system stiffness matrix
1- Unknowns:
𝑛 = 3(3) − 6 = 3
0
𝑄̅ = [0]
0

2- Element Model :
𝑓 𝑖 = 𝑘. 𝑑𝑖 ,
𝛽 0 0 −𝛽 0 0
0 12 6𝐿 0 −12 6𝐿
0 6𝐿 4𝐿2 0 −6𝐿 2𝐿2 𝐸𝐼 𝐴𝐿2
𝑘= 𝛼 𝛼 = 3 ,𝛽 =
−𝛽 0 0 𝛽 0 0 𝐿 𝐼
0 −12 −6𝐿 0 12 −6𝐿
[ 0 6𝐿 2𝐿2 0 −6𝐿 4𝐿2 ]
Benha University Department: Civil Engineering
Benha Faculty of Engineering Program: Structural

Semester: 1 (2020-2021) Subject: Structural Analysis


Examiner:Prof. Ashraf Abourayan Code: CE1311
Total Marks: 75

3- System Model:
- M-code:

0 1
0 2
0 3
𝑀=
1 0
2 0
[3 0]

- 𝐸𝑙𝑒𝑚𝑒𝑚𝑡 1 ∶ 𝜃 = 90, 𝑐1 = 0, 𝑐2 = 1

𝑔1 𝑔2 −𝑔4 12 0 6𝐿
𝑘 (1) = 𝛼 [ 𝑔2 𝑔3 −𝑔5 ] = 𝛼 [ 0 𝛽 0 ]
−𝑔4 −𝑔5 𝑔6 6𝐿 0 4𝐿2

- 𝐸𝑙𝑒𝑚𝑒𝑚𝑡 2 ∶ 𝜃 = 0, 𝑐1 = 1, 𝑐2 = 0

𝑔1 𝑔2 𝑔4 𝛽 0 0
- 𝑘 (2)
= 𝛼 [𝑔2 𝑔3 𝑔5 ] = 𝛼 [0 12 6𝐿 ]
𝑔4 𝑔5 𝑔6 0 6𝐿 4𝐿2

2 𝛽 + 12 0 6𝐿
𝑘 𝑠 = ∑ 𝑘𝑖 = 𝛼 [ 0 12 + 𝛽 6𝐿 ]
𝑖=1 6𝐿 6𝐿 8𝐿2
Benha University Department: Civil Engineering
Benha Faculty of Engineering Program: Structural

Semester: 1 (2020-2021) Subject: Structural Analysis


Examiner:Prof. Ashraf Abourayan Code: CE1311
Total Marks: 75

b- The system model equation:

0 −1 0 0 0 0 0 −10 0
1 0 0 0 0 0 10 0 0
0 0 1 0 0 0 6.67 6.67 0
𝐹1̂ = Ԉ1𝑇 𝑓 1̂ = =
0 0 0 0 −1 0 0 −10 1
0 0 0 1 0 0 10 0 2
[0 0 0 0 0 1] [−6.67] [−6.67] 3
−10
̂)
(1
⇒𝐹 =[ 0 ]
−6.67
0 1
10 2
0
̂
2 ̂
2 20 3 ̂)
(2
𝐹 =𝑓 = ⇒ 𝐹 = [10]
0 0
20
10 0
[−20] 0

−10 0 −10
𝑄̂ = [ 0 ] + [10] = [ 10 ]
−6.67 20 13.33

0 −10 10
𝑄 = 𝑄̅ − 𝑄̂ = [0] − [ 10 ] = [ −10 ]
0 13.33 −13.33
Benha University Department: Civil Engineering
Benha Faculty of Engineering Program: Structural

Semester: 1 (2020-2021) Subject: Structural Analysis


Examiner:Prof. Ashraf Abourayan Code: CE1311
Total Marks: 75

4- 𝑆𝑜𝑙𝑢𝑡𝑖𝑜𝑛

10 𝛽 + 12 0 6𝐿 𝑞1
[ −10 ] = 𝛼 [ 0 12 + 𝛽 6𝐿 ] [𝑞2 ]
−13.33 6𝐿 6𝐿 8𝐿2 𝑞3

1.5647 −0.0587 −0.3865


𝑞1 = , 𝑞2 = , 𝑞3 =
𝛼 𝛼 𝛼

5- Element model
Element1 :
𝐷11
1
0 0
𝐷2 0 0
1
𝐷 0 0
𝐷1 = 3 1 ⇒ → 𝑑1 = 𝛬1 𝐷1
𝐷4 1 𝑞1

𝐷51 2 𝑞2
3 [𝑞3 ]
[𝐷61 ]
0 1 0 0 0 0 0 0
−1 0 0 0 0 0 0 0
0 0 1 0 0 0 0 0
𝑑1 = =
0 0 0 0 1 0 𝑞1 𝑞2
0 0 0 −1 0 0 𝑞2 −𝑞1
[ 0 0 0 0 0 1] [𝑞3 ] [ 𝑞3 ]
Benha University Department: Civil Engineering
Benha Faculty of Engineering Program: Structural

Semester: 1 (2020-2021) Subject: Structural Analysis


Examiner:Prof. Ashraf Abourayan Code: CE1311
Total Marks: 75

𝛽 0 0 −𝛽 0 0 0
0 12 6𝐿 0 −12 6𝐿 0
1 1 1 0 6𝐿 4𝐿2 0 −6𝐿 2𝐿2 0
𝑓 = 𝑘 𝑑 =𝛼
−𝛽 0 0 𝛽 0 0 𝑞2
0 −12 −6𝐿 0 12 −6𝐿 −𝑞1
[0 6𝐿 2𝐿 2
0 −6𝐿 4𝐿2 ] [ 𝑞3 ]
0.02
9.50
25.18
=
−0.02
−9.50
[ 12.82 ]
0 0.02 0.02
10 9.50 19.50
6.67 25.18 31.85
𝑓 1̅ = 𝑓 1̂ + 𝑓 2 = + =
0 −0.02 −0.02
10 −9.50 0.50
[−6.67] [ 12.82 ] [ 6.15 ]

Element 2:
𝐷1 2 𝑞1 𝑞1
1
𝐷2 2 2 𝑞2 𝑞2
2
𝐷 3 𝑞3 𝑞
𝐷2 = 3 2 ⇒ → 𝑑2 = 𝐷2 = 3
𝐷4 0 0 0
𝐷5 2 0 0 0
2 0
[0] [0]
[𝐷6 ]
Benha University Department: Civil Engineering
Benha Faculty of Engineering Program: Structural

Semester: 1 (2020-2021) Subject: Structural Analysis


Examiner:Prof. Ashraf Abourayan Code: CE1311
Total Marks: 75

𝛽 0 0 −𝛽 0 0 𝑞1
0 12 6𝐿 0 −12 6𝐿 𝑞2
0 6𝐿 4𝐿 2
0 −6𝐿 2𝐿 2
𝑞3
𝑓 2 = 𝑘 2 𝑑2 = 𝛼
−𝛽 0 0 𝛽 0 0 0
0 −12 −6𝐿 0 12 −6𝐿 0
2] [
[ 0 6𝐿 2𝐿 2
0 −6𝐿 4𝐿 0]
0.50
−10.0
−26.14
=
−0.50
10.00
[−13.78]
0 0.50 0.50
10 −10.0 0
20 −26.14 −6.14
𝑓 2̅ = 𝑓 2̂ + 𝑓 2 = + =
0 −0.50 −0.50
10 10.00 20
[−20] [−13.78] [ 6.27 ]

6- Joint Forces:
0 −1 0 0.02 19.50
̅ 1 1𝑇 ̅ 1
𝑝1 = 𝐹𝑎 = ƛ 𝑓𝑎 = [1 0 0] [19.50] = [ 0.02 ]
0 0 1 31.85 31.85

0 −1 0 −0.02 0.5 0
̅ 1
̅ 2 1𝑇 ̅ 1 ̅ 2
𝑝2 = 𝐹𝑏 + 𝐹𝑎 = ƛ 𝑓𝑏 + 𝑓𝑎 = [1 0 0] [ 0.50 ] + [ 0 ] ≅ [0]
0 0 1 6.15 −6.14 0

2
−0.5
2
̅ ̅
𝑝3 = 𝐹𝑏 = 𝑓𝑏 = [ 20 ]
6.27

You might also like