ばねます日記を読みました １
私はＭＡＴＬＡＢをよくつかうのですが、ＳＣＩＬＡＢはＭＡＴＬＡＢと異なる部分も多いようです。
ＳＣＩＬＡＢを使用したブログでばねます日記というのを見つけてあります。
だれが何のために書いたのかという情報は無いようです。
しかし、参考文献の筆頭に挙げられているのが．．．
カヤバ工業の「自動車のサスペンション」だったりします (#^_^#)
いちばん基本的なのは１自由度振動系のシミュレーションでしょう。
ばねます日記ではＳＣＩＬＡＢに付属するＸＣＯＳというアプリケーションを使ってシミュレーションを行っています。
ＭＡＴＬＡＢに付属するＳＩＭＵＬＩＮＫと同様に、ＸＣＯＳもブロック線図を描いて使用するソフトウェアですがＳＩＭＵＬＩＮＫほど使いやすくはないと、私は思っています。
というか、ほとんど使った経験がありません。
ばねます日記には以下のように記されています。
周波数を徐々に変化させた正弦波をファイルから読み込むというのは、
とりあえず、難易度が高いと感じます ( ￣ロ￣)ﾉ
慣れていないソフトなので確認しながら進んでみます。
まずは．．．
正弦波をオシロスコープに入れます。
再生ボタンのような三角ボタンがあったので押してみると、シミュレーションが実行されました。
振幅０．２４にはしておいて．．．周波数１Ｈｚの結果。
つぎに１０Ｈｚ
正しいと思います ( ￣ロ￣)ﾉ
つぎに、微分いきます。
意外に苦労していたりします (￣－￣)
線の分岐が難しいのですが．．．
まず、線を出すときは入力側からドラッグなんですが．．．
曲がった線にしたい場合はマウスボタンを解除すると角ができるような．．．
とりあえず、線の分岐のところまで来たら右クリックにすると、垂直とか水平とかのメニューがでてくるので．．．「最適」にしとけば大丈夫のような気がします (￣－￣)
で、微分の結果１Ｈｚ
サインを微分したらコサインになっているので正しいと思います。
１０Ｈｚの場合
サインなのかコサインなのかはわかりにくくなりましたが振幅が１０倍になっているので．．．正しいと思います。
多少ですが自信もついてきたので、１自由度振動系のブロック図を描きます。
さほど苦労しませんでした。
位置と速度をフィードバックする線がありますが、右クリックの「最適」にするとうまく描けませんでした。
線を選択すると角の部分が黄緑色の□になるので、その□の部分をつまんでドラッグするように意識するとうまくいくような気がしました （´Д｀；）
ところで、ＳＣＩＬＡＢは６４ビット版と３２ビット版の両方をインストールしました。
それで、この段階で６４ビット版から３２ビット版に切り替えています。
上記のサインコサインのあたりはまでは６４ビット版のＳＣＩＬＡＢ６．１．０でした。
変数はワークスペースの方に入れておきました。
質量３００ｋｇ、ばね定数３０００Ｎ／ｍ、減衰係数３０００Ｎｓ／ｍ
．．．ここで間違いに気が付きました。
ばねます日記では入力信号の振幅を０．０２４ｍとしていました。
さきほどは０．２４ｍの振幅になっていたので、ここで修正します。
シミュレーション結果１Ｈｚ → 0.32
１０Ｈｚ → 3.2
さて、ここでも．．．いくつかの間違いに気がつきました。
まず、周波数の設定がＨｚでなくｒａｄ／ｓでした （；´Д｀）
つまり、ここまでの１Ｈｚと１０Ｈｚは０．３２Ｈｚと３．２Ｈｚの間違いです。
ところで、ばねます日記は何を言いたいのか？と考えて (￣－￣)
共振点と減衰を言いたいのだと感じました ( ・∀・)
そこで、あらためて、１．４Ｈｚと１０Ｈｚを確認します。
数値計算の刻み幅が０．１秒と大きかったので０．００１秒（１ＫＨｚ）に修正しました。
３０秒間の計算は長いので５秒間にします。
共振点付近の１．４Ｈｚ
振幅が０．０３５ｍより少し小さいぐらいだと思うので良いと思います。
１０Ｈｚの場合
細かい数値は読み取れませんが、振幅が１／６ぐらいまで落ちてきているということで．．．納得できる結果ではないかと思います v(￣ー￣)v