Macro Deodatis 1996 Wind Simulation Points: Difference between revisions

From LabRPS Documentation
Jump to navigation Jump to search
No edit summary
No edit summary
 
(4 intermediate revisions by the same user not shown)
Line 1: Line 1:
{{Macro
{{Macro
|Name=Macro Deodatis 1996 Wind Simulation Points
|Name=Macro Deodatis 1996 Wind Simulation Points
|Icon=Macro_Deodatis_1996_Wind_Simulation_Points.png
|Description=This macro create simulation point according to G. Deodatis 1996.
|Description=This macro create simulation point according to G. Deodatis 1996.
|Author=Koffi Daniel
|Author=Koffi Daniel
Line 7: Line 6:
|Date=2024-04-15
|Date=2024-04-15
|RPSVersion=All
|RPSVersion=All
|Download=[https://wiki.freecad.org/images/f/f5/Macro_3D_Parametric_Curve.png ToolBar Icon]
|SeeAlso=[[Macro_Deodatis_1996_Mean_Wind_Profile|Macro Deodatis 1996 Mean Wind Profile]]
|SeeAlso=[[Macro_Deodatis_1996_Mean_Wind_Profile|Macro Deodatis 1996 Mean Wind Profile]]
}}
}}
This macro create simulation point according to G. Deodatis 1996. There are three simulation points created with WindLabFeatureLocationDistribution feature.
This macro create simulation point according to G. Deodatis 1996. There are three simulation points created with WindLabFeatureLocationDistribution feature.
[[File:Deodatis1996WindSimulationPoints.png|600px]]
[[File:Deodatis1996WindSimulationPoints.png|1024px|center]]
 
{{Caption|Simulation Point Distribution}}
{{Caption|Simulation Point Distribution}}


Line 21: Line 20:
[https://github.com/LabRPS/LabRPS-macros/blob/master/WindLab/Deodatis1996/Deodatis1996_LocationDistribution.RPSMacro Get the code here!]
[https://github.com/LabRPS/LabRPS-macros/blob/master/WindLab/Deodatis1996/Deodatis1996_LocationDistribution.RPSMacro Get the code here!]


'''Macro Deodatis 1996 Wind Simulation Points.FCMacro'''
'''Macro Deodatis 1996 Wind Simulation Points.RPSMacro'''


{{MacroCode|code=
{{MacroCode|code=
Line 34: Line 33:
import WindLabObjects
import WindLabObjects


#install the WindLab plugin with all its plugged features
class Deodatis1996LocationDistribution():
installResuslt = WindLab.installPlugin("WindLabPlugin")
    def __init__(self):
        #install the WindLab plugin with all its plugged features
        installResuslt = WindLab.installPlugin("WindLabPlugin")


doc = LabRPS.newDocument()
        # create new document
        doc = LabRPS.newDocument()


# create new simulation with default name "Simulation"
        # create new simulation with default name "Simulation" in the above document
sim = WindLabObjects.makeSimulation(doc)
        sim = WindLabObjects.makeSimulation(doc)
      
      
# check if the simulation is successfully created
        # check if the simulation is successfully created
if not sim:
        if not sim:
    LabRPS.Console.PrintError("The simulation does not exist.\n")
            LabRPS.Console.PrintError("The simulation does not exist.\n")
            return None
 
        # create new location distribution in the simulation called "Simulation"
        loc = WindLabObjects.makeFeature("SimulationPoints", "Simulation", "Location Distribution",  "General Distribution")
      
      
# set simulation parameters
        # check if the created location distribution feature is null
sim.NumberOfFrequency = 2048
        if not loc:
sim.MaxFrequency = "0.64 Hz" # 4 rad/s;
            LabRPS.Console.PrintError("Error on creating the location distribution.\n")
sim.FrequencyIncrement = "0.00031 Hz" # 0.00195 rad/s;
            return None
sim.TimeIncrement = "0.785 s"
sim.NumberOfTimeIncrements = 9651


# create new location distribution in the simulation called "Simulation"
        # set the coordinates of the simulation points
loc = WindLabObjects.makeFeature("SimulationPoints", "Simulation", "Location Distribution", "General Distribution")
        v1 = vec(0, 0, 35)
   
        v2 = vec(0, 0, 40)
# check if the created location distribution feature is null
        v3 = vec(0, 0, 140)
if not loc:
        loc.Locations = [v1, v2, v3]
    LabRPS.Console.PrintError("Error on creating the location distribution.\n")


v1 = vec(0, 0, 35)
        # run the method "computeLocationCoordinateMatrixP3" of the feature and store the result.
v2 = vec(0, 0, 40)
        locRes = sim.computeLocationCoordinateMatrixP3()
v3 = vec(0, 0, 140)
loc.Locations = [v1, v2, v3]


# run the method "computeLocationCoordinateMatrixP3" of the feature and store the result.
        if LabRPS.GuiUp:
# Note that here, this method will be called for the active location distribution feature
            import WindLabGui
locRes = sim.computeLocationCoordinateMatrixP3()
            #show the location coordinate in Alphaplot
            WindLabGui.WindLabPyTool.showArray(sim.getSimulationData().numberOfSpatialPosition, 4, locRes)
            loc.ViewObject.Document.activeView().viewAxonometric()
            loc.ViewObject.Document.activeView().fitAll()


if LabRPS.GuiUp:
            import numpy
    import WindLabGui
            import matplotlib.pyplot as plt
    #show the location coordinate in Alphaplot
            arr = numpy.asarray(locRes)
    WindLabGui.WindLabPythonPluginUtilities.showArray(sim.getSimulationData().numberOfSpatialPosition, 4, locRes)
            fig, ax = plt.subplots()
    loc.ViewObject.Document.activeView().viewAxonometric()
            fig.suptitle('Location Distribution')
    loc.ViewObject.Document.activeView().fitAll()
            ax.plot(arr[:,0], arr[:,3], 'o-')
            ax.set_xlabel('Simulation Point Number')
            ax.set_ylabel('Height (m)')
            ax.grid(True)
            plt.show()
   
        doc.recompute()


    import numpy
Deodatis1996LocationDistribution()
    import matplotlib.pyplot as plt
    arr = numpy.asarray(locRes)
    fig, ax = plt.subplots()
    fig.suptitle('Location Distribution')
    ax.plot(arr[:,0], arr[:,3], 'o-')
    ax.set_xlabel('Simulation Point Number')
    ax.set_ylabel('Height (m)')
    ax.grid(True)
    plt.show()
 
doc.recompute()


}}
}}


{{clear}}
{{clear}}

Latest revision as of 13:44, 16 April 2024

Macro Deodatis 1996 Wind Simulation Points.png Macro Deodatis 1996 Wind Simulation Points

Description
This macro create simulation point according to G. Deodatis 1996.

Macro version: 1.0
Last modified: 2024-04-15
LabRPS version: All
Author: Koffi Daniel
Author
Koffi Daniel
Download
None
Links
Macro Version
1.0
Date last modified
2024-04-15
LabRPS Version(s)
All
Default shortcut
None
See also
Macro Deodatis 1996 Mean Wind Profile

This macro create simulation point according to G. Deodatis 1996. There are three simulation points created with WindLabFeatureLocationDistribution feature.

Deodatis1996WindSimulationPoints.png

Simulation Point Distribution

Points Heights‎

Script

You can find the source code on the following Github repository: Get the code here!

Macro Deodatis 1996 Wind Simulation Points.RPSMacro

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

import LabRPS
from LabRPS import Vector as vec

import WindLab
import WindLabObjects

class Deodatis1996LocationDistribution():
    def __init__(self):
        #install the WindLab plugin with all its plugged features
        installResuslt = WindLab.installPlugin("WindLabPlugin")

        # create new document
        doc = LabRPS.newDocument()

        # create new simulation with default name "Simulation" in the above document
        sim = WindLabObjects.makeSimulation(doc)
    
        # check if the simulation is successfully created
        if not sim:
            LabRPS.Console.PrintError("The simulation does not exist.\n")
            return None

        # create new location distribution in the simulation called "Simulation"
        loc = WindLabObjects.makeFeature("SimulationPoints", "Simulation", "Location Distribution",  "General Distribution")
    
        # check if the created location distribution feature is null
        if not loc:
            LabRPS.Console.PrintError("Error on creating the location distribution.\n")
            return None

        # set the coordinates of the simulation points
        v1 = vec(0, 0, 35)
        v2 = vec(0, 0, 40)
        v3 = vec(0, 0, 140)
        loc.Locations = [v1, v2, v3]

        # run the method "computeLocationCoordinateMatrixP3" of the feature and store the result.
        locRes = sim.computeLocationCoordinateMatrixP3()

        if LabRPS.GuiUp:
            import WindLabGui
            #show the location coordinate in Alphaplot
            WindLabGui.WindLabPyTool.showArray(sim.getSimulationData().numberOfSpatialPosition, 4, locRes)
            loc.ViewObject.Document.activeView().viewAxonometric()
            loc.ViewObject.Document.activeView().fitAll()

            import numpy
            import matplotlib.pyplot as plt
            arr = numpy.asarray(locRes)
            fig, ax = plt.subplots()
            fig.suptitle('Location Distribution')
            ax.plot(arr[:,0], arr[:,3], 'o-')
            ax.set_xlabel('Simulation Point Number')
            ax.set_ylabel('Height (m)')
            ax.grid(True)
            plt.show()
    
        doc.recompute()

Deodatis1996LocationDistribution()