• Proceedings of the KSR Conference
• /
• 2001.10a
• /
• pp.500-507
• /
• 2001

In this study, Interaction among each element was construed by the unit of new stiffness matrix to analyze the finite element about Japanese precast concrete slab track and improved slab track. Dynamic analysis which is assumed a static analysis and a trainload on the transverse and the longitudinal load of the train into a series periodic function was performed by using the common program. And then, the difference of the movement between an improved section and an existing structure type was realized. Longitudinal static analysis indicated that the stress of the improved section is smaller than that of the protrusion of the existing slab track. And static and dynamic analysis on transverse load showed a little decrease of the displacement on new slab track. But the dynamic analysis result showed that new track system was considerably decreased by 30% compared with the existing Japanese slab track.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.4
• /
• pp.439-450
• /
• 2003

Flat plate system has been adopted in many buildings constructed recently because of the advantage of reduced floor heights to meet the economical and architectural demands. Structural engineers commonly use the effective beam width model(EBWM) in practical engineering for the analysis of flat plate structures. However, in many cases, when it is difficult to use the EBWM, it is necessary to use a refined finite element model for an accurate analysis. But it would take significant amount of computational time and memory if the entire building structure was subdivided with finer meshes. An efficient analytical method is proposed in this study to obtain accurate results in significantly reduced computational time. The proposed method employs super elements developed using the matrix condensation technique and fictitious beams are used in the development of super elements to enforce the compatibility at the interfaces of super elements. The stiffness degradation of flat plate system considered in the EBWM was taken into account by reducing the elastic modulus of floor slabs in this study. Static and dynamic analyses of example structures were performed and the efficiency and accuracy of the proposed method were verified by comparing the results with those of the refined finite element model and the EBWM.

• Structural Engineering and Mechanics
• /
• v.38 no.1
• /
• pp.27-37
• /
• 2011

This paper presents the results of a study on the capability of nonlinear quasi-static finite element modelling in simulating the hysteretic behaviour of CFRP and GFRP-retrofitted RC exterior beam-column joints under cyclic loads. Four specimens including two plain and two CFRP/GFRP-strengthened beam-column joints tested by Mahini and Ronagh (2004) and other researchers are modelled using ANSYS. Concrete in compression is defined by the modified Hognestad model and anisotropic multi-linear model is employed for modelling the stress-strain relations in reinforcing bars while anisotropic plasticity is considered for the FRP composite. Both concrete and FRP are modelled using solid elements whereas space link elements are used for steel bars considering a perfect bond between materials. A step by step load increment procedure to simulate the cyclic loading regime employed in the testing. An automatically reforming stiffness matrix strategy is used in order to simulate the actual seismic performance of the RC concrete after cracking, steel yielding and concrete crushing during the push and pull loading cycles. The results show that the hysteretic simulation for all specimens is satisfactory and therefore suggest that the numerical model can be used as an inexpensive tool to design of FRP-strengthened RC beam-column joints under cyclic loads.

• Journal of computational fluids engineering
• /
• v.10 no.2
• /
• pp.38-47
• /
• 2005

Substructuring methods are often used in finite element structural analyses. In this study a multi-level substructuring(MLSS) algorithm is developed and proposed as a possible candidate for finite element fluid solvers. The present algorithm consists of four stages such as a gathering, a condensing, a solving and a scattering stage. At each level, a predetermined number of elements are gathered and condensed to form an element of higher level. At the highest level, each sub-domain consists of only one super-element. Thus, the inversion process of a stiffness matrix associated with internal degrees of freedom of each sub-domain has been replaced by a sequential static condensation of gathered element matrices. The global algebraic system arising from the assembly of each sub-domain matrices is solved using a well-known iterative solver such as the conjugare gradient(CG) or the conjugate gradient squared(CGS) method. A time comparison with CG has been performed on a 2-D Poisson problem. With one domain the computing time by MLSS is comparable with that by CG up to about 260,000 d.o.f. For 263,169 d.o.f using 8 x 8 sub-domains, the time by MLSS is reduced to a value less than $30\%$ of that by CG. The lid-driven cavity problem has been solved for Re = 3200 using the element interpolation degree(Deg.) up to cubic. in this case, preconditioning techniques usually accompanied by iterative solvers are not needed. Finite element formulation for the incompressible flow has been stabilized by a modified residual procedure proposed by Ilinca et al.[9].

• 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.

• Journal of the Korean Society for Railway
• /
• v.2 no.3
• /
• pp.36-45
• /
• 1999

The basic idea of cable-stayed girder bridges is the utilization of high strength cables to provide intermediate supports for the bridge girder so that the girder can span a much longer distance. In the cable-stayed bridge, the cables exhibit nonlinear behavior because of the change in sag, due to the dead weight of the cable, which occurs with changing tension in the cable resulting from the movement of the end points of the cable as the bridge is loaded. Techniques required for the static analysis of cable-stayed bridges has been developed by many researchers. However, little work has been done on the dynamic analysis of such structures. To investigate the characteristics of the dynamic response of long-span cable-stayed bridges due to various dynamic loadings likes moving traffic loads. two different 3-D cable-stayed bridge models are considered in this study. Two models are exactly the same in structural configurations but different in finite element discretization. Modal analysis is conducted using the deformed dead-load tangent stiffness matrix. A new concept was presented by using divided a cable into several elements in order to study the effect of the cable vibration (both in-plane and swinging) on the overall bridge dynamics. The result of this study demonstrates the importance of cable vibration on the overall bridge dynamics.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.1A
• /
• pp.31-43
• /
• 2009

In this study, the numerical method is presented, which can consider the various train types and can solve the equations of motion for a vehicle-bridge interaction analysis by non-iteration procedure through formulating the coupled equations of motion. The coupled equations of motion for the vehicle-bridge interaction are solved by the Newmark ${\beta}$ of a direct integration method, and by composing the effective stiffness matrix and the effective force vector according to a analysis step, those can be solved with the same manner of the solving procedure of equilibrium equations in static analysis. Also, the effective stiffness matrix is reconstructed by the Skyline method for increasing the analysis effectiveness. The Cholesky's matrix decomposition scheme is applied to the analysis procedure for minimizing the numerical errors that can be generated in directly calculating the inverse matrix. The equations of motion for the conventional trains are derived, and the numerical models of the conventional trains are idealized by a set of linear springs and dashpots with 16 degrees of freedom. The bridge models are simplified by the 3 dimensional space frame element which is based on the Euler-Bernoulli theory. The rail irregularities of vertical and lateral directions are generated by the PSD functions of the Federal Railroad Administration (FRA). The results of the vehicle-bridge interaction analysis are verified by the experimental results for the railway plate girder bridges of a span length with 12 m, 18 m, and the experimental and analytical data are applied to the low pass filtering scheme, and the basis frequency of the filtering is a 2 times of the 1st fundamental frequency of a bridge bending.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.6 no.2
• /
• pp.121-131
• /
• 1986

The stiffness matrix of circularly curved frame elements including the warping effects is formulated by the solutions of vlasov's differential equations, and the procedure for the elastic static analysis of curved girder systems by the displacement method is presented. The validity of this method has been demonstrated by comparing the analysis results with other solutions. And if the tangential lines of the two frame element axes connected at any nodal point coincide, the transformation to the global coordinate system can be omitted when we analyze the structures consisting of circularly curved elements. The theory introduced in this thesis can be applied with sufficient accuracy to the structures built up with horizontally circular curved frame elements which have closed or open cross sections and are symmetric to the axis perpendicular to the plane of the curvature, such as prestressed concrete box girder bridges.

• Structural Engineering and Mechanics
• /
• v.13 no.5
• /
• pp.557-568
• /
• 2002

For the evaluation of the capability of a tubular member of an offshore structure to absorb the collision energy, a simple method can be employed for the collision analysis without performing the detailed analysis. The most common simple method is the rigid-plastic method. However, in this method any characteristics for horizontal movement and rotation at the ends of the corresponding tubular member are not included. In a real structural system of an offshore structure, tubular members sustain a certain degree of elastic support from the adjacent structure. End fixity has influences in the behaviors of a tubular member. Three-dimensional FEM analysis can include the effect of end fixity fully, however in viewpoints of the inherent computational complexities of the 3-D approach, this is not the recommendable analysis at the initial design stage. In this paper, influence of end fixity on the behaviors of a tubular member is investigated, through a new approach and other approaches. A new analysis approach that includes the flexibility of the boundary points of the member is developed here. The flexibility at the ends of a tubular element is extracted using the rational reduction of the modeling characteristics. The property reduction is based on the static condensation of the related global stiffness matrix of a model to end nodal points of the tubular element. The load-displacement relation at the collision point of the tubular member with and without the end flexibility is obtained and compared. The new method lies between the rigid-plastic method and the 3-demensional analysis. It is self-evident that the rigid-plastic method gives high strengthening membrane effect of the member during global deformation, resulting in a steeper slope than the present method. On the while, full 3-D analysis gives less strengthening membrane effect on the member, resulting in a slow going load-displacement curve. Comparison of the load-displacement curves by the new approach with those by conventional methods gives the figures of the influence of end fixity on post-yielding behaviors of the relevant tubular member. One of the main contributions of this investigation is the development of an analytical rational procedure to figure out the post-yielding behaviors of a tubular member in offshore structures.

• Journal of Korean Society of Steel Construction
• /
• v.10 no.1 s.34
• /
• pp.47-61
• /
• 1998

The behavior of horizontally curved I-girder bridges is complex because the flexural and torsional behavior of curved girders are coupled due to their initial curvature. Also, the behavior is affected by cross beams. To investigate the behavior of horizontally curved I-girder bridges, it is necessary to consider curved girders with cross beams. In order to perform free vibration analyses of horizontally curved I-girder bridges, a finite element formulation is presented here and a finite element analysis program is developed. The formulation that is presented here consists of curved and straight beam elements, including the warping degree of freedom. Based on the theory of thin-walled curved beams, the shape functions of the curved beam elements are derived from homogeneous solutions of the static equilibrium equations. Third-order hermits polynomials are used to form the shape functions of the straight beam elements. In the finite element analysis program, global stiffness and mass matrix are composed, based on the Cartesian coordinate system. The Gupta method is used to efficiently solve the eigenvalue problem. Comparing the results of several examples here with those of previous studies, the formulation presented is verified. The validity of the program developed is shown by comparing results with those analyzed by the shell element.