Trusses 192022 1 Element Formulation by Virtual Work
Trusses 1/9/2022 1
Element Formulation by Virtual Work u Use virtual work to derive element stiffness matrix based on assumed displacements – Principle of virtual work states that if a general structure that is in equilibrium with its applied forces deforms due to a set of small compatible virtual displacements, the virtual work done is equal to its virtual strain energy of internal stresses. 1/9/2022 2
u At element level, d. Ue = d. We – d. Ue = virtual strain energy of internal stresses – d. We = virtual work of external forces acting through virtual displacements 1/9/2022 3
u We now assume a simple displacement function to define the displacement of every material point in the element. u Usually use low order polynomials u Here u = a 1 + a 2 x – u is axial displacement – a 1, a 2 are constants to be determined – x is local coordinate along member 1/9/2022 4
u The constants are found by imposing the known nodal displacements ui, uj at nodes i and j ui = a 1 + a 2 xi uj = a 1 + a 2 xj u ui, uj are nodal displacements u xi, xj are nodal coordinates 1/9/2022 5
u letting xi = 0, xj = L, we get – a 1 = ui – a 2 = (uj-ui)/L u We can write – [N] = matrix of element shape functions or interpolation functions – {d} = nodal displacements 1/9/2022 6
N 1=1 Variation of N 1 N 2=1 Variation of N 2 1/9/2022 7
u Strain is given by u where [B] is a matrix relating strain to nodal displacement (matrix of derivatives of shape function) 1/9/2022 8
u Now s = E(e- eo )= E[B]{d}-E eo – Stress and strain are constant in a member u Define internal virtual strain energy for a set of virtual displacements {dd} to be 1/9/2022 9
- de = virtual strain - s = stress level at equilibrium – d. V = volume u Virtual work of nodal forces is T d. We = {dd} {f} u Then, virtual work is given by 1/9/2022 10
u Substituting u Canceling and rearranging gives T {dd} gives [k]{d}={F} where For thermal problem 1/9/2022 11
u for a truss we get u this formulation method also applies to 2 -d and 3 -d elements 1/9/2022 12
Procedure for Direct Stiffness Method (Displacement Method) 1. Discretize into finite elements, Identify nodes, elements and number them in order. 2. Develop element stiffness matrices [Ke] for all the elements. 3. Assemble element stiffness matrices to get the global stiffness matrix ([KG] =S [Ke]). The size of of global stiffness matrix = total d. o. f of the structure including at boundary nodes. Assembly is done by matching element displacement with global displacements. Also develop appropriate force vector (by adding element force vectors) such that equation of the type [KG] {u}={F} is obtained. 1/9/2022 13
Procedure for Direct Stiffness Method 4. Apply kinematic boundary conditions. Without applying boundary conditions, [KG] will be singular. (minimum number of boundary conditions required is to arrest ‘Rigid Body’ displacements). 5. Solve for unknown displacements {u} ( {u}= [KG] – 1{F}). 6. Once displacements are determined find (a) reactions by picking up appropriate rows from the equation {F}=[KG] {u}, (b) Find element forces {f}=[Ke] {ue}, (c) Element stresses given by {se}= [D][B]{ue}. 1/9/2022 14
1/9/2022 15
2 A, L, E 1/9/2022 16
1/9/2022 17
Reactions { 1/9/2022 18
Element Forces f 2 f 1 A, L, E 2 A, L, E 1/9/2022 19
1/9/2022 20
1/9/2022 21
1/9/2022 22
u 1 u 2 1/9/2022 u 3 23
1/9/2022 24
1/9/2022 25
1/9/2022 26
1/9/2022 27
Direct Element Formulation u truss element acts like 1 -d spring – l >> transverse dimensions – pinned connection to other members (only axial loading). – usually constant cross section and modulus of elasticity 1/9/2022 28
» A = cross section area » E = modulus of elasticity » L = length 1/9/2022 29
u Assume displacements are much smaller than overall geometry – vertical displacements of horizontal member produce no vertical force u Stiffness matrix is written in local element coordinates aligned along element axis u want stiffness matrix for arbitrary orientation 1/9/2022 30
u rotate coordinate systems using rotation matrix [R] u displacement components in global coordinates are related to displacement components in local coordinates by {d’}=[R]{d} – {d} = displacement in global coordinates – {d’} = displacement in local element coordinates 1/9/2022 31
1/9/2022 32
1/9/2022 33
u start with member on x axis, element equations are or {k’}{d’}={f’} u Note that y equations are all zero 1/9/2022 34
1/9/2022 35
At node i A similar matrix can be obtained at node j 1/9/2022 36
u Matrix [R] is: 1/9/2022 37
u Similarly , force components are related by {f’} = [R]{f} u Local force displacement relation is [k’]{d’} = {f’} u global force displacement relation is [k][R]{d} = [R]{f} -1 T u using fact that [R] = [R] , we get T [R] [k][R]{d} = {f} 1/9/2022 38
u then [k] = stiffness matrix in global T coordinates is [R] [k’][R] 1/9/2022 39
u Structure equation is [k] {D} = {F} – [k] = structure stiffness matrix – {D} = nodal displacement vector – {F} = applied load vector 1/9/2022 40
1/9/2022 41
2 E le m e n t C o o r d in a t e i- n o d e y x 1 1 2 0 3 0 C o o r d in a te 0 2 L s in 4 5 L e n g th j- n o d e 2 2 x L c o s 4 5 L s in 4 5 1/9/2022 C y s L c o s 4 5 s in 4 5 L c os 4 5 -sin 4 5 42
1/9/2022 43
1/9/2022 44
1/9/2022 45
1/9/2022 46
1/9/2022 47
1/9/2022 48
1/9/2022 49
1/9/2022 50
v’ 3 1/9/2022 v 3 u’ 3 q u 3 51
1/9/2022 52
1/9/2022 53
1/9/2022 54
Finite Element Model u usually use existing codes to solve problems u user responsible for – creating the model – executing the program – interpreting the results 1/9/2022 55
u arrangement of nodes and elements is known as the mesh u plan to make the mesh model the structure as accurately as possible 1/9/2022 56
u for a truss – each member is modeled as 1 truss element – truss members or elements are connected at nodes – node connections behave like pin joints – truss element behaves in exact agreement with assumptions – no need to divide a member into more than 1 element 1/9/2022 57
– such subdivision will cause execution to fail » due to zero stiffness against lateral force at the node connection where 2 members are in axial alignment 1/9/2022 58
1/9/2022 59
u there is geometric symmetry – often possible to reduce the size of problem by using symmetry – need loading symmetry as well 1/9/2022 60
1/9/2022 61
u Fig. 3 -5 and 3 -6 show symmetric loads and the reduced model – need to impose extra conditions along the line of symmetry » displacement constraints: nodes along the line of symmetry must always move along that line » changed loads: the load at the line of symmetry is split in two 1/9/2022 62
Computer input assistance ua preprocessor is used to assist user input u required inputs are – data to locate nodes in space – definition of elements by node numbers – type of analysis to be done – material properties – displacement conditions – applied loads 1/9/2022 63
u interactive preprocessors are preferable – you can see each node as it is created – elements are displayed as they are created – symbols are given for displacement and load conditions – usually allow mesh generation by replication or interpolation of an existing mesh – allow inserting nodes along lines – allow entering a grid by minimum and maximum positions plus a grid spacing 1/9/2022 64
u truss element consists of 2 node numbers that connect to form element u other information for truss is – modulus of elasticity – cross sectional area u data can form a material table u assign element data by reference to the table 1/9/2022 65
u boundary or displacement conditions are set by selecting a node and setting its displacement U do not over constrain a structure by prescribing zero displacements where there is no physical support 1/9/2022 66
u loading conditions are set by selecting nodes and specifying force or moment components U check model carefully at this point 1/9/2022 67
Analysis Step u mostly transparent to user u small truss models have enough accuracy and performance for an accurate solution u a large model has a large number of elements and nodes 1/9/2022 68
u numerical solution may not be accurate if there are full matrices u get better accuracy if the nonzero terms are close to the diagonal – reduces the number of operations and round off error (banded matrix) 1/9/2022 69
u in FE model, element or node numbering can affect bandwidth – good numbering pattern can minimize bandwidth – different methods based on node or element numbering – to minimize, plan numbering pattern so nodes that connect through an element have their equations assembled close together 1/9/2022 70
u In Fig. 3 -7, node numbers are considered, X’s show nonzero terms 1/9/2022 71
u In Fig. 3 -8, node numbers are considered 1/9/2022 72
u many programs have bandwidth or wavefront minimizers available u most programs will keep original numbering for display but use the minimized number scheme 1/9/2022 73
u numerical algorithms, numerical range of the computer affect solution u relative stiffness of members can influence results – problems when members of high and low stiffness connect – can exceed precision of computer – physical situation is usually undesirable 1/9/2022 74
u Approximation error for truss is zero u Most common error messages (errors) come from – incorrect definition of elements – incorrect application of displacement boundary conditions 1/9/2022 75
– may get non-positive definite structure stiffness matrix from not enough boundary conditions to prevent rigid body motion » two elements connect in-line Þ zero lateral stiffness » truss structure not kinematically stable (linkage) 1/9/2022 76
u next look at stress components – in continua, stress components are related to averaged quantities at the nodes – trusses have a stress in each member (not easy to plot) u truss model is exact so it does not usually need refinement 1/9/2022 77
Output Processing and Evaluation u Get numerical results with input data followed by all nodal displacements and element stresses u first graphic to look at is the deformed shape of the structure – nodal displacements are exaggerated to show structure deformation – check to ensure model behaves as expected 1/9/2022 78
u linear elastic analysis, failure is by – overstressing – buckling (have to find members with significant compression and use Euler's buckling equation) 1/9/2022 79
Final Remarks u few situations where a truss element is the right element for modeling behavior 1/9/2022 80
- Slides: 80