wawi 0.0.8__py3-none-any.whl → 0.0.11__py3-none-any.whl

wawi/model/_screening.py

@@ -0,0 +1,124 @@
+import numpy as np
+import json
+from ._hydro import Seastate
+from pathlib import Path
+
+'''
+SCREENING SUBMODULE
+'''
+
+class ScreeningCase:
+    def __init__(self, seastate, parvar, independent=True, name=None):
+        
+        self.name = name
+        self.seastate = seastate
+        self.assign_parvar(parvar)
+        self.independent = independent
+        self.combos = self.get_parameter_space()
+        self.ix = -1
+        
+        if not self.independent:
+            sz_prev = None
+            for key in self.parvar:
+                sz = len(self.parvar[key])
+                
+                if sz_prev is not None and sz!=sz_prev:
+                    raise ValueError('If dependent parameter arrays are requested, they must have the same length!')
+                    
+                sz_prev = sz*1
+        
+        
+    def assign_parvar(self, parvar):
+        self.parvar = dict()
+        
+        for key in parvar:
+            if type(parvar[key]) is str:
+                self.parvar[key] = eval(parvar[key])
+            else:
+                self.parvar[key] = np.array(parvar[key])
+                
+        # Convert angles
+        conversions = {'theta0': np.pi/180.0, 'thetaU': np.pi/180.0}
+        
+        for key in self.parvar:
+            if key in conversions:
+                self.parvar[key] = self.parvar[key]*conversions[key]
+                
+            
+    def get_parameter_space(self):
+        pars = [self.parvar[k] for k in self.parvar]
+        keys = [k for k in self.parvar if k]
+        
+        if self.independent:
+            combos = np.array(np.meshgrid(*pars)).reshape(len(keys),-1).T
+        else:
+            combos = np.vstack(pars).T
+            
+        combo_dicts = [dict(zip(keys, combo)) for combo in combos]
+        return combo_dicts
+
+    @property
+    def n(self):
+        if self.independent:
+            return np.prod([len(v) for v in self.parvar.values()])
+        else:
+            return len(list(self.parvar.values())[0])
+
+    # Alternative constructor
+    @classmethod
+    def from_json(cls, json_file, **kwargs):
+        with open(json_file, 'r') as fileobj:
+            data = json.load(fileobj)
+            
+        seastate = Seastate.from_json(data['seastate'], **kwargs)           
+    
+        # Update options if provided (to enable overriding options from screening setup)
+        if 'options' in data:
+            options = data['options']
+        else:
+            options = {}
+            
+        if 'pontoon_options' in data:
+            pontoon_options = data['pontoon_options']
+        else:
+            pontoon_options = {}
+            
+        seastate.options.update(**options)
+        seastate.pontoon_options.update(**pontoon_options)
+        
+        parvar = data['parvar']
+        if 'independent' in data:
+            independent = data['independent']
+        else:
+            independent = True
+            
+        return cls(seastate, parvar, independent=independent, name=Path(json_file).stem)
+    
+    def get_combo(self):
+        return self.combos[self.ix]
+    
+    def get_next_combo(self):
+        self.iterate_ix()
+        combo = self.combos[self.ix]
+        return combo
+            
+    def iterate_seastate(self):
+        combo = self.get_next_combo()
+        if combo is not None:
+            for key in combo:
+                setattr(self.seastate, key, combo[key])
+            
+        return self.seastate
+    
+    def get_next_seastate(self):
+        self.iterate_seastate()            
+        return self.seastate
+      
+    def iterate_ix(self):
+        if self.ix == (self.n-1):
+            self.ix = 0     #reset
+        else:
+           self.ix += 1
+        
+    def reset_ix(self):
+        self.ix = 0

wawi/time_domain.py

@@ -52,6 +52,51 @@ def simulate_mdof(S, omega, fs=None, tmax=None, reg_factor=None, zero_limit=1e-1
                   phase_angles=None, return_phase_angles=False, interpolation_kind='linear', 
                   component_scaling=None, print_status=False):
     
+    '''
+    Simulate time series from given cross-spectral density matrix.
+
+    Parameters
+    ----------------
+    S : float
+        cross-spectral density matrix (Ndofs x Ndofs x Nfreqs) as complex numpy array
+    omega : float
+        numpy array defining frequencies in rad/s
+    fs : float, optional
+        sampling frequency in Hz; if not given, the value is defined as two times the maximum of `omega`
+    tmax : float, optional
+        duration of realization in seconds
+    reg_factor : float, optional
+        to help the Cholesky decomposition to achieve a numerical solution, 
+        a diagonal matrix with the norm of all frequency components of the matrix S scaled by the given factor
+        - if no value is given (--> None imposed), no regularization is conducted; 
+        used as input to function `spectrum_to_process` which decomposes the spectral density 
+    zero_limit : float, optional
+        frequency components where the norm of S is below this limit is set to zero
+        if no value is given (--> None imposed), no zero limit is applied; 
+        used as input to function `spectrum_to_process` which decomposes the spectral density 
+    phase_angles : float, optional
+        if previously conducted simulations are to be recreated, they can by inputting the phase angles
+        (to get the phase angles from simulation, see `return_phase_angles` input parameter)
+    return_phase_angles : bool, default=False
+        whether or not to return phase angle enabling recreation of the simulation
+    interpolation_kind : optional, default='linear'
+        interpolation type used for interpolation (to ensure that we get defined duration and sampling) prior to ifft
+    component_scaling : float, optional
+        values to use for scaling of components prior to decomposition (and thereafter rescale afterwards);
+        can help if some components are much smaller than others; size should match number of DOFs in S
+    print_status : False, optional
+        whether or not to print status messages
+
+    Returns
+    ----------------
+    p : float
+        time history 
+    t : float
+        numpy array with time axis values corresponding to `p`
+    alpha : float
+        phase angles, only returned if `return_phase_angles` is True
+    '''
+    
     if omega[0] !=0:
         omega = np.insert(omega, 0, 0)
         S = np.dstack([0*S[:,:,0], S])

wawi/wind.py

@@ -6,7 +6,7 @@ from scipy.special import jv as besselj, yv as bessely
 from .general import rodrot, blkdiag
 from .plot import plot_ads
 
-conv_text=''' 
+conv_text = r''' 
 -----------------------------------------------------
 |                                                   |
 | ~ ~ ~~~ ~  ~~ ~ /^^^^^^^^^^^^\ 88ooo...  .  .  .  |

{wawi-0.0.8.dist-info → wawi-0.0.11.dist-info}/METADATA

@@ -1,6 +1,6 @@
-Metadata-Version: 2.2
+Metadata-Version: 2.4
 Name: wawi
-Version: 0.0.8
+Version: 0.0.11
 Summary: WAve and WInd response prediction
 Author-email: "Knut A. Kvåle" <knut.a.kvale@ntnu.no>, Ole Øiseth <ole.oiseth@ntnu.no>, Aksel Fenerci <aksel.fenerci@ntnu.no>, Øivind Wiig Petersen <oyvind.w.petersen@ntnu.no>
 License: MIT License
@@ -43,6 +43,7 @@ Requires-Dist: trame
 Requires-Dist: ipywidgets
 Requires-Dist: pyvistaqt
 Requires-Dist: beefpy
+Dynamic: license-file
 
 ![WAWI logo](https://raw.githubusercontent.com/knutankv/wawi/main/wawi-logo-animated.svg)
 =======================
@@ -133,6 +134,13 @@ Examples
 =======================
 Examples are provided as Jupyter Notebooks in the [examples folder](https://github.com/knutankv/wawi/tree/main/examples).
 
+The examples are structured in the following folders based on their topic:
+
+* **0 Model generation and setup** - *Describing how a model is defined, either using input files together with the function `wawi.io.import_folder` or using the classes available in the `wawi.model` module directly (in the latter case, the open source Python library BEEF is used to create the required parameters directly in the notebook).*
+* **1 Modal analysis** - *Describing how to set up iterative (and incremental) modal analyses to represent the contributions from aerodynamics and hydrodynamices. Also, an example showing how to set up a multi-modal flutter analysis is given.*
+* **2 Response prediction** - *Describing how to conduct response analyses using WAWI. This includes assigning wind states, sea states, and the necessary commands to run a frequency domain analysis. Furthermore, more advanced topics such as wave-current interaction, inhomogeneous waves and stochastic linearization to represent quadratic drag damping are given in separate examples. Three models are considered: (i) a simple curved floating bridge, (ii) a single beam, (iii) a suspension bridge.*
+* **3 Software interfacing** - *Describing how to export necessary data from other software (limited to Abaqus for now) to construct a WAWI model.*
+
 References
 =======================
 

wawi-0.0.11.dist-info/RECORD

@@ -0,0 +1,38 @@
+examples/3 Software interfacing/Abaqus model export/bergsoysund-export.py,sha256=-1jrDNmevOcB1IV0iLOiW0z-QiNkOv8GkGk-7K8CcK0,2606
+examples/3 Software interfacing/Abaqus model export/test.py,sha256=kAmol3XZQ2jkO2qujEQrqx38Q24LbzrYZzqE5ng5-E8,2188
+tests/test_IABSE_step11a.py,sha256=guRsaP9gpERvSp0Ui3hO-WH764BcIe6Qx8vC4QVTAoU,8013
+tests/test_IABSE_step11c.py,sha256=U_NcL6BAmnQG7MSjkr8KYk9hH8ZPdZgg17pM9ugUBzo,10412
+tests/test_IABSE_step2a.py,sha256=L5AyTOJZ_sPmcmW_4h-sQfRwvVW71LZkXc1yz33DABM,10203
+tests/test_wind.py,sha256=r4rx6f8Tn-u4fn4gGPBMrL_JJQ2SVHGQiQ9sQuMQCPo,1707
+wawi/__init__.py,sha256=rH8H6W3AN1iZb8rU9K_8k-mycRSvQNxFq97Biy0-hzw,183
+wawi/fe.py,sha256=22QKI1GlfsG7o_TpFXaKJfzmbO2_2zdIMaoJmaIZdmY,4001
+wawi/general.py,sha256=xHRoDkcchrL6Y7pTUqGsjBHh8YBLvX9dYcNXXCjQap8,13787
+wawi/identification.py,sha256=bVB6EVRR6J39OO9ckuzNJ6f0FwIo4cLqfYgrsIN89TE,1748
+wawi/io.py,sha256=cr8XLF7bwpTCxpZ84YjOyWjc5uuu-JjGv8vY9UKH-HI,25327
+wawi/modal.py,sha256=WjNGFsk0C3tYRy18Q9WNRCatmGJtq1JSv0WrkGV02Eo,20247
+wawi/plot.py,sha256=jllJcjZxTBqjzBoT4k9jLXVUnie8oqNr8371IJvCd3c,19791
+wawi/prob.py,sha256=0nCdKdwkNf4M6sHyCZuYlt06gD0NmqRNfl4KesgySWA,215
+wawi/random.py,sha256=MHPpyTlRJSJFkCmeTAmw4Q5K1BPoFVb0Nxg0jDhkuIM,871
+wawi/signal.py,sha256=9HJs7VUhXOccuYPo12A0IUVoBIAJ2e_9F3rL-q3JuP4,1179
+wawi/structural.py,sha256=t25ohH4uBbzUJ7Hqn_kUfYhxcikZkRp8da-9dn7aEbw,8341
+wawi/time_domain.py,sha256=q8-H2xceP-2BmtOfbRQqYhD1JSb0z7jGq-dL_MzAX40,6122
+wawi/tools.py,sha256=-hFBvf0qK4AMn2MQRhrOitDMMMKm2QuRkVfbPBefEkQ,332
+wawi/wave.py,sha256=hhKg3KhKMBhOoCI7g8PFOGUbYgVhDyGmnYdBEdJ8mkY,16064
+wawi/wind.py,sha256=VHy07IbrCJxlqR0Z5O5WRc9YE4jb-MqFbULPiXhj0NA,38211
+wawi/wind_code.py,sha256=8OKLPpqc84obNNKBoYb2NunKjcn6a3i_pAWpIFEwg4Q,223
+wawi/ext/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+wawi/ext/abq.py,sha256=Wo-Zzb6b3lr_z7cNKTW-NdYJuVeAwZgwzlGhF_5Jym0,8905
+wawi/ext/ansys.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+wawi/ext/orcaflex.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+wawi/ext/sofistik.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+wawi/model/__init__.py,sha256=u6B-dugP76LBkOUVMs0KAoQ81PeRHwFL0M8MbNeAaJA,218
+wawi/model/_aero.py,sha256=tvWMjMoVi9ZVBc3NZ_s-wSofbEk1cc5jf5s-Swv2RxQ,9337
+wawi/model/_dry.py,sha256=KpmFlSinoY6DrSyc3Y0M8w1-cCC7VCFep-uzcqZsHz4,3940
+wawi/model/_hydro.py,sha256=1tUj19OgvEjItNtBxXuPHf0zWtKv_ioTrOcgl9Ov3dg,26960
+wawi/model/_model.py,sha256=34qblCqGpIKwbsa38CJ3a8KbHYyX8LDWOFQhdRwK98Q,52504
+wawi/model/_screening.py,sha256=NRYkKq928z2lqMSUTpbQLls04td_9R_4dhkjU3Gv1oQ,3716
+wawi-0.0.11.dist-info/licenses/LICENSE,sha256=bH1aWhrNbbPLrYnVFRaoYYzcUr-figHjry-kGB7Tc54,1076
+wawi-0.0.11.dist-info/METADATA,sha256=TPaFfAkztMCRaw4F8U5JT6XqOK2IjY3orqmBmguLoaY,6530
+wawi-0.0.11.dist-info/WHEEL,sha256=CmyFI0kx5cdEMTLiONQRbGQwjIoR1aIYB7eCAQ4KPJ0,91
+wawi-0.0.11.dist-info/top_level.txt,sha256=Nk5G_ZwgZRCb9ZMWZdr1M3QIskX6kCnlqeMl67N3zg8,20
+wawi-0.0.11.dist-info/RECORD,,

{wawi-0.0.8.dist-info → wawi-0.0.11.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (75.8.0)
+Generator: setuptools (78.1.0)
 Root-Is-Purelib: true
 Tag: py3-none-any
 

wawi-0.0.11.dist-info/top_level.txt

@@ -0,0 +1,3 @@
+examples
+tests
+wawi

wawi-0.0.8.dist-info/RECORD

@@ -1,21 +0,0 @@
-wawi/__init__.py,sha256=oZOxto9PG8WxUfEE02tcFY0o5-yLRR8lK8WXH30zxDE,158
-wawi/fe.py,sha256=22QKI1GlfsG7o_TpFXaKJfzmbO2_2zdIMaoJmaIZdmY,4001
-wawi/general.py,sha256=xHRoDkcchrL6Y7pTUqGsjBHh8YBLvX9dYcNXXCjQap8,13787
-wawi/identification.py,sha256=bVB6EVRR6J39OO9ckuzNJ6f0FwIo4cLqfYgrsIN89TE,1748
-wawi/io.py,sha256=NrziwlebaZUEdu8nWKNXKExHNZZKMzJ0gSHciBTJ4As,25323
-wawi/modal.py,sha256=WjNGFsk0C3tYRy18Q9WNRCatmGJtq1JSv0WrkGV02Eo,20247
-wawi/plot.py,sha256=jllJcjZxTBqjzBoT4k9jLXVUnie8oqNr8371IJvCd3c,19791
-wawi/prob.py,sha256=0nCdKdwkNf4M6sHyCZuYlt06gD0NmqRNfl4KesgySWA,215
-wawi/random.py,sha256=MHPpyTlRJSJFkCmeTAmw4Q5K1BPoFVb0Nxg0jDhkuIM,871
-wawi/signal.py,sha256=9HJs7VUhXOccuYPo12A0IUVoBIAJ2e_9F3rL-q3JuP4,1179
-wawi/structural.py,sha256=t25ohH4uBbzUJ7Hqn_kUfYhxcikZkRp8da-9dn7aEbw,8341
-wawi/time_domain.py,sha256=Oe-jviwDtBgpSmA7ZVmKEqQ5tdvsekXwOakYO1qUsN4,3841
-wawi/tools.py,sha256=-hFBvf0qK4AMn2MQRhrOitDMMMKm2QuRkVfbPBefEkQ,332
-wawi/wave.py,sha256=hhKg3KhKMBhOoCI7g8PFOGUbYgVhDyGmnYdBEdJ8mkY,16064
-wawi/wind.py,sha256=1rUqiEdMaUrhprs87TBdX24l1uDCOINaYkBChMLYWso,38208
-wawi/wind_code.py,sha256=8OKLPpqc84obNNKBoYb2NunKjcn6a3i_pAWpIFEwg4Q,223
-wawi-0.0.8.dist-info/LICENSE,sha256=bH1aWhrNbbPLrYnVFRaoYYzcUr-figHjry-kGB7Tc54,1076
-wawi-0.0.8.dist-info/METADATA,sha256=nlNip0I5cZlcT3lgLf5lXTysEeQ2y8nSVMapDTvoXlA,5197
-wawi-0.0.8.dist-info/WHEEL,sha256=In9FTNxeP60KnTkGw7wk6mJPYd_dQSjEZmXdBdMCI-8,91
-wawi-0.0.8.dist-info/top_level.txt,sha256=sE2NH_xVXnDKTGeIYzeX5IyU_j5vYDLl32v5uLiOMT4,5
-wawi-0.0.8.dist-info/RECORD,,

wawi-0.0.8.dist-info/top_level.txt

@@ -1 +0,0 @@
-wawi