• The Journal of the Acoustical Society of Korea
• /
• v.31 no.1
• /
• pp.19-28
• /
• 2012

The time-domain solution from the Kirchhoff integral equation for an exterior problem is not unique at certain eigen-frequencies associated with the fictitious internal modes as happening in frequency-domain analysis. One of the solution methods is the CHIEF (Combined Helmholtz Integral Equation Formulation) approach, which is based on employing additional zero-pressure constraints at some interior points inside the body. Although this method has been widely used in frequency-domain boundary element method due to its simplicity, it was not used in time-domain analysis. In this work, the CHIEF approach is formulated appropriately for time-domain acoustic boundary element method by constraining the unknown surface pressure distribution at the current time, which was obtained by setting the pressure at the interior point to be zero considering the shortest retarded time between boundary nodes and interior point. Sound radiation of a pulsating sphere was used as a test example. By applying the CHIEF method, the low-order fictitious modes could be damped down satisfactorily, thus solving the non-uniqueness problem. However, it was observed that the instability due to high-order fictitious modes, which were beyond the effective frequency, was increased.

• Journal of Navigation and Port Research
• /
• v.26 no.4
• /
• pp.473-479
• /
• 2002

The estimation of the volume of a pit excavation is often required in many surveying, soil mechanics, highway applications and transportation engineering situations. The calculation of earthwork plays a major role in plan or design of many civil engineering projects such as seashore reclamation, and thus it has become very important to improve the accuracy of earthwork calculation. In this paper the spot height method, proposed formulas(A, B, C), and chen and Line method are compared with the volumes of the pits in these examples. And we proposed an algorithm of finding a terrain surface with the free boundary conditions and both direction spline method drawback, i.e., the modeling curves form peak points at the joints. To avoid this drawback, the cubic spline polynomial was chosen as the methematical model of the new method. From the characteristics of the cubic spline polynomial, the modeling curve of the new method was smooth and matched the ground profile well. As a result of this study, algorithm of proposed three methods to estimate pit excavation volume provided a better accuracy than spot height, chamber, chen and Lin method. And the mathematical model mentioned makes is thought to give a maximum acccuracy in estimating the volume of a pit excavation.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.60-68
• /
• 1991

The updated Lagrangian Finite Element Method is introduced to analyse rigid body-fluid impact problem which is characterized by incompressible Navier-Stokes equations and impact-contact conditions between free surface and rigid body. For the convenience of numerical computation, velocity fields are splinted into vicous and pressure parts, and then the governing equations and boundary conditions are decomposed in accordance with the decomposition. However, Viscous stresses acting an the solid boundaries are neglected on the assumption that very small velocity gradients may occur during extremely small time interval of the impact. Four coded quadrilateral elements are used to discretize the space domain and the fully explicit time-marching algorithm is employed with a reasonably small time step. At the beginning of each time step, contact velocity of the rigid body is computed from the momentum balance between the body and the fluid. The velocity field is then computed to satisfy the discretized equations of motions and incompressibility and contact constraints as well as an exact free surface boundary condition. At the end of each time step, the fluid domain is updated from the velocity field. In the present time stepping numerical analysis, behaviour of the free surface near the body can be observed without any difficulty which is very important in the water impact problem. The applicability of the algorithm is illustrated by a wedge type falling body problem. The numerical solutions for time-varying pressure distributions and impact loadings acting ion the surface are obtained.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.3
• /
• pp.1-13
• /
• 1998

In the present paper, the hydrodynamic characteristics of three dimensional hydrofoils moving with a constant speed below the free surface using a higher-order boundary element method based on 9-node Lagrangian curvilinear elements are investigated. A bi-quadratic spline scheme is employed to improve the numerical results on the free surface. To validate the present scheme, the calculated results are compared with the analytic solutions for a submerged sphere and a spheroid showing a good agreement. For the validation of the hydrofoil study, the computed lift and drag of a hydrofoil having $NACA64_{1}A412$ section with aspect ratio(A.R.) of 4 are compared with the experimental data by Wadlin et al.[28]. The comparison covers a number of variations of angle of attack and submergence depth. Then, using an A.R. hydrofoil with NACA0012 section, the free surface on the lift and drag are investigated and these are compared with the previous results. The wave elevations and patterns created by the aforementioned submerged bodies are also investigated with Froude numbers and submergences.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.11
• /
• pp.205-212
• /
• 2018

In order to overcome shortcomings of commercial analysis program for solving certain geotechnical problems, finite element method adopting isoparametric quadrilateral element was selected as a tool for analyzing soil behavior and calculating process was programmed. Two examples were considered in order to verify reliability of the developed program. One of the two examples is the case of acting isotropic confining pressure on finite element and the other is the case of acting shear stress on the sides of the finite element. Isoparametric quadrilateral element was considered as the finite element and displacements in the element can be expressed by node displacements and shape functions in the considered element. Calculating process for determining strain which is defined by derivatives using global coordinates was coded using the Jacobian and the natural coordinates. Four point Gauss rule was adopted to convert double integral which defines stiffness of the element into numerical integration. As a result of executing analysis of the finite element under isotropic confining pressure, calculated stress corresponding to four Gauss points and center of the element were equal to the confining pressure. In addition, according to the analyzed results for the element under shear stress, horizontal stresses and vertical stresses were varied with positions in the element and the magnitudes and distribution pattern of the stresses were thought to be rational.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.13 no.3 s.55
• /
• pp.164-172
• /
• 2009

Anchors and saddles are used to have sufficient geometrical rigidity and make target configuration of main cable of suspension bridge. Neglecting the sliding effect at saddles, points at them have been idealized as fixed nodes in lots of former studies. In general, sliding effects are reported to show significant structural behaviors of main cable and cause to the different responses of bridges. During early erection steps of the suspension bridge, especially, the sliding effect occurs easily because there is large difference of cable tension between main and side span in removing set-back ropes or not applying set-back. This study presents the finite element analysis considering to cable sliding effect and shows the comparison of differences between sliding and non-sliding at election sequence. The analysis of sliding between main cables and saddles needs to obtain more realistic responses because the analysis result can represent unfavorably different responses of bridges. Moreover, the sliding analysis method and results in this study can be used to basic criteria in engineering design and construction steps.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.19 no.3 s.73
• /
• pp.283-293
• /
• 2006

Recently, various types of structural systems for skyscrapers are studied as the height and size of the building structures rapidly increase due to social and economical needs. The mega frame system among them, which is the structural system developed recently, is known as a suitable structural system for skyscrapers because this structural system has sufficient stiffness against the lateral forces by combination of mega members which consist of many columns and girders. Since the mega frame structure has significant numbers of elements and nodes, it takes tremendous times and computer memories to analyze and design the structures. Therefore, the exclusive structural analysis program for mega frame system is developed to reduce the efforts and time required for the analysis and design of mega frame structure. To this end, an efficient modelling technique using the characteristics of mega frame structures and an efficient analytical model, which uses a few DOFs selected by the user using the matrix condensation method, are developed in this study. Static and dynamic analyses are conducted using an example structure. The effectiveness and accuracy of the developed program we verified by the comparison between the results of the proposed method and the conventional method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.3
• /
• pp.385-396
• /
• 1999

In the process of structural design for large-scale structures with several thousands of degrees of freedom, a plethora of structural calculations with large amount of data storage are required to obtain the forces and displacements of the members. However, current computational environment with single microprocessor such as a personal computer or a workstation is not capable of generating a high-level of efficiency in structural analysis and design process for large-scale structures. In this paper, a high-performance parallel computing system interconnected by a network of personal computers is proposed for an efficient structural analysis. Two distributed structural analysis algorithms are developed in the form of distributed or parallel preconditioned conjugate gradient (DPCG) method. To enhance the performance of the developed distributed structural analysis algorithms, the number of communications and the size of data to be communicated are minimized. These algorithms are applied to the structural analyses of three large space structures as well as a 144-story tube-in-tube framed structure.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.3
• /
• pp.37-44
• /
• 2004

This paper deals with the performance analysis of LIPCA(Light-weight Piezo-Composite actuator) including the material nun-linearity of the embedded 3203HD PZT wafer. For this analysis, we used a piezo-shell element code based on a nine-node assumed strain shell element formulation. The material non-linearity was implemented in the formulation due to a large observed discrepancy between the measured displacement and the computed actuation displacement based on the linear analysis. An experimentally extracted piezo-strain function of the PZT wafer and incremental formulation were incorporated into the linear finite element code to improve the accuracy of the estimated actuation displacement of the LIPCA. The non-linear piezo-shell program was used to predict the non-linear performance of the LIPCA. The simulated actuation displacement from the non-linear code showed much better agreement with the measured data.

• Geotechnical Engineering
• /
• v.9 no.1
• /
• pp.21-30
• /
• 1993

To know the importance of soil paramters which are used to estimate the deformation and porepressure of soil, the sensitivity for soil parameters in elastic analysis is analyzed. Using the consolidation teat results of several cohesive soils, soil parameters are estimated by back analysis method, and from the parameters the deformations and porepressures of the soil are estimated by elastic analysis, In elastic analysis for soil-deformation and porepressure, the sensitivity for the Young's modulus is large, and the esimation of Young's modulus is more important in pro- portion to the size of stress. Using the measured results during initial short period in small stress, the soil parameters can be correctly estimated by back analysis method. To decrease the iteration number in back analysis and to get the better paramters, the initial measurements in more nodes are required and the more accurate initial measurements are required.