Skip to content

work on example 8, where anba4 fails #41

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
ptrbortolotti merged 2 commits into from
Jan 8, 2026
Merged

work on example 8, where anba4 fails #41

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
b41f555
work on example 8, where anba4 fails
ptrbortolotti Jan 8, 2026
690d40f
Changing the custom mesh material number in the example since a mater…
JustinPorter88 Jan 8, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion SONATA/utl/converter_WT.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -270,7 +270,7 @@ def converter_WT(blade, cs_pos, byml, materials, mesh_resolution):

if 'fiber_orientation' in sec.keys():
set_interp = PchipInterpolator(sec['fiber_orientation']['grid'], sec['fiber_orientation']['values'])
tmp2[i]['segments'][0]['layup'][id_layer]['orientation'] = float(set_interp(x[i]) * 180 / np.pi) # added float
tmp2[i]['segments'][0]['layup'][id_layer]['orientation'] = float(set_interp(x[i])) # added float
else:
tmp2[i]['segments'][0]['layup'][id_layer]['orientation'] = 0.

Expand Down
81 changes: 33 additions & 48 deletions examples/8_fiber_orient/8_box_beam.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,5 +1,7 @@
name: Box_Beam
description: Simple rectangular beam
windIO_version: '2.0.1'
assembly:
number_of_blades: 1
components:
blade:
reference_axis:
Expand All @@ -15,15 +17,18 @@ components:
outer_shape:
chord:
grid: [0.0, 1.0]
values: [1., 1.]
values: [2., 2.]
twist:
grid: [0.0, 1.0]
values: [0.0, 0.0]
# Can use this to verify the SONATA applying twist option:
# values: [0.17453292519943295, 0.17453292519943295]
section_offset_y:
grid: [0.0, 1.0]
values: [0., 0.]
values: [1., 1.]
rthick:
grid: [0.0, 1.0]
values: [0.5, 0.5]
airfoils:
- name: rectangular
spanwise_position: 0.0
Expand All @@ -41,10 +46,10 @@ components:
material: CarbonUD
fiber_orientation:
grid: [0.0, 1.0]
values: [0.2618, 0.2618]
values: [15.0, 15.0]
thickness:
grid: [0.0, 1.0]
values: [0.01, 0.01]
values: [0.001, 0.001]
start_nd_arc:
anchor:
name: TE
Expand All @@ -57,10 +62,10 @@ components:
material: CarbonUD
fiber_orientation:
grid: [0.0, 1.0]
values: [0.2618, 0.2618]
values: [15.0, 15.0]
thickness:
grid: [0.0, 1.0]
values: [0.01, 0.01]
values: [0.001, 0.001]
start_nd_arc:
anchor:
name: TE
Expand All @@ -73,10 +78,10 @@ components:
material: CarbonUD
fiber_orientation:
grid: [0.0, 1.0]
values: [0.2618, 0.2618]
values: [15.0, 15.0]
thickness:
grid: [0.0, 1.0]
values: [0.01, 0.01]
values: [0.001, 0.001]
start_nd_arc:
anchor:
name: TE
Expand All @@ -89,10 +94,10 @@ components:
material: CarbonUD
fiber_orientation:
grid: [0.0, 1.0]
values: [0.2618, 0.2618]
values: [15.0, 15.0]
thickness:
grid: [0.0, 1.0]
values: [0.01, 0.01]
values: [0.001, 0.001]
start_nd_arc:
anchor:
name: TE
Expand All @@ -105,10 +110,10 @@ components:
material: CarbonUD
fiber_orientation:
grid: [0.0, 1.0]
values: [0.2618, 0.2618]
values: [15.0, 15.0]
thickness:
grid: [0.0, 1.0]
values: [0.01, 0.01]
values: [0.001, 0.001]
start_nd_arc:
anchor:
name: TE
Expand All @@ -121,10 +126,10 @@ components:
material: CarbonUD
fiber_orientation:
grid: [0.0, 1.0]
values: [0.2618, 0.2618]
values: [15.0, 15.0]
thickness:
grid: [0.0, 1.0]
values: [0.01, 0.01]
values: [0.001, 0.001]
start_nd_arc:
anchor:
name: TE
Expand All @@ -137,10 +142,10 @@ components:
material: CarbonUD
fiber_orientation:
grid: [0.0, 1.0]
values: [0.2618, 0.2618]
values: [15.0, 15.0]
thickness:
grid: [0.0, 1.0]
values: [0.01, 0.01]
values: [0.001, 0.001]
start_nd_arc:
anchor:
name: TE
Expand All @@ -153,10 +158,10 @@ components:
material: CarbonUD
fiber_orientation:
grid: [0.0, 1.0]
values: [0.2618, 0.2618]
values: [15.0, 15.0]
thickness:
grid: [0.0, 1.0]
values: [0.01, 0.01]
values: [0.001, 0.001]
start_nd_arc:
anchor:
name: TE
Expand All @@ -169,10 +174,10 @@ components:
material: CarbonUD
fiber_orientation:
grid: [0.0, 1.0]
values: [0.2618, 0.2618]
values: [15.0, 15.0]
thickness:
grid: [0.0, 1.0]
values: [0.01, 0.01]
values: [0.001, 0.001]
start_nd_arc:
anchor:
name: TE
Expand All @@ -185,10 +190,10 @@ components:
material: CarbonUD
fiber_orientation:
grid: [0.0, 1.0]
values: [0.2618, 0.2618]
values: [15.0, 15.0]
thickness:
grid: [0.0, 1.0]
values: [0.01, 0.01]
values: [0.001, 0.001]
start_nd_arc:
anchor:
name: TE
Expand All @@ -208,43 +213,24 @@ components:
airfoils:
- name: rectangular
coordinates: # follow s-coordinate along the outer profile of the structure
x: [1.0, 1.0, -1.0, -1.0, 1.0, 1.0] # in SI units - mm - box beam analogous to Popescu/Hodges 2000 paper p. 551/552 & Hodges 1992 paper
y: [0.0, 0.5, 0.5, -0.5, -0.5, 0.0]
x: [1.0, 1.0, 0.0, 0.0, 1.0, 1.0] # in SI units - mm - box beam analogous to Popescu/Hodges 2000 paper p. 551/552 & Hodges 1992 paper
y: [0.0, 0.25, 0.25, -0.25, -0.25, 0.0]
rthick: 0.5000
aerodynamic_center: 0.5000
polars: # dummy 2pi polars (not used in example)
- configuration: Default
re_sets:
- re: 3.00E+06
c_m:
cm:
grid: [-3.141592653589793, 3.141592653589793]
values: [0.0, 0.0]
c_l:
cl:
grid: [-3.141592653589793, 3.141592653589793]
values: [0.0, 0.0]
c_d:
cd:
grid: [-3.141592653589793, 3.141592653589793]
values: [0.5, 0.5]
materials:
- name: steel
description: Steel of the tower and monopile ASTM A572 Grade 50
source: http://www.matweb.com/search/DataSheet.aspx?MatGUID=9ced5dc901c54bd1aef19403d0385d7f
orth: 0
rho: 7800
alpha: 0.0
E: 200000000000.0
nu: 0.3
G: 79300000000.0
GIc: 0
GIIc: 0
alp0: 0
Xt: 450000000.0
Xc: 450000000.0
S: 0
Xy: 345000000.0
m: 3
A: 35534648443.719765
unit_cost: 0.7
- name: CarbonUD
E: [114500000000.0, 8390000000.0, 8390000000.0]
G: [5990000000.0, 5990000000.0, 5990000000.0]
Expand All @@ -265,4 +251,3 @@ materials:
waste: 0.05
fiber_density: 1800.
area_density_dry: 1.000
component_id: 4
2 changes: 1 addition & 1 deletion examples/8_fiber_orient/8_sonata_orientation.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -84,7 +84,7 @@
cells = np.array([[1, 2, 3, 0]])

# check job.materials
MatID = 2*np.ones(1)
MatID = np.ones(1)

job.blade_custom_mesh(nodes, cells, MatID, theta_3=15)

Expand Down
Loading