luchete80 / WeldForm

Weakly-Compressible Smoothed Particle Hydrodynamics Parallel Solver for Elasto Plastic and thermal coupled Solid Mechanics
GNU General Public License v3.0
9 stars 5 forks source link

Add a output ERROR note when intial gap is large since solver diverges #313

Open luchete80 opened 1 year ago

luchete80 commented 1 year ago

If contact compression example is like below works ok:

"RigidBodies": [{ "type": "Plane", "zoneId": 10, "start": [-0.25, -0.25, 0.5975], "partSide": 20, "dim": [0.5, 0.5, 0.0], "translation": [1.0, 0.0, 0.0], "scale": [1, 1, 1] }],

but, if this is changed to:

"start": [-0.25, -0.25, 0.5974], i:e, 1e-4 difference, it diverges with initial stresses huge


Total CPU time: 0.17774 Calculation Times Accel: 0.09%, Density: 0.062%, Stress: 0.1%, Energy: 0.017%, Contact: 0.034%, Nb: -5.5%, Update: 0.0056%, Output No. 2 at 1.6e-06 has been generated Current Time Step = 1.6e-06 Max plastic strain: 0in particle11399504 Max Displacements (No Cont Surf): 3 [ 0 0 0 ] Contact Force Sum 2.35251e+08, Reaction Sum 2.35251e+08 Contact Friction Work 0 External Forces Work 0 Int Energy: 0, Kin Energy: 0


Total CPU time: 0.653614 Calculation Times Accel: 0.17%, Density: 0.12%, Stress: 0.2%, Energy: 0.028%, Contact: 0.014%, Nb: -1.5%, Update: 0.011%, Output No. 3 at 1.1e-05 has been generated Current Time Step = 1.6e-06 Max plastic strain: 0.0055in particle11123 Plastic Work 3.5e+03 Max Displacements (No Cont Surf): 3 [ 0.0002329480.0002329480.000592889 ] Contact Force Sum 0, Reaction Sum 0 Contact Friction Work 0 External Forces Work 0 Int Energy: 7815.29, Kin Energy: 102.691

luchete80 commented 1 year ago

But with initial gap larger, contact works ok:


Total CPU time: 0.176164 Calculation Times Accel: 0.085%, Density: 0.062%, Stress: 0.1%, Energy: 0.017%, Contact: 0.028%, Nb: -5.5%, Update: 0.0057%, Output No. 2 at 1.6e-06 has been generated Current Time Step = 1.6e-06 Max plastic strain: 0in particle11399504 Max Displacements (No Cont Surf): 3 [ 0 0 0 ] Contact Force Sum 0.000742735, Reaction Sum 0.000742735 Contact Friction Work 0 External Forces Work 0 Int Energy: 0, Kin Energy: 0


Total CPU time: 0.613463 Calculation Times Accel: 0.17%, Density: 0.12%, Stress: 0.22%, Energy: 0.029%, Contact: 0.015%, Nb: -1.6%, Update: 0.0065%, Output No. 3 at 1.1e-05 has been generated Current Time Step = 1.6e-06 Max plastic strain: 0in particle11123 Max Displacements (No Cont Surf): 3 [ 2.06778e-062.06778e-068.28172e-06 ] Contact Force Sum 331301, Reaction Sum 219817 Contact Friction Work 0 External Forces Work 11.1035 Int Energy: 2.59655, Kin Energy: 0.219819

luchete80 commented 1 year ago

TODO: Add a ERROR note when intial gap is large