Skip to content

test for fix to STR_CID for solids #48

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.

Sign up for GitHub

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

Merged
victorkemp merged 1 commit into from
Oct 24, 2026
Merged

test for fix to STR_CID for solids #48

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 4 additions & 0 deletions cases_new.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2,6 +2,10 @@
# Put new test case definitions at the top so they run first and you can see the results easier.
# ==============================================================================================

# Bug of wrong transform for solid stress transformation
pth; hexa_rotated_cordm_stress.bdf; SC/1/SOLIDSTRESSES/EID/1-7/CORNER/0/XX;; -1; 5e-5%
pth; hexa_rotated_cordm_stress.bdf; SC/1/SOLIDSTRESSES/EID/1-7/CORNER/0/YY;; 0; 1e-13

# Bug of using only the first 4 nodes of CHEXA for thermal stress
pth; chexa_linear_temperature_stress.bdf; SC/1/SOLIDSTRESSES/EID/1,2/CORNER/0/ZZ;; 0; 1e-6
pth; chexa_linear_temperature_stress.bdf; SC/1/SOLIDSTRESSES/EID/1,2/CORNER/1-4/YY;; -2e9; 5e-5%
Expand Down
305 changes: 305 additions & 0 deletions decks/hexa_rotated_cordm_stress.bdf
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,305 @@
$ Generated by Mecway 35 for Mystran
$-----------------------EXECUTIVE CONTROL SECTION-----------------------
SOL 101
CEND
$-------------------------CASE CONTROL SECTION--------------------------
DISPLACEMENT = ALL
SPCFORCES = ALL
MPCFORCES = ALL
STRESS(CORNER) = ALL
STRAIN(FIBER, CORNER) = ALL
SPC = 1
LOAD = 1
BEGIN BULK
$---------------------------BULK DATA SECTION---------------------------
PARAM, POST,-1
GRID 1 0 4.0+0 -5.0-1 0.0+0 0
GRID 2 0 5.0+0 -5.0-1 0.0+0 0
GRID 3 0 5.0+0 5.0-1 0.0+0 0
GRID 4 0 4.0+0 5.0-1 0.0+0 0
GRID 5 0 4.0+0 -5.0-1 1.0+0 0
GRID 6 0 5.0+0 -5.0-1 1.0+0 0
GRID 7 0 5.0+0 5.0-1 1.0+0 0
GRID 8 0 4.0+0 5.0-1 1.0+0 0
GRID* 9 0 3.99311282771+0 5.52313267266-1*
* 0.0+0 0
GRID* 10 0 4.959038654+0 8.11132312368-1*
* 0.0+0 0
GRID* 11 0 4.70021960889+0 1.77705813866+0*
* 0.0+0 0
GRID* 12 0 3.73429378261+0 1.51823909355+0*
* 0.0+0 0
GRID* 13 0 3.99311282771+0 5.52313267266-1*
* 1.0+0 0
GRID* 14 0 4.959038654+0 8.11132312368-1*
* 1.0+0 0
GRID* 15 0 4.70021960889+0 1.77705813866+0*
* 1.0+0 0
GRID* 16 0 3.73429378261+0 1.51823909355+0*
* 1.0+0 0
GRID* 17 0 3.71410161514+0 1.56698729811+0*
* 0.0+0 0
GRID* 18 0 4.58012701892+0 2.06698729811+0*
* 0.0+0 0
GRID* 19 0 4.08012701892+0 2.93301270189+0*
* 0.0+0 0
GRID* 20 0 3.21410161514+0 2.43301270189+0*
* 0.0+0 0
GRID* 21 0 3.71410161514+0 1.56698729811+0*
* 1.0+0 0
GRID* 22 0 4.58012701892+0 2.06698729811+0*
* 1.0+0 0
GRID* 23 0 4.08012701892+0 2.93301270189+0*
* 1.0+0 0
GRID* 24 0 3.21410161514+0 2.43301270189+0*
* 1.0+0 0
GRID* 25 0 3.18198051534+0 2.47487373415+0*
* 0.0+0 0
GRID* 26 0 3.88908729653+0 3.18198051534+0*
* 0.0+0 0
GRID* 27 0 3.18198051534+0 3.88908729653+0*
* 0.0+0 0
GRID* 28 0 2.47487373415+0 3.18198051534+0*
* 0.0+0 0
GRID* 29 0 3.18198051534+0 2.47487373415+0*
* 1.0+0 0
GRID* 30 0 3.88908729653+0 3.18198051534+0*
* 1.0+0 0
GRID* 31 0 3.18198051534+0 3.88908729653+0*
* 1.0+0 0
GRID* 32 0 2.47487373415+0 3.18198051534+0*
* 1.0+0 0
GRID* 33 0 2.43301270189+0 3.21410161514+0*
* 0.0+0 0
GRID* 34 0 2.93301270189+0 4.08012701892+0*
* 0.0+0 0
GRID* 35 0 2.06698729811+0 4.58012701892+0*
* 0.0+0 0
GRID* 36 0 1.56698729811+0 3.71410161514+0*
* 0.0+0 0
GRID* 37 0 2.43301270189+0 3.21410161514+0*
* 1.0+0 0
GRID* 38 0 2.93301270189+0 4.08012701892+0*
* 1.0+0 0
GRID* 39 0 2.06698729811+0 4.58012701892+0*
* 1.0+0 0
GRID* 40 0 1.56698729811+0 3.71410161514+0*
* 1.0+0 0
GRID* 41 0 1.51823909355+0 3.73429378261+0*
* 0.0+0 0
GRID* 42 0 1.77705813866+0 4.70021960889+0*
* 0.0+0 0
GRID* 43 0 8.11132312368-1 4.959038654+0*
* 0.0+0 0
GRID* 44 0 5.52313267266-1 3.99311282771+0*
* 0.0+0 0
GRID* 45 0 1.51823909355+0 3.73429378261+0*
* 1.0+0 0
GRID* 46 0 1.77705813866+0 4.70021960889+0*
* 1.0+0 0
GRID* 47 0 8.11132312368-1 4.959038654+0*
* 1.0+0 0
GRID* 48 0 5.52313267266-1 3.99311282771+0*
* 1.0+0 0
GRID 49 0 5.0-1 4.0+0 0.0+0 0
GRID 50 0 5.0-1 5.0+0 0.0+0 0
GRID 51 0 -5.0-1 5.0+0 0.0+0 0
GRID 52 0 -5.0-1 4.0+0 0.0+0 0
GRID 53 0 5.0-1 4.0+0 1.0+0 0
GRID 54 0 5.0-1 5.0+0 1.0+0 0
GRID 55 0 -5.0-1 5.0+0 1.0+0 0
GRID 56 0 -5.0-1 4.0+0 1.0+0 0
$--1---><--2---><--3---><--4---><--5---><--6---><--7---><--8---><--9---><--10-->
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Material $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
MAT1 1 1.0+11 5.0+10 0.0+0 0.0+0 0.0+0 0.0+0
$--1---><--2---><--3---><--4---><--5---><--6---><--7---><--8---><--9---><--10-->
SPC 1 1 1 0.0+0
SPC 1 1 2 0.0+0
SPC 1 1 3 0.0+0
SPC 1 4 1 0.0+0
SPC 1 4 2 0.0+0
SPC 1 4 3 0.0+0
SPC 1 5 1 0.0+0
SPC 1 5 2 0.0+0
SPC 1 5 3 0.0+0
SPC 1 8 1 0.0+0
SPC 1 8 2 0.0+0
SPC 1 8 3 0.0+0
SPC 1 9 1 0.0+0
SPC 1 9 2 0.0+0
SPC 1 9 3 0.0+0
SPC 1 12 1 0.0+0
SPC 1 12 2 0.0+0
SPC 1 12 3 0.0+0
SPC 1 13 1 0.0+0
SPC 1 13 2 0.0+0
SPC 1 13 3 0.0+0
SPC 1 16 1 0.0+0
SPC 1 16 2 0.0+0
SPC 1 16 3 0.0+0
SPC 1 17 1 0.0+0
SPC 1 17 2 0.0+0
SPC 1 17 3 0.0+0
SPC 1 20 1 0.0+0
SPC 1 20 2 0.0+0
SPC 1 20 3 0.0+0
SPC 1 21 1 0.0+0
SPC 1 21 2 0.0+0
SPC 1 21 3 0.0+0
SPC 1 24 1 0.0+0
SPC 1 24 2 0.0+0
SPC 1 24 3 0.0+0
SPC 1 25 1 0.0+0
SPC 1 25 2 0.0+0
SPC 1 25 3 0.0+0
SPC 1 28 1 0.0+0
SPC 1 28 2 0.0+0
SPC 1 28 3 0.0+0
SPC 1 29 1 0.0+0
SPC 1 29 2 0.0+0
SPC 1 29 3 0.0+0
SPC 1 32 1 0.0+0
SPC 1 32 2 0.0+0
SPC 1 32 3 0.0+0
SPC 1 33 1 0.0+0
SPC 1 33 2 0.0+0
SPC 1 33 3 0.0+0
SPC 1 36 1 0.0+0
SPC 1 36 2 0.0+0
SPC 1 36 3 0.0+0
SPC 1 37 1 0.0+0
SPC 1 37 2 0.0+0
SPC 1 37 3 0.0+0
SPC 1 40 1 0.0+0
SPC 1 40 2 0.0+0
SPC 1 40 3 0.0+0
SPC 1 41 1 0.0+0
SPC 1 41 2 0.0+0
SPC 1 41 3 0.0+0
SPC 1 44 1 0.0+0
SPC 1 44 2 0.0+0
SPC 1 44 3 0.0+0
SPC 1 45 1 0.0+0
SPC 1 45 2 0.0+0
SPC 1 45 3 0.0+0
SPC 1 48 1 0.0+0
SPC 1 48 2 0.0+0
SPC 1 48 3 0.0+0
SPC 1 49 1 0.0+0
SPC 1 49 2 0.0+0
SPC 1 49 3 0.0+0
SPC 1 52 1 0.0+0
SPC 1 52 2 0.0+0
SPC 1 52 3 0.0+0
SPC 1 53 1 0.0+0
SPC 1 53 2 0.0+0
SPC 1 53 3 0.0+0
SPC 1 56 1 0.0+0
SPC 1 56 2 0.0+0
SPC 1 56 3 0.0+0
$--1---><--2---><--3---><--4---><--5---><--6---><--7---><--8---><--9---><--10-->
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Default $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
PSOLID 1 1 0 2 1
CHEXA 1 1 1 2 3 4 5 6+
+ 7 8
CORD2R* 2 0 0.0+0 0.0+0*
* 0.0+0 0.0+0 0.0+0 9.99999999999-1*
* 9.65925826289-1 2.58819045102-1 0.0+0 *
*
PSOLID 3 1 2 2 1
CHEXA 2 3 9 10 11 12 13 14+
+ 15 16
CORD2R* 4 0 0.0+0 0.0+0*
* 0.0+0 0.0+0 0.0+0 1.0+0*
* 8.66025403784-1 5.0-1 0.0+0 *
*
PSOLID 5 1 4 2 1
CHEXA 3 5 17 18 19 20 21 22+
+ 23 24
CORD2R* 6 0 0.0+0 0.0+0*
* 0.0+0 0.0+0 0.0+0 1.0+0*
* 7.07106781186-1 7.07106781186-1 0.0+0 *
*
PSOLID 7 1 6 2 1
CHEXA 4 7 25 26 27 28 29 30+
+ 31 32
CORD2R* 8 0 0.0+0 0.0+0*
* 0.0+0 0.0+0 0.0+0 1.0+0*
* 5.0-1 8.66025403784-1 0.0+0 *
*
PSOLID 9 1 8 2 1
CHEXA 5 9 33 34 35 36 37 38+
+ 39 40
CORD2R* 10 0 0.0+0 0.0+0*
* 0.0+0 0.0+0 0.0+0 9.99999999999-1*
* 2.58819045102-1 9.65925826289-1 0.0+0 *
*
PSOLID 11 1 10 2 1
CHEXA 6 11 41 42 43 44 45 46+
+ 47 48
CORD2R 12 0 0.0+0 0.0+0 0.0+0 0.0+0 0.0+0 1.0+0+
+ 0.0+0 1.0+0 0.0+0
PSOLID 13 1 12 2 1
CHEXA 7 13 49 50 51 52 53 54+
+ 55 56
$--1---><--2---><--3---><--4---><--5---><--6---><--7---><--8---><--9---><--10-->
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ pressure $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
FORCE* 2 26 0 1.0+0*
* -1.76776695297-1-1.76776695297-1 5.3619654621-17
FORCE* 2 27 0 1.0+0*
* -1.76776695297-1-1.76776695297-1-1.4367343154-17
FORCE* 2 30 0 1.0+0*
* -1.76776695297-1-1.76776695297-1 5.3619654621-17
FORCE* 2 31 0 1.0+0*
* -1.76776695297-1-1.76776695297-1-1.4367343154-17
FORCE* 2 18 0 1.0+0*
* -2.16506350946-1 -1.25-1 8.0123445266-18
FORCE* 2 19 0 1.0+0*
* -2.16506350946-1 -1.25-1 2.9902476861-17
FORCE* 2 22 0 1.0+0*
* -2.16506350946-1 -1.25-1 8.0123445266-18
FORCE* 2 23 0 1.0+0*
* -2.16506350946-1 -1.25-1 2.9902476861-17
FORCE* 2 2 0 1.0+0*
* -2.5-1 0.0+0 1.1730886562-17
FORCE* 2 3 0 1.0+0*
* -2.5-1 0.0+0 4.3780264669-17
FORCE* 2 6 0 1.0+0*
* -2.5-1 0.0+0 1.1730886562-17
FORCE* 2 7 0 1.0+0*
* -2.5-1 0.0+0 4.3780264669-17
FORCE* 2 10 0 1.0+0*
* -2.41481456572-1-6.47047612756-2-9.4769576625-17
FORCE* 2 11 0 1.0+0*
* -2.41481456572-1-6.47047612756-2-6.0701633556-17
FORCE* 2 14 0 1.0+0*
* -2.41481456572-1-6.47047612756-2-9.4769576625-17
FORCE* 2 15 0 1.0+0*
* -2.41481456572-1-6.47047612756-2-6.0701633556-17
FORCE* 2 34 0 1.0+0*
* -1.25-1-2.16506350946-1 2.9327216406-18
FORCE* 2 35 0 1.0+0*
* -1.25-1-2.16506350946-1 1.0945066167-17
FORCE* 2 38 0 1.0+0*
* -1.25-1-2.16506350946-1 2.9327216406-18
FORCE* 2 39 0 1.0+0*
* -1.25-1-2.16506350946-1 1.0945066167-17
FORCE* 2 42 0 1.0+0*
* -6.47047612756-2-2.41481456572-1-8.9483515584-17
FORCE* 2 43 0 1.0+0*
* -6.47047612756-2-2.41481456572-1-4.0973785181-17
FORCE* 2 46 0 1.0+0*
* -6.47047612756-2-2.41481456572-1-8.9483515584-17
FORCE* 2 47 0 1.0+0*
* -6.47047612756-2-2.41481456572-1-4.0973785181-17
FORCE* 2 50 0 1.0+0*
* 0.0+0 -2.5-1 1.1730886562-17
FORCE* 2 51 0 1.0+0*
* 0.0+0 -2.5-1 4.3780264669-17
FORCE* 2 54 0 1.0+0*
* 0.0+0 -2.5-1 1.1730886562-17
FORCE* 2 55 0 1.0+0*
* 0.0+0 -2.5-1 4.3780264669-17
$--1---><--2---><--3---><--4---><--5---><--6---><--7---><--8---><--9---><--10-->
LOAD 1 1.0+0 1.0+0 2
ENDDATA