A comparative study of fiber-matrix interactions by using micromechanical models

The purpose of this study is to describe the interfacial interactions in terms of stress distributions on short fibers in fiber-matrix unit-cell models. The fiber and matrix are subjected to tensile loading. The study consists of three main parts. First, fiber-matrix cell segments are modeled using a 3D finite-element analysis (FEA) with ANSYS. Three different finite-element geometrical unit-cell models are generated in order to simulate the Cox analytical model: a fiber-matrix combination, a single fiber, and a single matrix element. The second part contains the results of 3D FE analyses, which are applied to the Cox formulations by using a computer program developed. In the last part, the analytical solutions for distributions of normal and shear stresses are investigated. Cox 2D linear elasticity solutions, together with finite-element ones, are presented in detail in graphs. The interfacial interactions between the fibers and matrix are also discussed considering the relative changes in the distributions of normal and shear stresses. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 44, No. 4, pp. 505–520, July–August, 2008.     

The influence of imperfections on the creep behavior of woven polymer composites at elevated temperatures

The dynamic behavior of compression molded polymer/woven graphite fiber composites at elevated temperatures is investigated analytically. This is performed with the objective of predicting the initiation of catastrophic failure that may occur after prolonged usage of the material at these temperatures. Special attention is paid to the behavior of the voids present in them where the failure may occur. The polymer matrix is modeled as a 4-parameter model (Maxwell-Voigt combination) (Govindarajan et al., in: Advances in Computer-Aided Engineering, ASME, 1994) while the composite structure is modeled using the fiber undulation model (Ishikawa and Chou, J. Mater. Sci. 17, pp. 3211–3220, 1982). The relation between the polymer properties and the ambient temperature is modeled after Arhenius' relation (Govindarajan et al., 1994; Ferry, Viscoelastic Properties of Polymers, Wiley, New York, 1961). The multiple phases in the matrix are taken into account through Eshelby's theory (Proc. Royal Soc. London A 241, pp. 376–396, 1957) and its extension for multiple occurrences of the same phase (Tanden and Weng, Polymer Composites 5, pp. 327–333, 1984; Weng, Internat. J. Eng. Sci. 22 (7), pp. 845–856, 1984) which assumes an ellipsoidal shape for inclusions. The resulting elastic equations are transformed into the time domain using Laplace transformation and the correspondence principle (Govindarajan et al., 1994; Wang and Weng, ASME J. Appl. Mech, 1992). All the voids are considered to be prolate ellipsoids with the 1-axis being the axis of symmetry. The distribution of voids is assumed to be of a Gaussian form with respect to the aspect ratio. The response of the composite under creep condition (constant load) has been simulated. Relations between the applied stress and the stresses in the matrix/void phase are also supplied, so that the influence of the voids may be characterized. The model is then applied to simulate the behavior of an epoxy/woven graphite composite to obtain the numerical results.     

Investigation of the mechanical behavior of two-component granular composites in terms of structural models

A structural model is suggested for elastomers filled with particles of two fractions — with diameters exceeding 10 μm and submicronic ones. In each fraction, the particle diameter varies randomly, but between the fractions, the average particle diameter differs by several orders of magnitude. It is assumed that the small particles, together with the matrix, behave as a homogeneous medium relative to the large ones. By using this model, the mechanical behavior of composites based on elastomers filled with different volume contents of solid particles is investigated. __________ Translated from Mekhanika Kompozitnykh Materialov, Vol. 43, No. 2, pp. 191–200, March–April, 2007.     

Application of the semy-analytical micromechanical methods for optimization of the elastic response of metal-ceramic composites

Metal-ceramic composites produced by melt infiltration of the ceramic preform were studied. The ceramic preform of these materials is manufactured through freeze-casting process and the microstructure can be presented as distribution of the lamellar domains. Micromechanical models were used for the calculation of the effective elastic properties of the domains. The minimum compliance problem was solved for sample subjected to the four-point bending test. The optimal orientation and volume fracture of the micro constituents were identified using different models. The difference between the initial and optimized design was analyzed. The convergence of the provided optimization procedure was very rapid. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical method of mixed finite volume-modified upwind fractional step difference for three-dimensional semiconductor device transient behavior problems

Transient behavior of three-dimensional semiconductor device with heat conduction is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditions. The electric potential is defined by an elliptic equation and it appears in the following three equations via the electric field intensity. The electron concentration and the hole concentration are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite volume element approximation, keeping physical conservation law, is used to get numerical values of the electric potential and the accuracy is improved one order. Two concentrations and the heat conduction are computed by a fractional step method combined with second-order upwind differences. This method can overcome numerical oscillation, dispersion and decreases computational complexity. Then a three-dimensional problem is solved by computing three successive one-dimensional problems where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L² norm is derived by using prior estimate theory and other special techniques of partial differential equations. This type of mass-conservative parallel method is important and is most valuable in numerical analysis and application of semiconductor device.     

Characterization of the classical domains by Caratheodory volume elements and Eisenman-Kobayashi volume elements

We give a sufficient and necessary condition that a domain is biholomorphic to the classical domain.     

On a Representative Volume in the Micromechanics of Particulate Composites

A rectangle filled with closely packed spheres of random size and properties is considered as a micromechanical model of a two-phase particulate composite. A numerical simulation is used to determine the effective mechanical properties of the assembly and their scatter as a function of the number of spheres. It is shown that, in a system with relatively small number of particles (up to 300), the scatter of Young's modulus decreases with the system size. However, the rate of the scatter decrease becomes smaller with growing size of the system, so that the convergence to zero most likely takes place at infinity.     

有界对称域的特征

We give a necessary and sufficient condition for a domain to be biholomorphic to a bounded symmetric domain.     

Investigation of the dynamic characteristics of composites

The dynamic characteristics E and G and the damping capacity δ* of polymeric resins and composites based on glass, carbon, and boron fibers have been investigated. It is shown that the mechanical losses are correlated with the modulus of elasticity and the shear modulus of the resin and composites with different types of reinforcement. The vibrational strength of various structural materials is estimated.     

Analysis of elastomeric composites based on fiber-reinforced systems 3. Two-directional composites

Results of investigation of deformation of elastomeric composite materials with a two-directional reinforcement scheme are presented. The study is performed on the basis of a structural macroscopic theory. The matrix of the composites analyzed is of a poorly compressible material. The fibers of both reinforcing systems are simulated as compressible bodies. Dependences of the parameters of tensile and shear strains on the strain values for different geometries of fiber arrangement are obtained.State Metallurgical Academy of Ukraine, Dnepropetrovsk, Ukraine. Translated from Mekhanika Kompozitnykh Materialov, Vol. 35, No. 4, pp. 479–492, May–June, 1999.     

Constitutive model of ferroelectric composites with a viscoelastic and dielectric relaxation matrix I

Based on micromechanics and Laplace transformation, a constitutive model of ferroelectric composites with a linear elastic and linear dielectric matrix is developed and extended to the ferroelectric composites with a viscoelastic and dielectric relaxation matrix. Thus, a constitutive model for ferroelectric composites with a viscoelastic and dielectric relaxation matrix has been set up Project supported by the National Natural Science Foundation of China (Grant No. 19891180).     

3D modeling of crack growth in electro-magneto-thermo-elastic coupled viscoplastic multiphase composites

This work presents a time-domain hypersingular integral equation (TD-HIE) method for modeling 3D crack growth in electro-magneto-thermo-elastic coupled viscoplastic multiphase composites (EMTE-CVP-MCs) under extended incremental loads rate through intricate theoretical analysis and numerical simulations. Using Green’s functions, the extended general incremental displacement rate solutions are obtained by time-domain boundary element method. Three-dimensional arbitrary crack growth problem in EMTE-CVP-MCs is reduced to solving a set of TD-HIEs coupled with boundary integral equations, in which the unknown functions are the extended incremental displacement discontinuities gradient. Then, the behavior of the extended incremental displacement discontinuities gradient around the crack front terminating at the interface is analyzed by the time-domain main-part analysis method of TD-HIE. Also, analytical solutions of the extended singular incremental stresses gradient and extended incremental integral near the crack fronts in EMTE-CVP-MCs are provided. In addition, a numerical method of the TD-HIE for a 3D crack subjected to extended incremental loads rate is put forward with the extended incremental displacement discontinuities gradient approximated by the product of time-domain basic density functions and polynomials. Finally, examples are presented to demonstrate the application of the proposed method.     

Measuring and modeling the thermal conductivities of three-dimensionally woven fabric composites

The effect of a three-dimensional fiber reinforcement on the out-of-plane thermal conductivity of composite materials is investigated. Composite preforms with different fibers in the thickness direction were fabricated. After in fusion by using a vacuum-assisted resin transfer molding process, their through-thickness thermal conductivities were evaluated. The measured thermal conductivities showed a significant increase compared with those of a typical laminated composite. Although the through-thickness thermal conductivity of the samples increased with through-thickness fiber volume fraction, its values did not match those predicted by the simple rule of mixtures. By using finite-element models to better under stand the behavior of the composite material, improvements in an existing analytical model were performed to predict the effective thermal conductivity as a function of material properties and in-contact thermal properties of the composite. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 45, No. 2, pp. 241–254, March–April, 2009.     

Bounds and estimates on the effective properties for nonlinear composites

In this paper we derive lower bounds and upper bounds on the effective properties for nonlinear heterogeneous systems. The key result to obtain these bounds is to derive a variational principle, which generalizes the variational principle by P. Ponte Castaneda from 1992. In general, when the Ponte Castaneda variational principle is used one only gets either a lower or an upper bound depending on the growth conditions. In this paper we overcome this problem by using our new variational principle together with the bounds presented by Lukkassen, Persson and Wall in 1995. Moreover, we also present some examples where the bounds are so tight that they may be used as a good estimate of the effective behavior.     

Comparative micromechanical analysis of nonelastic behavior of unidirectional plastics with tetragonal and hexagonal structures

A mathematical model for determining the effective elastic properties and describing the processes of inelastic deformation and damage accumulation of unidirectional fiber-reinforced composites with tetragonal and hexagonal structures is developed. A comparative analysis of the effective elastic moduli of glass, boron, organic, and carbon unidirectional plastics shows that, if the fiber volume fraction does not exceed 0.5, the effective elastic properties calculated by the models presented give closely related results. The calculation results for nonlinear fields of deformation and failure are presented and the limiting strength surfaces of fibrous glass plastics with hexagonal and tetragonal structures are obtained for different transverse loading paths. It is found that the structure of a composite affects significantly its strength properties.Submitted to the 11th International Conference on Mechanics of Composite Materials (Riga, June 11–15, 2000).Perm' State Technical University, Perm', Russia. Translated from Mekhanika Kompozitnykh Materialov, Vol. 36, No. 4, pp. 455–464, July–August, 2000.     

Formation of a transition layer on the fillers of polymer composites

Based on a plane model of composites, the effect of a transition layer on the elastic modulus E_c of the composites is analyzed in the case where, under the action of a load, the transition layer is formed both on the side of matrix and filler. In evaluating E_c, it is assumed that the elastic modulus in the layer grows linearly from the elastic modulus of matrix to that of filler, but pores in the filler are impermeable to matrix macromolecules. Analytic relation ships are found which allow one to determine the volume fractions of the transition layer on the side of matrix and filler if the experimental elastic modulus of the composite is known. These relationships are used to find the magnitude of the layer in epoxy composites with various fillers and to evaluate its effect on the compressive elastic modulus of the composites. __________ Translated from Mekhanika Kompozitnykh Materialov, Vol. 42, No. 5, pp. 693–700, September–October, 2006.