• Ocean Systems Engineering
• /
• v.4 no.4
• /
• pp.327-342
• /
• 2014

The tension leg platform (TLP) is one of the compliant structures which are generally used for deep water oil exploration. With respect to the horizontal degrees of freedom, it behaves like a floating structure moored by vertical tethers which are pretension due to the excess buoyancy of the platform, whereas with respect to the vertical degrees of freedom, it is stiff and resembles a fixed structure and is not allowed to float freely. In the current study, a numerical study for square TLP using modified Morison equation was carried out in the time domain with water particle kinematics using Airy's linear wave theory to investigate the effect of changing the tether tension force on the stiffness matrix of TLP's, the dynamic behavior of TLP's; and on the fatigue stresses in the cables. The effect was investigated for different parameters of the hydrodynamic forces such as wave periods, and wave heights. The numerical study takes into consideration the effect of coupling between various degrees of freedom. The stiffness of the TLP was derived from a combination of hydrostatic restoring forces and restoring forces due to cables. Nonlinear equation was solved using Newmark's beta integration method. Only uni-directional waves in the surge direction was considered in the analysis. It was found that for short wave periods (i.e., 10 sec.), the surge response consisted of small amplitude oscillations about a displaced position that is significantly dependent on tether tension force, wave height; whereas for longer wave periods, the surge response showed high amplitude oscillations that is significantly dependent on wave height, and that special attention should be given to tethers fatigue because of their high tensile static and dynamic stress.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.23 no.3
• /
• pp.240-248
• /
• 2013

Torsional vibrations, such as the gear rattle of the manual transmission in vehicle systems, are correlated with the firing stroke from the engine. These vibro-impacts can be examined based upon linear time-invariant analysis. In order to understand the gear dynamics, a specific manual transmission with a front-engine front-wheel drive configuration is investigated. A method to reduce the degrees of freedom is suggested based upon the eigensolutions and frequency response functions, which will lead to the development of an efficient matrix size. The dynamic characteristics of single- and dual-mass flywheels are then compared. The effect of the dual-mass flywheel is investigated based upon the mobility analysis, which will lead to understanding of the concepts for avoiding vibro-impacts. A linear time-invariant system model is examined by employing the effective clutch stiffness from a two-stage clutch damper. Thus, the relationship between the dynamic characteristics and the clutch damper can be predicted by assuming a combination of different stage stiffness levels.

• Journal of the Korea Concrete Institute
• /
• v.25 no.5
• /
• pp.529-537
• /
• 2013

The aim of the current study is to develop a nonlinear isoparametric layered frame finite element (FE) analysis of FRP strengthened reinforced concrete (RC) beam or column members by a force-based FE formulation. In sections, concrete is modeled in the triaxial stress-strain relationship state and the FRP sheet is modeled as layered composite materials in two-dimension. The element stiffness matrix derived by the force-based FE has the force-interpolation functions without assuming the displacement shape functions. A lateral load test of RC column strengthened by GFRP sheets was analyzed by the developed force-based FE model. From comparative studies of the experimental and analysis results, it was shown to compare with the stiffness FE method that the force-based FE analysis could give more accurate predictions in the overall lateral load-deflection response as well as in nonlinear deformations and damages in the column plastic hinge region.

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

• Bulletin of the Society of Naval Architects of Korea
• /
• v.22 no.1
• /
• pp.21-30
• /
• 1985

This paper deals whth the buckling as well as postbuckling analysis of axisymmertric shells taking the initial deflection effects into account. Incremental equilibrium equations, based on the principle of virtual work, were derived by the finite element method, the successive step-by-step Newton-Raphson iterative technique was adopted. To define the transition pattern of postbuckling behavior from the prebuckling state more accurately, a simple solution method was developed, i.e. the critical load was calculated by the load extrapolation method with the determinant of tangent stiffness matrix and the equilibrium configuration in the immediate postbuckling stage was obtained by perturbation scheme and eigenvalue analysis. Degenerated isoparametric shell elements were used to analyse the axisymmetric shell of revolution. And by the method developed in this paper, the computer program applicable to the nonlinear analysis of both thin and moderately thick shells was constructed. To verify the capabilities and accuracies of the present solution method, the computed results were compared with the results of analytical solutions. These results coincided fairly well in both the small deflection and large deflection ranges. Various numerical analyses were done to show the effect of initial deflection and shape of shells on buckling load and postbuckling behavior. Futhermore, corrected directions of applied loads at every increment steps were used to determine the actual effects of large deflection in non-conservative load systems such as hydrostatic pressure load. The following conclusions can be obtained. (1) The method described in this paper was found to be both economic and effective in calculating buckling load and postbuckling behavior of shell structure. (2) Buckling and postbuckling behavior of spherical caps is critically dependent upon their geometric configuration, i.e. the shape of spherical cap and quantities of the initial deflection. (3) In the analysis of large deflection problems of shells by the incremental method, corrections of the applied load directions are needed at every incremental step to compensate the follower force effects.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.4
• /
• pp.202-213
• /
• 1992

For the analysis and design of ship structures the generalized slope deflection method(GSDM) taking account of axial elongation effect as well as the bending and shearing deformation is developed. Using the span point concept, the existing slope deflection method is easy to transform the variable section to the equivalent uniform one under the bending moment and the shear force, but it is difficult to analyze the web frame with inclined members because the axial deformation effect is not considered. In the present method, the equilibrium conditions including all force components(i.e. axial force, shear force, bending moment) are formulated at the both ends of the variable section beam, such that the usual space frame stiffness equation which can be solved easily by the matrix method is derived. The accuracy and applicability of the present method is demonstrated by analyzing the ship web frame structures.

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

• Structural Engineering and Mechanics
• /
• v.12 no.6
• /
• pp.657-668
• /
• 2001

Fiber reinforced cementitious composites are nowadays widely applied in civil engineering. The postcracking performance of this material depends on the interaction between a steel fiber, which is obliquely across a crack, and its surrounding matrix. While the partly debonded steel fiber is subjected to pulling out from the matrix and simultaneously subjected to transverse force, it may be modelled as a Bernoulli-Euler beam partly supported on an elastic foundation with non-linearly varying modulus. The fiber bridging the crack may be cut into two parts to simplify the problem (Leung and Li 1992). To obtain the transverse displacement at the cut end of the fiber (Fig. 1), it is convenient to directly solve the corresponding differential equation. At the first glance, it is a classical beam on foundation problem. However, the differential equation is not analytically solvable due to the non-linear distribution of the foundation stiffness. Moreover, since the second order deformation effect is included, the boundary conditions become complex and hence conventional numerical tools such as the spline or difference methods may not be sufficient. In this study, moment equilibrium is the basis for formulation of the fundamental differential equation for the beam (Timoshenko 1956). For the cantilever part of the beam, direct integration is performed. For the non-linearly supported part, a transformation is carried out to reduce the higher order differential equation into one order simultaneous equations. The Runge-Kutta technique is employed for the solution within the boundary domain. Finally, multi-dimensional optimization approaches are carefully tested and applied to find the boundary values that are of interest. The numerical solution procedure is demonstrated to be stable and convergent.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.6
• /
• pp.1299-1308
• /
• 1994

This paper is considered on the dynamic behavior and the dynamic impact coefficient on the cable-stayed bridge under the vehicle load. The method of static analysis, that is, the transfer matrix method is used to get influence values about displacements, section forces of girder and cable forces. Gotten influence values were used as basic data to analyse dynamic behavior. This paper used the transfer matrix method because it is relatively simpler than the finite element method, and calculating speed of computer is very fast and the precision of computation is high. In the process of dynamic analysis, the uncoupled equation of motion is derived from simultaneous equation of the motion of cable-stayed bridge and vehicle travelling by using mode shape, which was borne from system of undamped free vibration. The solution of the uncoupled equation of motion, that is, time history of response of deflections, velocity and acceleration on reference coordinate system, is found by Newmark-${\beta}$ method, a kind of direct integral method. After the time history of dynamic response was gotten, and it was transfered to the time history of dynamic response of cable-stayed bridge by linear transformation of coordinates. As a result of this numerical analysis, effect of dynamic behavior for cable-stayed bridge under the vehicle load has varied depending on parameter of design, that is, the ratio of span, the ratio of main span length, tower height, the flexural rigidity of longitudinal girder, the flexural rigidity of tower, and the cable stiffness, investigated. Very good agreements with the existing solution in the literature are shown for the uncracked plate as well as the cracked plate.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.5
• /
• pp.475-482
• /
• 2010

In this paper, an efficient two-level parallel domain decomposition algorithm is suggested to solve large-DOF structural problems. Each subdomain is composed of the coarse problem and local problem. In the coarse problem, displacements at coarse nodes are computed by the iterative method that does not need to assemble a stiffness matrix for the whole coarse problem. Then displacements at local nodes are computed by Multi-Frontal Sparse Solver. A parallel version of PCG(Preconditioned Conjugate Gradient Method) is developed to solve the coarse problem iteratively, which minimizes the data communication amount between processors to increase the possible problem DOF size while maintaining the computational efficiency. The test results show that the suggested algorithm provides scalability on computing performance and an efficient approach to solve large-DOF structural problems.