CV312 Lab 5
CV312 Lab 5
Aim
The aim of this software lab is to determine the stress intensity factors and analyze the plate with
the central crack.
Procedure
In this software lab the Ansys software has been used and the first step to create the file and name
as “CentralCrack”. Followed by defining the units, element types and the material properties. The
5 key points were created and those point are shown below;
Key point X Y
number
1 0 0
2 0.02 0
3 0.1 0
4 0.1 0.1
5 0 0.1
Then the line segments were drawn by the appropriate key points to create a shape shown below;
Figure 1
MM/CV 312 Solid Mechanics 2 – Lab 5
Moving on, the line 5 had the spacing ratio of 0.2 and with it the discretizing of the lines 3,4,5 was
done. Then at the crack tip the concentration key point was created. The first-row radius was
0.0025 for the key point 2 and with regards to key point the number of elements around
circumference was 8 and for the midside node position ‘skewed 1/4pt’ was selected. Moreover,
the lines with the boundary conditions and load of -100e6 was applied to create the area. The
meshing and solving of the model were done. Then by using the nodes created through meshing
the crack-face path was defined and named as K1. The 3 nodes and activated was used to define
the local crack-tip coordinate system. Lastly, the KCALC was used to determine the Mode-I stress
intensity Factor.
Figure 2
Loading
The applied stress has created a pressure on the loading of the plates whereby the load pressure
value is -100e6.
Figure 3
MM/CV 312 Solid Mechanics 2 – Lab 5
Freedom case
As the displacement was symmetrical with the area the boundary conditions were defined as
structurally.
Figure 4
MM/CV 312 Solid Mechanics 2 – Lab 5
Figure 5
The theoretical value for the stress intensity factor shown below.
𝑎
𝜂=
𝑏
0.02
=
0.2
= 𝟎. 𝟏
1
𝐶 = (1 − 0.1𝜂2 + 0.96𝜂4 )√
cos 𝜋𝜂
1
= (1 − 0.1(0.1)2 + 0.96(0.1)4 )√
cos 𝜋(0.1)
= 𝟏. 𝟎𝟐𝟒𝟒𝟖
𝐾𝑝𝑖𝑘𝑒𝑙𝑦 = 𝐶𝜎√𝜋𝑎
= 𝟐𝟓. 𝟔𝟖𝟎𝑴𝑷𝒂√𝒎
MM/CV 312 Solid Mechanics 2 – Lab 5
Hence, the above calculation was done using the W. D Pikely equation. The percentage of accuracy
of the software analysis is shown below.
𝐾𝑝𝑖𝑘𝑒𝑙𝑦 − 𝐾𝑎𝑛𝑠𝑦𝑠
𝜖= × 100
𝐾𝑝𝑖𝑘𝑒𝑙𝑦
25.680 − 26.6533
= × 100
25.680
= 𝟑. 𝟕%
When reference to the model and the final result the deformed length is decreased as it reaches the
end of the crack. The defining of the material can be the cause of this cracks for example its elastic
properties. The ability of a material to resist the propagation of cracks can be represented fracture
toughness [3]. The stress intensity factor describes the stress state at a crack tip, is related to the
rate of crack growth, and is used to establish failure criteria due to fracture [2]. It is one of the most
fundamental and useful parameters in all of fracture mechanics [1]. The figure below shows the
crack reformed under the tensile stress applied where the crack is uplifted by a crack.
Figure 6
MM/CV 312 Solid Mechanics 2 – Lab 5
Conclusion
All in all, the stress intensity factor of the cracked plate was determined by the ANSYS software.
Where the lab manual and the model that was created in ANSYS had different node numbers which
is shown in figure 4. To clarify this error the results were compared theoretically by using W. D
Pikely equation and the percentage of accuracy of the software analysis was 3.7%.
References
[1] B. McGinty, "Stress Intensity Factor", Fracturemechanics.org, 2021. [Online]. Available:
https://www.fracturemechanics.org/sif.html. [Accessed: 29- May- 2021].
[2] R. Ribeiro and M. Hill, "A benchmark fracture mechanics solution for a two-dimensional
eigenstrain problem considering residual stress, the stress intensity factor, and
superposition", Engineering Fracture Mechanics, vol. 163, pp. 313-326, 2016. Available:
10.1016/j.engfracmech.2016.06.007.
[3] Z. Zhang, "More accurate stress intensity factor derived by finite element analysis for the ISRM
suggested rock fracture toughness specimen—CCNBD", ResearchGate, 2021. [Online].
Available:
https://www.researchgate.net/publication/240362736_More_accurate_stress_intensity_factor_der
ived_by_finite_element_analysis_for_the_ISRM_suggested_rock_fracture_toughness_specimen-
CCNBD. [Accessed: 29- May- 2021].