E=210*10^6;
I=1*10^-4;
A=1*10^-2;
k 1=
10^5 *
[3.5000 0 0 -3.5000 0 0
0 0.0117 0.0350 0 -0.0117 0.0350
0 0.0350 0.1400 0 -0.0350 0.0700
-3.5000 0 0 3.5000 0 0
0 -0.0117 -0.0350 0 0.0117 -0.0350
0 0.0350 0.0700 0 -0.0350 0.1400]
T 1=
0 1 0 0 0 0
-1 0 0 0 0 0
0 0 1 0 0 0
0 0 0 0 1 0
0 0 0 -1 0 0
0 0 0 0 0 1
Kele1 =
10^5 *
0.0117 0 -0.0350 -0.0117 0 -0.0350
0 3.5000 0 0 -3.5000 0
-0.0350 0 0.1400 0.0350 0 0.0700
-0.0117 0 0.0350 0.0117 0 0.0350
0 -3.5000 0 0 3.5000 0
-0.0350 0 0.0700 0.0350 0 0.1400
k2(local) =
10^5 *
3.7123 0 0 -3.7123 0 0
0 0.0139 0.0394 0 -0.0139 0.0394
0 0.0394 0.1485 0 -0.0394 0.0742
-3.7123 0 0 3.7123 0 0
0 -0.0139 -0.0394 0 0.0139 -0.0394
0 0.0394 0.0742 0 -0.0394 0.1485
T 2=
0.7071 0.7071 0 0 0 0
-0.7071 0.7071 0 0 0 0
0 0 1.0000 0 0 0
0 0 0 0.7071 0.7071 0
0 0 0 -0.7071 0.7071 0
0 0 0 0 0 1.0000
Kele(2) =
10^5 *
1.8631 1.8492 -0.0278 -1.8631 -1.8492 -0.0278
1.8492 1.8631 0.0278 -1.8492 -1.8631 0.0278
-0.0278 0.0278 0.1485 0.0278 -0.0278 0.0742
-1.8631 -1.8492 0.0278 1.8631 1.8492 0.0278
-1.8492 -1.8631 -0.0278 1.8492 1.8631 -0.0278
-0.0278 0.0278 0.0742 0.0278 -0.0278 0.1485
k 3=
10^5 *
3.7123 0 0 -3.7123 0 0
0 0.0139 0.0394 0 -0.0139 0.0394
0 0.0394 0.1485 0 -0.0394 0.0742
-3.7123 0 0 3.7123 0 0
0 -0.0139 -0.0394 0 0.0139 -0.0394
0 0.0394 0.0742 0 -0.0394 0.1485
T3 =
0.7071 -0.7071 0 0 0 0
0.7071 0.7071 0 0 0 0
0 0 1.0000 0 0 0
0 0 0 0.7071 -0.7071 0
0 0 0 0.7071 0.7071 0
0 0 0 0 0 1.0000
Kele(3) =
10^5 *
1.8631 -1.8492 0.0278 -1.8631 1.8492 0.0278
-1.8492 1.8631 0.0278 1.8492 -1.8631 0.0278
0.0278 0.0278 0.1485 -0.0278 -0.0278 0.0742
-1.8631 1.8492 -0.0278 1.8631 -1.8492 -0.0278
1.8492 -1.8631 -0.0278 -1.8492 1.8631 -0.0278
0.0278 0.0278 0.0742 -0.0278 -0.0278 0.1485
k 4=
10^5 *
3.5000 0 0 -3.5000 0 0
0 0.0117 0.0350 0 -0.0117 0.0350
0 0.0350 0.1400 0 -0.0350 0.0700
-3.5000 0 0 3.5000 0 0
0 -0.0117 -0.0350 0 0.0117 -0.0350
0 0.0350 0.0700 0 -0.0350 0.1400
T 4=
0 1 0 0 0 0
-1 0 0 0 0 0
0 0 1 0 0 0
0 0 0 0 1 0
0 0 0 -1 0 0
0 0 0 0 0 1
Kele(4) =
10^5*
0.0117 0 -0.0350 -0.0117 0 -0.0350
0 3.5000 0 0 -3.5000 0
-0.0350 0 0.1400 0.0350 0 0.0700
-0.0117 0 0.0350 0.0117 0 0.0350
0 -3.5000 0 0 3.5000 0
-0.0350 0 0.0700 0.0350 0 0.1400
Kglobal=10^5 *
0.0117 0 -0.0350 -0.0117 0 -0.0350 0 0 0 0 0 0 0 0 0
0 3.5000 0 0 -3.5000 0 0 0 0 0 0 0 0 0 0
-0.0350 0 0.1400 0.0350 0 0.0700 0 0 0 0 0 0 0 0 0
-0.0117 0 0.0350 1.8748 1.8492 0.0072 -1.8631 -1.8492 -0.0278 0 0 0 0 0 0
0 -3.5000 0 1.8492 5.3631 0.0278 -1.8492 -1.8631 0.0278 0 0 0 0 0 0
-0.0350 0 0.0700 0.0072 0.0278 0.2885 0.0278 -0.0278 0.0742 0 0 0 0 0 0
0 0 0 -1.8631 -1.8492 0.0278 3.7262 0 0.0557 -1.8631 1.8492 0.0278 0 0 0
0 0 0 -1.8492 -1.8631 -0.0278 0 3.7262 0 1.8492 -1.8631 0.0278 0 0 0
0 0 0 -0.0278 0.0278 0.0742 0.0557 0 0.2970 -0.0278 -0.0278 0.0742 0 0 0
0 0 0 0 0 0 -1.8631 1.8492 -0.0278 1.8748 -1.8492 0.0072 -0.0117 0 0.0350
0 0 0 0 0 0 1.8492 -1.8631 -0.0278 -1.8492 5.3631 -0.0278 0 -3.5000 0
0 0 0 0 0 0 0.0278 0.0278 0.0742 0.0072 -0.0278 0.2885 -0.0350 0 0.0700
0 0 0 0 0 0 0 0 0 -0.0117 0 -0.0350 0.0117 0 -0.0350
0 0 0 0 0 0 0 0 0 0 -3.5000 0 0 3.5000 0
0 0 0 0 0 0 0 0 0 0.0350 0 0.0700 -0.0350 0 0.1400
F+Ff=Kglobal*d
F=Nodal forces
Ff=lump vector
d=nodal displacement
d=
0
0
0
-0.0608
-0.0002
0.0078
-0.0530
-0.0083
-0.0075
-0.0449
-0.0001
0.0117
0
0
0
F =
43.6514
72.0911
-158.1780
5.0000
0.0000
0.0000
0.0000
-60.0000
0.0000
-0.0000
0.0000
0.0000
11.3486
47.9089
-75.0930
Consider the frame in Fig. 1, the node and element numbers as well as the material and geometrica...
Analyse the beam shown in Figure 4 using the stiffiness method. Node D is fixed and node 2 and 3 are rollers. A uniform distributed load of 1 kN/m is acting on member 1 . And a load of 10 kN is acting at the middle of member2. EI is constant for all members a) Identify the force vector of the structure; [4 marks] b) Identify the displacement vector of the structure; [2 marks] c) Determine the stiffness matrices of...
Consider a three-node quadratic element in one dimension with unequally spaced nodes (Fig. 2). (a) Obtain the B matrix. (b) Consider an element with x1=0, x2=1/5 and x3=1 Evaluate strain & in terms of u2 and u (u0), and check what happens when s approaches 0. (c) If you evaluate K by one-point quadrature using B EADB for same coordinates as in (b) and constrain node 1 (i.e. 0), is K invertible?
The plane frame in fig. 1 has columns of length h pinned at the base. The beam is welded to the columns and is of length 2l. Equal loads P are applied at the mid-span of both the beam and the right column. Find the collapse load PC in terms of MP , where MP is the plastic moment of the columns and beam. l is 3 m and h is 5m. The plane frame in fig. 1 has...
MEE456 Homework Consider the beam shown below: 550 Analyze this beam using a 3-node, 2-element model (hand calculations or E-200 G Pa and i=1 x 10" m'. Report on a summary page: a) Sketch a diagram of the beam showing the nodes, elements and loading (FE model). b) Present the displacements, rotations, reactions and, element forces and moments. c) Sketch a free body diagram of the beam showing the original loading and the reactions. d) Sketch a free body diagram...
Element 1 is a steel bar that has a circular cross-section with a radius of 30 mm. Element 2 is an aluminum bar that has a circular cross-section with a radius of 50 mm. Element 3 is a steel bar that has a circular cross- section with a radius of 60 mm. Assume for steel, the moduļus, E, is 2.0E11 Pa, and the density, p, is 7800 kg/m3. Assume for aluminum, E-7.0E10 Pa, and p-2700 kg/m3. The rigid, massless rod...
Section 1: Finite Element Derivation and Validation In this section of the report you will develop your own Finite Element method for 1-dimensional axial loading. The governing equation for displacement, u is Poisson's Equation: อั1 where E is the modulus of elasticity, A(a) is the cross-sectional area as a function of length, and q(x) is the loading distribution as a function of length. The weak form of this equation with 0 1. Starting from the weak form of the governing...
a. Draw a complete free-body diagram for the beam shown in Fig. 1. Must show and label x-y coordinate system, support reactions, loads, and dimensions b. Derive formulas for all support reactions using conditions of equilibrium. Must show all work. c. Using the student variables shown below, calculate beam support reactions and draw shear and bending moment diagrams. To receive full credit, support reactions and loads must be shown on a free-body diagram. Each shear and bending moment diagram must...
What is the reaction Ax to the nearest 0.01 kN for the following frame? El are constant and the same for all members. Consider only bending energy. W = 67 kN/m w -5 m-+ 10 m ---5 m 5 m А B 422 CHAPTER 10 ANALYSIS OF STATICALLY INDETERMINATE STRUCTURES BY THE FORCE METHOD 10.5 Force Method of Analysis: Frames The force method is very useful for solving problems involving statically indeterminate frames that have a single story and unusual...
QUESTION 1 [25 marks A frame loaded with a uniformly distributed load at Member AB and point load at Member BC and joint B. It has pinned supports A and C, while joint B is fixed connected, as can be seen in Figure 1. Take E-200 GPa. a) Using the slope-deflection method, calculate the moments and illustrate the bending moment diagram. [15 marks) b) Then calculate the shear forces and sketch the shear force diagram. [10 marks) 22 KN 10...
1. The cantilever beam AB of length L shown in Fig.(a carries a uniformly distributed load of intensity wo, which includes the weight of the beam. Equation of the elastic curve is given by the following differential equation Af _irer where M is the bending moment of the beam. it is also given that of the beam and I modulus of Inertia. where E- Modulus of the elasticity v is the deflection of the beam. Compute the maximum displacement of...