Inhomogeneities can be beneficial to mechanical properties, such as toughness, if distributed through a material effectively. Such inhomogeneities come in an array of formats, from grain boundaries and ceramic precipitations in engineered metal alloys to ceramic particulates in mammalian bone and crustaceans. In these cases, the organization of inhomogeneities is described by an optimized microscale architecture that is integral to macroscale toughness. Toughness represents the ability of materials systems to inhibit the initiation and propagation of cracks from defects and is often the determining factor of overall material performance, especially for more brittle systems under cyclic loading such as bone and structural composites. Specifically, toughness is the energy that a material can absorb prior to and during fracture and includes the initiation toughness, which inhibits the starting of a crack from a void or defect, the propagation toughness, which inhibits the growth of a nanoscopic crack, and the total toughness, which inhibits the complete failure or breaking of a sample or part. In all cases, the toughness is derived from both configurational system properties such as the loading and the sample or part geometry and inherent materials properties such as the fracture toughness.
Crack deflection is often used as evidence of substantive fracture toughness in both natural and synthetic composites. The associated microstructures of these tough composites are often highly anisotropic, exhibiting both easy and hard directions for crack propagation. Cracks generally deflect away from hard directions to propagate instead in easy directions. Allowing a crack to deflect toward easy directions ties the composite's overall fracture toughness to its weak axis. Instead, enhanced fracture toughness can be realized when cracks are forced to propagate along the hard directions. Composite manufacturing using 3D magnetic printing techniques can enable the creation of detailed microstructures that can locally tune the fracture toughness anisotropy. The technology described herein includes a family of architectures reminiscent of a herringbone pattern, in which propagating cracks are effectively trapped along the spine between two regions. Along this spine, a crack can be inhibited from deflecting and forced to propagate along a relatively hard direction, resulting in a significant boost in initial and overall fracture toughness. Toughened herringbone patterns can be expanded to herringbone mosaics and an approach is described for patterning composites to provide isotropic and enhanced fracture toughness. Both numerical and experimental modalities to systematically test and characterize this new class of composite architectures are described.
Other aspects, features, and embodiments include the following:
1. A composite material comprising:
a volume of a matrix material and a plurality of particles disposed within the matrix material, the particles having a geometric shape with a longest dimension and a shortest dimension different from the longest dimension; and
the volume has an anisotropic structure in at least two dimensions comprising:
a volume of a polymer matrix material and a plurality of reinforcing particles disposed within the matrix material, the particles having a platelet shape with a longest dimension and a shortest dimension different from the longest dimension, wherein the longest dimension is at least 200 nm and the particles have an aspect ratio of the longest dimension to the shortest dimension of at least 2; and
the volume has an anisotropic structure providing increased fracture toughness along a selected direction and an anisotropic fracture energy of at least 3, the anisotropic structure comprising:
(a) determining a pattern of an anisotropic structure of a composite material comprising a volume of a matrix material and a plurality of particles disposed within the matrix material, the particles having a geometric shape with a longest dimension and a shortest dimension different from the longest dimension, the anisotropic structure comprising:
(b) providing a precursor matrix material and the plurality of particles disposed within the precursor matrix material; and
(c) consolidating the precursor matrix material to form the composite material with the plurality of particles disposed within the matrix material according to the pattern.
40. The method of 39, further comprising determining the non-orthogonal angles between the particles within the adjacent volume portions and the interface between the adjacent volume portions as a function of fracture energy anisotropy of the composite material.
41. The method of any of 39-40, wherein the non-orthogonal angles are determined to provide a fracture energy anisotropy of at least 3.
42. The method of any of 39-41, wherein at least a subset of the plurality of the volume portions are arranged in the mosaic pattern to be adjacent to at least three other volume portions.
43. The method of any of 39-42, wherein the orientation of the particles in at least some of the adjacent volume portions are mirror symmetric with respect to the interface.
44. The method of any of 39-43, wherein the orientation of the particles in the adjacent volume portions are arranged to form a herringbone pattern in a plane normal to the interface, the interface forming a spine of the herringbone pattern.
45. The method of any of 39-44, wherein the non-orthogonal angle is selected to range from 20°±1° to 80°±1° with respect to the interface, the non-orthogonal angle is selected to range from 30°±1° to 75°±1° with respect to the interface, or the non-orthogonal angle is selected to range from 45°±1° to 70°±1° with respect to the interface.
46. The method of any of 39-45, wherein the volume comprises a layer having opposed surfaces and a thickness between the opposed surfaces, and the volume portions are arranged in a polygonal pattern within the layer.
47. The method of any of 39-46, wherein the polygonal pattern is a hexagonal pattern.
48. The method of any of 39-47, wherein the orientation of the particles in the adjacent volume portions of the hexagonal pattern is selected to differ by 30°±1° or by 60°±1°.
49. The method of any of 39-48, wherein the volume portions each have an area in a plane normal to the interface that at least ten times greater than the shortest dimension of the particles within the volume portions.
50. The method of any of 39-49, wherein the volume further comprises a plurality of layers, each layer having opposed surfaces and a thickness between the opposed surfaces, the layers arranged in a stack in a thickness direction, and the volume portions are arranged in a polygonal pattern within each layer.
51. The method of any of 39-50, wherein the polygonal pattern within each layer is a hexagonal pattern, and the hexagonal pattern within each layer is rotated in a plane parallel to the opposed surfaces with respect to the hexagonal pattern in an adjacent layer.
52. The method of any of 39-51, wherein the polygonal pattern within each layer is rotated with respect to the polygonal pattern in an adjacent layer to provide a Bouligand structure.
53. The method of any of 39-52, wherein each of the volume portions is a Wigner-Seitz unit cell.
54. The method of any of 39-53, wherein the anisotropic structure comprises a rotating herringbone structure, rows of stacked herringbone structures, or a Bouligand structure.
55. The method of any of 39-54, wherein the anisotropic structure of the volume is anisotropic in two dimensions or in three dimensions.
56. The method of any of 39-55, wherein the geometric shape of the particles is a platelet, flake, filing, fiber, rod, or whisker.
57. The method of any of 39-56, wherein the particles have a longest dimension ranging from 200 nm to 1000 μm.
58. The method of any of 39-57, wherein the particles have a longest dimension ranging from 1 μm to 20 μm.
59. The 58 method of any of 39-58, wherein the particles have an aspect ratio of the longest dimension to the shortest dimension of at least 2.
60. The method of any of 39-59, wherein the particles have an aspect ratio of the longest dimension to the shortest dimension ranging from 2 to 200.
61. The method of any of 39-60, wherein the particles are magnetically responsive particles comprising at least in part a magnetic material.
62. The method of any of 39-61, wherein the magnetic material comprises a ferromagnetic material, a paramagnetic material, a superparamagnetic material, iron oxide, iron, cobalt, nickel, an iron alloy, a cobalt alloy, or a nickel alloy.
63. The method of any of 39-62, wherein the magnetic material comprises particles, microbeads, nanoparticles, filings, fibers, flakes, rods, whiskers, or platelets.
64. The method of any of 39-63, wherein the magnetically responsive particles comprise a non-magnetic material coupled with the magnetic material.
65. The method of any of 39-64, wherein the non-magnetic material comprises aluminum oxide, calcium phosphate, copper, glass, calcium sulfate, nylon, polystyrene, or silicon carbide.
66. The method of any of 39-65, wherein the non-magnetic material comprises discontinuous fibers, rods, platelets, flakes, or whiskers.
67. The method of any of 39-66, wherein the non-magnetic material is coated with the magnetic material.
68. The method of any of 39-67, wherein the matrix material is a polymer and the particles comprise aluminum oxide platelets coated with a magnetically responsive material.
69. The method of any of 39-68, wherein the matrix material is a thermoset polymer, polyurethane, silicone, acrylic, acrylonitrile butadiene styrene, nylon, poly lactic acid, polycarbonate, polyethylene, polypropylene, polystyrene, polyvinyl chloride, polytetrafluoroethylene, polymethylmethacrylate, collagen, gelatin, or alginate.
70. A method of making the composite material of any of 1-38.
, the fracture energy anisotropy. The dashed line shows the optimal tilt angle for platelets αΓ to produce the toughest composite. “X” marks show parameters used for phase-field simulations shown in
=4 corresponding to experimental samples.
=4 (see also
=4, 6.
The technology described herein relates to composite materials of a matrix material reinforced with particles oriented in a variety of herringbone architectures to increase fracture toughness. In some embodiments, the composite material is a polymeric composite reinforced with alumina micro-platelets. The herringbone architectures can result in remarkable toughening behavior, described herein using linear elastic fracture mechanics (LEFM) theoretical analysis as well as numerical and experimental studies. The technology highlights how for herringbone architectures below a critical angle, the material can effectively trap developing cracks along a cleavage plane of the architecture, driving the crack in a direction of high fracture energy (toughness). These herringbone patterns can be extended to enhance fracture toughness in the second in-plane dimension through patterning herringbones in two-dimensional mosaic patterns. Patterning of composites in a herringbone architecture can provide a straightforward way to boost the overall material toughness.
The experimental measurements, described below, showed that near the optimal angle, initial fracture toughness demonstrated substantial improvements including average increases of 185% and 549% in GC compared to the monolithic samples of αΓ=π/2 and αΓ=0, respectively. The experimental observations combined with the numerical simulations, also highlight that the optimal structure needs to be chosen not too close to the theoretical LEFM optimal angle due to their augmented susceptibility to T-stress instabilities.
The composite materials with the herringbone patterns described herein can provide an increase in toughness 3 to 5 times greater than that of composite materials with conventional filler alignments. The technology can be applied to 3D composite printing to effectively toughen parts independent of local geometries and applied stresses. The technology can significantly increase the initial fracture resistance of a polymer system reinforced with filler. In some embodiments, the technology can provide increased fracture toughness isotropically, such that the direction or orientation of the pattern relative to external load and geometry does not matter significantly. The technology can be advantageous since toughness can be a difficult performance metric to increase.
The technology can be widely applicable and transferrable. The technology can be applied in stereolithography apparatus (SLA) 3D Printing of fiber reinforced composites. The technology can be applied in carbon fiber composites industries that employ milled carbon fiber as an additional toughening agent. It can be applied to fiber glass systems, carbon fiber systems, and ceramic matrix composite systems. The technology can be applied in injection molding flows for filled systems. The technology can be applied to provide isotropic toughening of SLA printed composites and toughening of interlaminate materials. The technology can be applied to enhance the toughness in discontinuous fiber reinforced composites, in some case by three times or more. The technology can be advantageous, because toughness can be a difficult performance metric to increase.
By way of further explanation, a variety of mechanisms have been proposed to increase the initial resistance of a polymeric system to the formation of a fracture, as well as its capability to impede the fracture's propagation through the material. Many of these efforts have taken the form of adding filler particles possessing a higher mechanical strength and stiffness to the weaker bulk polymer, increasing its resistance to tensile or mode I fracturing through the transference of mechanical stress from the polymer to the dispersed particles within the composite. A general understanding is that orienting fillers parallel with the principle stress direction is the best way to increase material toughness. The technology described herein provides a further increase in toughness by introducing conflicting alignment patterns of filler within the composite. These patterns range from herringbone patterns for one dimensional (1D) toughening to more complex patterns including herringbone weaves for two dimensional (2D) and three dimensional (3D) toughening. Such patterns can be fabricated with 3D magnetic printing techniques, although a variety of fabrication techniques can be used.
Specifically, the orientation of microscale reinforcement plays a role in the anisotropic fracture toughness of such structures. For example, shear-aligned discontinuous fibers in glass-filled nylon formed by injection molding and weave-aligned continuous carbon fiber epoxy can demonstrate toughness in the direction perpendicular to the reinforcement particles, tending to deflect cracks away from this axis based on maximizing the energy release rate. Therefore, deflected cracks no longer propagate perpendicular to the reinforcing fibers, meaning that the microstructure is not toughening to its full potential.
An ability to finely tune composite microstructures can enable tuning macroscopic part performance through the local control of discontinuous fibers. Complex structures can be possible by a variety of manufacturing techniques, allowing for significant increases in the resolution of sample features and precision control over the filler materials. Selective irradiation of voxel-specific areas, combined with a magnetic alignment system and fillers pre-doped with ultra-high magnetic response iron oxide nanoparticles, as in Martin et. al., can be used to construct 2D and 3D structures containing multiple regions of differing alignments. (Martin, J. J.; Fiore, B. E.; Erb, R. M., Designing bioinspired composite reinforcement architectures via 3D magnetic printing. Nature Communications 2015, 6 (1), 8641.) Other techniques, such as those employed by Ferrand et al., have achieved high-precision alignments via the use of a unique form of slip-casting, in which a fluid suspension of particles can be deposited into a pre-designed, porous mold, utilizing capillary forces to remove the fluid and magnetic fields to control the particle alignment. (Le Ferrand, H.; Bouville, F.; Niebel, T. P.; Studart, A. R., Magnetically assisted slip casting of bioinspired heterogeneous composites. Nature Materials 2015, 14 (11), 1172-1179.)
In the technology described herein, three-dimensional (3D) magnetic printing techniques are used to fabricate thermoset polymers reinforced with discontinuous ceramic particles ordered in highly-oriented “herringbone” architectures. The angular contrast in this structure creates an optically antagonistic pattern termed “herringbone” from the angled skeletal structure of the herring fish.
In some embodiments, a composite material can comprise a volume of a matrix material and a plurality of particles disposed within the matrix material. The particles can have a geometric shape with a longest dimension and a shortest dimension different from the longest dimension. The volume can have an anisotropic structure comprising a plurality of volume portions, adjacent volume portions abutting along an interface comprising a surface forming a common boundary between the adjacent volume portions. The particles within any pair of the adjacent volume portions can be disposed in orientations in which the geometric shape with the longest dimension is disposed at a non-orthogonal angle with respect to the interface, with the orientation of the particles in one of the pair of the adjacent volume portions differing from the orientation of the particles in the of the pair of the adjacent volume portions. The orientations of the particles in the pair of the adjacent volume portions can be selected so that a crack propagating on a crack propagation path through the one or the other of the pair of the adjacent volume portions stops at or deflects to propagate along the interface. This anisotropic structure can be extended in two dimensional and three dimensional directions. The volume portions can have a dimension that is at least 10 times a dimension of the particles.
In some embodiments, at least a subset of the plurality of the volume portions are arranged in the mosaic pattern to be adjacent to at least three other volume portions. In some embodiments, the non-orthogonal angles between the particles within the adjacent volume portions and the interface between the adjacent volume portions are a function of a fracture energy anisotropy of the composite material. In some embodiments, the non-orthogonal angle of each of the volume portions is selected such that a fracture energy anisotropy of the composite material is equal to or greater than 3. In some embodiments, the anisotropic structure provides a reduction of a driving force of a crack tip regardless of crack initiation location or direction as the crack tip propagates through the volume. In some embodiments, the anisotropic structure provides a mechanical property tuned along a selected direction. In some embodiments, the anisotropic structure provides compressive loading reinforcement, tensile loading reinforcement, or both compressive loading and tensile loading reinforcement. In some embodiments, the anisotropic structure provides increased fracture toughness along the interface.
More particularly,
This architecture can be extended into two-dimensions and three dimensions through tiling the herringbone pattern into various mosaic patterns. Referring to
In some embodiments, a herringbone alignment architecture can employ a Bouligand structure, which can provide effective reinforcing techniques for certain loading conditions. In particular, Bouligand structures can be altered to accommodate a series of cells or volume portions, each with reinforcing particles aligned in a particular orientation. The interfaces between each cell or volume portion can effectively pin a crack or drive propagation of the crack along the interface. In some embodiments, the cells can be arranged with the herringbone orientations rotating in a descending or ascending manner. In some embodiments, a composite material can be provided in stacked layers, each layer composed of rows of stacked herringbone architecture. In some embodiments, the architectures can combine both compressive and tensile loading reinforcement.
A further example geometry for a 3D pattern is that of a Wigner-Seitz cell 42, illustrated in
In order to calculate the fracture energy required to initiate a fracture in this “pinning” architecture, the angle at which the crack is anticipated to propagate immediately following initiation (θ*) is first determined. This angle, generally referred to as the “kinking angle,” is determined by the competition between the rate of energy released by the crack tip G(θ) (i.e., the energy of the bonds severed to produce a new surface) and the fracture energy of the crack tip itself, Γ(θ). This initial direction, therefore, is whichever “kink” angle produces the greatest decrease in the energy dissipated by the fracture. The resulting crack's fracture energy is a function of both this kinking/deflection angle and the material's so-called fracture energy anisotropy, , a value determined by the ratio of the fracture energy for platelets aligned parallel to the direction of Mode I applied stress
to the fracture energy of platelets aligned perpendicular to Mode I applied stress (˜0 radians), expressed as
These methods for solving these equations can be expanded to determine the approximate fracture energy values for a variety of particle alignment angles (±αΓ) ranging from parallel to perpendicular global alignments. In the case of the herringbone architecture, these calculations account for the effects of two separate, though symmetric, alignment architectures, with the total initial fracture energy being a function of the system's estimated kinking angle, particle alignment angle, and the fracture energy anisotropy. This allows for the approximation of the energy released by the fracture initiation Gc, as well as estimating a range of fracture energy values for the full gamut of possible particle alignment angles.
The volume percentage of particles to matrix material can be selected to achieve a desired fracture energy anisotropy. The volume percentage can depend on the materials selected and the size of the particles. Particle size can range from nanoparticles to microparticles to macroparticles. As use herein, nanoparticles include particles having a longest dimension in the range of 1 to 1000 nm (10−9 m to 10−6 m). Microparticles include particles having a longest dimension in the range of 1 to 1000 μm (10−6 m to 10−3 m). Macroparticles include particles having a longest dimension above 10−3 m.
In one example, a macro-scale anisotropic fracture energy emerges for a structure of several centimeters in size that is greater than 3000 times a particle size. In one example, using Al2O3 particles having a dimension of 7.5 μm dispersed in a polymer resin matrix material formed from EBECRYL 230 urethane acrylate and isobornyl acrylate, a volume percentage of ˜4% was found to achieve an anisotropy of 3.1.
The alignment architectures, and the methods used to produce them, can have applications in the study of mechanical principles and the production of further methods of creating common commercial products of superior robustness and longevity.
The results from experiments, discussed further below, point to a variety of additional possibilities in the realms of material strengthening. The ability to align filler particles into precise, complex internal architectures has shown an ability not only to elevate a polymer's initial resistance to fracture above that of monolithically aligned or randomly aligned systems, but to create patterns that can overcome the normally anisotropic nature of this strengthening effect.
The full potential of the principles uncovered during these herringbone architecture studies can be achieved through the design of further sample geometries which can apply the herringbone strengthening along all possible degrees of freedom.
The composite material can be fabricated from a variety of matrix materials and particles. In some embodiments, the composite material can be a polymer matrix material and the particles can be reinforcing particles of a ceramic or other material. In some embodiments, the reinforcing particles can be aluminum oxide platelets.
In some embodiments, the matrix material can be a thermoplastic polymer, polyurethane, silicone, acrylic, ABS, nylon, PLA, polycarbonate, polyethylene, polypropylene, polystyrene, polyvinyl chloride, polytetrafluoroethylene. In some embodiments, the matrix material can be collagen, gelatin, or alginate. In some embodiments, the matrix material can be a combinations of materials.
In some embodiments, the particles can comprise platelets, filings, fibers, flakes, rods, whisker, microbeads, or nanoparticles. In some embodiments, the particles can be anisotropic in shape in at least one dimension. In some embodiments, the particles can have a longest dimension ranged from 200 nm to 1000 μm. In some embodiments, the particles can have an aspect ratio ranging from 2 to 200.
In some embodiments, the particles can be a magnetically responsive material. In some embodiments, the magnetically responsive material can be a ferromagnetic material, a paramagnetic material, a superparamagnetic material, iron oxide, iron, cobalt, nickel, an iron alloy, a cobalt alloy, or a nickel alloy. In some embodiments, the magnetically responsive particles can comprise a non-magnetic material coupled with a magnetic material. In some embodiments, the non-magnetic material can be aluminum oxide, calcium phosphate, copper, glass, calcium sulfate, nylon, polystyrene, or silicon carbide. In some embodiments, the non-magnetic material can comprise discontinuous fibers, rods, platelets, flakes, or whiskers. In some embodiments, the non-magnetic material is coated with the magnetic material.
A variety of three-dimensional manufacturing techniques can be used to fabricate the composite materials described herein. Suitable manufacturing techniques include, without limitation, additive manufacturing (such as three-dimensional printing techniques), molding, and tape casting, In some embodiments, a peel-off stereolithographic printer system, capable of selective voxel irradiation and three dimensional magnetic alignment, can be used to create such structures. This permits the creation not only of additional three dimensional structures, but customized interior filler alignments which can provide finely tuned mechanical properties along any axis desired. Such manufacturing techniques provide for rapid prototyping and fabrication of structures incorporating herringbone and other patterns and geometries that provide a variety of internal interfaces and inhomogeneities and that can incorporate a variety of mechanical and material properties. Further suitable three-dimensional manufacturing techniques are described in WO 2015/188175, entitled “Additive Manufacturing of Discontinuous Fiber Composites Using Magnetic Fields,” and WO 2017/100271, entitled “Direct Write Three-Dimensional Printing of Aligned Composite Materials,” both of which are incorporated by reference herein.
This technology provides polymeric composites reinforced with alumina micro-platelets oriented in monolithic and herringbone architectures. It was demonstrated that the herringbone architectures resulted in remarkable toughening behavior using LEFM theoretical analysis as well as numerical and experimental studies. The technology highlighted how for herringbone architectures below a critical angle, the samples could effectively trap developing cracks along the cleavage plane of the architecture, driving the crack in a direction of high fracture energy (toughness). The experimental measurements showed that near the optimal angle, initial fracture toughness demonstrated substantial improvements including average increases of 185% and 549% in GC compared to the monolithic samples of αΓ=π/2 and αΓ=0, respectively. The experimental observations combined with the numerical simulations, also highlighted that an optimal structure can be chosen not too close to the theoretical LEFM optimal angle due to their augmented susceptibility to T-stress instabilities. These herringbone patterns can be extended to enhance fracture toughness in the second in-plane dimension through patterning herringbones in two-dimensional mosaic. Patterning of composites in a herringbone architecture can provide a straightforward way to boost the overall material toughness. Further, herringbone reinforcement as described herein can be achieved using techniques other than the 3D magnetic printing technique as described herein.
The technology described herein can provide conflicting alignment patterns that have been investigated and demonstrated a 3-5× toughness improvement over other conventional filler alignments. The technology can be applied to 3D composite printing to effectively toughen parts independent of local geometries and applied stresses. The technology has been confirmed to significantly increase the initial fracture resistance of a polymer system reinforced with filler. The technology can perform well isotropically; direction/orientation of pattern relative to external load and geometry does not matter significantly.
The crack path in an anisotropic solid is controlled by the competition between the energy release rate from the crack tip and the fracture energy. As such, in most simple anisotropic structures (for high enough anisotropy), the crack driven against the tough direction deflects. The technology herein provides a strategy for designing a tough architecture by forcing the fracture to propagate along the tough path. The effective manufacturing of such hierarchical composites is demanding both in terms of the ability to control a microstructure's orientation and its chemistry.
The herringbone architecture provides a minimal design by which the crack is trapped at the boundary of two mirror-symmetric platelet orientations. This is achieved by linear elastic fracture mechanics (LEFM) calculations where the platelet angles are optimized to enforce propagation in the toughest direction possible without deflection. To calculate the platelets' orientation, the deflection angle of the crack propagating at the mirror-symmetry plane of the herringbone architecture is calculated. It is assumed that the deflection angle is governed by maximum energy dissipation such that
θ*=argmax(G(θ)−Γ(θ)) (1)
It is assumed that the elastic properties of the herringbone composite remains isotropic and use closed-form expressions to calculate the energy release rate from stress-intensity-factors as a function of the kink angle where pure Mode-I loading is assumed. Furthermore, it is assumed that the fracture energy anisotropy is of the form:
Γ(θ;)=
√{square root over (sin2(θ)+
−2 cos2(θ)}) (2)
where θ is the angle of the crack with respect to the platelets' tangent, is the maximum fracture toughness attainable at fracture energy anisotropy
, and
=Γ(π/2)/Γ(0) is the fracture energy anisotropy.
Using Eq. 2, the fracture energy can be constructed for a Herringbone architecture where the platelets are rotated by ∓αΓ (counterclockwise for θ>0 and clockwise for θ<0) from the loading direction
ΓHB(θ;,αΓ)=Γ(|θ|+αΓ,
) (3)
The kink angle can be found for given fracture energy anisotropy and for all platelet orientations π/4≤αΓ≤π/2, where αΓ=π/2 corresponds to a monolithic perpendicular architecture with platelets parallel to the loading direction and αΓ=π/4 corresponds to herringbone structure with platelets at oriented at ∓π/4. The kink angle can then be calculated by varying the load G(0)/ between 0 and 1 and finding the first instance of propagation where G(θ*)=ΓHB(θ;
,αΓ).
≤3.1, the cracks do not kink. For higher fracture energy anisotropy, the cracks kink in the perpendicular architecture, but if the platelets are rotated backward, the kinking is eventually suppressed, and the crack is trapped at the mirror-symmetry boundary. This can be intuitively understood since in the herringbone architecture, θ=0 is a cleavage angle; thus for the crack to kink, it has to propagate in a tougher direction.
at different
values. Here to compare values at different fracture energy anisotropies, perfect toughening is assumed where the maximum fracture energy at given anisotropy scales as
=
Γmin. Simply put this means that the fracture energy for propagation along the platelets stays constant and increases linearly with A for propagation perpendicular to the platelets. As expected, the maximum dissipation through facture is attained when the crack is trapped at the cleavage angle. It should be noted, however, that toughening never reaches the maximum value
Γmin since the trapped crack effectively propagates against an ensemble of randomly tilted platelets at ±αΓ.
The theoretical findings can be verified using phase-field simulations. The phase-field method for fracture couples the displacement field to a smeared fracture field. The fracture field varies smoothly between ϕ=1 in pristine material to ϕ=0 at the cracks. For a domain ⊂n, first introduce the free energy
as the sum of the bulk stored energy
and the surface energy S required for creation of new cracks Γ⊂Ω:
(u,ϕ):=
(u,ϕ)+S(ϕ) (4)
The bulk elastic energy can be written as
(u,Γ):=½∫Ωg(ϕ)CijklϵijϵkldV (5)
where ϵij=(ui,j+uj,i)/2 is the linearized strain tensor and Cijkl=λδijδkl+μ(δikjδjl+δilδjk) is the elastic constitutive tensor and λ=Ev/(1−v2) and μ=E/(2(1+v)) are Lame parameters for elastic modulus E and Poisson ratio v for plane-stresses.
Furthermore, the fracture energy can be written as
where ξ is the process zone size, and Cw=∫01√{square root over (w(ϕ))}dϕ is a scaling constant. To recover the Herringbone fracture energy, the anisotropy tensor D is set as
is the rotation matrix. In this simulation, the KKL model is used, i.e., g(φ)=3ϕ4−4ϕ3, w(φ)=1−g(ϕ). (Murakami, Y., Material defects as the basis of fatigue design. International Journal of Fatigue 2012, 41, 2-10.)
Finally, the fracture propagation is posed as a variational problem where at each time step the displacement field u and the phase-field φ are sought, such that they minimize the approximate free energy:
u
t* ,ϕt*=argmin (u,φ) ∀u∈At, ϕt≥ϕt−1 (9)
where At is the admissible set of displacement satisfying the Dirichlet boundary conditions u=ût on ∂DΩ. The Euler-Lagrange equations derived variationally from Eq. (9) are discretized using the Galerkin finite element method and solved using distributed data structures provided by libMesh and linear algebra solvers in PETSc. (Balay, S.; Abhyankar, S.; Adams, M. F.; Brown, J.; Brune, P.; Buschelman, K.; Dalcin, L.; Dener, A.; Eijkhout, V.; Gropp, W. D.; Karpeyev, D.; Kaushik, D.; Knepley, M. G.; May, D. A.; McInnes, L. C.; Mills, R. T.; Munson, T.; Rupp, K.; Sanan, P.; Smith, B. F.; Zampini, S.; Zhang, H.; Zhang, H. PETSc users manual. Technical Report ANL-95/11; 2019.)
The phase-field simulations were performed in a sample geometry shown in =4 (corresponding to the experimental samples) propagating in the herringbone structure. As predicted by the LEFM results (
=4 the crack is trapped at the symmetry line, αΓ=9π/20 offers the best initial resistance to fracture. This verifies the closed-form LEFM analysis where it can be seen that αΓ=9π/20 is closer to optimal tilt angle line (dashed line in
Finally, the importance of the optimal tilt angle can be highlighted by comparing the load-displacement curves for the herringbone architecture for different fracture energy anisotropies. As shown in =4 can perform as well as the tougher composite in at a less optimal architecture (αΓ=π/4).
To investigate the anticipated mechanical enhancement of herringbone architectures, a series of samples were 3D magnetic printed with acrylate-based photo-polymers reinforced by 5 vol % alumina micro-platelets (corresponding to a fracture energy anisotropy of =4). Specifically, herringbone architectures were printed with two alignment zones below and above the spine with respective orientation angles of αΓ=±0, αΓ=±π/4, αΓ=±π/3, αΓ=±7π/18, and αΓ=±π/2, as shown in
Samples were tested under tension in an Instron Universal tester with the tension grips clamping the samples behind the notch (hatched region in
For the two monolithic sample groupings tested, the αΓ=±π/2 samples exhibited a 43% increase in stiffness and a 127% increase in GC compared to the αΓ=±0 samples. This monolith comparison is consistent with both the theoretical predictions as well as conventional guidelines for “optimal” composite design in which reinforcement is supposed to be oriented orthogonal to the notch. In contrast, herringbone architectures show a surprising increase in GC over both monolithic orientations. In the case of the αΓ=π/3 herringbone architecture, the samples demonstrated an improvement to GC with average increases of 185% and 549% in GC compared to the monolithic samples of θ=π/2 and θ=0, respectively, shown in ≃4 can be predicted for the experimental samples. For
≃4, the LEFM analysis predicts an optimal tilt angle for the herringbone architecture αΓ≃0.47π compared to αΓ≃0.33π observed experimentally. This suggests that the precise form of the fracture energy anisotropy is the main source of this discrepancy between the theoretical LEFM predictions and the experimental results. A simple two-fold fracture energy anisotropy function was used in the LEFM analysis and phase-field simulations. While this simple form may be different from the true fracture energy anisotropy of these composites, the theoretical framework, including the existence of an optimal herringbone architecture remains valid as experimentally demonstrated.
Other herringbone architectures showed similar trends. Subtle gains were found in the αΓ=±π/4 herringbone architecture in which samples demonstrated an average increase of 9% and 148% in GC compared to the monolithic samples of αΓ=±π/2 and αΓ=±0, respectively. All cracks were found to propagate straight along the symmetry line (cleavage plane) of the herringbone architecture. As the cracks were pinned in this direction, the samples exhibited large areas under the normalized force versus displacement curves. Experimental samples with αΓ=±7π/18 herringbone architectures (that sit close to the optimal tilt angle value predicted by the theoretical LEFM analysis) demonstrated crack kinking. As previously highlighted by the phase-field simulations, architectures close to the optimal theoretical tilt angle are more sensitive to T-stress instability and kinking. While in the numerical simulations kinking only occurred after propagation, heterogeneities in the experimental sample can result in expedited kinking as experimentally observed. Heterogeneities can include local inconsistencies in the alignment or concentrations of the ceramic particles along the symmetry line of the herringbone pattern. In the experiments, these αΓ=±7π/18 samples on average showed a much noisier and reduced fracture energy Gc compared to the other herringbone architectures highlighting the reduced stability of the crack path. This noise may have caused most samples in this group to behave like the monolithic sample of αΓ=±7π/2. Nevertheless, the αΓ=±7π/18 herringbone architectures still demonstrate enhanced total toughness, taken as the area under the curve relative to their monolithic counterparts.
As demonstrated through these experiments and predicted by the theoretical LEFM results, the herringbone architectures demonstrate the ability to significantly improve the toughness of the printed specimens. GC and the areas under the normalized force versus displacement are both strong functions of the center-line alignment angles of the herringbone architecture with experimental values peaking around αΓ=±π/4. Further reduction of experimental heterogeneities should result in pushing the optimal herringbone alignment closer to the numerical predictions, effectively increasing further the GC of herringbone composites.
Leveraging the remarkable increases to fracture toughness in herringbones, these principles can be applied to make composites that exhibit more isotropic enhanced fracture toughness. In the simple herringbone patterns of
Preparation of magnetized alumina micro-platelets: The creation of precise alignment architectures made from aluminum oxide required a high degree of precision in particle alignment. A variety of techniques to achieve this effect exist, but a great degree of precision and customizability can be produced through the application of external, alternating magnetic fields. In order to make these particles susceptible to these fields, it was necessary to coat them with a material with a high degree of magnetic response. Commercially available iron oxide nanoparticles provided this responsiveness (often referred to as “Ultra High Magnetic Response” particles), displaying an exceptional level of paramagnetic response to fields induced by external solenoid devices, allowing for the aluminum oxide to achieve a 3D alignment with fields alternating polarity at a sufficiently high frequency. To coat the aluminum oxide platelets with the iron oxide, a suspension of 10 grams of platelets was placed in a suspension of ˜1000 mL of deionized water. This suspension was placed on a hot plate with a rotating magnetic field with a stir bar placed within, rotating at ˜600-700 rpms to produce a mixing vortex within the suspension. This was followed by the addition of 375 micro-liters of iron oxide ferrofluid (nanoparticles suspended in water) to the mixing alumina. This solution was then allowed to mix for ˜24 hours, allowing the iron oxide particles to adsorb to the surface of the aluminum oxide platelets. This is followed by a filtering process, in which the solid platelets were separated from the water suspension, and a ˜12 hour exposure to 90° C. in a convection oven to evaporate any remaining water.
Preparation of photocurable resins: Samples were prepared using a 1:9 mass ratio of polyeurethane and isobornyl acrylate, mixed with a 2% mass percentage of Phenylbis (2,4,6-trimethylbenzoyl) phosphine oxide and 1-Hydroxycyclohexyl phenyl ketone photoinitiates. This mixture was hand mixed, followed by 2 minutes of tip sonication to break up any material aggregates, and 20 minutes of bath sonicated to disperse the initiate powders within the polyeurethane-IBOA system, finalizing the resin which forms the bulk matrix of the composite polymer. Once the pure polymer resin preparation was finished, a 5% volume fraction of the aluminum oxide platelets was added relative to the resin and again tip and bath sonicated for similar times to break up aluminum oxide aggregates and disperse them uniformly within the system.
3D magnetic printing and pre-notching of samples: Following the final mixture of the composite resin system, the material was inserted between two glass substrates. The substrates were suspended over a single, static Z-axis solenoid and between two X-Y axis alternating field solenoids, allowing for alignment of the super-paramagnetic aluminum oxide platelets while the bulk resin was still in a liquid phase. Following 15 seconds of magnetic alignment, a projector connected to an external laptop displayed an image projecting the desired section of the sample under production in ˜405 nanometer wavelength UV light. This initiated material crosslinking in the area under UV exposure, transitioning the desired section of the sample into a solid with the desired filler alignment. This step was repeated for each of the sample sections for the remaining desired alignments, creating a single, continuous fracture coupon with several different areas of uniform alignment. This was followed by a cleaning of the surface with diluted isopropyl alcohol to remove non-polymerized or partially polymerized resin from the surface of the test sample. Following the final cleaning step, the sample was exposed to a higher-power projector for 30 minutes (15 minutes each side) in order to complete the system polymerization, reducing sample ductility. They were then placed in an oven at 90° C. for 2 hours in order to reduce stresses in the sample. 20 minutes into this exposure, the samples were temporarily removed to extend the pre-engineered notch by 1 mm in order to produce a consistent starting location for the sample fracture for all tensile tests. After the 2 hours elapsed, the samples were removed from the oven and allowed to drop below their glass transition temperature, returning to their original ductility and ready to be fractured.
Mode I fracture tests: Once sample preparation was completed, fracture testing was carried out using an Instron Tensile Tester. Sample ends were held in place using Instron clamps, with one clamp remaining stationary while the other moved upwards from its zeroed position at a pre-selected strain rate. In order to limit the effects of out-of-plane deformation which might cause fracture of differing modes besides that of Mode I, the samples were placed between a pair of glass slides, separated by a set of spacers to ensure the slides and samples were not in contact initially, and held in place by an external clamp. In addition, to ensure that contact between the samples and slides during testing would not affect the results due to frictional concerns, the interior of the slides was covered with soap-water to reduce the force due to the sample dragging across the glass surface.
Prior ex-situ studies have demonstrated the strain-rate dependence of polymeric systems behavior, with lower strain rate/higher temperature tensile conditions producing viscoelastic response and higher strain rate/lower temperature conditions producing more brittle system responses. As a result, samples were exposed to tensile force conditions at strain rates of 20 mm/min, which prior in-situ lab tests had shown to be sufficient to produce non-viscoelastic mechanical results. The test was continued until an initiated fracture propagated throughout the entire sample and was manually discontinued. Fracture results were then output as a set of force versus displacement data, which was later normalized using the average thickness of the sample under test. Additionally, tests were video recorded to aid in determining which portions of the force versus displacement curves were associated with different fracture behaviors during testing. This allowed to ensure the determination of exactly what point of the curve could be said to be associated with the fracture initiation and where this fracturing behavior could be said to end.
Determination of G: To calculate the critical energy release rate Gc, elastic finite element analysis was performed of the experimental specimen geometry (
Simple dimensional analysis showed that the energy release rate from the tip of a crack assuming LEFM, for displacement controlled, can be written as:
where E is the elastic modulus of the material, L0 is a characteristic dimension used to adimensionalize the sample for the computer simulation (here, take L0=10 mm), δ is the applied displacement, and
Furthermore, the elastic modulus of samples E=1.13 Gpa was obtained by fitting the numerical stiffness of the samples to the experimental stiffness obtained load-displacement curves. Combined with the polynomial fit to constant
G=(3.07−0.7√{square root over (Lc)})δ2 (10)
Determination of initial fracture energy: Assuming the Griffith criterion G=Gc the fracture energy was calculated using Eq. 10 as the critical energy release rate for the critical load δc at onset of fracture propagation for the initial pre-notch Lc≃1 mm.
As used herein, “consisting essentially of” allows the inclusion of materials or steps that do not materially affect the basic and novel characteristics of the claim. Any recitation herein of the term “comprising,” particularly in a description of components of a composition or in a description of elements of a device, can be exchanged with “consisting essentially of” or “consisting of.”
To the extent that the appended claims have been drafted without multiple dependencies, this has been done only to accommodate formal requirements in jurisdictions that do not allow such multiple dependencies. It should be noted that all possible combinations of features that would be implied by rendering the claims multiply dependent are explicitly envisaged and should be considered part of the invention.
The present technology has been described in conjunction with certain preferred embodiments and aspects. It is to be understood that the technology is not limited to the exact details of construction, operation, exact materials or embodiments or aspects shown and described, and that various modifications, substitution of equivalents, alterations to the compositions, and other changes to the embodiments and aspects disclosed herein will be apparent to one of skill in the art.
This application claims priority under 35 U.S.C. § 119(e) of U.S. Provisional Application No. 62/985,966, filed on 6 Mar. 2020, entitled “Ceramic-Reinforced Polymer Composites with Herringbone Architecture,” and U.S. Provisional Application No. 63/124,187, filed on 11 Dec. 2020, entitled “Ceramic-Reinforced Polymer Composites with Herringbone Architecture,” the disclosures of which are hereby incorporated by reference.
This invention was made with government support under Grant Number 1536354 awarded by the National Science Foundation. The government has certain rights in the invention.
Number | Date | Country | |
---|---|---|---|
62985966 | Mar 2020 | US | |
63124187 | Dec 2020 | US |