Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in / Register
Toggle navigation
Menu
The Gitlab server is succesfully updated to version 15.1.0
Open sidebar
Frederik Zahle
SEAMRotor
Commits
3ddd6d68
Commit
3ddd6d68
authored
Oct 02, 2015
by
Frederik Zahle
Browse files
fixed error in norm
parent
d99cb462
Changes
2
Hide whitespace changes
Inline
Side-by-side
src/SEAMRotor/SEAMRotor.py
View file @
3ddd6d68
...
...
@@ -135,7 +135,7 @@ class SEAMBladeStructure(Component):
# Calculating thickness flapwise fatigue
tfat_flap
=
np
.
zeros
(
self
.
Nsections
)
for
i
in
range
(
self
.
Nsections
):
norm
=
self
.
solve_tfat_flap
(
0.01
,
C
[
i
],
thick
[
i
],
M
ex
t_flap
[
i
],
1.
)
norm
=
self
.
solve_tfat_flap
(
0.01
,
C
[
i
],
thick
[
i
],
M
fa
t_flap
[
i
],
1.
)
res
=
minimize
(
self
.
solve_tfat_flap
,
0.01
,
args
=
(
C
[
i
],
thick
[
i
],
Mfat_flap
[
i
],
norm
),
bounds
=
[(
1.e-6
,
0.5
)],
method
=
'SLSQP'
,
tol
=
1.e-8
)
tfat_flap
[
i
]
=
res
[
'x'
]
if
not
res
[
'success'
]:
print
'WARNING solve_tfat_flap'
,
i
,
res
...
...
@@ -143,7 +143,7 @@ class SEAMBladeStructure(Component):
# Calculating thickness edgewise fatigue
tfat_edge
=
np
.
zeros
(
self
.
Nsections
)
for
i
in
range
(
self
.
Nsections
):
norm
=
self
.
solve_tfat_edge
(
0.01
,
C
[
i
],
thick
[
i
],
M
ex
t_flap
[
i
],
1.
)
norm
=
self
.
solve_tfat_edge
(
0.01
,
C
[
i
],
thick
[
i
],
M
fa
t_flap
[
i
],
1.
)
res
=
minimize
(
self
.
solve_tfat_edge
,
0.01
,
args
=
(
C
[
i
],
thick
[
i
],
Mfat_edge
[
i
],
norm
),
bounds
=
[(
1.e-6
,
0.5
)],
method
=
'SLSQP'
,
tol
=
1.e-8
)
tfat_edge
[
i
]
=
res
[
'x'
]
if
not
res
[
'success'
]:
print
'WARNING solve_tfat_edge'
,
i
,
res
...
...
src/SEAMRotor/test/test_SEAMRotor.py
View file @
3ddd6d68
...
...
@@ -8,20 +8,20 @@ class SEAMRotorTestCase(unittest.TestCase):
def
setUp
(
self
):
pass
def
tearDown
(
self
):
pass
# add some tests here...
#def test_SEAMRotor(self):
#pass
if
__name__
==
"__main__"
:
# unittest.main()
top
=
SEAMBladeStructure
()
top
.
Nsections
=
21
.
top
.
Nsections
=
21
top
.
Neq
=
1e7
top
.
WohlerExpFlap
=
10.0
top
.
PMtarget
=
1.0
...
...
@@ -94,6 +94,3 @@ if __name__ == "__main__":
plt
.
grid
()
plt
.
show
()
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment