Skip to content

BUSH forces: the global-to-basic transform of the node forces used one component three times - #326

Open
Valladarex wants to merge 1 commit into
MystranSolver:devfrom
Valladarex:pr/bush-forces
Open

Valladarex wants to merge 1 commit into
MystranSolver:devfrom
Valladarex:pr/bush-forces

Conversation

@Valladarex

Copy link
Copy Markdown

Fixes #325.

One commit, in Source/LK9/L92/ELEMENT_RECOVERY_SUPPORT.f03.

  • Each grid's forces and moments are transformed with all three components.
  • Each grid writes only its own 6 terms.

Test: MYSTRAN_Validation MystranSolver/MYSTRAN_Validation#39, cbush_rotated_cd_forces.bdf: a CBUSH between coincident grids in a rotated CD, the first fixed, with a load in basic of (30, 50, 80, 6, 0, 8) on the second.

  • The six BUSH forces must equal the load.
  • On dev they are (13.7, 25.5, 39.0, 2.5, -5.8, 3.5). With this PR they equal the load.

Also checked on dev:

  • The first grid in basic with the second rotated gives the same wrong forces.
  • The first rotated with the second in basic gives the right ones, so that case is not added as a test.

Validation suite: with both PRs, 1 failure of 3,261 checks, the same one current dev has without them ("mitc4 thermal strain free rectangle cylinder"). The 6 new checks pass.

…e component three times

In CALC_ELEM_NODE_FORCES (BUSH), the node forces of a grid with a non-basic CD were taken to basic as
T0G(J,1)*PEG(J) + T0G(J,2)*PEG(J) + T0G(J,3)*PEG(J): the same component three times instead of the three
components, and when grid B was in basic the copy for it overwrote all 12 terms, grid A's too. The BUSH
engineering forces were then wrong whenever the BUSH's second grid has a rotated output system (with the first
grid rotated and the second in basic they came out right).

Each grid's forces and moments are now transformed with all three components, and each grid writes only its
own 6 terms.

Test (MYSTRAN_Validation): cbush_rotated_cd_forces.bdf, a CBUSH between coincident grids in a rotated CD, the
first fixed, a load (30, 50, 80, 6, 0, 8) in basic on the second: the BUSH forces must equal the load. Current
dev gives (13.7, 25.5, 39.0, 2.5, -5.8, 3.5).
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.

1 participant