num4mechaequ 0.0.3 → 0.1.2

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: aa1f15732e36a08d7b119cc9d638d27ea13182fd1f4a36bac3b4e823996fcffd
-  data.tar.gz: 486c9ee96ffdf6db1a83e23f68bc146b7ea7af7b7298ec13215b71f06078cb91
+  metadata.gz: ece1c1f8ac4951ed3e9225bc0139f6a332f4ecc22d271cdcb149cd191e28840f
+  data.tar.gz: e9579a5fe2e05b269f0e0aac9c232588bcbe2660268098087311337b694ecd08
 SHA512:
-  metadata.gz: 4e6936f8107a121a98a2d3f010a99219b8a09bb1f1d769b748300953a42a9d4a6f9780f68a33a838122a88e94a0e962c25d96b87675b15abcf025f3385e09af1
-  data.tar.gz: d1886b86214b6670ce2c5698b1ca4f639a59f5d7cc884be83385b302425d40da6640e57637a47dbe35edaa1a126cdc23da6cb59bcb697ee2499d705d4bd9c88b
+  metadata.gz: 6c943c0d56cdb23ba471236f6e772c70f20f0bc96303d9a3138ce7247e27c4ad48a52e63047ad055a2fb03522bdf99eeed6bc4b09b3d104ce19a5276046a8656
+  data.tar.gz: 2428c482848f27714fe833dd741347aa30e29f37191faecf35984a28bdab2c38a4c96787a4d1e7469afb94d927792322570769775bb20c44de330800c250cd7a

data/CHANGELOG.md

@@ -2,6 +2,30 @@
 
 ## Unreleased
 
+## [0.1.2] - 2022-07-13
+
+### Added
+- add function of UCM
+- add function of pendulumMotion
+
+### Fixed
+- change order of parameter of SHM
+
+  SHM(k, m, t, h0, v0) -> SHM(m, k, t, h0, v0)
+
+## [0.1.1] - 2022-07-10
+
+### Changed
+- change of function name
+  * motiEqu4freeFallEqu -> freeFallMotion
+  * springFreqEqu       -> SHM
+
+### Fixed
+- fix appointed time. (from 0 to n)
+
+### Added
+- add function of projectileMotion
+
 ## [0.0.3] - 2022-07-03
 
 ### Changed

data/lib/num4mechaequ.rb

@@ -4,39 +4,55 @@ require 'num4simdiff'
 # 数値計算による力学方程式の解法するライブラリ
 #
 module Num4MechaEquLib
-    @g = 9.80665                       # 重力加速度
-    @t = 0.001                         # 時間刻み
+    @g = 9.80665                       # 重力加速度(m/s2)
+    @dt = 0.001                        # 時間刻み
     @w = 0.0
+    # 単振動
     @springFunc = Proc.new do | n, yi |
         f = []
         f[0] = yi[1]                   # hの傾き値
-        f[1] = @w * yi[0] - @g# vの傾き値
+        f[1] = - @w * @w * yi[0]       # vの傾き値
         next f
     end
+    # 自由落下運動
     @motionFunc = Proc.new do | n, yi |
         f = []
         f[0] = yi[1]                   # hの傾き値
-        f[1] = - @g                    # vの傾き値
+        f[1] = - @g + @w * yi[1]       # vの傾き値
+        next f
+    end
+    # 放物運動
+    @projectileXFunc = Proc.new do | n, yi |
+        f = []
+        f[0] = yi[1]                   # hの傾き値
+        f[1] = 0                       # vの傾き値        
+        next f
+    end
+    @projectileYFunc = Proc.new do | n, yi |
+        f = []
+        f[0] = yi[1]                   # hの傾き値
+        f[1] = - @g                    # vの傾き値        
         next f
     end
     class << self
 
         #
-        # バネの固有振動
-        # @overload springFreqEqu(k, m, h0, v0)
-        #   @param [double] k バネ定数
+        # 単振動(simple harmonic motion)
+        # @overload SHM(m, k, t, h0, v0)
         #   @param [double] m 重りの重さ
+        #   @param [double] k バネ定数
+        #   @param [double] t 時間
         #   @param [double] h0 初期位置値
         #   @param [double] v0 初期速度
-        #   @return [hash[]] 0秒から1秒までの位置(h)と速度(v)の値
+        #   @return [hash[]] 0秒からt秒までの位置(h)と速度(v)の値
         #
-        def springFreqEqu(k, m, h0, v0)
+        def SHM(m, k, t, h0, v0)
             @w = Math.sqrt((k / m))
             hvt = []
             yi_1 = []
             yi = [h0, v0]
-            0.step(1, @t) { |x|
-              yi_1 = Num4SimDiffLib.rungeKuttaMethod(yi, @t, @springFunc)
+            0.step(t, @dt) { |x|
+              yi_1 = Num4SimDiffLib.rungeKuttaMethod(yi, @dt, @springFunc)
               hvt.push({"t" => x, 
                         "h" => yi_1[0], "v" => m * yi_1[1]})
               yi = yi_1
@@ -45,18 +61,21 @@ module Num4MechaEquLib
         end
         #
         # 自由落下による運動方程式(空気抵抗有り)
-        # @overload motiEqu4freeFallEqu(m, h0, v0)
+        # @overload freeFallMotion(m, c, t, h0, v0)
         #   @param [double] m 重りの重さ
+        #   @param [double] c 空気抵抗
+        #   @param [double] t 時間
         #   @param [double] h0 初期位置値
         #   @param [double] v0 初期速度
-        #   @return [hash[]] 0秒から1秒までの位置(h)と速度(v)の値
+        #   @return [hash[]] 0秒からt秒までの位置(h)と速度(v)の値
         #
-        def motiEqu4freeFallEqu(m, h0, v0)
+        def freeFallMotion(m, c, t, h0, v0)
+            @w = c / m
             hvt = []
             yi_1 = []
             yi = [h0, v0]
-            0.step(1, @t) { |x|
-              yi_1 = Num4SimDiffLib.rungeKuttaMethod(yi, @t, @motionFunc)
+            0.step(t, @dt) { |x|
+              yi_1 = Num4SimDiffLib.rungeKuttaMethod(yi, @dt, @motionFunc)
               hvt.push({"t" => x, 
                         "h" => yi_1[0], "v" => m * yi_1[1]})
               yi = yi_1
@@ -64,6 +83,73 @@ module Num4MechaEquLib
             return hvt
 
         end
+        #
+        # 放物運動
+        # @overload projectileMotion(m, theta, t, h0, v0)
+        #   @param [double] m     重りの重さ
+        #   @param [double] theta 角度(ラジアン指定)
+        #   @param [double] t 時間
+        #   @param [double] h0 初期位置値
+        #   @param [double] v0 初期速度
+        #   @return [hash[]] 0秒からt秒までの位置(h)と速度(v)の値
+        #
+        def projectileMotion(m, theta, t, h0, v0)
+            hvt = []
+            yxi_1 = []
+            yyi_1 = []
+            yxi = [h0, v0 * Math.cos(theta)]
+            yyi = [h0, v0 * Math.sin(theta)]
+            0.step(t, @dt) { |x|
+              yxi_1 = Num4SimDiffLib.rungeKuttaMethod(yxi, @dt, @projectileXFunc)
+              yyi_1 = Num4SimDiffLib.rungeKuttaMethod(yyi, @dt, @projectileYFunc)
+              hvt.push({"t" => x, 
+                        "x" => {"h" => yxi_1[0], "v" => m * yxi_1[1]},
+                        "y" => {"h" => yyi_1[0], "v" => m * yyi_1[1]},
+                       }
+                      )
+              yxi = yxi_1
+              yyi = yyi_1
+            }
+            return hvt
+        end
+        #
+        # 等速円運動(Uniform Circular motion)
+        # @overload UCM(m, r, w, t)
+        #   @param [double] m 重りの重さ
+        #   @param [double] r 半径
+        #   @param [double] w 角速度
+        #   @param [double] t 時間
+        #   @return [hash[]] 0秒からt秒までの位置(h)と速度(v)の値
+        def UCM(m, r, w, t)
+            hvt = []
+            h = []
+            v = []
+            0.step(t, @dt) { |x|
+              sinwt = Math.sin(w * x)
+              coswt = Math.cos(w * x)
+              h = [r * coswt, r * sinwt]
+              v = [r * w * -sinwt, r * w * coswt]
+              a = [-r * w * w * coswt, -r * w * w * sinwt]
+              hvt.push({"t" => x, 
+                        "x" => {"h" => h[0], "v" => m * v[0]},
+                        "y" => {"h" => h[1], "v" => m * v[1]},
+                       }
+                      )
+            }
+            return hvt
+        end
+        #
+        # 振り子運動
+        # @overload pendulumMotion(m, l, t, h0, v0)
+        #   @param [double] m 重りの重さ
+        #   @param [double] l 糸の長さ
+        #   @param [double] t 時間
+        #   @param [double] h0 初期位置値
+        #   @param [double] v0 初期速度
+        #   @return [hash[]] 0秒からt秒までの位置(h)と速度(v)の値
+        def pendulumMotion(m, l, t, h0, v0)
+            return SHM(m, l, t, h0, v0)
+        end 
     end
 end
 

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: num4mechaequ
 version: !ruby/object:Gem::Version
-  version: 0.0.3
+  version: 0.1.2
 platform: ruby
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2023-07-03 00:00:00.000000000 Z
+date: 2023-07-10 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: num4simdiff
@@ -45,7 +45,7 @@ licenses:
 - MIT
 metadata:
   changelog_uri: http://github.com/siranovel/num4mechaequation/blob/main/CHANGELOG.md
-  documentation_uri: https://rubydoc.info/gems/num4mechaequ/0.0.3
+  documentation_uri: https://rubydoc.info/gems/num4mechaequ/0.1.2
   homepage_uri: http://github.com/siranovel/num4mechaequation
   wiki_uri: https://github.com/siranovel/mydocs/tree/main/num4phyequation/num4mechaequation
 post_install_message: