WindLab Tutorial GeorgeDeodatis 1996 Example SimulationPoints: Difference between revisions

Latest revision as of 19:20, 17 January 2025

Topic
Tutorial
Wind Simulation Points (George Deodatis, 1996)
Level
Beginner
Time to complete
10 minutes
Authors
Koffi Daniel
LabRPS version
0.1.0 or above
Example files
None
See also
None
Contents 1 Introduction 2 Requirements 3 Plugin Installation 4 Create New Feature 5 Input Points Coordinates 6 Showing Results 7 Scripting 8 What next?

Introduction

Deodatis (1996) provides a comprehensive method for simulating ergodic wind velocity time histories by modeling the wind as a stationary random process with well-defined statistical properties. His approach is particularly focused on multivariate wind simulations (i.e., for multiple correlated wind components at different locations in space. These locations where the wind velocity components are considered in space are called Simulation Point in WindLab. This example is meant to show how a wind simulation points in LabRPS's WindLab Workbench looks like in the LabRPS interface and how their coordinates can be visualized. Every computation in LabRPS needs a RPS Feature. Here we rather need a WindLab Feature (RPS Feature) which should be provided by a plugin.

Requirements

A compatible version of LabRPS designated in the tutorial overview.
Use the Help → About LabRPS to see the version of LabRPS installed.
No external software is needed for the computation of the locations coordinates as well as for visualizing the results.
Install the plugin that will provide the feature for the computation of the locations coordinates.

Plugin Installation

The first step is to start LabRPS by double-cliking its icon. Then active the WindLab workbench as shown in the following picture. Note that plugins are loaded according to the active workbench. If you do not activate the WindLab workbench, no wind velocity related plugin will be loaded.

According to Deodatis (1996), the simulation points are three in total and their distribution in space does not follow any particular distribution in space. To compute such simulation distribution, the General Distribution feature from the WindLab Plugin can be used. Let's install the feature first. Go to Tools → RPS Features..., select the WindLabPlugin in the list and click Install as shown in the following picture.

The Feature installation wizard will be launched. Please install the feautre by following the steps as shown in the following pictures:

Create New Feature

Please follow the following steps to create new WindLab simulation and add new Location distribution feature to it.

Press the New button to create new document or
- Select the File → New option from the menu or
- Use the keyboard shortcut: Ctrl+N.
Press the New Simulation button to create new WindLab simulation or
- Select the WindLab → New Simulation option from the menu.
- This will display the WindLab simulation creation task dialog as follows:

Press the Create Feature button to create new Location Distribution feature or
- Select the WindLab → Create Feature option from the menu.

Input Points Coordinates

We have created a locations distribution (simulation points) feature of type General Distribution. Now, we have to use the created feature to compute the coordinates of the simulation points. In the tree view, click on the simulationPoints feature, all the properties of this feature will be shown in the property editor. In the property editor scroll down and find a property named Locations. Click the ... near that property value to show a the vectors dialog that will allow you to input the coordinates of the three points: Point1(0,0,35), Point2(0,0,40) and Point3(0,0,140). Before inputting the points coordinates, you should first click the Table button to expand the points list. Note that value input in the LabRPS has default unit depending on the unit system you are using. Here, we are using the default unit system which is the Standard(mm/kg/s/degree). However, the coordinates are kept in metre because the vector dialog does not convert the values. That is why we keep the coordinates in metre as points coordinates outputs are expected in metre in LabRPS.

Showing Results

Now that everything is ready, you can display the simulation points in table form. Right click on the feature in the tree view and choose Compute Location Array in the context menu. The points will be shown as follows.

Scripting

you can manipulate the created object easily through scripting as follows:

# -*- coding: utf-8 -*-
# (c) Koffi Daniel 2024

import WindLab
import GeneralToolsGui
import WindLabObjects
from LabRPS import Vector as vec
import LabRPS
import numpy
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

def compute():
    installResuslt = WindLab.installPlugin("WindLabPlugin")

    doc = LabRPS.ActiveDocument
    if not doc:
       doc = LabRPS.newDocument()

    # create WindLab simulation called "Simulation"
    sim = WindLabObjects.makeSimulation(doc, "Simulation")
    
    # check if the simulation does really exist
    if not sim:
       LabRPS.Console.PrintError("The simulation does not exist.\n")
       return None 
  
    featureType = "General Distribution"
    featureGroup = "Location Distribution"
  
    # create the feature and add it to the existing simulation (you may refer to the WindLab Workbench page in 
    # case you don't understand the next line)
    genSimPoints= WindLabObjects.makeFeature("MyNewFeature", sim.Name, featureType, featureGroup)
    
    # check if the created feature is good
    if not genSimPoints:
       LabRPS.Console.PrintError("Error on creating the uniform points feature.\n")
       return None

    # create the simulation points by their coordinates
    v1 = vec(0, 0, 35)
    v2 = vec(0, 0, 40)
    v3 = vec(0, 0, 140)

    # add the points to the locations
    genSimPoints.Locations = [v1, v2, v3]

    # compute the simulation points coordinates. WindLab will internally use the "genSimPoints" feature
    simPoints = sim.computeLocationCoordinateMatrixP3()

    # now you can convert the coordinate matrix to numpy array and use it for any other purposes
    arr = numpy.asarray(simPoints)

    # Example 3D points
    x = arr[:,1]
    y = arr[:,2]
    z = arr[:,3]

    # you can also show the result in a table, pass False as last argument to the function to ask 
    # LabRPS to only show the data without plotting them
    GeneralToolsGui.GeneralToolsPyTool.showArray(sim.getSimulationData().numberOfSpatialPosition, 4, simPoints, False)

    # Create a figure
    fig = plt.figure()

    # Add 3D axes
    ax = fig.add_subplot(111, projection='3d')

    # Plot points
    ax.scatter(x, y, z, color='blue')

    # Hide all axes and labels
    ax.set_axis_off()

    # Set the title
    ax.set_title('3D Plotting of Points')

    # Show the plot
    plt.show()

compute()

you can find the entire scripting for this tutorial here.

What next?

We are now finished with the basic workflow for the WindLab Workbench feature creation.
You are now prepared to do the second WindLab tutorial.
We will create mean wind profile that will use the simulation points coordinates created in this tutorial.