search

dankelley / mooring

Design and analysis of oceanographic moorings
https://dankelley.github.io/mooring
4 stars 1 forks source link

problem with 1" neutral wire #40

Closed dankelley closed 3 years ago

dankelley commented 3 years ago

author: kelley date: 2021-06-21 output: "reprex::reprex\_document" title: satin-urial_reprex.R

library(mooring)
a <- anchor(depth=110)
W <- wire('1in Neutral Wire',length=100)
w <- wire('1/2in Neutral Wire',length=100)
f <- float("HMB 20")
m <- discretise(mooring(a,w,f))
M <- discretise(mooring(a,W,f))
par(mfrow=c(1,2))
plot(knockdown(m,0.5),which="knockdown",title="thin wire")
plot(knockdown(M,0.5),which="knockdown",title="thick wire (should deform more)")

w$CD
#> [1] 1.5
W$CD
#> [1] 0

Created on 2021-06-21 by the reprex package (v2.0.0)

dankelley commented 3 years ago

The problem is that the CD value is zero. I'll fix that tomorrow. I'll also look into whether the radius ratio should be 2 -- perhaps the name is not quite right (this seems very likely) or there is a rounding problem in Dewey's matlab file (also not too unlikely).

dankelley commented 3 years ago

Fixed in "main" commit a6a7628fdd8626ae1150084879ee526d5a0d93e2. I wrote code to detect zero CD in Dewey's file (it had 2 hits). I replaced them with median values in the category.