support.f90

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

        SUBROUTINE Support
        USE CtrlData
        USE MeshData
        USE SolvData
        USE GlobData
        USE FactorData

        ktype =-1
        DO ipoin = 1, npoin
          DO idofn = 1, mdofn
            IF( lfixs( idofn, ipoin ) .GT. 0 ) THEN
              IF( ktype .NE. lfixs( idofn, ipoin ) - 1 ) THEN
                ktype = lfixs( idofn, ipoin ) - 1
                CALL Factor( ktype )
              END IF
              idofs = ( ipoin - 1 ) * mdofn + idofn
              sprts( idofn, ipoin ) = vfixs( idofn, ipoin ) * fctor
              IF( nincc .NE. 0 )      sprts( idofn, ipoin ) =                  &
              sprts( idofn, ipoin ) - disps( idofs, kload )
            END IF
          END DO
        END DO
        RETURN
        END