num4mechaequ 0.1.2 → 0.1.4

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: ece1c1f8ac4951ed3e9225bc0139f6a332f4ecc22d271cdcb149cd191e28840f
-  data.tar.gz: e9579a5fe2e05b269f0e0aac9c232588bcbe2660268098087311337b694ecd08
+  metadata.gz: 33f7db372d5c6d97cce0dfd47d9040d10178f575bc6394b84322079d96f605e6
+  data.tar.gz: a4ef6a9911a6ae1654c095def96600d283bf69e73657c0f39c8fbeb879e42090
 SHA512:
-  metadata.gz: 6c943c0d56cdb23ba471236f6e772c70f20f0bc96303d9a3138ce7247e27c4ad48a52e63047ad055a2fb03522bdf99eeed6bc4b09b3d104ce19a5276046a8656
-  data.tar.gz: 2428c482848f27714fe833dd741347aa30e29f37191faecf35984a28bdab2c38a4c96787a4d1e7469afb94d927792322570769775bb20c44de330800c250cd7a
+  metadata.gz: 014cabe8283d052239d19853b1591e2e981b2e1cc7c97a98d58d4d47439fc26393ec7f991d8c99628a8a5d448e6f86aa735b515ea5a40c816d339804497cb49b
+  data.tar.gz: de07a537f9ec64ee9046a6c84e758b3949916bc679fb851711ec00d01850b442bb53a6a46a85448cf62d3cb85c54002d4ce4fb6667b002828429f6a61d37c2cf

data/CHANGELOG.md

@@ -2,6 +2,17 @@
 
 ## Unreleased
 
+## [0.1.4] - 2023-10-18
+
+### add
+- add sample of yardoc
+
+## [0.1.3] - 2022-07-20
+
+### Added
+- add function of DHM
+- add function of forcedOscillation
+
 ## [0.1.2] - 2022-07-13
 
 ### Added

data/lib/num4mechaequ.rb

@@ -7,8 +7,10 @@ module Num4MechaEquLib
     @g = 9.80665                       # 重力加速度(m/s2)
     @dt = 0.001                        # 時間刻み
     @w = 0.0
+    @l = 0.0
+    @ft = 0.0                          # 外力
     # 単振動
-    @springFunc = Proc.new do | n, yi |
+    @SHMFunc = Proc.new do | n, yi |
         f = []
         f[0] = yi[1]                   # hの傾き値
         f[1] = - @w * @w * yi[0]       # vの傾き値
@@ -18,7 +20,7 @@ module Num4MechaEquLib
     @motionFunc = Proc.new do | n, yi |
         f = []
         f[0] = yi[1]                   # hの傾き値
-        f[1] = - @g + @w * yi[1]       # vの傾き値
+        f[1] = @w * yi[1] - @g         # vの傾き値
         next f
     end
     # 放物運動
@@ -34,10 +36,24 @@ module Num4MechaEquLib
         f[1] = - @g                    # vの傾き値        
         next f
     end
+    # 減衰振動
+    @DHMFunc = Proc.new do | n, yi |
+        f = []
+        f[0] = yi[1]
+        f[1] = - @l * yi[1] - @w * @w * yi[0] 
+        next f
+    end
+    # 強制振動
+    @forcedFunc = Proc.new do | n, yi |
+        f = []
+        f[0] = yi[1]
+        f[1] = - @w * @w * yi[0] + @ft
+        next f
+    end 
     class << self
 
-        #
         # 単振動(simple harmonic motion)
+        #
         # @overload SHM(m, k, t, h0, v0)
         #   @param [double] m 重りの重さ
         #   @param [double] k バネ定数
@@ -45,6 +61,12 @@ module Num4MechaEquLib
         #   @param [double] h0 初期位置値
         #   @param [double] v0 初期速度
         #   @return [hash[]] 0秒からt秒までの位置(h)と速度(v)の値
+        # @example
+        #   m = 1.0
+        #   k = 1.0
+        #   h0 = 1.0
+        #   v0 = 0.0
+        #   yi_1 = Num4MechaEquLib.SHM(m, k, 2, h0, v0)
         #
         def SHM(m, k, t, h0, v0)
             @w = Math.sqrt((k / m))
@@ -52,15 +74,15 @@ module Num4MechaEquLib
             yi_1 = []
             yi = [h0, v0]
             0.step(t, @dt) { |x|
-              yi_1 = Num4SimDiffLib.rungeKuttaMethod(yi, @dt, @springFunc)
+              yi_1 = Num4SimDiffLib.rungeKuttaMethod(yi, @dt, @SHMFunc)
               hvt.push({"t" => x, 
                         "h" => yi_1[0], "v" => m * yi_1[1]})
               yi = yi_1
             }
             return hvt
         end
-        #
         # 自由落下による運動方程式(空気抵抗有り)
+        #
         # @overload freeFallMotion(m, c, t, h0, v0)
         #   @param [double] m 重りの重さ
         #   @param [double] c 空気抵抗
@@ -68,6 +90,12 @@ module Num4MechaEquLib
         #   @param [double] h0 初期位置値
         #   @param [double] v0 初期速度
         #   @return [hash[]] 0秒からt秒までの位置(h)と速度(v)の値
+        # @example
+        #   m = 1.0
+        #   c = 0.0
+        #   h0 = 1.0
+        #   v0 = 0.0
+        #   yi_1 = Num4MechaEquLib.freeFallMotion(m, c, 2, h0, v0)
         #
         def freeFallMotion(m, c, t, h0, v0)
             @w = c / m
@@ -83,8 +111,8 @@ module Num4MechaEquLib
             return hvt
 
         end
-        #
         # 放物運動
+        #
         # @overload projectileMotion(m, theta, t, h0, v0)
         #   @param [double] m     重りの重さ
         #   @param [double] theta 角度(ラジアン指定)
@@ -92,6 +120,10 @@ module Num4MechaEquLib
         #   @param [double] h0 初期位置値
         #   @param [double] v0 初期速度
         #   @return [hash[]] 0秒からt秒までの位置(h)と速度(v)の値
+        # @example
+        #   m = 1.0
+        #   theta = 5 * Math::PI / 180.0
+        #   yi_1 = Num4MechaEquLib.projectileMotion(m, theta, 2, 0.0, 3.0)
         #
         def projectileMotion(m, theta, t, h0, v0)
             hvt = []
@@ -112,14 +144,20 @@ module Num4MechaEquLib
             }
             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)の値
+        # @example
+        #   m = 1.0
+        #   r  = 3
+        #   w  = 2
+        #   yi_1 = Num4MechaEquLib.UCM(m, r, w, 2)
+        #
         def UCM(m, r, w, t)
             hvt = []
             h = []
@@ -138,8 +176,8 @@ module Num4MechaEquLib
             }
             return hvt
         end
-        #
         # 振り子運動
+        #
         # @overload pendulumMotion(m, l, t, h0, v0)
         #   @param [double] m 重りの重さ
         #   @param [double] l 糸の長さ
@@ -147,9 +185,80 @@ module Num4MechaEquLib
         #   @param [double] h0 初期位置値
         #   @param [double] v0 初期速度
         #   @return [hash[]] 0秒からt秒までの位置(h)と速度(v)の値
+        # @example
+        #   m = 1.0
+        #   l  = 10
+        #   h0 = 1.0
+        #   v0 = 0.0
+        #   yi_1 = Num4MechaEquLib.pendulumMotion(m, l, 2, h0, v0)
+        #
         def pendulumMotion(m, l, t, h0, v0)
             return SHM(m, l, t, h0, v0)
-        end 
+        end
+        # 減衰振動（damped harmonic motion)
+        #
+        # @overload DHM(m, k, b, t, h0, v0)
+        #   @param [double] m 重りの重さ
+        #   @param [double] k 比例定数
+        #   @param [double] b 抵抗力の比例定数
+        #   @param [double] t 時間
+        #   @param [double] h0 初期位置値
+        #   @param [double] v0 初期速度
+        #   @return [hash[]] 0秒からt秒までの位置(h)と速度(v)の値
+        # @example
+        #   m = 1.0
+        #   k = 1.0
+        #   h0 = 1.0
+        #   v0 = 0.0
+        #   yi_1 = Num4MechaEquLib.DHM(m, k, 0.6, 2, h0, v0)
+        #
+        def DHM(m, k, b, t, h0, v0)
+            @l = 2 * m * b
+            @w = Math.sqrt((k / m))
+            hvt = []
+            yi_1 = []
+            yi = [h0, v0]
+            0.step(t, @dt) { |x|
+              yi_1 = Num4SimDiffLib.rungeKuttaMethod(yi, @dt, @DHMFunc)
+              hvt.push({"t" => x, 
+                        "h" => yi_1[0], "v" => m * yi_1[1]})
+              yi = yi_1
+            }
+            return hvt
+        end
+        # 強制振動
+        #
+        # @overload forcedOscillation(m, k, w0, w, t, h0, v0)
+        #   @param [double] m 重りの重さ
+        #   @param [double] k 比例定数
+        #   @param [double] w0 振動系の固有角振動
+        #   @param [double] w 外力の角振動数
+        #   @param [double] t 時間
+        #   @param [double] h0 初期位置値
+        #   @param [double] v0 初期速度
+        #   @return [hash[]] 0秒からt秒までの位置(h)と速度(v)の値
+        # @example
+        #   m = 1.0
+        #   k = 1.0
+        #   f0 = 2.5
+        #   h0 = 1.0
+        #   v0 = 0.0
+        #   yi_1 = Num4MechaEquLib.forcedOscillation(m, k, 0.6, f0, 2, h0, v0)
+        #
+        def forcedOscillation(m, k, w0, w, t, h0, v0)
+            @w = Math.sqrt((k / m))
+            hvt = []
+            yi_1 = []
+            yi = [h0, v0]
+            0.step(t, @dt) { |x|
+              @ft = w / m * Math.cos(w * x)
+              yi_1 = Num4SimDiffLib.rungeKuttaMethod(yi, @dt, @forcedFunc)
+              hvt.push({"t" => x, 
+                        "h" => yi_1[0], "v" => m * yi_1[1]})
+              yi = yi_1
+            }
+            return hvt
+        end
     end
 end
 

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: num4mechaequ
 version: !ruby/object:Gem::Version
-  version: 0.1.2
+  version: 0.1.4
 platform: ruby
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2023-07-10 00:00:00.000000000 Z
+date: 2023-10-18 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.1.2
+  documentation_uri: https://rubydoc.info/gems/num4mechaequ/0.1.4
   homepage_uri: http://github.com/siranovel/num4mechaequation
   wiki_uri: https://github.com/siranovel/mydocs/tree/main/num4phyequation/num4mechaequation
 post_install_message: