本稿では，砂の力学モデルとしてのひずみ空間での多重せん断モデル(Iai＆Ozutsumi, 2005)に基づく線形弾性体の構成式を提案し，幾何学的非線形性を考慮した大変形解析(有限ひずみ解析)に必要な定式化を示す。本定式化では，体積ひずみと等方圧力，および，仮想単純モデルにおけるせん断ひずみとせん断応力との間に線形関係を仮定している。定式化を基に，有限要素法による数値解析に必要な離散化を行い，Total Lagrangian法，および，Updated Lagrangian法の両者によるプログラミングを実施することで得られた要素解析例を併せて示す。