• Smart Structures and Systems
• /
• v.13 no.4
• /
• pp.587-614
• /
• 2014

In recent years the interest in online monitoring of lightweight structures with ultrasonic guided waves is steadily growing. Especially the aircraft industry is a driving force in the development of structural health monitoring (SHM) systems. In order to optimally design SHM systems powerful and efficient numerical simulation tools to predict the behaviour of ultrasonic elastic waves in thin-walled structures are required. It has been shown that in real industrial applications, such as airplane wings or fuselages, conventional linear and quadratic pure displacement finite elements commonly used to model ultrasonic elastic waves quickly reach their limits. The required mesh density, to obtain good quality solutions, results in enormous computational costs when solving the wave propagation problem in the time domain. To resolve this problem different possibilities are available. Analytical methods and higher order finite element method approaches (HO-FEM), like p-FEM, spectral elements, spectral analysis and isogeometric analysis, are among them. Although analytical approaches offer fast and accurate results, they are limited to rather simple geometries. On the other hand, the application of higher order finite element schemes is a computationally demanding task. The drawbacks of both methods can be circumvented if regions of complex geometry are modelled using a HO-FEM approach while the response of the remaining structure is computed utilizing an analytical approach. The objective of the paper is to present an efficient method to couple different HO-FEM schemes with an analytical description of an undisturbed region. Using this hybrid formulation the numerical effort can be drastically reduced. The functionality of the proposed scheme is demonstrated by studying the propagation of ultrasonic guided waves in plates, excited by a piezoelectric patch actuator. The actuator is modelled utilizing higher order coupled field finite elements, whereas the homogenous, isotropic plate is described analytically. The results of this "semi-analytical" approach highlight the opportunities to reduce the numerical effort if closed-form solutions are partially available.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.12
• /
• pp.2121-2133
• /
• 1997

A finite element code is developed for 3-D self-contact problems, using continuum elements with a SRI(Selective Reduced Integration) scheme to prevent locking phenomenon by the incompressibility of rubber. Contact treatment is carried out in two ways : using the displacement constraints in case of rigid contact ; and imposing the same contact forces on two contact boundaries in case of self-contact. The finite element code developed is applied to the deformation analysis of C.V.joint boots which maintain lubrication conditions and protect the C.V.joint assembly from impact and dust. The boot accompanies large rotation depending on the rotation of the wheel axis and leading to the self-contact phenomena of the boot bellows. Since this contact phenomenon causes wear of the product and has great influence on the endurance life of the product, it is indispensable to carry out stress analysis of the rubber boots. In case of self-contact, various methods for determining contact forces have been suggested with an appropriate contact formulation. Especially, the types of penetration in self-contact are modularized to accelerate conputation with a contact algorithm.

• Journal of Power System Engineering
• /
• v.16 no.5
• /
• pp.40-46
• /
• 2012

For free surface flow problem, a high-order spectral/boundary element method is adapted as an efficient numerical tool. This method is one of the most efficient numerical methods by which the nonlinear gravity waves can be simulated and hydrodynamic forces also can be calculated in time domain. In this method, the velocity potential is expressed as the sum of surface potential and body potential. Then, surface potential is solved by using the high-order spectral method and body potential is solved by using the high-order boundary element method. Using the combination of these two methods, the free surface flow problems of a submerged moving body are solved in time domain. In the present study, lifting surface theory is added to the former work to include effects of lift force. Therefore, a new formulation for the basic mathematical theory is introduced to contain the lift body in calculation.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.2
• /
• pp.183-192
• /
• 2004

New dynamic analysis procedures lot the vertically drilled sea water pile are suggested and demonstrated by the typical design Problem. Pile structure submerged in the sea water as well as forces by the ocean waves and tidal currents are modeled and formulated by finite element method. To obtain wave forces for the finite element equation, Airy's wave theory is tested and selected among others. Lateral lifting forces induced by the vortex shedding of current flow is simply based on the harmonic function with the Strouhal frequency and lifting coefficient. Natural frequencies and frequency responses for the pile are calculated by NASTRAN using the results of the formulation. Dynamic displacement and stress results obtained by these procedures are shown to be applicable to predict the dynamic behaviors of the ocean pile by the wave and lifting forces as a preliminary design analysis.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.1
• /
• pp.24-34
• /
• 1990

In this paper, the formulation of a new simplified finite element is made to analyze the ultimate strength of damaged tubular members subjected to combined axial force and end moment. A damaged tubular member that has the bending deformation and the local dent is modeled by beam elements. Tangent elastic stiffness matrix of a beam element which contains the effect of the geometric nonlinearity is derived by using the updated Lagrangian approach. Here the contribution of the stiffness in the dented area is neglected since its resistance against the external loads is considered to be small. A fully plastic interaction curve of the element under combined loads taking account of the local dent effect is selected as a yielding criterion at each nodal point. Also tangent elasto-plastic stiffness matrix of the element is formulated by plastic node method. Comparison with the present solution and the existing experimental results is made showing that the present method gives quite an accurate solution.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.4
• /
• pp.307-314
• /
• 2009

The present paper briefly describes the space frame element and the fundamental strategies in computational elastic bifurcation theory of geometrically nonlinear, single load parameter conservative elastic spatial structures. A method for large deformation(rotation) analysis of space frame is based on an eulerian formulation, which takes into consideration the effects of large joint translations and rotations with finite deformation(rotation). The local member force-deformation relationships are based on the beam-column approach, and the change in member chord lengths caused by axial strain and flexural bowing are taken into account. and the derived geometric stiffness matrix is unsymmetric because of the fact that finite rotations are not commutative under addition. To detect the singular point such as bifurcation point, an iterative pin-pointing algorithm is proposed. And the path switching mode for bifurcation path is based on the non-negative eigen-value and it's corresponding eigen-vector. Some numerical examples for bifurcation analysis are carried out for a plane frame, plane circular arch and space dome structures are described.

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

The hydrodynamic problem is treated here when a pressurized bag is submerged partially into water and the end points of it oscillate. SES(Surface Effect Ship)has a bag filled with pressurized air at the stern in order to prevent the air leakage, and the pitch motion of SES is largely affected by the hydrodynamic force of the bag. The shape of a bag can be determined with the pressure difference between inside and outside. Once the hydrodynamic pressure is given, the shape of a bag can be obtained, however in order to calculate the hydrodynamic pressure we should know the shape change of the bag, and vice versa. Therefore the type of boundary condition on the surface of a bag is a moving boundary like a free surface boundary. In this paper, the formulation of this problem was done and linearized. The calculation scheme for the radiation problem of an oscillating bag is shown in comparison with the case that the bag is treated as rigid body. The hydrodynamic forces are calculated for various values of the pressure inside the bag and the submerged depth.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.1
• /
• pp.95-110
• /
• 2010

A numerical model to simulate bond-slip behavior of composite beam bridges is introduced in this paper. Assuming a linear bond stress-slip relation along the interface between the slab and girder, the slip behavior is implemented into a finite element formulation. Adopting the introduced model, the slip behavior can be taken account even in a beam element which is composed of both end nodes only. Governing equation of the slip behavior, based on the linear partial interaction theory, can be determined from the force equilibrium and a constant curvature distribution across the section of a composite beam. Since the governing equation for the slip behavior requires the moment values at both end nodes, the piecewise linear distribution of the constant bending moment in an element is assumed. Analysis results by the model are compared with numerical results and experimental values, and load-displacement relations of composite beams were then evaluated to verify the validity of the proposed model.

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

• Geomechanics and Engineering
• /
• v.6 no.6
• /
• pp.517-530
• /
• 2014

Seismic ground response analysis is one of the most important issues in geotechnical earthquake engineering. Conventional seismic site response and free field analysis of layered soils does not consider the effect of surcharge mass which may be present on the top layer. Surcharge mass may develop extra inertial force to the soil and, hence, significantly affect on the results of seismic ground response analysis. Methods of analysis of ground response may also be categorized into time domain and frequency domain concepts. Simplicity in developing analytical relations and accuracy in considering soil dynamic properties dependency to loading frequency are benefits of frequency domain analysis. In this part of the paper, seismic ground response is analyzed using transfer function method for soil layers considering surcharge mass on the top layer. Equation of motion, wave equation, is solved using amended boundary conditions which effectively take the impact of surcharge mass into account. A computer program is developed by MATLAB software based on the solution method developed for wave equation. Layered soils subjected to earthquake loading were numerically studied and solved especially by the computer program developed in this research. Results obtained were compared with those given by DEEP SOIL computer program. Such comparison showed the accuracy of the program developed in this study. Also in this part, the effects of geometrical and mechanical properties of soil layers and especially the impact of surcharge mass on transfer function are investigated using the current approach and the program developed. The efficiency and accuracy of the method developed here is shown through some worked examples and through comparison of the results obtained here with those given by other approaches. Discussions on the results obtained are presented throughout in this part.