• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.1
• /
• pp.65-73
• /
• 1987

An efficient plate bending finite element has been developed using higher-order inplane displacement profiles of the plate. The 6-noded, 21-d.o.f. triangular element including shear deformation effect has been derived from the plate-like continuum by the Galerkin's weighted residual method. Square plate examples were tested with selected element meshes and several aspect ratios for their static behavior under uniformly distributed load. The result of the example tests indicated consistently good performance of the present higher-order plate bending element in comparison with the thin and thick plate solution and other existing finite element solutions.

• Proceedings of the KIEE Conference
• /
• 1992.07b
• /
• pp.649-652
• /
• 1992

ECR(Electron Cyclotron Resonance) occure at ${\omega}_c$=${\omega}$, ${\omega}_c$:electron cycltron frequency, ${\omega}$:electromagnetic wave frequency. ECR system have several merit, 1) power transefer efficiency 2) low neutral gas pressure (below 1 mTorr) 3) high plasma density($10^{12}$ $cm^{-3}$). It is applicated variously in the field of semiconductor and new materials as the manufacturing equipment. Magnetic field in ECR system contruct resonance layer (${\omega}$=2.45GHz, $B_z$=875 Gauss) and control plasma. Plasma is almost generated at resonance layer. If the distance between substrate and resonance layer is short, uniformity of plasma is related with profile of resonance layer. Plasma have the property "Cold in Field", so directonality of magnetic field is one of the control factors of anisotropic etching. In this study, we calculate B field and flux line distribution, optimize geometry and submagnet current and improve of magnetic field directionality (99.9%) near substrate. For the purpose of calculation, vector potential A(r,z) and magnetic field B(r,z), green function and numerical integration is used. Object function for submagnet optimization is magnetic field directionality on the substrate and Powell method is used as optimization skim.

• Computers and Concrete
• /
• v.11 no.1
• /
• pp.51-61
• /
• 2013

This paper describes earthquake response of the Arhavi Highway Tunnel its geometrical properties, 3D finite element model and the linear time history analyses under a huge ground motion considering soil-structure interaction. The Arhavi Highway Tunnel is one of the tallest tunnels constructed in the Black Sea region of Turkey as part of the Coast Road Project. The tunnel has two tubes and each of them is about 1000 m tall. In the study, lineartime history analyses of the tunnel are performed applying north-south, east-west and up accelerations components of 1992 Erzincan, Turkey ground motion. In the time history analyses, Rayleigh damping coefficients are calculated using main natural frequency obtained from modal analysis. Element matrices are computed using the Gauss numerical integration technique. The Newmark method is used in the solution of the equation of motion. Because of needed too much memory for the analyses, the first 10 second of the ground motions, which is the most effective duration, is taken into account in calculations. The results obtained 3D finite element model are presented. In addition, the displacement and stress results are observed to be allowable level of the concrete material during the earthquakes.

• Geomechanics and Engineering
• /
• v.4 no.4
• /
• pp.263-279
• /
• 2012

The present paper deals with the analysis of combined footings resting on geosynthetic reinforced granular fill overlying stone column improved poor soil. An attempt has been made to study the influence of inclusion of geosynthetic layer on the deflection of the footing. The footing has been idealized as a beam having finite flexural rigidity. Granular fill layer has been represented by Pasternak shear layer and stone columns and poor soil have been represented by nonlinear Winkler springs. Nonlinear behavior of granular fill layer, stone columns and the poor soil has been considered by means of hyperbolic stress strain relationships. Governing differential equations for the soil-foundation system have been derived and solution has been obtained employing finite difference scheme by means of iterative Gauss Elimination method. Results of a detailed parametric study have been presented, for a footing supporting typically five columns, in non-dimensional form in respect of deflection with and without geosynthetic inclusion. Geosynthetic layer has been found to significantly reduce the deflection of the footing which has been quantified by means of parametric study.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.3
• /
• pp.599-606
• /
• 1988

Structures of interfacial waves in nearly-horizontal countercurrent stratified air-water flow have been measured by means of a needle contact method. Based upon a statistical analysis for the liquid film distribution, statistical properties of the waves such as mean film thickness, mean wave amplitude and rms value of the wave fluctuation have been calculated. The results show that the film distribution can be described by a Gaussian probability density function for the three-dimensional wave regime. It is also indicated that the mean film thick ness and the rms value of the wave fluctuation increase as gas and liquid flow rates are increased in countercurrent two-phase flow. The dimensionless intensity of the wave fluctuation may be regarded as a function of the Froude number and the dimensionless mean film thickness.

• Research in Mathematical Education
• /
• v.8 no.2
• /
• pp.107-114
• /
• 2004

University teachers of linear algebra often feel annoyed and disarmed when faced with the inability of their students to cope with concepts that they consider to be very simple. Usually, they lay the blame on the impossibility for the students to use geometrical intuition or the lack of practice in basic logic and set theory. J.-L. Dorier [(2002): Teaching Linear Algebra at University. In: T. Li (Ed.), Proceedings of the International Congress of Mathematicians (Beijing: August 20-28, 2002), Vol. III: Invited Lectures (pp. 875-884). Beijing: Higher Education Press] mentioned that the situation could not be improved substantially with the teaching of Cartesian geometry or/and logic and set theory prior to the linear algebra. In East Asian countries, science-orientated mathematics curricula of the high schools consist of calculus with many other materials. To understand differential and integral calculus efficiently or for other reasons, students have to learn a lot of content (and concepts) in linear algebra, such as ordered pairs, n-tuple numbers, planar and spatial coordinates, vectors, polynomials, matrices, etc., from an early age. The content of linear algebra is spread out from grades 7 to 12. When the high school teachers teach the content of linear algebra, however, they do not concern much about the concepts of content. With small effort, teachers can help the students to build concepts of vocabularies and languages of linear algebra.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.1 no.2
• /
• pp.22-31
• /
• 1997

We developed a general purpose program for the analysis of flows in the combustor with recirculating flow regime and simulated the flows. The program uses non-staggered grids based on finite volume method and the primitive variables are cartesian velocities. The combustion model is irreversible one step reaction with infinite chemistry The Favre averaged governing equations are considered and the clipped gaussian distribution is considered as a probability density function of the conserved scalar. We calculated turbulent diffusion flame with recirculating flow regime. Simulation shows two recirculating regions like experimental results. Velocity, turbulent kinetic energy, temperature and concentration distribution in simulation agree well with experimental data.

• Structural Engineering and Mechanics
• /
• v.18 no.6
• /
• pp.731-744
• /
• 2004

Modified virtual crack closure integral (MVCCI) technique has become very popular for computation of strain energy release rate (SERR) and stress intensity factor (SIF) for 2-D crack problems. The objective of this paper is to propose a numerical integration procedure for MVCCI so as to generalize the technique and make its application much wider. This new procedure called as numerically integrated MVCCI (NI-MVCCI) will remove the dependence of MVCCI equations on the type of finite element employed in the basic stress analysis. Numerical studies on fracture analysis of 2-D crack (mode I and II) problems have been conducted by employing 4-noded, 8-noded (regular & quarter-point), 9-noded and 12-noded finite elements. For non-singular (regular) elements at crack tip, NI-MVCCI technique generates the same results as MVCCI, but the advantage for higher order regular and singular elements is that complex equations for MVCCI need not be derived. Gauss numerical integration rule to be employed for 8-noded singular (quarter-point) element for accurate computation of SERR and SIF has been recommended based on the numerical studies.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.10
• /
• pp.2640-2652
• /
• 1994

The stiffness of 6-node isotropic element is stiffer than that of 8-node isotropic element of same configuration. This phenomenon was called 'Relative Stiffness Stiffening Phenomenon'. In this paper, an equation of sampling point modification which correct this phenomenon was derived for the composite plate, as well as an equation for an isotropic plate. The relative stiffness stiffening phenomena of an isotropic plate element could be corrected by modifying Gauss sampling points in the numerical integration of stiffness matrix. This technique could also be successfully applied to the static analyses of composite plate modeled by the 3-dimensional 16-node elements. We predicted theoretical errors of stiffness versus the number of layers that result from the reduction of numerical integration order. These errors coincide very well with the actual errors of stiffness. Therefore, we can choose full integration of reduced integration based upon the permissible error criterion and the number of layers by using the thoretically predicted error.

• Journal of Magnetics
• /
• v.18 no.2
• /
• pp.173-177
• /
• 2013

Cantilever-detected high-frequency electron spin resonance (ESR) is a powerful method of sub-terahertz and terahertz ESR spectroscopy for a tiny magnetic sample at low temperature. In this technique, a small magnetization change associated with ESR transition is detected as deflection of a sample-mounted cantilever. So far, we have succeeded in ESR detection at 370 GHz using a commercial piezoresistive microcantilever. The spin sensitivity was estimated to ${\sim}10^{12}$ spins/gauss. In order to further increase the sensitivity, we adopt a fiber-optic-based detection system using a Fabry-Perot interferometer in place of piezoresistive system. Fabry-Perot cavity is formed between an optical-fiber end and microcantilever surface, and a change in the interference signal, corresponding to the cantilever deflection, is sensitively detected. This system is suitable for low-temperature and high-magnetic-field experiments because of its compact setup and less heat dissipation. In this study, performance of Fabry-Perot interferometer is evaluated, and its application to cantilever-detected ESR measurement is described.