With MSC, transverse shear element force (per unit length) is correct (1.0e+03) but in Mystran, it's 0.0.
Input:
$ Generated by Mecway 29
$-----------------------EXECUTIVE CONTROL SECTION-----------------------
SOL 101
CEND
$-------------------------CASE CONTROL SECTION--------------------------
DISPLACEMENT = ALL
SPCFORCES = ALL
MPCFORCES = ALL
STRESS(CORNER) = ALL
STRAIN(CORNER) = ALL
FORCE(CORNER) = ALL
SPC = 1
LOAD = 1
BEGIN BULK
$---------------------------BULK DATA SECTION---------------------------
GRID* 1 0 0.0+00 0.0+00*
* 0.0+00 0
GRID* 2 0 1.0-03 0.0+00*
* 0.0+00 0
GRID* 3 0 1.0-03 1.0-03*
* 0.0+00 0
GRID* 4 0 0.0+00 1.0-03*
* 0.0+00 0
$--1---><------2-------><------3-------><------4-------><------5-------><--6--->
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Material $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
MAT1* 1 2.0+11 1.0+11 0.0+00*
* 0.0+00 0.0+00 0.0+00
$--1---><------2-------><------3-------><------4-------><------5-------><--6--->
SPC* 1 1 1 0.0+00*
*
SPC* 1 1 2 0.0+00*
*
SPC* 1 1 3 0.0+00*
*
SPC* 1 1 4 0.0+00*
*
SPC* 1 1 5 0.0+00*
*
SPC* 1 1 6 0.0+00*
*
SPC* 1 2 1 0.0+00*
*
SPC* 1 2 2 0.0+00*
*
SPC* 1 2 3 0.0+00*
*
SPC* 1 2 4 0.0+00*
*
SPC* 1 2 5 0.0+00*
*
SPC* 1 2 6 0.0+00*
*
$--1---><------2-------><------3-------><------4-------><------5-------><--6--->
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Default $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
PSHELL* 1 1 1.0-04 1*
* 0.0+00
CQUAD4* 1 1 1 2*
* 3 4 0.0+00 0.0+00
$--1---><------2-------><------3-------><------4-------><------5-------><--6--->
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ force $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
FORCE* 1 3 0 1.0+00*
* 0.0+00 0.0+00 5.0-01
FORCE* 1 4 0 1.0+00*
* 0.0+00 0.0+00 5.0-01
$--1---><------2-------><------3-------><------4-------><------5-------><--6--->
ENDDATA
MSC output:
F O R C E S I N Q U A D R I L A T E R A L E L E M E N T S ( Q U A D 4 ) OPTION = BILIN
ELEMENT - MEMBRANE FORCES - - BENDING MOMENTS - - TRANSVERSE SHEAR FORCES -
ID GRID-ID FX FY FXY MX MY MXY QX QY
0 1 CEN/4 0.0 0.0 0.0 -8.535399E-17 5.000000E-01 -4.368062E-16 -2.273737E-13 1.000000E+03
1 0.0 0.0 0.0 9.097149E-33 5.000000E-01 -4.368062E-16 -2.273737E-13 1.000000E+03
2 0.0 0.0 0.0 4.548574E-33 5.000000E-01 -4.368062E-16 -2.273737E-13 1.000000E+03
3 0.0 0.0 0.0 -1.707080E-16 5.000000E-01 -4.506840E-16 -2.273737E-13 1.000000E+03
4 0.0 0.0 0.0 -1.707080E-16 5.000000E-01 -4.368062E-16 -2.273737E-13 1.000000E+03
Mystran output:
E L E M E N T E N G I N E E R I N G F O R C E S
F O R E L E M E N T T Y P E Q U A D 4
Element Location N o r m a l F o r c e s M o m e n t s T r a n s v e r s e
ID S h e a r F o r c e s
Nxx Nyy Nxy Mxx Myy Mxy Qx Qy
1 0.000000E+00 0.000000E+00 0.000000E+00 1.035782E-01 5.000000E-01 -1.221245E-16 0.000000E+00 0.000000E+00
GRD 1 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00 5.000000E-01 -5.178908E-02 0.000000E+00 0.000000E+00
GRD 2 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00 5.000000E-01 5.178908E-02 0.000000E+00 0.000000E+00
GRD 3 0.000000E+00 0.000000E+00 0.000000E+00 2.071563E-01 5.000000E-01 5.178908E-02 0.000000E+00 0.000000E+00
GRD 4 0.000000E+00 0.000000E+00 0.000000E+00 2.071563E-01 5.000000E-01 -5.178908E-02 0.000000E+00 0.000000E+00
With MSC, transverse shear element force (per unit length) is correct (1.0e+03) but in Mystran, it's 0.0.
Input:
MSC output:
Mystran output: