Cadquery for wind turbine blade

[kkpal3]

Hi, I'm trying to use cadquery to build a wind turbine blade. I'm defining the blade sections by sections along the spanwise direction. Each section has a specific airfoil shape, twist and chord length. I have successfully build the airfoils sections along the span by inputting coordinates and closing them by polyline() function.

After all the sections are built, I define a path along the spanwise direction to do the sweep() function. However, I always get errors like "Standard_ConstructionError: BRepFIll:: profiles are inconsistent", or "Standard_NoSuchObject:BRepFill::SameNumberByPolarMethod failed", even though I use the same number of coordinates for each airfoil section! I believe this problem is caused by the different airfoil shape along the span, as I can successfully build the blade using the same airfoil for all spanwise locations.

Does anyone have suggestions how to solve this problem? Many thanks!

[jmwright]

I think that @winksaville has done some work along these lines.

Potentially related discussion: https://groups.google.com/g/cadquery/c/_7JI1Ffl2N4/m/LbVVfSC3CAAJ

There might be some potentially useful info here, although loft is being used instead of sweep.

Is there a simple example you can post of what you're trying to accomplish?

[kkpal3]

Thank you Jeremy! I will have a look the references you mentioned.
Meanwhile, here is a sketch of what I want to do, and a test case is also attached.

testCase.zip

I have information of the blade stored in "BladeInfo.dat". It specifies the the following:
BlSpan: spanwise location of the airfoil section
BlTwist: local twist angle of the airfoil section
BlChord: local chord length of the airfoil
BlAFID: the index of airfoil type. The coordinates of different airfoils are stored in files named nacaxxxx.txt. They have the same number of points for different airfoils.

For now, I can build the airfoil sections in cad-query. I wanted to construct the blade from these sectional data. If I use the same type of airfoils, i.e., the last column of "BladeInfo.dat" (BlAFID) having the same value, the sweep operation is successful.
However, when I use different values of BlAFID, errors as I mentioned previously occurs!

[adam-urbanczyk]

I can run your code without any issue using the master version of CQ / CQ-editor. Which version are you using?

[marcus7070]

I can run your code without any issue using the master version of CQ / CQ-editor.

@adam-urbanczyk just to be sure, the example as @kkpal3 provided it does run, but when you modify the provided BladeInfo.dat and change one of the numbers in the last column to eg. 2, then it doesn't run. Is that the case for you?

When I change the last profile from 1 to 2, I get about 5 minutes of 100% cpu usage then Standard_ConstructionError: BRepFill :: profiles are inconsistent.

[kkpal3]

Thank you @adam-urbanczyk and @marcus7070 for the discussions, and sorry for not being clear about my questions. Yes, I meant when I use a different type of airfoil along the spanwise direction, for example, change the last profile from 1 to 2, it will give me errors.

I just checked my code again and found that I made a mistake!! I was not using the same number of points for the airfoil type 1, naca 6412. Once I discard this type and use the rest of the three airfoils for the test, cad-query seems to work well!!

[adam-urbanczyk]

Good to hear @kkpal3 . You might be also interested ins using spline instead of polyline and note that one profile is misisng or the spine should not start at (0,0,0). You can also just used loft and not specify any extrusion path at all.

[kkpal3]

Thanks @adam-urbanczyk for the suggestions.
I investigated my codes a little bit more, and found that the sweep operation still failed when there is a large change in the sectional geometry, for example, from a circle to a thin airfoil, using the same number of points for both geometries.
Anyway, I will try what you suggested and report more thorough results later.

[adam-urbanczyk]

Did you mange to get some interesting result @kkpal3 ? I'm really curious.

[kkpal3]

Sorry for the delay @adam-urbanczyk. I haven't got a satisfactory wind blade model yet. I will explain what I have tried so far. Let me know if you have any suggestions!
The code is here.
dtuTest.zip

Firstly, there is always problem when I try to loft over more than 2 blade sections whose shape changes drastically, say from circle to airfoil, see the fig below.

Secondly, even with similar airfoil shapes, there seems to be wiggles over a large number of sections, which is undesirable, see the fig below.

Thirdly, sweep seems to perform inferior to loft, so I did not try it much.

The best I can get so far is to loft segment-by-segment, as loft with two any TWO sections always work nicely. However, when using a lot number of sections, this approach also results in error (this error can be reproduced when using "BladeInfo.dat" instead of "BladeInfo2.dat", the latter is an abridged version of the former). The problem with this approach is that the surface of the blade is no longer smooth, as the slope at intersecting airfoil sections of two loft process is not ensured to match.

[parnoldx]

For wings a Gordon algorithm would be a better fit imo. Not sure if it can be made easily available but here they ported it to python from tigl c++ tpaviot/pythonocc-core#409

[adam-urbanczyk]

@kkpal3 I could not get this to work either. One way to slightly improve the situation is to use BRepOffsetAPI_ThruSections.SetMaxDegree from OCC (not exposed currently). Probably you'll need to investigate using Gordon Surface for this as mentioned above. It is not implemented in CQ (PR would be very welcome). You could use https://github.com/DLR-SC/tigl or https://github.com/tomate44/CurvesWB to check if it solves the issue.

[kkpal3]

Thanks for the comments. I just found a piece of code based on open cascade here. It seems to work nicely.
I'm very new to open cascade, but do you think it's possible to implement it in cadquery?