<content> CreateMesh,γ </content>

24648be5-4b77-4c5c-9289-0c4715e7d06a 6 4 5 false false 0 decimal

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

Schatz linkage

Pin,CreateMesh,γ a Ω1 1000 a cos a sin - 0 a sin a cos 33 mat : a Ω2 a cos 0 a sin 010 a sin - 0 a cos 33 mat : a Ωz a cos a sin - 0 a sin a cos 000133 mat : 13 mat u v C 0.08 v sin * 0.08 v cos * 0.03 u * 31 mat : j C 0 π 02 π * 88 CreateMesh eval : L 3 sqrt : O1 L 0013 mat : s 1 : 15 mat u v C1 0.02 v sin * 0.02 v cos * 0.3 u * 31 mat : r C1 0 π 02 π * 88 CreateMesh eval : 12 mat u v C2 0.02 v sin * 0.02 v cos * 0.5 - (1 0.25 u * - (* 31 mat : m C2 0 π 02 π * 88 CreateMesh eval : 12 mat pin 0 s - 00 s 023 mat : h L 0 0.5 - 00 0.5 - 23 mat : γ 0.8232 - 0.4194 - 0.3827 0.5677 0.6187 - 0.543 0.009 0.6643 0.7474 33 mat 11 * : 13 mat 51 line

<content> construction of a cylinder by two points on its axis </content>

Construction of a cylinder by two points on its axis axis t el cil v1 axis t el 1 row axis t el 2 row - : ev v1 v1 norme / : 12 mat d ev 2 el (2^ev 3 el (2^+ sqrt eval : α2 ev 3 el d 1015 -^+ / acos eval : 12 mat β2 d acos eval : ψ ev 1 el 0 ≥ (ev 2 el 0 ≥ (& β2 - Ω2 α2 Ω1 * eval ev 1 el 0 ≤ (ev 2 el 0 ≥ (& β2 Ω2 α2 Ω1 * eval ev 1 el 0 ≤ (ev 2 el 0 ≤ (& β2 Ω2 α2 - Ω1 * eval β2 - Ω2 α2 - Ω1 * eval if if if : el ψ * 51 line :

<content> Draghilev Method </content>

u v C1 0.02 v sin * 0.02 v cos * 0.3 u * 31 mat : M1 C1 0 π 02 π * 88 CreateMesh eval : 12 mat

Draghilev Method tmax D a b c d e S a b c d e submatrix : n f rows 1 + : k 1 n range : u k el x k el : 13 mat b 1 - (f u Jacobian 1 n 1 - n n S * : A f u Jacobian 1 n 1 - 1 n 1 - S : Δ 1 el b A 1 n 1 - 2 n 1 - S augment det - : Δ n 1 - el A 1 n 1 - 1 n 2 - S b augment det - : Δ n el A det - : 12 mat k 2 n 2 - (range : Δ k el A 1 n 1 - 1 k 1 - S b A 1 n 1 - k 1 + n 1 - S augment det - : 41 line k 1 n range : D2 k el Δ k el Δ 1 el 2^Δ 2 el (2^+ sqrt / : t x D1 D2 : 13 mat X0 00110116 mat transpose x f f : al_nleqsolve : X0 0 tmax N t x D1 rkfixed 81 line :

f x 1 el (2^x 2 el (2^+ x 3 el (2^+ 1 - x 1 el x 4 el - (2^x 2 el x 5 el - (2^+ x 3 el x 6 el - (2^+ 1 - x 4 el L - (2^x 5 el (2^+ x 6 el (2^+ 1 - x 3 el x 6 el * 12 / - L x 4 el + 2 / x 5 el 2 / x 6 el 2 / 1 L 001 x 1 el x 2 el x 3 el 1 x 4 el x 5 el x 6 el 144 mat det 51 mat :

tmax 3.4 : Δt 0.05 : N tmax Δt / : B1 tmax D 1 N 2 f rows 2 + submatrix : 14 mat x1 B1 1 col : x2 B1 2 col : x3 B1 3 col : x4 B1 4 col : x5 B1 5 col : x6 B1 6 col : 16 mat a x1 x2 x3 augment eval : c x4 x5 x6 augment eval : 12 mat 31 line

:Estimated length of middle link : [L= 1]

k 1 N range :

L x1 x4 - (2^vectorize x2 x5 - (2^vectorize + x3 x6 - (2^vectorize + :

Error of changing the length of the middle link in one cycle:

ε L 1 - :

Error of changing the length of the middle

k Δt * ε 10 6 ^ * augment

ε 1.4226 109 -^* - 3.2267 109 -^* - 9.9705 109 -^* - 2.2575 108 -^* - 4.2899 108 -^* - 7.4238 108 -^* - 1.2219 107 -^* - 1.9587 107 -^* - 3.0825 107 -^* - 4.6843 107 -^* - 6.4494 107 -^* - 7.2065 107 -^* - 6.4674 107 -^* - 5.0234 107 -^* - 3.3529 107 -^* - 151 mat