Fix Nan in Blade Stiffened Shell Constitutive Model

smdogroup / tacs

Finite-element library for analysis and adjoint-based gradient evaluation

Apache License 2.0

108 stars 75 forks source link

Closed sean-engelstad closed 7 months ago

sean-engelstad commented 7 months ago

Getting nans in the bucklingEnvelopeFailureSens while running the Mach wing in Funtofem
N1 << N1crit leads to root becomes 0 due to machine precision. Add small tolerance to prevent divide by zero.

sean-engelstad commented 7 months ago

@A-CGray - I made the changes you requested

sean-engelstad commented 7 months ago

Had to test it out again => switched from max(., .) to if statement