Skip to content

OP2 solid corner stress is center stress duplicated #246

@victorkemp

Description

@victorkemp

Problem should be in WRITE_ELEM_STRESSES.f90:

           ! TODO: we repeat the center node N times because the corner results have not been calculated
           WRITE(OP2) (EID_OUT_ARRAY(I,1)*10+DEVICE_CODE, CID, CEN_WORD, NNODES-1,             &
                        ! grid_id
                        ! 21
                       (GID_OUT_ARRAY(I,J),                                                    &
                        ! oxx             txy                s1                  a1  a2  a3
                       REAL(OGEL(I,1),4), REAL(OGEL(I,4),4), REAL(OGEL(I,9), 4), 0., 0., 0.,   &
                        ! p                   ovm
                       REAL(OGEL(I,12),4), REAL(OGEL(I,7),4),  &
                        ! syy             tyz                s2                  b1  b2  b3
                       REAL(OGEL(I,2),4), REAL(OGEL(I,5),4), REAL(OGEL(I,10),4), 0., 0., 0.,   &
                        ! szz             txz                s3                  c1  c2  c3
                       REAL(OGEL(I,3),4), REAL(OGEL(I,6),4), REAL(OGEL(I,11),4), 0., 0., 0.,   &
                       J=1,NNODES), I=1,NUM)

J is not used to index OGEL, so it writes the same values for each iteration of the J loop. OGEL does contain corner stress which is used correctly for F06 in the code below that.

Metadata

Metadata

Assignees

No one assigned

    Labels

    bugSomething isn't working

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions