A self-regularized approach for deriving the free-free flexibility and stiffness matrices

Using the SVD, we make the singular matrix become the nonsingular bordered matrix.The constants and constraints of the bordered matrix can be explained.We derive the free-free flexibility/stiffness matrices from the bordered matrices.The equilibrium/compatibility of the unnatural force/displacement can be tested. Motivated by Fichera's idea for regularizing the rank-deficiency model, we derive the free-free flexibility matrices by inverting the bordered stiffness matrix. The singular stiffness matrix of a free-free structure is expanded to a bordered matrix by adding n slack variables, where n is the nullity of the singular stiffness matrix. Besides, the corresponding n constraints are accompanied to result in a nonsingular matrix. The constraints filter out the homogeneous solution for the regularized solution. By inverting the nonsingular matrix, we can obtain the free-free flexibility matrix from the submatrices. The value of the extra degree of freedom shows the role of no solution (nonzero case) or infinite solution (zero case) with respect to the loading vector. After constructing the bordered system, the equilibrium of the specified force and the compatibility of the specified displacement can be tested according the zero slack variable. Similarly, the free-free flexibility matrix is obtained from the free-free stiffness matrix. Finally, four examples, a rod with symmetric stiffness, a plane truss, a beam and a bar with unsymmetric stiffness, were demonstrated to see the validity of the present formulation.

[1]  Alphose Zingoni,et al.  Group-theoretic insights on the vibration of symmetric structures in engineering , 2014, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[2]  T. Healey,et al.  Exact block diagonalization of large eigenvalue problems for structures with symmetry , 1991 .

[3]  Carlos A. Felippa,et al.  The construction of free–free flexibility matrices for multilevel structural analysis , 2002 .

[4]  S. Mukherjee,et al.  Elimination of rigid body modes from discretized boundary integral equations , 1998 .

[5]  Decomposition of symmetric mass–spring vibrating systems using groups, graphs and linear algebra , 2006 .

[6]  A. Zingoni A group-theoretic finite-difference formulation for plate eigenvalue problems , 2012 .

[7]  Shyh-Rong Kuo,et al.  Regularization methods for ill-conditioned system of the integral equation of the first kind with the logarithmic kernel , 2014 .

[8]  E. Riks An incremental approach to the solution of snapping and buckling problems , 1979 .

[9]  Carlos A. Felippa,et al.  The construction of free–free flexibility matrices as generalized stiffness inverses , 1998 .

[10]  A. Eriksson USING EIGENVECTOR PROJECTIONS TO IMPROVE CONVERGENCE IN NON-LINEAR FINITE ELEMENT EQUILIBRIUM ITERATIONS. , 1987 .

[11]  R. Vodička,et al.  On the removal of the non‐uniqueness in the solution of elastostatic problems by symmetric Galerkin BEM , 2006 .

[12]  Jeng-Tzong Chen,et al.  Derivation of stiffness and flexibility for rods and beams by using dual integral equations , 2008 .

[13]  P. Wriggers Nonlinear finite element analysis of solids and structures , 1998 .

[14]  R. Hartwig Singular Value Decomposition and the Moore–Penrose Inverse of Bordered Matrices , 1976 .

[15]  Federico París,et al.  Note on the removal of rigid body motions in the solution of elastostatic traction boundary value problems by SGBEM , 2006 .

[16]  Carlos A. Felippa,et al.  A direct flexibility method , 1997 .

[17]  Adi Ben-Israel,et al.  Generalized inverses: theory and applications , 1974 .

[18]  V. Mantič,et al.  ON THE REMOVAL OF RIGID BODY MOTIONS IN THE SOLUTION OF ELASTOSTATIC PROBLEMS BY DIRECT BEM , 1996 .