日本包装学会誌VblL3ﾉVbL4ag94）

一ﾉ職論文

衝撃モデルによる倒立剛体の転倒限界
西原主計＊横山孝之率＊OverturnThreshoIdsforTheInvertedRigidBody
toShockModels

ＫａｚｕｅＮＩＳＨＩＨＡＲＡ.，ＴａｋａｙｕｋｉＹＯＫＯＹＡＭＡ｡。

Tomaintainthesafetyofpr℃ducts，ｉｔｉｓｉｍｐｏＵｔａｎｔｔｏｕｎｄｅ応tandthecimumstanceln
whichtheligidbodytumove■andwiththisknowledge，todesignthecontainerhardto
upsetevenintheshoCkenvlronment,Ｉｎｔｈｉｓｐａｐｅｎａｔｆｉｌ石tusmgsometypicalshockmodels，computersimulationswerCachievedtogiveovelturnth1℃sholdaccelelationsonthemverted

ngidbody，ｇｅｔｔｍｇｔｗｏｓｉｍｐｌｅｅｑｕａｔｉｏｎｓｆＯｒｔｈｅｔｈ妃sholdaccelelRationsThetwoequations
werechamctelizedinthecaseoｆｔｈｅｔｉｍｅｍｔｅｇｌｎｌａ1℃ａｏｆｓｈｏｃｋｗａｖｅ，ｔｈａｔｉｓｂｃａｓｅｏｆｔｈｅａｒｅａ≧１ａｎｄｃａｓｅｏｆｔｈｅａｒｅａ＜１，bothofwhichgavegoodappmximationｓｔｏｔｈｅｏriginal
diffe定ntialequationinawiderlangeofnon-dimensionalfi℃quencyQ・
Secondly，theshockexpelimentiscamiedoutfOrasteelbaUovertumingapplatuswith
vanouscushioningpads・Ｉｔｗａｓｅｖｉｄｅｎｔｔｈａｔｔｈｅｔｅｎｄｅｎｃｙｏｆｔheexpe㎡mentalth1℃sholdaccelerationscomcidedwiththatofaboveequati０，s・SooverturnmgthTesholdswillbe

presumedwidelyifthedurationtimeaｎｄｉｔｓｉｎｔｅｇｍｌａ1℃ａｏｆａｃｔｕｍｌｓｈｏｃｋａ尼partlyassumed．

Keywords：Overturnofrigidbody,Ovel･tumofcontainer,OverturnthlCshold,TurnoverbFall，
UpseLShockmodel

製品の安全性を確保するためには、物体の転倒はどういう場合に起るかを知ることと、その知見により転倒しにくい容器を設計していくことが重要である。

そこで本研究では、まず、各種衝撃波形のモデルを用い、剛体の転倒限界を与える計算機シミュレーシ
ョンを行い､二つの簡便な転倒限界式を導いた。その簡易転倒限界式は､衝撃面積≧１又はく１の場合に分けられ、無次元化振動数Ｑの広範囲にわたって、元の転倒微分方程式に対しよい近似を与えた。
次に、鋼球の転勤装置と各種緩衝材による衝撃実験を行った。その実験による衝撃転倒限界加速度の特性は､上に求めた簡易転倒限界式の特性傾向に合致した。よって､衝撃持続時間･衝撃面liKiが幾つ力､推測できれば、上の簡易転倒限界式から転倒限界が推定できることが分かった。

キーワード：容器の転倒、倒立剛体、転勤、倒れ、衝撃転倒。神奈川工科大学（〒243-02神奈川県厚木市下荻野1030)：KanagawalnstituteofTechnology,1030,Shimo-Ogino，Atsugi-shi,Kanagawa,243-02.．二国機械工業（株）（〒213神奈川県川崎市高津区久地843-5）：Nikuni-Kikai
KogyoCo.,843-5,Kuji,Takatsu-ku,Kawasaki-shi,Kanagawa,２１３－２４９－

衝浄ピモ完ルによる釘、;W汝の転倒限界

1．まえがき2.運動方程式

包装物輸送において、これまで輸送振動、水平方向の加速度を受けるＦｉｇ．１のような落下衝撃､回転衝撃など､主として大振動･大転倒系を考える。同図の（a）と（b）は同一衝撃に関する動力学的研究がなされ、緩衝包方程式系であるが､以下Fig.１（a)すなわち、装技術は進展してきた。しかし、輸送過程・倒立剛体の質量を、､重心高さをh､上面幅を消費者環境において、個々の製品の倒れ、躍2a､底面幅を2ｄとし、水平方向の加速度ｕにり・衝突などの動力学的研究は地味で数が少より転倒する系について調べる。重心Ｏからないように思われる。

底面端Ｐまでの距離をr､重心を通る鉛直軸と製品の中には瓶類、缶類のごとく脊の高いOPとの傾き角を０、その静止状態での傾き角ものがあり、また、積み重ねたとき倒れ易すを０．、重心回りの回転角をａ（すべり回転成ぐなるものがある。加わる振動・衝撃が小さ分を含む)、重心の相対速度をｖ、Ｐ点での剛くても製品は躍り、他と衝突し、転倒する場

体のすべり速度をv`、抗力をＲ、摩擦力をＦ、合が見受けられる。転倒自体が直ちに製品の動摩擦係数を似とし、図の矢印の方向を正に品質を損うものではないが、容器を傷付けたとる。０点、Ｐ点での慣性モーメントをそれり、場合によっては破損し安全性にかかわるぞれL､ＩＤとする。Ｐ点において剛体の傾きがことが予想される。したがって、流通過程に静止状態から増加する方向に①をとる。重心おいても消費者環境においても、製品の安全Ｏが床点Ｐと垂直になるときを転倒限界とす性を最低限確保するためには、少なくとも物る。

体の転倒はどういう場合に起るかを知ることと、その知見により転倒しにくい容器を設計Ｒしていくことが重要である。

本研究では、大局的な転倒特性を得るた
Ｖｙ

め、まず製品及び床を剛体と考え、先に著者らが導いた転倒近似方程式Ｕ２)のうちすべり・躍りのない状態式に対し、計算機シミュレー

Ｆションにより、各種衝撃波形のモデルを印加して物体の転倒限界を調べ、簡便な転倒限界。ＩｒＴＩｂ

式を導びいた。併せて各種緩衝材による鋼球（a）containerの衝撃転勤実験と比較し検討した。

（b）SteelballFigllnvertedrigidbodies

物体の転倒問題は、包装輸送面では今後の容器設計、陳列・積み重ねなど、家庭では収運動方程式を簡単ｲﾋするため次の仮定を設納、取り扱いなど、また安全面ではプロパンける。

ガスボンベ、それらの転倒検知センサなど3M》（１）Ｐ点で剛体の飛び跳ねはない

多くの分野にまたがる基本問題である。

（２）運動は左右半面で対称である－２５０－

日本包鍵学会誌ＶｂＬ３Ⅳ０．‘α994）（３）０，００，９などの２乗オーダは省略できる

ピックアップを台車に載せ、衝突面に各種形（４）摩擦は乾性摩擦則に従い、かつ,LA＜＜状・厚さの発泡ウレタンフォーム（以下スポ００とするンジと略す)、発泡スチロール、エアキャヅ（５）端点Ｐでの曲率半径は考えない

プ、ゴムシート、丸めたチリ紙及びそれらのここで、すべり速度はvb＝ｖ－ｒａであり、組み合わせによる緩衝材を置き、台車をゴムvs≠Ｏで動摩擦、ｖｂ＝Ｏで静摩擦である。

紐で引っ張り衝突させる。鋼球の転がり上が０＝（汀／2）_β－の、ｔａｎβ＝ｈ／d、ｒ＝り点を見ながら、加速度波形とそのピーク値､/~TIzT~５万である。また、の、０，００、似はをオシロスコープで読み取りＦｉｇ．２を得た。

小さいとして、転倒方程式は換算時間領域Tれら)、回転軸が他端に移ることによる速度変０等による減衰を含んだ変軸速度係数で定義さ．
」●●

ｌかび直前の値を表す。ｅは弾性変形､熱､音伝播

－８●
●

４

ただし、添字の正負は０．に到達した直後及●

８ど'+＝ｅ駐２（□）ロ．［①８ｍ己一・二ｍ①亘］Ｅロゼｇｏ。i鋪P洲｝①１において以下のように表される2)。

」篭|:鏡ａ．
●●、●

0５０１００１５０２００
Durationtime,Ｔ［ｍｓ］

化が主体をなす。時間尺度をＴ＝〃t､′を丁に関する微分、①イーｇ／r、Ｑ＝し／①｡と置Fig2Expe｢imentonshockove『turning

いて、ど＝①／００、初期荷重定数ｂ＝１－’／００，ｐ＝ｕｏ／goo、ｕ･を加速度振幅、ｊFig.２から、衝撃加速度Ｇ(Ｇは重力加速度(「）を最大１及び持続時間兀なる外力関数、を1単位とする）は､実時間領域での衝撃持続

例えば､半正弦波ではしを加振振動数として、時間Ｔが小さいときは緩衝材の材質、形状に‘（て）＝sｉｎ丁である。変位の最大ｌｆｌｍより１Ｇ程度ばらつくが、Ｔが大きくなるにつ＝ｂのときを転倒限界とする。システムは振れて緩衝材に依存せず、ばらつきは小さくな

動数比Ｑと転倒感度p／ｂで特徴付けられる。っている。全体の傾向として、転倒限界が－

なお、式展開は［付録］を参照。

つの単調な曲線の回りにばらついているように見られるし、また、剛体の転倒限界が衝撃

3.実験

持続時間又はその時間積分など、材質による衝撃関数に依存しているようにも見られる。

Ｆｉｇ．１（b）において、鋼球、＝２８１９、そ

材質的にはａはスポンジをちぎったものかの半径ｒ＝0.951ｃｍ、アルミ台の溝半径。＝又はそれとエアキャヅプの重ね合わせ、ｂは0.34ｃｍ、摩擦係数（鉄・アルミ）α＝0.0016）

スポンジ単体又はチリ紙を丸めたもの、ｃはとし、静止状態の傾き角０｡＝0365rad、ｂ＝エアキャップ単体、又は発泡スチロールとス0.997,⑩｡＝32.1である。この装置と加速度ポンジの重ね合わせ、ｄはゴムシート又はそ－２５１－

溺薄gぞ完ルによるβym5V）敏の転倒腹深め、第２章に得られているころがり.転倒方程式!）を整備し、設計などに応用するため、より見通しのよい簡易転倒限界式を導くこととする。まず、換算時間領域の式（１）において、倒れ限界は、
ど＝ｂ、ｒ’＝0て＝て、(2)

ただし、実時間領域のＴは換算時間領域にお(a）shockwave

O5G／div,２０，s／div（exp・No.1）
symbol▼：ｓｔａｒｔａｎｄｓｔｏｐｏｆｉｍｐａｃｔ

いてて、＝汀である。このときの式（１）を以下のような範囲に分ける。すなわち、(1)Ｑ＜＜１のとき：

式（1）は－ど±ｂ＝Ｐ⑩（て）となりど＝ｂのオーダであるから、ｐはＯに近づく。

(2)Ｑ＜１でＱど"＝どの範囲：
このときはＰ＝ｂとなる。数値計算でもこのような範囲がある。

(3)Ｑが１の近傍ないしそれ以上で十分大
きい範囲：

式（１）を直接積分し、(b）apparatus
l：polyurethanefoamcushion

ｗｌｉ鰄十I罪±b川

２：SteelbaUandaccelerometerinthebox

３：triggerequipment
４：guide一prm(r)…。（３）Fig3Shockwaveandapparatus

れと発泡スチロールの重ね合わせ、ｅは発泡

ここに、式（３）の右辺の積分値、すなわちＰの係数を衝撃面積ｓと定義する。左辺第１項

スチロール単体などである。

は零、したがって､左辺第２項は倒れ限界におなお、Fig.３は衝撃値2.6G、持続時間Ｔ＝54,sのときの加速度波形と緩衝材である。

いてＱとｓの関数Ｃ（Ｑ,ｓ）となっていると考えられる。そこで、その加速度波形は正弦波ｎ乗型のように見えps＝Ｃ（Ｑ，ｓ）

るが、実際に衝突している時間はこれより２倍程度長い。

(4)と置く。上の範囲(1)、（2）は特定できないので、（３）によって式（４）に従う簡易転倒限

4.数値シミュレーション界式を求める。
Ｆｉｇ．４に今回シミュレーションに用いた衝ここでは、剛体の転倒限界を数式化するた撃波形のモデル、ＴａｂＩｅｌにそれらの衝撃面－２５２－日本包装学会誌 VbL3jVb． ４口994)1２
３４５６７８９１１

100

5０

Ｔ３Ｔｍ＝汀

(a） shortｊｅｒｋｗｉｔｈｓｍａｎａｎｄｌｏｎｇｏｆｆｓｅｔ
ｆｏｒＮｏ､１ｔｏ７ａｎｄ３１ｔｏ３３０

１０１て２０５２１５０て１ｒｏ2０
１ｑ

Ｃ、□旨冒吾目①⑪日ゴゼロシ。

Ｓ２

0．５１２５１０２０
Non-dimensionalfrequency,Ｑ

Fig.５ove｢turnth『eshoIdaccelerationscalculatedp
simulatedandexperimented

so1idune：calculatedfromeq.(1)●：egmerimentedand▽：ｓｉｍｕｌａｔｅｄｂｙｅｑｓ（５）ａｎｄ（６）Ｏ７ｚ(b）largeimplusewithshortrisetime

TabIelNormalizedshockmcdeI

ｆｏｒＮｏＢｔｏｌｌ

Fig.４ＭｏｄｅＩｓｏｆｎｏ｢maIizedshockacceleration

ｕ９７５３３３７７７５６７－２５３－７７７に対しては第１近似．

ノノノ算されていることを見出した。そこで、ｓ＜１３３３り、ｓ≧１ではその因子に（s／Ｑ）の因子が加６５４３２２２

かに、ｓ＜１では（Ｑ／s）の因子が支配的であ

幻れれれれれれ兀兀兀兀れれ幻

次に、Fig.５からps、ｓ、Ｑの関係を取り出したところFig.６となった。Ｆｉｇ．６から大ま１９７５３３３を省略した。

刎幻呵れれれが００００汀兀兀
５４３２

特性はそれぞれ番号５，６，７に近いので作図●●■ひ●●●■●●■●の■ｌに対応している。なお、番号３１，３２，３３の

唖鈍感型”趣”加珂０噸卿到樒

た結果がFig.５である。図中の番号はＴａｂＩｅ００００００１１１２３０００に印加し、Runge-Kutta-Gill法で計算し１２３４５６７８９ｍｎｍ犯鋼積sを示す。まず、その衝撃くり（て）を式（１）Ｎｏ．ｓて】西ｑｎｃｏｍｍｅｎｔｓ
0.01
0.01
0.01
0.01

0.01
0.05
0.2
４sine4２

sine2

１sme

●rectangular
0.01

simiIarto50.01

similaｒｔｏ６

0.01

ＳｉｍｉＩａＴｔｏ７

__■陣
ニヨロご⑩、両の騏切:Ｏ￣

ＩＩ＝｜⑫１１－⑰⑩Ｗ』一価堂ｒベ一叶溌凹肖壷Characteristicvalue9ps／Ｑ

回蠅綱。垣・胸Ｓ珊露寓叶乳（］）「ｅ蒲野洲一宮、片。ベ翠鰄アコ②．、丙刺鼓ｒ汁（●巴）・糾騨壷伸寡Ｓ富コトｒべ再診噌勿ⅢでＩ唖山Ｙ邑

議Ｓ噸卯］診凹ゾ②ｅ宙。」－乳Ｓ鐵圏。冑針が、作堂卑己永・鞭Ｆ〈ｍ小作蝿加］ゾ凹診⑭Ｓ

燕博丙醸Ｊペゾ乢口鷺｝戟椰昔鋼苦叶「響⑩叶片ｊ作０鉢萱Ｊくう叡卉Ｊ汽柳削伽・糾露颪詞一℃・吋汰年副一畔鋪璽茸Ｓ茸圃も蕊憾鴬嶺Ｓ

憾う娑蕪制剣艸斡（さくゾコロ・切号⑩弄珍針靱
繭珂州彗⑩Ｓ鰯ラ苞酬騨習昔劉印昔州苦Ｓ爵

露丙宍Ｊべ醗醗茸Ｓ淵鬮萱蕪冊剣脚酬貯Ｊ丙衙Ｊ汁・淵⑩号ご蝋蕊茸興蕊艤諭黙覇亜目費巾〈働戸醗耀団認め費升叫〈歓びｄ繍璽醗耀

言偽岡已ごぐう腓ヴ汁与“図Ｓ所司聖牟行乢口憧｝腓営が。

「百・②苧⑩升頚潴悶通う汁鐡繕蕪甑覇亜目津ヰペグ小巻墹含ｅ②作］肌ｐい］Ｃ側Ｓご崗什醗脚診蝋卯（」ｌ①）Ｓ詮診（⑪ｌ屋）Ｓ詮丙ｅ議興刀垣・ｍｅ院彗苧ｅ逹勿Ⅲマーピ山疋匿

ペゲゾ郡昔苧⑩学識腓昔叡ｐ什菖冊眺営汁ｂが。ｒザＦ“『Ｓ割轡雪ご粥鰍哉汁脈ラトｒ

Ⅱ『ｇ、這醗ＨＰ〈翻戦冒空、昇饗議轍燗斗

蹄叩く、腸汽頓汎“ぜＩ乢刈〈電｝閉軸’・汁議宥黙劔持ｓ籔菫制ｒくち碑ご調剣Ｓ舗璽酎（⑤）

什、穿索鐸四８コ晋一の円ｅ舟Ｊ一価転ぶ苦汁・

耀崗再慧詳凹需ｄ断叡・

鉢サゴ園最高ｅ醜男ごいグ温蔚Ｓ作向いｎ

テＦＪＪ伊診Ｍ①丙醸帥一氏ｒ汁ⅡＯＳ鍬璽扇

マーピ山里帆胖戸州Ｓ餅室叶鮪乳一ｍｕへ超爵帥言で汁齢軍卦鮪乳ごぶ／牢掴藻珀勿憎ｓⅡＪＳ璽帥乳啼蚊紳がｅ萱嵐雌剣針叡・州のＣ亡四『－．コ⑭（⑪）回。。（⑤）○・の⑪⑬。。⑪『、ＢＧ：回牌。．。←。

－い、ムー蕊艤繍璽扇細僻蠅鰯佇汁睦即・凹鮒書Ｓ蹄。識罫

萱ヨ謡甸猷が佇畑苦含が。

小、作戦昌藷ゾサＪ診用昨鶴啼江傍寺蔵鰡理

稗“ｎｓ尚壺挑汽胖戸豐寡ｓ齢窒扇瓠啼苔瓠乳餅塘う汁・醗耀回蕗費竿職制眺伽輌恥汽】駒（四℃己『○涕・酉冒○回ロロ『◎〆。

一●－G°●四￣￣。。⑪哺崗卿「⑪丘一⑮酌。。①薫一、｜の．一切。｛四百℃『。×一『ゴ四｛⑩こ昔鹿舜⑩衝う。

昔①餅欝恥斗叡汁与丙一蜂“囲制醗遼も腓竺斗刈、ミテｒ汁（『召）・卵彗⑩ｅ尚置挑管伊勺醗刺竿鰄ｒ談「一垣．、訂１．＋厚十十野牛（函）勵璃巾Ⅵ￥両氏⑬忠已壗嚢坪ｓ玲塞罰髄LA）．＠．①。①｜昌一。。⑩画ヨ。。○つ。⑩四．Ｑｂ●

ロⅡ四○＋餅‐いトー●一脚②日本包装学会霧Vbk3jVu‘α994）

計に資するため、衝撃に対し剛体の転倒限界●●

０＋=言ｅ８－を与えるシミュレーションを行い、かつ、数(a3）ａｔＯ＝８０種の緩衝材による鋼球の衝撃転勤実験を行

添字の正負は０．に到達した直後及び直前

い、以下の結論を得た。の値を表す。仮定の(4)、（３）により式（１）転倒方程式を数値解析し、衝撃面積≧(al）は1及びく１の場合に分けて簡易転倒限界式を導いた。それらは無次元化振動数Ｑの広ｍＺｄ＝－ｍﾕｭ（cosO＋Llsme）－ｍｇ（sinO-lucosO）（a4）

範囲にわたって、元の微分方程式に対しよここに｡と０の関係はい近似を与える。

Ｍ＋０１＝００，ｄｏ＞０（a5）（２）鋼球の転勤装置と各種緩衝材による衝撃実験において、その衝撃転倒限界値の特であるから、０を①に置き換え、ｅ＞Ｏで６ａ限界式から転倒限界が推定できる。

ｊ

持続時間・衝撃面i積が分かれば、簡易転fi3}くｓｉｎＯ＝ｓiｎ（００－の－の＋Jucos（０．－．－匹）
cos8-cos（８０－⑦－匹）
一旦Ｓｍ（Ｏｏ－ｄ－Ｊｕ）

界式の特性傾向に合致した。よって、衝撃、－－トー』

性は、シミュレーションによる簡易転倒限の近似を与える。本文のパラメータにより次の近似方程式

おわりにＱ２。.＋Sｍ（０ｏ－ｄ－ｊＵ）＝

本論では剛体の転倒限界の内すべり・躍りｕ”のない場合に限定したが、ほかにも、容器底

一丁CＯＳ（０．－．－回）（a7）

形状、積み重ね、内容物の揺れ、台の材質なとなる。さらに本文の無次元化と、仮定

ど転倒に関係する因子は多いと思われる。今(2)を用いて式（１）を得る。

後の研究の発展が期待される。
＜引用文献＞[付録］

１）西原主計、精密機械、４４(2)，１８５（1978）

Ｆｉｇ．１の剛体の重心Ｏの運動方程式は、

２）西原主計、奥田欲、和田充雄、渡辺健朗、製品

:二二潟期二:鴇鵠｝（aｎ

重心０回りの回転運動方程式は、●●

10α＝－ＦI.

科学研究所報告、第87号、５１（1979）
３）西原主計、渡辺健朗、奥田欲、和田充雄、製品科学研究所報告、第80号、１（1977）
４）西原主計､精密工学会誌、５９(2),３２３（1993）
５）木村隼、飯田汲事、地震、６(3)，１２５（1934）(a2）６）日本機械学会、“機械工学便覧"、第５章、ｐ３－３５（1968）

衝突による減衰力は方程式に入れないで次（原稿受付1994年３月１６日）の条件で与える。

（審査受理1994年８月１０日）－２５５－