Skip to content

[StructuralMechanicsApplication] added calculation of external and internal forces in timoshenko beams - #14837

Merged
aronnoordam merged 4 commits into
masterfrom
struct/external_internal_force_beams
Oct 2, 2026
Merged

aronnoordam merged 4 commits into
masterfrom
struct/external_internal_force_beams

Conversation

@aronnoordam

Copy link
Copy Markdown
Member

added calculation of external and internal forces in timoshenko beams, which can be usefull in specific solving strategies and convergence criteria

@aronnoordam
aronnoordam marked this pull request as ready for review October 1, 2026 14:14
@aronnoordam
aronnoordam requested a review from a team as a code owner October 1, 2026 14:14
BoundedMatrix<double, 3, 3> T;
// From global to local, need to transpose the matrix
noalias(T) = trans(GetConsistentFrenetSerretMatrix3D(r_geom));
// From global to local (the rows of T are the local axes)

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

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

Hi, I solved an issue related to this not long ago. Did I do something wrong? is this a new correction?

@AlejandroCornejo AlejandroCornejo Oct 1, 2026 •

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

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

Maybe I forgot to also correct the NodalValues method... My bad guys

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

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

Am I right on this @aronnoordam ?

Copy link
Copy Markdown
Member Author

Choose a reason for hiding this comment

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

yes I saw #14680 this fix, but it was not fixed in the getNodalValuesVector

@AlejandroCornejo

Copy link
Copy Markdown
Member

In general I'm good with the additions. Just a comment (take it or not), to avoid duplications of code maybe we could implement a CalculateInternalAndExternalForcesVector(bool, Internal, bool External) and compute both inside with the flags. Then in the API CalculateRightHandSide call the other one with both bools as true. What do you think?

@aronnoordam

Copy link
Copy Markdown
Member Author

@AlejandroCornejo makes sense to me, I've made the change. Note that the internal forces are now multiplied with -1, such the the calculation of the RHS is correct

@AlejandroCornejo

Copy link
Copy Markdown
Member

Thank you @aronnoordam for the additions and the fix.

@aronnoordam
aronnoordam merged commit 048e64a into master Oct 2, 2026
10 checks passed
@aronnoordam
aronnoordam deleted the struct/external_internal_force_beams branch October 2, 2026 11:02
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants