search

luchete80 / WeldFormGPU

CUDA implementation of Weakly-Compressible Smoothed Particle Hydrodynamics for Elasto Plastic and thermal coupled Mechanics
GNU General Public License v3.0
2 stars 2 forks source link

Mech example not moving #62

Closed luchete80 closed 2 years ago

luchete80 commented 2 years ago

Veloc is not updated in second time step in gpu case (text says accel but it is vel):

Particle 1282Acc: 0, 0, 0 Particle 1283Acc: 0, 0, 0 Particle 1284Acc: 0, 0, 0 Particle 1285Acc: 0, 0, 0 Particle 1286Acc: 0, 0, 0 Particle 1287Acc: 0, 0, 0 Particle 1288Acc: 0, 0, 0 Particle 1289Acc: 0, 0, 0 Particle 1290Acc: 0, 0, 0 Particle 1291Acc: 0, 0, 0 Particle 1292Acc: 0, 0, 0 Particle 1293Acc: 0, 0, 0 Particle 1294Acc: 0, 0, 0 Particle 1295Acc: 0, 0, 0 Particle 1296Acc: 0, 0, 0 Particle 1297Acc: 0, 0, 0 Particle 1298Acc: 0, 0, 0 Particle 1299Acc: 0, 0, 0 Particle 1300Acc: 0, 0, 0 Particle 1301Acc: 0, 0, 0 Particle 1302Acc: 0, 0, 0 Particle 1303Acc: 0, 0, 0 Particle 1304Acc: 0, 0, 0 Particle 1305Acc: 0, 0, 0 Particle 1306Acc: 0, 0, 0 Particle 1307Acc: 0, 0, 0 Particle 1308Acc: 0, 0, 0 Particle 1309Acc: 0, 0, 0 Particle 1310Acc: 0, 0, 0 Particle 1311Acc: 0, 0, 0 Particle 1312Acc: 0, 0, 0 Particle 1313Acc: 0, 0, 0 Particle 1314Acc: 0, 0, 0 Particle 1315Acc: 0, 0, 0 Particle 1316Acc: 0, 0, 0 Particle 1317Acc: 0, 0, 0 Particle 1318Acc: 0, 0, 0 Particle 1319Acc: 0, 0, 0

whereas in CPU version, boundary particles has different acc

rticle 1250, ID 1, Acc 3 [ 0 0 0 ] Particle 1251, ID 1, Acc 3 [ 0 0 0 ] Particle 1252, ID 1, Acc 3 [ 0 0 0 ] Particle 1253, ID 1, Acc 3 [ 0 0 0 ] Particle 1254, ID 1, Acc 3 [ 0 0 0 ] Particle 1255, ID 1, Acc 3 [ 0 0 0 ] Particle 1256, ID 1, Acc 3 [ 0 0 0 ] Particle 1257, ID 1, Acc 3 [ 0 0 0 ] Particle 1258, ID 1, Acc 3 [ 0 0 0 ] Particle 1259, ID 1, Acc 3 [ 0 0 0 ] Particle 1260, ID 3, Acc 3 [ 0 0 -1 ] Particle 1261, ID 3, Acc 3 [ 0 0 -1 ] Particle 1262, ID 3, Acc 3 [ 0 0 -1 ] Particle 1263, ID 3, Acc 3 [ 0 0 -1 ] Particle 1264, ID 3, Acc 3 [ 0 0 -1 ] Particle 1265, ID 3, Acc 3 [ 0 0 -1 ] Particle 1266, ID 3, Acc 3 [ 0 0 -1 ] Particle 1267, ID 3, Acc 3 [ 0 0 -1 ]

luchete80 commented 2 years ago

If artif viscosity are null, accelerations appear at third step

At second time step accel is diff than 0 in non bc particles, in gpu are null:

Particle 1249Acc: 0, 0, 0 Particle 1250Vel: 0, 0, 0 Particle 1250Acc: 0, 0, 0 Particle 1251Vel: 0, 0, 0 Particle 1251Acc: 0, 0, 0 Particle 1252Vel: 0, 0, 0 Particle 1252Acc: 0, 0, 0 Particle 1253Vel: 0, 0, 0 Particle 1253Acc: 0, 0, 0 Particle 1254Vel: 0, 0, 0 Particle 1254Acc: 0, 0, 0 Particle 1255Vel: 0, 0, 0 Particle 1255Acc: 0, 0, 0 Particle 1256Vel: 0, 0, 0 Particle 1256Acc: 0, 0, 0 Particle 1257Vel: 0, 0, 0 Particle 1257Acc: 0, 0, 0 Particle 1258Vel: 0, 0, 0 Particle 1258Acc: 0, 0, 0 Particle 1259Vel: 0, 0, 0 Particle 1259Acc: 0, 0, 0 Particle 1260Vel: 0, 0, -1 Particle 1260Acc: 0, 0, 0 Particle 1261Vel: 0, 0, -1 Particle 1261Acc: 0, 0, 0 Particle 1262Vel: 0, 0, -1 Particle 1262Acc: 0, 0, 0 Particle 1263Vel: 0, 0, -1 Particle 1263Acc: 0, 0, 0 Particle 1264Vel: 0, 0, -1

CPU

Particle 1252, ID 1, Acc 3 [ 1924.17 1924.17 -72492.1 ] Vel3 [ 0.0013849 0.0013849-0.0521753 ] Particle 1253, ID 1, Acc 3 [ 11019.2 7375.34 -56124.7 ] Vel3 [ 0.007930970.00530832-0.0403951 ] Particle 1254, ID 1, Acc 3 [ -7375.34 11019.2 -56124.7 ] Vel3 [ -0.005308320.00793097-0.0403951 ] Particle 1255, ID 1, Acc 3 [ -102.22 11019.2 -63397.4 ] Vel3 [ -7.35714e-050.00793097-0.0456296 ] Particle 1256, ID 1, Acc 3 [ 1.13687e-12 11019.2 -63448.5 ] Vel3 [ 8.18248e-190.00793097-0.0456664 ] Particle 1257, ID 1, Acc 3 [ 2.27374e-13 11019.2 -63448.5 ] Vel3 [ 1.6365e-190.00793097-0.0456664 ] Particle 1258, ID 1, Acc 3 [ 102.22 11019.2 -63397.4 ] Vel3 [ 7.35714e-050.00793097-0.0456296 ] Particle 1259, ID 1, Acc 3 [ 7375.34 11019.2 -56124.7 ] Vel3 [ 0.005308320.00793097-0.0403951 ] Particle 1260, ID 3, Acc 3 [ 0 0 0 ] Vel3 [ 0 0 -1 ] Particle 1261, ID 3, Acc 3 [ 0 0 0 ] Vel3 [ 0 0 -1 ]

luchete80 commented 2 years ago

Problem is at the calculus of strain rate which is used to calculate sigma:

CPU output

Output No. 3 at 1.43948e-06 has been generated Current Time Step = 1.43948e-06 Total CPU time: 0.0429409 , Nb: 0, StAcc: 0, PrAcc: 0.016Ls Acc: 0.015Contact: 0Msh: 0, BC: 0, mv: 0 Max plastic strain: 0in particle0 Steps count in this interval: 1Total Step count2 Total Nb search time in this interval: 0Average Nb search time in this interval: 0 Avg Neighbour Count42 Particle 1250 Sigma (3,3): [ -721795 0-8.02029e-12 ] [ 0 -721795 1592.08 ] [ -8.02029e-12 1592.08-1.46545e+06 ]

Particle 1250 StrRate (3,3): [ 0 0-2.15106e-16 ] [ 0 0 0.0426999 ] [ -2.15106e-16 0.0426999 -19.9451 ]

luchete80 commented 2 years ago

There are not nb set. Only velocity are set

Time, i,j,vab: 2.6525e-311 0 0 0.000000 -1.000000 0.000000 Time, i,j,vab: 2.6525e-311 0 0 0.000000 0.000000 0.000000 Time, i,j,vab: 2.7777e-311 0 0 0.000000 1.000000 0.000000

luchete80 commented 2 years ago

In the acceleration calculus in step 2, strain rate is still null but sigma i and j are not zero. (Writing file is before move and calculate stress and strain)

Max plastic strain: 0in particle0 Steps count in this interval: 1Total Step count2 Total Nb search time in this interval: 0Average Nb search time in this interval: 0 Avg Neighbour Count42 Particle 1250 Sigma (3,3): [ -721795 0-8.02029e-12 ] [ 0 -721795 1592.08 ] [ -8.02029e-12 1592.08-1.46545e+06 ]

Particle 1250 StrRate (3,3): [ 0 0-2.15106e-16 ] [ 0 0 0.0426999 ] [ -2.15106e-16 0.0426999 -19.9451 ]

------------------------ END STEP--------------------------- Time, i,j,vab: 2.8790e-006 1196 1250 0.000000 0.000000 0.000000 Strain Rate (3,3): [ -0 0 -0 ] [ 0 0 0 ] [ -0 0 -0 ]

Particle 1250, Time 2.87896e-06, i, j ,Sigmai 1196; 1250, (3,3): [ -29845.6 0 0 ] [ 0 -29845.6 796.257 ] [ 0 796.257 -60595.3 ] , (3,3): [ -721795 0-8.02029e-12 ] [ 0 -721795 1592.08 ] [ -8.02029e-12 1592.08-1.46545e+06 ]

but in GPU program, they are null

Time 1.4395e-06 Particle 1250 Sigmazz 0.000000 --------------------------- END STEP -------------------------- Time, i,j,vab: 2.8790e-06 1250 1196 0.000000 0.000000 0.000000 Strain Rate 0.000000 0.000000 0.000000 Particle 1250 Time 2.8790e-06, Sigmaizz 0.000000 , Sigmajzz 0.000000 Time, i,j,vab: 2.8790e-06 1250 1189 0.000000 0.000000 0.000000

luchete80 commented 2 years ago

Original

Output No. 3 at 1.43948e-06 has been generated Current Time Step = 1.43948e-06 Total CPU time: 0.0405519 , Nb: 0, StAcc: 0, PrAcc: 0Ls Acc: 0.031Contact: 0Msh: 0, BC: 0, mv: 0 Max plastic strain: 0in particle0 Steps count in this interval: 1Total Step count2 Total Nb search time in this interval: 0Average Nb search time in this interval: 0 Avg Neighbour Count42 Particle 1250 Pressure 969682 Particle 1250 ShearStress (3,3): [ 247886 08.27901e-12 ] [ 0 247886 1592.08 ] [ 8.27901e-12 1592.08 -495772 ]

Particle 1250 Sigma (3,3): [ -721795 08.27901e-12 ] [ 0 -721795 1592.08 ] [ 8.27901e-12 1592.08-1.46545e+06 ]

Particle 1250 StrRate (3,3): [ 0 02.22045e-16 ] [ 0 0 0.0426999 ] [ 2.22045e-16 0.0426999 -19.9451 ]

GPU

Strain Rate 0.000000 0.000000 0.000000 Particle 1250 Time 1.4395e-06, Sigmaizz 0.000000 , Sigmajzz 0.000000 Time 1.4395e-06 Particle 1250, presure 0.000000 , ShearStresszz 0.000000 Sigmazz 0.000000 --------------------------- END STEP -------------------------- Time, i,j,vab: 2.8790e-06 1250 1196 0.000000 0.000000 0.000000

luchete80 commented 2 years ago

Cs is not initialized

Particle 1250 Time 1.4395e-06, Sigmaizz 0.000000 , Sigmajzz 0.000000 PresEq[i], Cs[i], P0[i],Density[i], RefDensity[i]: 0.000000 0.000000 0.000000 1000.038759 0.000000 Time 1.4395e-06 Particle 1250, presure 0.000000 , ShearStresszz 0.000000 Sigmazz 0.000000

luchete80 commented 2 years ago

Shear stress is null due to null Strain rate in previous step (in particle 1250):

Correct result is (without Arif Visc)

Particle 1250 ShearStress (3,3): [ 247886 0-8.02029e-12 ] [ 0 247886 1592.08 ] [ -8.02029e-12 1592.08 -495772 ]

Particle 1250 Sigma (3,3): [ -721795 0-8.02029e-12 ] [ 0 -721795 1592.08 ] [ -8.02029e-12 1592.08-1.46545e+06 ]

Particle 1250 StrRate (3,3): [ 0 0-2.15106e-16 ] [ 0 0 0.0426999 ] [ -2.15106e-16 0.0426999 -19.9451 ]

------------------------ END STEP--------------------------- Time, i,j,vab: 2.8790e-006 1192 1250 0.000000 0.000000 0.000000

luchete80 commented 2 years ago

Corrected an incrementation error 2bda0d9

luchete80 commented 2 years ago

Updated calculations

CPU

Particle 1250 Pressure 969682 Particle 1250 ShearStress (3,3): [ 247886 0-8.02029e-12 ] [ 0 247886 1592.08 ] [ -8.02029e-12 1592.08 -495772 ]

Particle 1250 Sigma (3,3): [ -721795 0-8.02029e-12 ] [ 0 -721795 1592.08 ] [ -8.02029e-12 1592.08-1.46545e+06 ]

Particle 1250 StrRate (3,3): [ 0 0-2.15106e-16 ] [ 0 0 0.0426999 ] [ -2.15106e-16 0.0426999 -19.9451 ]

------------------------ END STEP--------------------------- Time, i,j,vab: 2.8790e-006 1192 1250 0.000000 0.000000 0.000000 Strain Rate (3,3):

GPU

Particle 1250 Time 1.4395e-06, Sigmaizz 0.000000 , Sigmajzz 0.000000 Time, i,j,vab: 1.4395e-06 1250 1192 0.000000 0.000000 0.000000 Strain Rate i 0.000000 0.000000 0.000000 Particle 1250 Time 1.4395e-06, Sigmaizz 0.000000 , Sigmajzz 0.000000 TOTAL (SUM) Strain Rate part 1250 0.000000 0.000000 0.000000 PresEq[i], Cs[i], P0[i],Density[i], RefDensity[i]: 0.000000 5001.814486 0.000000 2700.038759 2700.000000 Time 1.4395e-06 Particle 1250, presure 969681.519080 , ShearStresszz 0.000000 Sigmazz 0.000000

luchete80 commented 2 years ago

Corrected some mult error 3d587bb

luchete80 commented 2 years ago

Still not adding to Strain rate

Strain Rate 1316 1250 -0.000000 0.000000 -25426.539424 Time, i,j,vab, xij, GK: 1.4395e-06 1250 1233 0.000000 0.000000 0.000000 0.030000 0.060000 0.000000 -39542.145757 Time, i,j,vab, xij, GK: 1.4395e-06 1198 1250 0.000000 0.000000 0.000000 0.030000 0.030000 -0.030000 -847508.127410 Time, i,j,vab, xij, GK: 1.4395e-06 1310 1250 0.000000 0.000000 -1.000000 0.000000 0.000000 0.029999 -5922770.102976 StrainRate_zz tempt 1316 1250 -25426.539424 0.000000