Proceedings Article | 29 March 2019
KEYWORDS: Composites, Interfaces, Polymers, Solid modeling, Atmospheric modeling, Solids, Particles, Atmospheric particles, Data modeling, Mechanics
Polymer composites of particulate smart materials are increasingly relevant due to emerging trends in using additive manufacturing of replacement parts in automotive and aerospace applications. The mechanical properties of the polymer matrix is sufficiently well understood from decades of research and the mechanical properties of the interface between the particulate (smart material) and matrix (polymer) phase is not very well defined. The challenges associated with understanding the interface requires a multidisciplinary approach and expertise and varies dramatically between various material combinations. In this article, we establish a robust approach to experimentally determine the linear interface modulus using cohesive zone model and demonstrate its application for mechanoluminescent elastomeric composite.
The phenomenon of light emission induced by any mechanical action is termed mechanoluminescence and is classified based on the applied load regime as fracto, plastico, and elastico-mechanoluminescence (EML). EML is repeatable as the applied load is within elastic limits and induced strains are recoverable. Repetitiveness paves way for utilization in several applications including, but not limited to, stress sensing, stress visualization, fatigue monitoring, damage detection and failure prevention. Previously published work on EML materials have mostly been exploratory in nature with motivation to fabricate brighter EML materials and characterize/model their mechanism of emission. Through such work, doped zinc sulfide phosphors and doped strontium aluminate phosphors have been reported to emit the brightest emission till date. Characterization studies generally involve testing composites with EML crystals impregnated in a polymeric matrix, or EML thin films grown on various substrates. With both approaches, the focus is generally placed on correlating light output to various applied macroscopic inputs. However, there are no known efforts focusing on interfacial adhesion of EML crystals with the surrounding matrix or the substrate. Preliminary experimental and numerical research work on elastico-mechanoluminescent (EML) composites revealed that the interfacial adhesion between particle and matrix played a significant role in stress transfer and EML emission. Contemporary models for particle-reinforced composites focus on estimating macroscopic composite properties [1, 2 and 3]. No known model focusses on stresses imparted to the functional filler particles at the interface. Since EML emission depends on stress on the filler EML particles and not the macroscopic applied stress, it becomes necessary to develop models for stress transfer between the matrix and the particle at the interface. Estimation of stress transfer at the interface requires knowledge of stress distribution in the elastomer matrix. Hence, the objective of the models to be developed for elastomer composites is to capture variations in macroscopic stress of the matrix due to softening of matrix and degradation of interface with fatigue.
EML-Elastomer Composites
Elastomers are highly non-linear elastic materials which undergo significant softening during first few actuation cycles (Mullins effect) [4]. Mechanical behavior is rate-dependent exhibiting hysteresis during cyclic loading [5], and fatigue analysis of this complex system requires knowledge about loading history, composite composition and filler properties [6]. Debonding of filler-matrix interface is also an important non-linearity that cannot be neglected. While focusing attention on damage and fatigue of filled elastomers, micromechanics based and finite element based approaches have been adopted in combination by researchers to estimate macroscopic mechanical properties of the composite [7, 8]. The micromechanics approach to constitutive modeling of particulate composites generally builds upon the pioneering works of J.D. Eshelby in 1957 and T. Mori and K. Tanaka in 1973 [9, 10]. J.D. Eshelby derived analytical solutions to problems involving ellipsoidal elastic inclusions in infinite elastic bodies which are in general referred to as Eshelby’s solutions. Mori and Tanaka modeled composites with N number of filler particles as ellipsoidal inclusions in perfect adherence with matrix. Extension of the Mori-Tanaka model to elastomeric composites requires consideration of nonlinear hyperelastic and viscoelastic behavior of matrix as well as particle debonding. The non-linearity from hyperelastic behavior of matrix is generally accounted for by finite strain deformations and strain energy based formulations (for example, Ogden’s model), while particle debonding is modeled through non-linear cohesive laws arrived at from experimental data. Cohesive law is the relationship between traction forces and displacements at the interface [11]. Viscoelastic damage and hysteresis have to be accounted for separately, by modifying appropriate energy functions to include viscous dissipation determined from experimental data [12, 13, and 14].
In this article, we establish a robust approach to experimentally determine the linear interface modulus using cohesive zone model and demonstrate its application for mechanoluminescent elastomeric composite. The experimental design uses photoluminescence property of phosphors (Osram Sylvania GG45 - ZnS:Cu) to determine damage propagation and digital image correlation (DIC) to determine stress and strain in the polymer composite at the micron scale.A sample measuring 15mm x 5mm x 1.5mm is prepared using PDMS and GG45 and a notch is introduced to the sample. A microscope is modified to pull the sample apart at a constant speed and the process is observed using a microscope objective and a 4K video camera. The photoluminescence of dispersed GG45 particles is adjusted by mixing black ink to the composite to control halo effect and subsequent processing in image processing. Stress and strain computed from NICORR running in MATLAB is coupled with the opening displacement to obtain a measure of stiffness in Stage-I, II and III stiffness in a sample of varying extents of initial damage. The stiffness measured from linear softening CZM response is used to obtain the interface stiffness that can be subsequently applied in finite element models for modeling the overall property of composite. We demonstrate that an undamaged sample has the highest stiffness. With the formation of micron scale damage in the polymer composite at the interface, the interface stiffness decreases and reduces the maximum load that can be applied to the material in subsequent cycles.
References:
[1] Allan Zhong, Wolfgang G. Knauss, X. (2000). Effects of particle interaction and size variation on damage evolution in filled elastomers. Mechanics of composite materials and structures, 7(1), 35-53.
[2] Yeh, J. R. (1992). The effect of interface on the transverse properties of composites. International journal of solids and structures, 29(20), 2493-2502.
[3] Nakamura, Y., Yamaguchi, M., Okubo, M., & Matsumoto, T. (1992). Effect of particle size on the fracture toughness of epoxy resin filled with spherical silica. Polymer, 33(16), 3415-3426.
[4] Mullins, L., and N. R. Tobin. (1957). "Theoretical model for the elastic behavior of filler-reinforced vulcanized rubbers." Rubber Chemistry and Technology 30.2: 555-571.
[5] Ravichandran, G., and C. T. Liu. (1995). "Modeling constitutive behavior of particulate composites undergoing damage." International Journal of Solids and Structures 32.6-7: 979-990.
[6] Dorfmann, A., and Ray W. Ogden. (2004). "A constitutive model for the Mullins effect with permanent set in particle-reinforced rubber." International Journal of Solids and Structures 41.7:1855-1878.Eshelby, John D. (1957). "The determination of the elastic field of an ellipsoidal inclusion, and related problems." Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences. Vol. 241. No. 1226. The Royal Society.
[7] Tan, H., Huang, Y., Liu, C., Ravichandran, G., Inglis, H. M., & Geubelle, P. H. (2007). The uniaxial tension of particulate composite materials with nonlinear interface debonding. International Journal of Solids and Structures, 44(6), 1809-1822.
[8] Toulemonde, P. A., Diani, J., Gilormini, P., & Desgardin, N. (2016). On the account of a cohesive interface for modeling the behavior until break of highly filled elastomers. Mechanics of Materials, 93, 124-133.
[9] Eshelby, John D. (1957). "The determination of the elastic field of an ellipsoidal inclusion, and related problems." Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences. Vol. 241. No. 1226. The Royal Society.
[10] Mori, Tanaka, and K. Tanaka. (1973). "Average stress in matrix and average elastic energy of materials with misfitting inclusions." Acta Metallurgica 21.5:571-574.
[11] Needleman, Alan. (1987). "A continuum model for void nucleation by inclusion debonding." Journal of applied mechanics 54.3:525-531.
[12] Ravichandran, G., & Liu, C. T. (1995). Modeling constitutive behavior of particulate composites undergoing damage. International Journal of Solids and Structures, 32(6-7), 979-990.
[13] Yang, B. J., Kim, B. R., & Lee, H. K. (2012). Micromechanics-based viscoelastic damage model for particle-reinforced polymeric composites. Acta Mechanica, 223(6), 1307.
[14] Jung, G. D., Youn, S. K., & Kim, B. K. (2000). A three-dimensional nonlinear viscoelastic constitutive model of solid propellant. International Journal of Solids and Structures, 37(34), 4715-4732.
