Skip to content
Snippets Groups Projects
Commit e8fdc33b authored by Malo Rosemeier's avatar Malo Rosemeier
Browse files

added docu, clean up

parent 780a250d
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
.gitignore
+ 1
− 2
View file @ e8fdc33b
......@@ -5,5 +5,4 @@
*.settings
*~
becas_*
*.buc
*.ldc
\ No newline at end of file
finstrip_*
\ No newline at end of file
README.md
+ 28
− 7
View file @ e8fdc33b
......@@ -6,9 +6,13 @@ FINSTRIPWrapper
FINSTRIPWrapper wraps WMC's [FINSTRIP](https://wmc.eu/bladknik.php) (FINite STRIP) buckling tool into Python and OpenMDAO enabling batch runs for optimization purposes.
# Dependencies and supported versions
* Python 2.7.x
* Numpy
* FINSTRIP as of FOCUS 6.3
* OpenMDAO 1.7.3
* FUSED-Wind dev
* PGL - Parametric Geometry Library
* BECAS (optional)
# OS
......@@ -53,10 +57,10 @@ sudo apt-get install wine1.8-amd64
winetricks
```
* install Sentinel HASP at least version 7.55, download from https://sentinelcustomer.gemalto.com/sentineldownloads/
* install Sentinel HASP (at least version 7.55), download from https://sentinelcustomer.gemalto.com/sentineldownloads/
* Sentinel HASP,LDK Script for Linux 32-bit Compatibility Packages
* Sentinel HASP LDK Ubuntu, DEB Run-time Installer
* add license server to `/etc/hasplm/hasplm.ini`
* add license server to `/etc/hasplm/hasplm.ini`, e. g.:
```bash
serveraddr = 10.16.51.84
```
......@@ -68,7 +72,7 @@ sudo service aksusbd start
sudo service aksusbd status
```
* check if license server is reachable via
* check if license dongle is reachable via
http://localhost:1947/_int_/devices.html
* add FINSTRIP executable to wine's PATH
......@@ -85,7 +89,7 @@ PATH = z:\\\home\\\mrosemeier\\\finstrip\\\
wine finstrip
```
* to ommit popup windows in batch install xvfb:
* to ommit popup windows in batch run install xvfb:
```bash
sudo apt-get install xvfb
```
......@@ -95,12 +99,29 @@ sudo apt-get install xvfb
xvfb-run -a wine finstrip test.buc
```
* clone FISNTRIPWrapper:
* clone PGL - Parametric Geometry Library:
```bash
$ cd git
$ git clone ssh://git@github.com:OpenMDAO/OpenMDAO.git
$ git checkout 1.7.3
$ git clone git@gitlab.windenergy.dtu.dk:frza/PGL.git
$ git checkout master
```
* clone FINSTRIPWrapper:
```bash
$ cd git
$ git clone git@gitlab.iwes.fraunhofer.de:bdt/FINSTRIPWrapper.git
$ git checkout master
```
* in case you want to calculate cross section properties for FINSTRIP with [BECAS](http://www.becas.dtu.dk/) install BECAS and clone BECASWrapper:
```bash
$ cd git
$ git clone git@gitlab.windenergy.dtu.dk:HAWTOpt2/BECASWrapper.git
$ git checkout master
```
## Framework
......
finstrip_wrapper/core/blade.py
+ 51
− 23
View file @ e8fdc33b
import numpy as np
from finstrip_wrapper.core.main import FINSTRIPWrapper, status
from finstrip_wrapper.core.main import FINSTRIPWrapper
from PGL.components.airfoil import AirfoilShape
class FINSTRIPWrapperBlade(FINSTRIPWrapper):
''' Base class for all blade models, i.e. BladeShell and BladeBeam
'''
A basic wrapper for FINSTRIP that provides three primary functionalities:
file input generation, execution and file output reading.
This class extends FINSTRIPWrapper doing blade specific setup
The following keys should be passed as dictionary, otherwise the standard
values are used.
keys
----
name: str
job name
elsize_factor: float
element length as fraction of arc length of airfoil surface
core_mats: list
list of material name strings to be identified as core
spanpos: float
span wise position of cross section
use_becas_stiffness: bool
use BECAS cross section properties if True
'''
def __init__(self, **kwargs):
......@@ -13,9 +32,14 @@ class FINSTRIPWrapperBlade(FINSTRIPWrapper):
# data base output name
self.name = 'blade'
self.core_mats = [] # list of materials to be identified as core
# element length as fraction of arc length of airfoil surface
self.elsize_factor = 0.015
# list of material name strings to be identified as core
self.core_mats = []
self.spanpos = 0 # position along span of analysed cross-section
# use BECAS cross section properties or not
self.use_becas_stiffness = False
for k, w in kwargs.iteritems():
......@@ -28,7 +52,7 @@ class FINSTRIPWrapperBlade(FINSTRIPWrapper):
self._init_finstrip_input_files_names()
def init_dicts(self, cs2d):
# dictionaries written by FINSTRIPBladeStructure according to FUSED-Wind format
# dictionaries written by FINSTRIPBladeStability according to FUSED-Wind format
# NOTE: absolut dimensions are required!
self.cs2d = cs2d
......@@ -79,7 +103,9 @@ class FINSTRIPWrapperBlade(FINSTRIPWrapper):
self._write_finstrip_analysis_files()
def _det_maxlength_element(self):
# calculate arc length of airfoil surface
''' calculate element length based on fraction of arc length of airfoil
surface
'''
af = AirfoilShape()
af.initialize(points=self.cs2d['coords'])
......@@ -138,7 +164,7 @@ class FINSTRIPWrapperBlade(FINSTRIPWrapper):
G13 = matprops[i, 7] * 1E-6 # to N/mm*2
G23 = matprops[i, 8] * 1E-6 # to N/mm*2
rho = matprops[i, 9] * 1E-9 # to kg/mm*3
# deterimne material type
# determine material type
if label in self.core_mats:
# mat_type = 'core'
_write_core(
......@@ -199,30 +225,32 @@ class FINSTRIPWrapperBlade(FINSTRIPWrapper):
I_xx = self.csprops_ref[7]
I_yy = self.csprops_ref[8]
I_xy = self.csprops_ref[9]
nu = self.csprops_ref[17] + np.pi / 2
# TODO: check if sign is correct in I_xy!
nu = self.csprops_ref[17]
# TODO: check if sign is correctly taken for I_xy!
# A = self.csprops_ref[15]
# Gere 12-30
theta = 0.5 * np.arctan(-2 * I_xy / (I_xx - I_yy))
theta = 0.5 * np.arctan(-2 * I_xy / (I_xx - I_yy)) + np.pi / 2
# Gere 12-32
R = np.sqrt((0.5 * (I_xx - I_yy))**2 + I_xy**2)
#R = np.sqrt((0.5 * (I_xx - I_yy))**2 + I_xy**2)
# Gere 12-31a
theta1 = 0.5 * np.arccos((I_xx - I_yy) / (2 * R))
#theta1 = 0.5 * np.arccos((I_xx - I_yy) / (2 * R))
# Gere 12-31b
theta2 = 0.5 * np.arcsin(-I_xy / R)
#theta2 = 0.5 * np.arcsin(-I_xy / R)
# if not np.sign(theta) == np.sign(nu):
# raise RuntimeError('BECAS input not consistent')
# check the output format of BECASWrapper
config = {}
config['BECASWrapper'] = {}
config['BECASWrapper']['hawc2_FPM'] = True
from becas_wrapper.becas_bladestructure import set_cs_size
cs_size, _ = set_cs_size(config)
if cs_size == len(self.csprops):
#hawc2_FPM = True
AE = self.csprops[20] * 1E-3 # to N/mm
......@@ -241,7 +269,7 @@ class FINSTRIPWrapperBlade(FINSTRIPWrapper):
(x_e, y_e, AE, EI_xx, EI_yy, EI_xy) + n)
fid = open(self.file_input, 'w')
for i, line in enumerate(out_str):
for line in out_str:
fid.write(line)
fid.close()
......@@ -258,12 +286,11 @@ class FINSTRIPWrapperBlade(FINSTRIPWrapper):
my = self.csloads[case, 4]
mz = self.csloads[case, 5]
n = '\n'
# write header
out_str.append(
'loads, %s, %i, %g, %g, %g, %g, %g' % (self.spanpos, fx, fy, fz, mx, my, mz) + n)
fid = open(self.file_load, 'w')
for i, line in enumerate(out_str):
for line in out_str:
fid.write(line)
fid.close()
......@@ -301,14 +328,15 @@ class FINSTRIPWrapperBlade(FINSTRIPWrapper):
# x-u p1, y-u p1, x-u p2, y-u p2
for case in range(self.ncase):
for mode in range(self.nmode):
self.mode_shape[case, mode, :, :] = np.loadtxt('R' +
'%05i' % (self.spanpos) +
'_LC' +
'%03d' % (case) +
'_mode'
'%03d' % (mode) +
self.txt_sfx,
skiprows=2)
# convert to m
self.mode_shape[case, mode, :, :] = 1E-3 * np.loadtxt('R' +
'%05i' % (self.spanpos) +
'_LC' +
'%03d' % (case) +
'_mode'
'%03d' % (mode) +
self.txt_sfx,
skiprows=2)
idxs = []
# identify surface and web nodes
......
finstrip_wrapper/core/main.py
+ 19
− 40
View file @ e8fdc33b
......@@ -7,44 +7,33 @@ import sys
class FINSTRIPWrapper(object):
"""
A basic wrapper of FINSTRIP that provides two primary functionalities:
model generation, meshing and solving
parameters
----------
finstrip_work_dir: Str
ANSYS working directory
job_name: str
standard: 'blade'
ANSYS job name
file_input: str
standard: 'runFINSTRIP.inp'
FINSTRIP input file name
number_cores: int
number of cores for ANSYS execution
workspace_memory: int
ANSYS workspace memory, units in MByte
database_memory: int
standards: 100
ANSYS database memory, units in MByte
"""
'''
A basic wrapper for FINSTRIP that provides three primary functionalities:
file input generation, execution and file output reading.
The following keys should be passed as dictionary, otherwise the standard
values are used.
keys
----
nmode: int
number of buckling modes to be processed
lambda_min: float
minimum half wavelength in m
lambda_max: float
maximum half wavelength in m
'''
def __init__(self, **kwargs):
super(FINSTRIPWrapper, self).__init__()
# element length as fraction of arc length of airfoil surface
self.elsize_factor = 0.015
# number of buckling modes to be processed
self.nmode = 21
# finstrip program parameters
self.lambda_min = 0.01 # minimum half wavelength in m
self.lambda_max = 15.126 # maximum half wavelength in m
self.nmode = 21
#self.ncoord = 79
#self.ncase = 1
# flags
self.debug_mode = True
self.dry_run = False
......@@ -61,21 +50,11 @@ class FINSTRIPWrapper(object):
# input file names
self.buc_sfx = '.buc'
self.lcd_sfx = '.ldc'
# self.fdb_prfx = 'db' # why no "." ??
self.txt_sfx = '.txt'
self.file_input = self.name + self.buc_sfx
self.file_load = self.name + self.lcd_sfx
#self.file_output = self.name + self.fcdb_prfx
#self.file_mat = 'MAT' + self.fd_sfx
#self.file_failmat = 'FAILMAT' + self.fd_sfx
#self.file_et = 'ET' + self.fd_sfx
#self.file_mesh = 'MESH' + self.fi_prfx
#self.file_variables = 'VARS' + self.fi_prfx
#self.file_vars_analysis = 'VARS_ANALYSIS' + self.fi_prfx
# output file names
self.file_rst_eigen = self.name + '_result_status' + self.txt_sfx
self.file_rst_loadcase = self.name + '_result_status' + self.txt_sfx
def compute(self):
''' Runs FINSTRIP - setup/execution and post processing
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment

Datenschutz und Impressum