第一工業大学研究報告第28号（2016）pp.31－34

31

第一工業大学研究報告
第 28 号 (2016)

しらす斜面崩壊警報システムにおけるしらす斜面崩壊警報システムにおけるひずみ測定の精度に関する研究

ひずみ測定の精度に関する研究
山田

山田
1

2

第一工業大学

第一工業大学3

猛矢1・福重彩夏2・岡林巧3

猛矢 1 ・ 福重

彩夏 2 ・ 岡林

巧3

情報電子システム工学科（〒 899-4395 鹿児島県霧島市国分中央 1-10-2）
E-mail: t-yamada@daiichi-koudai.ac.jp
情報電子システム工学科卒業（〒 899-4395 鹿児島県霧島市国分中央 1-10-2）

第一工業大学

自然環境工学科（〒 899-4395 鹿児島県霧島市国分中央 1-10-2）
E-mail: t-okabayashi@daiichi-koudai.ac.jp

A study of the accuracy of the strain measurement
on Shirasu slope failure warning system
Takeshi YAMADA，Ayaka FUKUSHIGE，Takumi OKABAYASHI
Department of Information and Electronics Systems Engineering, Daiichi Institute of Technology（1-10-2 Kokubu Chuo, Kirishima City, Kagoshima 899-4395, Japan）
E-mail: t-yamada@daiichi-koudai.ac.jp
Department of Environmental Engineering, Daiichi Institute of Technology（1-10-2 Kokubu Chuo, Kirishima City, Kagoshima 899-4395, Japan）

E-mail: t-okabayashi@daiichi-koudai.ac.jp

Abstract : This paper is a study on the accuracy of the strain measurement of Shirasu slope failure warning
system produced last year. It was discussed accuracy of the device by comparing the experimental results withthe strain measuring device was fabricated and the theoretical value. As a result, it was found that the causeof the error is due to the fluctuation of the acceleration sensor itself and the numerical calculation algorithm.
In addition, the error was found to be suppressed to 0.04 m.
Keywords : Shirasu slope failure, accuracy of the strain measurement
1. はじめに
本研究は，昨年度作製した「しらす斜面崩壊警報システム」[1] のひずみ測定装置の精度についての研究である。現段階において本システムの斜面崩壊予測は，ひ定装置で実験を行い，理論値との比較を行ったうえで，どこまで精度を高められるかを考察する。

2. ひずみ測定装置の概要

ずみ測定のみで行っている。そのため，ひずみ測定のまず昨年度作製したひずみ測定装置の概要について精度を正確に知ることは，斜面崩壊の予測を行う上で記述する。本装置は，Raspberry Pi と呼ばれる超小型欠かすことができない。本研究は，作製したひずみ測コンピュータ（図 1）と ST マイクロ社の加速度セン

第一工業大学研究報告

32

サー lis3dh （図 2）を利用し作成された。Raspberry

第28号（2016）

なり大きな誤差が含まれる結果となった。総時間とは移動させた時間である。1 秒間移動させたが，装置で計測された時間は平均 1.01 秒となった。また計算処理
時間の関係で加速度の取得時間間隔は平均 0.145 秒となった。

変位 < r >図 1: Raspberry Pi

図 2: 加速度センサー

1.38 m

表 1: 実験結果
総時間 加速度の取得間隔 ∆t

1.01 s0.145 sPi に加速度センサーを搭載することで，自身の変位が計測可能となる。具体的には，i 回目に x 方向の加速度ax (i) が得られた時，速さの変化量∆vx (i) = ax (i)∆t

次に混入する誤差について，理論的に考察する。混(1)から（ただし，∆t は加速度の取得時間間隔），そのときの変位∆x(i) = vx (i − 1)∆t +vx (i) =i


1ax (i)∆t22∆vx (j)4. 混入誤差
入する誤差の要因としては，加速度センサー自身の持つ揺らぎによるものと計算アルゴリズムによるものの2 つが考えられる。そこで，これら 2 つの要因によりどの程度誤差が混入するかを考える。

(2)

4.1 加速度センサー自身の持つ揺らぎ(3)

加速度センサーの取得する値は，真の値のまわりでを計算する。これを足し合わせることで変位が求まる。

せた状態での x 方向の取得加速度（10000 回分）の分揺らいでいる。図 3 は加速度センサーを水平に静止さj=1式 (1)，(2)，(3) を用いて加速度を n 回取得した時の変位を計算すると

x(n) =n布である。この分布は平均 0，分散 σx2 = 0.00184 のガウス分布となった。今回我々が用いた数値計算は式 (6)で変位を計算するため，この程度の揺らぎのある加速∆x(i)(4)∆t2   2aj + ai
=2 i=1j=1ni−1(5)ぎは，足し合わせる数が多くなるほど大きくなる。
毎回取得する加速度の母分散 σx2 は変わらないので，加速度を n 回取得した時の分散 σx2n は，式 (6) よりとなる。これをさらに整理しながら計算するとn∆t2  2i − 1 aix(n) =2 i=1度の値を足し合わせていくことになる。そのため計算により得られた変位も揺らぎを持っており，その揺ら

i=1n ∆t2 σx2n = n(6)となる。y ，z 方向も同様に計算することで，自身の変位（ r = x2 + y 2 + z 2 ）が得られる。
これを斜面に設置することで，ひずみが生じた際，本装置の変位によりひずみの大きさ，また成長速度を得 ∆t2=n差でひずみを測定できるか調べた。まず本装置を水平面上で，1 秒間に 1m 移動させるという実験を 1000 回行った。表 1 が実験結果であり，変位 < r > は変位
r = x2 + y 2 の 1000 回の平均だが，結果は 1.38m と· n22σx2

いるとσx2n =大きくなるときに斜面が崩壊すると予測し警報を出す。

我々はこのひずみ測定装置を使って，どの程度の誤2i=1σx2となる。この式に 1 秒あたりの取得回数 n =ることができる。我々は，ひずみの成長速度が急激に

3. 実験22(2i − 1)

1 2
σ4∆t x(7)(8)1を用∆t(9)となる。図 6 は，実験結果（ σx2 = 0.00184，∆t = 0.145 ）を代入したものである。グラフを見てわかるように，∆tが 0 に近づく（取得回数が多い）と σx2n が大きくなることがわかる。

実験では，装置を水平面上で移動させているので，x方向の議論だけでは不十分である。そこで y 方向も考x2n + yn2 ）を求める。σr2 を求2r=慮した σr（ただし，

山田・福重・岡林：しらす斜面崩壊警報システムにおけるひずみ測定の精度に関する研究

0.010.008

P (a)

0.006

σx2n0.004

0.002

000.20.4

0.60.81∆t
図 4: 分散 σx2nax
図 3: 加速度 ax の分布

 x 2nσxn+ y 2nσy nである。つまり 1 回当たりの x 方向の誤差 (n+1)∆t1v(t)dt∆Ex = v(n∆t)∆t + a(n∆t)(∆t)2 −2n∆t(∆t)2  2a (n + 1)∆t − a(n∆t)(17)
2となり，1 回当たりこれだけの誤差が含まれることになる。式 (17) より (∆t)2 に比例した誤差が混入してくることがわかる。図 6 は，実験から得られる 2a (n +1)∆t − a(n∆t) = 3.56 を代入したグラフである。

めるにあたり，まずr2 =傾き a(n∆t)
Dv(t)

(10)

ynxn，はそれぞれ平均 0，分散 1 のガσxn σy nウス分布に従うので，r2 は自由度 2 の χ2 分布に従う。

を考える。

AまたB
σx2n  σy2n(11)(12)σx2n r2∆Ex 0.40.2C0n∆t(n + 1)∆t
図 5: 混入誤差

4−π 2σxnσr2 =2(13)となり，式 (9) を代入すると

σr2 =4−π 2σ8∆t x(14)となる。この式に実験で得られた ∆t = 0.145，σx2 =0.00184 を代入すると σr2 = 0.00136 となる。
4.2 アルゴリズムにより混入する誤差

次に数値計算アルゴリズムにより混入する誤差について考える。今回採用した計算は，取得した加速度 a1から単純に等加速度直線運動の式 x = v0 t + at2 に代2入し求めたものである。つまり n 回目に取得した加速度 a(n∆t) を用いて

∆x = v(n∆t)∆t +1a(n∆t)(∆t)22(15)とし，微小変位 ∆x を計算した。しかし，この計算では図 5 の台形 ABCD の面積を求めており，真の変位ではない。真の変位は∆x = (n+1)∆tn∆tv(t)dt(16)0.20.40.60.81∆t
図 6: 誤差 ∆Ex算した時の混入誤差2ているので，変数変換を行うことで r の分散 σr2 を求めることができる。実際に σr2 を計算すると

01 回当たりにこれだけの誤差が混入するので，n 回計となる。つまり r が自由度 2 の χ 分布に従うとわかっ210.80.6と近似すると

r=33

Ex =n

∆Ei(18)i=1(19)
 n∆Ex∆t  2a (n + 1)∆t − a(n∆t)(20)=2となる。図 7 は，実験から得られた 2a (n + 1)∆t −a(n∆t) = 3.56 を代入したものである。式 (20) に
1および実験から得られた ∆t = 0.145，2a (n +n=∆t
1)∆t − a(n∆t) = 3.56 を代入すると Ex = 0.261 となる。

y 方向も同様に誤差が混入する仮定すると，誤差 Eは図 8 のようになり，計算すると
(21)E = Ex2 + Ey2
 0.368(22)となる。

5. 実験値と理論値の比較
5.1 実験値と理論値の比較
表 2 が実験値と理論値をまとめたものである。実験
では装置を 1m 移動させたが，測定結果（< r >）は1.38m と大きなずれがあった。しかし，数値計算により混入する誤差が E = 0.368m，また加速度センサー
自身の持つ揺らぎの分散 σr2 = 0.00136（σr = 0.0369）

第一工業大学研究報告34により混入する誤差であることが，実験値と理論値の1Ey0.8

Ex 0.4比較によりわかった。さらにこの装置においては，計E0.6

算上ではあるが誤差を 0.04m まで抑えられることもわかった。

y

ひずみ測定は「斜面崩壊警報システム」において，斜0.20第28号（2016）

00.2

0.4

0.60.8面の崩壊を予測する非常に重要な要素である。今のま1∆tx

図 7: 誤差 ExExまでは誤差の混入が大きく実用化は難しい。アルゴリ

図 8: 誤差 E

ズムの工夫や，より高性能な加速度センサーを用いるということを考えると，この結果は妥当なものであり，誤差の要因はこの 2 つであったといえる。

<r>

表 2: 実験値と理論値<t>∆tσr21.381.010.1450.00136などしてひずみ測定の精度を上げることが今後の課題である。
謝辞

本研究は，第一工業大学研究開発助成金の支援によσrE

0.03690.368

り実施された。ここに記して深く感謝の意を表する。
参考文献[1] 山田猛矢，岡林巧: しらす斜面崩壊警報防災システ

5.2 誤差を抑える
そこで，計算アルゴリズムを工夫することにより，どれだけ誤差の混入を抑えることができるか考える。これまで見てきたように，∆t を小さくすれば誤差 E を抑えることができるが，∆t が小さくなることにより取得回数が増えるため分散 σr2 は大きくなってしまう。では，誤差をできるだけ抑えるためには ∆t をどのように取ればよいだろうか。
図 9 が σr2 ，E および 2 つを足し合わせた σr2 + E のグラフである。この σr2 +E を最小にする ∆t が，我々の使用した装置，行った数値計算において，分散を抑えつつ誤差を最小にする ∆t である。その値は ∆t = 0.00846である。数値計算のアルゴリズムを工夫し，仮にこの値が実現できれば，混入する誤差は 0.04m 程度に抑えられる。

σr2 + EE
σr2∆t
図 9: 誤差を抑える

6. まとめと今後の課題

昨年度作製したひずみ測定装置を使って変位測定の実験を行ったところ，平均で 0.368m の誤差が混入することがわかった。またその要因が，加速度センサー自身の持つ揺らぎによるものと数値計算アルゴリズム

ムにおける情報収集無線ネットワークシステムの構
築，第一工業大学研究報告 第 27 号

[2] 福重彩夏: 斜面崩壊警報システムにおけるひずみ測定に関する研究，第一工業大学 卒業論文 2015