structure.f90

边界元程序,供力学工作者参考,希望对大家有所帮助

        SUBROUTINE UpdateStructure
        USE CtrlData
        USE ElmtData
        USE GlobData
        iswth = mswth * 2 + 7
        DO ielem = 1, nelem
          CALL ElmOpt( ielem, iswth )
          IF( nerrc .GT. 0 ) RETURN
        END DO
        DO idofs = 1, ndofs
          disph( idofs, kload ) = disph( idofs, kload ) +                      &
          disps( idofs, kload )
          disps( idofs, kload ) = 0.0D0
        END DO
        END