diff options
| author | Marius Peter <blendoit@gmail.com> | 2019-06-20 18:02:31 -0700 |
|---|---|---|
| committer | Marius Peter <blendoit@gmail.com> | 2019-06-20 18:02:31 -0700 |
| commit | 4ee52cf24db8148c54f02022517a06b662509f92 (patch) | |
| tree | 9caf8b934cff255090466951c4e438c97bb7df77 /creator.py | |
| parent | d80486b906e68380afdcf3a3693eaa3e81b38a9c (diff) | |
elliptical and rectangular lift distributions
Diffstat (limited to 'creator.py')
| -rw-r--r-- | creator.py | 20 |
1 files changed, 10 insertions, 10 deletions
@@ -154,12 +154,12 @@ class Airfoil(Coordinates): """ y_c = float() if 0 <= x < p_c: - y_c = (m / (p**2)) * (2 * p * (x / self.chord) - - (x / self.chord)**2) + y_c = (m / (p ** 2)) * (2 * p * (x / self.chord) + - (x / self.chord) ** 2) elif p_c <= x <= self.chord: - y_c = (m / ((1 - p)**2)) * ((1 - 2 * p) - + 2 * p * (x / self.chord) - - (x / self.chord)**2) + y_c = (m / ((1 - p) ** 2)) * ((1 - 2 * p) + + 2 * p * (x / self.chord) + - (x / self.chord) ** 2) return (y_c * self.chord) def get_thickness(x): @@ -169,17 +169,17 @@ class Airfoil(Coordinates): y_t = 5 * t * self.chord * ( +0.2969 * sqrt(x / self.chord) - 0.1260 * (x / self.chord) - - 0.3516 * (x / self.chord)**2 - + 0.2843 * (x / self.chord)**3 - - 0.1015 * (x / self.chord)**4) + - 0.3516 * (x / self.chord) ** 2 + + 0.2843 * (x / self.chord) ** 3 + - 0.1015 * (x / self.chord) ** 4) return y_t def get_theta(x): dy_c = float() if 0 <= x < p_c: - dy_c = ((2 * m) / p**2) * (p - x / self.chord) + dy_c = ((2 * m) / p ** 2) * (p - x / self.chord) elif p_c <= x <= self.chord: - dy_c = (2 * m) / ((1 - p)**2) * (p - x / self.chord) + dy_c = (2 * m) / ((1 - p) ** 2) * (p - x / self.chord) theta = atan(dy_c) return theta |