• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 1997.10a
• /
• pp.151-158
• /
• 1997

The main concern of numerical dynamic analysis of large structures is to find an acceptable solution with fewer mode shapes and less computational efforts. component mode method utilizes substructure technique to reduce the degrss of freedom but have a disadvantage to not consider the dynamic characteristics of loads. Ritz Vector method consider the load characteristics but requires many integrations and errors are accumulated. In this study, to prove the effectiveness of component mode method, Lanczos algorithm are introduced. To prove the effectiveness of this method, example structures areanalyzed and the results are compared with SAP90.

• Structural Engineering and Mechanics
• /
• v.81 no.4
• /
• pp.411-428
• /
• 2022

This paper presents the Campbell diagram analysis of the rotordynamic system using the full order model (FOM) and the reduced order model (ROM) techniques to determine the critical speeds, identify the stability and reduce the computational time. Due to the spin-speed-dependent matrices (e.g., centrifugal stiffening matrix), several model order reduction (MOR) techniques may be considered, such as the modal superposition (MS) method and the Krylov subspace-based MOR techniques (e.g., Ritz vector (RV), quasi-static Ritz vector (QSRV), multifrequency quasi-static Ritz vector (MQSRV), multifrequency/ multi-spin-speed quasi-static Ritz vector (MMQSRV) and the combined Ritz vector & modal superposition (RV+MS) methods). The proposed MMQSRV method in this study is extended from the MQSRV method by incorporating the rotational-speed-dependent stiffness matrices into the Krylov subspace during the MOR process. Thus, the objective of this note is to respond to the question of whether to use the MS method or the Krylov subspace-based MOR technique in establishing the Campbell diagram of the shaft-disc-blade assembly systems in three-dimensional (3D) finite element analysis (FEA). The Campbell diagrams produced by the FOM and various MOR methods are presented and discussed thoroughly by computing the norm of relative errors (ER). It is found that the RV and the MS methods are dominant at low and high rotating speeds, respectively. More precisely, as the spinning velocity becomes large, the calculated ER produced by the RV method is significantly increased; in contrast, the ER produced by the MS method is smaller and more consistent. From a computational point of view, the MORs have substantially reduced the time computing considerably compared to the FOM. Additionally, the verification of the 3D FE rotordynamic model is also provided and found to be in close agreement with the existing solutions.

• Computational Structural Engineering
• /
• v.9 no.2
• /
• pp.115-120
• /
• 1996

The main concern of numerical dynamic analysis of large structures is to find an acceptable solution with fewer mode shapes and less computational efforts. Component mode method utilizes substructure technique to reduce the degree of freedom but have a disadvantage to not consider the dynamic characteristics of loads. Ritz Vector method consider the load characteristics but requires many integrations and errors are accumulated. In this study, to improve the effectiveness of component mode method, Lanczos algorithm is introduced. To prove the effectiveness of this method, example structure are analyzed and the results are compared with SAP90.

• Steel and Composite Structures
• /
• v.29 no.6
• /
• pp.785-802
• /
• 2018

In this paper, two different computational methods, called Rayleigh-Ritz and collocation are developed to estimate the ultimate strength of composite plates. Progressive damage behavior of moderately thick composite laminated plates is studied under in-plane compressive load and uniform lateral pressure. The formulations of both methods are based on the concept of the principle of minimum potential energy. First order shear deformation theory and the assumption of large deflections are used to develop the equilibrium equations of laminated plates. Therefore, Newton-Raphson technique will be used to solve the obtained system of nonlinear algebraic equations. In Rayleigh-Ritz method, two degradation models called complete and region degradation models are used to estimate the degradation zone around the failure location. In the second method, a new energy based collocation technique is introduced in which the domain of the plate is discretized into the Legendre-Gauss-Lobatto points. In this new method, in addition to the two previous models, the new model named node degradation model will also be used in which the material properties of the area just around the failed node are reduced. To predict the failure location, Hashin failure criteria have been used and the corresponding material properties of the failed zone are reduced instantaneously. Approximation of the displacement fields is performed by suitable harmonic functions in the Rayleigh-Ritz method and by Legendre basis functions (LBFs) in the second method. Finally, the results will be calculated and discussions will be conducted on the methods.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.6
• /
• pp.1363-1370
• /
• 1995

A nonlinear dynamic modeling method for simply supported structures undergoing large overall motion is suggested. The modeling method employs Rayleigh-Ritz mode technique and Von Karman nonlinear strain measures. Numerical study shows that the suggested modeling method provides qualitatively different results from those of the Classical Linear Cartesian modeling method. Especially, natural frequency variations and residual deformation due to membrane strain effects are observed in the numerical results obtained by the suggested modeling method.

• Structural Engineering and Mechanics
• /
• v.68 no.5
• /
• pp.603-619
• /
• 2018

In this research, an effective computational technique is carried out for nonlinear and post-buckling analyses of cracked imperfect composite plates. The laminated plates are assumed to be moderately thick so that the analysis can be carried out based on the first-order shear deformation theory. Geometric non-linearity is introduced in the way of von-Karman assumptions for the strain-displacement equations. The Ritz technique is applied using Legendre polynomials for the primary variable approximations. The crack is modeled by partitioning the entire domain of the plates into several sub-plates and therefore the plate decomposition technique is implemented in this research. The penalty technique is used for imposing the interface continuity between the sub-plates. Different out-of-plane essential boundary conditions such as clamp, simply support or free conditions will be assumed in this research by defining the relevant displacement functions. For in-plane boundary conditions, lateral expansions of the unloaded edges are completely free while the loaded edges are assumed to move straight but restricted to move laterally. With the formulation presented here, the plates can be subjected to biaxial compressive loads, therefore a sensitivity analysis is performed with respect to the applied load direction, along the parallel or perpendicular to the crack axis. The integrals of potential energy are numerically computed using Gauss-Lobatto quadrature formulas to get adequate accuracy. Then, the obtained non-linear system of equations is solved by the Newton-Raphson method. Finally, the results are presented to show the influence of crack length, various locations of crack, load direction, boundary conditions and different values of initial imperfection on nonlinear and post-buckling behavior of laminates.

• Advances in nano research
• /
• v.13 no.4
• /
• pp.393-406
• /
• 2022

In the present article, functionally graded small-scaled plates based on modified strain gradient theory (MSGT) are studied for analyzing the nonlinear bending and post-buckling responses. Von-Karman's assumptions are applied to incorporate geometric nonlinearity and the first-order shear deformation theory is used to model the plates. Modified strain gradient theory includes three length scale parameters and is reduced to the modified couple stress theory (MCST) and the classical theory (CT) if two or all three length scale parameters become zero, respectively. The Ritz method with Legendre polynomials are used to approximate the unknown displacement fields. The solution is found by the minimization of the total potential energy and the well-known Newton-Raphson technique is used to solve the nonlinear system of equations. In addition, numerical results for the functionally graded small-scaled plates are obtained and the effects of different boundary conditions, material gradient index, thickness to length scale parameter and length to thickness ratio of the plates on nonlinear bending and post-buckling responses are investigated and discussed.

• Computational Structural Engineering
• /
• v.3 no.3
• /
• pp.125-131
• /
• 1990

A finite element technique for the time dependent large structural problems is presented. It is based on the error estimation for the bases of solution spaces. An a-posteriori energy norm of residual error serves as the error indicator. Mode shapes which are calculated by scaling the Ritz vectors are applied to discretize the continuous spatial domain. Finally, the performance of the proposed methods is demonstrated by solving simple examples.

• Structural Engineering and Mechanics
• /
• v.53 no.2
• /
• pp.337-354
• /
• 2015

In this article, free vibration of functionally graded (FG) elliptic plates subjected to various classical boundary conditions has been investigated. Literature review reveals no study has been performed based on functionally graded elliptic plates till date. The mechanical kinematic relations are considered based on classical plate theory. Rayleigh-Ritz technique is used to obtain the generalized eigenvalue problem. The material properties of the FG plate are assumed to vary along thickness direction of the constituents according to power-law form. Trial functions denoting the displacement components are expressed in simple algebraic polynomial forms which can handle any edge support. The objective is to study the effect of geometric configurations and gradation of constituent volume fractions on the natural frequencies. New results for frequency parameters are incorporated after performing a test of convergence. A comparison study is carried out with existing literature for validation in special cases. Three-dimensional mode shapes for circular and elliptic FG plates are also presented with various boundary conditions at the edges.

• Interaction and multiscale mechanics
• /
• v.1 no.2
• /
• pp.211-230
• /
• 2008

Owing to the growing size of the eigenvalue problem and the growing number of eigenvalues desired, solution methods of iterative nature are becoming more popular than ever, which however suffer from low efficiency and lack of proper convergence criteria. In this paper, three efficient iterative eigenvalue algorithms are considered, i.e., subspace iteration method, iterative Ritz vector method and iterative Lanczos method based on the cell sparse fast solver and loop-unrolling. They are examined under the mode error criterion, i.e., the ratio of the out-of-balance nodal forces and the maximum elastic nodal point forces. Averagely speaking, the iterative Ritz vector method is the most efficient one among the three. Based on the mode error convergence criteria, the eigenvalue solvers are shown to be more stable than those based on eigenvalues only. Compared with ANSYS's subspace iteration and block Lanczos approaches, the subspace iteration presented here appears to be more efficient, while the Lanczos approach has roughly equal efficiency. The methods proposed are robust and efficient. Large size tests show that the improvement in terms of CPU time and storage is tremendous. Also reported is an aggressive shifting technique for the subspace iteration method, based on the mode error convergence criteria. A backward technique is introduced when the shift is not located in the right region. The efficiency of such a technique was demonstrated in the numerical tests.