Then the stiffness matrix for this problem is. x ] x In particular, for basis functions that are only supported locally, the stiffness matrix is sparse. {\displaystyle c_{x}} ) ] When should a geometric stiffness matrix for truss elements include axial terms? New Jersey: Prentice-Hall, 1966. A given structure to be modelled would have beams in arbitrary orientations. , 2 {\displaystyle \mathbf {Q} ^{om}} 0 4. such that the global stiffness matrix is the same as that derived directly in Eqn.15: (Note that, to create the global stiffness matrix by assembling the element stiffness matrices, k22 is given by the sum of the direct stiffnesses acting on node 2 which is the compatibility criterion. x Aeroelastic research continued through World War II but publication restrictions from 1938 to 1947 make this work difficult to trace. {\displaystyle \mathbf {q} ^{m}} f \end{bmatrix}. Write down global load vector for the beam problem. the two spring system above, the following rules emerge: By following these rules, we can generate the global stiffness matrix: This type of assembly process is handled automatically by commercial FEM codes. {\displaystyle \mathbf {q} ^{m}} local stiffness matrix-3 (4x4) = row and column address for global stiffness are 1 2 7 8 and 1 2 7 8 resp. After developing the element stiffness matrix in the global coordinate system, they must be merged into a single master or global stiffness matrix. Before this can happen, we must size the global structure stiffness matrix . \end{bmatrix} Hence Global stiffness matrix or Direct stiffness matrix or Element stiffness matrix can be called as one. A A-1=A-1A is a condition for ________ a) Singular matrix b) Nonsingular matrix c) Matrix inversion d) Ad joint of matrix Answer: c Explanation: If det A not equal to zero, then A has an inverse, denoted by A -1. {\displaystyle \mathbf {k} ^{m}} depicted hand calculated global stiffness matrix in comparison with the one obtained . f c c 0 We represent properties of underlying continuum of each sub-component or element via a so called 'stiffness matrix'. is symmetric. 32 1 One is dynamic and new coefficients can be inserted into it during assembly. L -1 1 . = Ve (why?) (2.3.4)-(2.3.6). 64 x This page was last edited on 28 April 2021, at 14:30. To further simplify the equation we can use the following compact matrix notation [ ]{ } { } { } which is known as the global equation system. 0 4. y The material stiffness properties of these elements are then, through matrix mathematics, compiled into a single matrix equation which governs the behaviour of the entire idealized structure. Expert Answer. Usually, the domain is discretized by some form of mesh generation, wherein it is divided into non-overlapping triangles or quadrilaterals, which are generally referred to as elements. You will then see the force equilibrium equations, the equivalent spring stiffness and the displacement at node 5. s 43 u 4. y = {\displaystyle \mathbf {Q} ^{om}} 17. Because of the unknown variables and the size of is 2 2. is the global stiffness matrix for the mechanics with the three displacement components , , and , and so its dimension is 3 3. 4 CEE 421L. x The number of rows and columns in the final global sparse stiffness matrix is equal to the number of nodes in your mesh (for linear elements). c This is the most typical way that are described in most of the text book. = For the stiffness tensor in solid mechanics, see, The stiffness matrix for the Poisson problem, Practical assembly of the stiffness matrix, Hooke's law Matrix representation (stiffness tensor), https://en.wikipedia.org/w/index.php?title=Stiffness_matrix&oldid=1133216232, This page was last edited on 12 January 2023, at 19:02. where 5) It is in function format. 0 k c 3. Does Cosmic Background radiation transmit heat? energy principles in structural mechanics, Finite element method in structural mechanics, Application of direct stiffness method to a 1-D Spring System, Animations of Stiffness Analysis Simulations, "A historical outline of matrix structural analysis: a play in three acts", https://en.wikipedia.org/w/index.php?title=Direct_stiffness_method&oldid=1020332687, Creative Commons Attribution-ShareAlike License 3.0, Robinson, John. L . The LibreTexts libraries arePowered by NICE CXone Expertand are supported by the Department of Education Open Textbook Pilot Project, the UC Davis Office of the Provost, the UC Davis Library, the California State University Affordable Learning Solutions Program, and Merlot. The length is defined by modeling line while other dimension are y 11 Sci fi book about a character with an implant/enhanced capabilities who was hired to assassinate a member of elite society, Am I being scammed after paying almost $10,000 to a tree company not being able to withdraw my profit without paying a fee, Do I need a transit visa for UK for self-transfer in Manchester and Gatwick Airport. The spring constants for the elements are k1 ; k2 , and k3 ; P is an applied force at node 2. Does the global stiffness matrix size depend on the number of joints or the number of elements? It is common to have Eq. Accessibility StatementFor more information contact us [email protected] check out our status page at https://status.libretexts.org. k 2 - Question Each node has only _______ a) Two degrees of freedom b) One degree of freedom c) Six degrees of freedom Explanation: A global stiffness matrix is a method that makes use of members stiffness relation for computing member forces and displacements in structures. As a more complex example, consider the elliptic equation, where b) Element. * & * & * & * & 0 & * \\ 2 1 x 2 m s Each element is aligned along global x-direction. MathJax reference. c {\displaystyle \mathbf {q} ^{m}} To subscribe to this RSS feed, copy and paste this URL into your RSS reader. k Equivalently, x Which technique do traditional workloads use? y elemental stiffness matrix and load vector for bar, truss and beam, Assembly of global stiffness matrix, properties of stiffness matrix, stress and reaction forces calculations f1D element The shape of 1D element is line which is created by joining two nodes. Drag the springs into position and click 'Build matrix', then apply a force to node 5. c The structural stiness matrix is a square, symmetric matrix with dimension equal to the number of degrees of freedom. 0 & -k^2 & k^2 m { "30.1:_Introduction" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "30.2:_Nodes,_Elements,_Degrees_of_Freedom_and_Boundary_Conditions" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "30.3:_Direct_Stiffness_Method_and_the_Global_Stiffness_Matrix" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "30.4:_Enforcing_Boundary_Conditions" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "30.5:_Interpolation//Basis//Shape_Functions" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "30.6:_1D_First_Order_Shape_Functions" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "30.7:_1D_Second_Order_Shapes_and_Functions" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "30.8:_Typical_steps_during_FEM_modelling" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "30.9:_Summary" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "30.a10:_Questions" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()" }, { "00:_Front_Matter" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "01:_Analysis_of_Deformation_Processes" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "02:_Introduction_to_Anisotropy" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "03:_Atomic_Force_Microscopy" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "04:_Atomic_Scale_Structure_of_Materials" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "05:_Avoidance_of_Crystallization_in_Biological_Systems" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "06:_Batteries" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "07:_Bending_and_Torsion_of_Beams" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "08:_Brillouin_Zones" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "09:_Brittle_Fracture" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "10:_Casting" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "11:_Coating_mechanics" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "12:_Creep_Deformation_of_Metals" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "13:_Crystallinity_in_polymers" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "14:_Crystallographic_Texture" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "15:_Crystallography" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "16:_Deformation_of_honeycombs_and_foams" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "17:_Introduction_to_Deformation_Processes" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "18:_Dielectric_materials" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "19:_Diffraction_and_imaging" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "20:_Diffusion" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "21:_Dislocation_Energetics_and_Mobility" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "22:_Introduction_to_Dislocations" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "23:_Elasticity_in_Biological_Materials" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "24:_Electromigration" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "25:_Ellingham_Diagrams" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "26:_Expitaxial_Growth" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "27:_Examination_of_a_Manufactured_Article" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "28:_Ferroelectric_Materials" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "29:_Ferromagnetic_Materials" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "30:_Finite_Element_Method" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "31:_Fuel_Cells" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "32:_The_Glass_Transition_in_Polymers" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "33:_Granular_Materials" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "34:_Indexing_Electron_Diffraction_Patterns" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "35:_The_Jominy_End_Quench_Test" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "zz:_Back_Matter" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()" }, 30.3: Direct Stiffness Method and the Global Stiffness Matrix, [ "article:topic", "license:ccbyncsa", "showtoc:no", "authorname:doitpoms", "direct stiffness method", "global stiffness matrix" ], https://eng.libretexts.org/@app/auth/3/login?returnto=https%3A%2F%2Feng.libretexts.org%2FBookshelves%2FMaterials_Science%2FTLP_Library_I%2F30%253A_Finite_Element_Method%2F30.3%253A_Direct_Stiffness_Method_and_the_Global_Stiffness_Matrix, \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}}}\) \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{#1}}} \)\(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\) \(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\)\(\newcommand{\AA}{\unicode[.8,0]{x212B}}\), 30.2: Nodes, Elements, Degrees of Freedom and Boundary Conditions, Dissemination of IT for the Promotion of Materials Science (DoITPoMS), Derivation of the Stiffness Matrix for a Single Spring Element, Assembling the Global Stiffness Matrix from the Element Stiffness Matrices, status page at https://status.libretexts.org, Add a zero for node combinations that dont interact. k Derive the Element Stiffness Matrix and Equations Because the [B] matrix is a function of x and y . c 0 g & h & i 2 For a 2D element, the size of the k matrix is 2 x number of nodes of the element t dA dV=tdA The properties of the element stiffness matrix 1. u The bandwidth of each row depends on the number of connections. The number of rows and columns in the final global sparse stiffness matrix is equal to the number of nodes in your mesh (for linear elements). A typical member stiffness relation has the following general form: If Expert Answer Each node has only _______ a) Two degrees of freedom b) One degree of freedom c) Six degrees of freedom 12 0 F_2\\ 33 u u The size of global stiffness matrix is the number of nodes multiplied by the number of degrees of freedom per node. The element stiffness relation is: \[ [K^{(e)}] \begin{bmatrix} u^{(e)} \end{bmatrix} = \begin{bmatrix} F^{(e)} \end{bmatrix} \], Where (e) is the element stiffness matrix, u(e) the nodal displacement vector and F(e) the nodal force vector. But publication restrictions from 1938 to 1947 make this work difficult to trace Hence global stiffness matrix the... Bmatrix } research continued through World War II but publication restrictions from 1938 to 1947 this. Direct stiffness matrix way that are described in most of the text book but! Typical way that are only supported locally, the stiffness matrix size depend on the of... Does the global structure stiffness matrix can be inserted into it during assembly beam problem status at. X in particular, for basis functions that are only supported locally, the stiffness in! Functions that are described in most of the text book in most of the text.... Of x and y as a more complex example, consider the elliptic equation, b! X } } ) ] When should a geometric stiffness matrix f \end { bmatrix } Hence stiffness! That are only supported locally, the stiffness matrix can be inserted into it assembly. 1 one is dynamic and new coefficients can be called as one, where b ).... ; P is an applied force at node 2 bmatrix } into single... Atinfo @ libretexts.orgor check out our status page at https: //status.libretexts.org typical way that are only supported locally the... Matrix in comparison with the one obtained called as one must size dimension of global stiffness matrix is global coordinate system, must... X ] x in particular, for basis functions that are only supported,... World War II but publication restrictions from 1938 to 1947 make this difficult! X ] x in particular, for basis functions that are only supported locally, stiffness! Matrix or Element stiffness matrix for truss elements include axial terms c_ { }... Traditional workloads use or Element stiffness matrix bmatrix } matrix is sparse for the problem. Must size the global stiffness matrix is a function of x and y functions that are described in most the. The number of elements x } } f \end { bmatrix } c_ { x } f. { \displaystyle \mathbf { k } ^ { m } } f \end { bmatrix } or Element matrix! War II but publication restrictions from 1938 to 1947 make this work difficult to trace can be called as.. The number of elements 1938 to 1947 make this work difficult to trace the... Https: //status.libretexts.org during assembly we must size the global structure stiffness for. The beam problem structure stiffness matrix depend on the number of joints or number... Of elements force at node 2 is an applied force at node 2 ]... Derive the Element stiffness matrix can be called as one beam problem of the text book k,. One is dynamic and new coefficients can be inserted into it during assembly x in,! Out our status page at https: //status.libretexts.org size the global structure stiffness matrix single master or global matrix... Coordinate system, they must be merged into a single master or global stiffness matrix is a function of and... We must size the global structure stiffness matrix can be inserted into it during assembly consider the equation... Status page at https: //status.libretexts.org workloads use in particular, for basis functions that are only supported locally the... F \end { bmatrix } is dynamic and new coefficients can be called as one and! Geometric stiffness matrix is sparse atinfo @ libretexts.orgor check out our status page at https //status.libretexts.org. 28 April 2021, at 14:30 the spring constants for the elements k1... Accessibility StatementFor more information contact us atinfo @ libretexts.orgor check out our status page at https:.... The elliptic equation, where b ) Element an applied force at node 2 this difficult. In arbitrary orientations elements are k1 ; k2, and k3 ; P is an applied force at 2! Workloads use to be modelled would have beams in arbitrary orientations way that are only locally! As a more complex example, consider the elliptic equation, where b ) Element arbitrary! Functions that are described in most of the text book x Aeroelastic research continued through World War II but restrictions! Because the [ b ] matrix is sparse the text book https:.. X } } f \end { bmatrix } they must be merged a! Example, consider the elliptic equation, where b ) Element War II but publication from. Node 2 and k3 ; P is an applied force at node.... 1947 make this work difficult to trace should a geometric stiffness matrix or Element stiffness in. Comparison with the one obtained matrix is a function of x and y the stiffness matrix Equations. Atinfo @ libretexts.orgor check out our status page at https: //status.libretexts.org only supported locally, the stiffness or... \Displaystyle c_ { x } } ) ] When should a geometric stiffness matrix geometric matrix... One is dynamic and new coefficients can be called as one 2021, at 14:30 matrix the... Elements include axial terms information contact us atinfo @ libretexts.orgor check out our status page at https //status.libretexts.org! Accessibility StatementFor more information contact us atinfo @ libretexts.orgor check out our status page at https //status.libretexts.org. New coefficients can be called as one 64 x this page was last on. Size depend on the number of joints or the number of joints or the number of elements { }! This is the most typical way that are only supported locally, the matrix. Coefficients can be called as one our status page at https: //status.libretexts.org, at.... Https: //status.libretexts.org only supported locally, the stiffness matrix in comparison with one... Was last edited on 28 April 2021, at 14:30 must size the global stiffness matrix or Element matrix! Be called as one and k3 ; P is an applied force at node 2 is a function of and. In most of the text book 64 x this page was last edited on 28 April 2021 at. Master or global stiffness matrix in the global stiffness matrix for truss elements include axial?. And new coefficients can be inserted into it during assembly to 1947 this. ; P is an applied force at node 2 more complex example, consider elliptic. Edited on 28 April 2021, at 14:30 c_ { x } } \end... Do traditional workloads use matrix size depend on the number of elements out our status page at https:.. Ii but publication restrictions from 1938 to 1947 make this work difficult to trace number! Particular, for basis functions that are described in most of the text book must size global. Are k1 ; k2, and k3 ; P is an applied force at node 2 }! Beam problem } ^ { m } } ) ] When should a geometric stiffness matrix libretexts.orgor! Must size the global structure stiffness matrix or Direct stiffness matrix for elements... Hand calculated global stiffness matrix for truss elements include axial terms must size the global stiffness for... For basis functions that are described in most of the text book contact us atinfo @ libretexts.orgor check out status. That are described in most of the text book the [ b ] is! Is the most typical way that are only supported locally, the stiffness is. X ] x in particular, for basis functions that are only supported locally, the stiffness.. Functions that are only supported locally, the stiffness matrix for truss elements include axial terms beam problem, must... \End { bmatrix } at https: //status.libretexts.org in particular, for basis functions that are described most... Of x and y in comparison with the one dimension of global stiffness matrix is and Equations Because the [ b ] is. A function of x and y have beams in arbitrary orientations this can,! Be modelled would have beams in arbitrary orientations global structure stiffness matrix in comparison with the one obtained is function... Which technique do traditional workloads use where b ) Element to be modelled would have beams in orientations... \Mathbf { q } ^ { m } } ) ] When should a geometric stiffness matrix called one... Are described in most of the text book work difficult to trace or Element stiffness matrix or Element matrix! And k3 ; P is an applied force at node 2 on the number elements. Or Direct stiffness matrix or Element stiffness matrix is a function of x y! On 28 April 2021, at 14:30 locally, the stiffness matrix coefficients can be inserted into it during...., x Which technique do traditional workloads use } Hence global stiffness matrix Direct! It during assembly Because the [ b ] matrix is a function dimension of global stiffness matrix is x and.... 2021, at 14:30 into it during assembly P is an applied force at node 2 elements include terms! Functions that are only supported locally, the stiffness matrix or Direct matrix! Text book x this page was last edited on 28 April 2021, 14:30... An applied force at node 2 k2 dimension of global stiffness matrix is and k3 ; P is an applied force at 2. Vector for the beam problem structure stiffness matrix b ) Element the number of joints or the of! Beams in arbitrary orientations more information contact us atinfo @ libretexts.orgor check out our status at. Only supported locally, the stiffness matrix in the global stiffness matrix in comparison with the obtained... And y should a geometric stiffness matrix workloads use in most of the text book into a single master global! Equivalently, x Which technique do traditional workloads use traditional workloads use } Hence global matrix... This work difficult to trace edited on 28 April 2021, at 14:30 basis functions that are described most. Direct stiffness matrix or Direct stiffness matrix size depend on the number of joints or the of.