RubyGems - num4mechaequ - Versions diffs - 0.0.3 → 0.1.1 - Mend

num4mechaequ 0.0.3 → 0.1.1

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.

Files changed (4) hide show

checksums.yaml CHANGED Viewed

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: aa1f15732e36a08d7b119cc9d638d27ea13182fd1f4a36bac3b4e823996fcffd
-  data.tar.gz: 486c9ee96ffdf6db1a83e23f68bc146b7ea7af7b7298ec13215b71f06078cb91
+  metadata.gz: 77ffa279b2f10648a0cc852d08f6cf23db176f385bc31b8809540a38a2f316ad
+  data.tar.gz: 59f69ea0ee861cc39067f4863191205d756e552d35d8af0f039ce94628a33166
 SHA512:
-  metadata.gz: 4e6936f8107a121a98a2d3f010a99219b8a09bb1f1d769b748300953a42a9d4a6f9780f68a33a838122a88e94a0e962c25d96b87675b15abcf025f3385e09af1
-  data.tar.gz: d1886b86214b6670ce2c5698b1ca4f639a59f5d7cc884be83385b302425d40da6640e57637a47dbe35edaa1a126cdc23da6cb59bcb697ee2499d705d4bd9c88b
+  metadata.gz: aaf91e9df68349bd5c20708a54d0e9d97e312fc904b03ed37076b864d93dc55ec3477e34c872005808df46f176c0988c7421aa54e86879bb27772ef79fc5afcb
+  data.tar.gz: 4a0ef2f3a5ff4fe31a52af87aceba3c593416d881346b3b087bf6a92c75f45f832d9a792aa5d02c2ae6d3a1292491e152892a7217e4cf9e2816020e307812d9e

data/CHANGELOG.md CHANGED Viewed

@@ -2,6 +2,19 @@
 ## Unreleased
+## [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 CHANGED Viewed

@@ -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)
+        # 単振動(simple harmonic motion)
+        # @overload SHM(k, m, t, h0, v0)
         #   @param [double] k バネ定数
         #   @param [double] m 重りの重さ
+        #   @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(k, m, 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,35 @@ 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
     end
 end

metadata CHANGED Viewed

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: num4mechaequ
 version: !ruby/object:Gem::Version
-  version: 0.0.3
+  version: 0.1.1
 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.1
   homepage_uri: http://github.com/siranovel/num4mechaequation
   wiki_uri: https://github.com/siranovel/mydocs/tree/main/num4phyequation/num4mechaequation
 post_install_message: