An efficient method of model order reduction is proposed for the dynamic computation of a flexible multibody system undergoing both large overall motions and large deformations. The system is initially modeled by using the nonlinear finite elements of absolute nodal coordinate formulation and then locally linearized at a series of quasi-static equilibrium configurations according to the given accuracy in dynamic computation. By using the Craig-Bampton method, the reduced model is established by projecting the incremental displacements of the locally linearized system onto a set of local modal bases at the quasi-static equilibrium configuration accordingly. Afterwards, the initial conditions for the dynamic computation for the reduced model via the generalized-α integrator can be determined from the modal bases. The analysis of computation complexity is also performed. Hence, the proposed method gives time-varying and dimension-varying modal bases to elaborate the efficient model reduction. Finally, three examples are presented to validate the accuracy and efficiency of the proposed method.
KEYWORDSabsolute nodal coordinate formulation, Craig-Bampton method, flexible multibody dynamics, model order reduction Nomenclature: C, vector of constraints; d, vector of augmented external forces; e, vector of generalized coordinates; E, Young's modulus; I, moment of inertia; h, time step; K, global stiffness matrix;K, augmented global stiffness matrix; M, global mass matrix;M, augmented global mass matrix; n, total number of degrees of freedom; n c , number of constraint equations; n e , number of finite elements; n m , number of dynamic equations; n p , number of degrees of freedom per node; n s , revolution per minute; n t , number of truncated modes; n t 1 , initial number of truncated modes; n t max , maximal number of truncated modes; q, vector of augmented modal coordinates; Q, vector of external forces; r, vector of global position; t, time; t N , moment of the Nth local linearization; T, transformation matrix;ũ, vector of augmented coordinates; U, strain energy; XYZ, global frame of coordinates; X i Y i Z i , the ith body frame of coordinates; C , matrix of constrained modes; V , matrix of vibration modes; , vector of modal coordinates of internal part; , vector of Lagrange multipliers; j , the jth eigenvalue; , diagonal matrix of eigenvalues; Ω, configuration region; , vector of augmented acceleration-like variables; s , rotation angle; s , angular velocity; Δ max , threshold of local linearization; , specified value of mode selection; Δ RMS , root of mean squared error. Superscripts and subscripts: () B , symbol related to boundary part; () I , symbol related to internal part; () BI , () IB , symbol related to coupling parts; () e , the Jacobian with respect to vector e; () N , symbol corresponding to the Nth local linearization;() k , symbol corresponding to the kth time step; () d , reserved quantity; () r , residual quantity. Overbar: (^), quantity in the reduced model; (˜), quantity in the non-reduced model; (...