• Earthquakes and Structures
• /
• v.5 no.2
• /
• pp.161-205
• /
• 2013

We study seismically induced, anti-plane strain wave motion in a non-homogeneous geological region containing tunnels. Two different scenarios are considered: (a) The first models two tunnels in a finite geological region embedded within a laterally inhomogeneous, layered geological profile containing a seismic source. For this case, labelled as the first boundary-value problem (BVP 1), an efficient hybrid technique comprising the finite difference method (FDM) and the boundary element method (BEM) is developed and applied. Since the later method is based on the frequency-dependent fundamental solution of elastodynamics, the hybrid technique is defined in the frequency domain. Then, an inverse fast Fourier transformation (FFT) is used to recover time histories; (b) The second models a finite region with two tunnels, is embedded in a homogeneous half-plane, and is subjected to incident, time-harmonic SH-waves. This case, labelled as the second boundary-value problem (BVP 2), considers complex soil properties such as anisotropy, continuous inhomogeneity and poroelasticity. The computational approach is now the BEM alone, since solution of the surrounding half plane by the FDM is unnecessary. In sum, the hybrid FDM-BEM technique is able to quantify dependence of the signals that develop at the free surface to the following key parameters: seismic source properties and heterogeneous structure of the wave path (the FDM component) and near-surface geological deposits containing discontinuities in the form of tunnels (the BEM component). Finally, the hybrid technique is used for evaluating the seismic wave field that develops within a key geological cross-section of the Metro construction project in Thessaloniki, Greece, which includes the important Roman-era historical monument of Rotunda dating from the 3rd century A.D.

• Journal of the Korean Operations Research and Management Science Society
• /
• v.7 no.1
• /
• pp.11-16
• /
• 1982

A method is described for the solution of a linearly constrained program with parametric nonlinear objective function. The algorithm proposed in this paper may be regarded as an extension of the simplex method for parametric linear programming. Namely, it specifies the basis at each stage such that feasibility ana optimality of the original problem are satisfied by the optimal solution of the reduced parametric problem involving only nonbasic variables. It is shown that under appropriate assumptions the algorithm is finite. Parametric procedures are also indicated for solving each reduced parametric problem by maintaining the Kuhn-Tucker conditions as the parameter value varies.

• Journal of the Korean Data and Information Science Society
• /
• v.5 no.1
• /
• pp.11-20
• /
• 1994

In recent years, nonparametric kernel estimation of regresion function are abundant and widely applicable to many areas of statistics. Most of modern researches concerned with the fixed global bandwidth selection which can be used in the estimation of regression function with all the same value for all x. In this paper, we propose a method for selecting locally varing bandwidth based on bootstrap method in kernel estimation of fixed design regression. Performance of proposed bandwidth selection method for finite sample case is conducted via Monte Carlo simulation study.

• Journal of Theoretical and Applied Mechanics
• /
• v.3 no.1
• /
• pp.45-65
• /
• 2002

A constitutive model on oorthotropic thermo-elasto-viscoplasticity for fiber-reinforced composite materials Is illustrated, and their thermomechanical responses are predicted with the fully-coupled finite element formulation. The unmixing-mixing scheme can be adopted with the multipartite matrix method as the constitutive model. Basic assumptions based upon the composite micromechanics are postulated, and the strain components of thermal expansion due to temperature change are included In the formulation. Also. more than two sets of mechanical variables, which represent the deformation states of multipartite matrix can be introduced arbitrarily. In particular, the unmixing-mixing scheme can be used with any well-known isotropic viscoplastic theory of the matrix material. The scheme unnecessitates the complex processes for developing an orthotropic viscoplastic theory. The governing equations based on fully-coupled thermomechanics are derived with constitutive arrangement by the unmixing-mixing concept. By considering some auxiliary conditions, the Initial-boundary value problem Is completely set up. As a tool of numerical analyses, the finite element method Is used with isoparametric Interpolation fer the displacement and the temperature fields. The equation of mutton and the energy conservation equation are spatially discretized, and then the time marching techniques such as the Newmark method and the Crank-Nicolson technique are applied. To solve the ultimate nonlinear simultaneous equations, a successive iteration algorithm is constructed with subincrementing technique. As a numerical study, a series of analyses are performed with the main focus on the thermomechanical coupling effect in composite materials. The progress of viscoplastic deformation, the stress-strain relation, and the temperature History are careful1y examined when composite laminates are subjected to repeated cyclic loading.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.1
• /
• pp.61-65
• /
• 2016

Composites are materials that are widely used in industries such as automobile and aircraft. The composite material is required as a material for using in a high temperature environment as well as acting as a high pressure environment like the nozzle part of the ship. It is important to know the properties of composites. Result values obtained substituting the properties of matrix and fiber numerically have an large error compared with experimental value. In this study we utilize CASADsolver EDISON program for using Finite Element Method. Properties by substituting the fiber and Matrix properties of the composite material properties were compared with those measured in the experiment and calculated by the empirical properties.

• Transactions of The Korea Fluid Power Systems Society
• /
• v.8 no.1
• /
• pp.1-9
• /
• 2011

The main design factors which effect on operating speed of linear actuator for valve operation are mass of plunger, electromagnetic motive force, inductance, and return spring, and these factors are not independent but related with each other in view point of design and electromagnetic theory. It is impossible to increase the operating speed by only change the value of any one design factor. The change of any one value results in change of any value related it in various design factors. This paper presents a speed increasing method of linear actuator using a solenoid design method by some governing equations which are composed of electromagnetic theory and empirical knowledge and permanent magnets as assistant material, and proved the propriety by experiments.

• Publications of The Korean Astronomical Society
• /
• v.15 no.spc2
• /
• pp.65-71
• /
• 2000

We have investigated the numerical methods to calculate model atmosphere for the analysis of spectral lines emitted from the sun and stars. Basic equations used in our calculations are radiative transfer, statistical equilibrium and charge-particle conservations. Transfer equation has been solved to get emitting spectral line profile as an initial value problem using Adams-Bashforth-Moulton method with accuracy as high as 12th order. And we have calculated above non linear differential equations simultaneously as a boundary value problem by finite difference method of 3 points approximation through Feautrier elimination scheme. It is found that all computing programs coded by above numerical methods work successfully for our model atmosphere.

• Transactions of the Korean Society of Machine Tool Engineers
• /
• v.10 no.5
• /
• pp.85-95
• /
• 2001

This study proposed a technique for inverse problem, linear approximation of contact position and loading in single and double meshing of spiral bevel gear , using 2-dimension model considered near the tooth by root stress. Determine root stress is carried out far the gear tooth by finite element method and boundary element method. Boundary element discretization near contact point is carefully performed to keep high computational accuracy. And from those estimated results, the comparing estimate value with boundary element method value was discussed.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.369-376
• /
• 2002

The stochastic analysis of semi-infinite domain is presented using the weighted integral method, which is improved to include the higher order terms in expanding the displacement vector. To improve the weighted integral method, the Lagrangian remainder is taken into account in the expansion of the status variable with respect to the mean value of the random variables. In the resulting formulae only the 'proportionality coefficients' are introduced in the resulting equation, therefore no additional computation time and memory requirement is needed. The equations are applied in analyzing the semi-infinite domain. The results obtained by the improved weighted integral method are reasonable and are in good agreement with those of the Monte Carlo simulation. To model the semi-infinite domain, the Bettess's infinite element is adopted, where the theoretical decomposition of the strain-displacement matrix to calculate the deviatoric stiffness of the semi-infinite domains is introduced. The calculated value of mean and the covariance of the displacement are revealed to be larger than those given by the finite domain assumptions which is thought to be rational and should be considered in the design of structures on semi-infinite domains.

• Structural Engineering and Mechanics
• /
• v.68 no.6
• /
• pp.713-720
• /
• 2018

In this paper a new numerical method to determine the elastic curve of the simply supported beams of variable cross-section is demonstrated. In general case it needs to solve linear or small nonlinear second order differential equations with prescribed boundary conditions. For numerical solution the initial values of the slope and the deflection of the end cross-section of the beam is necessary. For obtaining the initial values a lively procedure is developed: it is a special application of the shooting method because boundary value problems can be transformed into initial value problems. As a result of these transformations the initial values of the differential equations are obtained with high accuracy. Procedure is applied for calculating of elastic curve of a simply supported beam of variable cross-section. Results of these numerical procedures, analytical solution of the linearized version and finite element method are compared. It is proved that the suggested procedure yields technically accurate results.