Next Article in Journal
Bone Regeneration Using Rat-Derived Dedifferentiated Fat Cells Combined with Activated Platelet-Rich Plasma
Next Article in Special Issue
Quantification Approaches for Fatigue Crack Resistance of Thermoplastic Tape Layered Composites with Multiple Delaminations
Previous Article in Journal
The Significance of Coherent Transformation on Grain Refinement and Consequent Enhancement in Toughness
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fracture of Thin-Walled Polyoxymethylene Bulk Specimens in Modes I and III

1
Institute for Mechanics and Materials, Technische Hochschule Mittelhessen, 35390 Gießen, Germany
2
Materials Science and Testing of Polymers, Montanuniversität Leoben, 8700 Leoben, Austria
3
Polymer Competence Center Leoben GmbH, Roseggerstraße 12, 8700 Leoben, Austria
*
Author to whom correspondence should be addressed.
Materials 2020, 13(22), 5096; https://doi.org/10.3390/ma13225096
Submission received: 16 October 2020 / Revised: 4 November 2020 / Accepted: 7 November 2020 / Published: 12 November 2020
(This article belongs to the Special Issue Fracture Mechanics Investigation of Polymeric Materials)

Abstract

:
Thin-walled polymeric components are used in many applications. Hence, knowledge about their fracture behavior in bulk is beneficial in practice. Within this study, the double cantilever beam (DCB) and out-of-plane double cantilever beam (ODCB) tests are enhanced to enable the testing of such bulk specimens in mode I and mode III on the basis of the J-integral. This paper then presents and discusses the experimental results following the investigation of a semicrystalline polymer (polyoxymethylen) under quasi-static load conditions. From the experiments, fracture energies of similar magnitude in both mode I and mode III were determined. In mode III, pop-in fracture was observed. Furthermore, the fracture surfaces were investigated regarding the mode I and mode III dominant crack growth mechanisms, based on the morphology of the tested material. For specimens tested in mode I, no signs of plastic deformation were observed, and the fracture surface appears flat. In mode III, some samples display a twisted fracture surface (twisting angle close to 45 ° ), which indicates local mode I crack growth. A transfer of the presented methodology to other (more ductile) polymeric materials is deemed possible without further restrictions. In addition, the presented setup potentially enables an investigation of polymeric bulk specimens in mixed mode I+III.

1. Introduction

Polymers are used in a wide range of applications, from daily use (packaging, bottles, etc.) to mechanical components (pipes, rollers, gearwheels, etc.). Unfortunately, compared to other materials such as metals or ceramics, there is still a large variety of open questions that remain to be addressed regarding the design of polymeric assemblies. To guarantee the safety of polymeric components, it is beneficial to possess good knowledge of their fracture mechanical properties. Especially when used as tube materials, information about their behavior as thin-walled components is also necessary, as fracture mechanical tests are generally conducted under plain strain conditions and, hence, greater specimen widths are tested.
Due to its material properties such as high stiffness, dimensional stability and fatigue resistance, polyoxymethylene (POM) is currently used in many technical applications [1,2,3]. While the fracture mechanical properties of POM are already available in mode I and mixed mode I+III under monotonic loading and cyclic fatigue tests [4,5,6,7,8], to the authors’ knowledge, little is known about its fracture mechanical properties in pure mode III. This paper addresses the need for the determination of fracture mechanical characteristics of POM for thin-walled components, so far lacking scientific literature.
Generally, only a few methods are available for determining the fracture energy of polymeric materials in pure mode III, of which many have drawbacks either in terms of experimental effort, costs or unwanted contributions in other modes that cannot be deemed negligible [9]. Furthermore, previous studies and therein proposed experimental setups have exclusively focused on linear elastic fracture mechanics, limiting their field of application to brittle and quasi-brittle materials. Hence, a number of questions regarding the mode III fracture of ductile polymers remain to be addressed.
To determine the fracture mechanical properties of POM and possibly other polymers in both pure modes I and III, we propose experimental setups based on the double cantilever beam (DCB) and out-of-plane double cantilever beam (ODCB) tests evaluated by using the J-integral. DCB tests with this method of data reduction have historically mostly been used for the experimental investigation of adhesives or composites, see [10,11,12,13] for example, whereas the ODCB test was exclusively used for the investigation of adhesives thus far [14]. The proposed modifications to the DCB and ODCB tests offer the possibility to determine the fracture energy in pure modes I and III on the basis of the J-integral. Hence, non-linear fracture behavior can also be characterized with the herein proposed setups, implying that the presented experiments are also used for materials that are not necessarily brittle or quasi-brittle. Conveniently, for the special case of linear elastic and brittle or quasi-brittle materials, the obtained J-integral measurements can be directly converted to approximate the fracture toughness values K I c and K I I I c .
This paper begins by examining the theoretical background of the experimental determination of the J-integral in pure modes I and III. We will then describe the modifications made to the DCB and ODCB tests, which enabled the determination of the energy release rate (ERR) of a POM homopolymer (POM-H) in modes I and III. Then, after the determination of the ERR in modes I and III, a comparison with relevant literature is sought to classify the generated results in mode I and III. After an examination of the fracture surfaces and methodological critique, we will then present a summary of the most important results.

2. Theoretical Background

The ERR G is central to the field of fracture mechanics. It can generally be defined as the decrease in potential energy W per increase in fracture surface area, yielding
G = Δ W b Δ a
for plane problems. Here, a denotes the crack length, and b is the out-of-plane thickness of the body. For a crack to propagate, the ERR must equal the critical ERR G = G c .
Next to the ERR, the J-integral according to Rice [15]
J = S W d y t i δ i x d S
represents an alternative approach to determine the release of energy during fracture. Here, t i are the components of the traction vector, δ i are the components of the displacement vector along an arbitrary path S containing the crack tip in counter-clockwise direction, and x is the direction of crack propagation (see Figure 1 for reference). This two-dimensional integral is path independent. It should be noted that the J-integral, according to its initial definition, strictly only applies to materials for which a strain energy density exists, implying that it may only be used for materials that behave in a hyperelastic manner. One should also note that the strain energy density must not explicitly depend on the x-coordinate for J to be path independent. However, the J-integral was also found to be path independent for small scale yielding and even during hardening at the crack tip under a monotonically increasing load [16,17]. An important key factor resulting from this path independence is that contributions from external loads are in equilibrium with the value of J at the crack tip if the path is chosen around the outer bounds of the tested specimen, allowing the determination of J at the crack tip from the outer loading conditions. Furthermore, the J-integral and the ERR are equivalent if the crack grows straight ahead, the deformation at the crack tip is largely linear elastic, and only small-scale yielding is present [18].
In the following, we will briefly present the necessary equations for the evaluation of the conducted experiments in modes I and III and state the most important findings of prior studies. For more detailed derivations, we would like to refer to relevant literature, e.g., [10,11,14,19,20].
Consider a specimen with the corresponding external loads as displayed in Figure 2. Under these given forces F and moments M, the J-integral yields
J = F y θ 1 + θ 2 b + M y κ b + M x , up 2 + M x , low 2 2 b 1 μ I t + M z 2 2 b 1 E I z
with the measured rotational angle at the points of load introduction θ i and the curvature κ . E and μ denote the elastic Young’s and shear modulus of the specimen. I z and I t are the second moment of area around the z-axis and the torsional second moment of area around the x-axis, respectively. By deconstructing the above equation, one obtains the contributions to J in the individual modes:
J I = F y θ 1 + θ 2 b
is the contribution in pure mode I,
J I I I = M y κ b
is the pure mode III contribution,
J I * = M x , up 2 + M x , low 2 2 b 1 μ I t
is an unintended contribution to J by an “out-of-plane mode I”—loading due to the specimen twisting, and
J I + I I = M z 2 2 b 1 E I z
is an unintended contribution in modes I and II due to the finite width of the specimen.
The instantaneous experimental determination of the ERR based on the J-integral is relatively unsophisticated in mode I, as only force and rotational angle at the points of load introduction have to be measured [10]. To obtain the curvature κ in mode III, an increased experimental effort and the use of additional measuring equipment are required. Provided that the load is introduced to the crack tip by linear-elastically deforming beams, this experimental effort can be reduced by using the Irwin–Kies Equation
G I I I = M y 2 2 b d C I I I d a
with the change of compliance d C I I I in shear direction, and the change of the lever arm d a during crack propagation [14]. Loh and Marzi were able to establish a connection between the two evaluation methods from Equations (5) and (8) using Bernoulli beam theory, and determined that J I I I = G I I I for this evaluation method [14]. Their beam theory approach yields
κ = M y 2 d C I I I d a ,
offering the possibility to significantly reduce the experimental effort, if d C I I I / d a is determined beforehand.
Furthermore, it should be mentioned that the unintended contributions due to the mode III load were found to be negligible in the studies of Loh and Marzi [14,20,21], in which adhesives were investigated. Because of the changed specimen geometry we use within our study, it must be investigated to which extent these unintended contributions are present at the point of material failure with our proposed mode III setup to determine the “purity” of the mode III fracture.

3. Materials and Methods

3.1. Tested Material—Polyoxymethylene Homopolymer

The investigated material is a POM-H (Delrin 111PF) from DuPont (DuPont de Nemours, Wilmington, NC, USA) supplied as tubular granules, with a nominal diameter of 4 mm and a nominal height of 2 mm. POM-H shows good mechanical properties (Young’s Modulus E around 3500 MPa and Poisson’s ratio ν around 0.42) [1,22,23,24], high ductility down to 0 ° C, high abrasion resistance and a low friction coefficient. Furthermore, the material is able to resist high amounts of constant loading or fatigue loading. Therefore, POM-H is used in more advanced applications with higher requirements concerning the material properties. POM-H is categorized as a technical thermoplastic material. In preliminary investigations, a glass transition temperature for the tested POM-H of −64 ° C was obtained using dynamic mechanical thermal analysis.

3.2. Specimen Preparation

The geometry of the tested POM-H specimens is displayed in Figure 3. The specimens were milled and cut from plates (dimensions approx. 200 × 200 × 5.5 mm), which were produced by compression molding (Hydrostat 300, Schwabenthan, Germany) with an immersion edge tool (TT-260, Tool-Temp). We want to note that extreme caution was taken during the milling and cutting processes, to ensure minimum temperature yield. For a better overview, the different stages of processing as well as the used tools and milling parameters are disclosed in Appendix A.
Prior to experimental investigation, the specimens were notched via pushing with a thin razor blade. The depth of the notches lay in the range between 1.0 ± 0.4 mm (mean and standard deviation). We want to mention that notching proved to be very difficult because of the brittle material behavior, leading to a larger standard deviation of the achieved notch length.
In the mode III investigation, aluminum reinforcements with a length of 160 mm and the cross-section area displayed in Figure 4 are connected to the POM-H specimens. The aluminum reinforcements were milled from an aluminum alloy (AlZn5.5MgCu, material grade number 3.4365) by a professional supplier (Feiler GmbH, Ehringshausen, Germany). The aluminum reinforcements were added to avoid energy dissipation outside of the crack tip due to plastic deformations in the lever arms caused by the introduction of an external moment. Furthermore, as the aluminum reinforcements are very much stiffer than the POM-H specimens and behave linear-elastically, it is ensured that possible influences of material nonlinearities are repressed so that the evaluation of the mode III ERR remains possible with Equation (8).
The POM-H specimens and the reinforcements were joined using thermal clamping. For this purpose, the specimens were cooled to a temperature of −60 ° C with an ultra-deep freezer Herafreeze Basic (Thermo Electron LED GmbH, Langenselbold, Germany) and, hence, subjected to thermal shrinkage. After two to three hours of cooling, the POM-H specimens were removed from the industrial freezer and inserted into the aluminum reinforcements. Then, steel gauge tape with a thickness of 0.20 mm is inserted into the small gaps between specimen and reinforcement on each side. Through thermal expansion of the POM-H specimens, the specimens and the aluminum reinforcements are, hence, uptight at room temperature. It should be noted that the glass transition temperature of POM-H is close to being reached during the cooling procedure, which could, although deemed very unlikely, influence the fracture behavior observed in the experiments. For reasons of space, the investigation of this influencing factor is postponed to possible future works. For a better overview, the specimens in the mode I and mode III investigation are shown in Figure 5.
Overall, 29 specimens were produced (eleven for the mode I investigation, 18 for the mode III investigation).

3.3. Experimental Setups and Test Evaluation

In Figure 6, the test setup as realized in a biaxial tension-torsional servo-hydraulic test machine of type MTS Landmark Bionix (MTS Systems, Eden Prairie, MN, USA) is displayed. For both the mode I and the mode III setup, a distance between the point of load introduction and crack tip of 70 mm was maintained. It should be noted that for every tested specimen, the width of the specimen at the initial crack tip was measured with an optical microscope prior to evaluation. All tests were performed under laboratory conditions (20–23 ° C, approx. 40% RH).
For the mode I tests, the specimen is clamped directly into the test machine without attaching reinforcements. The angles at the points of load introduction θ 1 and θ 2 are measured with two high resolution rotary incremental shaft encoders. A cross head velocity of v T = 0.05 mm/s was selected during the mode I investigation. To determine the critical ERR in mode I, eleven specimens were tested. The tests were terminated after the first brittle crack extension for being able to measure the generated crack length after the experiment.
Next to the determination of the mode I ERR using the J-integral, the modified beam theory approach according to ASTM standard D5528 [25] is also used to determine the ERR. In this method of data reduction, the ERR is calculated from
G I = 3 2 F δ b ( a + Δ )
with the correction term Δ , which is derived from a regression between the cube root of compliance and crack length. Prior to the analysis, the displacement measurement was corrected with the experimentally obtained compliance of the test setup (561.6 Nmm).
In mode III, the specimen is loaded with an angle-rate of α ˙ = 0.05 deg/s. To achieve pure mode III and avoid lateral forces, the setup contains two orthogonal linear slides below the bottom clamping device. Furthermore, axial forces are controlled to be naught by the control system of the test machine, leading to a vanishing mode-I-contribution. To directly obtain the mode III ERR from the J-integral, the curvature of three specimens is measured using two strain-gauges in half-bridge circuit at the external surface of the aluminum reinforcements. With the measured averaged strain ε and the width c of the aluminum parts, the curvature is computed with κ = 2 ε / c . Thereupon, d C I I I / d a is determined from the given measured values in order to enable the calculation of the ERR via Equation (8), which leads to a reduction of experimental requirements, as the external moment is solely needed for evaluation. At last, using 15 more specimens, the critical mode III ERR of POM-H is determined.

3.4. Crack Length Determination and Fracture Surface Analysis

The produced crack length during the experiment in mode I was determined after testing. Therefore, the crack path of the mode I specimens, which were not fully fractured, were covered with a black ink to mark the crack advancement. Afterwards, the specimens were completely broken and the crack length was obtained as the length of the covered area via a digital calliper (Kellner & Kunz AG, Wels, Austria), with a measurement accuracy of 0.02 mm.
To obtain more information about the fracture process, a detailed analysis of the fracture surfaces was conducted. Therefore, an optical microscope SZX12 (Olympus Life Science Europe GmbH, Hamburg, Germany) was used. The investigation of the fracture surface provides valuable information about the fracture mechanisms and crack path formations in modes I and III.

4. Results and Discussion

4.1. Experimental Results in Mode I

In Figure 7, the measured force-displacement curves are displayed. As shown, the relationship between force and displacement is almost completely linear up until the point of a sudden, brittle fracture. Furthermore, the cube root of the specimen compliance is displayed over the crack length with the 5% confidence bands of the performed regression, from which a correction factor Δ of 37.66 mm can be derived. Using the modified beam theory approach according to Equation (10) yields a critical ERR of (6.58 ± 1.09) kJ/m2.
However, the confidence bands suggest that the determination of Δ leads to a considerably large uncertainty in the calculation of G I c in this case. To roughly quantify this uncertainty, the critical ERR and its standard deviation depending on the obtained correction factor Δ is displayed in Figure 8. At the 5% confidence bands of Δ , ERRs between (5.18 ± 0.85) kJ/m2 and (8.13 ± 1.34) kJ/m2 are obtained. It is therefore argued that the calibration of the corrected crack length in this method of data reduction may be a strong source of error, which would be eliminated in the calculation of the J-integral via Equation (4).
The measured mode I J-integral is shown over the cross head displacement δ in Figure 9. From the conducted experiments, a value for J I c of (8.84 ± 1.39) kJ/m2 can be obtained. On the right side of Figure 9, both the maximum values of J I as well as the values at crack arrest J I , rest are displayed over the generated crack surface. Interestingly, these results imply that J I c increases linearly with the generated crack surface, whereas J I , rest linearly decreases. In Table 1, the parameters of both of these linear regressions are summarized. Both regression lines just narrowly fail to meet in the intersection with the ordinate axis. However, within the confidence bands of the regression analyses, it is possible that an intersection of both lines exists on the ordinate axis. One may hypothesize that this intersection is the minimum possible value of J able to cause crack propagation. Considering the overlaps between both confidence bands, the intersection is most likely located in the range between 4.56 and 6.19 kJ/m2.
In practical application, these results can be advantageous: Firstly, knowing the relationships between J I c and J I , rest with the generated crack surface allows a rough determination of the fracture energy of a thin-walled component (with a similar geometry to the tested specimens) in hindsight, by simply measuring the area of the fracture surface. Secondly, the intersection between the regression lines is a useful parameter for the design of a thin-walled polymeric component, as it represents a lower limit to the crack driving force. We suggest that further research is performed in these areas to determine whether the obtained results also remain valid for larger sample sizes, other specimen geometries, and different polymeric materials.

4.2. Experimental Results in Mode III

4.2.1. Determination of d C I I I / d a

Within this study, three specimens are used to experimentally determine the constant d C I I I / d a . In Figure 10, the measured moment M is displayed over the rotational angle of the biaxial testing machine α .
As can be seen in Equation (9), d C I I I / d a can be obtained from the relationship 2 κ / M y . Here, d C I I I / d a is determined from a linear regression between 2 κ and M y up until the measured maximum of J I I I . The value for d C I I I / d a is then derived from the slope of the regression line. The measured values of curvature κ on the outside of the reinforcements and the bending moment M y as well as the performed regression are also displayed in Figure 10. The slope of the regression line equates to (6.75 ± 0.01) × 10−9 1/Nmm2 (±5% confidence interval of regression slope). It can be observed that the measured values are captured by the regression line with great accuracy (correlation coefficient R 2 = 0.999 ).
For verification purposes, d C I I I / d a is also determined analytically with d C I I I / d a = 2 / ( E I y ) , the aluminum’s elastic modulus of E = ( 70 ± 1 ) GPa and the specimen holders area moment of inertia I y = 4116.66 mm4. From this, a value of ( 6.940 ± 0.101 ) × 10 9 1/Nmm2 can be derived, which fits with the experimentally obtained d C I I I / d a , given that the compliance of the aluminum reinforcements itself should be greater than the compliance of the holder with an inserted POM specimen. Because of the relatively good agreement between experimentally measured and analytically determined d C I I I / d a , the evaluation of the unintended contributions according to Equations (6) and (7) is performed with the moments of inertia of the reinforcements.

4.2.2. Critical Mode III ERR

The experimental results with 15 additional specimens are displayed in Figure 11. From this, a critical mode III ERR of (7.59 ± 1.19) kJ/m2 can be derived. We note from the diagram that prior to critical failure all samples show at least one significant, sudden drop in J I I I . This so-called pop-in phenomenon has been mostly observed for steels or weldments in mode I testing. Studies found that pop-ins can both result from a local unstable crack growth that is then stabilized by the surrounding material [26] or by the formation of cracks perpendicular to the plane of the initial pre-crack [27]. To our knowledge, pop-in in pure mode III fracture of polymers is observed for the very first time within this study.
In Figure 12 the ratio between the values of J I I I at which pop-ins occurred and J I I I c is displayed. Here, pop-in ratios between 53% and 93% can be observed. For the given sample size, the results do not allow any statement as to whether the pop-in-ratio is dependent on the measured J I I I c . The occurrence of pop-ins can be critical for the structural integrity of a thin-walled POM-H component loaded in shear as cracks could grow prior to critical failure, potentially weakening thin-walled components to a significant extent. Because pop-ins were already observed shortly above 50% of J I I I c , it is indicated that shear loads could be very much critical for a thin-walled POM-H component.
The measured unintended contributions normalized to the current measured mode III ERR are displayed over J I I I / J I I I c in Figure 13. One should note that the contribution of J I * is considerably larger than the measured mode III ERR at the start of the test. This may be partly attributed to the fact that the breakaway force of the linear slides must first be overcome at the beginning of the test. The load history is therefore not to be considered as a pure mode III loading process. However, it can also be observed that the unintended contributions are, in fact, negligible at the points of fracture. At fracture, J I * takes up (3.02 ± 1.01) % of J I I I and tends further towards naught, whereas J I + I I only takes up (0.04 ± 0.04)‰. This means that, according to the measurements, the fracture process can be considered as a pure mode III fracture.

4.3. Approximate Determination of the Stress Intensity Factors K I c and K I I I c

As a side effect, for the special case of linear elastic, isotropic, and brittle or quasi-brittle materials, the stress intensity factors in mode I and mode III can directly be related to the critical values of G (and J) in a given loading mode. In mode I, the stress intensity factor K I c is calculated with
K I c = E J I c = E G I c
with E = E in plane stress and E = E / ( 1 ν 2 ) in plane strain. Unfortunately, it is unclear at this point, whether the crack tip was loaded in plane stress or plane strain due to the necessary addition of the side-grooves to the specimen. Although the thin geometry of the specimens should induce a plane stress state, the addition of side-grooves is known to induce stress triaxialities that could lead to a state of plain strain [28].
In mode III, i.e., for antiplane shear, the stress intensity factor K I I I c can be derived with
K I I I c = 2 μ J I I I c = 2 μ G I I I c .
As shown in Section 4.1, a critical ERR J I c of (8.84 ± 1.39) kJ/m2 was determined from the experiments performed in this study, which roughly equates to a mode I stress intensity factor K 1 c of 5.6 MPa m1/2 (plane stress) or K I c of 6.1 MPa m1/2 (plane strain), respectively. In mode III, a fracture energy of (7.59 ± 1.19) kJ/m2 was obtained, which equates to a K I I I c of about 4.3 MPa m1/2.

4.4. Summary and Discussion of the Obtained Fracture Mechanical Properties

In this section, we want to shortly summarize and discuss the obtained fracture mechanical properties displayed in Table 2.
In the mode I investigation, a relatively large discrepancy between the experimentally obtained mode I G I c and J I c was found. The authors believe that this discrepancy is due to a large uncertainty in the determination of the corrected crack length Δ from the crack length measurement and specimen compliance. We want to point out that using the J-integral method of evaluation offers the possibility to determine the fracture energy of a polymeric bulk specimen without measuring crack length or the elastic properties of the material. As only force and rotational angle have to be measured for the determination of the mode I J-integral, a simultaneous determination of the fracture energy during the experiment is possible. This allows controlling the experiments on specific values of J, which poses an interesting topic for future research.
Regarding the obtained mode III J-integral, we want to point out two important factors: Firstly, we want to emphasize that the evaluation of the J-integral with the measured curvature κ might be necessary in some cases. With the aluminum reinforcements, a constant d C I I I / d a could be ensured in this study, rendering J I I I and G I I I equivalent. However, this may change if the reinforcements are removed. Furthermore, if the material was more ductile, a determination of J I I I with the curvature κ may prove to be more accurate. Secondly, it should also be mentioned that J I I I can be determined with both evaluation methods ( κ or d C I I I / d a ) instantaneously during the experiment. This, as in mode I testing, allows controlling the experiments on specific values of J in future studies.
A review of literature provided mode I stress intensity factors K I c for POM-H between 2.5 and 6.9 MPa m1/2 [4,6,7]. Thus, the obtained mode I stress intensity factor between 5.6 and 6.1 MPa m1/2 lies within the range of reference measurements from literature. We want to emphasize that our results match with the plane strain fracture toughness of the same material determined with compact tension specimens [6,7]. Unfortunately, we were unable to better investigate the stress state at the crack tip within the scope of this work, and we cannot report the mode I stress intensity factor at a greater accuracy. To the authors’ knowledge, no comparative values are available under pure mode III. Hence, one of the highlights of our study is that the determination of the pure mode III fracture toughness of POM was made possible with our setup for the first time. However, future studies on the topic are suggested in order to verify the determined mode III fracture energy and fracture toughness.

4.5. Investigation of the Fracture Surfaces

4.5.1. Mode I

An optical analysis of the fracture surface was conducted in this study to gain more information about ongoing crack growth mechanisms during testing. This is a common method to investigate the crack growth process after testing and to identify changes in the latter with the surface structure. An overview picture of the fracture of a specimen tested in mode I and a light microscope image of the fracture surface close to the initial notch are shown in Figure 14.
The macroscopically observed fracture surface of monotonically loaded mode I specimens is flat (see overview picture in Figure 14), with three different areas on the fracture surface (marked with (1) to (3) on the images in Figure 14). The first area represents the pre-notch, which was generated via a razor blade before testing. Area (2) marks the produced crack growth during testing and area (3) was generated after testing to determine the ligament length.
The quality of the pre-notch (sharp crack tip without any signs of deformation) has a major influence on the resulting fracture mechanical parameters. Therefore, it is of high interest to guarantee a similar notch quality for all investigated specimens. For the used POM-H specimens, the notching procedure was found to be challenging because of the rather brittle material behavior. This difficulty manifests itself as small imperfections at the crack tip, as shown in Figure 14. Hence, small variations of the produced pre-notch led to an increased standard deviation of the determined fracture mechanical parameters.
No indications of plastic deformation, which would have been visible as intensive white zones on the fracture surfaces, were found in the investigation. Furthermore, as already observed during the experimental investigation, the specimens fractured in a brittle manner. This observation is in contrast to previous monotonic mode I tests on the same material class [6], in which an intensive white zone was found whose formation is attributed to crazes and micro-voids within the material.
In general, deformation rate, specimen thickness, and the added side-grooves also have a large influence on the size of the plastic zone. Based on the fracture surface investigation and the brittle fracture without any indications of plastic deformations, it is suggested that the specimens fractured in a plane strain state. This is also supported by the addition of the side-grooves, which lead to an increased triaxiality along the crack front. However, this finding suggests that further investigations are required to gain a better understanding of the influence of specimen thickness and groove shape on the fracture behavior of POM-H.

4.5.2. Mode III

To examine the fracture process of the mode III samples, a fracture surface analysis via an optical microscope was conducted similar to the mode I procedure. An overview of two representative fracture surfaces loaded in mode III is shown in Figure 15. Interestingly, two different types of mode III fracture surfaces were observed. A larger group of the mode III specimens displayed extensive twisting and crack plane deflection (Figure 15a), whereas in some cases, a less deflected crack flank was observed (Figure 15b).
The twisting of the crack during propagation is known as mode I branching. Here, the crack deflects under mixed-mode conditions to a local mode I loading [29,30]. Such mode I branches were also observed in a previous study on the same material in mode I/III fatigue tests [8]. Furthermore, the occurrence of shear lips and slant crack growth is quite common for thin specimens subjected to mode III [31,32,33]. As the measurements obtained in Section 4.2.2 demonstrated that the fracture was initiated in pure mode III and the fracture surfaces show typical signs of mode III fracture and the subsequent mixed-mode crack propagation, we can conclude that our adaption to the ODCB test allows the determination of the mode III fracture energy without further restrictions.
After a comparison with the results shown in Figure 12, a dependency of the macroscopically observed angle of the fracture surface and the pop-in ratio is suggested. A small ratio of J I I I , pop - in / J I I I c seems to lead to a more flat and straight fracture, as shown in Figure 15b. This means that the macroscopically twisted and deflected fracture surfaces may be related to a higher pop-in ratio. However, to prove this assumption, further investigations are necessary. Especially a determination of the twisting angle at the crack tip using optical methods could be of huge benefit. Apart from the implied change in pop-in ratio, no connection between other experimental results and the twisting angle of the fracture surfaces could be found within the framework of this study.
Akin to the mode I fracture, no indications of plastic deformation were observed on the mode III fracture surfaces. Figure 16 presents a detailed picture of the fracture surfaces of a macroscopically twisted specimen (Figure 16a) and a specimen with a less deflected fracture surface (Figure 16b). The pre-notch is marked as area (1) and the mode III crack growth is marked as area (2).

4.6. Advantages, Limitations and Research Proposals

Before summarizing the most important results obtained within our paper, we would like to discuss the advantages and limitations of our work and propose topics for future studies.
We are aware that our research may have limitations: As we mainly tried to enable the testing of thin-walled polymeric components, we were unable to examine some possible influencing factors in more detail, unfortunately. It is generally well known that both the thickness of the specimen as well as the position of the initial crack tip may have a large influence on the fracture behavior. Furthermore, the stress state at the crack tip is likely influenced by the geometry of the grooves on the side of the specimen. As previously mentioned, the absence of signs of plastic deformation at the crack tip suggests a fracture in plane strain, although the slender specimen geometry should lead to a plane stress state. Hence, the impact of the grooves’ geometry is unclear at this point. In this case, a simulative study using finite element analysis should provide important insights, but as we primarily focused on the experimental setup and the methods of evaluation, we have refrained from performing simulations at this point in time. Such detailed analyses would not have been possible within the framework of this study without neglecting other important areas of our investigations. We propose that further research be undertaken in these areas.
In our investigation, we had to rely on thermal clamping of the specimens. Currently, we also cannot certainly rule out that the cooling procedure has influenced the overall fracture behavior of the specimen in mode III. Furthermore, one could argue that the addition of the aluminum reinforcements influences the load introduction to the crack tip, possibly impacting the measured results. We are aware of this limitation and, hence, propose that further research is undertaken to improve the specimen geometry. An alternative can be a change of specimen thickness at the lever arms, which may allow omitting the aluminum reinforcements in the first place.
As discussed above, the instantaneous determination of J during the experiment in both modes I and III can be used to control the experiments on J. This enables testing in mixed-mode I+III under constant mixed-mode ratios, which can help to better understand the fracture behavior of polymeric bulk specimens. A revision of the sample geometry will therefore be required so that the aluminum holders become obsolete. In future studies, this should be investigated in more detail.
We also want to emphasize that the proposed setups should also be tested with other polymeric materials, as we exclusively focused on POM-H in our research. The selection of POM for this study was mainly due to its high relevance among engineering plastics and our prior knowledge of the material. Another advantage of POM is its limited dependency on strain rate in the elastic range [34], which allows the strain energy density to be considered as rate-independent in good approximation.
One of the obvious advantages of our methodology is that using the J-integral allows for the investigation of ductile materials. As a nice side effect, in case of a brittle or quasi-brittle failure, the stress intensity factors in modes I and III can be determined. Furthermore, we found evidence to suggest that our mode III setup enables pure mode III testing of various kinds of polymers, as our setup allows for a precise experimental determination of the contributions to the fracture process in other fracture modes.

5. Conclusions

Within the limitations of our study, the following conclusions can be drawn:
  • Our proposed methodology provides the possibility to determine the energy release rate of a polymeric material in pure mode I on the basis of the J-integral. It also allows us to measure the energy release rate in pure mode III, which is not yet possible with setups found within literature. For isotropic, brittle or quasi-brittle materials, the results can approximately be converted into the stress intensity factors in mode I and mode III.
  • The observed pop-in fracture in mode III could be crucial for the structural integrity of thin-walled POM-H components. Crack growth prior to reaching the critical energy release rate can significantly weaken a structure when loaded in shear.
  • The analysis of the fracture surfaces showed no signs of plastic deformation close to the initial notch in both modes I and III. The mode I specimens displayed a macroscopically flat fracture surface, whereas the mode III samples showed a deflection of the crack plane. The twisting of the crack path was attributed to a mixed-mode crack propagation and requires further investigation.
Additionally, we suggest that further research should be undertaken in the following areas:
  • Future studies should address the applicability of the presented test setups to other materials, especially more ductile polymers.
  • A simulative study of the experimental design could provide information on the validity of the test setups proposed in this study. Furthermore, an investigation of the influence of thermal clamping and the stress state at the crack tip using finite element analysis is suggested.
  • The influences of specimen thickness, initial crack position, groove geometry, and influences due to notching have not been investigated within this study. We suggest that further studies focus on this, as their influence on the fracture behavior determined from the given test setup is not yet foreseeable.
  • The test setups can be modified and “superimposed” to achieve a mixed mode I/III load. For this, it may be necessary to revise the clamping procedure.
  • The fracture surface investigation should be expanded to obtain more information about the influence of pop-ins and the deflection of the crack-plane in the mode III tests.

Author Contributions

Conceptualization, P.S., A.G. and S.M.; methodology, P.S.; investigation, P.S. and A.G.; formal analysis, P.S.; data curation, P.S.; writing—original draft preparation, P.S. and A.G.; writing—review and editing, P.S.; visualization, P.S.; supervision, M.B. and S.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This article is part of P. Schrader’s doctoral thesis at the Doctoral Center for Engineering Sciences of the Research Campus of Central Hessen under the supervision of the Justus-Liebig-University Gießen in cooperation with the University of Applied Sciences of Central Hessen (Technische Hochschule Mittelhessen). The authors want to thank Bernd Haar for the production of the POM plates.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
POMPolyoxymethylene
DCBDouble cantilever beam
ODCBOut-of-plane-loaded double cantilever beam
ERREnergy release rate

Appendix A. Production of the POM-H Specimens

For a better overview, we want to disclose the manufacturing steps and parameters for the production of the POM-H specimens in greater detail:
  • Step 1: Compression molding tubular granules into plates of approx. 200 × 200 × 5.5 mm; Hydrostat 300 (Schwabenthan, Germany) with an immersion edge tool (TT-260, Tool-Temp)
  • Step 2: Cutting the pressed plates into a rectangular shape of size 160 × 25 mm; Bäuerle tablesaw KSW-7 (Riston Werkzeug GmbH, Fellbach, Germany)
  • Step 3: Milling the rectangles to correct thickness; universal milling machine Deckel FFP3 (Friedrich Deckel AG, München, Germany) with HSS face-milling cutter ø60 mm at 500 rpm
  • Step 4: Milling the slit between the lever arms; universal milling machine Deckel FFP3 (Friedrich Deckel AG, München, Germany) with HSS prismatic disk cutter ø250 mm at 80 rpm
  • Step 5: Cutting the sidegrooves; Bäuerle tablesaw KSW-7 (Riston Werkzeug GmbH, Fellbach, Germany) with prismatic sawblade

References

  1. Berer, M.; Major, Z. Characterization of the global deformation behaviour of engineering plastics rolls. Int. J. Mech. Mater. Des. 2010, 6, 1–9. [Google Scholar] [CrossRef]
  2. Berer, M.; Major, Z. Characterisation of the Local Deformation Behaviour of Engineering Plastics Rolls. Strain 2012, 48, 225–234. [Google Scholar] [CrossRef]
  3. Berer, M.; Mitev, I.; Pinter, G. Finite element study of mode I crack opening effects in compression-loaded cracked cylinders. Eng. Fract. Mech. 2017, 175, 1–14. [Google Scholar] [CrossRef]
  4. Hashemi, S.; Elmes, P.; Sandford, S. Hybrid effects on mechanical properties of polyoxymethylene. Polym. Eng. Sci. 1997, 37, 45–58. [Google Scholar] [CrossRef]
  5. Berer, M.; Pinter, G. Determination of crack growth kinetics in non-reinforced semi-crystalline thermoplastics using the linear elastic fracture mechanics (LEFM) approach. Polym. Test. 2013, 32, 870–879. [Google Scholar] [CrossRef]
  6. Berer, M.; Pinter, G.; Feuchter, M. Fracture mechanical analysis of two commercial polyoxymethylene homopolymer resins. J. Appl. Polym. Sci. 2014, 131, n/a–n/a. [Google Scholar] [CrossRef]
  7. Ramoa, B.; Berer, M.; Schwaiger, M.; Pinter, G. Effect of cyclic fatigue on the fracture toughness of Polyoxymethylene. J. Physics Conf. Ser. 2017, 843, 012052. [Google Scholar] [CrossRef] [Green Version]
  8. Gosch, A.; Berer, M.; Hutař, P.; Slávik, O.; Vojtek, T.; Arbeiter, F.J.; Pinter, G. Mixed Mode I/III fatigue fracture characterization of Polyoxymethylene. Int. J. Fatigue 2020, 130, 105269. [Google Scholar] [CrossRef]
  9. Aliha, M.; Bahmani, A.; Akhondi, S. Determination of mode III fracture toughness for different materials using a new designed test configuration. Mater. Des. 2015, 86, 863–871. [Google Scholar] [CrossRef]
  10. Anthony, J.; Paris, P. Instantaneous evaluation of J and C. Int. J. Fract. 1988, 38, R19–R21. [Google Scholar] [CrossRef]
  11. Olsson, P.; Stigh, U. On the determination of the constitutive properties of thin interphase layers—An exact inverse solution. Int. J. Fract. 1989, 41, R71–R76. [Google Scholar] [CrossRef]
  12. Sørensen, B.F.; Jacobsen, T.K. Determination of cohesive laws by the J integral approach. Eng. Fract. Mech. 2003, 70, 1841–1858. [Google Scholar] [CrossRef]
  13. Gunderson, J.D.; Brueck, J.F.; Paris, A.J. Alternative test method for interlaminar fracture toughness of composites. Int. J. Fract. 2007, 143, 273–276. [Google Scholar] [CrossRef]
  14. Loh, L.; Marzi, S. An Out-of-plane Loaded Double Cantilever Beam (ODCB) test to measure the critical energy release rate in mode III of adhesive joints: Special issue on joint design. Int. J. Adhes. Adhes. 2018, 83, 24–30. [Google Scholar] [CrossRef]
  15. Rice, J.R. A Path Independent Integral and the Approximate Analysis of Strain Concentration by Notches and Cracks. J. Appl. Mech. 1968, 35, 379–386. [Google Scholar] [CrossRef] [Green Version]
  16. Hutchinson, J.W. Singular behaviour at the end of a tensile crack in a hardening material. J. Mech. Phys. Solids 1968, 16, 13–31. [Google Scholar] [CrossRef]
  17. Rice, J.R.; Rosengren, G.F. Plane strain deformation near a crack tip in a power-law hardening material. J. Mech. Phys. Solids 1968, 16, 1–12. [Google Scholar] [CrossRef]
  18. Zehnder, A.T. Fracture Mechanics; Lecture Notes in Applied and Computational Mechanics; Springer Science+Business Media: London, UK; New York, NY, USA, 2012; Volume 62. [Google Scholar]
  19. Loh, L.; Marzi, S. Mixed-mode I+III tests on hyperelastic adhesive joints at prescribed mode-mixity. Int. J. Adhes. Adhes. 2018, 85, 113–122. [Google Scholar] [CrossRef]
  20. Loh, L.; Marzi, S. A Mixed-Mode Controlled DCB test on adhesive joints loaded in a combination of modes I and III. Procedia Struct. Integr. 2018, 13, 1318–1323. [Google Scholar] [CrossRef]
  21. Loh, L.; Marzi, S. A novel experimental methodology to identify fracture envelopes and cohesive laws in mixed-mode I + III. Eng. Fract. Mech. 2019, 214, 304–319. [Google Scholar] [CrossRef]
  22. Tscharnuter, D.; Gastl, S.; Pinter, G. Modeling of the nonlinear viscoelasticity of polyoxymethylene in tension and compression. Int. J. Eng. Sci. 2012, 60, 37–52. [Google Scholar] [CrossRef]
  23. Brandstetter, D. Überprüfung eines Zugkriechprüfstands auf Genauigkeit und Reproduzierbarkeit. Bachelor’s Thesis, Montanuniversität Leoben, Leoben, Austria, 2014. [Google Scholar]
  24. Usman, K. Mechanical Behaviour of POM in Dependence on Molar Mass, Temperature and Loading Regime. Master’s Thesis, Technische Universität Chemnitz, Chemnitz, Germany, 2019. [Google Scholar]
  25. D30 Committee. ASTM D5528–13, Test Method for Mode I Interlaminar Fracture Toughness of Unidirectional Fiber-Reinforced Polymer Matrix Composites; ASTM Internationa: West Conshohocken, PA, USA, 2013. [Google Scholar] [CrossRef]
  26. Berejnoi, C.; Perez Ipiña, J.E.; Llorente, C.L. Reproducibility of pop-ins in laboratory testing of welded joints. Mater. Res. 2000, 3, 139–146. [Google Scholar] [CrossRef] [Green Version]
  27. Pisarski, H.G.; Hammond, R.; Watt, K. Significance of Splits and Pop-Ins Observed During Fracture Toughness Testing of Line Pipe Steel. In Proceedings of the 2008 7th International Pipeline Conference, Calgary, AB, Canada, 29 September–3October 2008; Volume 3, pp. 473–482. [Google Scholar] [CrossRef]
  28. Anderson, T.L. Fracture Mechanics: Fundamentals and Applications, 4th ed.; Chapman and Hall/CRC: Boca Raton, FL, USA, 2017. [Google Scholar]
  29. Pokluda, J.; Pippan, R.; Vojtek, T.; Hohenwarter, A. Near-threshold behaviour of shear-mode fatigue cracks in metallic materials. Fatigue Fract. Eng. Mater. Struct. 2014, 37, 232–254. [Google Scholar] [CrossRef]
  30. Vojtek, T.; Pokluda, J.; Šandera, P.; Horníková, J.; Hohenwarter, A.; Pippan, R. Analysis of fatigue crack propagation under mixed mode II+III in ARMCO iron. Int. J. Fatigue 2015, 76, 47–52. [Google Scholar] [CrossRef]
  31. Pook, L.P. Crack profiles and corner point singularities. Fatigue Fract. Eng. Mater. Struct. 2000, 23, 141–150. [Google Scholar] [CrossRef]
  32. Pook, L.P. Five decades of crack path research. Eng. Fract. Mech. 2010, 77, 1619–1630. [Google Scholar] [CrossRef]
  33. Pook, L.P. The linear elastic analysis of cracked bodies and crack paths. Theor. Appl. Fract. Mech. 2015, 79, 34–50. [Google Scholar] [CrossRef]
  34. Plummer, C.J.G.; Béguelin, P.; Kausch, H.H. The temperature and strain-rate dependence of mechanical properties in polyoxymethylene. Polym. Eng. Sci. 1995, 35, 1300–1312. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of the line J-integral around a notch for a plane problem.
Figure 1. Schematic representation of the line J-integral around a notch for a plane problem.
Materials 13 05096 g001
Figure 2. Schematic drawing of a loaded specimen. Measured forces and moments are displayed in blue color; measured or prescribed rotations are displayed in red.
Figure 2. Schematic drawing of a loaded specimen. Measured forces and moments are displayed in blue color; measured or prescribed rotations are displayed in red.
Materials 13 05096 g002
Figure 3. Geometry of the tested polyoxymethylene homopolymer (POM-H) specimens (in mm). The specimens were notched via pushing with a thin razor blade (notch depth ≈ 1 mm).
Figure 3. Geometry of the tested polyoxymethylene homopolymer (POM-H) specimens (in mm). The specimens were notched via pushing with a thin razor blade (notch depth ≈ 1 mm).
Materials 13 05096 g003
Figure 4. Geometry of the aluminum specimen holder (in mm). The nominal geometrical moments of inertia for the given geometry equate to I y = 4116.66 mm4, I z = 3343.02 mm4, and I t = 7459.68 mm4 around the area’s centroid.
Figure 4. Geometry of the aluminum specimen holder (in mm). The nominal geometrical moments of inertia for the given geometry equate to I y = 4116.66 mm4, I z = 3343.02 mm4, and I t = 7459.68 mm4 around the area’s centroid.
Materials 13 05096 g004
Figure 5. Specimen in mode I investigation (left) and mode III investigation (right). In mode I, the specimen is clamped directly into the test setup whereas in mode III, aluminum specimen holders are added.
Figure 5. Specimen in mode I investigation (left) and mode III investigation (right). In mode I, the specimen is clamped directly into the test setup whereas in mode III, aluminum specimen holders are added.
Materials 13 05096 g005
Figure 6. Experimental setup with POM-H specimens. Mode I setup (left) and mode III setup (right).
Figure 6. Experimental setup with POM-H specimens. Mode I setup (left) and mode III setup (right).
Materials 13 05096 g006
Figure 7. Measured force-displacement curves in mode I (left) and cube root of compliance over crack length with 5% confidence bands (right). Eleven specimens were tested.
Figure 7. Measured force-displacement curves in mode I (left) and cube root of compliance over crack length with 5% confidence bands (right). Eleven specimens were tested.
Materials 13 05096 g007
Figure 8. Influence of the correction factor Δ on the energy release rate (ERR) obtained from modified beam theory.
Figure 8. Influence of the correction factor Δ on the energy release rate (ERR) obtained from modified beam theory.
Materials 13 05096 g008
Figure 9. Measured values of the mode I J-integral (left) and obtained values for J I c and J I , rest over the generated crack surface (right).
Figure 9. Measured values of the mode I J-integral (left) and obtained values for J I c and J I , rest over the generated crack surface (right).
Materials 13 05096 g009
Figure 10. Measured moment M y over the rotational angle α (left) and linear regression between the curvature 2 κ and the measured moment M y (right) obtained with three specimens.
Figure 10. Measured moment M y over the rotational angle α (left) and linear regression between the curvature 2 κ and the measured moment M y (right) obtained with three specimens.
Materials 13 05096 g010
Figure 11. Experimental results of the mode III investigation. The shown mode III ERR of 15 specimens was determined from the Irwin-Kies Equation.
Figure 11. Experimental results of the mode III investigation. The shown mode III ERR of 15 specimens was determined from the Irwin-Kies Equation.
Materials 13 05096 g011
Figure 12. Ratio between the observed pop-in values of J I I I and J I I I c .
Figure 12. Ratio between the observed pop-in values of J I I I and J I I I c .
Materials 13 05096 g012
Figure 13. Artificial contributions determined within the mode III investigation. The unintended out-of-plane mode I load J I * (left) and the mode I/II contribution J I + I I (right) are displayed normalized on J I I I c over J I I I / J I I I c .
Figure 13. Artificial contributions determined within the mode III investigation. The unintended out-of-plane mode I load J I * (left) and the mode I/II contribution J I + I I (right) are displayed normalized on J I I I c over J I I I / J I I I c .
Materials 13 05096 g013
Figure 14. Fracture surface images from a monotonically loaded mode I specimen analyzed using a light microscope (16× magnification)—three areas of different crack areas were observed: pre-notch before testing (1), the crack growth area (2) and the area produced during the separation of the two specimen halves (3).
Figure 14. Fracture surface images from a monotonically loaded mode I specimen analyzed using a light microscope (16× magnification)—three areas of different crack areas were observed: pre-notch before testing (1), the crack growth area (2) and the area produced during the separation of the two specimen halves (3).
Materials 13 05096 g014
Figure 15. Fracture surface images from monotonically loaded mode III specimens (side view)—two different types of surface structure were observed: twisted crack planes along the fracture surface (a) and a nearly straight crack plane with little deflections (b).
Figure 15. Fracture surface images from monotonically loaded mode III specimens (side view)—two different types of surface structure were observed: twisted crack planes along the fracture surface (a) and a nearly straight crack plane with little deflections (b).
Materials 13 05096 g015
Figure 16. Fracture surface images from two mode III specimens analyzed using a light microscope (16× magnification). Image of a deflected fracture surface (a) and a flat fracture surface (b)—two different crack areas were observed on the fractured specimen half: pre-notch before testing (1) and area of crack growth during testing (2).
Figure 16. Fracture surface images from two mode III specimens analyzed using a light microscope (16× magnification). Image of a deflected fracture surface (a) and a flat fracture surface (b)—two different crack areas were observed on the fractured specimen half: pre-notch before testing (1) and area of crack growth during testing (2).
Materials 13 05096 g016
Table 1. Function parameters and correlation coefficients of the linear regressions between J I c and J I , test vs. crack length.
Table 1. Function parameters and correlation coefficients of the linear regressions between J I c and J I , test vs. crack length.
SlopeIntercept R 2
[(kJ/m2)/mm2][kJ/m2][-]
J I c 3.249 × 10−26.120.583
J I , rest −2.286 × 10−25.370.783
Table 2. Summary of the most important experimentally obtained parameters.
Table 2. Summary of the most important experimentally obtained parameters.
Mode IMode III
G c from modified beam theory(6.58 ± 1.09) kJ/m2n.a.
critical value of J(8.84 ± 1.39) kJ/m2(7.59 ± 1.19) kJ/m2
approx. fracture toughness K5.6 MPa m1/2 (plane stress)4.3 MPa m1/2
6.1 MPa m1/2 (plane strain)
min. of J to cause crack propagation4.56–6.19 kJ/m2n.a.
pop-in rationot observed53–93%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Schrader, P.; Gosch, A.; Berer, M.; Marzi, S. Fracture of Thin-Walled Polyoxymethylene Bulk Specimens in Modes I and III. Materials 2020, 13, 5096. https://doi.org/10.3390/ma13225096

AMA Style

Schrader P, Gosch A, Berer M, Marzi S. Fracture of Thin-Walled Polyoxymethylene Bulk Specimens in Modes I and III. Materials. 2020; 13(22):5096. https://doi.org/10.3390/ma13225096

Chicago/Turabian Style

Schrader, Peer, Anja Gosch, Michael Berer, and Stephan Marzi. 2020. "Fracture of Thin-Walled Polyoxymethylene Bulk Specimens in Modes I and III" Materials 13, no. 22: 5096. https://doi.org/10.3390/ma13225096

APA Style

Schrader, P., Gosch, A., Berer, M., & Marzi, S. (2020). Fracture of Thin-Walled Polyoxymethylene Bulk Specimens in Modes I and III. Materials, 13(22), 5096. https://doi.org/10.3390/ma13225096

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop