Three-dimensional finite element analysis of composites with coated spherical inclusions

A three-dimensional finite element analysis has been used to determine the internal stresses in a three-phase composite. The stresses have been determined for a variety of interphase properties, the thicknesses of the interphase and the volume fractions of particles. Young's modulus has been calculated from a knowledge of these stresses and the applied deformation. The calculations show that stress distributions in the matrix and the mechanical properties are sensitive to the interphase property in the three-phase composites. The interfacial stresses in the three-dimensional analysis are in agreement with results obtained by an axisymmetric analysis. The predicted bulk modulus in three-dimensional analysis agrees well with the theoretical solution obtained by Qui and Weng, but it presents a great divergence from that in axisymmetric analyses. An investigation indicates that this divergence may be caused by the difference in the unit cell structure between two models. A comparison of the numerically predicted bulk and shear modulus for two-phase composites with the theoretical results indicates that the three-dimensional analysis gives quite satisfactory results.

Numerical simulation of axisymmetric unsteady incompressible flow by a vorticity-velocity method

A new numerical method for solving the axisymmetric unsteady incompressible Navier-Stokes equations using vorticity-velocity variables and a staggered grid is presented. The solution is advanced in time with an explicit two-stage Runge-Kutta method. At each stage a vector Poisson equation for velocity is solved. Some important aspects of staggering of the variable location, divergence-free correction to the velocity held by means of a suitably chosen scalar potential and numerical treatment of the vorticity boundary condition are examined. The axisymmetric spherical Couette flow between two concentric differentially rotating spheres is computed as an initial value problem. Comparison of the computational results using a staggered grid with those using a non-staggered grid shows that the staggered grid is superior to the non-staggered grid. The computed scenario of the transition from zero-vortex to two-vortex flow at moderate Reynolds number agrees with that simulated using a pseudospectral method, thus validating the temporal accuracy of our method.

A new method for evaluation of stress intensities for interface cracks

A new method is presented for calculating the values of K-I and K-II in the elasticity solution at the tip of an interface crack. The method is based on an evaluation of the J-integral by the virtual crack extension method. Expressions for calculating K-I and K-II by using the displacements and the stiffness derivative of the finite element solution and asymptotic crack tip displacements are derived. The method is shown to produce very accurate solutions even with coarse element mesh.

Experimental-study of free-surface oscillations of a liquid bridge by optical diagnostics

A non-contact optical method, consisting of a projecting grating technique for the relative measurement of a surface, and a technique of absolute measurement at a fixed point on the surface, are applied to measure the free surface vibration in a liquid bridge of half floating zone with small typical scale of a few of mm for emphasizing the thermocapillary effect in comparison with the effect of buoyancy. The radii variations in both longitudinal and azimuthal directions are obtained, and, then, the feature of surface wave could be analyzed in detail. The results show that there are values of principal oscillatory frequencies at different positions of free surface. The amplitudes of surface waves in longitudinal and azimuthal directions are several mum and several tenths of mum in order of magnitude. The phase of two-dimensional surface waves is different at different height for fixed cross section or at different azimuthal angle for fixed height. The wave features are discussed for the cases of typical parameter ranges.

Equilibrium of isolated magnetostatic flux structure

Two-dimensional magnetostatic models of flux structure confined in stratified atmosphere are discussed in the present paper. The magnetic field in the flux structure is assumed to be force-free at the first step. Numerical solutions for this nonlinear free boundary problem are obtained by finite element method. Results show clearly the relation between the inside fields and outside pressure, especially the influence of atmospheric pressure distribution on the flux structure.

Free surface vibration in oscillatory convection of half floating zone

Projecting an orthographical grating mask (20pl/mm) on the surface of a small liquid bridge and receiving the reflected distortion image, one can calculate out reversely the shape of free surface of a liquid bridge. In this way we measured the surface shape of a small floating zone and the two-dimensional deformation of its vibration. The mechanism of thermocapillary oscillatory convection and the three-dimensional variation of the free surface are revealed experimentally. The principle for space experiment has been studied in our laboratory.

A finite-element analysis of viscoelastically damped sandwich plates

A finite element analysis associated with an asymptotic solution method for the harmonic flexural vibration of viscoelastically damped unsymmetrical sandwich plates is given. The element formulation is based on generalization of the discrete Kirchhoff theory (DKT) element formulation. The results obtained with the first order approximation of the asymptotic solution presented here are the same as those obtained by means of the modal strain energy (MSE) method. By taking more terms of the asymptotic solution, with successive calculations and use of the Padé approximants method, accuracy can be improved. The finite element computation has been verified by comparison with an analytical exact solution for rectangular plates with simply supported edges. Results for the same plates with clamped edges are also presented.

Research on Bearing Capacity Calculation Method of Large Diameter Cylinder in Soft Clay

A large diameter cylinder inserted in soils is a new type of engineering structures used in offshore and port engineering. The mechanism of its bearing capacity and the analysis of its stability are important to its design and applications. In this paper, the finite element method is used to analyze the reacting forces of the soft soil foundation on the structure under the wave action. A simplified method is proposed, based on the plastic limit method, for the safety and stability analysis. Our analysis shows that the assumptions made in this paper and the mechanism used are reasonable, and the results obtained are appropriate. The calculation method is very efficient and can be used to evaluate main parameters of the structure in its preliminary designs.

Development of a Cell Size Relaxed Scheme for the Direct Simulation Monte Carlo Method

The conventional direct simulation Monte Carlo (DSMC) method has a strong restriction on the cell size because simulated particles are selected randomly within the cell for collisions. Cells with size larger than the molecular mean free path are generally not allowed in correct DSMC simulations. However, the cell-size induced numerical error can be controlled if the gradients of flow properties are properly involved during collisions. In this study, a large cell DSMC scheme is proposed to relax the cell size restriction. The scheme is applied to simulate several test problems and promising results are obtained even when the cell size is greater than 10 mean free paths of gas molecules. However, it is still necessary, of course, that the cell size be small with respect to the flow field structures that must be resolved.

A video method for quantifying size distribution, density, and three-dimensional spatial structure of reef fish spawning aggregations

There is a clear need to develop fisheries independent methods to quantify individual sizes, density, and three dimensional characteristics of reef fish spawning aggregations for use in population assessments and to provide critical baseline data on reproductive life history of exploited populations. We designed, constructed, calibrated, and applied an underwater stereo-video system to estimate individual sizes and three dimensional (3D) positions of Nassau grouper (Epinephelus striatus) at a spawning aggregation site located on a reef promontory on the western edge of Little Cayman Island, Cayman Islands, BWI, on 23 January 2003. The system consists of two free-running camcorders mounted on a meter-long bar and supported by a SCUBA diver. Paired video “stills” were captured, and nose and tail of individual fish observed in the field of view of both cameras were digitized using image analysis software. Conversion of these two dimensional screen coordinates to 3D coordinates was achieved through a matrix inversion algorithm and calibration data. Our estimate of mean total length (58.5 cm, n = 29) was in close agreement with estimated lengths from a hydroacoustic survey and from direct measures of fish size using visual census techniques. We discovered a possible bias in length measures using the video method, most likely arising from some fish orientations that were not perpendicular with respect to the optical axis of the camera system. We observed 40 individuals occupying a volume of 33.3 m3, resulting in a concentration of 1.2 individuals m–3 with a mean (SD) nearest neighbor distance of 70.0 (29.7) cm. We promote the use of roving diver stereo-videography as a method to assess the size distribution, density, and 3D spatial structure of fish spawning aggregations.

A method for estimating the rate of shedding of tags from yellowfin tuna

ENGLISH: The present paper describes a new method for estimating the shedding rate of tags. The method utilizes not only data on tagging and recovery of fish marked with two tags but also data from those marked with one. One important advantage of the new technique is that the estimates of the shedding rates are free from distortion caused by variations in fishing intensity during the total recovery period. The idea of this method appears to be implicit in a short note by Gulland (1963). This technique has been applied to the data obtained by the Inter-American Tropical Tuna Commission in a tagging cruise off the west coast of southern Baja California, during June 1963, at which time both single and double-tagged yellowfin tuna were released. Details of the tagging procedure and equipment have been described by Fink (1965b). The results presented in the present paper are for yellowfin tuna tagged with dart tags. Estimates of shedding should be made separately for each species investigated and also for each type of tag used, since these rates may be variable and often unexpectedly high (Springer and McErlean 1961, Chadwick 1963). SPANISH: El presente estudio describe un nuevo método para estimar las tasas del desprendimiento de marcas. El método emplea no solamente los datos sobre la marcación y recobro de peces marcados con dos marcas, pero también datos de los peces marcados con una marca. Una ventaja importante de la nueva técnica, es que las estimaciones de las tasas de desprendimiento son libres de alteración, causada por las variaciones en la intensidad de pesca durante el período total de recobro. La idea de este método parece ser implícita en un breve apunte por Gulland (1963). Esta técnica se ha aplicado a los datos obtenidos por la Comisión Interamericana del Atún Tropical, en un crucero de marcación efectuado frente a la costa occidental al sur de Baja California, en junio de 1963, tiempo en el cual fueron liberados atunes aleta amarilla marcados tanto con una como con dos marcas. Los detalles del procedimiento de la marcación y del equipo usado han sido descritos por Fink (1965b). Los resultados presentados en este estudio, pertenecen al atún aleta amarilla marcado con marcas de dardo. Las estimaciones del desprendimiento deben efectuarse separadamente para cada especie que ha sido investigada y también para cada tipo de marca usado, ya que estas tasas pueden ser vaiables, y a menudo inesperadamente altas (Springer y McErlean 1961, Chadwick 1963). (PDF contains 20 pages.)

An egg production method for estimating spawning biomass of pelagic fish: Application to the northern anchovy, Engraulis mordax

Fishery scientists engaged in estimating the size of free-swimming populations have never had a technique available to them whereby all the parameters could be estimated from a resource survey and where no parameter values need to be assumed. Recognizing the need for a technique of this kind, the staff of the Coastal Fisheries Resources Division of the Southwest Fisheries Center (SWFC) devised an egg production method for anchovy biomass assessment. Previously, anchovy biomass was estimated by approximate methods derived from a long-time series and anchovy larval abundance, which required about 5 ma of shiptime each year to integrate the area under a seasonal spawning curve. One major assumption used in the larval abundance census method is that there is constant proportionality between larval numbers and spawning biomass. This has now proved to be erroneous. (PDF file contains 105 pages.)

A simple method to simulate shrinkage-induced cracking in cement-based composites by lattice-type modeling

A new numerical procedure is proposed to investigate cracking behaviors induced by mismatch between the matrix phase and aggregates due to matrix shrinkage in cement-based composites. This kind of failure processes is simplified in this investigation as a purely spontaneous mechanical problem, therefore, one main difficulty during simulating the phenomenon lies that no explicit external load serves as the drive to propel development of this physical process. As a result, it is different from classical mechanical problems and seems hard to be solved by using directly the classical finite element method (FEM), a typical kind of "load -> medium -> response" procedures. As a solution, the actual mismatch deformation field is decomposed into two virtual fields, both of which can be obtained by the classical FEM. Then the actual response is obtained by adding together the two virtual displacement fields based on the principle of superposition. Then, critical elements are detected successively by the event-by-event technique. The micro-structure of composites is implemented by employing the generalized beam (GB) lattice model. Numerical examples are given to show the effectiveness of the method, and detailed discussions are conducted on influences of material properties.

A surface energy model and application to mechanical behavior analysis of single crystals at sub-micron scale

Size effects of mechanical behaviors of materials are referred to the variation of the mechanical behavior due to the sample sizes changing from macroscale to micro-/nanoscales. At the micro-/nanoscale, since sample has a relatively high specific surface area (SSA) (ratio of surface area to volume), the surface although it is often neglected at the macroscale, becomes prominent in governing the energy effect, although it is often neglected at the macroscale, becomes prominent in governing the mechanical behavior. In the present research, a continuum model considering the surface energy effect is developed through introducing the surface energy to total potential energy. Simultaneously, a corresponding finite element method is developed. The model is used to analyze the axial equilibrium strain problem for a Cu nanowire at the external loading-free state. As another application of the model, from dimensional analysis, the size effects of uniform compression tests on the microscale cylinder specimens for Ni and Au single crystals are analyzed and compared with experiments in literatures. (C) 2009 Elsevier B.V. All rights reserved.

An Instrumented Indentation Method for Young's Modulus Measurement with Accuracy Estimation

The relationships between indentation responses and Young's modulus of an indented material were investigated by employing dimensional analysis and finite element method. Three representative tip bluntness geometries were introduced to describe the shape of a real Berkovich indenter. It was demonstrated that for each of these bluntness geometries, a set of approximate indentation relationships correlating the ratio of nominal hardness/reduced Young's modulus H (n) /E (r) and the ratio of elastic work/total work W (e)/W can be derived. Consequently, a method for Young's modulus measurement combined with its accuracy estimation was established on basis of these relationships. The effectiveness of this approach was verified by performing nanoindentation tests on S45C carbon steel and 6061 aluminum alloy and microindentation tests on aluminum single crystal, GCr15 bearing steel and fused silica.