
メタマテリアルを均質化してみた
Homogenizing Mechanical Metamaterial
あるメタマテリアル構造を効率的にシミュレーションするため,均一なシェル要素に置き換えます。
In order to simulate a metamaterial structure efficiently, it is replaced by a uniform shell element.

山村悟史
Satoshi Yamamura
2023,09,08 2023,09,08ミッドソールに3Dプリントされた特殊な構造が埋め込まれたアディダスのシューズ「4DFWD」は踏み込み力を推進力に変換する機能を持ちます。本記事ではこの機能を圧縮をせん断に変換するメカニカルメタマテリアルとして解釈し、その設計方法を解説します。本記事では圧縮をせん断に変換するマクロ物性について材料力学の観点からその特性を紐解きます。
The adidas "4DFWD" shoe, which has a special 3D-printed structure embedded in the midsole, has the function of converting stepping force into propulsive force. This article interprets this function as a mechanical metamaterial that converts compression into shear and explains how to design it. This article explains the macro physical properties that convert compression into shear from the viewpoint of material mechanics.
須藤海
Kai Suto
2022,11,04 2022,11,04ミッドソールに3Dプリントされた特殊な構造が埋め込まれたアディダスのシューズ「4DFWD」は踏み込み力を推進力に変換する機能を持ちます。本記事ではこの機能を圧縮をせん断に変換するメカニカルメタマテリアルとして解釈し、その設計方法を解説します。
本記事は以下の3部構成のうち1つ目について解説します。
(2) 圧縮をせん断に変換するメカニカルメタマテリアルの設計
(3) 踏み込み力を推進力に変換するミッドソールの設計
アディダスのシューズ「4DFWD」は踏み込み力を推進力に変換し、ランニング時のブレーキを抑制する効果を持ちます。この力の変換は3DPされたミッドソールのラティスと呼ばれる複雑形状によって実現されています。ラティスはユニットセルと呼ばれる単位構造がXYZ方向に繰り返される周期構造を持っており、このユニットセルが持つ弾性率などの力学特性をある種の材料特性とみなすことができます。ラティスの形状を工夫することで、通常の材料では実現できないような特性、例えばこのミッドソールにおいては踏み込み力を推進力に変換する特性、すなわち「圧縮をせん断に変換する」材料特性を設計することができます。このように幾何的形状によって特殊な材料物性が実現された材料を「メカニカルメタマテリアル (mechanical metamaterials)」と呼びます。
[引用:ADIDAS 4D FWD. Innovation in shock absorption! Spoiler: Hyperbola! (Youtube)]
それでは4DFWDのミッドソールが持つ「圧縮をせん断に変換」する材料特性とは何か?について解説します。 今回ミッドソールを真横からみることで、2次元の問題として考えます。 2次元の材料が変形する自由度は以下となります。
\varepsilon_x
:X方向の伸びひずみ\varepsilon_y
:Y方向の伸びひずみ\gamma_{xy}=\gamma_{yx}
:せん断ひずみこれらをまとめてひずみベクトル\boldsymbol\varepsilon = [\varepsilon_x, \varepsilon_y, \gamma_{xy}]^\mathrm{T}
と表記します。
上記材料の変形に対応する応力は以下の3つです。
\sigma_x
:X方向の引張応力\sigma_y
:Y方向の引張応力\tau_{xy}=\tau_{yx}
:せん断応力これらをまとめて応力ベクトル\boldsymbol\sigma = [\sigma_x, \sigma_y, \tau_{xy}]^\mathrm{T}
と表記します。
次に材料に応力が発生した時にどのようにひずむのか?という以下のような関係式を考えます。
\boldsymbol\varepsilon = \mathbf{S} \boldsymbol\sigma
\tag{1-1}
これは成分表示すると以下のようになります。
\left[
\begin{array}{c}
\varepsilon_x \\\\
\varepsilon_y \\\\
\gamma_{xy}
\end{array}
\right]
=\left[
\begin{array}{ccc}
S_{11} & S_{12} & S_{13} \\\\
S_{21} & S_{22} & S_{23} \\\\
S_{31} & S_{32} & S_{33}
\end{array}
\right]
\left[
\begin{array}{c}
\sigma_x \\\\
\sigma_y \\\\
\tau_{xy}
\end{array}
\right]
\tag{1-2}
この関係式を2次元線形弾性材料の構成則と呼び、行列\mathbf{S}
をコンプライアンス行列と呼びます。
材料が均質に分布した等方性材料ではコンプライアンス行列はヤング率(縦弾性係数)E
、せん断率(横弾性係数)G
、ポアソン比\nu
を用いて以下のように表されます。
\mathbf{S} =
\left[
\begin{array}{ccc}
\frac{1}{E} & -\frac{\nu}{E} & 0 \\\\
-\frac{\nu}{E} & \frac{1}{E} & 0 \\\\
0 & 0 & \frac{1}{G}
\end{array}
\right]
\tag{1-3}
等方性材料の場合、せん断率G
はヤング率E
とポアソン比\nu
を用いて以下のように表されます。
G = \frac{E}{2(1+\nu)}
\tag{1-4}
それではこの関係式上で圧縮がせん断に変換されるとはどういうことかを説明します。
上面に圧縮力を与えるということはY方向に応力-\sigma_y
を与えることに対応します。
引張方向が正となっているので圧縮方向は負になることにご注意ください。
このとき応力ベクトルは\boldsymbol\sigma=[0, -\sigma_y, 0]^\mathrm{T}
となり生じるせん断ひずみ\gamma_{xy}
は式\text{(1-2)}
より以下の式で表されます。
\gamma_{xy} = S_{32} (-\sigma_y)\tag{2-1}
したがって圧縮応力をせん断ひずみに変換するためにはS_{32}
が0でないような材料特性を設計する必要があります。
先程紹介した等方性材料の場合、S_{32}=0
となるためせん断ひずみは発生しません。
それでは材料を回転させてみるとどうでしょうか?
材料を角度\theta
回転させることはコンプライアンス行列に右から応力の変換行列\mathbf{T}_\sigma(\theta)
、左からひずみの変換行列\mathbf{T}_\varepsilon(\theta)
の逆行列をかけることに対応し、新たなコンプライアンス行列を\mathbf{S}(\theta)
とすると以下のようになります。
\mathbf{S}(\theta)=\mathbf{T}_{\varepsilon}(\theta)^\mathrm{-1}\mathbf{S}\mathbf{T}_{\sigma}(\theta)\tag{2-2}
ここで変換行列\mathbf{T}_{\sigma}(\theta)
, \mathbf{T}_{\varepsilon}(\theta)
は以下であらわされます。
\begin{align}
\mathbf{T}_\sigma(\theta)&=
\left[
\begin{array}{ccc}
\cos^2{\theta} & \sin^2{\theta} & 2\cos{\theta}\sin{\theta} \\\\
\sin^2{\theta} & \cos^2{\theta} & -2\cos{\theta}\sin{\theta} \\\\
-\cos{\theta}\sin{\theta} & \cos{\theta}\sin{\theta} & \cos^2{\theta} - \sin^2{\theta}
\end{array}
\right]\tag{2-3}\\
\mathbf{T}_\varepsilon(\theta)&=
\left[
\begin{array}{ccc}
\cos^2{\theta} & \sin^2{\theta} & \cos{\theta}\sin{\theta} \\\\
\sin^2{\theta} & \cos^2{\theta} & -\cos{\theta}\sin{\theta} \\\\
-2\cos{\theta}\sin{\theta} & 2\cos{\theta}\sin{\theta} & \cos^2{\theta} - \sin^2{\theta}
\end{array}
\right]\tag{2-4}
\end{align}
これを計算すると圧縮応力をせん断ひずみに変換するS(\theta)_{32}
は以下となります。
S(\theta)_{32} = \sin{4\theta}\cdot\frac{1}{4}\left(\frac{1}{G} - \frac{2(1+\nu)}{E}\right)\tag{2-5}
等方性材料の場合式\text{(1-4)}
が成り立つので回転させてもS(\theta)_{32}=0
となります。
よって等方材料の場合圧縮応力をせん断ひずみに変換することはできないことがわかりました。
逆に式(\text{1-4)}
が成り立たない特殊な材料物性を設計することができれば、式\text{(2-5)}
からS(\theta)_{32} \neq 0
が成り立ち、圧縮応力をせん断ひずみに変換することが可能となります。このような線形弾性物性を持つ材料のクラスを正方異方性材料(square anisotropic material)と呼びます。
以下に等方性材料と正方異方性材料の回転角度毎の方向ヤング率(directional Young's modulus)を示しました。 等方性材料の場合、どの方向でも一定のヤング率を持ちますが、正方異方性材料は角度毎に異なるヤング率を示します。角度によって異なるヤング率を持つ材料を使うことで「圧縮応力をせん断ひずみに変換する」などの不思議な力学的現象を作り出すことができるのです。
「圧縮がせん断に変換される」とはどういうことか?を材料力学を用いて説明しました。 等方性材料は上記特性を持たないことを示し、特殊な材料物性を持つ正方異方性材料が必要であることがわかりました。
次回、正方異方性物性を持つユニットセル構造を紹介し、具体的にどのような変形をするか解説します。
Nature Architects株式会社
Nature Architects, Inc.
代表取締役:大嶋 泰介
CEO:Taisuke Ohsima
所在地:東京都港区赤坂8-4-7 アパートメントカーム3D
ADDRESS:Apartment Calm 3D, 8-4-7, Akasaka, Tokyo, 107-0052, Japan
設立:2017年5月
Funded:May, 2017
URL:https://nature-architects.com/
あるメタマテリアル構造を効率的にシミュレーションするため,均一なシェル要素に置き換えます。
In order to simulate a metamaterial structure efficiently, it is replaced by a uniform shell element.
山村悟史
Satoshi Yamamura
2023,09,08 2023,09,08Karamba3dは3DCADのRhinocerosのGrasshopper上で動く構造解析プラグインです。 Karamba3dはAnsysのような有名な商用CAEではないのでその性能について気になりますが、 こういった情報はあまりないため、この記事では簡単なモデルを対象に確認していきます。
Karamba3d is a structural analysis plug-in that runs on top of 3D CAD Rhinoceros Grasshopper. Karamba3d is not a well-known commercial CAE like Ansys, so I am curious about its performance, However, as there is not much information on this, this article will check it on a simple model.
夏目大彰
Hiroaki Natsume
2023,06,09 2023,06,09