SEARCH · Search NASA

Results for “virtual crack closure technique”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2021

An Enhanced Virtual Crack Closure Technique for Stress Intensity Factor Calculation along Arbitrary Crack Fronts and the Application in Hydraulic Fracturing Simulation

The virtual crack closure technique (VCCT) is widely used for calculating energy release rates along crack fronts and modeling the propagation of cracks in solid materials. Although the VCCT formulation for smooth crack fronts has been sufficiently addressed in the literature, the application of VCCT to a nonsmoothed crack front with sharp corners warrants further investigation. In this study, we describe an enhanced VCCT to calculate energy release rates and stress intensity factors for cracks with arbitrary shapes in 3D domains discretized on structured grids. The formulations of the enhanced VCCT were developed and implemented into a multiphysics simulation environment capable of simulating crack propagation in the framework of linear elastic fracture mechanics. Comparisons with existing analytical/numerical solutions and other VCCT approaches were performed to verify the enhanced VCCT in terms of SIF calculation along nonsmoothed crack fronts. We then applied the enhanced VCCT to a hydraulically driven penny-shaped fracture problem to further demonstrate its capability to simulate nonsmoothed fracture propagation.

02 PETROLEUM↗

NASA NTRS2019

A New Approach to Model Delamination Growth in Fatigue Using the Virtual Crack Closure Technique Without Re-Meshing

A modeling approach is proposed to simulate delamination propagation in fatigue that combines the Virtual Crack Closure Technique (VCCT) with a progressive nodal release strategy. The progressive nodal release alleviates the artificial stress concentrations found when using the VCCT with instantaneous release to model 3D delamination without re-meshing. This enables crack shapes that do not conform to the underlying mesh to be readily simulated. The progressive release is formulated such that no artificial growth or healing occurs as a result of a change in loading conditions. Furthermore, expressions for calculating energy release rates at partially released nodes are proposed and applied to arbitrarily shaped delamination fronts. The shape of the delamination front is considered explicitly in the energy release rate calculations, and is determined by considering the normalized crack positions associated with the partial nodal release. Contrary to previous implementations, the local delamination propagation direction is not assumed to follow mesh lines, but instead is computed as part of the iterative procedure proposed to determine the maximum energy release along the crack front. The results obtained suggest the approach can accurately simulate Mode I, Mixed Mode I/II and Mode II fatigue delamination growth and be used to simulate growth of arbitrarily shaped cracks that do not conform to the underlying mesh.

N V De Carvalho↗

NASA NTRS2006

Influence of Finite Element Software on Energy Release Rates Computed Using the Virtual Crack Closure Technique

Strain energy release rates were computed along straight delamination fronts of Double Cantilever Beam, End-Notched Flexure and Single Leg Bending specimens using the Virtual Crack Closure Technique (VCCT). Th e results were based on finite element analyses using ABAQUS# and ANSYS# and were calculated from the finite element results using the same post-processing routine to assure a consistent procedure. Mixed-mode strain energy release rates obtained from post-processing finite elem ent results were in good agreement for all element types used and all specimens modeled. Compared to previous studies, the models made of s olid twenty-node hexahedral elements and solid eight-node incompatible mode elements yielded excellent results. For both codes, models made of standard brick elements and elements with reduced integration did not correctly capture the distribution of the energy release rate acr oss the width of the specimens for the models chosen. The results suggested that element types with similar formulation yield matching results independent of the finite element software used. For comparison, m ixed-mode strain energy release rates were also calculated within ABAQUS#/Standard using the VCCT for ABAQUS# add on. For all specimens mod eled, mixed-mode strain energy release rates obtained from ABAQUS# finite element results using post-processing were almost identical to re sults calculated using the VCCT for ABAQUS# add on.

Krueger, Ronald↗

NASA NTRS2018

Combining Progressive Nodal Release with the Virtual Crack Closure Technique to Model Fatigue Delamination Growth Without Re-Meshing

The present work summarizes an approach to model mixed-mode 3D fatigue crack growth using the Virtual Crack Closure Technique (VCCT) without requiring re-meshing. It is demonstrated that the proposed approach can be used to simulate crack shapes that do not conform to the underlying mesh. The proposed approach relies solely on Paris Law characterization data to model delamination growth. Growth is determined as a post-processing step at the end of each increment, and hence no convergence issues associated with the progressive nodal release are encountered. This approach can be readily applied using standard solid element formulations and is implemented via an interface user element in Abaqus/Standard.

Carvalho, N. V. De↗

NASA NTRS2013

Floating Node Method and Virtual Crack Closure Technique for Modeling Matrix Cracking-Delamination Interaction

A novel approach is proposed for high-fidelity modeling of progressive damage and failure in composite materials that combines the Floating Node Method (FNM) and the Virtual Crack Closure Technique (VCCT) to represent multiple interacting failure mechanisms in a mesh-independent fashion. In this study, the approach is applied to the modeling of delamination migration in cross-ply tape laminates. Delamination, matrix cracking, and migration are all modeled using fracture mechanics based failure and migration criteria. The methodology proposed shows very good qualitative and quantitative agreement with experiments.

DeCarvalho, N. V.↗

NASA NTRS2013

Floating Node Method and Virtual Crack Closure Technique for Modeling Matrix Cracking-Delamination Migration

DeCarvalho, Nelson V.↗

NASA NTRS2016

Modeling Fatigue Damage Onset and Progression in Composites Using an Element-Based Virtual Crack Closure Technique Combined With the Floating Node Method

A new methodology is proposed to model the onset and propagation of matrix cracks and delaminations in carbon-epoxy composites subject to fatigue loading. An extended interface element, based on the Floating Node Method, is developed to represent delaminations and matrix cracks explicitly in a mesh independent fashion. Crack propagation is determined using an element-based Virtual Crack Closure Technique approach to determine mixed-mode energy release rates, and the Paris-Law relationship to obtain crack growth rate. Crack onset is determined using a stressbased onset criterion coupled with a stress vs. cycle curve and Palmgren-Miner rule to account for fatigue damage accumulation. The approach is implemented in Abaqus/Standard® via the user subroutine functionality. Verification exercises are performed to assess the accuracy and correct implementation of the approach. Finally, it was demonstrated that this approach captured the differences in failure morphology in fatigue for two laminates of identical stiffness, but with layups containing θdeg plies that were either stacked in a single group, or distributed through the laminate thickness.

De Carvalho, Nelson V.↗

NASA NTRS1987

Strain-energy release rate analysis of a laminate with a postbuckled delamination

The objectives are to present the derivation of the new virtual crack closure technique, evaluate the accuracy of the technique, and finally to present the results of a limited parametric study of laminates with a postbuckled delamination. Although the new virtual crack closure technique is general, only homogeneous, isotropic laminates were analyzed. This was to eliminate the variation of flexural stiffness with orientation, which occurs even for quasi-isotropic laminates. This made it easier to identify the effect of geometrical parameters on G. The new virtual crack closure technique is derived. Then the specimen configurations are described. Next, the stress analyses is discussed. Finally, the virtual crack closure technique is evaluated and then used to calculate the distribution of G along the delamination front of several laminates with a postbuckled delamination.

Whitcomb, John D.↗

NASA NTRS1989

Strain-energy release rate analysis of plates with postbuckled delaminations

Whitcomb, John D.↗

NASA NTRS1988

A virtual crack-closure technique for calculating stress intensity factors for cracked three dimensional bodies

A three-dimensional virtual crack-closure technique is presented which calculates the strain energy release rates and the stress intensity factors using only nodal forces and displacements from a standard finite element analysis. The technique is an extension of the Rybicki-Kanninen (1977) method, and it assumes that any continuous function can be approximated by a finite number of straight line segments. Results obtained by the method for surface cracked plates with and without notches agree favorably with previous results.

Shivakumar, K. N.↗

NASA NTRSDate not supplied

Development of a Combined Cohesive and Virtual Crack-Closure Technique Approach to Represent R-Curves

Resistance curve (R-curve) effects due to fiber bridging, crack migration, and other blunting mechanisms are common in composite laminates. These mechanisms can dramatically increase the resistance to crack propagation but can be challenging to simulate. Delamination with R-curves can be analyzed using cohesive zone modeling (CZM) or the virtual crack closure technique (VCCT). Large fracture process zones can be simulated with CZM, but they require highly refined meshes. Coarser meshes can be used with VCCT, but this method is only applicable to small fracture process zones with R-curves defined as functions of position. Therefore, a technique with the computational efficiency of VCCT and the natural ability of cohesive elements to represent large fracture process zones is desirable. An approach is proposed that starts as a CZM in which the cohesive traction separation law (TSL) is separated into high-strength (HS) and low-strength (LS) components. The HS part, which is responsible for the mesh requirements of the CZM analysis, is replaced by VCCT. The combined method is evaluated by analyzing the response of double cantilevered beam specimens. Two specimen layup configurations, [0/90/90/0] 3s and [0/90/90/0] 9s , are evaluated. The R-curve response from each configuration is determined and applied to the constitutive properties of models built with CZM and with VCCT. The results and computational efficiency of the CZM, VCCT, and combined approaches are compared. The results are indicative that a combined cohesive/VCCT approach can enable progressive failure analyses to retain the computational efficiency of VCCT with the ability of the cohesive elements to capture R-curve effects.

Virtual Crack-Closure Technique↗

NASA NTRS1993

Analysis of local delaminations caused by angle ply matrix cracks

Two different families of graphite/epoxy laminates with similar layups but different stacking sequences, (0,theta,-theta) sub s and (-theta/theta/0) sub s were analyzed using three-dimensional finite element analysis for theta = 15 and 30 degrees. Delaminations were modeled in the -theta/theta interface, bounded by a matrix crack and the stress free edge. The total strain energy release rate, G, along the delamination front was computed using three different techniques: the virtual crack closure technique (VCCT), the equivalent domain Integral (EDI) technique, and a global energy balance technique. The opening fracture mode component of the strain energy release rate, Gl, along the delamination front was also computed for various delamination lengths using VCCT. The effect of residual thermal and moisture stresses on G was evaluated.

Salpekar, Satish A.↗

NASA NTRS2016

Mesh Convergence Requirements for Composite Damage Models

The ability of the finite element method to accurately represent the response of objects with intricate geometry and loading renders the finite element method as an extremely versatile analysis technique for structural analysis. Finite element analysis is routinely used in industry to calculate deflections, stress concentrations, natural frequencies, buckling loads, and much more. The method works by discretizing complex problems into smaller, simpler approximations that are valid over small uniform domains. For common analyses, the maximum size of the elements that can be used is often be determined by experience. However, to verify the quality of a solution, analyses with several levels of mesh refinement should be performed to ensure that the solution has converged. In recent years, the finite element method has been used to calculate the resistance of structures, and in particular that of composite structures. A number of techniques such as cohesive zone modeling, the virtual crack closure technique, and continuum damage modeling have emerged that can be used to predict cracking, delaminations, fiber failure, and other composite damage modes that lead to structural collapse. However, damage models present mesh refinement requirements that are not well understood. In this presentation, we examine different mesh refinement issues related to the representation of damage in composite materials. Damage process zone sizes and their corresponding mesh requirements will be discussed. The difficulties of modeling discontinuities and the associated need for regularization techniques will be illustrated, and some unexpected element size constraints will be presented. Finally, some of the difficulties in constructing models of composite structures capable of predicting transverse matrix cracking will be discussed. It will be shown that to predict the initiation and propagation of transverse matrix cracks, their density, and their saturation may require models that are significantly more refined than those that have been contemplated in the past.

Davila, Carlos G.↗

NASA NTRS2020

DEVELOPMENT AND APPLICATION OF BENCHMARK EXAMPLES FOR MIXED-MODE I/II FATIGUE GROWTH PREDICTIONS

The development of benchmark examples for assessing growth prediction capabilities is presented and demonstrated for a commercial code. The examples are based on finite element models of Mixed-Mode Bending (MMB) specimens for three different mode ratios. The examples are independent of the analysis software used and allow the assessment of the automated delamination propagation, onset and growth prediction capabilities in commercial finite element codes based on the virtual crack closure technique (VCCT). First, based on existing results for delamination propagation under static loading, new benchmark examples for delamination growth under fatigue loading were created. Second, starting from an initially straight front, the delamination was allowed to grow under fatigue loading. The number of cycles to delamination onset and the number of cycles during delamination growth for each growth increment were obtained from the automated analysis and compared to the benchmark examples. Good agreement between the results obtained from the growth analysis and the benchmark results could be achieved by selecting the appropriate input parameters. The benchmarking procedure proved valuable by highlighting the issues associated with choosing the input parameters of the particular implementation. Overall the results from automatic propagation are encouraging, however, selecting the appropriate input parameters was not straightforward and often required an iterative procedure.

Composites↗

NASA NTRSDate not supplied

Recent Advances in Discrete Crack Modeling Applied to Laminated Composites with Emphasis on: Floating Node Method, VCCT and Cohesive Zone Modeling

The present talk will provide an overview of the work performed during the Advanced Composites Project (ACP) on the development, and verification and validation of the Floating Node Method (FNM) as well as the Virtual Crack Closure Technique (VCCT) and cohesive zone modeling (CZM). The FNM is a finite element-based technique to represent crack networks. The complex nature of matrix crack-delamination interactions observed in unidirectional (UD) tape laminates suggests that a methodology such as the FNM method may be required to be able to accurately simulate damage progression in these laminates. Simulating crack onset and growth within the context of the FNM relies on techniques such as VCCT and CZM. The talk is organized to provide, via select examples, an overview of the breadth of the Verification & Validation (V&V) exercises performed during the ACP, and how these challenged the state-of-the art and guided further developments in discrete crack modeling, while helping to establish confidence in the progress made and map the challenges ahead. The performance of the VCCT and CZM individually, and in combination with the FNM can be assessed through verification exercises. These exercises typically consist of a comparison of simulation results to known numerical or analytical solutions. Verification is key to identify implementation issues and limitations that, otherwise, may remain undetected and cloud any subsequent validation efforts. Indeed, a subset of these numerical exercises led to further developments of the VCCT and the FNM method as will be illustrated. Before embarking on the subsequent validation of the framework, it is critical to have adequate characterization data. However, the testing campaign conducted revealed material responses that challenged the state-of-the-art and required further developments. The developments in CZM technology associated with the modeling of the responses of hybrid interfaces (fabric/UD) will be given as an example. Finally, the talk will conclude with a summary of the validation exercises performed under quasi-static and fatigue loadings, highlighting some of the key achievements, outstanding challenges and lessons learned.

finite elements↗

NASA NTRS2002

Analysis of Composite Skin-Stiffener Debond Specimens Using Volume Elements and a Shell/3D Modeling Technique

The debonding of a skin/stringer specimen subjected to tension was studied using three-dimensional volume element modeling and computational fracture mechanics. Mixed mode strain energy release rates were calculated from finite element results using the virtual crack closure technique. The simulations revealed an increase in total energy release rate in the immediate vicinity of the free edges of the specimen. Correlation of the computed mixed-mode strain energy release rates along the delamination front contour with a two-dimensional mixed-mode interlaminar fracture criterion suggested that in spite of peak total energy release rates at the free edge the delamination would not advance at the edges first. The qualitative prediction of the shape of the delamination front was confirmed by X-ray photographs of a specimen taken during testing. The good correlation between prediction based on analysis and experiment demonstrated the efficiency of a mixed-mode failure analysis for the investigation of skin/stiffener separation due to delamination in the adherents. The application of a shell/3D modeling technique for the simulation of skin/stringer debond in a specimen subjected to three-point bending is also demonstrated. The global structure was modeled with shell elements. A local three-dimensional model, extending to about three specimen thicknesses on either side of the delamination front was used to capture the details of the damaged section. Computed total strain energy release rates and mixed-mode ratios obtained from shell/3D simulations were in good agreement with results obtained from full solid models. The good correlations of the results demonstrated the effectiveness of the shell/3D modeling technique for the investigation of skin/stiffener separation due to delamination in the adherents.

Krueger, Ronald↗

NASA NTRS2001

A Shell/3D Modeling Technique for the Analyses of Delaminated Composite Laminates

A shell/3D modeling technique was developed for which a local three-dimensional solid finite element model is used only in the immediate vicinity of the delamination front. The goal was to combine the accuracy of the full three-dimensional solution with the computational efficiency of a plate or shell finite element model. Multi-point constraints provided a kinematically compatible interface between the local three-dimensional model and the global structural model which has been meshed with plate or shell finite elements. Double Cantilever Beam (DCB), End Notched Flexure (ENF), and Single Leg Bending (SLB) specimens were modeled using the shell/3D technique to study the feasibility for pure mode I (DCB), mode II (ENF) and mixed mode I/II (SLB) cases. Mixed mode strain energy release rate distributions were computed across the width of the specimens using the virtual crack closure technique. Specimens with a unidirectional layup and with a multidirectional layup where the delamination is located between two non-zero degree plies were simulated. For a local three-dimensional model, extending to a minimum of about three specimen thicknesses on either side of the delamination front, the results were in good agreement with mixed mode strain energy release rates obtained from computations where the entire specimen had been modeled with solid elements. For large built-up composite structures modeled with plate elements, the shell/3D modeling technique offers a great potential for reducing the model size, since only a relatively small section in the vicinity of the delamination front needs to be modeled with solid elements.

Krueger, Ronald↗

NASA NTRS1994

Separation of crack extension modes in composite delamination problems

This work concerns fracture mechanics modeling of composite delamination problems. In order to predict delamination resistance, an applied stress intensity factor, K, or energy release rate, G, must be compared to a mode-dependent critical value of K or G from experiment. In the interfacial fracture analysis of most applications and some tests, the mode of crack extension is not uniquely defined. It is instead a function of distance from the crack tip due to the oscillating singularity existing at the tip. In this work, a consistent method is presented of extracting crack extension modes in such cases. In particular, use of the virtual crack closure technique (VCCT) to extract modes of crack extension is studied for cases of a crack along the interface between two in-plane orthotropic materials. Modes of crack extension extracted from oscillatory analyses using VCCT are a function of the virtual crack extension length, delta. Most existing efforts to obtain delta-independent modes of crack extension involve changing the analysis in order to eliminate its oscillatory nature. One such method involves changing one or more properties of the layers to make the oscillatory exponent parameter, epsilon, equal zero. Standardized application of this method would require consistent criteria for identifying which properties can be altered without changing the physical aspects of the problem. Another method involves inserting a thin homogeneous layer (typically referred to as a resin interlayer) along the interface and placing the crack within it. The drawbacks of this method are that it requires increased modeling effort and introduces the thickness of the interlayer as an additional length parameter. The approach presented here does not attempt to alter the interfacial fracture analysis to eliminate its oscillatory behavior. Instead, the argument is made that the oscillatory behavior is non-physical and that if its effects were separated from VCCT quantities, then consistent, delta-independent modes of crack extension could be defined. Knowledge of the near-tip fields in a planar orthotropic material interfacial fracture analysis is used to determine the explicit delta dependence of VCCT parameters. Once this delta dependence is determined, energy release rates are defined with this delta dependence factored out. This modified VCCT method is applied to results from two finite element test cases. It is shown that, as predicted, delta-independent modes of crack extension result. The modified VCCT approach shows potential as a consistent method of extracting crack extension modes. It uses the same information from a finite element analysis (i.e., nodal forces and displacements) as the traditional VCCT method does. The A-independent modes extracted using the modified VCCT approach can also be used as guides to test the convergence of finite element solutions.

Beuth, Jack L.↗