LabRPS at 16:39, 17 January 2025

2025-01-17T16:39:35Z

← Older revision		Revision as of 09:39, 17 January 2025
Line 23:		Line 23:
	'''Macro_Deodatis_1996_Simulation_Method.RPSMacro'''		'''Macro_Deodatis_1996_Simulation_Method.RPSMacro'''

	{{~~MacroCode~~\|code=		{{Code\|code=

	# -- coding: utf-8 --		# -- coding: utf-8 --

LabRPS at 14:02, 8 January 2025

2025-01-08T14:02:16Z

← Older revision		Revision as of 07:02, 8 January 2025
Line 1:		Line 1:
	{{Macro		{{Macro
	\|Name=Macro Deodatis 1996 ~~Wind~~ Simulation ~~Points~~		\|Name=Macro Deodatis 1996 Simulation Method
	\|Description=This macro ~~create simulation point~~ according to G. Deodatis 1996.		\|Description=This macro allows to simulate random wind velocity according to G. Deodatis (1996) simulation method.
	\|Author=LabRPS		\|Author=LabRPS
	\|Version=1.0		\|Version=1.0
	\|Date=~~2024~~-04-15		\|Date=2025-01-18
	\|RPSVersion=All		\|RPSVersion=All
	\|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]]

LabRPS at 13:50, 8 January 2025

2025-01-08T13:50:12Z

← Older revision		Revision as of 06:50, 8 January 2025
Line 8:		Line 8:
	\|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 allows to simulate random wind velocity according to G. Deodatis (1996) simulation method presented in the paper "Simulation of Ergodic Multivariate Stochastic Processes" which outlines an approach for simulating multivariate stochastic processes that are ergodic. Ergodicity in this context refers to the property that time averages converge to ensemble averages for a given process. The model focuses on methods that can generate multivariate realizations of such processes while ensuring that the statistical properties of the simulated data reflect the underlying stochastic processes.		This macro allows to simulate random wind velocity according to G. Deodatis (1996) simulation method presented in the paper [https://doi.org/10.1061/(ASCE)0733-9399(1996)122:8(778) Simulation of Ergodic Multivariate Stochastic Processes] which outlines an approach for simulating multivariate stochastic processes that are ergodic. Ergodicity in this context refers to the property that time averages converge to ensemble averages for a given process. The model focuses on methods that can generate multivariate realizations of such processes while ensuring that the statistical properties of the simulated data reflect the underlying stochastic processes.

	[[File:Deodatis1996SimulationMethodData.png\|1024px\|center]]		[[File:Deodatis1996SimulationMethodData.png\|1024px\|center]]

LabRPS: Created page with "{{Macro |Name=Macro Deodatis 1996 Wind Simulation Points |Description=This macro create simulation point according to G. Deodatis 1996. |Author=LabRPS |Version=1.0 |Date=2024-04-15 |RPSVersion=All |SeeAlso=Macro Deodatis 1996 Mean Wind Profile }} This macro allows to simulate random wind velocity according to G. Deodatis (1996) simulation method presented in the paper "Simulation of Ergodic Multivariate Stochastic Processes" whic..."

2025-01-08T13:47:36Z

Created page with "{{Macro |Name=Macro Deodatis 1996 Wind Simulation Points |Description=This macro create simulation point according to G. Deodatis 1996. |Author=LabRPS |Version=1.0 |Date=2024-04-15 |RPSVersion=All |SeeAlso=Macro Deodatis 1996 Mean Wind Profile }} This macro allows to simulate random wind velocity according to G. Deodatis (1996) simulation method presented in the paper "Simulation of Ergodic Multivariate Stochastic Processes" whic..."

New page

{{Macro
|Name=Macro Deodatis 1996 Wind Simulation Points
|Description=This macro create simulation point according to G. Deodatis 1996.
|Author=LabRPS
|Version=1.0
|Date=2024-04-15
|RPSVersion=All
|SeeAlso=[[Macro_Deodatis_1996_Mean_Wind_Profile|Macro Deodatis 1996 Mean Wind Profile]]
}}
This macro allows to simulate random wind velocity according to G. Deodatis (1996) simulation method presented in the paper "Simulation of Ergodic Multivariate Stochastic Processes" which outlines an approach for simulating multivariate stochastic processes that are ergodic. Ergodicity in this context refers to the property that time averages converge to ensemble averages for a given process. The model focuses on methods that can generate multivariate realizations of such processes while ensuring that the statistical properties of the simulated data reflect the underlying stochastic processes.

[[File:Deodatis1996SimulationMethodData.png|1024px|center]]

{{Caption|Simulated wind velocity data}}

[[File:Deodatis1996SimulationMethodPlot.png|600px|1024px|center]]

==Script==

You can find the source code on the following Github repository:
[https://github.com/LabRPS/LabRPS-macros/blob/master/WindLab/Deodatis1996/Deodatis1996_SimulationMethod.RPSMacro Get the code here!]

'''Macro_Deodatis_1996_Simulation_Method.RPSMacro'''

{{MacroCode|code=

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

import LabRPS
from LabRPS import Vector as vec

import WindLab
import WindLabObjects

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

__Name__ = 'Deodatis1996SimulationMethod'
__Comment__ = 'Simulate random velocity according to Deodatis 1996'
__Author__ = 'Koffi Daniel'
__Version__ = '0.1.0'
__Date__ = '08/01/2025'
__License__ = ''
__Web__ = ''
__Wiki__ = 'https://wiki.labrps.com/Macro_Deodatis_1996_Simulation_Method'
__Icon__ = ''
__Help__ = ''
__Status__ = ''
__Requires__ = ''
__Contact__ = ''
__Communication__ = ''
__Files__ = ''

import LabRPS
import WindLab
import WindLabObjects
from LabRPS import Vector as vec
import time

def simulate():
# Plugin
installResuslt = WindLab.installPlugin("WindLabPlugin")
if not installResuslt:
LabRPS.Console.PrintError("The installation the WindLabPlugin has failed.\n")
return None

# Document
doc = LabRPS.newDocument()

# Simulation
sim = WindLabObjects.makeSimulation(doc, "Simulation")
if not sim:
LabRPS.Console.PrintError("The simulation does not exist.\n")
return None

# set simulation parameters
sim.NumberOfFrequency = 2048
sim.MaxFrequency = "0.64 Hz" # 4 rad/s;
sim.FrequencyIncrement = "0.00031 Hz" # 0.00195 rad/s;
sim.TimeIncrement = "0.785 s"
sim.NumberOfTimeIncrements = 9651

# Simulation points
loc = WindLabObjects.makeFeature("SimulationPoints", "Simulation", "General Distribution", "Location Distribution")
if not loc:
LabRPS.Console.PrintError("Error on creating the location distribution.\n")
return None

v1 = vec(0, 0, 35)
v2 = vec(0, 0, 40)
v3 = vec(0, 0, 140)
loc.Locations = [v1, v2, v3]

# Mean wind Profile
mean = WindLabObjects.makeFeature("MeanSpeed", "Simulation", "Logarithmic Law Profile", "Mean Wind Profile")
if not mean:
LabRPS.Console.PrintError("The creation of the mean wind profile was not successuful.\n")
return None

mean.TerrainRoughness = '0.001266 m'
mean.ShearVelocity = '1.76 m/s'

# Frequencies
frequency = WindLabObjects.makeFeature("Frequencies", "Simulation", "Double Index Frequency Discretization", "Frequency Distribution")
if not frequency:
LabRPS.Console.PrintError("Error on creating the frequency distribution.\n")
return None

# Spectrum
spectrum = WindLabObjects.makeFeature("Spectrum", "Simulation", "Kaimal Along Wind Spectrum", "Along Wind Spectrum")
if not spectrum:
LabRPS.Console.PrintError("Error on creating the spectrum model.\n")
return None

# Coherence
coherence = WindLabObjects.makeFeature("CoherenceFunction", "Simulation", "Davenport Coherence Function", "Coherence Function")
if not coherence:
LabRPS.Console.PrintError("The creation of the coherence was not successuful.\n")
return None
coherence.ExponentialDecayCz = 10

# Spectrum decomposition
spectrumD = WindLabObjects.makeFeature("SpectrumDecomposition", "Simulation", "Cholesky Decomposition", "Spectrum Decomposition Method")
if not spectrumD:
LabRPS.Console.PrintError("Error on creating the spectrum decomposition method.\n")
return None

# Random phase
randomness = WindLabObjects.makeFeature("RandomPhases", "Simulation", "Uniform Random Phases", "Randomness Provider")
if not randomness:
LabRPS.Console.PrintError("The creation of the randomness provider was not successuful.\n")
return None

# Simulation method
simMethod = WindLabObjects.makeFeature("SimulationMethod", "Simulation", "Deodatis 1996", "Simulation Method")
if not simMethod:
LabRPS.Console.PrintError("Error on creating the simulation method.\n")
return None

# Run simulation and output the first(0) sample
# store starting time
begin = time.time()
velocities = sim.simulate(0)
# store end time
end = time.time()
LabRPS.Console.PrintMessage(f"Total runtime of the simulaltion is {end - begin} seconds\n")

if LabRPS.GuiUp:
import WindLabGui
import GeneralToolsGui
WindLabGui.setActiveSimulation(sim)
GeneralToolsGui.GeneralToolsPyTool.showArray(sim.getSimulationData().numberOfTimeIncrements, sim.getSimulationData().numberOfSpatialPosition + 1, velocities, True)

simulate()

}}

{{clear}}

Macro Deodatis 1996 Simulation Method - Revision history

LabRPS at 16:39, 17 January 2025

LabRPS at 14:02, 8 January 2025

LabRPS at 13:50, 8 January 2025