Macro Deodatis 1996 Wind Simulation Points - Revision history

LabRPS: /* Script */

2025-01-18T02:59:13Z

Script

← Older revision		Revision as of 19:59, 17 January 2025
Line 23:		Line 23:

	{{Code\|code=		{{Code\|code=

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

	import WindLab		import WindLab
	import GeneralToolsGui		import GeneralToolsGui

LabRPS at 16:15, 17 January 2025

2025-01-17T16:15:43Z

← Older revision		Revision as of 09:15, 17 January 2025
Line 22:		Line 22:
	'''Macro_Deodatis_1996_Wind_Simulation_Points.RPSMacro'''		'''Macro_Deodatis_1996_Wind_Simulation_Points.RPSMacro'''

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

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

	~~import LabRPS~~
	~~from LabRPS import Vector as vec~~

	import WindLab		import WindLab
			import GeneralToolsGui
	import WindLabObjects		import WindLabObjects
			from LabRPS import Vector as vec
			import LabRPS
			import numpy

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

	~~# create new document~~		doc = LabRPS.ActiveDocument
	doc = LabRPS.newDocument()		if not doc:
			doc = LabRPS.newDocument()

	# create ~~new~~ simulation ~~with default name~~ "Simulation" ~~in the above document~~		# create WindLab simulation called "Simulation"
	sim = WindLabObjects.makeSimulation(doc)		sim = WindLabObjects.makeSimulation(doc, "Simulation")

			# check if the simulation does really exist
			if not sim:
			LabRPS.Console.PrintError("The simulation does not exist.\n")
			# abord the computation

			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 ~~simulation~~ is ~~successfully created~~		# check if the created feature is good
	if not ~~sim~~:		if not genSimPoints:
	LabRPS.Console.PrintError("~~The simulation does not exist~~.\n")		LabRPS.Console.PrintError("Error on creating the uniform points feature.\n")
	~~return None~~		# abord the computation

	# create ~~new location distribution in~~ the simulation ~~called "Simulation"~~		# create the simulation points by their coordinates
	~~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:~~
	~~LabRPS.Console.PrintError~~(~~"Error on creating the location distribution.\n"~~)
	~~return None~~

	# ~~set~~ the ~~coordinates of~~ the ~~simulation points~~		# add the points to the locations
	~~v1 = vec(0, 0, 35000)~~		genSimPoints.Locations = [v1, v2, v3]
	~~v2 = vec(0, 0, 40000)~~
	~~v3 = vec(0, 0, 140000)~~
	~~loc~~.Locations = [v1, v2, v3]

	# ~~run~~ the ~~method~~ "~~computeLocationCoordinateMatrixP3~~" ~~of the~~ feature ~~and store the result.~~		# compute the simulation points coordinates. WindLab will internally use the "genSimPoints" feature
	~~locRes~~ = sim.computeLocationCoordinateMatrixP3()		simPoints = sim.computeLocationCoordinateMatrixP3()

	~~if LabRPS.GuiUp:~~		# now you can convert the coordinate matrix to numpy array and use it for any other purposes
	~~import WindLabGui~~		arr = numpy.asarray(simPoints)
	#~~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~~		# you can also show the result in a table, pass False as last argument to the function to ask
	~~import matplotlib.pyplot as plt~~		# LabRPS to only show the data without plotting them
	~~arr = numpy~~.~~asarray(locRes)~~		GeneralToolsGui.GeneralToolsPyTool.showArray(sim.getSimulationData().numberOfSpatialPosition, 4, simPoints, False)
	~~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~~()		compute()

	}}		}}

	{{clear}}		{{clear}}

LabRPS at 21:58, 15 November 2024

2024-11-15T21:58:54Z

← Older revision		Revision as of 14:58, 15 November 2024
Line 20:		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~~.RPSMacro'''		'''Macro_Deodatis_1996_Wind_Simulation_Points.RPSMacro'''

	{{MacroCode\|code=		{{MacroCode\|code=

LabRPS at 21:52, 15 November 2024

2024-11-15T21:52:30Z

← Older revision		Revision as of 14:52, 15 November 2024
Line 2:		Line 2:
	\|Name=Macro Deodatis 1996 Wind Simulation Points		\|Name=Macro Deodatis 1996 Wind Simulation Points
	\|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=LabRPS
	\|Version=1.0		\|Version=1.0
	\|Date=2024-04-15		\|Date=2024-04-15
Line 58:		Line 58:

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

LabRPS at 21:44, 16 April 2024

2024-04-16T21:44:56Z

← Older revision		Revision as of 14:44, 16 April 2024
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

LabRPS at 20:56, 16 April 2024

2024-04-16T20:56:44Z

← Older revision		Revision as of 13:56, 16 April 2024
Line 34:		Line 34:
	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"		# create new location distribution in the simulation called "Simulation"
	loc = WindLabObjects.makeFeature("SimulationPoints", "Simulation", "Location Distribution", "General Distribution")		loc = WindLabObjects.makeFeature("SimulationPoints", "Simulation", "Location Distribution", "General Distribution")

	# check if the created location distribution feature is null		# check if the created location distribution feature is null
	if not loc:		if not loc:
	LabRPS.Console.PrintError("Error on creating the location distribution.\n")		LabRPS.Console.PrintError("Error on creating the location distribution.\n")
			return None

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

	# run the method "computeLocationCoordinateMatrixP3" of the feature and store the result.		# run the method "computeLocationCoordinateMatrixP3" of the feature and store the result.
	~~# Note that here, this method will be called for the active location distribution feature in case they are more than one~~		locRes = sim.computeLocationCoordinateMatrixP3()
	locRes = sim.computeLocationCoordinateMatrixP3()

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

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

	doc.recompute()		doc.recompute()

			Deodatis1996LocationDistribution()

	}}		}}

	{{clear}}		{{clear}}

LabRPS at 20:44, 16 April 2024

2024-04-16T20:44:58Z

← Older revision		Revision as of 13:44, 16 April 2024
Line 7:		Line 7:
	\|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]]
	}}		}}
Line 22:		Line 21:
	[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=

LabRPS at 19:52, 16 April 2024

2024-04-16T19:52:55Z

← Older revision		Revision as of 12:52, 16 April 2024
Line 11:		Line 11:
	}}		}}
	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}}

LabRPS at 18:18, 16 April 2024

2024-04-16T18:18:33Z

← Older revision		Revision as of 11:18, 16 April 2024
Line 45:		Line 45:
	if not sim:		if not sim:
	LabRPS.Console.PrintError("The simulation does not exist.\n")		LabRPS.Console.PrintError("The simulation does not exist.\n")

	~~# 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~~

	# create new location distribution in the simulation called "Simulation"		# create new location distribution in the simulation called "Simulation"
Line 66:		Line 59:

	# run the method "computeLocationCoordinateMatrixP3" of the feature and store the result.		# run the method "computeLocationCoordinateMatrixP3" of the feature and store the result.
	# Note that here, this method will be called for the active location distribution feature		# Note that here, this method will be called for the active location distribution feature in case they are more than one
	locRes = sim.computeLocationCoordinateMatrixP3()		locRes = sim.computeLocationCoordinateMatrixP3()

LabRPS at 18:14, 16 April 2024

2024-04-16T18:14:44Z

← Older revision		Revision as of 11:14, 16 April 2024
Line 10:		Line 10:
	\|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\|600px]]