MultiPhysicsVault/.raw/MITC공부/MITC공부_001.md

1. MITC shell element

• 3차원 솔리드 형상으로부터 쉘형상을 표현 (유한요소 정식화가 다른 쉘요소에 비해 간단)
• 쉘 이론을 사용하지 않고 3차원 응력, 변형률을 사용하여 쉘을 표현할 수 있다.
  임의의 형상에 대한 두꺼운 쉘과 얇은 쉘 모두 적용 가능
  Locking을 방지하기 위해 횡방향 전단 변형률에 보간법 사용

2. Kinematics

flowchart

graph TD
    A["0 Ω"] --> B["0 v_n²"]
    B --> C["t u"]
    C --> D["t + Δt u"]
    D --> E["4"]
    E --> F["3"]
    F --> G["2"]
    G --> H["1"]
    H --> I["t + Δt v_n²"]
    I --> J["t + Δt Ω"]
    J --> K["t + Δt u"]
    K --> L["t + Δt v_n²"]
    L --> M["t + Δt u"]
    M --> N["t + Δt v_n²"]
    N --> O["t + Δt u"]
    O --> P["t + Δt v_n²"]
    P --> Q["t + Δt u"]
    Q --> R["t + Δt v_n²"]
    R --> S["t + Δt u"]
    S --> T["t + Δt v_n²"]
    T --> U["t + Δt u"]
    U --> V["t + Δt v_n²"]
    V --> W["t + Δt u"]
    W --> X["t + Δt v_n²"]
    X --> Y["t + Δt u"]
    Y --> Z["t + Δt v_n²"]
    Z --> AA["t + Δt u"]
    AA --> AB["t + Δt v_n²"]
    AB --> AC["t + Δt u"]
    AC --> AD["t + Δt v_n²"]
    AD --> AE["t + Δt u"]
    AE --> AF["t + Δt v_n²"]
    AF --> AG["t + Δt u"]
    AG --> AH["t + Δt v_n²"]
    AH --> AI["t + Δt u"]
    AI --> AJ["t + Δt v_n²"]
    AJ --> AK["t + Δt u"]
    AK --> AL["t + Δt v_n²"]
    AL --> AM["t + Δt u"]
    AM --> AN["t + Δt v_n²"]
    AN --> AO["t + Δt u"]
    AO --> AP["t + Δt v_n²"]
    AP --> AQ["t + Δt u"]
    AQ --> AR["t + Δt v_n²"]
    AR --> AS["t + Δt u"]
    AS --> AT["t + Δt v_n²"]
    AT --> AU["t + Δt u"]
    AU --> AV["t + Δt v_n²"]
    AV --> AW["t + Δt u"]
    AW --> AX["t + Δt v_n²"]
    AX --> AY["t + Δt u"]
    AY --> AZ["t + Δt v_n²"]
    AZ --> BA["t + Δt u"]
    BA --> BB["t + Δt v_n²"]
    BB --> BC["t + Δt u"]
    BC --> BD["t + Δt v_n²"]
    BD --> BE["t + Δt u"]
    BE --> BF["t + Δt v_n²"]
    BF --> BG["t + Δt u"]
    BG --> BH["t + Δt v_n²"]
    BH --> BI["t + Δt u"]
    BI --> BJ["t + Δt v_n²"]
    BJ --> BK["t + Δt u"]
    BK --> BL["t + Δt v_n²"]
    BL --> BM["t + Δt u"]
    BM --> BN["t + Δt v_n²"]
    BN --> BO["t + Δt u"]
    BO --> BP["t + Δt v_n²"]
    BP --> BQ["t + Δt u"]
    BQ --> BR["t + Δt v_n²"]
    BR --> BS["t + Δt u"]
    BS --> BT["t + Δt v_n²"]
    BT --> BU["t + Δt u"]
    BU --> BV["t + Δt v_n²"]
    BV --> BW["t + Δt u"]
    BW --> BX["t + Δt v_n²"]
    BX --> BY["t + Δt u"]
    BY --> BZ["t + Δt v_n²"]
    BZ --> CA["t + Δt u"]
    CA --> CB["t + Δt v_n²"]
    CB --> CC["t + Δt u"]
    CC --> CD["t + Δt v_n²"]
    CD --> CE["t + Δt u"]
    CE --> CF["t + Δt v_n²"]
    CF --> CG["t + Δt u"]
    CG --> CH["t + Δt v_n²"]
    CH --> CI["t + Δt u"]
    CI --> CJ["t + Δt v_n²"]
    CJ --> CK["t + Δt u"]
    CK --> CL["t + Δt v_n²"]
    CL --> CM["t + Δt u"]
    CM --> CN["t + Δt v_n²"]
    CN --> CO["t + Δt u"]
    CO --> CP["t + Δt v_n²"]
    CP --> CQ["t + Δt u"]
    CQ --> CR["t + Δt v_n²"]
    CR --> CS["t + Δt u"]
    CS --> CT["t + Δt v_n²"]
    CT --> CU["t + Δt u"]
    CU --> CV["t + Δt v_n²"]
    CV --> CW["t + Δt u"]
    CW --> CX["t + Δt v_n²"]
    CX --> CY["t + Δt u"]
    CY --> CZ["t + Δt v_n²"]
    CZ --> DA["t + Δt u"]
    DA --> DB["t + Δt v_n²"]
    DB --> DC["t + Δt u"]
    DC --> DD["t + Δt v_n²"]
    DD --> DE["t + Δt u"]
    DE --> DF["t + Δt v_n²"]
    DF --> DG["t + Δt u"]
    DG --> DH["t + Δt v_n²"]
    DH --> DI["t + Δt u"]
    DI --> DJ["t + Δt v_n²"]
    DJ --> DK["t + Δt u"]
    DK --> DL["t + Δt v_n²"]
    DL --> DM["t + Δt u"]
    DM --> DN["t + Δt v_n²"]
    DN --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    D --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    Do --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    EO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    O --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    AO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    BO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    OD --> DO
    OD --> DO
    OD --> DO
    OD --> DO
    OD --> DO
    OD --> DO
    OD --> DO
    OD --> DO
    OD --> DO
    OD --> DO
    OD --> DO
    OD --> DO
    OD --> DO
    OD --> DO
    OD --> DO
    OD --> DO
    OD --> DO
    OD --> DO
    OD --> DO
    OD --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DN --> DO
    DN --> DO
    DN --> DO
    DN --> DO
    DN --> DO
    DN --> DO
    DN --> DO
    DN --> DO
    DN --> DO
    DN --> DO
    DN --> DO
    DN --> DO
    DN --> DO
    DN --> DO
    DN --> DO
    DN --> DO
    DN --> DO
    DN --> DO
    DN --> DO
    DN --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO
    DO --> DO

text_image

mid surface t+Δt E₃ t+Δt Vₙ² α₁² t+Δt E₂ α₂² t+Δt V₂² t+Δt E₁ h₂ t+Δt x

Shell의 초기 위치 벡터는 다음과 같이 shape function으로 나타낼 수 있다.

{ } ^ { 0 } \mathbf { X } = \sum _ { i = 1 } ^ { 4 } \phi _ { i } \left( \xi ^ { 1 } , \xi ^ { 2 } \right) { } ^ { 0 } \mathbf { X } _ { i } + \frac { \xi ^ { 3 } } { 2 } \sum _ { i = 1 } ^ { 4 } h _ { i } \phi _ { i } \left( \xi ^ { 1 } , \xi ^ { 2 } \right) { } ^ { 0 } \mathbf { V } _ { n } ^ { i }

마찬가지로 시간이 t, t + \Delta t 일 때 위치벡터는 다음과 같다.

{ } ^ { t } \mathbf { x } = \sum _ { i = 1 } ^ { 4 } \phi _ { i } \left( \xi ^ { 1 } , \xi ^ { 2 } \right) { } ^ { t } \mathbf { x } _ { i } + \frac { \xi ^ { 3 } } { 2 } \sum _ { i = 1 } ^ { 4 } h _ { i } \phi _ { i } \left( \xi ^ { 1 } , \xi ^ { 2 } \right) { } ^ { t } \mathbf { V } _ { n } ^ { i }

{ } ^ { t + \Delta t } \mathbf { x } = \sum _ { i = 1 } ^ { 4 } \phi _ { i } \left( \xi ^ { 1 } , \xi ^ { 2 } \right) ^ { t + \Delta t } \mathbf { x } _ { i } + \frac { \xi ^ { 3 } } { 2 } \sum _ { i = 1 } ^ { 4 } h _ { i } \phi _ { i } \left( \xi ^ { 1 } , \xi ^ { 2 } \right) ^ { t + \Delta t } \mathbf { V } _ { n } ^ { i }

시간이 t일 때와 t+\Delta t 일 때 변위 벡터는 다음과 같이 계산 할 수 있다.

\begin{array}{l} { } ^ { t } \mathbf { u } = { } ^ { t } \mathbf { x } - { } ^ { 0 } \mathbf { X } \\ = \sum_ {i = 1} ^ {4} \phi_ {i} (\xi^ {1}, \xi^ {2}) ^ {t} \mathbf {x} _ {i} + \frac {\xi^ {3}}{2} \sum_ {i = 1} ^ {4} h _ {i} \phi_ {i} (\xi^ {1}, \xi^ {2}) ^ {t} \mathbf {V} _ {n} ^ {i} - \sum_ {i = 1} ^ {4} \phi_ {i} (\xi^ {1}, \xi^ {2}) ^ {0} \mathbf {X} _ {i} - \frac {\xi^ {3}}{2} \sum_ {i = 1} ^ {4} h _ {i} \phi_ {i} (\xi^ {1}, \xi^ {2}) ^ {0} \mathbf {V} _ {n} ^ {i} \\ = \sum_ {i = 1} ^ {4} \phi_ {i} \left(^ {t} \mathbf {x} _ {i} - ^ {0} \mathbf {X} _ {i}\right) + \frac {\xi^ {3}}{2} \sum_ {i = 1} ^ {4} h _ {i} \phi_ {i} \left(^ {t} \mathbf {V} _ {n} ^ {i} - ^ {0} \mathbf {V} _ {n} ^ {i}\right) \\ { } ^ { t + \Delta t } \mathbf { u } = { } ^ { t + \Delta t } \mathbf { x } - { } ^ { 0 } \mathbf { X } \\ = \sum_ {i = 1} ^ {4} \phi_ {i} (\xi^ {1}, \xi^ {2}) ^ {t + \Delta t} \mathbf {x} _ {i} + \frac {\xi^ {3}}{2} \sum_ {i = 1} ^ {4} h _ {i} \phi_ {i} (\xi^ {1}, \xi^ {2}) ^ {t + \Delta t} \mathbf {V} _ {n} ^ {i} - \sum_ {i = 1} ^ {4} \phi_ {i} (\xi^ {1}, \xi^ {2}) ^ {0} \mathbf {X} _ {i} - \frac {\xi^ {3}}{2} \sum_ {i = 1} ^ {4} h _ {i} \phi_ {i} (\xi^ {1}, \xi^ {2}) ^ {0} \mathbf {V} _ {n} ^ {i} \\ = \sum_ {i = 1} ^ {4} \phi_ {i} \left(^ {t + \Delta t} \mathbf {x} _ {i} - ^ {0} \mathbf {X} _ {i}\right) + \frac {\xi^ {3}}{2} \sum_ {i = 1} ^ {4} h _ {i} \phi_ {i} \left(^ {t + \Delta t} \mathbf {V} _ {n} ^ {i} - ^ {0} \mathbf {V} _ {n} ^ {i}\right) \\ \end{array}

\begin{array}{l} { } ^ { t + \Delta t } \mathbf { u } = { } ^ { t + \Delta t } \mathbf { x } - { } ^ { 0 } \mathbf { X } \\ = \sum_ {i = 1} ^ {4} \phi_ {i} (\xi^ {1}, \xi^ {2}) ^ {t + \Delta t} \mathbf {x} _ {i} + \frac {\xi^ {3}}{2} \sum_ {i = 1} ^ {4} h _ {i} \phi_ {i} (\xi^ {1}, \xi^ {2}) ^ {t + \Delta t} \mathbf {V} _ {n} ^ {i} - \sum_ {i = 1} ^ {4} \phi_ {i} (\xi^ {1}, \xi^ {2}) ^ {0} \mathbf {X} _ {i} - \frac {\xi^ {3}}{2} \sum_ {i = 1} ^ {4} h _ {i} \phi_ {i} (\xi^ {1}, \xi^ {2}) ^ {0} \mathbf {V} _ {n} ^ {i} \\ = \sum_ {i = 1} ^ {4} \phi_ {i} \left(^ {t + \Delta t} \mathbf {x} _ {i} - ^ {0} \mathbf {X} _ {i}\right) + \frac {\xi^ {3}}{2} \sum_ {i = 1} ^ {4} h _ {i} \phi_ {i} \left(^ {t + \Delta t} \mathbf {V} _ {n} ^ {i} - ^ {0} \mathbf {V} _ {n} ^ {i}\right) \\ \end{array}

따라서 시간 t와 t + \Delta t 사이의 incremental displacement는 다음과 같이 나타낼 수 있다.

\begin{array}{l} \Delta^ {t} \mathbf {u} = ^ {t + \Delta t} \mathbf {u} - ^ {t} \mathbf {u} \\ = \sum_ {i = 1} ^ {4} \phi_ {i} \left(^ {t + \Delta t} \mathbf {x} _ {i} - ^ {0} \mathbf {X} _ {i}\right) + \frac {\xi^ {3}}{2} \sum_ {i = 1} ^ {4} h _ {i} \phi_ {i} \left(^ {t + \Delta t} \mathbf {V} _ {n} ^ {i} - ^ {0} \mathbf {V} _ {n} ^ {i}\right) - \sum_ {i = 1} ^ {4} \phi_ {i} \left(^ {t} \mathbf {x} _ {i} - ^ {0} \mathbf {X} _ {i}\right) - \frac {\xi^ {3}}{2} \sum_ {i = 1} ^ {4} h _ {i} \phi_ {i} \left(^ {t} \mathbf {V} _ {n} ^ {i} - ^ {0} \mathbf {V} _ {n} ^ {i}\right) \\ = \sum_ {i = 1} ^ {4} \phi_ {i} \left(^ {t + \Delta t} \mathbf {x} _ {i} - ^ {t} \mathbf {x} _ {i}\right) + \frac {\xi^ {3}}{2} \sum_ {i = 1} ^ {4} h _ {i} \phi_ {i} \left(^ {t + \Delta t} \mathbf {V} _ {n} ^ {i} - ^ {t} \mathbf {V} _ {n} ^ {i}\right) \\ = \sum_ {i = 1} ^ {4} \phi_ {i} \Delta^ {t} \mathbf {u} _ {i} + \frac {\xi^ {3}}{2} \sum_ {i = 1} ^ {4} h _ {i} \phi_ {i} \left(- \alpha_ {1} ^ {i} {} ^ {t} \mathbf {V} _ {2} ^ {i} + \alpha_ {2} ^ {i} {} ^ {t} \mathbf {V} _ {1} ^ {i}\right) \\ = \sum_ {i = 1} ^ {4} \phi_ {i} \Delta^ {t} \mathbf {u} _ {i} + \frac {\xi^ {3}}{2} \sum_ {i = 1} ^ {4} h _ {i} \phi_ {i} \Delta^ {t} \mathbf {V} _ {n} ^ {i} \\ \end{array}

위치 벡터와 변위 벡터를 matrix form으로 나타내면

^ 0 \mathbf {X} = ^ {0} \mathbf {N} ^ {0} \mathbf {X} _ {n}

= \left[ \begin{array}{l l l l} \phi_ {1} & \phi_ {2} & \phi_ {3} & \phi_ {4} \end{array} \right] \left[ \begin{array}{l} ^ {0} \mathbf {X} _ {1} \\ ^ {0} \mathbf {X} _ {2} \\ ^ {0} \mathbf {X} _ {3} \\ ^ {0} \mathbf {X} _ {4} \end{array} \right] + \left[ \begin{array}{c c c c} \frac {\xi^ {3}}{2} h _ {1} \phi_ {1} & \frac {\xi^ {3}}{2} h _ {2} \phi_ {2} & \frac {\xi^ {3}}{2} h _ {3} \phi_ {3} & \frac {\xi^ {3}}{2} h _ {4} \phi_ {4} \end{array} \right] \left[ \begin{array}{l} ^ {0} \mathbf {V} _ {n} ^ {1} \\ ^ {0} \mathbf {V} _ {n} ^ {2} \\ ^ {0} \mathbf {V} _ {n} ^ {3} \\ ^ {0} \mathbf {V} _ {n} ^ {4} \end{array} \right]

= \left[ \begin{array}{c c c c c c c c} \phi_ {1} & \frac {\xi^ {3}}{2} h _ {1} \phi_ {1} & \phi_ {2} & \frac {\xi^ {3}}{2} h _ {2} \phi_ {2} & \phi_ {3} & \frac {\xi^ {3}}{2} h _ {3} \phi_ {3} & \phi_ {4} & \frac {\xi^ {3}}{2} h _ {4} \phi_ {4} \end{array} \right] \left[ \begin{array}{c} ^ 0 \mathbf {X} _ {1} \\ ^ 0 \mathbf {V} _ {n} ^ {1} \\ ^ 0 \mathbf {X} _ {2} \\ ^ 0 \mathbf {V} _ {n} ^ {2} \\ ^ 0 \mathbf {X} _ {3} \\ ^ 0 \mathbf {V} _ {n} ^ {3} \\ ^ 0 \mathbf {X} _ {4} \\ ^ 0 \mathbf {V} _ {n} ^ {4} \end{array} \right]

{ } ^ { t } \mathbf { x } = { } ^ { 0 } \mathbf { N } ^ { t } \mathbf { x } _ { n }

= \left[ \begin{array}{c c c c c c c c c} \phi_ {1} & \frac {\xi^ {3}}{2} h _ {1} \phi_ {1} & \phi_ {2} & \frac {\xi^ {3}}{2} h _ {2} \phi_ {2} & \phi_ {3} & \frac {\xi^ {3}}{2} h _ {3} \phi_ {3} & \phi_ {4} & \frac {\xi^ {3}}{2} h _ {4} \phi_ {4} \end{array} \right] \left[ \begin{array}{c} ^ {t} \mathbf {x} _ {1} \\ ^ {t} \mathbf {V} _ {n} ^ {1} \\ ^ {t} \mathbf {x} _ {2} \\ ^ {t} \mathbf {V} _ {n} ^ {2} \\ ^ {t} \mathbf {x} _ {3} \\ ^ {t} \mathbf {V} _ {n} ^ {3} \\ ^ {t} \mathbf {x} _ {4} \\ ^ {t} \mathbf {V} _ {n} ^ {4} \end{array} \right]

\mathbf {\Lambda} ^ {t + \Delta t} \mathbf {x} = \mathbf {\Lambda} ^ {0} \mathbf {N} ^ {t + \Delta t} \mathbf {x} _ {n}

= \left[ \begin{array}{c c c c c c c c c} \phi_ {1} & \frac {\xi^ {3}}{2} h _ {1} \phi_ {1} & \phi_ {2} & \frac {\xi^ {3}}{2} h _ {2} \phi_ {2} & \phi_ {3} & \frac {\xi^ {3}}{2} h _ {3} \phi_ {3} & \phi_ {4} & \frac {\xi^ {3}}{2} h _ {4} \phi_ {4} \end{array} \right] \left[ \begin{array}{c} t + \Delta t \mathbf {X} _ {1} \\ t + \Delta t \mathbf {V} _ {n} ^ {1} \\ t + \Delta t \mathbf {X} _ {2} \\ t + \Delta t \mathbf {V} _ {n} ^ {2} \\ t + \Delta t \mathbf {X} _ {3} \\ t + \Delta t \mathbf {V} _ {n} ^ {3} \\ t + \Delta t \mathbf {X} _ {4} \\ t + \Delta t \mathbf {V} _ {n} ^ {4} \end{array} \right]

\Delta^ {t} \mathbf {x} = ^ {0} \mathbf {N} \Delta^ {t} \mathbf {x} _ {n}

= \left[ \begin{array}{c c c c c c c c} \phi_ {1} & \frac {\xi^ {3}}{2} h _ {1} \phi_ {1} & \phi_ {2} & \frac {\xi^ {3}}{2} h _ {2} \phi_ {2} & \phi_ {3} & \frac {\xi^ {3}}{2} h _ {3} \phi_ {3} & \phi_ {4} & \frac {\xi^ {3}}{2} h _ {4} \phi_ {4} \end{array} \right] \left[ \begin{array}{l} \Delta^ {t} \mathbf {x} _ {1} \\ \Delta^ {t} \mathbf {V} _ {n} ^ {1} \\ \Delta^ {t} \mathbf {x} _ {2} \\ \Delta^ {t} \mathbf {V} _ {n} ^ {2} \\ \Delta^ {t} \mathbf {x} _ {3} \\ \Delta^ {t} \mathbf {V} _ {n} ^ {3} \\ \Delta^ {t} \mathbf {x} _ {4} \\ \Delta^ {t} \mathbf {V} _ {n} ^ {4} \end{array} \right]

변위 벡터도 마찬가지로 나타낼 수 있다.

{ } ^ { t } \mathbf { u } = ^ { 0 } \mathbf { N } \left( { } ^ { t } \mathbf { x } _ { n } - ^ { 0 } \mathbf { X } _ { n } \right)

= \left[ \begin{array}{c c c c c c c c} \phi_ {1} & \frac {\xi^ {3}}{2} h _ {1} \phi_ {1} & \phi_ {2} & \frac {\xi^ {3}}{2} h _ {2} \phi_ {2} & \phi_ {3} & \frac {\xi^ {3}}{2} h _ {3} \phi_ {3} & \phi_ {4} & \frac {\xi^ {3}}{2} h _ {4} \phi_ {4} \end{array} \right] \left[ \begin{array}{c} ^ {t} \mathbf {x} _ {1} - ^ {0} \mathbf {X} _ {1} \\ ^ {t} \mathbf {V} _ {n} ^ {1} - ^ {0} \mathbf {V} _ {n} ^ {1} \\ ^ {t} \mathbf {x} _ {2} - ^ {0} \mathbf {X} _ {2} \\ ^ {t} \mathbf {V} _ {n} ^ {2} - ^ {0} \mathbf {V} _ {n} ^ {2} \\ ^ {t} \mathbf {x} _ {3} - ^ {0} \mathbf {X} _ {3} \\ ^ {t} \mathbf {V} _ {n} ^ {3} - ^ {0} \mathbf {V} _ {n} ^ {3} \\ ^ {t} \mathbf {x} _ {4} - ^ {0} \mathbf {X} _ {4} \\ ^ {t} \mathbf {V} _ {n} ^ {4} - ^ {0} \mathbf {V} _ {n} ^ {4} \end{array} \right]

t + \Delta t _ {\mathbf {u}}

= \left[ \begin{array}{c c c c} \phi_ {1} & \phi_ {2} & \phi_ {3} & \phi_ {4} \end{array} \right] \left[ \begin{array}{c} ^ {t + \Delta t} \mathbf {x} _ {1} - ^ {0} \mathbf {X} _ {1} \\ ^ {t + \Delta t} \mathbf {x} _ {2} - ^ {0} \mathbf {X} _ {2} \\ ^ {t + \Delta t} \mathbf {x} _ {3} - ^ {0} \mathbf {X} _ {3} \\ ^ {t + \Delta t} \mathbf {x} _ {4} - ^ {0} \mathbf {X} _ {4} \end{array} \right] + \left[ \begin{array}{c c c c} \frac {\xi^ {3}}{2} h _ {1} \phi_ {1} & \frac {\xi^ {3}}{2} h _ {2} \phi_ {2} & \frac {\xi^ {3}}{2} h _ {3} \phi_ {3} & \frac {\xi^ {3}}{2} h _ {4} \phi_ {4} \end{array} \right] \left[ \begin{array}{c} ^ {t + \Delta t} \mathbf {V} _ {n} ^ {1} - ^ {0} \mathbf {V} _ {n} ^ {1} \\ ^ {t + \Delta t} \mathbf {V} _ {n} ^ {2} - ^ {0} \mathbf {V} _ {n} ^ {2} \\ ^ {t + \Delta t} \mathbf {V} _ {n} ^ {3} - ^ {0} \mathbf {V} _ {n} ^ {3} \\ ^ {t + \Delta t} \mathbf {V} _ {n} ^ {4} - ^ {0} \mathbf {V} _ {n} ^ {4} \end{array} \right]

= \left[ \begin{array}{c c c c} \phi_ {1} & \phi_ {2} & \phi_ {3} & \phi_ {4} \end{array} \right] \left[ \begin{array}{c} ^ {t + \Delta t} \mathbf {u} _ {1} \\ ^ {t + \Delta t} \mathbf {u} _ {2} \\ ^ {t + \Delta t} \mathbf {u} _ {3} \\ ^ {t + \Delta t} \mathbf {u} _ {4} \end{array} \right] + \left[ \begin{array}{c c c c} \frac {\xi^ {3}}{2} h _ {1} \phi_ {1} & \frac {\xi^ {3}}{2} h _ {2} \phi_ {2} & \frac {\xi^ {3}}{2} h _ {3} \phi_ {3} & \frac {\xi^ {3}}{2} h _ {4} \phi_ {4} \end{array} \right] \left[ \begin{array}{c} \Delta^ {t} \mathbf {V} _ {n} ^ {1} + \Delta^ {0} \mathbf {V} _ {n} ^ {1} \\ \Delta^ {t} \mathbf {V} _ {n} ^ {2} + \Delta^ {0} \mathbf {V} _ {n} ^ {2} \\ \Delta^ {t} \mathbf {V} _ {n} ^ {3} + \Delta^ {0} \mathbf {V} _ {n} ^ {3} \\ \Delta^ {t} \mathbf {V} _ {n} ^ {4} + \Delta^ {0} \mathbf {V} _ {n} ^ {4} \end{array} \right]

\left(^ {t + \Delta t} \mathbf {V} _ {n} ^ {1} - ^ {0} \mathbf {V} _ {n} ^ {1} = ^ {t + \Delta t} \mathbf {V} _ {n} ^ {1} - ^ {t} \mathbf {V} _ {n} ^ {1} + ^ {t} \mathbf {V} _ {n} ^ {1} - ^ {0} \mathbf {V} _ {n} ^ {1} = \Delta^ {t} \mathbf {V} _ {n} ^ {1} + \Delta^ {0} \mathbf {V} _ {n} ^ {1}\right)

= \left[ \begin{array}{l l l l} \phi_ {1} & \phi_ {2} & \phi_ {3} & \phi_ {4} \end{array} \right] \left[ \begin{array}{c} ^ {t + \Delta t} \mathbf {u} _ {1} \\ ^ {t + \Delta t} \mathbf {u} _ {2} \\ ^ {t + \Delta t} \mathbf {u} _ {3} \\ ^ {t + \Delta t} \mathbf {u} _ {4} \end{array} \right] + \left[ \begin{array}{c c c c} \frac {\xi^ {3}}{2} h _ {1} \phi_ {1} & \frac {\xi^ {3}}{2} h _ {2} \phi_ {2} & \frac {\xi^ {3}}{2} h _ {3} \phi_ {3} & \frac {\xi^ {3}}{2} h _ {4} \phi_ {4} \end{array} \right] \left[ \begin{array}{c} - \alpha_ {1} ^ {1 t} \mathbf {V} _ {2} ^ {1} + \alpha_ {2} ^ {1 t} \mathbf {V} _ {1} ^ {1} \\ - \alpha_ {1} ^ {2 t} \mathbf {V} _ {2} ^ {2} + \alpha_ {2} ^ {2 t} \mathbf {V} _ {1} ^ {2} \\ - \alpha_ {1} ^ {3 t} \mathbf {V} _ {2} ^ {3} + \alpha_ {2} ^ {3 t} \mathbf {V} _ {1} ^ {3} \\ - \alpha_ {1} ^ {4 t} \mathbf {V} _ {2} ^ {4} + \alpha_ {2} ^ {4 t} \mathbf {V} _ {1} ^ {4} \end{array} \right]

+ \left[ \begin{array}{c c c c} \frac {\xi^ {3}}{2} h _ {1} \phi_ {1} & \frac {\xi^ {3}}{2} h _ {2} \phi_ {2} & \frac {\xi^ {3}}{2} h _ {3} \phi_ {3} & \frac {\xi^ {3}}{2} h _ {4} \phi_ {4} \end{array} \right] \left[ \begin{array}{c} \Delta^ {0} \mathbf {V} _ {n} ^ {1} \\ \Delta^ {0} \mathbf {V} _ {n} ^ {2} \\ \Delta^ {0} \mathbf {V} _ {n} ^ {3} \\ \Delta^ {0} \mathbf {V} _ {n} ^ {4} \end{array} \right]

\Rightarrow^ {t + \Delta t} \mathbf {u} = ^ {t} \mathbf {N} ^ {t + \Delta t} \mathbf {u} _ {n} + ^ {0} \tilde {\mathbf {N}} \Delta^ {0} \tilde {\mathbf {X}} ^ {n}

= \left[ \begin{array}{c c c c c c c c c c c c c c c} \phi_ {1} & - \frac {\xi^ {3}}{2} h _ {1} \phi_ {1} \mathbf {v} _ {2} ^ {1} & \frac {\xi^ {3}}{2} h _ {1} \phi_ {1} \mathbf {v} _ {1} ^ {1} & \phi_ {2} & - \frac {\xi^ {3}}{2} h _ {2} \phi_ {2} \mathbf {v} _ {2} ^ {2} & \frac {\xi^ {3}}{2} h _ {2} \phi_ {2} \mathbf {v} _ {1} ^ {2} & \phi_ {3} & - \frac {\xi^ {3}}{2} h _ {3} \phi_ {3} \mathbf {v} _ {2} ^ {3} & \frac {\xi^ {3}}{2} h _ {3} \phi_ {3} \mathbf {v} _ {1} ^ {3} & \phi_ {4} & - \frac {\xi^ {3}}{2} h _ {4} \phi_ {4} \mathbf {v} _ {2} ^ {4} & \frac {\xi^ {3}}{2} h _ {4} \phi_ {4} \mathbf {v} _ {1} ^ {4} \end{array} \right] \left[ \begin{array}{c} t + \Delta t \mathbf {u} _ {1} \\ \alpha_ {1} ^ {1} \\ \alpha_ {2} ^ {1} \\ t + \Delta t \mathbf {u} _ {2} \\ \alpha_ {1} ^ {2} \\ \alpha_ {2} ^ {2} \\ t + \Delta t \mathbf {u} _ {3} \\ \alpha_ {1} ^ {3} \\ \alpha_ {2} ^ {3} \\ t + \Delta t \mathbf {u} _ {4} \\ \alpha_ {1} ^ {4} \\ \alpha_ {2} ^ {4} \end{array} \right]

+ \left[ \begin{array}{c c c c} \frac {\xi^ {3}}{2} h _ {1} \phi_ {1} & \frac {\xi^ {3}}{2} h _ {2} \phi_ {2} & \frac {\xi^ {3}}{2} h _ {3} \phi_ {3} & \frac {\xi^ {3}}{2} h _ {4} \phi_ {4} \end{array} \right] \left[ \begin{array}{c} \Delta^ {0} \mathbf {v} _ {n} ^ {1} \\ \Delta^ {0} \mathbf {v} _ {n} ^ {2} \\ \Delta^ {0} \mathbf {v} _ {n} ^ {3} \\ \Delta^ {0} \mathbf {v} _ {n} ^ {4} \end{array} \right]

3. FE Formulation

현재 형상(t+Δt)에서의 평형방정식은 다음과 같다.

\nabla_ {X} \cdot^ {t + \Delta t} \boldsymbol {\sigma} + \rho^ {t + \Delta t} \mathbf {f} = \rho^ {t + \Delta t} \ddot {\mathbf {u}}

가상일 원리를 적용하면

\begin{array}{l} \int_ {V} \delta^ {t + \Delta t} \mathbf {u} \cdot \left(\nabla_ {X} \cdot^ {t + \Delta t} \boldsymbol {\sigma} + \rho^ {t + \Delta t} \mathbf {f}\right) d V = \int_ {V} \delta^ {t + \Delta t} \mathbf {u} \cdot \rho^ {t + \Delta t} \ddot {\mathbf {u}} d V \\ \Rightarrow \int_ {V} \delta u _ {i} \left(\frac {\partial \sigma_ {i j}}{\partial x _ {j}} + \rho f _ {i}\right) d V = \int_ {V} \delta u _ {i} \cdot \rho \ddot {u} _ {i} d V \quad \left(\delta u _ {i} \frac {\partial \sigma_ {i j}}{\partial x _ {j}} = \frac {\partial \left(\delta u _ {i} \sigma_ {i j}\right)}{\partial x _ {j}} - \frac {\partial \delta u _ {i}}{\partial x _ {j}} \sigma_ {i j}\right) \\ \Rightarrow \int_ {V} \left\{\frac {\partial \left(\delta u _ {i} \sigma_ {i j}\right)}{\partial x _ {j}} - \frac {\partial \delta u _ {i}}{\partial x _ {j}} \sigma_ {i j} + \delta u _ {i} \rho f _ {i} \right\} d V = \int_ {V} \delta u _ {i} \cdot \rho \ddot {u} _ {i} d V \quad \rightarrow \text { Gauss's   divergence   theorem } \\ \end{array}

\Rightarrow \int_ {\partial V} \delta u _ {i} \sigma_ {i j} n _ {j} d A + \int_ {V} \left\{- \frac {\partial \delta u _ {i}}{\partial x _ {j}} \sigma_ {i j} + \delta u _ {i} \rho f _ {i} \right\} d V = \int_ {V} \delta u _ {i} \cdot \rho \ddot {u} _ {i} d V \quad \left(\text { geometric   B   .   C   瓦   natural   B   .   C }\right)

\Rightarrow \int_ {\partial V _ {g}} \delta u _ {i} \sigma_ {i j} n _ {j} d A + \int_ {\partial V _ {m}} \delta u _ {i} \bar {t} _ {i} d A + \int_ {V} \left\{- \frac {\partial \delta u _ {i}}{\partial x _ {j}} \sigma_ {i j} + \delta u _ {i} \rho f _ {i} \right\} d V = \int_ {V} \delta u _ {i} \cdot \rho \ddot {u} _ {i} d V

\left(\int_ {\partial V _ {g}} \delta u _ {i} \sigma_ {i j} n _ {j} d A = 0, \text {   기하학적   경계조건에서   가상변위가   } 0\right)

\Rightarrow \int_ {\partial V _ {m}} \delta u _ {i} \bar {t} _ {i} d A + \int_ {V} \left\{- \frac {\partial \delta u _ {i}}{\partial x _ {j}} \sigma_ {i j} + \delta u _ {i} \rho f _ {i} \right\} d V = \int_ {V} \delta u _ {i} \cdot \rho \ddot {u} _ {i} d V

\left(\frac {\partial \delta u _ {i}}{\partial x _ {j}} \sigma_ {i j} = \frac {1}{2} \left(\frac {\partial \delta u _ {i}}{\partial x _ {j}} + \frac {\partial \delta u _ {i}}{\partial x _ {j}}\right) \sigma_ {i j} = \delta \varepsilon_ {i j} \sigma_ {i j}\right)

\Rightarrow \int_ {\partial V _ {m}} \delta u _ {i} \bar {t} _ {i} d A + \int_ {V} \delta u _ {i} \rho f _ {i} d V = \int_ {V} \delta u _ {i} \cdot \rho \ddot {u} _ {i} d V + \int_ {V} \delta \varepsilon_ {i j} \sigma_ {i j} d V

비선형 해석을 위해 하중을 조금씩 증가시켜 물체의 변형을 순차적으로 구해 나가며 이러한 반복계산에 있어 기준이 되는 물체의 형상을 설정하는 방법에는 크게 Total Lagrange formulation과 Updated Lagrange formulation이 있다. Total Lagrange formulation은 변형과 관련된 변수들을 초기형상을 이용하여 정의하고 Updated Lagrange formulation은 변수들을 현재 변형된 형상으로 정의한다.

\int_ {V} \delta u _ {i} \cdot \rho \ddot {u} _ {i} d V + \int_ {V} \delta \varepsilon_ {i j} \sigma_ {i j} d V = \int_ {\partial V _ {m}} \delta u _ {i} \overline {{t}} _ {i} d A + \int_ {V} \delta u _ {i} \rho f _ {i} d V

\left( \begin{array}{l} \int_ {V} (\bullet) d V = \int_ {V _ {0}} (\bullet) \det (\mathbf {F}) d V _ {0} = \int_ {V _ {0}} (\bullet) J d V _ {0} \\ \int_ {\partial V} (\bullet) \mathbf {n} d A = \int_ {\partial V} (\bullet) J \mathbf {F} ^ {- T} \tilde {\mathbf {n}} d A (N a n s o n ^ {\prime} s f o r m u l a) \\ \int_ {V _ {0}} \delta \varepsilon_ {i j} \sigma_ {i j} d V _ {0} = \int_ {V _ {0}} \delta F _ {i j} P _ {i j} d V _ {0} = \int_ {V _ {0}} \delta E _ {i j} S _ {i j} d V _ {0} \end{array} \right)

\Rightarrow \int_ {V _ {0}} \delta u _ {i} \cdot \rho J \ddot {u} _ {i} d V _ {0} + \int_ {V _ {0}} \delta E _ {i j} S _ {i j} d V _ {0} = \int_ {\partial V _ {0 m}} \delta u _ {i} \sigma_ {i j} J \left[ F ^ {- T} \right] _ {k j} \tilde {n} _ {k} d A _ {0} + \int_ {V _ {0}} \delta u _ {i} \rho J f _ {i} d V _ {0}

\Rightarrow \int_ {V _ {0}} \delta u _ {i} \cdot \rho_ {0} \ddot {u} _ {i} d V _ {0} + \int_ {V _ {0}} \delta E _ {i j} S _ {i j} d V _ {0} = \int_ {\partial V _ {0 m}} \delta u _ {i} \sigma_ {i j} J \left[ F ^ {- T} \right] _ {k j} \tilde {n} _ {k} d A _ {0} + \int_ {V _ {0}} \delta u _ {i} \rho_ {0} f _ {i} d V _ {0}

\Rightarrow \int_ {V _ {0}} \delta^ {t + \Delta t} \mathbf {u} \cdot \rho_ {0} ^ {t + \Delta t} \ddot {\mathbf {u}} d V _ {0} + \int_ {V _ {0}} \delta^ {t + \Delta t} _ {0} \mathbf {E} \cdot {} _ {0} ^ {t + \Delta t} \mathbf {S} d V _ {0} = \int_ {\partial V _ {0 m}} \delta^ {t + \Delta t} \mathbf {u} J ^ {t + \Delta t} \boldsymbol {\sigma} F ^ {- T t + \Delta t} \tilde {\mathbf {n}} d A _ {0} + \int_ {V _ {0}} \delta^ {t + \Delta t} \mathbf {u} \rho_ {0} ^ {t + \Delta t} \mathbf {f} d V _ {0}

\left(J \boldsymbol {\sigma} F ^ {- T} = \mathbf {P} = \mathbf {F S}\right)

위 식을 정리하면 다음과 같다.

\int_ {V _ {0}} \delta^ {t + \Delta t} \mathbf {u} \cdot \rho_ {0} ^ {t + \Delta t} \ddot {\mathbf {u}} d V _ {0} + \int_ {V _ {0}} \delta^ {t + \Delta t} _ {0} \mathbf {E} \cdot {} _ {0} ^ {t + \Delta t} \mathbf {S} d V _ {0} = \int_ {\partial V _ {0 m}} \delta^ {t + \Delta t} \mathbf {u} ^ {t + \Delta t} _ {0} \mathbf {F} ^ {t + \Delta t} _ {0} \mathbf {S} ^ {t + \Delta t} \tilde {\mathbf {n}} d A _ {0} + \int_ {V _ {0}} \delta^ {t + \Delta t} \mathbf {u} \rho_ {0} ^ {t + \Delta t} \mathbf {f} d V _ {0}

\Leftrightarrow \int_ {V _ {0}} \delta^ {t + \Delta t} \mathbf {u} \cdot \rho_ {0} ^ {t + \Delta t} \ddot {\mathbf {u}} d V _ {0} + \int_ {V _ {0}} \delta_ {0} ^ {t + \Delta t} \mathbf {F}: _ {0} ^ {t + \Delta t} \mathbf {P} d V _ {0} = \int_ {\partial V _ {0 m}} \delta^ {t + \Delta t} \mathbf {u} _ {0} ^ {t + \Delta t} \mathbf {P} ^ {t + \Delta t} \tilde {\mathbf {n}} d A _ {0} + \int_ {V _ {0}} \delta^ {t + \Delta t} \mathbf {u} \rho_ {0} ^ {t + \Delta t} \mathbf {f} d V _ {0}

E 는 Green-Lagrange strain, S 는 2 ^{nd} PK stress, P 는 1 ^{st} PK stress, F 는 deformation gradient를 나타낸다. 여기에서는 Green-Lagrange strain과 2 ^{nd} PK stress를 사용한다.

\int_ {V _ {0}} \delta^ {t + \Delta t} \mathbf {u} \cdot \rho_ {0} ^ {t + \Delta t} \ddot {\mathbf {u}} d V _ {0} + \int_ {V _ {0}} \delta_ {0} ^ {t + \Delta t} \mathbf {E} \cdot {} _ {0} ^ {t + \Delta t} \mathbf {S} d V _ {0} = \int_ {\partial V _ {0 m}} \delta^ {t + \Delta t} \mathbf {u} _ {0} ^ {t + \Delta t} \mathbf {F} _ {0} ^ {t + \Delta t} \mathbf {S} ^ {t + \Delta t} \tilde {\mathbf {n}} d A _ {0} + \int_ {V _ {0}} \delta^ {t + \Delta t} \mathbf {u} \rho_ {0} ^ {t + \Delta t} \mathbf {f} d V _ {0}

먼저 좌변을 살펴보면 첫 번째 항을 다음과 같이 정리 할 수 있다.

\int_ {V _ {0}} \delta^ {t + \Delta t} \mathbf {u} \cdot \rho^ {t + \Delta t} \ddot {\mathbf {u}} J d V _ {0} = \int_ {V _ {0}} \delta \left(^ {t} \mathbf {N} ^ {t + \Delta t} \mathbf {u} _ {n} + ^ {0} \tilde {\mathbf {N}} \Delta^ {0} \tilde {\mathbf {X}} _ {n}\right) \cdot \rho^ {t} \mathbf {N} ^ {t + \Delta t} \ddot {\mathbf {u}} J d V _ {0}

= \int_ {V _ {0}} \delta \left(^ {t} \mathbf {N} ^ {t + \Delta t} \mathbf {u} _ {n}\right) \cdot \rho \left(^ {t} \mathbf {N} ^ {t + \Delta t} \ddot {\mathbf {u}}\right) J d V _ {0}

= \left[ \delta^ {t + \Delta t} \mathbf {u} _ {n} \right] ^ {T} \int_ {V _ {0}} \rho \left[ ^ {t} \mathbf {N} \right] ^ {T} \left[ ^ {t} \mathbf {N} \right] J d V _ {0} \left[ ^ {t + \Delta t} \ddot {\mathbf {u}} _ {n} \right]

두 번째 항도 마찬가지로 다음과 같이 정리 할 수 있다. 먼저 Green-Lagrange strain은 다음과 같이 나타낼 수 있다.

\begin{array}{l} { } _ { 0 } ^ { t + \Delta t } \mathbf { E } = \frac { 1 } { 2 } \Big ( { } _ { 0 } ^ { t + \Delta t } \mathbf { g } _ { i } \cdot { } _ { 0 } ^ { t + \Delta t } \mathbf { g } _ { j } - { } _ { 0 } \mathbf { G } _ { i } \cdot { } _ { 0 } \mathbf { G } _ { j } \Big ) \Big ( { } _ { 0 } \mathbf { G } ^ { i } \otimes _ { 0 } \mathbf { G } ^ { j } \Big ) \\ = \frac {1}{2} \left(\frac {\partial_ {0} ^ {t + \Delta t} \mathbf {x}}{\partial \xi^ {i}} \cdot \frac {\partial_ {0} ^ {t + \Delta t} \mathbf {x}}{\partial \xi^ {j}} - \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {j}} \cdot \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {j}}\right) \left(_ {0} \mathbf {G} ^ {i} \otimes_ {0} \mathbf {G} ^ {j}\right) \\ = \frac {1}{2} \left(\frac {\partial \left(^ {0} \mathbf {X} + ^ {t + \Delta t} \mathbf {u}\right)}{\partial \xi^ {i}} \cdot \frac {\partial \left(^ {0} \mathbf {X} + ^ {t + \Delta t} \mathbf {u}\right)}{\partial \xi^ {j}} - \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {j}} \cdot \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {j}}\right) \left(_ {0} \mathbf {G} ^ {i} \otimes_ {0} \mathbf {G} ^ {j}\right) \\ = \frac {1}{2} \left(\frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {j}} \cdot \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {j}} + \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {i}} \cdot \frac {\partial^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {j}} + \frac {\partial^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {j}} + \frac {\partial^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {j}} - \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {j}} \cdot \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {j}}\right) \left(_ {0} \mathbf {G} ^ {i} \otimes_ {0} \mathbf {G} ^ {j}\right) \\ = \frac {1}{2} \left(\frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {i}} \cdot \frac {\partial \left(^ {t} \mathbf {u} + \Delta^ {t} \mathbf {u}\right)}{\partial \xi^ {j}} + \frac {\partial \left(^ {t} \mathbf {u} + \Delta^ {t} \mathbf {u}\right)}{\partial \xi^ {i}} \cdot \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {j}} + \frac {\partial \left(^ {t} \mathbf {u} + \Delta^ {t} \mathbf {u}\right)}{\partial \xi^ {i}} \cdot \frac {\partial \left(^ {t} \mathbf {u} + \Delta^ {t} \mathbf {u}\right)}{\partial \xi^ {j}}\right) \left(_ {0} \mathbf {G} ^ {i} \otimes_ {0} \mathbf {G} ^ {j}\right) \\ \end{array}

정리하면 Green-Lagrange strain을 Δu에 대한 상수 term, 선형 term, 비선형 term으로 나눌 수 있다.

{ } _ { 0 } ^ { t + \Delta t } \mathbf { E } = { } _ { 0 } ^ { t + \Delta t } \mathbf { E } _ { 0 } + { } _ { 0 } ^ { t + \Delta t } \mathbf { E } _ { C } + { } _ { 0 } ^ { t + \Delta t } \mathbf { E } _ { L }

\left( \begin{array}{l} ^ {t + \Delta t} _ {0} \mathbf {E} _ {0} = \frac {1}{2} \left(\frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {i}} \cdot \frac {\partial^ {t} \mathbf {u}}{\partial \xi^ {j}} + \frac {\partial^ {t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {j}} + \frac {\partial^ {t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial^ {t} \mathbf {u}}{\partial \xi^ {j}}\right) \left(_ {0} \mathbf {G} ^ {i} \otimes_ {0} \mathbf {G} ^ {j}\right) \\ ^ {t + \Delta t} _ {0} \mathbf {E} _ {C} = \frac {1}{2} \left(\frac {\partial^ {t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial \Delta^ {t} \mathbf {u}}{\partial \xi^ {j}} + \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {i}} \cdot \frac {\partial \Delta^ {t} \mathbf {u}}{\partial \xi^ {j}} + \frac {\partial \Delta^ {t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {j}} + \frac {\partial \Delta^ {t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial^ {t} \mathbf {u}}{\partial \xi^ {j}}\right) \left(_ {0} \mathbf {G} ^ {i} \otimes_ {0} \mathbf {G} ^ {j}\right) \\ ^ {t + \Delta t} _ {0} \mathbf {E} ^ {L} = \frac {1}{2} \left(\frac {\partial \Delta^ {t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial \Delta^ {t} \mathbf {u}}{\partial \xi^ {j}}\right) \left(_ {0} \mathbf {G} ^ {i} \otimes_ {0} \mathbf {G} ^ {j}\right) \end{array} \right)

와 같다. 두 번째 항을 위의 표현으로 나타내면

\int_ {V _ {0}} \delta^ {t + \Delta t} _ {0} \mathbf {E}: ^ {t + \Delta t} _ {0} \mathbf {S} d V _ {0} = \int_ {V _ {0}} \delta^ {t + \Delta t} _ {0} \mathbf {E}: ^ {t + \Delta t} _ {0} \mathbf {C}: ^ {t + \Delta t} _ {0} \mathbf {E} d V _ {0}

= \int_ {V _ {0}} \left(\underbrace {\delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {0}} _ {0} + \delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {C} + \delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {L}\right): _ {0} ^ {t + \Delta t} \mathbf {C}: \left(_ {0} ^ {t + \Delta t} \mathbf {E} _ {0} + _ {0} ^ {t + \Delta t} \mathbf {E} _ {C} + _ {0} ^ {t + \Delta t} \mathbf {E} _ {L}\right) d V _ {0}

= \int_ {V _ {0}} \underbrace {\delta^ {t + \Delta t} \mathbf {E} _ {C} : _ {0} ^ {t + \Delta t} \mathbf {C} : _ {0} ^ {t + \Delta t} \mathbf {E} _ {0}} _ {\text {constant term}} d V _ {0} + \int_ {V _ {0}} \underbrace {\delta^ {t + \Delta t} \mathbf {E} _ {C} : _ {0} ^ {t + \Delta t} \mathbf {C} : _ {0} ^ {t + \Delta t} \mathbf {E} _ {C} + \delta^ {t + \Delta t} \mathbf {E} _ {L} : _ {0} ^ {t + \Delta t} \mathbf {C} : _ {0} ^ {t + \Delta t} \mathbf {E} _ {0}} _ {\text {linear term}} d V _ {0}

+ \int_ {V _ {0}} \underbrace {\delta^ {t + \Delta t} \mathbf {E} _ {C} : {} _ {0} ^ {t + \Delta t} \mathbf {C} : {} _ {0} ^ {t + \Delta t} \mathbf {E} _ {L} + \delta^ {t + \Delta t} \mathbf {E} _ {L} : {} _ {0} ^ {t + \Delta t} \mathbf {C} : {} _ {0} ^ {t + \Delta t} \mathbf {E} _ {C} + \delta^ {t + \Delta t} \mathbf {E} _ {L} : {} _ {0} ^ {t + \Delta t} \mathbf {C} : {} _ {0} ^ {t + \Delta t} \mathbf {E} _ {L}} _ {\text {high order term}} d V _ {0}

와 같다. 이후 선형화 시키고(High order term 무시) component형태로 나타내면 아래와 같다.

\int_ {V _ {0}} \underbrace {\delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {C} : _ {0} ^ {t + \Delta t} \mathbf {C} : _ {0} ^ {t + \Delta t} \mathbf {E} _ {0}} _ {\text {constant term}} d V _ {0} + \int_ {V _ {0}} \underbrace {\delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {C} : _ {0} ^ {t + \Delta t} \mathbf {C} : _ {0} ^ {t + \Delta t} \mathbf {E} _ {C} + \delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {L} : _ {0} ^ {t + \Delta t} \mathbf {C} : _ {0} ^ {t + \Delta t} \mathbf {E} _ {0}} _ {\text {linear term}} d V _ {0}

= \int_ {V _ {0}} \underbrace {\left[ \delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {C} \right] _ {i j} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {C} \right] ^ {i j k l} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {E} _ {0} \right] _ {k l}} _ {\text { constant   term }} d V _ {0}

+ \int_ {V _ {0}} \underbrace {\left[ \delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {C} \right] _ {i j} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {C} \right] ^ {i j k l} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {E} _ {C} \right] _ {k l} + \left[ \delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {L} \right] _ {i j} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {C} \right] ^ {i j k l} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {E} _ {0} \right] _ {k l}} _ {\text {linear term}} d V _ {0}

= \int_ {V _ {0}} \underbrace {\left[ \delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {C} \right] _ {i j} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {S} _ {0} \right] ^ {i j}} _ {\text {constant term}} d V _ {0} + \int_ {V _ {0}} \underbrace {\left[ \delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {C} \right] _ {i j} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {S} _ {C} \right] ^ {i j} + \left[ \delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {L} \right] _ {i j} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {S} _ {0} \right] ^ {i j}} _ {\text {linear term}} d V _ {0}

Green-Lagrange strain의 변분을 구해보면 다음과 같다. 변분은 incremental displacement에만 적용되는데 이는 incremental displacement가 정의되지 않았기 때문이다. 따라서 \delta^{t+\Delta t}\mathbf{u}=\delta\left(^{t}\mathbf{u}+\Delta^{t}\mathbf{u}\right)=\delta\Delta^{t}\mathbf{u} 이고 \delta^{t+\Delta t}_{0}\mathbf{E}_{0}=0 이 성립한다.

\delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {C} = \frac {1}{2} \left(\frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {i}} \cdot \frac {\partial \delta^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {j}} + \frac {\partial \delta^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {j}} + \frac {\partial \delta^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {j}} + \frac {\partial^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial \delta^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {j}}\right) \left(_ {0} \mathbf {G} ^ {i} \otimes_ {0} \mathbf {G} ^ {j}\right)

\delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {L} = \frac {1}{2} \left(\frac {\partial \delta^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial \Delta_ {0} \mathbf {u}}{\partial \xi^ {j}} + \frac {\partial \Delta_ {0} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial \delta^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {j}}\right) \left(_ {0} \mathbf {G} ^ {i} \otimes_ {0} \mathbf {G} ^ {j}\right)

Component form으로 나타내면

\begin{array}{l} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {E} _ {0} \right] _ {i j} = \frac {1}{2} \left(\frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {i}} \cdot \frac {\partial^ {t} \mathbf {u}}{\partial \xi^ {j}} + \frac {\partial^ {t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {j}} + \frac {\partial^ {t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial^ {t} \mathbf {u}}{\partial \xi^ {j}}\right) \\ = \left[ ^ {0} \mathbf {X} _ {n} \right] ^ {T} \frac {1}{2} \left\{\left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {j}} \right] + \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {j}} \right] \right\} \left[ ^ {t} \mathbf {x} _ {n} - ^ {0} \mathbf {X} _ {n} \right] \\ + \frac {1}{2} \left[ ^ {t} \mathbf {x} _ {n} - ^ {0} \mathbf {X} _ {n} \right] ^ {T} \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {j}} \right] \left[ ^ {t} \mathbf {x} _ {n} - ^ {0} \mathbf {X} _ {n} \right] \\ = \left[ ^ {0} \mathbf {X} _ {n} \right] ^ {T} \frac {1}{2} \left\{\left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {j}} \right] + \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {j}} \right] \right\} \left[ ^ {t} \mathbf {x} _ {n} - ^ {0} \mathbf {X} _ {n} \right] \\ + \frac {1}{2} \left[ ^ {t} \mathbf {x} _ {n} - ^ {0} \mathbf {X} _ {n} \right] ^ {T} \frac {1}{2} \left\{\left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {j}} \right] + \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {j}} \right] \right\} \left[ ^ {t} \mathbf {x} _ {n} - ^ {0} \mathbf {X} _ {n} \right] \\ = \frac {1}{2} \left[ ^ {t} \mathbf {x} _ {n} + ^ {0} \mathbf {X} _ {n} \right] ^ {T} \underbrace {\frac {1}{2} \left\{\left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {j}} \right] + \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {j}} \right] \right\}} _ {[ \mathbf {e} ] _ {i j}} \left[ ^ {t} \mathbf {x} _ {n} - ^ {0} \mathbf {X} _ {n} \right] \\ \end{array}

\left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {E} _ {C} \right] _ {i j} = \frac {1}{2} \left(\frac {\partial^ {t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial \Delta^ {t} \mathbf {u}}{\partial \xi^ {j}} + \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {i}} \cdot \frac {\partial \Delta^ {t} \mathbf {u}}{\partial \xi^ {j}} + \frac {\partial \Delta^ {t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {j}} + \frac {\partial \Delta^ {t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial^ {t} \mathbf {u}}{\partial \xi^ {j}}\right)

= \left[ ^ {0} \mathbf {X} _ {n} \right] ^ {T} \frac {1}{2} \left\{\left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {t} \mathbf {N}}{\partial \xi^ {j}} \right] + \left[ \frac {\partial^ {t} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {j}} \right] \right\} \left[ \Delta^ {t} \mathbf {u} _ {n} \right]

+ \left[ ^ {t} \mathbf {x} _ {n} - ^ {0} \mathbf {X} _ {n} \right] ^ {T} \frac {1}{2} \left\{\left[ \frac {\partial^ {t} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {j}} \right] + \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {t} \mathbf {N}}{\partial \xi^ {j}} \right] \right\} \left[ \Delta^ {t} \mathbf {u} _ {n} \right]

= \left[ ^ {t} \mathbf {x} _ {n} \right] ^ {T} \underbrace {\frac {1}{2} \left\{\left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {t} \mathbf {N}}{\partial \xi^ {j}} \right] + \left[ \frac {\partial^ {t} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {j}} \right] \right\}} _ {[ \mathbf {a} ] _ {i j}} \left[ \Delta^ {t} \mathbf {u} _ {n} \right]

\begin{array}{l} \left[ \begin{array}{c} ^ {t + \Delta t} _ {0} \mathbf {E} _ {L} \end{array} \right] _ {i j} = \frac {1}{2} \left(\frac {\partial \Delta^ {t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial \Delta^ {t} \mathbf {u}}{\partial \xi^ {j}}\right) \\ = \left[ \Delta^ {t} \mathbf {u} _ {n} \right] ^ {T} \frac {1}{2} \left[ \frac {\partial^ {t} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {t} \mathbf {N}}{\partial \xi^ {j}} \right] \left[ \Delta^ {t} \mathbf {u} _ {n} \right] \\ \end{array}

\left[ \delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {C} \right] _ {i j} = \frac {1}{2} \left(\frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {i}} \cdot \frac {\partial \delta^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {j}} + \frac {\partial \delta^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {j}} + \frac {\partial \delta^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {j}} + \frac {\partial^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial \delta^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {j}}\right)

= \frac {1}{2} \left( \begin{array}{c} \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {i}} \cdot \frac {\partial \delta \left(^ {t} \mathbf {N} ^ {t + \Delta t} \mathbf {u} _ {n} + ^ {0} \tilde {\mathbf {N}} \Delta^ {0} \tilde {\mathbf {X}} _ {n}\right)}{\partial \xi^ {j}} + \frac {\partial \delta \left(^ {t} \mathbf {N} ^ {t + \Delta t} \mathbf {u} _ {n} + ^ {0} \tilde {\mathbf {N}} \Delta^ {0} \tilde {\mathbf {X}} _ {n}\right)}{\partial \xi^ {i}} \cdot \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {j}} \\ + \frac {\partial \left(^ {t} \mathbf {N} ^ {t + \Delta t} \mathbf {u} _ {n} + ^ {0} \tilde {\mathbf {N}} \Delta^ {0} \tilde {\mathbf {X}} _ {n}\right)}{\partial \xi^ {i}} \cdot \frac {\partial^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {j}} + \frac {\partial^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial \left(^ {t} \mathbf {N} ^ {t + \Delta t} \mathbf {u} _ {n} + ^ {0} \tilde {\mathbf {N}} \Delta^ {0} \tilde {\mathbf {X}} _ {n}\right)}{\partial \xi^ {j}} \end{array} \right)

= \frac {1}{2} \left(\frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {i}} \cdot \frac {\partial \delta^ {t} \mathbf {N} ^ {t + \Delta t} \mathbf {u} _ {n}}{\partial \xi^ {j}} + \frac {\partial \delta^ {t} \mathbf {N} ^ {t + \Delta t} \mathbf {u} _ {n}}{\partial \xi^ {i}} \cdot \frac {\partial^ {0} \mathbf {X}}{\partial \xi^ {j}} + \frac {\partial^ {t} \mathbf {N} ^ {t + \Delta t} \mathbf {u} _ {n}}{\partial \xi^ {i}} \cdot \frac {\partial^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {j}} + \frac {\partial^ {t + \Delta t} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial^ {t} \mathbf {N} ^ {t + \Delta t} \mathbf {u} _ {n}}{\partial \xi^ {j}}\right)

= \left[ \boldsymbol {\delta} ^ {t + \Delta t} \mathbf {u} _ {n} \right] ^ {T} \frac {1}{2} \left\{\left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {t} \mathbf {N}}{\partial \xi^ {j}} \right] + \left[ \frac {\partial^ {t} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {j}} \right] \right\} \left[ ^ {0} \mathbf {X} _ {n} \right]

+ \left[ \delta^ {t + \Delta t} \mathbf {u} _ {n} \right] ^ {T} \frac {1}{2} \left\{\left[ \frac {\partial^ {t} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {j}} \right] + \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {t} \mathbf {N}}{\partial \xi^ {j}} \right] \right\} \left[ ^ {t} \mathbf {x} _ {n} - ^ {0} \mathbf {X} _ {n} \right]

= \left[ \delta^ {t + \Delta t} \mathbf {u} _ {n} \right] ^ {T} \underbrace {\frac {1}{2} \left\{\left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {t} \mathbf {N}}{\partial \xi^ {j}} \right] + \left[ \frac {\partial^ {t} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \left[ \frac {\partial^ {0} \mathbf {N}}{\partial \xi^ {j}} \right] \right\}} _ {[ \mathbf {a} ] _ {i j}} \left[ ^ {t} \mathbf {x} _ {n} \right]

\left[ \delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {L} \right] _ {i j} = \frac {1}{2} \left(\frac {\partial \delta \left(^ {t} \mathbf {N} ^ {t + \Delta t} \mathbf {u} _ {n} + ^ {0} \tilde {\mathbf {N}} \Delta^ {0} \tilde {\mathbf {X}} _ {n}\right)}{\partial \xi^ {i}} \cdot \frac {\partial \Delta_ {0} \mathbf {u}}{\partial \xi^ {j}} + \frac {\partial \Delta_ {0} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial \delta \left(^ {t} \mathbf {N} ^ {t + \Delta t} \mathbf {u} _ {n} + ^ {0} \tilde {\mathbf {N}} \Delta^ {0} \tilde {\mathbf {X}} _ {n}\right)}{\partial \xi^ {j}}\right)

= \frac {1}{2} \left(\frac {\partial \delta^ {t} \mathbf {N} ^ {t + \Delta t} \mathbf {u} _ {n}}{\partial \xi^ {i}} \cdot \frac {\partial \Delta_ {0} \mathbf {u}}{\partial \xi^ {j}} + \frac {\partial \Delta_ {0} \mathbf {u}}{\partial \xi^ {i}} \cdot \frac {\partial \delta^ {t} \mathbf {N} ^ {t + \Delta t} \mathbf {u} _ {n}}{\partial \xi^ {j}}\right)

= \left[ \delta^ {t + \Delta t} \mathbf {u} _ {n} \right] ^ {T} \underbrace {\frac {1}{2} \left(\left[ \frac {\partial^ {t} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \cdot \left[ \frac {\partial^ {t} \mathbf {N}}{\partial \xi^ {j}} \right] + \left[ \frac {\partial^ {t} \mathbf {N}}{\partial \xi^ {i}} \right] ^ {T} \cdot \left[ \frac {\partial^ {t} \mathbf {N}}{\partial \xi^ {j}} \right]\right)} _ {[ \mathbf {c} ] _ {i j}} \left[ \Delta^ {t} \mathbf {u} _ {n} \right]

다시 가상일 항으로 돌아와서 위에 구한 Green-Lagrange strain을 대입하면

\int_ {V _ {0}} \underbrace {\left[ \delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {C} \right] _ {i j} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {S} _ {0} \right] ^ {i j}} _ {\text {constant term}} d V _ {0} + \int_ {V _ {0}} \underbrace {\left[ \delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {C} \right] _ {i j} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {S} _ {C} \right] ^ {i j} + \left[ \delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {L} \right] _ {i j} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {S} _ {0} \right] ^ {i j}} _ {\text {linear term}} d V _ {0}

\left( \begin{array}{l} \left[ \delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {C} \right] _ {i j} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {S} _ {0} \right] ^ {i j} = \left[ \delta^ {t + \Delta t} \mathbf {u} _ {n} \right] ^ {T} \left[ \mathbf {a} \right] _ {i j} \left[ \begin{array}{c} t \\ \mathbf {x} _ {n} \end{array} \right] \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {C} \right] ^ {i j k l} \frac {1}{2} \left[ \begin{array}{c} t \\ \mathbf {x} _ {n} + ^ {0} \mathbf {X} _ {n} \end{array} \right] ^ {T} \left[ \mathbf {e} \right] _ {k l} \left[ \begin{array}{c} t \\ \mathbf {x} _ {n} - ^ {0} \mathbf {X} _ {n} \end{array} \right] \\ \left[ \delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {C} \right] _ {i j} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {S} _ {C} \right] ^ {i j} = \left[ \delta^ {t + \Delta t} \mathbf {u} _ {n} \right] ^ {T} \left[ \mathbf {a} \right] _ {i j} \left[ \begin{array}{c} t \\ \mathbf {x} _ {n} \end{array} \right] \left[ \begin{array}{c} t + {\Delta t} \\ 0 \end{array} \mathbf {C} \right] ^ {i j k l} \left[ \begin{array}{c} t \\ \mathbf {x} _ {n} \end{array} \right] ^ {T} \left[ \mathbf {a} \right] _ {k l} \left[ \begin{array}{c} \Delta^ {t} \mathbf {u} _ {n} \end{array} \right] \\ \left[ \delta_ {0} ^ {t + \Delta t} \mathbf {E} _ {L} \right] _ {i j} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {S} _ {0} \right] ^ {i j} = \left[ \delta^ {t + \Delta t} \mathbf {u} _ {n} \right] ^ {T} \left[ \mathbf {c} \right] _ {i j} \left[ \begin{array}{c} \Delta^ {t} \mathbf {u} _ {n} \end{array} \right] \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {C} \right] ^ {i j k l} \frac {1}{2} \left[ \begin{array}{c} t \\ \mathbf {x} _ {n} + ^ {0} \mathbf {X} _ {n} \end{array} \right] ^ {T} \left[ \mathbf {e} \right] _ {k l} \left[ \begin{array}{c} t \\ \mathbf {x} _ {n} - ^ {0} \mathbf {X} _ {n} \end{array} \right] \end{array} \right)

= \left[ \delta^ {t + \Delta t} \mathbf {u} _ {n} \right] ^ {T} \left\{ \begin{array}{l} \left(\int_ {V _ {0}} \left[ \mathbf {a} \right] _ {i j} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {S} _ {0} \right] ^ {i j} d V _ {0}\right) \left[ \begin{array}{c} t \\ \mathbf {x} _ {n} \end{array} \right] \\ + \left(\int_ {V _ {0}} \left[ \mathbf {a} \right] _ {i j} \left[ \begin{array}{c} t \\ \mathbf {x} _ {n} \end{array} \right] \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {C} \right] ^ {i j k l} \left[ \begin{array}{c} t \\ \mathbf {x} _ {n} \end{array} \right] ^ {T} \left[ \mathbf {a} \right] _ {k l} + \left[ \mathbf {c} \right] _ {i j} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {S} _ {0} \right] ^ {i j} d V _ {0}\right) \left[ \Delta^ {t} \mathbf {u} _ {n} \right] \end{array} \right\}

와 같다. 이제 우변의 항들을 정리해보면 아래와 같다. 여기서 body force에 대한 영향은 무시한다.

\begin{array}{l} \int_ {\partial V _ {0 m}} \delta^ {t + \Delta t} \mathbf {u} _ {0} ^ {t + \Delta t} \mathbf {F} _ {0} ^ {t + \Delta t} \mathbf {S} ^ {t + \Delta t} \tilde {\mathbf {n}} d A _ {0} = \int_ {\partial V _ {0 m}} \delta \left(^ {t} \mathbf {N} ^ {t + \Delta t} \mathbf {u} _ {n} + ^ {0} \tilde {\mathbf {N}} \Delta^ {0} \tilde {\mathbf {X}} _ {n}\right) _ {0} ^ {t + \Delta t} \mathbf {F} _ {0} ^ {t + \Delta t} \mathbf {S} ^ {t + \Delta t} \tilde {\mathbf {n}} d A _ {0} \\ = \int_ {\partial V _ {0 m}} \delta \left(^ {t} \mathbf {N} ^ {t + \Delta t} \mathbf {u} _ {n}\right) _ {0} ^ {t + \Delta t} \mathbf {F} _ {0} ^ {t + \Delta t} \mathbf {S} ^ {t + \Delta t} \tilde {\mathbf {n}} d A _ {0} \\ = \left[ \delta^ {t + \Delta t} \mathbf {u} _ {n} \right] ^ {T} \int_ {\partial V _ {0 m}} \left[ ^ {t} \mathbf {N} \right] _ {0} ^ {T t + \Delta t} \mathbf {F} _ {0} ^ {t + \Delta t} \mathbf {S} ^ {t + \Delta t} \tilde {\mathbf {n}} d A _ {0} \\ = \left[ \delta^ {t + \Delta t} \mathbf {u} _ {n} \right] ^ {T} \int_ {\partial V _ {0 m}} \left[ ^ {t} \mathbf {N} \right] ^ {T} [ \mathbf {t} ] d A _ {0} \\ \end{array}

따라서 가상변위를 지워 모든 식을 정리하면

\begin{array}{l} \int_ {V _ {0}} \delta^ {t + \Delta t} \mathbf {u} \cdot \rho_ {0} ^ {t + \Delta t} \ddot {\mathbf {u}} d V _ {0} + \int_ {V _ {0}} \delta^ {t + \Delta t} _ {0} \mathbf {E}: ^ {t + \Delta t} _ {0} \mathbf {S} d V _ {0} = \int_ {\partial V _ {0 m}} \delta^ {t + \Delta t} \mathbf {u} ^ {t + \Delta t} _ {0} \mathbf {F} ^ {t + \Delta t} _ {0} \mathbf {S} ^ {t + \Delta t} \tilde {\mathbf {n}} d A _ {0} + \int_ {V _ {0}} \delta^ {t + \Delta t} \mathbf {u} \rho_ {0} ^ {t + \Delta t} \mathbf {f} d V _ {0} \\ \Rightarrow \underbrace {\int_ {V _ {0}} \rho \left[ ^ {t} \mathbf {N} \right] ^ {T} \left[ ^ {t} \mathbf {N} \right] J d V _ {0}} _ {\mathbf {M}} ^ {t + \Delta t} \ddot {\mathbf {u}} + \underbrace {\int_ {V _ {0}} \left[ \mathbf {a} \right] _ {i j} \left[ ^ {t + \Delta t} {} _ {0} \mathbf {S} _ {0} \right] ^ {i j} d V _ {0} {} ^ {t} \mathbf {x}} _ {\mathbf {f} _ {\text {int}}} \\ + \underbrace {\int_ {V _ {0}} \left[ \mathbf {a} \right] _ {i j} \left[ ^ {t} \mathbf {x} _ {n} \right] \left[ ^ {t + \Delta t} _ {0} \mathbf {C} \right] ^ {i j k l} \left[ ^ {t} \mathbf {x} _ {n} \right] ^ {T} \left[ \mathbf {a} \right] _ {k l} + \left[ \mathbf {c} \right] _ {i j} \left[ ^ {t + \Delta t} _ {0} \mathbf {S} _ {0} \right] ^ {i j} d V _ {0}} _ {\mathbf {K} _ {t}} \Delta^ {t} \mathbf {u} = \underbrace {\int_ {\partial V _ {0 m}} \left[ ^ {t} \mathbf {N} \right] ^ {T} [ \mathbf {t} ] d A _ {0}} _ {\mathbf {P} _ {d i s t}} + \mathbf {P} _ {c o n} \\ \Rightarrow \mathbf {M} ^ {t + \Delta t} \ddot {\mathbf {u}} + \mathbf {K} _ {t} \Delta^ {t} \mathbf {u} = \mathbf {P} _ {\text { dist }} + \mathbf {P} _ {\text { con }} - \mathbf {f} _ {\text { int }} \\ \end{array}

와 같이 정리할 수 있다. 여기서 M은 mass matrix, K_{t} 는 tangent stiffness matrix, P_{dist} 는 분포하중에 의한 힘, P_{con} 는 집중하중, f_{int} 는 변형에 의한 힘을 나타낸다.

4. Constitutive matrix

Plane stress 가정을 사용하는 구성행렬은 다음과 같다.

\left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {C} \right] _ {x ^ {1} x ^ {2} x ^ {3}} = \frac {E}{1 - \nu^ {2}} \left[ \begin{array}{c c c c c c} 1 & \nu & 0 & 0 & 0 & 0 \\ n & 1 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & \kappa \frac {1 - \nu}{2} & 0 & 0 \\ 0 & 0 & 0 & 0 & \kappa \frac {1 - \nu}{2} & 0 \\ 0 & 0 & 0 & 0 & 0 & \kappa \frac {1 - \nu}{2} \end{array} \right]

위의 구성 행렬은 local Cartesian coordinate에서 정의되었기 때문에 transformation matrix를 이용하여 natural coordinate로 바꾸어 줄 수 있다.

\left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {C} \right] _ {\xi^ {1} \xi^ {2} \xi^ {3}} = \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {T} \right] ^ {T} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {C} \right] _ {x ^ {1} x ^ {2} x ^ {3}} \left[ \begin{array}{c} t + \Delta t \\ 0 \end{array} \mathbf {T} \right]

이 때 전단보정계수 \kappa 는 \frac{5}{6} 를 사용하였다.