MYSTRAN
MYSTRAN copied to clipboard
MYSTRANsolver
Rotated coordinates with SOL5 is wrong
This issue affects CBAR as well as CQUAD. Previously, it caused a crash and Bruno fixed that but the solution is still wrong.
This one has buckling factor = 5.367289E+01
SOL 105
CEND
DISPLACEMENT = ALL
FORCE = ALL
SPC = 1
SUBCASE 1
LOAD = 1
SUBCASE 2
METHOD = 1
BEGIN BULK
GRID* 1 0 0. 0.*
* 0. 0
GRID* 2 0 0.424264069 0.*
* -0.424264069 0
GRID* 3 0 0.212132034 0.*
* -0.212132034 1
$--1---><------2-------><------3-------><------4-------><------5-------><--6--->
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Material $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
MAT1* 1 2.+11 0.*
* 0. 0. 0.
$--1---><------2-------><------3-------><------4-------><------5-------><--6--->
CORD2R* 1 0 0. 0.*
* 0. -0.707106781 0. 0.707106781*
* 0.707106781 0. 0.707106781 *
*
SPC* 1 1 1 0.*
*
SPC* 1 1 2 0.*
*
SPC* 1 1 3 0.*
*
SPC* 1 1 4 0.*
*
SPC* 1 1 5 0.*
*
SPC* 1 1 6 0.*
*
SPC* 1 3 1 0.*
*
$--1---><------2-------><------3-------><------4-------><------5-------><--6--->
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Default $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
PBAR* 2 1 0.0002 6.66666667-09*
* 1.66666667-09 4.57760417-09 0.
CBAR* 1 2 1 3*
* 0. 1. 0.
CBAR* 2 2 3 2*
* 0. 1. 0.
$--1---><------2-------><------3-------><------4-------><------5-------><--6--->
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ force $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
FORCE* 1 2 0 1.*
* -1000. 0. 1000.
$--1---><------2-------><------3-------><------4-------><------5-------><--6--->
EIGRL* 1 1*
* 0 1. MASS
ENDDATA
This one should be equivalent but is aligned with the x-axis and doesn't have the CORD2R. Buckling factor = 4.111297E+00 which I think is correct.
SOL 105
CEND
DISPLACEMENT = ALL
FORCE = ALL
SPC = 1
SUBCASE 1
LOAD = 1
SUBCASE 2
METHOD = 1
BEGIN BULK
GRID* 1 0 0. 0.*
* 0. 0
GRID* 2 0 0.6 0.*
* 5.69645242-16 0
GRID* 3 0 0.3 0.*
* 1.83690953-17 0
$--1---><------2-------><------3-------><------4-------><------5-------><--6--->
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Material $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
MAT1* 1 2.+11 0.*
* 0. 0. 0.
$--1---><------2-------><------3-------><------4-------><------5-------><--6--->
SPC* 1 1 1 0.*
*
SPC* 1 1 2 0.*
*
SPC* 1 1 3 0.*
*
SPC* 1 1 4 0.*
*
SPC* 1 1 5 0.*
*
SPC* 1 1 6 0.*
*
SPC* 1 3 3 0.*
*
$--1---><------2-------><------3-------><------4-------><------5-------><--6--->
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Default $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
PBAR* 1 1 0.0002 6.66666667-09*
* 1.66666667-09 4.57760417-09 0.
CBAR* 1 1 1 3*
* 0. 1. 0.
CBAR* 2 1 3 2*
* 0. 1. 0.
$--1---><------2-------><------3-------><------4-------><------5-------><--6--->
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ force $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
FORCE* 1 2 0 1.*
* -1414.21356 0. 0.
$--1---><------2-------><------3-------><------4-------><------5-------><--6--->
EIGRL* 1 1*
* 0 1. MASS
ENDDATA
This bug also affects static (SOL 101) engineering forces (FORCE case control entry) for at least BAR elements.