On numerical implementation of an isotropic elastic-viscoplastic constitutive model for bulk metallic glasses

An efficient algorithm within the finite deformation framework is developed for finite element implementation of a recently proposed isotropic, Mohr-Coulomb type material model, which captures the elastic-viscoplastic, pressure sensitive and plastically dilatant response of bulk metallic glasses. The constitutive equations are first reformulated and implemented using an implicit numerical integration procedure based on the backward Euler method. The resulting system of nonlinear algebraic equations is solved by the Newton-Raphson procedure. This is achieved by developing the principal space return mapping technique for the present model which involves simultaneous shearing and dilatation on multiple potential slip systems. The complete stress update algorithm is presented and the expressions for viscoplastic consistent tangent moduli are derived. The stress update scheme and the viscoplastic consistent tangent are implemented in the commercial finite element code ABAQUS/Standard. The accuracy and performance of the numerical implementation are verified by considering several benchmark examples, which includes a simulation of multiple shear bands in a 3D prismatic bar under uniaxial compression.

Influence of an additional discrete elastic support on the stability of a Leipholz column

The stability characteristics of a conservatively loaded structure are expected to improve if additional supports are provided to the structure. The same, however, may not be said of a non-conservatively loaded structure; several factors, such as the location and stiffness of supports, type of structure and loading, have a significant influence on the stability characteristics. The influence of an arbitrarily located elastic support on the stability characteristics of a Leipholz column is examined.

Analysis of reflection characteristics of a normal incidence plane wave on resonant sound absorbers: A finite element approach

Resonant sound absorbers are used widely as anechoic coatings in underwater applications. In this paper a finite element scheme based on the Galerkin technique is used to analyze the reflection characteristics of the resonant absorber when insonified by a normal incidence plane wave. A waveguide theory coupled with an impedance matching condition in the fluid is used to model the problem. It is shown in this paper that the fluid medium encompassing the absorber can be modeled as an elastic medium with equivalent Lamé constants. Quarter symmetry conditions within the periodic unit cell are exploited. The finite element results are compared with analytical results, and with results published elsewhere in the literature. It is shown in the process that meshing of the fluid domain can be obviated if the transmission coefficients or reflection coefficients only are desired as is often the case. Finally, some design curves for thin resonant absorbers with water closure are presented in this paper.

Dynamic response of a beam on elastic foundation of finite depth under a moving force

In this paper, dynamic response of an infinitely long beam resting on a foundation of finite depth, under a moving force is studied. The effect of foundation inertia is included in the analysis by modelling the foundation as a series of closely spaced axially vibrating rods of finite depth, fixed at the bottom and connected to the beam at the top. Viscous damping in the beam and foundation is included in the analysis. Steady state response of the beam-foundation system is obtained. Detailed numerical results are presented to study the effect of various parameters such as foundation mass, velocity of the moving load, damping and axial force on the beam. It is shown that foundation inertia can considerably reduce the critical velocity and can also amplify the beam response.

Examination of Stress-Fields In Elastic Plastic Indentation

A block of high-purity copper was indented by a 120-degrees diamond-tipped cone. Strain gauges were placed on the surface to measure the radial strains at different surface locations, during loading as well as unloading. The competence of three stress fields proposed for elastic-plastic indentation is assessed by comparing the predicted surface radial strains with those experimentally observed.

Finite element estimation of elastic interlaminar stresses in laminates

Low interlaminar strength and the consequent possibility of interlaminar failures in composite laminates demand an examination of interlaminar stresses and/or strains to ensure their satisfactory performance. As a first approximation, these stresses can be obtained from thickness-wise integration of ply equilibrium equations using in-plane stresses from the classical laminated plate theory. Implementation of this approach in the finite element form requires evaluation of third and fourth order derivatives of the displacement functions in an element. Hence, a high precision element developed by Jayachandrabose and Kirkhope (1985) is used here and the required derivatives are obtained in two ways. (i) from direct differentiation of element shape functions; and (ii) by adapting a finite difference technique applied to the nodal strains and curvatures obtained from the finite element analysis. Numerical results obtained for a three-layered symmetric and a two-layered asymmetric laminate show that the second scheme is quite effective compared to the first scheme particularly for the case of asymmetric laminates.

Coupled electrostatic-elastic analysis for topology optimization using material interpolation

In this paper, we present a novel analytical formulation for the coupled partial differential equations governing electrostatically actuated constrained elastic structures of inhomogeneous material composition. We also present a computationally efficient numerical framework for solving the coupled equations over a reference domain with a fixed finite-element mesh. This serves two purposes: (i) a series of problems with varying geometries and piece-wise homogeneous and/or inhomogeneous material distribution can be solved with a single pre-processing step, (ii) topology optimization methods can be easily implemented by interpolating the material at each point in the reference domain from a void to a dielectric or a conductor. This is attained by considering the steady-state electrical current conduction equation with a `leaky capacitor' model instead of the usual electrostatic equation. This formulation is amenable for both static and transient problems in the elastic domain coupled with the quasi-electrostatic electric field. The procedure is numerically implemented on the COMSOL Multiphysics (R) platform using the weak variational form of the governing equations. Examples have been presented to show the accuracy and versatility of the scheme. The accuracy of the scheme is validated for the special case of piece-wise homogeneous material in the limit of the leaky-capacitor model approaching the ideal case.

Glasslike elastic properties in the omega-beta alloys

We have measured the internal friction and speed of sound in several polycrystalline alloys, using compound torsional oscillators at frequencies between 60 kHz and 100 kHz and temperatures between 50 mK and 100 K. By combining these data with existing elastic and thermal data on similar alloys, we find that those alloys which can undergo diffusionsless phase transitions, such as Ti:Nb, Ti:V, or Zr:Nb in certain ranges of composition have glasslike excitations, since they have elastic properties which agree in magnitude and temperature dependence with those of amorphous solids. By contrast, crystalline continuous solution alloys, such as Nb:Ta, or alloys with diffusive phase transitions, such as high-pressure quenched Al94Si6, have the same elastic properties as are known for crystals.

The gasification of wood-char spheres in CO2---N2 mixtures: analysis and experiments

The gasification of charcoal spheres in an atmosphere of carbon-dioxide-nitrogen mixture involving diffusion and reactions in the pores is modelled and the results are compared with experiments of Standish and Tanjung and those performed in the laboratory on wood-char spheres to determine the effects of diameter, density, gas composition and flow. The results indicate that the conversion time, t(c) approximately d1.03 for large particles (> 5 mm), departing substantially from the t(c) approximately d2 law valid for diffusion limited conditions. The computational studies indicate that the kinetic limit for the particle is below 100 mum. The conversion time varies inversely as the initial char density as expected in the model. Predictions from the model show that there is no significant change in conversion time up to 60% N2 consistent with the CO2-N2 experiments. The variation of diameter and density with time are predicted. The peculiar dependence of conversion time on flow velocity in the experiments is sought to be explained by opposing free and forced convection heat transfer and the attempt is only partly successful. The studies also indicate that the dependence on the CO concentration with low CO2 is significant, indicating the need for multistep reaction mechanism against the generally accepted single-step reaction.

Use of L matrix to obtain reliable ab initio force constants of polyatomic molecules: ethylene as a test system

Direct use of experimental eigenvalues of the vibrational secular equation on to the ab initio predicted eigenvector space is suggested as a means of obtaining a reliable set of intramolecular force constants. This method which we have termed RECOVES (recovery in the eigenvector space) is computationally simple and free from arbitrariness. The RECOVES force constants, by definition, reproduce the experimental vibrational frequencies of the parent molecule exactly. The ab initio calculations were carried out for ethylene as a test molecule and the force constants obtained by the present procedure also correctly predict the vibrational frequencies of the deuterated species. The RECOVES force constants for ethylene are compared with those obtained by using the SQM procedure.

Calculation of elastic pinning strength of semicoherent Y2BaCuO5 precipitates in YBa2Cu3O7

The pinning energy due to the elastic interaction of a semicoherent Y2BaCuO5 precipitate with the YBa2Cu3O7 matrix is computed. This is achieved by setting up dislocation arrays at the interface. The elastic stresses generated by such arrays are integrated over a fluxoid volume to obtain the energy. It is seen that this elastic interaction energy makes an additive contribution to the total J(c) value.

Nonlinear static and dynamic damped response of a composite rectangular plate resting on a two-parameter elastic foundation

Nonlinear static and dynamic response analyses of a clamped. rectangular composite plate resting on a two-parameter elastic foundation have been studied using von Karman's relations. Incorporating the material damping, the governing coupled, nonlinear partial differential equations are obtained for the plate under step pressure pulse load excitation. These equations have been solved by a one-term solution and by applying Galerkin's technique to the deflection equation. This yields an ordinary nonlinear differential equation in time. The nonlinear static solution is obtained by neglecting the time-dependent variables. Thc nonlinear dynamic damped response is obtained by applying the ultraspherical polynomial approximation (UPA) technique. The influences of foundation modulus, shear modulus, orthotropy, etc. upon the nonlinear static and dynamic responses have been presented.

Interaction analysis of strip footing resting on a non-homogeneous elastic medium

In this paper, a plane stress solution for the interaction analysis of strip footing resting on (i) a non-homogeneous elastic half-plane and (ii) a non-homogeneous elastic layer resting on a rigid stratum has been presented. The analysis has been done using a combined analytical and FEM method in which the discretization of the half-plane is not required and thereby minimizes the computational efforts considerably. The contact pressure distribution and the settlement profile for the selected cases of varying modulus half-plane, which has more relevance to foundation engineering, have been given. Experimental verification through a photoelastic method of stress analysis has been carried out for the case of footing on Gibson elastic half-plane, and the contact pressure distribution thus obtained has been compared with the theoretical results. Copyright (C) 1996 Elsevier Science Ltd