matlib.f90

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

!       ===============================================================
!       |iswth 调用功能开关
!       |      ＝＝ 1   读入材料参数并获取材料参数
!       |      ＝＝ 2   计算切线模量
!       |      ＝＝ 3   由应变增量计算应力增量
!       |      ＝＝ 4
!       |      ＝＝ 5
!       |      ＝＝ 6
!       |      ＝＝ 7
!       |      ＝＝ 8   计算历史量
!       ===============================================================

        SUBROUTINE MatLib( dmatx, strsg, strng, dstrs, dstrn, ietyp,           &
                           imtyp, imats, iposp, iposh, iswth, nstrg )
        USE CtrlData
        USE MeshData
        USE GlobData
        IMPLICIT DOUBLE PRECISION( a-h, o-z )
        DIMENSION strsg( nstrg ), dstrs( nstrg )
        DIMENSION strng( nstrg ), dstrn( nstrg ), dmatx( nstrg, nstrg )

        DO istrg = 1, nstrg
          DO jstrg = 1, nstrg
            dmatx( istrg, jstrg ) = 0.0D0
          END DO
        END DO

        SELECT CASE( imtyp )
        CASE( 1 )
          CALL Mat001( dmatx, strsg, strng, dstrs, dstrn, ietyp,           &
                       imats, iposh, iposp, iswth )
        CASE( 301 )
          CALL Mat301( dmatx, strsg, strng, dstrs, dstrn, ietyp,           &
                       imats, iposh, iposp, iswth, nstrg )
        END SELECT
        END