Dynamic fracture behavior of ceramics at elevated temperatures by caustics

The fracture behavior of partially stabilized zirconia (PSZ) and silicon-nitride ceramics (Si₃N₄) is investigated under dynamic loading at elevated temperatures up to 1200°C using the caustic method combined with an ultra high-speed camera. The values of the dynamic fracture toughnessK _Id and the crack-propagation fracture toughnessK _ID are obtained, and it is shown that a dynamic effect on these values is observed in PSZ but not in Si₃N₄. The dynamic crack arrest toughnessK _Ia is found to exist for PSZ.     

A simplified optical method for measuring residual stress by rapid cooling in thermosetting resin strip

This paper presents a simplified optical method for measuring the residual stresses by rapid cooling in thermosetting resin strips. First, the fundamental equations for calculating the residual stress from the residual birefringence were obtained by the linear photoviscoelastic theory. The specimens were then subjected to rapid cooling. After rapid cooling, the residual stress was measured by two methods, the simplified optical method mentioned above and the well-known layer-removal method. The effectiveness of the simplified optical method was discussed by comparing results of the two methods.     

A New Flexible Multibody Beam Element Based on the Absolute Nodal Coordinate Formulation Using the Global Shape Function and the Analytical Mode Shape Function

Several techniques for the reduced dimensionality of finite elementformulations were considered as component mode reduction methods in themiddle sixties. These techniques are widely used in flexiblemultibody simulations for solving small deformation problems. Theabsolute nodal coordinate formulation for solving large rotation anddeformation problems has been established as a full finite elementmethod instead of using similar kinds of reduction techniques. In thispaper, a reduced order absolute nodal coordinate formulation is newlyestablished by introducing the global beam shape function and theanalytical deformation modes as a full finite element. This formulationleads to a constant and symmetric mass matrix as the conventionalabsolute nodal coordinate formulation, and makes it possible to reducethe number of elements and system coordinates of the beam structurewhich undergoes large rotations and large deformations. Numericalexamples show that the excellent agreements between thepresent formulation and the conventional absolute nodal coordinateformulation using a large number of elements are examined. These results demonstratethat the present formulation has high accuracy in the sense that thepresent solutions are similar to the conventional ones with fewersystem coordinates, and high efficiency in computation.     

Numerical simulation of unsteady interaction of centrifugal impeller with its diffuser using Lagrangian discrete vortex method

In this study, an advanced Lagrangian vortex- boundary element method is applied to simulate the unsteady impeller-diffuser interactions in a diffuser pump not only for design but also for off-design considerations. In velocity calculations based on the Biot-Savart law we do not have to grid large portions of the flow field and the calculation points are concentrated in the regions where vorticity is present. Lagrangian representation of the evolving vorticity field is well suited to moving boundaries. An integral pressure equation shows that the pressure distribution can be estimated directly from the instantaneous velocity and vorticity field. The numerical results are compared with the experimental data and the comparisons show that the method used in this study can provide us insight into the complicated unsteady impeller-diffuser interaction phenomena in a diffuser pump.     

Thermal stresses around two parallel cracks in two bonded dissimilar elastic half-planes

Summary Thermal stresses around two parallel cracks in two bonded dissimilar elastic half-planes are determined. One of the cracks lies in the upper half-plane, while the other is in the lower half-plane. Uniform heat flow is assumed to be at right angles to the interface. Application of the Fourier transform technique reduces the problem to that of solving dual integral equations. To solve the equations, the difference of the crack surface temperature and those of the crack surface displacements are expanded in a series of functions which are automatically zero outside the cracks. The unknown coefficients in the series are solved by the Schmidt method. The stress intensity factors are calculated numerially for composite materials featuring a ceramic upper half-plane and a steel lower half-plane.

---

Wärmespannungen um zwei parallele Risse in zwei verbundenen, verschiedenen, elastischen Halbunendlichplatten

Übersicht Es werden die Wärmespannungen um zwei parallele Risse in zwei verbundenen, verschiedenen, elastischen Halbunendlichplatten bestimmt. Einer der beiden Risse liegt in der oberen Halbunendlichplatte, der andere in der unteren. Es wird angenommen, daß ein gleichmäßiger Wärmefluß senkrecht zur Grenzfläche erfolgt. Die Anwendung der Fourier-Transformation reduziert das Problem auf die Lösung dualer Integralgleichungen. Zur Lösung der Gleichungen werden die Temperatur-sowie die Verschiebungsdifferenzen an der Rißoberfläche in eine Reihe von Funktionen entwickelt, die außerhalb der Risse automatisch zu Null werden. Die unbekannten Koeffizienten dieser Reihe werden dann über das Schmidt-Verfahren bestimmt. Anschließend werden für Verbundmaterialien, bei denen die obere Halbunendlichplatte aus Keramik und die untere aus Stahl besteht, die Spannungsintensitätsfaktoren numerisch berechnet.

     

Multi-scale finite element analyses of sheet metals by using SEM-EBSD measured crystallographic RVE models

In this study, two multi-scale analyses codes are newly developed by combining a homogenization algorithm and an elastic/crystalline viscoplastic finite element (FE) method (Nakamachi, E., 1988. A finite element simulation of the sheet metal forming process. Int. J. Numer. Meth. Eng. 25, 283–292; Nakamachi, E., Dong, X., 1996. Elastic/crystalline viscoplastic finite element analysis of dynamic deformation of sheet metal. Int. J. Computer-Aided Eng. Software 13, 308–326; Nakamachi, E., Dong, X., 1997. Study of texture effect on sheet failure in a limit dome height test by using elastic/crystalline viscoplastic finite element analysis. J. Appl. Mech. Trans. ASME(E) 64, 519–524; Nakamachi, E., 1998. Elastic/crystalline viscoplastic finite element modeling based on hardening–softening evaluation equation. In: Proc. of the 6th NUMIFORM, pp. 315–321; Nakamachi, E., Hiraiwa, K., Morimoto, H., Harimoto, M., 2000a. Elastic/crystalline viscoplastic finite element analyses of single- and poly-crystal sheet deformations and their experimental verification. Int. J. Plasticity 16, 1419–1441; Nakamachi, E., Xie, C.L., Harimoto, M., 2000b. Drawability assessment of BCC steel sheet by using elastic/crystalline viscoplastic finite element analyses. Int. J. Mech. Sci. 43, 631–652); (1) a “semi-implicit” finite element (FE) code and (2) a “dynamic explicit” FE code. These were applied to predict the plastic strain induced yield loci and the formability of sheet metal in the macro scale, and simultaneously the crystal texture and hardening evolutions in the micro scale. The isotropic and kinematical hardening laws are employed in the crystalline plasticity constitutive equation. For the multi-scale structure, two-scales are considered. One is a microscopic polycrystal structure and the other a macroscopic elastic plastic continuum. We measure crystal morphologies by using the SEM-EBSD apparatus with a unit of about 3.8 μm voxel, and define a three dimensional (3D) representative volume element (RVE) for the micro polycrystal structure, which satisfy the periodicity condition of crystal orientation distribution. A “micro” finite element modeling technique is newly established to minimize the total number of finite elements in the micro scale. Next, the “semi-implicit” crystallographic homogenization FE code, which employs the SEM-EBSD measured RVE, is applied to the 99.9% pure-iron uni-axial tensile problem to predict the texture evolution and the subsequent yield loci in the various strain paths. These “semi implicit” results reveal that the plastic strain induced anisotropy in the micro and macro levels can be predicted by our FE analyses. The kinematical hardening law leads a distinct plastic strain induced anisotropy. Our “dynamic-explicit” FE code is applied to simulate the limit dome height (LDH) test problem of the mild steel DQSK, the high strength steel HSLA and the aluminum alloy AL6022 sheet metals, which were adopted as the NUMISHEET2005 Benchmark sheet metals (Smith, L.M., Pourboghrat, F., Yoon, J.-W., Stoughton, T.B., 2005. NUMISHEET2005. In: Proc. of 6th Int. Conf. Numerical Simulation of 3D Sheet Metal Forming Processes, PART A and B(Benchmark), pp. 409–451) to estimate formability. The “dynamic explicit” results reveal that the initial crystal orientation distribution has a large affects to a plastic strain induced texture and anisotropic hardening evolutions and sheet formability.     

Velocity statistics along the stagnation line of an axi- symmetric stagnating turbulent flow

 Velocity statistics along the stagnation line of an axi-symmetric wall stagnating turbulent flow are studied experimentally. A low turbulence, uniform air flow from a nozzle type air supply with an exit diameter of 50 mm stagnates at a wall located 50 mm downstream. A flow velocity is set to 3 m/s, 10 mm downstream from the exit of the air supply. Instantaneous values of streamwise and radial velocities are measured by laser-Doppler velocimetry. The turbulence level in the air flow is changed by use of turbulence generator. When the turbulence generator is not installed in the air supply, the mean velocity profile in the streamwise direction fits well with that of a laminar viscous flow with the rms value of velocity fluctuations low near the wall. With the turbulence generator installed, a significant turbulence structure appears near the wall. When the wall is approached, the rms value of velocity fluctuations in the streamwise direction decreases monotonically while the profile of the rms value in the radial direction reaches a maximum near the wall. The increase in the rms value of velocity fluctuations in the radial direction near the wall is attributed to the bi-modal histogram of the fluctuating velocity in the radial direction. Near the wall, the instantaneous stagnation streamline fluctuates and the probability of the mean location of the stagnation point reaches a maximum not at the stagnation line but on a circle around the stagnation line, resulting in the bi-modal histogram. Turbulence statistics, the rms value of velocity fluctuation and the turbulent kinetic energy, can be normalized successfully by similarity parameters based on the strain rate and the reference turbulent kinetic energy introduced by Champion and Libby. Received: 7 April 1995/Accepted: 27 September 1996     

Direct-contact heat transfer to a spherical liquid/vapor two-phase bubble trailing a wake

The paper deals with the heat transfer causing evaporation of liquid drops in a medium of an immiscible, less volatile liquid. Each drop turns into a two-phase bubble, consisting of a growing vapor phase and a reducing liquid phase, which continues to buoy up in the medium. The bubble is modeled as a sphere in which the yet-to-be vaporized liquid spreads over the rear surface while the rest is occupied by the heat-insulating vapor phase. The rear surface to serve as the effective heat transfer area is assumed to be covered with an axisymmetric wake instead of a boundary layer flow. The quasi-steady, overall heat transfer through the wake in the medium and the layer of the yet-to-be vaporized liquid in the bubble is predicted and compared with relevant experimental results.Die Untersuchung befaßt sich mit der durch die Verdampfung von Flüssigkeitstropfen in einem aus unvermischbarer, flüchtiger Flüssigkeit bestehendem, Medium verursachten Wärmeübertragung. Aus jedem Tropfen entsteht eine Zweiphasenblase, bestehend aus einer wachsenden Dampfphase und einer abnehmenden Flüssigkeitsphase, die in das Medium aufsteigt. Die Blase wird anhand eines Kugelmodells betrachtet, an dem sich die noch nicht verdampfte Flüssigkeit über die Rückseite verteilt und die restliche Fläche mit der wärmeisolierenden Dampfphase behaftet ist. Es wird angenommen, daß sich an der Rückseite als wirksamer Wärmeübergangsbereich ein achsensymmetrischer Nachlauf bildet. Es werden die quasi-gleichförmige Wärmeübertragung durch den Nachlauf in das Medium und durch die Schicht Lage der nicht verdampften Flüssigkeit in der Blase mathematisch berechnet und mit relevanten Versuchsergebnissen verglichen.     

Experimental study on a separate type thermosyphon

The passages of vapor flow and the returning liquid flow, are perfectly partitioned in a separate thermosyphon. Therefore the flooding limit can be eliminated, and practicability based on its construction is highly evaluated. In the present work, a container tube made of heat resisting glass, in which an electric heater is inserted, is selected as the heating section of the experimental equipment. Distilled water is used as the working fluid. The influence of the heater type, the diameter of evaporating section and the liquid fill charge on the heat transfer performance have been studied. The larger heat transfer coefficient is achieved in the case of the U type heater, the larger tube diameter and the less liquid fill charge of the evaporating section. The useful correlation equations of the heat transfer coefficient in the evaporator have been derived.

---

Experimentelle Untersuchung an einem Trennstrom-Thermosiphon

Zusammenfassung In einem Trennstrom-Thermosiphon sind die Strömungskanäle für die Dampfströmung und den Kondensatrücklauf vollkommen separiert. Die Flutgrenze kann deshalb unberücksichtigt bleiben. In dieser Untersuchung dient ein temperaturbeständiges Glasrohr mit innenliegender Heizung als Wärmequelle. Arbeitsmedium ist destilliertes Wasser. Untersucht wurden der Einfluß des Heizelementtyps, des Durchmessers der Verdampferstrecke und der Füllmenge auf das Wärmeüber-gangsverhalten. Das U-Typ Heizelement liefert bei größerem Rohrdurchmesser und geringerer Füllmenge in der Verdampferstrecke die höchsten Wärmeübergangskoeffizienten. Für diese wurden Korrelationsgleichungen aufgestellt.