- Home
- Documents
*Introduction to Finite Element Analysis Using SolidWorks ... Introduction to Finite Element Analysis*

prev

next

out of 22

View

249Download

16

Embed Size (px)

Introduction to Finite Element Analysis Using SolidWorks Simulation 2012Randy H. Shih

SDCP U B L I C A T I O N S

www.SDCpublications.comBetter Textbooks. Lower Prices.

Schroff Development Corporation

Visit the following websites to learn more about this book:

http://www.sdcpublications.com/Textbooks/Introduction-Finite-Element-Analysis-Using/ISBN/978-1-58503-704-9/http://www.amazon.com/gp/product/1585037044?ie=UTF8&tag=sdcpublications&linkCode=as2&camp=211189&creative=374929&creativeASIN=1585037044http://search.barnesandnoble.com/booksearch/isbninquiry.asp?EAN=1585037044http://books.google.com/books?vid=ISBN1585037044&printsec=frontcover

Introduction to Finite Element Analysis 2-1

Chapter 2 Truss Elements in Two-Dimensional Spaces

50 lbs

9 in.

12 in.

Perform 2D Coordinate Transformation. Expand the Direct Stiffness Method to

2D Trusses. Derive the general 2D element Stiffness

Matrix. Assemble the Global Stiffness Matrix for

2D Trusses. Solve 2D trusses using the Direct

Stiffness Method.

2-2 Introduction to Finite Element Analysis

Introduction This chapter presents the formulation of the direct stiffness method of truss elements in a two-dimensional space and the general procedure for solving two-dimensional truss structures using the direct stiffness method. The primary focus of this text is on the aspects of finite element analysis that are more important to the user than the programmer. However, for a user to utilize the software correctly and effectively, some understanding of the element formulation and computational aspects are also important. In this chapter, a two-dimensional truss structure consisting of two truss elements (as shown below) is used to illustrate the solution process of the direct stiffness method.

Truss Elements in Two-Dimensional Spaces As introduced in the previous chapter, the system equations (stiffness matrix) of a truss element can be represented using the system equations of a linear spring in one-dimensional space.

Free Body Diagram:

K

F

K = EA/L

F

K

F1 F2 Node 1 Node 2

+X1 +X2

Truss Elements in Two-Dimensional Spaces 2-3

The general force-displacement equations in matrix form:

F1 + K - K X1 F2 - K + K X2

For a truss element, K = EA/L

F1 + 1 - 1 X1 F2 - 1 + 1 X2

For truss members positioned in two-dimensional space, two coordinate systems are established:

1. The global coordinate system (X and Y axes) chosen to represent the entire structure.

2. The local coordinate system (X and Y axes) selected to align the X axis along the length of the element.

The force-displacement equations expressed in terms of components in the local XY coordinate system:

F1X + 1 - 1 X1 F2X - 1 + 1 X2

=

= EA L

F1

F2

Node 1

Node 2 X1

X2 Y2

Y1

K

= EA L

Y

X

Global X axis

Local X axes

2-4 Introduction to Finite Element Analysis

The above stiffness matrix (system equations in matrix form) can be expanded to incorporate the two force components at each node and the two displacement components at each node.

F1X +1 0 -1 0 X1 F1Y 0 0 0 0 Y1

F2X -1 0 +1 0 X2 F2Y 0 0 0 0 Y2 In regard to the expanded local stiffness matrix (system equations in matrix form):

1. It is always a square matrix. 2. It is always symmetrical for linear systems. 3. The diagonal elements are always positive or zero.

The above stiffness matrix, expressed in terms of the established 2D local coordinate system, represents a single truss element in a two-dimensional space. In a general structure, many elements are involved, and they would be oriented with different angles. The above stiffness matrix is a general form of a SINGLE element in a 2D local coordinate system. Imagine the number of coordinate systems involved for a 20-member structure. For the example that will be illustrated in the following sections, two local coordinate systems (one for each element) are needed for the truss structure shown below. The two local coordinate systems (X1Y1 and X2Y2) are aligned to the elements.

= EA L

Nodal Displacements (Local Coordinate System)

Force Components (Local Coordinate System)

X

Y

Element 1 Element 2

X1

X2

Y2

Y1

A

B

Truss Elements in Two-Dimensional Spaces 2-5

In order to solve the system equations of two-dimensional truss structures, it is necessary to assemble all elements stiffness matrices into a global stiffness matrix, with all the equations of the individual elements referring to a common global coordinate system. This requires the use of coordinate transformation equations applied to system equations for all elements in the structure. For a one-dimensional truss structure (illustrated in chapter 2), the local coordinate system coincides with the global coordinate system; therefore, no coordinate transformation is needed to assemble the global stiffness matrix (the stiffness matrix in the global coordinate system). In the next section, the coordinate transformation equations are derived for truss elements in two-dimensional spaces. Coordinate Transformation A vector, in a two-dimensional space, can be expressed in terms of any coordinate system set of unit vectors. For example,

Vector A can be expressed as: A = X i + Y j

Where i and j are unit vectors along the X and Y axes.

Magnitudes of X and Y can also be expressed as:

X = A cos () Y = A sin ()

Where X, Y and A are scalar quantities.

Therefore,

A = X i + Y j = A cos () i + A sin () j ---------- (1)

A

i

j

X

Y

2-6 Introduction to Finite Element Analysis

Next, establish a new unit vector (u) in the same direction as vector A.

Vector A can now be expressed as: A = A u ------------- (2)

Both equations (the above (1) and (2)) represent vector A:

A = A u = A cos () i + A sin () j

The unit vector u can now be expressed in terms of the original set of unit vectors i and j:

u = cos () i + sin () j

Now consider another vector B:

Vector B can be expressed as:

B = - X i + Y j

Where i and j are unit vectors along the X- and Y-axes.

Magnitudes of X and Y can also be expressed as components of the magnitude of the vector:

X = B sin () Y = B cos ()

Where X, Y and B are scalar quantities.

Therefore, B = - X i + Y j = - B sin () i + B cos () j ---------- (3)

A

i

j u

B

i

j

X

Y

Truss Elements in Two-Dimensional Spaces 2-7

Next, establish a new unit vector (v) along vector B.

Vector B can now be expressed as: B = B v ------------- (4)

Equations (3) and (4) represent vector B:

B = B v = - B sin () i + B cos () j

The unit vector v can now be expressed in terms of the original set of unit vectors i and j:

v = - sin () i + cos () j We have established the coordinate transformation equations that can be used to transform vectors from ij coordinates to the rotated uv coordinates. Coordinate Transformation Equations:

u = cos () i + sin () j v = - sin () i + cos () j In matrix form,

u cos () sin () i v - sin () cos () j

B

i

j v

i

j v

u

=

Direction cosines

2-8 Introduction to Finite Element Analysis

The above direction cosines allow us to transform vectors from the GLOBAL coordinates to the LOCAL coordinates. It is also necessary to be able to transform vectors from the LOCAL coordinates to the GLOBAL coordinates. Although it is possible to derive the LOCAL to GLOBAL transformation equations in a similar manner as demonstrated for the above equations, the MATRIX operations provide a slightly more elegant approach. The above equations can be represented symbolically as: {a} = [ l ] {b} where {a} and {b} are direction vectors, [l] is the direction cosines. Perform the matrix operations to derive the reverse transformation equations in terms of the above direction cosines:

{b} = [ ? ] {a} First, multiply by [ l ] 1 to remove the direction cosines from the right hand side of the original equation. {a} = [ l ] {b} [ l ] 1{a} = [ l ] 1 [ l ] {b} From matrix algebra, [ l ] 1 [ l ]= [ I ] and [ I ]{b} = {b} The equation can now be simplified as [ l ] 1{a} = {b} For linear statics analyses, the direction cosines is an orthogonal matrix and the inverse of the matrix is equal to the transpose of the matrix. [ l ] 1 = [ l ] T Therefore, the transformation equation can be expressed as: [ l ] T {a} = {b}

Truss Elements in Two-Dimensional Spaces 2-9

The transformation equations that enable us to transform any vector from a LOCAL coordinate system to the GLOBAL coordinate system become: LOCAL coordinates to the GLOBAL coordinates:

i cos () - sin ()