

February 23, 2022 at 10:26 amjosemishaelSubscriber
I am performing the static crack growth of a 3D object (a thin walled cylinder) with Semi Elliptical crack. I used SMART crack growth method for the static crack growth. I have to perform fatigue crack growth later.
I have few questions related with the postprocessing of my static crack growth simulation results and implementation of number of contours for obtaining the results for fatigue crack growth simulation.
 Does the semielliptical crack grows in both directions during crack growth, i.e. does the minor radius (a) and major radius (c) increases during the crack growth? If both of them increase, how can I obtain the crack extensions in both directions? The crack extension probe provides only one value per crack front node. Is it 'a' or 'c' for the new crack?
 Is there any possibility to visualize the semielliptical crack alone with the mesh around the crack tip? Something like as shown in the reference 'ANSYS Mechanical APDL Technology showcase Example problems' (Rectangular Block Geometry with Semicircular Surface Flaw and Sweep mesh around crack front, under Chapter  Evaluation of MixedMode Stress intensity Factors and Tstress for 3D Surface Flaws)
 During crack growth simulations, sometimes I received the following warning. ' Contour integration for crack 1 shows that some of contours are beyond the FE model and are to be modified. Please carefully evaluate the reported contour integral results. ' My understanding is that after crack growth, few contours may be outside the FE model and not able to considered in the fracture parameter integral calculations. Is there any way to reduce the number of contours as the simulation progresses? For example, during the beginning, I would like to consider 6 contours and as the crack progresses, the number of contours to be reduced into 4 or 5. This will be critical for the fatigue analysis of a thin walled structure.
I appreciate if someone can provide me some insights to resolve the problems.

March 3, 2022 at 6:18 amDavid WeedAnsys Employee
Thank you for your questions; see my responses below.
1) A semielliptical crack will grow in all directions, but keep in mind that every node on the crack front has a local coordinate system, where the:
xaxis is the direction of crack growth
yaxis defines the direction normal to the crack plane
zaxis defines the direction tangent to the crack plane
See the MAPDL Fracture Analysis Guide, 2.1.3. Step 3: Define Crack Information: https://ansyshelp.ansys.com/account/secured?returnurl=/Views/Secured/corp/v221/en/ans_frac/fraccalcprocedure.html
The crackextension direction (e1), crackplane normal (e2) at the tip, and the cracktangent (e3) at the tip together constitute an orthonormal local coordinate system at the crack tip.
A scalar value for the crack extension increment is sufficient to represent crack growth. In the case that you have mixedmode crack growth, in a fatigue scenario, the direction of the crack may change. In that case, if you wish to know the crack direction, you can use the formula in the MAPDL Fracture Analysis Guide, 3.1.3.1. Paris' Law:
As an aside, it may be helpful to see how the crack is setup using APDL commands and subsequently visualize the local coordinate systems for each node along the crack front. In Mechanical, if you open the ds.dat file (located in the Solver Files directory), you can see how the crack is setup using APDL commands:
cint,ctnc,NS_SECRACK_FRONT,_iAssistNode! Define crack tip node component
cint,ncon,3! Define number of contours
cint,normal,13,2! Define crack plane normal
You can read more about these commands and how they define the crack in the MAPDL Fracture Analysis Guide, 2.1.3.1. Define the CrackTip Node Component and CrackPlane Normal (CINT,CTNC): https://ansyshelp.ansys.com/account/secured?returnurl=/Views/Secured/corp/v221/en/ans_frac/fraccalcprocedure.html
To visualize the local coordinate systems, load the ds.dat into MAPDL, select the nodes along the crack front (using the Named Selection, NS_SECRACK_FRONT and the CMSEL command) and issue the command CINT,PLOT,1 (where '1' is the crack ID). See the plot of the nodes and then the local coordinate systems:
2) It appears that in the Tech Demo showcase, a separate crack mesh was created and brought into Mechanical. In contrast, when you use the native semielliptical crack object, you can't visualize the bodies separately in the same way as the graphic shows. However, if you have the mesh view toggled on and highlight the semielliptical crack object, you can see an outline of the contours overlaid onto the crack mesh:
3) During crack growth simulations, sometimes I received the following warning. ' Contour integration for crack 1 shows that some of contours are beyond the FE model and are to be modified. Please carefully evaluate the reported contour integral results. ' My understanding is that after crack growth, few contours may be outside the FE model and not able to considered in the fracture parameter integral calculations. Is there any way to reduce the number of contours as the simulation progresses? For example, during the beginning, I would like to consider 6 contours and as the crack progresses, the number of contours to be reduced into 4 or 5. This will be critical for the fatigue analysis of a thin walled structure.
One approach is to use the CINT,NCON command to adjust the contours during the analysis. However, for crack growth analysis, the suggested method is to use the CGROW,RMCONT,ESIZE command. It allows you to adjust the size of the mesh as the crack approaches the boundary. Essentially, with a smaller element size, the calculated crack extension increment will be smaller, which can help the fracture contours to remain inside of the FE domain. You may consider using a table to adjust the element size using the aforementioned command so that you can better control when the mesh will become more fine, particularly at the FE boundary. For information on that command, I would recommend reading this help section in the MAPDL Fracture Guide, 3.2.1.6. CrackGrowth Dynamic CrackIncrement Control: https://ansyshelp.ansys.com/account/secured?returnurl=/Views/Secured/corp/v221/en/ans_frac/fracSMART.html
I hope that this information is helpful to you.

March 3, 2022 at 10:14 amjosemishaelSubscriberThank you very much for the detailed explanation.

March 3, 2022 at 6:58 pmDavid WeedAnsys Employeeyou're welcome and let me know if you have any follow up questions.

April 3, 2022 at 7:38 pmTomCSubscriberI have a question in this vein. I'm running some models where cracks on the ~10 micron scale grow until failure on the 1000 microns scale.
It's not computationally efficient to run the same mesh required at 10 microns all the way out to 1000 microns. With Premeshed Cracks, I just run the geometry in steps and it works fairly well. With SemiElliptical cracks, I could do the same, but I'm not sure how to match the size of the minor and major radius between models. If I have a 20 um x 100 um crack that I run 10 substeps of SMART crack growth on, it seems inaccurate to say that the new, coarser mesh model will be 20+Crack Extension x 100+Crack Extension.
I'm not sure where to find the information I need to define a new, larger crack (with a correspondingly coarser mesh).

 You must be logged in to reply to this topic.
Earth Rescue – An Ansys Online Series
The climate crisis is here. But so is the human ingenuity to fight it. Earth Rescue reveals what visionary companies are doing today to engineer radical new ideas in the fight against climate change. Click here to watch the first episode.
Ansys Blog
Subscribe to the Ansys Blog to get great new content about the power of simulation delivered right to your email on a weekly basis. With content from Ansys experts, partners and customers you will learn about product development advances, thought leadership and trends and tips to better use Ansys tools. Sign up here.
 Saving & sharing of Working project files in .wbpz format
 An Unknown error occurred during solution. Check the Solver Output…..
 Understanding Force Convergence Solution Output
 Solver Pivot Warning in Beam Element Model
 Colors and Mesh Display
 How to calculate the residual stress on a coating by Vickers indentation?
 whether have the difference between using contact and target bodies
 What is the difference between bonded contact region and fixed joint
 The solver engine was unable to converge on a solution for the nonlinear problem as constrained.
 User manual

2630

2104

1327

1110

461
© 2023 Copyright ANSYS, Inc. All rights reserved.