Nonaxisymmetric bifurcation behaviours of composite tubes two different materials subjected to uniform radial shrinkage at the external surface have been investigated and compared with those of single tube. The effect of material parameters normalized with respect to those of outer tube upon the bifurcation point and corresponding mode has been clarified. The parameters substantially affect the bifurcation mode with long-wavelength so that the composite tube with low hardening exponent or with high yield stress of inner tube destabilizes the overall deformation of the tube. However surface type bifurcation, short-wavelength mode, shown on the traction-free inner surface is hardly affected by the material parameters. The surface type bifurcation completely depends on the material characteristics of inner tube and the bifurcation point of composite tube almost coincides with the of single tube.

In this paper research result of transient thermal stress analysis of power plant turbine rotors for life prediction under severs operating conditions is presented. Galerkin's recurrence scheme is used for numerical solution of discretized FEM equation of transient heat conduction equation. Boundary conditions for the equation and operating conditions are intensively investigated for accurate life prediction of turbine rotors in operation. A computer program for on-site application is developed and tested. Distribution of thermal stress in turbine rotors during various operating condition is analyzed with the program and it is found that the peak thermal stress appears during cold stage conditions at the first stage of high pressure rotors.

The method of analyzing flexible mechanisms with clearances was studied considering flexibility of beams in the mechanism using finite elements. Both ends of a beam were modeled as free following Dubowsky's impact pair model. Instead some force constraints were imposed at imposed at the connections between adjoining links. Coulomb model has been developed using dry frictions in place of tangential damping forces in the impact pair model and the contact compliance and damping coefficient approximated in a form of root function were used. As examples, impacts of a rigid ball in a cylinder, impact beam model and four-bar mechanisms made up of three flexible links with clearance connections were simulated numerically. The results from examples showed similar but a little bit smaller magnitude of impact forces compared with published studies.

In this paper, a new method is suggested to analyze impulsive stresses at loading poing of concentrated impact load under certain impact conditions determined by impact velocity, stiffness of plate and mass of impact body, etc. The impulsive stresses are analyzed by using the three dimensional dynamic theory of elasticity so as to analytically clarify the generation phenomenon of cone crack at the impact fracture of fragile materials (to be discussed if the second paper). The Lagrange's plate theory and Hertz's law of contact theory are used for the analysis of impact load, and the approximate equation of impact load is suggested to analyze the impulsive stresses at the impact point to decide the ranage of impact load factor. When impact load factors are over and under 0.263, approximate equations are suggested to be F(t)=Aexp(-Bt)sinCt and F(t)=Aexp(-bt) {1-exp(Ct)} respectively. Also, the inverse Laplace transformation is done by using the F.F.T.(fast fourier transform) algorithm. And in order to clarity the validity of stress analysis method, experiments on strain fluctuation at impact point are performed on a supported square glass plate. Finally, these analytical results are shown to be in close agreement with experimental results.

A three-noded, with three degree-of-freedom at each node, in-plane curved beam element is formulated and employed in eigen-analysis of constant curvature beam. The conventional quadratic shape functions used in a three noded C .deg. type curved beam element produce such an undesirable large stiffness that a significant error is introduced in displacements and stresses. These phenomena are called 'Stiffness Locking Phenomena', which result from spurious strain energy due to inappropriate assumptions on independent isoparametric quadratic interpolation functions. Stiffness locking phenomena can be alleviated by using modified interpolation functions which get rid of spurious constraints of conventional interpolation functions. Eigenvalues and their modes as well as displacements and stresses may be locked because they are related to stiffness. Using modified curved beam element in eigenvalue problem of cantilever and arch, the property and performance of modified curved beam element are examined by numerical experimentations. In these eigen-analyses, mass matrices are calculated by using both modified and unmodified curved beam element, are compared with theoretical solutions. These comparisons show that the performance of the modified curved beam element is better than that of the unmodified curved beam element.

In this study, in order to perform dynamic design of machine tools reasonably and effectively, a method was formulated to be applicable to the damped structures connected by joints having elasticity and damping by using substructure synthesis method. And a nonlinear solution method was proposed and it formulates the nonlinear parts by describing functions and uses the reducing transformation matrix by the substructure synthesis method. The results of frequency response analysis of a machine tool, where an NC lathe was partitioned by three parts of spindle, housing and bed-base part and the nonlinearity of bearing parts between spindle and housing was modelled, showed force dependency of the response.

In adaptive control, the lack of persistent and rich excitation causes the estimated parameters to drift, which degrade the performance of the system and may introduces instability to the system in a stochastic environment. To solve the problem of the parameter drift, the concept of single parameter adaptation is presented. For the parameter identification, a priori error is directly used for adaptation error. The structure of the controller is based upon the minimum variance control technique. The stability and robustness analysis is carried out by the sector stability theorem for the second order system. The computer simulation is performed to justify the theoretical analysis for the various cases.

Techniques which are able to predict and control dynamic characteristics, not affecting the vibrational characteristics on the modification of structural design, are being studied. As one of these techniques, experimental modal analysis is widely applied by many researchers. In this study, modal analysis is performed using transfer matrix method by a macro computer. The developed program would estimate the structural modal parameters precisely, and the validity of this program is certified by comparing with the experimental results of .GAMMA A. structure. Estimated modal parameters(natural frequency, vibrational mode, equivalent mass, etc.) are in accord with the experimental results. Also, the optimal location of the additive mass is determined by the evaluation of the vibrational mode and the equivalent mass. The relation between the additive mass and the equivalent mass is specified, and we come to know that the ratio of equivalent mass to additive mass alter linearly within the range of 20%.

Power transmission characteristics and speed ratio-torque load-axial force relationship for a metal V-belt CVT were investigated theoretically and experimentally. In the metal V-belt CVT drive, it was found that the power was transmitted by thrust force and that band tension should be greater than compression force between the metal blocks. The experimental results for the speed ratio-torque load-axial force relationship showed good agreement with those predicted using equations developed from theoretical considerations.

As a basic study to clarify the scoring resistance in lubricated sliding contact, the temperature rise on frictional surface was analyzed by theoretical method and the effects of various factors on the temperature rise were examined. On the basic of the results obtained theoretically, the practical equations to calculate the maximum average temperature of the contact surface were proposed which are applicable to sliding contact. Then, the effects of sliding velocity and oil temperature on the seizure behavior, and the relation between seizure and temperature rise were investigated. The following conclusions are deduced : The maximum average temperature rise and the other bulk temperature. The former is affected by the size of heat supply region and the sliding velocity, the latter is affected by heat transfer coefficient. Without regard to the operating condition such as sliding velocity, oil temperature and operating time at each load-step, the maximum average temperature just before seizure is nearly constant except in the region of lower velocity. Consequently, the maximum average temperature of the contact surface in boundary lubrication is a useful criterion to predict the scoring of sliding contact.

A data structure representing 3-D objects was designed for the personal computer. It is very simple to be used in the personal computers which have small memory and low speed. A homogeneous transformation for developing 3-D objects was derived. Each developing procedure consists of five transformations : one translational-three rotational-one translational. Developing a solid is a creative work. So the results of developing vary with the order of surfaces to be developed. One method to reduce the length of seam was considered. The programs used in this study were written in Pascal and Assembly and a modeller that generates 3-D primitives was included. This program is an interactive dual-screen system. While all the menus in Korean are displayed at the monochrome monitor, the development figures with projective views are drawn at the color monitor. The program has wide applications for plate works.

An interactive computer program to design a mold cavity with the proper rigging system has been developed. In addition to the pattern and the risers generated in part 1 of this work, the various components of the gating system are generated in complete three dimensional models by a rational approach. Then they are laid interactively by the user, and united together with the pattern and the risers to result in the three dimensional model of the mold assembly. Finally, the vents and the mold box are constructed following the user's interactive specification and then the mold cavity is completed in a three dimensional geometric model by subtraction the mold assembly and the vents from the mold box. The three dimensional model of a mold cavity is useful for many related applications such as the solidification simulation for mold evaluation and the NC tool path generation for mold production.

This study investigates the behavior of fatigue crack propagating between holes of holes filled with another materials. When holes of the holes filled with another materials are located symmetrically near a center crack, it is noted that the crack propagation rate is influenced by both the bonding force of the brazing part and the elastic modulus of another material. It is experimentally and analytically confirmed that the center crack stops when its tip reaches near the center line of the holes and a small crack is initiated from the boundaries of holes of the holes filled with another materials and it propagates to final fracture.

A method is proposed which can be used to obtain the fibesr content distribution of compression molded long fiber-reinforced thermoplastic sheet for nonisothermal state. The fiber is modelled to be a sphere. Once the one-dimensional unsteady state heat conduction equation in solved, the mean temperature in defined across the thickness direction. The viscosity of matrix is determined with the mean temperature. Using the obtained viscosity, two-dimensional sheet0like part compression molding is simulated with the finite element method. Comparison with experiments shows that the method accurately predicts the distribution.

This paper is studied to know corrosive wear mechanism of STS304 steel on atmospherical temperature against mating material as the same. The corrosive test was carried out by rubbing the annular surface of two test pieces in distilled water and NaCl aqueous solution. The corrosive wear mechanism was investigated by S.E.M. The experimental results show that there is one Lcr transferring from severe wear to mild wear on change of NaCl concentration and atmospherical temperature, and which is the other still remaining in server wear state. It was found that the critical sliding distance Lcr shorten with increasing NaCl concentration but it is longer with ascending atmospherical temperature and the mild wear state still continues under the condition of high generation rate and elimination rate of the corrosive product. Considering upon the result, the model of corrosive wear mechanism is proposed.

Elastic modulus of unidirectional short fiber composite has theoretically derived with the consideration of Poisson's ratios of matrix and fiber. Unidirectional short fiber composite is modeled as an aggregate of grains developed by Kerner. Under the assumption of extra strain at fiber ends, the strain distribution along the fiber's length is determined, and the elastic modulus is derived from this distribution. For the consideration of effects of Poisson's ratio, Kerner's results for particulate composites are adapted as boundary conditions. The effect of differences in Poisson's ratio of fiber and matrix on elastic modulus is studied. Proposed equation shows a good agreement with experimental data of Halpin and Tock, et al.

Acoustic emission signals due to leak from circular holes of 0.4, 1, 2 and 4mm diameter and rectangular slits of different geometry having the same cross section as 4mm diameter hole was studied both analytically and experimentally. Acoustic emission signals from a wide-band type transducer were transformed to digital signals through a digital oscilloscope, and $V_{rms}$ and frequency spectrum were obtained by processing digital signals. Relationships between acoustic parameters and fluid mechanical parameters were derived analytically. A quadrapole aerodynamic model was applied in the analysis of leak from the circular holes and $V_{rms}$ was found to be proportional to the root square of leak rate through the circular hole. A modified model based on dipole source mechanism and laminar equivalent diameter was applied in the analysis of leak signals from the rectangular slits. In the case of constant pressure, $V_{rms}$ increased as the laminar equivalent diameter of slit decreased. In the case of constant laminar equivalent diameter, however the result was similar to that for leak from the circular hole. The frequency spectra of leak signals shows the same frequency characteristics irrespective of the pressure difference.rence.

The experimental study investigates two aspects of counterflow sensible heat regenerator operation. First, it examines the regenerator performance in periodic steady state operation with spatially nonuniform inlet temperature in one of the fluid stream. Second, the study examines the transient response of a regenerator to a step change in the inlet temperature of one of the fluid streams. The effect of transient inlet temperatures is analyzed in terms of the response of the outlet fluid temperatures to a step change in temperature of one of the inlet fluid streams. The effect of temperature nonuniformities is analyzed in terms of the change of temperature nonuniformities is analyzed in terms of the change in steady state effectiveness due to a circumferential temperature distribution in one of the inlet fluid streams. an experimental analysis has been conducted using a counterflow, parallel passage, and rotary regenerator made from polyethylene film. Efficiencies follow similar trends with increasing matrix to fluid capacity rate ratio for the balanced and symmetric regenerator with nonuniform inlet temperature.

Steady state turbulent airflow and unsteady characteristics of generation, transportation, and recovery behavior of contaminate particles in the simplified 2 dimensional Vertical Laminar Flow (VLF) type clean room was numerically simulated using the low Reynolds number k-over bar.epsilon- turbulent model. Characteristics of airflow in VLF type clean room are greatly affected by the recirculation zone around working surface. The recirculation zone must be considered at the time of clean room design because the recirculation zone whose area increases with increment of inlet velocity exerts bad influence upon the performance of clean room in terms of particle contamination. The location of maximum particle concentration changes from the location of particle source to the recirculation zone, while averaged particle concentration is reduced exponentially with time. Recovery time of clean room with spontaneous particle generation source is inversely proportional to inlet velocity. We introduce nondimensionalized recovery time through the dimensional analysis, which can indicates the general performance of clean room with design structure change. It was identified that .tau. is independent of inlet velocity and background concentration. Therefore .tau. can be the simple factor to compare the different structure of clean room in terms of dynamic response to contamination and becomes larger with better structure of clean room.

Comprehensive numberical computations have been made for four turbulent swirling jets with and without recirculation to critically evaluate the accuracy and universality of several exising turbulence models as well as of the modified k-.epsilon. model proposed in the present study. A numerical scheme based on the full Navier-Stoke equations ha been developed and used for this purpose. Inlet conditions are given by experiments, whenever possible, to minimize the error due to incorrect initial conditions. The standard k-.epsilon. model performs well for the strongly swirling jets with recirculation while it underpredicts the influence of swirl for weakly swirling jets. Rodi's swirl correction and algebraic stress model do not exhibit universality for the swirling jets. The present modified k-.epsilon. model derived from algebraic stress model accounts for anisotropy and streamline curvature effect on turbulence. This model performs consistently better than others for all cases. It may be because these flows have a strong dependence of stresses on the local strain of the mean flow. The predictions of truculence intensities indicate that this model successfully reflect the curvature effect in swirling jets, i.e. the stabilizing and destabilizing effects of swirl on turbulence transport.

The interaction of forced convection-conduction with thermal radiation in laminar boundary layer over a plate fin is studied numerically. The analysis is based on complete solution whereby the heat conduction equation for the fin is solved simultaneously with the conservation equations for mass, momentum and energy in the fluid boundary layer adjacent to the fin. The fluid is a gray medium and diffusion(Rosseland) approximation is used to describe the radiative heat flux in the energy equation. The resulting boundary value problem are convection-conduction parameter N$_{c}$ and radiation-conduction parameter m, Prandtl number Pr. Numerical results are presented for gases with the Prandtl numbers of 0.7 & 5 with values of N$_{c}$ and M ranging from 0 to 10 respectively. The object of this study is to provide the first results on forced convection-radiation interaction in boundary layer flow over a semi-infinite flay plate which can be used for comparisons with future studies that will consider a more accurate expression for the radiative heat flux. The agreement of the results from the complete solution presented by E. M. Sparrow and those from this paper for the special case of M=0 is good. The overall rate of heat transfer from the fin considering radiative effect is higher than that from the fin neglecting radiative effect. The local heat transfer coefficient with radiative effect is higher than that without radiative effect. In the direction from tip to base, those coefficients decrease at first, attain minimum, and then increase. The larger values of N$_{c}$ M, Pr give rise to larger fin temperature variations and the fin temperature without radiative effect is always higher than that with radiative effect.

A rapid expansion of moist air of condensible gas through a supersonic nozzle gives rise to condensation of nonequilibrium and equilibrium processes. Because most of the effects of condensation on the flow are caused by process of nonequilibrium condensation, it is very important to know the onset. condition of nonequilibrium condensation. In the present study, the relation between the initial relative stagnation humidity and the onset Mach number, for the case of the similarity law suggested by Zierep and Lin. Furthermore, the present theoretical result number is compared with the experimental, numerical and other results.

Turbulent jet in a crossflow, issuing from a row of holes on a convex surface of 90 .deg. bend duct, is predicted by a 3-dimensional numerical method. The Cartesian coordinate system in adopted in upstream and downstream tangents and the cylindrical polar coordinate system in curved region. The Reynolds stresses and heat fluxes are obtained from a standard k-e model in the core region and van Driest model in the vicinity of the wall. The governing equations are discretized by a finite volume method and solutions are obtained by a locally elliptic calculation procedure. Pressure and convective terms are treated by SIMPLE algorithm and hybrid scheme respectively. A vortex initially induced by the injected jet has been built up due to the interaction with the secondary flow caused by pressure gradient and centrifugal force. The vortex structure has a strong influence on the wall cooling effectiveness. Another vortex like horseshoe is formed in the vicinity of the injection hole and its strength is getting weak as it moves downward.

A heat-flux meter has been designed and manufactured to measure the heat flux from the flame. A calibration method of the heat-flux meter by a calibration furnace has also been proposed. The k-type (Chromel-Alumel) thermocouple material has been used as the material for the beat-flux meter. The electormotive force (e.m.f.) from the K-type thermocouple is shown to be linearly proportional to the heat flux absorbed. The characteristics of the heat-flux meter become better as the radius of heat absorbing disk becomes larger and its thickness thinner.

It was known that the fracture incidences of offshore structure were mostly originated from the surface defects. Especially, in the case of the welded structures, since the welded region has some defects and incomplete beads which are apt to behave like the surface cracks, it has been necessary to evaluate the environmental effects on crack growth at HAZ for the design crack growth behaviour at multi-pall HAZ for SWS41 steel under free corrosion and cathodic protection(-0.9V vs Ag/Agcl) conditions. The results are summarized as follows ; (1) Crack growth rate of the as weld in air was faster than that of the parent and PWHT specimens over all .DELTA.K rang. (2) In free corrosion test, surface crack growth rate of the as welded was decreased in comparison with that of the parents. (3) In fatigue test under cathodic protection, cathodic electric potential(-0.9V vs Ag/Agcl) for the SWS41 steel parent was effective, while for the as welded ineffective. (4) There was a tendency that the exponent(m) of the Paris' equation was decreased in order of microhardness magnititude in air and under cathodic protection conditions and vise versa in free corrosion. (5) Fracture surface has dimples and ductile striations in air test, but transgranular cracks and brittle striations under cathodic protection test.