No changes between revisions
/Modules/Mechanical/WINDGAUGE01A/WINDGAUGE01A_small.png
Cannot display: file marked as a binary type.
svn:mime-type = application/octet-stream
Property changes:
Deleted: svn:mime-type
-application/octet-stream
\ No newline at end of property
/Modules/Mechanical/WINDGAUGE01A/SW/DataLogger.py
0,0 → 1,102
#!/usr/bin/python
 
# plot with >> plot 'last.txt' u 1:2 w l axes x1y1, 'last.txt' u 1:4 w l axes x1y2, 'last.txt' u 1:3
 
import os
import time
import datetime
import sys
import numpy as np
from gps import *
from pymlab import config
import threading
 
gpsd = None
 
class GpsPoller(threading.Thread):
def __init__(self):
threading.Thread.__init__(self)
global gpsd #bring it in scope
gpsd = gps(mode=WATCH_ENABLE)
self.current_value = None
self.running = True
def run(self):
global gpsd
while gpsp.running:
gpsd.next()
 
 
 
cfg = config.Config(
i2c = {
"port": 1,
},
bus = [
{
"name": "rps",
"type": "rps01",
},
],
)
 
 
cfg.initialize()
 
print "RPS01A logger"
sensor = cfg.get_device("rps")
 
 
try:
angles = np.zeros(5)
angles[4] = sensor.get_angle(verify = False)
time.sleep(0.01)
angles[3] = sensor.get_angle(verify = False)
time.sleep(0.01)
angles[2] = sensor.get_angle(verify = False)
time.sleep(0.01)
angles[1] = sensor.get_angle(verify = False)
n = 0
speed = 0
AVERAGING = 50
filen = 'log%0.0f.txt'%time.time()
f = open(filen,'w')
os.remove("last.txt")
os.symlink(filen, "last.txt")
gpsp = GpsPoller()
gpsp.start()
 
while True:
for i in range(AVERAGING):
time.sleep(0.01)
angles[0] = sensor.get_angle(verify = False)
if (angles[0] + n*360 - angles[1]) > 300:
n -= 1
angles[0] = angles[0] + n*360
 
elif (angles[0] + n*360 - angles[1]) < -300:
n += 1
angles[0] = angles[0] + n*360
 
else:
angles[0] = angles[0] + n*360
speed += (-angles[4] + 8*angles[3] - 8*angles[1] + angles[0])/12
angles = np.roll(angles, 1)
speed = speed/AVERAGING
 
g_spd = gpsd.fix.speed
print "W_Spd: %0.2f \t Angle: %0.2f \t G_Spd %0.2f" % (speed, angles[0], g_spd)
f.write("%0.2f %0.2f %0.2f %0.2f\r\n" %(time.time(), abs(speed), angles[0], g_spd))
f.flush()
 
except KeyboardInterrupt:
gpsp.running = False
gpsp.join()
sys.exit(0)
/Modules/Mechanical/WINDGAUGE01A/SW/wind_gauge.py
0,0 → 1,107
#!/usr/bin/python
 
# MLAB meteostation wind speed gauge with magnetic rotation sensor.
# This simple algorithm calculate difference between five time equidistant points during the rotation. The result is angular speed per time step.
# Size of time-step could be varied depending on expected wind speed range to measure.
# Algorithm should be expanded by Kalman filtering to minimize dependence on fast reading.
# The measuring principle could introduce time-stamped reading to increase precision of measurement. It could be possible because the readings are not exactly time equidistant in real Linux word.
 
#uncomment for debbug purposes
#import logging
#logging.basicConfig(level=logging.DEBUG)
 
import time
import datetime
import sys
import numpy as np
from pymlab import config
 
#### Script Arguments ###############################################
 
if len(sys.argv) != 2:
sys.stderr.write("Invalid number of arguments.\n")
sys.stderr.write("Usage: %s PORT ADDRESS\n" % (sys.argv[0], ))
sys.exit(1)
 
port = eval(sys.argv[1])
#### Sensor Configuration ###########################################
 
''''
cfg = config.Config(
i2c = {
"port": port,
},
 
bus = [
{
"type": "i2chub",
"address": 0x72,
"children": [
{"name": "encoder", "type": "rps01", "channel": 1, }
],
},
],
)
 
'''
cfg = config.Config(
i2c = {
"port": port,
},
bus = [
{
"name": "encoder",
"type": "rps01",
},
],
)
 
 
cfg.initialize()
 
print "RPS01A magnetic position sensor RPS01 readout example \r\n"
sensor = cfg.get_device("encoder")
 
print sensor.get_address()
print sensor.get_zero_position()
 
#### Data Logging ###################################################
 
try:
angles = np.zeros(5)
angles[4] = sensor.get_angle(verify = False)
time.sleep(0.01)
angles[3] = sensor.get_angle(verify = False)
time.sleep(0.01)
angles[2] = sensor.get_angle(verify = False)
time.sleep(0.01)
angles[1] = sensor.get_angle(verify = False)
n = 0
speed = 0
AVERAGING = 50
 
while True:
for i in range(AVERAGING):
time.sleep(0.01)
angles[0] = sensor.get_angle(verify = False)
if (angles[0] + n*360 - angles[1]) > 300:
n -= 1
angles[0] = angles[0] + n*360
 
elif (angles[0] + n*360 - angles[1]) < -300: # compute angular speed in backward direction.
n += 1
angles[0] = angles[0] + n*360
 
else:
angles[0] = angles[0] + n*360
speed += (-angles[4] + 8*angles[3] - 8*angles[1] + angles[0])/12
angles = np.roll(angles, 1)
 
speed = speed/AVERAGING # apply averaging on acummulated value.
print "Speed: %0.2f \t Total Angle: %0.2f \r\n" % (speed, angles[0])
 
except KeyboardInterrupt:
sys.exit(0)
/Modules/Mechanical/WINDGAUGE01A/SW/Data_analyser.ipynb
0,0 → 1,200
{
"metadata": {
"name": "",
"signature": "sha256:9453f2d297004717b7b0ab5bb80d8d05d3ff319f7e609fa9565e753aa36cef87"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": [
"import h5py\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 1
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"file = h5py.File('test_rps.hdf5', 'r') # 'r' means that hdf5 file is open in read-only mode\n",
"dataset = file['RPS01']"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 2
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"print dataset.value[4,2]"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"36.7822\n"
]
}
],
"prompt_number": 3
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"plt.plot( dataset.value[:,0], dataset.value[:,2], color = 'r')\n",
"#plt.xlim([1.45415117E9,1.4541513E9])\n",
"plt.show()"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 4
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"prev_val= dataset.value[0,2]\n",
"n = 0\n",
"angle = np.zeros((dataset.shape[0]))\n",
"for i in range(dataset.value.shape[0]):\n",
" if (dataset.value[i,2] - prev_val) > 300:\n",
" n -= 1\n",
" angle[i] = dataset.value[i,2] + n*360\n",
" prev_val = dataset.value[i,2]\n",
" elif -(dataset.value[i,2] - prev_val) > 300: # compute angular speed in backward direction.\n",
" n += 1\n",
" angle[i] = dataset.value[i,2] - n*360\n",
" prev_val = dataset.value[i,2]\n",
" else:\n",
" angle[i] = dataset.value[i,2] + n*360\n",
" prev_val = dataset.value[i,2]\n",
" "
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 6
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Five point difference numerical calculation. Source: http://mathfun528.blogspot.cz/2011/07/numerical-differentiation.html"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"angle_speed = np.zeros_like(angle)\n",
"\n",
"for i in range(2,angle.shape[0]-2):\n",
" angle_speed[i] = (-angle[i + 2] + 8*angle[i + 1] - 8*angle[i - 1] + angle[i - 2])/12"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 7
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"fig, ax1 = plt.subplots()\n",
"\n",
"ax2 = ax1.twinx()\n",
"ax1.set_xlabel('Sample #')\n",
"ax1.set_ylabel('Angle')\n",
"ax2.set_ylabel('Angular speed')\n",
"\n",
"ax1.plot(dataset.value[:,0], angle,'b',dataset.value[:,0], dataset.value[:,2],'r')\n",
"ax2.plot(dataset.value[:,0], angle_speed,'g')\n",
"\n",
"plt.show()"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 8
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"print angle"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"[ 3.64746094e+01 3.65405273e+01 3.77929688e+01 ..., -1.31529551e+05\n",
" -1.31530452e+05 -1.31528452e+05]\n"
]
}
],
"prompt_number": 10
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 4
},
{
"cell_type": "markdown",
"metadata": {},
"source": []
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 45
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 45
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"metadata": {},
"outputs": []
}
],
"metadata": {}
}
]
}
/Modules/Mechanical/WINDGAUGE01A/SW/test_rps.hdf5
Cannot display: file marked as a binary type.
svn:mime-type = application/octet-stream
Property changes:
Added: svn:mime-type
+application/octet-stream
\ No newline at end of property
/Modules/Mechanical/WINDGAUGE01A/SW/test_rps.log.hdf5
Cannot display: file marked as a binary type.
svn:mime-type = application/octet-stream
Property changes:
Added: svn:mime-type
+application/octet-stream
\ No newline at end of property
/Modules/Mechanical/WINDGAUGE01A/PrjInfo.txt
5,9 → 5,9
Mìøení smìru vìtru
 
[InfoLongDescription.en]
Printable variation sensor for measuring the wind direction. The detection is used three-axis magnetometer.
 
 
[InfoLongDescription.cs]
Tisknutelná varianta senzoru pro mìøení smìru vìtru. K detekci je využito tøíosého magnetometru.
 
[End]
/Modules/Mechanical/WINDGAUGE01A/CAD/SRC/WINDGAUGE01A_D01.scad
87,5 → 87,3
}
}
 
 
//WINDGAUGE01A_D01();
/Modules/Mechanical/WINDGAUGE01A/CAD/SRC/WINDGAUGE01A_D02.scad
15,16 → 15,9
union()
{
//výstuha
translate([-D01_sirka_tyce/2+S01_sila_materialu,0,0])
rotate ([0,-90,0])
linear_extrude (height = S01_sila_materialu, convexity = 10)
polygon(points=[[D02_vyska_uchytky/2,D01_delka_tyce+S01_prumer_vnitrni/2+D01_material_pred_zavitem+D02_sila_materialu-1],[0, D01_delka_tyce+S01_prumer_vnitrni/2+D01_material_pred_zavitem+D02_sila_materialu-1],[0,S01_prumer_vnitrni/2+D01_material_pred_zavitem],[D01_material_pod_zavitem,S01_prumer_vnitrni/2+D01_material_pred_zavitem]]);
//výstuha 2
translate([D01_sirka_tyce/2,0,0])
rotate ([0,-90,0])
linear_extrude (height = S01_sila_materialu, convexity = 10)
polygon(points=[[D02_vyska_uchytky/2,D01_delka_tyce+S01_prumer_vnitrni/2+D01_material_pred_zavitem+D02_sila_materialu-1],[0, D01_delka_tyce+S01_prumer_vnitrni/2+D01_material_pred_zavitem+D02_sila_materialu-1],[0,S01_prumer_vnitrni/2+D01_material_pred_zavitem],[D01_material_pod_zavitem,S01_prumer_vnitrni/2+D01_material_pred_zavitem]]);
//zavit
translate([0,0,D01_material_pod_zavitem])
56,7 → 49,7
//uchyt na sloupek
translate([0,D01_delka_tyce+(S01_prumer_vnitrni-S01_tolerance_zavit)/2+D01_material_pred_zavitem+D02_sila_materialu+(D02_prumer_obruby)/2,D02_vyska_uchytky/2])
rotate(a=[0,0,270])
WINDGAUGE01A_D01();
WINDGAUGE01A_D01_2();
 
}
 
/Modules/Mechanical/WINDGAUGE01A/CAD/assembly.scad
29,31 → 29,26
{
union()
{
//WINDGAUGE1A_D01 DRŽÁK
//WINDGAUGE02A_D01 DRŽÁK
//-------------------------------------------------------------
//-------------------------------------------------------------
 
 
color(barva_D01)
translate([0,90,-2*posunuti_dilu+D02_vyska_uchytky/2])
rotate(a=[0,0,90])
WINDGAUGE01A_D01_2();
translate([0,0,-2*posunuti_dilu])
WINDGAUGE01A_D01();
//WINDGAUGE01A_D02 DRŽÁK druhy dil
//WINDGAUGE02A_D02 DRŽÁK druhy dil
//-------------------------------------------------------------
//-------------------------------------------------------------
 
 
 
//WINDGAUGE02A_S01 stator velky dil
//-------------------------------------------------------------
color(barva_D01)
translate([0,0,-2*posunuti_dilu])
WINDGAUGE01A_D02();
//-------------------------------------------------------------
 
color(barva_S01)
61,8 → 56,8
rotate(a=[0,180,0])
WINDGAUGE01A_S01();
 
color(barva_S02)
WINDGAUGE01A_S02();
//color(barva_S02)
//WINDGAUGE01A_S02();
 
color(barva_S03)
translate([0,0,S01_vyska_spodniho_zavitu+S01_sila_drzaku_RJ11-0.3-4*posunuti_dilu])
80,7 → 75,7
rotate(a=[0,0,0])
WINDGAUGE01A_R02();
translate([-R02_sila_materialu_kridla/2,S01_prumer_vnitrni/2+10,S01_vyska+2*S01_sila_materialu+3*posunuti_dilu+6*posunuti_dilu+R04_zavit_vyska+7]) rotate(a=[90,0,90])
translate([-R02_sila_materialu_kridla/2,0,S01_vyska+2*S01_sila_materialu+3*posunuti_dilu+6*posunuti_dilu+R04_zavit_vyska+10]) rotate(a=[0,90,0])
WINDGAUGE01A_R05();
 
color(barva_R03)
/Modules/Mechanical/WINDGAUGE01A/CAD/configuration.scad
75,7 → 75,7
D01_sirka_tyce=10;
 
//WINDGAUGE02A_D02 //uchyt na sloup
D02_sila_materialu=5;
D02_sila_materialu=4;
D02_vyska_uchytky=40;
D02_prumer_obruby=27; // 3/4 trubka
D02_delka_celeho_uchytu=70;
/Modules/Mechanical/WINDGAUGE01A/CAD/print_data/P_WINDGAUGE01A_D01.scad
9,7 → 9,7
PI=3.141592;
rotate([0, 0, 0])
translate([0,0,D02_vyska_uchytky/2])
WINDGAUGE01A_D01_2();
WINDGAUGE01A_D01();
 
 
/Modules/Mechanical/WINDGAUGE01A/TODO.txt
1,0 → 0,0

prejmenovat slozku configuration na src.
 
Kryt musi mit presah nejspis az pod lopatky.
 
Mam pochybnosti o presnosti tisku a pevnosti stredove osy. Myslim, ze by zde mel byt nerezovy sroub o spravnem prurezu primo do loziska.
 
Asi by se mela prodlouzit delka lopatek na maximalni tisknutelny rozmer, protoze vetsi polomer lopatek bude vice tolerantni k promenlivemu odporu statoru vuci rotoru, napriklad namrzanim, navlhnutim atd.. 
 
Musí se pocitat se zamrzanim, mozna bude nutne umistit vyhrivani.