<content> <p> Golden ratio min/max </p> </content>

fb200cb2-2945-4542-ba13-fd5031213d92 304 4 5 false false 0 0 decimal

&[DATE] &[TIME] - &[FILENAME]

References:

[1] Timoshenko - Vibration Problems in Engineering pag 342

[2]

[3]

t0 1 time :

ρ 8000 kg m 3^/ * :

ν 0.298 :

E 198000 MPa * :

G E 21 ν + (* / :

κ_hollow 21 ν + (43 ν * + / * : 0.5304

Hollow cylinder cross section

Isect Asect r Isect Asect / (sqrt :

ε κ_hollow :

Isect mass a E Isect mass / * (sqrt :

ω Isect mass β ω 2^Isect mass a 2^/ (4 nthroot :

L Isect Asect α 1 ε / Isect Asect r L / (2^* E G / * :

L Isect Asect γ Isect Asect r L / (2^:

ω Isect Asect mass L A ω Isect mass β L * (4^L Isect Asect α L Isect Asect γ + (* :

ω Isect Asect mass L B ω Isect mass β L * (8^L Isect Asect α (L Isect Asect γ (- (2^* 4 ω Isect mass β L * (4^* (+ (sqrt :

ω Isect Asect mass L η 2 sqrt 2 L * / ω Isect Asect mass L A - ω Isect Asect mass L B + (sqrt * :

ω Isect Asect mass L θ 2 sqrt 2 L * / ω Isect Asect mass L A ω Isect Asect mass L B + (sqrt * :

ω Isect Asect mass L σ1 1 L / ω Isect Asect mass L θ L * L Isect Asect α ω Isect Asect mass L θ L * / ω Isect mass β L * (4^* - (* :

ω Isect Asect mass L σ2 1 L / ω Isect Asect mass L η L * L Isect Asect α ω Isect Asect mass L η L * / ω Isect mass β L * (4^* + (* :

Isect k E Isect * :

ω Isect Asect mass L σ11 ω Isect Asect mass L θ 3^L Isect Asect α L 2^* ω Isect mass β 4^* ω Isect Asect mass L θ * - L Isect Asect γ L 2^* ω Isect mass β 4^* ω Isect Asect mass L σ1 * - (:

ω Isect Asect mass L σ21 ω Isect Asect mass L η 3^L Isect Asect α L 2^* ω Isect mass β (4^* ω Isect Asect mass L η * + L Isect Asect γ L 2^* ω Isect mass β 4^* ω Isect Asect mass L σ1 * + (:

ω Isect Asect mass L σ12 ω Isect Asect mass L θ (2^L Isect Asect α L 2^* ω Isect mass β (4^* - (:

ω Isect Asect mass L σ21 ω Isect Asect mass L η (3^L Isect Asect α L 2^* ω Isect mass β (4^* ω Isect Asect mass L η * + L Isect Asect γ L 2^* ω Isect mass β (4^* ω Isect Asect mass L σ2 * + (:

ω Isect Asect mass L σ22 ω Isect Asect mass L η (2^L Isect Asect α L 2^* ω Isect mass β (4^* + (:

ω Isect Asect mass L φ11 1 ω Isect Asect mass L σ12 ω Isect Asect mass L σ22 + / ω Isect Asect mass L σ22 ω Isect Asect mass L θ L * cos * ω Isect Asect mass L σ12 ω Isect Asect mass L η L * cosh * + (* :

ω Isect Asect mass L φ12 1 ω Isect Asect mass L σ1 ω Isect Asect mass L σ21 * ω Isect Asect mass L σ2 ω Isect Asect mass L σ11 * + / ω Isect Asect mass L σ21 ω Isect Asect mass L θ L * (sin * ω Isect Asect mass L σ11 ω Isect Asect mass L η L * (sinh * + (* :

ω Isect Asect mass L φ13 1 Isect k (ω Isect Asect mass L σ1 ω Isect Asect mass L σ21 * ω Isect Asect mass L σ2 ω Isect Asect mass L σ11 * + (* / ω Isect Asect mass L σ2 - (ω Isect Asect mass L θ L * (sin * ω Isect Asect mass L σ1 ω Isect Asect mass L η L * (sinh * + (* :

ω Isect Asect mass L φ14 1 Isect k (ω Isect Asect mass L σ12 ω Isect Asect mass L σ22 + (* / ω Isect Asect mass L θ L * cos - ω Isect Asect mass L η L * (cosh + (* :

ω Isect Asect mass L φ21 1 ω Isect Asect mass L σ12 ω Isect Asect mass L σ22 + / ω Isect Asect mass L σ1 ω Isect Asect mass L σ22 ω Isect Asect mass L θ L * (sin * * ω Isect Asect mass L σ2 ω Isect Asect mass L σ12 * ω Isect Asect mass L η L * (sinh * - (* :

ω Isect Asect mass L φ22 1 ω Isect Asect mass L σ1 ω Isect Asect mass L σ21 * ω Isect Asect mass L σ2 ω Isect Asect mass L σ11 * + / ω Isect Asect mass L σ2 ω Isect Asect mass L σ21 ω Isect Asect mass L θ L * ((cos * * ω Isect Asect mass L σ2 ω Isect Asect mass L σ11 * ω Isect Asect mass L η L * ((cosh * + (* :

ω Isect Asect mass L φ23 ω Isect Asect mass L σ1 ω Isect Asect mass L σ2 * Isect k (ω Isect Asect mass L σ1 ω Isect Asect mass L σ21 * ω Isect Asect mass L σ2 ω Isect Asect mass L σ11 * + (* / ω Isect Asect mass L θ L * cos - ω Isect Asect mass L η L * (cosh + (* :

ω Isect Asect mass L φ24 1 Isect k (ω Isect Asect mass L σ12 ω Isect Asect mass L σ22 + (* / ω Isect Asect mass L σ1 ω Isect Asect mass L θ L * (sin * ω Isect Asect mass L σ2 ω Isect Asect mass L η L * (sinh * + (* :

ω Isect Asect mass L φ31 Isect k (ω Isect Asect mass L σ12 ω Isect Asect mass L σ22 + (/ ω Isect Asect mass L σ11 ω Isect Asect mass L σ22 * ω Isect Asect mass L θ L * (sin * ω Isect Asect mass L σ12 ω Isect Asect mass L σ21 * ω Isect Asect mass L η L * sinh * + (* :

ω Isect Asect mass L φ32 Isect k ω Isect Asect mass L σ11 ω Isect Asect mass L σ21 * * - (ω Isect Asect mass L σ1 ω Isect Asect mass L σ21 * ω Isect Asect mass L σ2 ω Isect Asect mass L σ11 * + (/ ω Isect Asect mass L θ L * cos ω Isect Asect mass L η L * (cosh - (* :

ω Isect Asect mass L φ33 1 ω Isect Asect mass L σ1 ω Isect Asect mass L σ21 * ω Isect Asect mass L σ2 ω Isect Asect mass L σ11 * + (/ ω Isect Asect mass L σ2 ω Isect Asect mass L σ21 * ω Isect Asect mass L θ L * cos * ω Isect Asect mass L σ1 ω Isect Asect mass L σ21 * ω Isect Asect mass L η L * ((cosh * + (* :

ω Isect Asect mass L φ34 1 ω Isect Asect mass L σ12 ω Isect Asect mass L σ22 + / ω Isect Asect mass L σ11 - (ω Isect Asect mass L θ L * sin * ω Isect Asect mass L σ21 ω Isect Asect mass L η L * (sinh * + (* :

ω Isect Asect mass L φ41 Isect k (ω Isect Asect mass L σ12 * ω Isect Asect mass L σ22 * - ω Isect Asect mass L σ12 ω Isect Asect mass L σ22 + (/ ω Isect Asect mass L θ L * cos ω Isect Asect mass L η L * (cosh - (* :

ω Isect Asect mass L φ42 Isect k (- ω Isect Asect mass L σ1 ω Isect Asect mass L σ21 * ω Isect Asect mass L σ2 ω Isect Asect mass L σ11 * + (/ ω Isect Asect mass L σ12 ω Isect Asect mass L σ21 ω Isect Asect mass L θ L * (sin * * ω Isect Asect mass L σ11 ω Isect Asect mass L σ22 * ω Isect Asect mass L η L * (sinh * - (* :

ω Isect Asect mass L φ43 1 ω Isect Asect mass L σ1 ω Isect Asect mass L σ21 * ω Isect Asect mass L σ2 ω Isect Asect mass L σ11 * + (/ ω Isect Asect mass L σ2 ω Isect Asect mass L σ12 ω Isect Asect mass L θ L * (sin * * ω Isect Asect mass L σ1 ω Isect Asect mass L σ22 * ω Isect Asect mass L η L * (sinh * + (* :

ω Isect Asect mass L φ44 1 ω Isect Asect mass L σ12 ω Isect Asect mass L σ22 + (/ ω Isect Asect mass L σ12 ω Isect Asect mass L θ L * ((cos * ω Isect Asect mass L σ22 ω Isect Asect mass L η L * ((cosh * + (* :

D1 0.40 m * :

D2 0.20 m * :

D3 0.14 m * :

d1 0.16 m * :

d2 0.10 m * :

d3 0.04 m * :

L1 2 m * :

L2 4 m * :

L3 2 m * :

D d f_Isect π 64 / D (4^d (4^- (* :

D d f_Asect π D 2 / (2^d 2 / (2^- (* :

ρ D d f_mass ρ π D 2 / (2^d 2 / (2^- (* * :

L L1 :

mass ρ D1 d1 f_mass : 844.4601

Isect D1 d1 f_Isect : 0.0012

Asect D1 d1 f_Asect : 0.1056

ω Φ1 ω Isect Asect mass L φ11 ω Isect Asect mass L φ12 ω Isect Asect mass L φ13 ω Isect Asect mass L φ14 ω Isect Asect mass L φ21 ω Isect Asect mass L φ22 ω Isect Asect mass L φ23 ω Isect Asect mass L φ24 ω Isect Asect mass L φ31 ω Isect Asect mass L φ32 ω Isect Asect mass L φ33 ω Isect Asect mass L φ34 ω Isect Asect mass L φ41 ω Isect Asect mass L φ42 ω Isect Asect mass L φ43 ω Isect Asect mass L φ44 44 mat :

L L2 :

mass ρ D2 d2 f_mass :

Isect D2 d2 f_Isect :

Asect D2 d2 f_Asect :

ω Φ2 ω Isect Asect mass L φ11 ω Isect Asect mass L φ12 ω Isect Asect mass L φ13 ω Isect Asect mass L φ14 ω Isect Asect mass L φ21 ω Isect Asect mass L φ22 ω Isect Asect mass L φ23 ω Isect Asect mass L φ24 ω Isect Asect mass L φ31 ω Isect Asect mass L φ32 ω Isect Asect mass L φ33 ω Isect Asect mass L φ34 ω Isect Asect mass L φ41 ω Isect Asect mass L φ42 ω Isect Asect mass L φ43 ω Isect Asect mass L φ44 44 mat :

L L3 :

mass ρ D3 d3 f_mass :

Isect D3 d3 f_Isect :

Asect D3 d3 f_Asect :

ω Φ3 ω Isect Asect mass L φ11 ω Isect Asect mass L φ12 ω Isect Asect mass L φ13 ω Isect Asect mass L φ14 ω Isect Asect mass L φ21 ω Isect Asect mass L φ22 ω Isect Asect mass L φ23 ω Isect Asect mass L φ24 ω Isect Asect mass L φ31 ω Isect Asect mass L φ32 ω Isect Asect mass L φ33 ω Isect Asect mass L φ34 ω Isect Asect mass L φ41 ω Isect Asect mass L φ42 ω Isect Asect mass L φ43 ω Isect Asect mass L φ44 44 mat :

ω Φ ω Φ3 ω Φ2 * ω Φ1 * :

ω detZ ω Φ 33 el ω Φ 44 el * ω Φ 34 el ω Φ 43 el * - :

Ω ω Hz * detZ ω (0 Hz * 60000 Hz * solve (1 * Hz * :

Ω #

13700 Hz * detZ 3.7884 105^* - 2.4706 1012 -^* i * -

32800 Hz * detZ 1.8988 106^* - 2.6708 1011 -^* i * +

55000 Hz * detZ 9.8205 105^* 1.4784 109 -^* i * -

Ω 1 el detZ #

Ω 2 el detZ #

Ω 3 el detZ #

13800 Hz * detZ 89249.5758

32900 Hz * detZ 1.1857 106^* 2.118 1010 -^* i * +

55500 Hz * detZ 2.5669 106^* 1.6664 1010 -^* i * +

s Z s detZ 108^/ :

64 Z #

roots x Z x 165 - 165 solve :

r roots roots 0 * augment :

plot x Z axes r 31 sys :

x y axes x y / - :

1 time t0 - 4206 s *

plot

roots #

1 time t0 - 13 s *

x f x Z :

on failure , increase TOL TOL 1012 -^:

rootsDichotomy local command line style x0 x1 TOL D a x0 : b x1 : 12 mat xn b a + (2 / : ε TOL : 12 mat n 0 : xn f abs ε > xn f a f * 0 ≤ b xn : a xn : if xn b a + (2 / : n n 1 + : 31 line while xn 21 line 41 line :

r 163 - 160 - TOL D 150 - 140 - TOL D 130 - 120 - TOL D 105 - 100 - TOL D 80 - 75 - TOL D 60 - 50 - TOL D 35 - 30 - TOL D 15 - 12 - TOL D 81 mat : 161.500000000000 - 145.000000000000 - 125.000000000000 - 102.500000000000 - 77.500000000000 - 55.000000000000 - 32.500000000000 - 13.500000000000 - 81 mat

all r r - reverse stack :

all 161.5 - 145 - 125 - 102.5 - 77.5 - 55 - 32.5 - 13.5 - 13.5 32.5 55 77.5 102.5 125145 161.5 161 mat

1 time t0 - 13 s *

all all all 0 * augment :

plot x Z axes all 31 sys :

plot

Golden-Ratio search [min/max]

<content> <p> Golden ratio min/max </p> </content>

Code adapted from Valery Ochov Mathcad 11 f a b test TOL Extrema a b G 32 / b * 12 / a * - 12 / 5 sqrt * a b - (* + : y1 y2 ≥ (max in range ≡ y1 y2 ≤ (min in range ≡ x1 a b G eval : y1 x1 f eval : x2 b a G eval : y2 x2 f eval : 41 line iter 0 : iter 1000 < a b - abs TOL ≤ break test b x2 : x2 x1 : y2 y1 : x1 a b G eval : y1 x1 f eval : 51 line a x1 : x1 x2 : y1 y2 : x2 b a G eval : y2 x2 f eval : 51 line if iter iter 1 + : 21 line if while a b + 2 / 71 line :

a b + 2 / iter 21 mat

y1 y2 ≥ (max in range ≡ y1 y2 ≤ (min in range ≡ 21 sys

x f x Z eval :

ε 1012 -^:

162 - 147 - 125 - 102 - 79 - 55 - 33 - 14 - 14335579102125147162161 mat

X f 160 - 150 - y1 y2 ≥ ε Extrema f 140 - 120 - y1 y2 ≤ ε Extrema f 120 - 100 - y1 y2 ≥ ε Extrema 31 mat :

X X f vectorize 12 mat #

1 time t0 - 26 s *