num4mechaequ 0.0.1 → 0.0.2

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 6a18f5e7ca589b774ab80c59c19c94c24e4901277064c2bbd1de066c2a5eae6c
-  data.tar.gz: 11ecea5860ece5096d7537a9602f2f1c398e6854e4e38db128b101cb708c56eb
+  metadata.gz: 2924e1bd5f10cca8551f34b8d0e360206085627ebb69603c49df86bcf0777675
+  data.tar.gz: 0d74c4d37202716b4e43123dc10cd65c2862256bdbae7a95750f5d475a9faa39
 SHA512:
-  metadata.gz: 0260dc768dc3130a2ffc1173ba2fa4bc87a0c95e4cd386273e328d158e8a79606ab8bdb174b6e967194a944f0b861fa51d87c990a58a02c836d09e4219157ddc
-  data.tar.gz: cced4fc0b159eaf0ab3e6acff661ff9d243dedb129508fe4402883466f33b6dc9b3e3feb09e270e735b594708cb5274cb63ac943ab6e47b901eb016dbaf3e088
+  metadata.gz: 5ce0e327685b176f86aa6f1fdca4d44c6c0144299cb27f6a312d44eece88258841380b54decc0b0bfaebd0fe373ffee2b42b4dc11cce863e0f81e0ea7096a9f4
+  data.tar.gz: 0c35ff365050615fec0cf67c1edc7981a2a6d944a4e00cabcda63f065f124c7b67e1b0ca2a5c8675a7a4c337a806850f3928351748162a15743ce7e69568903f

data/CHANGELOG.md

@@ -2,6 +2,11 @@
 
 ## Unreleased
 
+## [0.0.2] - 2022-06-08
+
+### Added
+- add function of motiEqu4freeFallEqu
+
 ## [0.0.1] - 2022-05-11
 
 ### Fixed

data/lib/num4mechaequ.rb

@@ -4,17 +4,23 @@ require 'num4simdiff'
 # 数値計算による力学方程式の解法するライブラリ
 #
 module Num4MechaEquLib
-    @g = 9.8                           # 重力定数
+    @g = 9.80665                       # 重力加速度
     @t = 0.001                         # 時間刻み
-    @k = 0.0                           # バネ定数
-    @m = 0.0                           # 重りの重さ
+    @w = 0.0
     @springFunc = Proc.new do | n, yi |
         f = []
-        f[0] = yi[1]
-        f[1] = - (@k / @m) * yi[0] - @g
+        f[0] = yi[1]                   # hの傾き値
+        f[1] = @w * yi[0] - @g# vの傾き値
+        next f
+    end
+    @motionFunc = 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)
@@ -25,18 +31,39 @@ module Num4MechaEquLib
         #   @return [hash[]] 0秒から1秒までの位置(h)と速度(v)の値
         #
         def springFreqEqu(k, m, h0, v0)
-            @k = k
-            @m = m
+            @w = Math.sqrt((k / m))
             hvt = []
             yi_1 = []
             yi = [h0, v0]
             0.step(1, @t) { |x|
               yi_1 = Num4SimDiffLib.rungeKuttaMethod(yi, @t, @springFunc)
-              hvt.push({"t" => x, "h" => yi_1[0], "v" => yi_1[1]})
+              hvt.push({"t" => x, 
+                        "h" => yi_1[0], "v" => m * yi_1[1]})
               yi = yi_1
             }
             return hvt
         end
+        #
+        # 自由落下による運動方程式(空気抵抗有り)
+        # @overload motiEqu4freeFallEqu(m, h0, v0)
+        #   @param [double] m 重りの重さ
+        #   @param [double] h0 初期位置値
+        #   @param [double] v0 初期速度
+        #   @return [hash[]] 0秒から1秒までの位置(h)と速度(v)の値
+        #
+        def motiEqu4freeFallEqu(m, h0, v0)
+            hvt = []
+            yi_1 = []
+            yi = [h0, v0]
+            0.step(1, @t) { |x|
+              yi_1 = Num4SimDiffLib.rungeKuttaMethod(yi, @t, @motionFunc)
+              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.0.1
+  version: 0.0.2
 platform: ruby
 authors:
 - siranovel
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2023-05-08 00:00:00.000000000 Z
+date: 2023-06-08 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.1
+  documentation_uri: https://rubydoc.info/gems/num4mechaequ/0.0.2
   homepage_uri: http://github.com/siranovel/num4mechaequation
 post_install_message: 
 rdoc_options: []