fix: correct the effective stress calculation for thermoporomechanics

[frankfeifan]

This PR aims to correct the current GEOS effective stress calculation such that it is consistent with Terzaghi definition. Detailed explanations on the issue and this change can be seen in the file attached.
THMEffStress_GEOS.pdf

• Correct the effective and total stress calculation as well as dTotalStress_dTemperature in PorousSolid for the thermal case
• Add an integrated test to verify the implementation
• Remove the thermal strain part from the avgPlasticStrain calculation. (assumed avgPlasticStrain to be purely mechanical deformation).
• Revise the THM tutorial benchmark problem accordingly (PR#4110)

Another PR(s) needed to fix following items

• Fix the porosity update for both full implicit case and fixed-stress case (see formulation in the last slide of the pdf above). It can be complex in terms of derivative calculations.
• Prepare tutorial examples with these newly added tests

[frankfeifan]

An integrated test is added to verify the implementation: https://github.com/GEOS-DEV/GEOS/blob/deaa1bda781b98215e21a7671b55756f434a5604/inputFiles/thermoPoromechanics/ThermoDruckerPrager_1DCooling_smoke.xml

The previous GEOS total stress calculation cannot match the analytical solution:
displacementPlots.pdf
stressPlots.pdf

With this fix, the analytical solution can be matched:
displacementPlots.pdf
stressPlots.pdf

[jhuang2601]

Great update, thanks!