사용자정의 지반모델 (User Supplied Soil Material) Home > GTS NX > Appendix > USSM(사용자정의 지반모델)

개요

• 사용자 정의 지반 재료 모델을 FORTRAN으로 구현하여, DLL로 빌드하면 GTS NX에서 사용할 수 있습니다.
• 사용자 정의 지반 재료 모델은 이전 스텝의 수렴된 응력과 상태 변수 및 현재 스텝의 변형률/시간 증분을 입력받아 현재 스텝의 응력과 상태 변수를 계산해야 합니다.

방법

서브루틴

    SUBROUTINE User_Mod (
    JOB, MODEL, UNDRAINED, STEPNUMBER, ITERNUMBER, ELEMNUMBER, 
    INTEGNUMBER, X, Y, Z, T0, DT, USRVAL, SIG0, PP0, 
    USRSTA0, DEPS, SE, BULKWATER, SIG, PP, USRSTA, PLASTICSTATUS, NUS, 
    NONSYM, STRESSDEP, TIMEDEP, TANGENT, 
    PROJECTDIRECTORY, PROJECTLENGTH, ERROR )

예제 의사 코드

1. 상태 변수 초기화

    USRSTA0(1) = (SIG0(1) + SIG0(2) + SIG0(3))/3.D0

2. 유효 응력 계산

    DO I = 1, 6

        DO J = 1, 6

           SIG(I) = SIG(I) + SE(I, J) * DEPS(J)

        ENDDO

    ENDDO

3. 유효 재료 강성 행렬 계산

    E = USRVAL(1)

    V = USRVAL(2)

    G = E/2/(1.0+V)

    FAC = 2*G/(1.0-2*V)

    VAL1 = FAC*(1-V)

    VAL2 = FAC*V



    SE = 0.0

    SE(1,1) = VAL1

    SE(1,2) = VAL2

    SE(1,3) = VAL2

    SE(2,1) = VAL2

    SE(2,2) = VAL1

    SE(2,3) = VAL2

    SE(3,1) = VAL2

    SE(3,2) = VAL2

    SE(3,3) = VAL1

    SE(4,4) = G

    SE(5,5) = G

    SE(6,6) = G

4. 상태 변수 개수 반환

    NUS = 1

5. 행렬 특성 반환(NONSYM, STRESSDEP, TIMEDEP, TANGENT)

    NONSYM = 1

    STRESSDEP = 1

    TIMEDEP = 1

    TANGENT = 1

보조 서브루틴

• dll 파일에 포함된 모델의 개수/이름과 파라미터 개수/이름/단위를 가져오는 함수.

  
      SUBROUTINE GETMODELCOUNT( NMODEL )
  
      IMPLICIT NONE
  
      INTEGER (KIND=4) NMODEL
  
      !DEC$ ATTRIBUTES DLLEXPORT, STDCALL, REFERENCE ::GETMODELCOUNT
  
      NMODEL = 1
  
      RETURN
  
  END ! GETMODELCOUNT
  
  SUBROUTINE GETMODELNAME( MODEL , MODELNAME )
  
      IMPLICIT NONE
  
      INTEGER, INTENT(IN) :: MODEL
  
      CHARACTER (LEN= * ) MODELNAME
  
      CHARACTER (LEN=255) TNAME
  
      INTEGER LT
  
      !DEC$ ATTRIBUTES DLLEXPORT, STDCALL, REFERENCE ::  GETMODELNAME
  
      SELECT CASE (MODEL)
  
      CASE(1)
  
          TNAME = 'ELASTIC'
  
          LT = LEN_TRIM(TNAME)
  
          MODELNAME= CHAR(LT) // TNAME(1:LT)
  
      END SELECT ! MODEL
  
      RETURN
  
  END ! GETMODELNAME
  
  SUBROUTINE GETPARAMCOUNT( MODEL , NVAL )
  
      IMPLICIT NONE
  
      INTEGER  MODEL, NVAL
  
      !DEC$ ATTRIBUTES DLLEXPORT, STDCALL, REFERENCE ::  GETPARAMCOUNT
  
      NVAL = 2
  
  END ! GETPARAMCOUNT
  
  SUBROUTINE GETPARAMNAME( MODEL, PARAM, PARAMNAME )
  
      IMPLICIT NONE
  
      INTEGER, INTENT(IN) :: MODEL
  
      INTEGER, INTENT(IN) :: PARAM
  
      CHARACTER (LEN= * ) PARAMNAME
  
      CHARACTER (LEN=255) TNAME
  
      INTEGER LT
  
      !DEC$ ATTRIBUTES DLLEXPORT, STDCALL, REFERENCE ::  GETPARAMNAME
  
      SELECT CASE (MODEL)
  
      CASE(1)
  
          SELECT CASE (PARAM)
  
          CASE (1)
  
          TNAME = 'E'
  
          CASE (2)
  
          TNAME = 'NU'
  
          END SELECT ! PARAM
  
      END SELECT ! MODEL
  
      LT = LEN_TRIM(TNAME)
  
      PARAMNAME= CHAR(LT) // TNAME(1:LT)
  
      RETURN
  
  END ! GETPARAMNAME
  
  SUBROUTINE GETPARAMUNIT( MODEL, PARAM, PARAMUNIT )
  
      IMPLICIT NONE
  
      INTEGER, INTENT(IN) :: MODEL
  
      INTEGER, INTENT(IN) :: PARAM
  
      CHARACTER (LEN= * ) PARAMUNIT
  
      CHARACTER (LEN=255) TNAME
  
      INTEGER LT
  
      !DEC$ ATTRIBUTES DLLEXPORT, STDCALL, REFERENCE ::  GETPARAMUNIT
  
      SELECT CASE (MODEL)
  
      CASE(1)
  
          SELECT CASE (PARAM)
  
          CASE (1)
  
          TNAME = 'F/L^2#'
  
          CASE (2)
  
          TNAME = ''
  
          END SELECT ! PARAM
  
      END SELECT ! MODEL
  
      LT = LEN_TRIM(TNAME)
  
      PARAMUNIT= CHAR(LT) // TNAME(1:LT)
  
      RETURN
  
  END ! GETPARAMUNIT