• Journal of Ocean Engineering and Technology
• /
• v.27 no.1
• /
• pp.85-92
• /
• 2013

This paper deals with the resonant type of a WEC (wave energy converter) and the determination method of its geometric parameters which were obtained to construct the robust and optimal structure, respectively. In detail, the optimization problem is formulated with the constraints composed of the response surfaces which stand for the resonance period(heave, pitch) and the meta center height of the buoy. Use of a signal-to-noise ratio calculated from normalized multi-objective results with the weight factor can help to select the robust design level. In order to get the sample data set, the motion responses of the power buoy were analyzed using the BEM (boundary element method)-based commercial code. Also, the optimization result is compared with a robust design for a feasibility study. Finally, the power efficiency of the WEC with the optimum design variables is estimated as the captured wave ratio resulting from absorbed power which mainly related to PTO (power take off) damping. It could be said that the resultant of the WEC design is the economical optimal design which satisfy the given constraints.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2003.11a
• /
• pp.244-251
• /
• 2003

A two-dimensional BEM code, $FRACOD^{2D}$, was applied to simulate fracture initiation and propagation processes in a rock pillar during an in situ heater test of a rock pillar planned at the $\"{A}sp\"{o}$ Underground Rock laboratory of SKB, in Southern Sweden. To take the advantage of conventional BEM for simulating fracturing processes, but without efforts for domain integral transformation, a hybrid approach is developed to simulate the fracturing processes in rock pillar under coupled thermo-mechanical loading. The code FRACOD was used for simulating the fracture initiation and propagation processes with its boundary tractions reflecting the effects of the initial and redistributed thermomechanical stresses in the domain of interest at multiple excavation and heating steps were produced by a special algorithm of stress inversion, based on resultant thermo-mechanical stress fields at each excavation and heat loading step by a FEM code without considering fracturing processes. This hybrid approach can take the advantages of both types of numerical methods and avoids their shortcomings for fracturing process simulation and domain effects, respectively. In this paper, we present the hybrid approach for the stress, displacements, and fracturing processes at sequential excavation and heating steps of the in situ heater test as a predictive modelling, the formulation of the fracturing models and the predictive results. Two sections of borehole depth, 0.5 m and 1.5 m below the tunnel floor are considered. The pillar area is modelled with the FRACOD and the stress field produced by excavation and heating is transferred with corresponding boundary stresses. From the modelling results, the degree of fracturing and damage are evaluated for 120 days of heating. Dominated shear fracturing in the vicinity of the central pillar was observed from the models at both sections, but spalled area appears to be limited. Based on the modelling results, a sensitivity study for the effect of pre-existing fractures in the vicinity of the holes is also conducted, and the initiation and evolution of EDZ around the deposition holes are investigated using this particular numerical technique.

• The Journal of the Acoustical Society of Korea
• /
• v.20 no.6
• /
• pp.94-103
• /
• 2001

The goal of using numerical methods in this study is two-fold: to replicate a set of measured, individualized HRTFs by a computer simulation, and also to visualise the resultant sound field around the head. Two methods can be wed: the Boundary Element Method (BEM) and the Infinite-Finite Element Method (IFEM). This paper presents the results of a preliminary study carried out on a KEMAR dummy-head, the geometry of which was captured with a high accuracy 3-D laser scanner and digitiser. The scanned computer model was converted to a few valid BEM and IFEM meshes with different polygon resolutions, enabling us to optimise the simulation for different frequency ranges. The results show a good agreement between simulations and measurements of the sound pressure at the blocked ear-canal of the dummy-head. The principle of reciprocity provides an effect method to simulate HRTF database. The BEM was also used to investigate the total sound field around the head, providing a tool to visualise the sound field for different arrangements of virtual acoustic imaging systems.

• Steel and Composite Structures
• /
• v.19 no.2
• /
• pp.369-384
• /
• 2015

This work presents a free vibration analysis of functionally graded metal-ceramic (FG) beams with considering porosities that may possibly occur inside the functionally graded materials (FGMs) during their fabrication. For this purpose, a simple displacement field based on higher order shear deformation theory is implemented. The proposed theory is based on the assumption that the transverse displacements consist of bending and shear components in which the bending components do not contribute toward shear forces and, likewise, the shear components do not contribute toward bending moments. The most interesting feature of this theory is that it accounts for a quadratic variation of the transverse shear strains across the thickness, and satisfies the zero traction boundary conditions on the top and bottom surfaces of the beam without using shear correction factors. In addition, it has strong similarities with Euler-Bernoulli beam theory in some aspects such as equations of motion, boundary conditions, and stress resultant expressions. The rule of mixture is modified to describe and approximate material properties of the FG beams with porosity phases. By employing the Hamilton's principle, governing equations of motion for coupled axial-shear-flexural response are determined. The validity of the present theory is investigated by comparing some of the present results with those of the first-order and the other higher-order theories reported in the literature. Illustrative examples are given also to show the effects of varying gradients, porosity volume fraction, aspect ratios, and thickness to length ratios on the free vibration of the FG beams.

• Geomechanics and Engineering
• /
• v.30 no.6
• /
• pp.525-538
• /
• 2022

This paper presents a more general model for T-shaped combined footings that support two columns aligned on a longitudinal axis and each column provides an axial load and two orthogonal moments. This model can be applied to the following conditions: (1) without restrictions on its sides, (2) a restricted side and (3) two opposite sides restricted. This model considers the linear soil pressure. The recently published works have been developed for a restricted side and for two opposite sides restricted by Luévanos-Rojas et al. (2018a, b). The current model considers the uniform pressure distribution because the position of the resultant force coincides with the center of gravity of the surface of the footing in contact with the soil in direction of the longitudinal axis where the columns are located. This paper shows three numerical examples. Example 1 is for a T-shaped combined footing with a limited side (one column is located on the property boundary). Example 2 is for a T-shaped combined footing with two limited opposite sides (the two columns are located on the property boundary). Example 3 is for a T-shaped combined footing with two limited opposite sides, one column is located in the center of the width of the upper flange (b₁/2=L₁), and other column is located at a distance half the width of the strip from the free end of the footing (b₂/2=b-L₁-L). The main advantage of this work over other works is that this model can be applied to T-shaped combined footings without restrictions on its sides, a restricted side and two opposite sides restricted. It also shows the deficiencies of the current model over the new model.

• Structural Engineering and Mechanics
• /
• v.58 no.3
• /
• pp.397-422
• /
• 2016

In the present work, a simple first-order shear deformation theory is developed and validated for a variety of numerical examples of the thermal buckling response of functionally graded sandwich plates with various boundary conditions. Contrary to the conventional first-order shear deformation theory, the present first-order shear deformation theory involves only four unknowns and has strong similarities with the classical plate theory in many aspects such as governing equations of motion, and stress resultant expressions. Material properties and thermal expansion coefficient of the sandwich plate faces are assumed to be graded in the thickness direction according to a simple power-law distribution in terms of the volume fractions of the constituents. The core layer is still homogeneous and made of an isotropic material. The thermal loads are considered as uniform, linear and non-linear temperature rises within the thickness direction. The results reveal that the volume fraction index, loading type and functionally graded layers thickness have significant influence on the thermal buckling of functionally graded sandwich plates. Moreover, numerical results prove that the present simple first-order shear deformation theory can achieve the same accuracy of the existing conventional first-order shear deformation theory which has more number of unknowns.

• Journal of Advanced Marine Engineering and Technology
• /
• v.8 no.1
• /
• pp.85-99
• /
• 1984

The stability of the two-phase flow in a heated channel is of great importance in the design and operation of the boilers and light water nuclear reactors, because it can cause flow oscillations and lead to a violation of thermal limits with resultant overheating of the channels and cladding. This paper presents a systematic evaluation to the variation effects of the basic four (4) dimensionless parameters in a homogeneous equilibrium model. The flow stability is examined on the ground of static characteristic curves. The complicated transfer function of flow dynamics which gives consideration to the transport lag of density wave is derived, and the transient flow stability is analysed by applying the Nyquist stability criterion in control engineering. The analysis results summed up as follows 1. The coolant flow becomes stable in large friction number and specific flow, while it is unstabale in small friction number and flow. 2. Large phase-change number and Froude number destabilize the two-phase flow, but small numbers stabilize it. The effect to variation of phase-change number is more dominant compared with Froude number. 3. The dynamic analysis is required to hold the sufficient safety of heated channels since only static results does not keep it. The special attention could be payed in the design and operation of heat engines, because the unstaable region exists within the stable boundary at small and middle phase-change number and Froude number.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.7
• /
• pp.1391-1397
• /
• 1992

Improved power law flow model was suggested for the calculation of wake flow characteristics around the three dimensional ship stern in case of the formation of bilge vortex in the direction of stern. In comparison with the power law and Coles flow model, the flow velocity calculated based on this study was delayed around the boundary of inner layer and outer layer in reverse flow. More accurate results was obtained with this improved power law flow model by the velocity calculation around ship stern. Accuracy was validated with the comparison of other calculation results and experimental datas.

• Journal of Institute of Control, Robotics and Systems
• /
• v.14 no.11
• /
• pp.1146-1154
• /
• 2008

This paper presents a method for the acceleration analysis of multi-legged walking robots in consideration of the frictional ground contact. This method is based on both unified dynamic equation for finding the acceleration of a robot's body and constraint equation for satisfying no-slip condition. After the dynamic equation representing relationship between actuator torques and body acceleration, is derived from the force and acceleration relationship between foot and body's gravity center, the constraint equation is formulated to reconfigure the maximum torque boundaries satisfying no-slip condition from given original actuator torque boundaries. From application of the reconfigured torques to the dynamic equation, interested acceleration boundaries are obtained. The approach based on above two equations, is adapted to the changes of degree-of-freedoms of legs as well as friction of ground. And the method provides the maximum translational and rotational acceleration boundaries of body's center that are achievable in every direction without occurring slipping at the contact points or saturating all actuators. Given the torque limits in infinite normsense, the resultant accelerations are derived as a polytope. From the proposed method, we obtained achievable acceleration boundaries of 4-legged and 6-legged walking robot system successfully.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.8 no.3
• /
• pp.224-231
• /
• 1984

Having found by means of a tension-shear test an optimal spot welding condition under which the maximum weld strength is to be brought forth, this study made an examination of behavior of fracture concerned with behavior of stress distribution, observed around the nugget periphery of the specimens prepared under the optimal conditions, with one point spot welded mild steel sheets. The resultant findings are as follows: (1)There remarkably exists an optimal spot welding condition to indicate the maximum weld strength, and fracture of the specimens spot welded under that condition occurs outside the nugget boundary. (2)An experiment on the basis of a photoelastic model reveals that the maximum stress is distributed along the center line of the steel plate width but occurs on the region corresponding to heat affected zone of spot welds. (3)Heat affected zone of spot welds consists of coarse grains with considerably low micro Vickers hardness value and of fine grains of high micro Vickers hardness value, and in this unbalanced structure weak region are represented in coarse grain region, where fracture is initiated and continues its propagation.