search

luchete80 / WeldForm

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

Thermal Expansion #51

Open luchete80 opened 3 years ago

luchete80 commented 3 years ago

image

luchete80 commented 3 years ago

Example is not increasing temperature:

Output No. 2 at 0 has been generated Current Time Step = 0.001728 Total CPU time: 0.471879, Neigbour search time: 0.148, Pr Accel Calc time: 0Las Acel Calc Time0.075 Max plastic strain: 0in particle15 Steps count in this interval: 2419020609Total Step count1 Total Nb search time in this interval: 0.148Average Nb search time in this interval: 6.11818e-11 Avg Neighbour Count805 Max T, Min T, Avg temps: 20, 20, 20

Output No. 3 at 0.100224 has been generated Current Time Step = 0.001728 Total CPU time: 11.707, Neigbour search time: 0.148, Pr Accel Calc time: 0.009Las Acel Calc Time2.651 Max plastic strain: 0in particle15 Steps count in this interval: 58Total Step count59 Total Nb search time in this interval: 0Average Nb search time in this interval: 0 Avg Neighbour Count805 Max T, Min T, Avg temps: 0, 1000, 500

Output No. 4 at 0.200448 has been generated Current Time Step = 0.001728 Total CPU time: 22.9655, Neigbour search time: 0.148, Pr Accel Calc time: 0.017Las Acel Calc Time5.245 Max plastic strain: 0in particle15 Steps count in this interval: 58Total Step count117 Total Nb search time in this interval: 0Average Nb search time in this interval: 0 Avg Neighbour Count805 Max T, Min T, Avg temps: 0, 1000, 500

luchete80 commented 3 years ago

Now, Now is working, but values are wrong:

Output No. 2 at 0 has been generated Current Time Step = 0.001728 Total time: 0.485095, Neigbour search time: 0, Accel Calc time: 0.077 Max, Min, Avg temps: 20, 20, 20 Max flux: 266.667 Max Displacements: 0.416667, 0.416667, 0.416667, Max dTdt: 266.667in particle: 2535 Max dTdt: 220.015in particle: 3150 Max dTdt: 179.707in particle: 0 Max dTdt: 151.113in particle: 3150 Max dTdt: 131.453in particle: 3150 Max dTdt: 117.565in particle: 3150 Max dTdt: 107.229in particle: 3150 Max dTdt: 99.104in particle: 3150 Max dTdt: 92.4466in particle: 3150 Max dTdt: 86.8532in particle: 3150 Max dTdt: 82.083in particle: 3360 Max dTdt: 77.9649in particle: 3150 Max dTdt: 74.3675in particle: 3360 Max dTdt: 71.1915in particle: 3150 Max dTdt: 68.3641in particle: 3360 Max dTdt: 65.8295in particle: 3360 Max dTdt: 63.5424in particle: 0 Max dTdt: 61.4652in particle: 0 Max dTdt: 59.5676in particle: 210

luchete80 commented 3 years ago

StartAcceleration(); ----->THIS RESETS STRAIN RATE!!!! PrimaryComputeAcceleration(); LastComputeAcceleration(); //--->>>CALCULATES Strain Rate BUT FROM ZERO! //AND CALCULATES ACCELERATION FROM SIGMA (PREV. IT)

CalcConvHeat(); CalcTempInc(); CalcThermalExpStrainRate(); //Add Thermal expansion Strain Rate Term
//...... BUT SIGMA IS ALREADY CALCULATED! Move(deltat); /// Uses ACCELERATION TO MOVE PARTICLE (calc u && v) //CALCULATE SIGMA FOM STRAIN RATE

Detail of Move()


inline void Particle::Move_Leapfrog(Mat3_t I, double dt)
{
    if (FirstStep) {
        Densitya = Density - dt/2.0*dDensity;
        va = v - dt/2.0*a;
    }
    Densityb = Densitya;
    Densitya += dt*dDensity;
    Density = (Densitya+Densityb)/2.0;
    vb = va;
    va += dt*a;
    v = (va + vb)/2.0;
    x += dt*va;

    Displacement += dt*va;

    Mat2Leapfrog(dt);
    if (FirstStep) FirstStep = false;

}
luchete80 commented 3 years ago

At 0.1seconds in one example, StrainRate is not null but sigma is When calculating Sigma, StrainRate is not null [ 0 0 0.466139 ]

StrainRate(3,3): [ 0.455175 0 0 ] [ 0 0.455175 0 ] [ 0 0 0.455175 ]

StrainRate(3,3): [ 0.445735 0 0 ] [ 0 0.445735 0 ] [ 0 0 0.445735 ]

StrainRate(3,3): [ 0.438239 0 0 ] [ 0 0.438239 0 ] [ 0 0 0.438239 ]

StrainRate(3,3): [ 0.433041 0 0 ] [ 0 0.433041 0 ] [ 0 0 0.433041 ]

StrainRate(3,3): [ 0.429959 0 0 ] [ 0 0.429959 0 ] [ 0 0 0.429959 ]

StrainRate(3,3): [ 0.54564 0 0 ] [ 0 0.54564 0 ] [ 0 0 0.54564 ]

StrainRate(3,3): [ 0.544136 0 0 ] [ 0 0.544136 0 ] [ 0 0 0.544136 ]

StrainRate(3,3): [ 0.540001 0 0 ] [ 0 0.540001 0 ]

luchete80 commented 3 years ago

At the begining:

StrainRate(3,3): [ 0.179446 0 0 ] [ 0 0.179446 0 ] [ 0 0 0.179446 ]

StrainRate(3,3): [ 0.0183256 0 0 ] [ 0 0.0183256 0 ] [ 0 0 0.0183256 ]

StrainRate(3,3): [ 0.00114861 0 0 ] [ 00.00114861 0 ] [ 0 00.00114861 ]

StrainRate(3,3): [ 5.28963e-05 0 0 ] [ 05.28963e-05 0 ] [ 0 05.28963e-05 ]

StrainRate(3,3): [ 9.05906e-07 0 0 ] [ 09.05906e-07 0 ] [ 0 09.05906e-07 ]

StrainRate(3,3): [ 6.34881e-09 0 0 ] [ 06.34881e-09 0 ] [ 0 06.34881e-09 ]

StrainRate(3,3): [ 1.57471e-11 0 0 ]