Skip to content

Mitc8 with stress/strain/force output#147

Merged
zchlrnr merged 13 commits intoMystranSolver:mainfrom
victorkemp:mitc8
Jun 30, 2025
Merged

Mitc8 with stress/strain/force output#147
zchlrnr merged 13 commits intoMystranSolver:mainfrom
victorkemp:mitc8

Conversation

@victorkemp
Copy link
Copy Markdown
Collaborator

This is hopefully MITC8 finished in a basic useable state. I've removed the DEBUG,250,1 requirement so it should work with normal input files using the new CQUAD8 element.

Description of CQUAD8

   1       2        3       4      5        6       7      8       9      10
<------><------><------><------><------><------><------><------><------><------>
CQUAD8    EID      PID     G1      G2      G3      G4      G5      G6
           G7       G8

PID - must be PSHELL ID, not PCOMP or PCOMP1.
G1-G8 are all required. No optional G5-G8 like MSC.

PSHELL referred to by PID:

  • Must have MID2 = MID1
  • MID3 can be different from MID1/2
  • TS/T can be used
  • Must have blank or default values for:
    • 12I/TM**3
    • NSM (MPA)
    • Z1
    • Z2
    • MID4

Material must be MAT1 (isotropic) with zero density.

Solution type must be statics (SOL 1)

Not compatible with:

  • Thermal load
  • Pressure load
  • Gravity load

Stress, strain and engineering force output:

  • Only F06, not OP2 or NEU.
  • In the element coordinate system which is the same as MSC's.
  • At corner grid points and center regardless of whether CENTER or CORNER is specified in case control, same as MSC.
  • The center is at the centroid of the 4 corner grid points and the center values are the averages of the corner values.

I've also made a small change to the engineering force format in F06 for all shells to accommodate corner values. The Element ID column is further left and the header includes a "Location" column that's only populated for CQUAD8 elements. Hopefully this won't break any parsers that use space delimiters.
Previous:

                                                       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         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  1.150000E+01  2.000000E+00  3.046240E-03  2.500000E-02 -1.396284E-03  2.486571E-03  1.000000E-01
                       2  0.000000E+00  1.450000E+01  2.000000E+00  4.733241E-04  7.500000E-02 -4.358266E-04 -6.376932E-04  1.000000E-01
                       3  0.000000E+00  5.500000E+00  2.000000E+00  4.733241E-04  7.500000E-02  4.358266E-04  6.376932E-04  1.000000E-01
                       4  0.000000E+00  8.500000E+00  2.000000E+00  3.046240E-03  2.500000E-02  1.396284E-03 -2.486571E-03  1.000000E-01
                          ------------- ------------- ------------- ------------- ------------- ------------- ------------- -------------

This PR:

                                                       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  1.150000E+01  2.000000E+00  3.445167E-03  2.500000E-02 -1.629156E-03  2.520876E-03  1.000000E-01
        2                 0.000000E+00  1.450000E+01  2.000000E+00  4.161607E-04  7.500000E-02 -4.253891E-04 -2.104274E-03  1.000000E-01
        3                 0.000000E+00  5.500000E+00  2.000000E+00  4.161607E-04  7.500000E-02  4.253891E-04  2.104274E-03  1.000000E-01
        4                 0.000000E+00  8.500000E+00  2.000000E+00  3.445167E-03  2.500000E-02  1.629156E-03 -2.520876E-03  1.000000E-01
                          ------------- ------------- ------------- ------------- ------------- ------------- ------------- -------------

@zchlrnr
Copy link
Copy Markdown
Collaborator

zchlrnr commented Jun 30, 2025

Watershed moment.

@zchlrnr zchlrnr merged commit 81213b5 into MystranSolver:main Jun 30, 2025
1 check passed
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.

2 participants