Dear community, I am using the following code for a research, but I did not understant how read and check the results in Paraview. Can someone help me? The code was made by Hirshikesh, S. Natarajan, R. K. Annabattula, E. Martinez-Paneda. from dolfin import * mesh = Mesh('mesh.xml') # Define Spac…

To make it easier, this is my code: from dolfin import * from mshr import * import matplotlib.pyplot as plt # Mesh vertices=[Point(8.5,0),Point(5.1,-0.9),Point(0,-2), Point(0,-140),Point(22,-140),Point(22,140),Point(0,140),Point(0,2),Point(5.1,0.9)] domain=Polygon(vertices) mesh=generate_mesh(doma…

You have intentionally created a mesh that excludes the starter crack, therefore there are no points that match your Crack(x) still. The following visually and numerically exemplifies this. class CrackClass(SubDomain): def inside(self, x, on_boundary): return abs(x[1])<1e-3 and x[0]<=8.…

Hey Sven, I am using the same code here and I have a question, if I may ask… If I were to change the loading to shear, do I just change these line of code: bctop = DirichletBC(W.sub(1), load, top) … Traction = dot(sigma(unew),n) fx = Traction[0]*ds(1) fname.write(str(assemble(fx)) + “\n”) I …

The correct mesh is that: vertices=[Point(8.5,0),Point(5.1,-0.9),Point(2.9,-0.9),Point(1.5,-3), Point(0,-2),Point(0,-140),Point(22,-140),Point(22,140), Point(0,140),Point(0,2),Point(1.5,3),Point(2.9,0.9),Point(5.1,0.9)] [image] but I can’t rework the Crack(x) to match the test’s standard. It is e…

If you provide the test standard that could help us provide specific advice. In lieu of that I can suggest the following. Consider if crack tip plasticity is relevant to your standard. Potentially refine the mesh along the x-axis and create a elastic-plastic zone around x=8.5 As per ASTM E1820, pl…

Hey @Sven ! Yes I’ll move the content here. My question: How do I modify the code for shear loading? Specific question: I have a working code for tension loading (vertical loading) (same mesh and code as the authors) and I can compute the reaction force by using Traction = sigma*n where n is the …

arianaqyp, It would help if you included the code from the past post for reference. Your previous comment [image] arianaqyp: I am expecting it to be the same range of reaction force, but I can’t seem to figure out how… is to be expected. You are applying a Displacement driven loading cond…

Hey @Sven ! Thank you for your quick reply! I am still confused on how to replicate something like this: [Screen Shot 2022-02-16 at 12.14.14 PM] [Screen Shot 2022-02-16 at 12.14.24 PM] [Screen Shot 2022-02-16 at 12.14.44 PM] This is their tension graph: [Screen Shot 2022-02-16 at 12.14.56 PM] …

[image] arianaqyp: Is this the right way?: What I understand is that I am choosing the x/y-component and loading it in that direction. Correct. You are accessing that dimension of the VectorFunctionSpace and prescribing a fixed displacement. That is fine. For the post’s reference @arianaqyp ’…

Phase Field Fracture

I/O

Sven February 16, 2022, 1:25pm 18

@arianaqyp I see you have already opened a post related to this question : shear loading horizontal loading on a unit mesh. I recommend that you close that post or move this content there to avoid duplication.

Topic		Replies	Views
Code for bend test using phase field fracture variational formulation	16	1044	April 29, 2022
Phasefield Crack doesnt evolve as it should dolfinx	0	104	May 7, 2025
Challenges in reproducing phase-field fracture benchmark: Unphysical crack trajectory and mesh refinement behavior dolfinx mesh , boundary-conditions	11	238	July 1, 2026
Phase field fracture in Dolfinx General	10	1172	March 30, 2026
Phase field model in dolfinx: did not converge	2	1482	December 28, 2021

Phase Field Fracture

Related topics