Real-time Computer Simulation of Three Dimensional Elastostatics using the Finite Point Method

Real-time simulation of deformable solids is essential for some applications such as biological organ simulations for surgical simulators. In this work, deformable solids are approximated to be linear elastic, and an easy and straight forward numerical technique, the Finite Point Method (FPM), is used to model three dimensional linear elastostatics. Graphics Processing Unit (GPU) is used to accelerate computations. Results show that the Finite Point Method, together with GPU, can compute three dimensional linear elastostatic responses of solids at rates suitable for real-time graphics, for solids represented by reasonable number of points.

A finite element method for compressible viscous fluid flows

This work presents a mixed three-dimensional finite element formulation for analyzing compressible viscous flows. The formulation is based on the primitive variables velocity, density, temperature and pressure. The goal of this work is to present a `stable' numerical formulation, and, thus, the interpolation functions for the field variables are chosen so as to satisfy the inf-sup conditions. An exact tangent stiffness matrix is derived for the formulation, which ensures a quadratic rate of convergence. The good performance of the proposed strategy is shown in a number of steady-state and transient problems where compressibility effects are important such as high Mach number flows, natural convection, Riemann problems, etc., and also on problems where the fluid can be treated as almost incompressible. Copyright (C) 2010 John Wiley & Sons, Ltd.

Some algorithms for discrete time queues with finite capacity

We consider a discrete time queue with finite capacity and i.i.d. and Markov modulated arrivals, Efficient algorithms are developed to calculate the moments and the distributions of the first time to overflow and the regeneration length, Results are extended to the multiserver queue. Some illustrative numerical examples are provided.

Stability and throughput analysis of unslotted CDMA-ALOHA with finite number of users and code sharing

We consider a system comprising a finite number of nodes, with infinite packet buffers, that use unslotted ALOHA with Code Division Multiple Access (CDMA) to share a channel for transmitting packetised data. We propose a simple model for packet transmission and retransmission at each node, and show that saturation throughput in this model yields a sufficient condition for the stability of the packet buffers; we interpret this as the capacity of the access method. We calculate and compare the capacities of CDMA-ALOHA (with and without code sharing) and TDMA-ALOHA; we also consider carrier sensing and collision detection versions of these protocols. In each case, saturation throughput can be obtained via analysis pf a continuous time Markov chain. Our results show how saturation throughput degrades with code-sharing. Finally, we also present some simulation results for mean packet delay. Our work is motivated by optical CDMA in which "chips" can be optically generated, and hence the achievable chip rate can exceed the achievable TDMA bit rate which is limited by electronics. Code sharing may be useful in the optical CDMA context as it reduces the number of optical correlators at the receivers. Our throughput results help to quantify by how much the CDMA chip rate should exceed the TDMA bit rate so that CDMA-ALOHA yields better capacity than TDMA-ALOHA.

The log delta log t method for the determination of the coefficient of consolidation

The behaviour of saturated soils undergoing consolidation is very complex, It may not follow Terzaghi's theory over the entire consolidation process, Different soils may behave in such a way as to fit into Terzaghi's theory over some specific stages of the consolidation process (percentage of consolidation), This may be one of the reasons for the difficulties faced by the existing curve-fitting procedures in obtaining the coefficient of consolidation, c(v). It has been shown that the slope of the initial linear portion of the theoretical log U-log T curve is constant over a wider range of degree of consolidation, U, when compared with the other methods in use, This initial well-defined straight line in the log U-log T plot intersects the U = 100% line at T = pi/4, which corresponds to U = 88.3%, The proposed log delta-log t method is based on this observation, which gives the value of c(v) through simple graphical construction, In the proposed method, which is more versatile, identification of the characteristic straight lines is very clear; the intersection of these lines is more precise and the method does not depend upon the initial compression for the determination of c(v).

Role of indenter angle on the plastic deformation underneath a sharp indenter and on representative strains: An experimental and numerical study

The subsurface microhardness mapping technique of Chaudhri was utilized to determine the shape, size and distribution of plastic strain underneath conical indenters of varying semi-apex angles, alpha (55 degrees, 65 degrees and 75 degrees). Results show that the elastic-plastic boundary under the indenters is elliptical in nature, contradicting the expanding cavity model, and the ellipticity increases with alpha. The maximum plastic strain immediately under the indenter was found to decrease with increasing alpha. Complementary finite-element analysis was conducted to examine the ability of simulations to capture the experimental observations. A comparison of computational and experimental results indicates that the plastic strain distributions as well as the maximum strains immediately beneath the indenter do not match, suggesting that simulation of sharp indentation requires further detailed studies for complete comprehension. Representative strains, epsilon(r), evaluated as the volume-average strains within the elastic-plastic boundary, decrease with increasing alpha and are in agreement with those estimated by using the dimensional analysis. (C) 2011 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Hidden order in serrated flow of metallic glasses

We report results of statistical and dynamic analysis of the serrated stress-time curves obtained from compressive constant strain-rate tests on two metallic glass samples with different ductility levels in an effort to extract hidden information in the seemingly irregular serrations. Two distinct types of dynamics are detected in these two alloy samples. The stress-strain curve corresponding to the less ductile Zr65Cu15Ni10Al10 alloy is shown to exhibit a finite correlation dimension and a positive Lyapunov exponent, suggesting that the underlying dynamics is chaotic. In contrast, for the more ductile Cu47.5Zr47.5Al5 alloy, the distributions of stress drop magnitudes and their time durations obey a power-law scaling reminiscent of a self-organized critical state. The exponents also satisfy the scaling relation compatible with self-organized criticality. Possible physical mechanisms contributing to the two distinct dynamic regimes are discussed by drawing on the analogy with the serrated yielding of crystalline samples. The analysis, together with some physical reasoning, suggests that plasticity in the less ductile sample can be attributed to stick-slip of a single shear band, while that of the more ductile sample could be attributed to the simultaneous nucleation of a large number of shear bands and their mutual interactions. (C) 2011 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Nonlinear analysis of adhesively bonded lap joints considering viscoplasticity in adhesives

This paper presents nonlinear finite element analysis of adhesively bonded joints considering the elastoviscoplastic constitutive model of the adhesive material and the finite rotation of the joint. Though the adherends have been assumed to be linearly elastic, the yielding of the adhesive is represented by a pressure sensitive modified von Mises yield function. The stress-strain relation of the adhesive is represented by the Ramberg-Osgood relation. Geometric nonlinearity due to finite rotation in the joint is accounted for using the Green-Lagrange strain tensor and the second Piola-Kirchhoff stress tensor in a total Lagrangian formulation. Critical time steps have been calculated based on the eigenvalues of the transition matrices of the viscoplastic model of the adhesive. Stability of the viscoplastic solution and time dependent behaviour of the joints are examined. A parametric study has been carried out with particular reference to peel and shear stress along the interface. Critical zones for failure of joints have been identified. The study is of significance in the design of lap joints as well as on the characterization of adhesive strength. (C) 1999 Elsevier Science Ltd. All rights reserved.

Numerical simulations of fracture initiation in ductile solids under mode I, dynamic loading

In this paper, an overview of some recent computational studies by the authors on ductile crack initiation under mode I, dynamic loading is presented. In these studies, a large deformation finite element procedure is employed along with the viscoplastic version of the Gurson constitutive model that accounts for the micro-mechanical processes of void nucleation, growth and coalescence. A three-point bend fracture specimen subjected to impact, and a single edge notched specimen loaded by a tensile stress pulse are analysed. Several loading rates are simulated by varying the impact speed or the rise time and magnitude of the stress pulse. A simple model involving a semi-circular notch with a pre-nucleated circular hole situated ahead of it is considered. The growth of the hole and its interaction with the notch tip, which leads to plastic strain and porosity localization in the ligament connecting them, is simulated. The role of strain-rate dependence on ductile crack initiation at high loading rates, and the specimen geometry effect on the variation of dynamic fracture toughness with loading rate are Investigated.

Permeability and consolidation behavior of fly ashes

A knowledge of permeability and consolidation is essential in a number of engineering problems such as settlement, seepage, and stability of the structures. Since fly ash is used very widely for several geotechnical applications, there is a need to understand its permeability and consolidation behavior. This paper presents a detailed study conducted on two Indian fly ashes. It brings out the role of chemical composition (free lime) on the permeability and consolidation behavior of fly ashes. It is found that the permeability values computed based on grain-size distribution agree well with those obtained based on test data.

Rapid method of consolidation testing

The determination of consolidation characteristics forms an important aspect in the design of foundations and other earth-retaining structures. The conventional consolidation test as originally proposed by Teaaghi takes considerable time (more than 15 days in highly compressible soils with low coefficient of consolidation) and effort. Any effort to reduce the duration of testing will be desirable from several considerations. In this paper, an attempt has been made to propose a rapid method of consolidation testing. In the proposed method, the next load increment is applied as soon as the necessary time required to identify the percent consolidation is reached and to evaluate the coefficient of consolidation by one of the popular curve-fitting procedures. The rectangular hyperbola method has been used to identify the percent consolidation reached after any load increment, and to determine the coefficient of consolidation, before making the next load increment. The time required to complete the test using the rapid consolidation method could be as low as 4-5 h compared with 1 or 2 weeks in the case of the conventional consolidation test.

Simplified seepage consolidation test for soft sediments

The consistency of very soft sediments prevents the conventional oedometer test from being applied to study their compressibility and permeability characteristics. The hydraulic consolidation test in existence requires sophisticated instrumentation and testing procedures. The present paper proposes a seepage-force-induced consolidation testing procedure for studying the compressibility and permeability behavior of soft sediments at low effective stress levels. The good agreement that has been observed between the results obtained from the proposed method and the conventional oedometer test at overlapping effective stress levels indicates that the proposed method can be used to study the compressibility and permeability characteristics of soft sediments at low effective stress levels satisfactorily.

Stress-strain Curves and Corresponding Structural Parameters in Mulberry and Non-mulberry Silk Fibers

Mulberry fiber (Bivoltine) and non-mulberry fiber (Tassar) were subjected to stress-strain studies and the corresponding samples were examined using wide angle X-ray scattering studies. Here we have two different characteristic stress-strain curves and this has been correlated with changes in crystallite shape ellipsoids in all the fibers. Exclusive crystal structure studies of Tassar fibers show interesting feature of transformation from antiparallel chains to parallel chains.