Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Bugfix/wrong stress stuff #43

Open
wants to merge 18 commits into
base: master
Choose a base branch
from
Open

Conversation

EleniSmyrniou
Copy link
Collaborator

@EleniSmyrniou EleniSmyrniou commented Feb 17, 2025

Added the following for the norm calc:
Tests for n_method=True and n_method=False:
In the first test, we set n_method=True so that the stress exponent is forced to 0.5. We then compute the expected normalized tip resistance (Qtn) and friction ratio (Fr) using the same formulas and verify that the output arrays match our expectations (taking into account the clamping).

In the second test, we let the iterative procedure compute n (which for these inputs converges to 1.0). The expected values for Qtn and Fr are then computed using the resulting n.

Added the following unit tests for the effective stresses:

test_stress_calc_no_pwp_effect

  • Scenario: The computed pore water pressure (pwp) is zero.
  • Expected Outcome: Effective stress equals total stress.

test_stress_calc_positive_pwp

  • Scenario: A positive pwp is computed that subtracts from the total stress.
  • Expected Outcome: Effective stress is reduced by the pwp amount.

test_stress_calc_negative_pwp_clipped

  • Scenario: The calculated pwp would be negative (indicating suction).
  • Expected Outcome: Negative pwp values are clipped to zero, so effective stress remains equal to total stress.

test_stress_calc_negative_effective_clipped

  • Scenario: The pore water pressure is so high that it would make the effective stress negative.
  • Expected Outcome: Negative effective stress values are clipped to zero.

Copy link

@bramvandeneijnden bramvandeneijnden left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

correct changes to pore water pressure and calculation of Qtn and Fr

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants