Anisotropic effective moduli of microcracked materials under antiplane loading

This study focuses on the prediction of the anisotropic effective elastic moduli of a solid containing microcracks with an arbitrary degree of alignment by using the generalized self-consistent method (GSCM). The effective elastic moduli pertaining to anti-plane shear deformation are discussed in detail. The undamaged solid can be isotropic as well as anisotropic. When the undamaged solid is isotropic, the GSCM can be realized exactly. When the undamaged solid is anisotropic it is difficult to provide an analytical solution for the crack opening displacement to be used in the GSCM, thus an approximation of the GSCM is pursued in this case. The explicit expressions of coupled nonlinear equations for the unknown effective moduli are obtained. The coupled nonlinear equations are easily solved through iteration.     

Basic relations of the theory of nonlinear dynamic deformation of high-strength polymer films

We present basic relations of the theory of nonlinear dynamic deformation of rigid (high-strength) polymeric films. On the basis of the experimental results, we formulate the determining relations for the elastic and plastic stages of deformation. We introduce a measure of the stress-strain state and obtain the equations of motion for arbitrary curvilinear orthogonal coordinates in the case of large displacements and low strains and squares of the angles of rotation. Timoshenko Institute of Mechanics, National Academy of Sciences of Ukraine, Kiev, Ukraine. Translated from Problemy Prochnosti, No. 6, pp. 127–133, November–December, 1999.     

Mathematical foundations of the immersed finite element method

In this paper, we propose an immersed solid system (ISS) method to efficiently treat the fluid–structure interaction (FSI) problems. Augmenting a fluid in the moving solid domain, we introduce a volumetric force to obtain the correct dynamics for both the fluid and the structure. We further define an Euler–Lagrange mapping to describe the motion of the immersed solid. A weak formulation (WF) is then constructed and shown to be equivalent to both the FSI and the ISS under certain regularity assumptions. The weak formulation (WF) may be computed numerically by an implicit algorithm with the finite element method, and the streamline upwind/Petrov–Galerkin method. Compared with the successful immersed boundary method (IBM) by Peskin and co-workers (J Comput Phys 160:705–719, 2000; Acta Numerica 11:479–517, 2002; SIAM J Sci Stat Comput 13(6):1361–1376, 1992) the ISS method applies to more general geometries with non-uniform grids and avoids the inaccuracy in approximating the Dirac delta function     

A combined dislocation—cohesive zone model for fracture in a confined ductile layer

The collective dislocation behavior near a crack tip in a ductile layer sandwiched between two brittle solids is analyzed via two-dimensional dislocation dynamics (DD) simulations that incorporate a cohesive zone (CZ) model. The cohesive crack tip is treated as part of a much larger finite crack confined in the ductile layer. The underlying boundary value problem is formulated with a set of boundary integral equations and numerically evaluated with a collocation method. The fracture energy of the layered composite material is shown to be strongly correlated with the layer thickness and is directly influenced by the cohesive strength of the ductile layer (Hsia KJ et al. (1994) J Mech Phys Solids 6 877–896).     

Composite materials whose phases have partially equal elastic moduli

Hill [J. Mech. Phys. Solids 11 (1963) 357, 12 (1964) 199] discovered that, regardless of its microstructure, a linearly elastic composite of two isotropic phases with identical shear moduli is isotropic and has the effective shear modulus equal to the phase ones. The present work generalizes this result to anisotropic phase composites by showing and exploiting the fact that uniform strain and stress fields exist in every composite whose phases have certain common elastic moduli. Precisely, a coordinate-free condition is given to characterize this specific class of elastic composites; an efficient algebraic method is elaborated to find the uniform strain and stress fields of such a composite and to obtain the structure of the effective elastic moduli in terms of the phase ones; sufficient microstructure-independent conditions are deduced for the orthogonal group symmetry of the effective elastic moduli. These results are applied to elastic composites consisting of isotropic, transversely isotropic and orthotropic phases.     

Elastic property prediction by finite element analysis with random distribution of materials for tungsten/silver composite

In the present numerical study, we introduce a finite element analysis for heterogeneous materials via a random distribution of materials to predict effective elastic properties. With this random distributing strategy, a large scale parametric analysis via finite element becomes feasible for the multi-phase heterogeneous solids. Taking a well-documented tungsten–silver bi-continuous material as an example, the numerical prediction provided here for the effective properties is checked by experimental testing data available in open publication. Discussions on the present finite element prediction and other approaches are also made by comparing with Hashin and Shtrikman (J Mech Phys Solids 11:127–140, 1963) bounds in the composite mechanics.     

Numerical study of the effect of nonlinear control on the behavior of a liquid drop in elongational flow with vorticity

The problem of controlling a liquid drop suspended in an arbitrary two-dimensional elongational flow with vorticity is revisited. Bentley and Leal (J Fluid Mech 137:219–240, 1986) kept the drop centroid at the stagnation point using a linear proportional control strategy in a four-roll-mill apparatus that projects the drop’s motion onto the stable flow direction of the stagnation point. A nonlinear strategy based on the Poincaré–Bendixson theory to ensure a periodic motion of the drop centroid inside a prescribed area around the stagnation point is proposed and studied. In addition, a detailed numerical study is presented to illustrate the effect of the control on the drop motion. The present strategy is effective, allowing for deformation and changes in the drop orientation by less than 1% for extreme flow conditions that cannot be achieved by a four-roll-mill setup.     

The effect of crack face contact on the anisotropic effective moduli of microcrack damaged media

The generalized self-consistent method (GSCM) in conjunction with a computational finite element method is used to calculate the anisotropic effective moduli of a medium containing damage consisting of microcracks with an arbitrary degree of alignment. Since cracks respond differently under different external loads, the moduli of the medium subjected to tension, compression and an initially stress-free state are evaluated and shown to be significantly different, which will further affect the wave speed inside the damaged media. There are four independent material moduli for a 2-D plane stress orthotropic medium in tension or compression, and seven independent material moduli for a 2-D plane stress orthotropic cracked medium, which is initially stress free. When friction exists, it further changes the effective moduli. Numerical methods are used to take into account crack face contact and friction. The wave slowness profiles for microcrack damaged media are plotted using the predicted effective material moduli.     

Analysis of in-plane transonically propagating interface crack with a finite contact zone

The report of Lambros and Rosakis [(1995) J Mech Phys Solids 43(2): 169–188] has focused attention on steady-state transonic interfacial crack growth accounting for the phenomenon of crack face contact in elastic/rigid bimaterial but could not handle issues relating to energy transmission across the interface. The present paper attempts to provide a complete explicit expression of the asymptotic fields induced by transonically propagating interfacial crack in elastic/elastic bimaterial for in-plane case. The energy distribution on the contact area, crack tip and two singular characteristic lines is analysed thoroughly and compared with the dynamic separated J-integrals. The length of the contact zone is also discussed briefly by establishing energy fracture criterion that satisfies contact condition. The two-dimensional in-plane asymptotic deformation field surrounding the contact area of a crack propagating transonically along an elastic/elastic bimaterial interface is observed and discussed thoroughly.     

An anisotropic elastic formulation for configurational forces in stress space

A new variational principle for an anisotropic elastic formulation in stress space is constructed, the Euler–Lagrange equations of which are the equations of compatibility (in terms of stress), the equilibrium equations and the traction boundary condition. Such a principle can be used to extend recently obtained configurational balance laws in stress space to the case of anisotropy.     

Numerical modeling of micro- and macro-behavior of viscoelastic porous materials

This paper presents a numerical method for solving the two-dimensional problem of a polygonal linear viscoelastic domain containing an arbitrary number of non-overlapping circular holes of arbitrary sizes. The solution of the problem is based on the use of the correspondence principle. The governing equation for the problem in the Laplace domain is a complex hypersingular boundary integral equation written in terms of the unknown transformed displacements on the boundaries of the holes and the exterior boundaries of the finite body. No specific physical model is involved in the governing equation, which means that the method is capable of handling a variety of viscoelastic models. A truncated complex Fourier series with coefficients dependent on the transform parameter is used to approximate the unknown transformed displacements on the boundaries of the holes. A truncated complex series of Chebyshev polynomials with coefficients dependent on the transform parameter is used to approximate the unknown transformed displacements on the straight boundaries of the finite body. A system of linear algebraic equations is formed using the overspecification method. The viscoelastic stresses and displacements are calculated through the viscoelastic analogs of the Kolosov–Muskhelishvili potentials, and an analytical inverse Laplace transform is used to provide the time domain solution. Using the concept of representative volume, the effective viscoelastic properties of an equivalent homogeneous material are then found directly from the corresponding constitutive equations for the average field values. Several examples are given to demonstrate the accuracy of the method. The results for the stresses and displacements are compared with the numerical solutions obtained by commercial finite element software (ANSYS). The results for the effective properties are compared with those obtained with the self-consistent and Mori–Tanaka schemes.     

Molecular dynamics simulations of HMX crystal polymorphs using a flexible molecule force field

Molecular dynamics simulations using a recently developed quantum chemistry-based atomistic force field [J. Phys. Chem. B, 103 (1999) 3570 ] were performed in order to obtain unit cell parameters, coefficients of thermal expansion, and heats of sublimation for the three pure crystal polymorphs of octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX). The predictions for β-, α-, and δ-HMX showed good agreement with the available experimental data. For the case of β-HMX, anisotropic sound speeds were calculated from the molecular dynamics simulation-predicted elastic coefficients and compared with recent Impulsive Stimulated Light Scattering (ISLS) sound speed measurements. The level of agreement is encouraging. This revised version was published online in June 2006 with corrections to the Cover Date.     

An alternating method based on the VNA solution for analysis of damaged solid containing arbitrarily distributed elliptical microcracks

This paper presents the development of an alternating method for the interaction analysis of arbitrary distributed numerous elliptical microcracks. The complete analytical solutions (VNA solutions) for a single elliptical crack in an infinite solid, subject to arbitrary crack-face tractions, are implemented in the present alternating method, together with the coordinate transformations for stress tensors. First, the present method is verified by solving the problems of two interacting cracks for which accurate numerical solutions have been obtained previously. Next, the present method demonstrates obtaining efficient and accurate solutions for the problems of many interacting elliptical cracks, which cannot be solved in a practical sense by the ordinary numerical methods such as the finite element method. Furthermore, damaged solids containing periodically distributed elliptical microcracks are analyzed by the present alternating method. The effective elastic moduli are evaluated for varying microcrack density. Detailed structures of the interactions in the damaged solids are visualized and clarified. This revised version was published online in July 2006 with corrections to the Cover Date.     

An Altgernating Method Based on the VNA Solution for Analysis of Damaged Solid Containing Arbitrarily Distributed Elliptical Microcracks

This paper presents the development of an alternating method for the interaction analysis of arbitrary distributed numerous elliptical microcracks. The complete analytical solutions (VNA solutions) for a single elliptical crack in an infinite solid, subject to arbitrary crack-face tractions, are implemented in the present alternating method, together with the coordinate transformations for stress tensors. First, the present method is verified by solving the problems of two interacting cracks for which accurate numerical solutions have been obtained previously. Next, the present method demonstrates obtaining efficient and accurate solutions for the problems of many interacting elliptical cracks, which cannot be solved in a practical sense by the ordinary numerical methods such as the finite element method. Furthermore, damaged solids containing periodically distributed elliptical microcracks are analyzed by the present alternating method. The effective elastic moduli are evaluated for varying microcrack density. Detailed structures of the interactions in the damaged solids are visualized and clarified. This revised version was published online in August 2006 with corrections to the Cover Date.     

Solitary Waves of Maxwell's Equations in Nonlinear Anisotropic Media

Abstract

The (1 + 1)-D solitary wave solutions of Maxwell's equations in nonlinearity induced anisotropic media (in liquids such as carbon disulphide, and in crystals, etc.) are investigated. We find that there is no arbitrarily linearly polarized (in the x-y plane perpendicular to the propagation direction z) soliton solution from Maxwell's equations except that with linear polarization either in alignment with or orthogonal to the geometric axis of the light induced refractive index change. This contradicts the prediction of the vector nonlinear Schroedinger equation (an approximation of Maxwell's equations) which yields soliton solutions with an arbitrary linear polarization. However, Maxwell's equations are found to admit stable elliptically polarized solitary wave solutions which reduce to the stable circularly polarized solitary wave solutions of the vector nonlinear Schroedinger equation when the induced refractive index change approaches zero.     

Efficiency and accuracy of fluid-structure interaction simulations using an implicit partitioned approach

An implicit partitioned arbitrary Lagrangian– Eulerian approach for fluid-structure interaction computations is considered. Enhancements of the coupled solution procedure by nonlinear multigrid techniques, an adaptive underrelaxation, and proper grid movement techniques are investigated.     

Investigation of macrocrack-microcrack interaction problems in anisotropic elastic solids – Part I: General solution to the problem and application of the J-integral

A general solution is derived for the plane problem of multiple microcracks near the tip in process zone of a semi-infinite macrocrack in an anisotropic elastic solid. The pseudo-traction method, addressed thoroughly in isotropic cases, is extended to anisotropic cases. A system of Fredholm integral equations, with difficulty in evaluation of the singular integrals, is solved by invoking the asymptotic behavior of the pseudo-tractions on the macrocrack faces. The interaction effect of the release of residual stresses due to near-tip microcracking is then evaluated. The J-integral analysis is also performed to give a consistency check of the solution and some useful conclusions. This revised version was published online in July 2006 with corrections to the Cover Date.     

Simulation of nonisothermal moisture transfer and stresses in wood in drying

We obtain a numerical solution for a two-dimensional problem of nonisothermal moisture transfer in the anisotropic structure of wood (lumber) in convective drying. On the basis of elasticity theory we determine internal thermal and moisture stresses. Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 71, No. 3, pp. 404–411, May–June, 1998.