シミュレーションの作り方 7



その7 伸び縮みするばねを作る。



  物理には「つるまきばね」がよく登場するが、のこぎり型やSin曲線型ではなくより実物に近いもの(右図>を考えてみた。

  ばねの基本仕様として、総巻数N=10回巻、半径R=2 のばねを、座標χ0=0~座標χ1=25の範囲に作ることにする。(右図参考)
  これを各500個のセルでχ、y座標値を表しグラフ化する。

 

 

  そのためにまず、ばねの基本仕様の値とχ、y座標の各500個の値表示する計算式を下図のように入力する。

  基本仕様以外の T,⊿t、t の要素は、χ軸上を+方向に(χ1-χ0)/(NT)で等速運動している中心の周りを、半径R、周期T、角速度=2π/T、で等速円運動している物体の軌跡を 「つるまきばね」の形状と考えて作くるために必要な要素として登場させたものです。(等速χ座標系での等速円運動物体の軌跡)

  L7セルの⊿t は、ばね全体を表示するのに必要な時間はN×Tでこれを500のセル値で表わすので各セル間の時間⊿t=NT/500となり、これが”=$L$5*$L$4/500”となります。
  N4セルの t は、前の時刻にこの⊿tを加えたものなので、”=N3+$L$7”となります。
  P3セルのχは、χ=χ0+(χ1-χ0)t/(NT)+RCos(πー2πt/Tーπ/4)なので、 ”=$L$2++($L$3-$L$2)*N3/($L$4*$L$5)+$L$6*COS(PI()*(1-2*N3/$L$5)-PI()/4)”となります。(
  Q3セルのyは、y=RSin(πー2πt/Tーπ/4)なので、”=$L$6*SIN(PI()*(1-2*N3/$L$5)-PI()/4)”となります。(角度部分は初期位相πーπ/4から開始し右回りのばねを作るため)
忘れていました。L3セルにとりあえず”=M3 ”と計算式を入れておいてください。次にばねを振動させるときにまた書き換えますが・・・・ 
  N4、P3、Q3セルは502までコピーしておきます。

  グラフ用のデータはできたので、右図のように、散布グラフを立ち上げ「系列」で「追加」をクリックしたあと、χ値、y値を範囲指定します。
   (図では3~50となっているが、3~502まででよい)

  このままでは各座標がマークで表示され「つるまきばね」には見えないので、次に系列1の書式設定で見やすくします。 

  右図のように、記号①~⑦の順で「系列1」の書式設定をし直します。

  ①グラフ上をクリックする。
  ②③系列1を選択する。(系列の書式設定画面の出し方はこれにこだわりません) 
  ④系列の書式設定を選ぶ。
  ⑤パターンタブで⑥線を太い実線に指定する。色は何でもよい
  ⑦マーカーは「なし」にチェックする。

これで静止したつるまきばねかできました。(右図は、ばねが見にくかったので、プロットエリアの背景色を白にしてあります。

  右図のようにばねの先端をχ1中心に、振幅 5、周期 0.1 で単振動させる仕掛けを作ります。

  細かい説明はダブりますので省略します。

  振動時刻の黄緑のスピンボタンを押してみてください。

  χ1の座標はP502とQ502なのでこの座標を追加系列のχ、y値に指定し、スタイルを●にすれば、おもりの単振動となります。
  縦方向にばねを作りたいときは、系列1のχとyのデータを反対に設定します。これでOK

  つぎは、任意の向きでばねを振動させてみましょう。

  ばねの向きを変えるには座標変換を行えばよいので、このχ、yの値を原点中心に反時計方向に回転した時の変換座標χ’、y’をS、T列に表示します。
  変換式はχ’=χCosθーySinθ、y’=χSinθ+yCosθなので、下図の様にS3には”=P3*COS($L$16)-Q3*SIN($L$16)

T3  には”=P3*SIN($L$16)+Q3*COS($L$16)”を記入します。コピーの関係で、$マークのの有無には注意してください。

  次に、座標の回転角度をスピンボタンを使って設定します。プロパティは、LinckedCell=M16、Max=360、Min=0です。また、L16セルにはra表示に直す式 "=M16*PI()/180" を入力します。
  続いて、バネ「系列1」のχ値、y値をX'及びy’に切り替えます。まず、グラフ上バネ「系列1」を左クリックします。現れるχ値、y値のDATA範囲の「色囲み線」をドラッグして、それぞれχ’とy’のところに設定をします。これで設定は完了しました。

  このままで回転させるとバネが回転と共にゆがみます。原因は縦横の目盛り設定の不統一によるものです。したがって、上図のようにばねの仕様値と縦横の目盛りの統一とグラフサイズを調整します。

この節で「軸の周りの回転」については一部ふれたので、当初予定していた回転については省略します。今日は4月26日馴れないHP作りはつかれます。
右の図をクリックしてみてください。 (bane1.xls)