Difference between revisions of "Imperfect Interfaces"
| Line 34: | Line 34: | ||
== Imperfect Interface Options in MPM == | == Imperfect Interface Options in MPM == | ||
== Imperfect | == Imperfect Interfaces in FEA == | ||
To use imperfect interfaces in [[NairnFEA]], use the following procedure: | To use imperfect interfaces in [[NairnFEA]], use the following procedure: | ||
Revision as of 23:15, 14 September 2013
Both NairnMPM and NairnFEA can model imperfect interfaces between phases in composite structures.
Imperfect Interface Theory
Modeling interfaces, which are often finite-thickness interphases, in composite materials is difficult. In numerical modeling, it would seemingly be straight-forward to discretize interphase zones and thereby explicitly model all effects. This approach has two problems. First, interphase zones may be much smaller than the bulk materials. Resolving both bulk materials and a thin interphase would require a highly refined model, which may exceed computational capacity. Second, interphase properties may be unknown and/or may vary within a transition zone from one material to another. Alternative methods for interphase modeling are needed to overcome these challenges. This need is especially important in nanocomposites because the amount of interphase per unit volume of reinforcement greatly exceeds the amount of interphase in composites with mi cron or larger reinforcement phases. As a consequence, interphases are expected to play a larger role, good or bad, in nanocomposite properties.
One way to model interphases is to abandon attempts for explicit modeling and instead replace 3D interphases with 2D interfaces[1] The interphase effects are reduced to modeling the response of 2D interfaces due to tractions normal and tangential to the interfacial surface, which can be modeled by interface traction laws. Elimination of 3D interphases removes the resolution problem. The use of interface traction laws always replaces numerous unknown and potentially unmeasurable interphase properties with a much smaller number of interface parameters. If interface traction laws can be determined, one can potentially model interphase effects well. Example of using this for analytical modeling of interface effects in composite materials are in Hashin (1991)[2][3] and Nairn and Liu (1997)[4] Implementation of imperfect interfaces in FEA and in MPM is described in Nairn (2007).[5] An alternative method for modeling imperfect interfaces in MPM using multimaterial MPM is described in Nairn (2013).[6]
Imperfect Interface Traction Laws
To model 3D interphase effects using imperfect interfaces, the 2D interface is allowed to develop displacement discontinuities. For an isotropic interphase, it suffices to resolve interfacial traction into normal and tangential tractions (Tn and Tt) and assume they are functions of normal and tangential displacement discontinuities ([un] and [un]) at the interface:
[math]\displaystyle{ T_n = f_n([u_n]) \qquad {\rm and } \qquad T_t = f_t([u_t]) }[/math]
Interfaces in composite materials may develop potential energy that is needed for effective property analysis.[7] For an elastic interface, interfacial potential energy is:
[math]\displaystyle{ \phi_i = \int_{S_i}\left(\int \vec T\cdot[d\vec u]\right) dS }[/math]
where Si is interfacial area. The simplest assumption for traction laws is that they are linear and elastic. The tractions become:
[math]\displaystyle{ T_n = D_n[u_n] \qquad {\rm and } \qquad T_t = D_t[u_t] }[/math]
and interfacial energy becomes
[math]\displaystyle{ \phi_i = {1\over2} \int_{S_i} \left( D_n[u_n]^2 + D_t[u_t]^2\right) dS }[/math]
Notice that all 3D interphase properties have been reduced to just two interface parameters, Dn and Dt. These interface stiffnesses range from zero, for a debonded interface with no tractions, to infinity, for a perfect interface with no displacement discontinuity.
Imperfect Interface Options in MPM
Imperfect Interfaces in FEA
To use imperfect interfaces in NairnFEA, use the following procedure:
- Define one or more imperfect interface materials and enter their interface parameters. Because NairnFEA is for linear elastic analysis only, the imperfect interface materials have to be linear elastic as well. In other words, the traction laws are linear with interface parameters Dn and Dt.
- Define two different, but overlapping paths that have the some number and arrangement of nodes. In other words, define the paths with the same number of nodes, having the same ratio of element sizes, and using keypoints that are distinct by have matching coordinates.
- Use the two paths in two areas for the two materials to be connected by an imperfect interface.
- Use the two paths in an interface element that is compatible with type of solid elements used in the areas.
Tractions in stresses
Energy in element energy
Figure of elements
References
- ↑ Z. Hashin, "Thermoelastic properties of fiber composites with imperfect interface," Mech. of Materials, 8, 333–348 (1990).
- ↑ Z. Hashin, "Composite materials with viscoelastic interphase: Creep and relaxation," Mech. of Materials, 11, 135–148 (1991).
- ↑ Z. Hashin, "Thermoelastic properties of particulate composites with imperfect interface," Journalof the Mechanics and Physics of Solids, 39, 745–762 (1991).
- ↑ J. A. Nairn and Y. C. Liu, "Stress transfer into a fragmented, anisotropic fiber through an imperfect interface," Int. J. Solids Structures, 34, 1255–1281 (1997).
- ↑ J. A. Nairn, "Numerical implementation of imperfect interfaces," Computational Materials Science, 40, 525–536 (2007).
- ↑ J.A. Nairn, "Modeling Imperfect Interfaces in the Material Point Method using Multimaterial Methods," Computer Modeling in Eng. & Sci., 92, 271-299 (2013).
- ↑ Z. Hashin, "Extremum principles for elastic heterogenous media with imperfect interfaces and their application to bounding of effective moduli." Journal of the Mechanics and Physics of Solids, 40, 767–781 (1992).