• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1990.10a
• /
• pp.28-33
• /
• 1990

This paper presents an efficient and accurate method for nonlinear analysis of frame structures by idealized structural unit method. The main idea behind the present method is to minimize the cost of the computational effort by reducing the number of unknowns. An explicit form of the tangential elastic stiffness matrix of the element is derived by using updated Lagrangian approach. An ultimate limit state of the element is judged on the basis of the formation of a plastic hinge mechanism. The elasto-plastic stiffness matrix and the post-ultimate stiffness matrix of the element are formulated by plastic node method. A comparison between the present method is very efficient and accurate because the computing time required is very small while giving the accurate solution.

• Structural Engineering and Mechanics
• /
• v.4 no.2
• /
• pp.139-148
• /
• 1996

Vector algorithms and the relative importance of the four basic modules (computation of element stiffness matrices, assembly of the global stiffness matrix, solution of the system of linear simultaneous equations, and calculation of stresses and strains) of a finite element computer program for inelastic analysis of reinforced concrete shells are presented. Performance of the vector program is compared with a scalar program. For a cooling tower problem, the speedup factor from the scalar to the vector program is 34 for the element stiffness matrices calculation, 25.3 for the assembly of global stiffness matrix, 27.5 for the equation solver, and 37.8 for stresses, strains and nodal forces computations on a Gray Y-MP. The overall speedup factor is 30.9. When the equation solver alone is vectorized, which is computationally the most intensive part of a finite element program, a speedup factor of only 1.9 is achieved. When the rest of the program is also vectorized, a large additional speedup factor of 15.9 is attained. Therefore, it is very important that all the modules in a nonlinear program are vectorized to gain the full potential of the supercomputers. The vector finite element computer program for inelastic analysis of RC shells with layered elements developed in the present study enabled us to perform mesh convergence studies. The vector program can be used for studying the ultimate behavior of RC shells and used as a design tool.

• Structural Engineering and Mechanics
• /
• v.33 no.3
• /
• pp.307-323
• /
• 2009

This article presents the differential system that governs the mechanical behaviour of a curved-beam element, with varying cross-section area, subjected to generalized load. This system is solved by an exact procedure or by the application of a new numerical recurrence scheme relating the internal forces and displacements at the two end-points of an increase in its centroid-line. This solution has a transfer matrix structure. Both the stiffness matrix and the equivalent load vector are obtained arranging the transfer matrix. New structural matrices have been defined, which permit to determine directly the unknown values of internal forces and displacements at the two supported ends of the curved-beam element. Examples are included for verification.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.11 no.7
• /
• pp.2661-2669
• /
• 2010

In this study, a finite element formulation based first-order shear deformation theory is developed for non-linear behaviors of laminated composite plates containing matrix cracking. The multi-directional stiffness degradation is developed for adopting the stiffness variation induced from matrix cracking, which is proposed by Duan and Yao. The matrix cracking can be expressed in terms of the variation of material properties, such as Young's modulus, shear modulus and Possion ratio of plates, and sequently it is possible to predict the variation of the local stiffness. Using the assumed natural strain method, the present shell element generates neither membrane nor shear locking behavior. Numerical examples demonstrate that the present element behaves quite satisfactorily either for the linear or geometrical nonlinear analysis of laminated composite plates. The results of laminated composite plates with matrix cracking may be the benchmark test for the non-linear analysis of damaged laminated composite plates.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.19 no.10
• /
• pp.1021-1027
• /
• 2009

Finite element models of dynamic systems can be updated in two stages. In the first stage, mass and stiffness matrices are updated neglecting damping. In the second stage, a damping matrix is estimated with the mass and stiffness matrices fixed. Methods to estimate a damping matrix for this purpose are proposed in this paper. For a system with proportional damping, a damping matrix is estimated using the modal parameters extracted from the measured responses and the modal matrix calculated from the mass and stiffness matrices from the first stage. For a system with non-proportional damping, a damping matrix is estimated from the impedance matrix which is the inverse of the FRF matrix. Only one low or one column of the FRF matrix is measured, and the remaining FRFs are synthesized to obtain a full FRF matrix. This procedure to obtain a full FRF matrix saves time and effort to measure FRFs.

• Structural Engineering and Mechanics
• /
• v.62 no.6
• /
• pp.759-769
• /
• 2017

The effect of central axial load on natural frequencies of various thin-walled beams, are investigated by some researchers using different methods such as finite element, transfer matrix and dynamic stiffness matrix methods. However, there are situations that the load will be off centre. This type of loading is called eccentric load. The effect of the eccentricity of axial load on the natural frequencies of asymmetric thin-walled beams is a subject that has not been investigated so far. In this paper, the mentioned effect is studied using exact dynamic stiffness matrix method. Flexure and torsion of the aforesaid thin-walled beam is based on the Bernoulli-Euler and Vlasov theories, respectively. Therefore, the intended thin-walled beam has flexural rigidity, saint-venant torsional rigidity and warping rigidity. In this paper, the Hamilton‟s principle is used for deriving governing partial differential equations of motion and force boundary conditions. Throughout the process, the uniform distribution of mass in the member is accounted for exactly and thus necessitates the solution of a transcendental eigenvalue problem. This is accomplished using the Wittrick-Williams algorithm. Finally, in order to verify the accuracy of the presented theory, the numerical solutions are given and compared with the results that are available in the literature and finite element solutions using ABAQUS software.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.192-199
• /
• 2002

The use of frequency-dependent dynamic stiffness matrix (or spectral element matrix) in structural dynamics may provide very accurate solutions, while it reduces the number of degrees-of-freedom to improve the computational efficiency and cost problems. Thus, this paper develops a spectral element model for the thin plates moving with constant speed under uniform in-plane tension. The concept of Kantorovich method is used in the frequency-domain to formulate the dynamic stiffness matrix. The present spectral element model is evaluated by comparing its solutions with the exact analytical solutions. The effects of moving speed and in-plane tension on the flexural wave dispersion characteristics and natural frequencies of the plate are numerically investigated.

• Structural Engineering and Mechanics
• /
• v.21 no.5
• /
• pp.539-551
• /
• 2005

Finite elements based on isoparametric formulation are known to suffer spurious stiffness properties and corresponding stress oscillations, even when care is taken to ensure that completeness and continuity requirements are enforced. This occurs frequently when the physics of the problem requires multiple strain components to be defined. This kind of error, commonly known as locking, can be circumvented by using reduced integration techniques to evaluate the element stiffness matrices instead of the full integration that is mathematically prescribed. However, the reduced integration technique itself can have a further drawback - rank deficiency, which physically implies that spurious energy modes (e.g., hourglass modes) are introduced because of reduced integration. Such instability in an existing stiffness matrix is generally detected by means of an eigenvalue test. In this paper we show that a knowledge of the dimension of the solution space spanned by the column vectors of the strain-displacement matrix can be used to identify the instabilities arising in an element due to reduced/selective integration techniques a priori, without having to complete the element stiffness matrix formulation and then test for zero eigenvalues.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.46 no.5
• /
• pp.61-68
• /
• 2004

Methodologies of the finite element and boundary element are combined to achieve an efficient and accurate analysis model of frame structure containing shear wall. This model analyzes the frame by employing the finite element method and the shear wall by boundary element method. This study is applicable to a specific situation, where the boundary element is surrounded by finite elements. By employing FE dominant method in which boundary stiffness matrix is transformed into finite element stiffness matrix, boundary element and finite element method are combined to analyze frame structure with walls.

• Journal of the Korean Society for Precision Engineering
• /
• v.27 no.2
• /
• pp.109-116
• /
• 2010

It is important to make a mechanical structure precisely and reasonably in predicting the dynamic characteristics, controlling the vibration, and designing the structure dynamics. In finite element analysis model updating is appropriate as the design parameter is used to analyze the dynamic system. The errors can be contained from the physical parameters and the element modeling. From the dynamic test, more precise dynamic characteristics can be obtained. In this paper, model updating algorithm is developed using frequency difference between experiment and calculation. Modal frequencies are obtained by experiment and finite element analysis for beams with various cross section and shapes which have added masses and holes in the middle. For plates with and without groove, experiment and analyses are carried out by applying free boundary conditions as well. Mass and stiffness matrices are updated by comparing test and analytical modal frequencies. The result shows that the updated frequencies become closer to the test frequencies in case that both matrices are updated. An improved analytical model is obtained by changing model parameters such that the discrepancy between test and finite element frequencies is minimized. For beam and plate models updating of mass and stiffness matrices can improve the dynamical behavior of the model by acting on the physical parameters such as masses and stiffness.