ES2874575T3 - Procedure and system for determining the aerodynamic resistance of a cyclist - Google Patents
Procedure and system for determining the aerodynamic resistance of a cyclist Download PDFInfo
- Publication number
- ES2874575T3 ES2874575T3 ES16741578T ES16741578T ES2874575T3 ES 2874575 T3 ES2874575 T3 ES 2874575T3 ES 16741578 T ES16741578 T ES 16741578T ES 16741578 T ES16741578 T ES 16741578T ES 2874575 T3 ES2874575 T3 ES 2874575T3
- Authority
- ES
- Spain
- Prior art keywords
- person
- tfin
- tinic
- instants
- vector
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/08—Aerodynamic models
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C3/00—Measuring distances in line of sight; Optical rangefinders
- G01C3/02—Details
- G01C3/06—Use of electric means to obtain final indication
- G01C3/08—Use of electric radiation detectors
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/06—Measuring arrangements specially adapted for aerodynamic testing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Fluid Mechanics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Electromagnetism (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Force Measurement Appropriate To Specific Purposes (AREA)
Abstract
Procedimiento de determinación de la resistencia aerodinámica de un objeto o de una persona que tiene un movimiento propio, caracterizado por que comprende por lo menos las etapas siguientes: - adquirir entre dos instantes tinic y tfin unas nubes de puntos 3D que corresponden a las posiciones en el espacio de diferentes puntos posicionados en diferentes lugares del objeto o de la persona, - determinar unos modelos 3D en función del tiempo, para cada instante entre tinic y tfin, a partir de los puntos de las nubes de puntos adquiridos en cada instante entre los instantes tinic y tfin, - simular numéricamente el desplazamiento de los sólidos descritos por las nubes de puntos que corresponden al objeto o a la persona en un fluido, para cada instante entre los instantes tinic y tfin, estando dicha simulación en particular parametrizada por la dirección y la velocidad de desplazamiento del fluido con respecto al sólido representadas por un vector D, - calcular la media de las fuerzas de resistencia aerodinámica ejercidas sobre el objeto o la persona entre dichos instantes tinic y tfin, a partir de dicha simulación numérica del desplazamiento de los sólidos.Procedure for determining the aerodynamic resistance of an object or a person that has its own movement, characterized in that it comprises at least the following stages: - acquire between two instants tinic and tfin some 3D point clouds corresponding to the positions in the space of different points positioned at different places on the object or person, - determine 3D models as a function of time, for each instant between tinic and tfin, from the points of the point clouds acquired at each instant between instants tinic and tfin, - numerically simulate the displacement of the solids described by the point clouds that correspond to the object or the person in a fluid, for each instant between the instants tinic and tfin, said simulation being in particular parameterized by the direction and the speed of displacement of the fluid with respect to the solid represented by a vector D, - calculate the average of the forces of res aerodynamic resistance exerted on the object or person between said instants tinic and tfin, from said numerical simulation of the displacement of the solids.
Description
DESCRIPCIÓNDESCRIPTION
Procedimiento y sistema de determinación de la resistencia aerodinámica de un ciclistaProcedure and system for determining the aerodynamic resistance of a cyclist
Campo técnicoTechnical field
La presente invención se refiere a un procedimiento y a un sistema de determinación de la resistencia aerodinámica de un objeto y/o de una persona en movimiento. Más particularmente, la invención se refiere a un procedimiento y a un sistema de determinación de la resistencia aerodinámica de un ciclista. Le permitirá mejorar sus rendimientos adaptando su posición sobre la bicicleta para minimizar dicha resistencia aerodinámica.The present invention relates to a method and a system for determining the aerodynamic resistance of an object and / or a person in motion. More particularly, the invention relates to a method and a system for determining the aerodynamic resistance of a cyclist. It will allow you to improve your performance by adapting your position on the bike to minimize said aerodynamic drag.
Estado de la técnicaState of the art
Es muy conocido que la resistencia aerodinámica está compuesta por dos fuerzas, es decir la fuerza de resistencia al avance y la fuerza de sustentación. En materia de ciclismo, por ejemplo, a la vista de las velocidades de desplazamiento, la fuerza de sustentación se considera insignificante de manera que se asimila por lo tanto la resistencia aerodinámica a la resistencia al avance.It is well known that aerodynamic drag is made up of two forces, that is, the drag force and the lift force. In cycling, for example, in view of travel speeds, lift force is considered negligible so that aerodynamic resistance is therefore assimilated to drag.
De manera bien conocida, la fuerza de resistencia al avance viene dada por la ecuación siguiente:In a well-known manner, the drag force is given by the following equation:
en la que Ap designa el área frontal proyectada del ciclista y de su bicicleta, Cd el coeficiente de resistencia al avance, p la densidad del aire y Vf la velocidad de flujo del fluido sobre el cuerpo del ciclista.in which Ap designates the projected frontal area of the cyclist and his bicycle, Cd the coefficient of resistance to advance, p the density of the air and Vf the speed of fluid flow over the body of the cyclist.
En el marco de un estudio aerodinámico, la cantidad que se busca minimizar es el área frontal efectiva dada por la ecuación siguiente:In the framework of an aerodynamic study, the quantity to be minimized is the effective frontal area given by the following equation:
á / t R* á / t R *
!> .’ f.Yf !> . ' f. and f
Con el fin de determinar el área frontal efectiva, es necesario conocer por lo tanto la fuerza de resistencia al avance Ra. Se observará que el área frontal efectiva no es una constante, sino que evoluciona ligeramente con la velocidad de desplazamiento.In order to determine the effective frontal area, it is therefore necessary to know the drag force Ra. It will be seen that the effective frontal area is not a constant, but changes slightly with the speed of travel.
Con el fin de evaluar el área frontal efectiva de un ciclista, se conocen bien varios procedimientos, estando estos diferentes procedimientos descritos en la publicación de Vincent Chabroux "Approche aérodynamique et biomécanique de l'amélioration des performances de cyclistes en course contre la montre" PhD thesis, Université de la Méditerranée, 2010 y en la publicación de F. Grappe. "Cyclisme et optimisation de la performance: sciences et méthodologie de l'entraínement." Sciences et pratiques du sport, De Boeck, 2005.In order to evaluate the effective frontal area of a cyclist, several procedures are well known, these different procedures being described in Vincent Chabroux's publication "Approche aerodynamique et biomécanique de l'amélioration des performances de cyclistes en course contre la montre" PhD thesis, Université de la Méditerranée, 2010 and in the publication of F. Grappe. "Cyclisme et optimization de la performance: sciences et méthodologie de l'entraínement." Sciences et pratiques du sport, De Boeck, 2005.
Esencialmente, estos procedimientos se pueden clasificar en cuatro grupos: el análisis de potencia mecánica, la medición por remolcado, el procedimiento de desaceleración y la medición en túnel aerodinámico.Essentially, these procedures can be classified into four groups: the mechanical power analysis, the towed measurement, the deceleration procedure, and the wind tunnel measurement.
El análisis de la potencia mecánica consiste en medir la potencia suministrada por el ciclista gracias a un potenciómetro a diferentes velocidades y después en evaluar la fuerza de resistencia al avance por medio de una regresión lineal a partir de dichas mediciones de la potencia suministrada. Este procedimiento es relativamente preciso pero necesita varios ensayos en condiciones reales.The mechanical power analysis consists of measuring the power supplied by the cyclist thanks to a potentiometer at different speeds and then evaluating the drag force by means of a linear regression from these measurements of the supplied power. This procedure is relatively accurate but requires several tests under real-world conditions.
La medición por remolcado consiste en medir las fuerzas de tracción remolcando a un ciclista sobre su bicicleta con un automóvil a velocidad constante. Este procedimiento es bastante limitativo y carece de precisión debido a las turbulencias provocadas por el vehículo remolcador.Towing measurement consists of measuring traction forces by towing a cyclist on his bicycle with a car at constant speed. This procedure is quite restrictive and lacks precision due to turbulence caused by the tow vehicle.
El procedimiento por desaceleración consiste en medir la desaceleración del ciclista a diferentes velocidades. De esta manera, utilizando las leyes de Newton, se calcula la fuerza de resistencia al avance. Este procedimiento es muy pesado de realizar y necesita que el ciclista esté inmóvil.The deceleration procedure consists of measuring the cyclist's deceleration at different speeds. In this way, using Newton's laws, the drag force is calculated. This procedure is very heavy to perform and requires the rider to be immobile.
Las mediciones en túnel aerodinámico consisten en generar un flujo de aire sobre el conjunto ciclista/bicicleta y en cuantificar las fuerzas de reacción en el suelo gracias a una plataforma de fuerza. Este procedimiento, aunque es el más preciso, también es el más costoso.Wind tunnel measurements consist of generating an air flow over the cyclist / bicycle assembly and quantifying the reaction forces on the ground thanks to a force platform. This procedure, while the most accurate, is also the most expensive.
Se observará que en términos de aerodinamismo, es muy conocido determinar unas características aerodinámicas y/o mecánicas, tales como la deformación de un objeto o la resistencia al avance de un objeto, por ejemplo, mediante unas pruebas en túnel aerodinámico. Es el caso, en particular, de los documentos US 7,997,130 y US 2007/095135. It will be noted that in terms of aerodynamics, it is well known to determine aerodynamic and / or mechanical characteristics, such as the deformation of an object or the drag of an object, for example, by means of wind tunnel tests. This is the case, in particular, of documents US 7,997,130 and US 2007/095135.
El documento US 7,997,130 describe un sistema y un procedimiento para medir la deformación de un objeto, tal como un avión de combate por ejemplo, posicionado en un túnel aerodinámico. El objeto es posicionado en el túnel aerodinámico y se registra un sistema de adquisición de una nube de puntos que representa por lo menos una superficie del objeto. El objeto se desplaza en el flujo del túnel aerodinámico y se adquieren asimismo unas nubes de puntos de la superficie del objeto durante dicho desplazamiento con el fin de determinar por medio de un sistema informático, que recibe los datos relativos a las nubes de puntos, por lo menos la posición del objeto, su orientación y la deformación de la superficie del objeto.Document US 7,997,130 describes a system and a method for measuring the deformation of an object, such as a fighter plane for example, positioned in a wind tunnel. The object is positioned in the wind tunnel and an acquisition system is recorded for a point cloud representing at least one surface of the object. The object moves in the flow of the wind tunnel and point clouds are also acquired from the surface of the object during said movement in order to determine by means of a computer system, which receives the data related to the point clouds, by at least the position of the object, its orientation and the deformation of the object's surface.
El documento US 2007/095135 describe un procedimiento de determinación de la resistencia al avance de una aeronave que comprende proporcionar un modelo de la aeronave, que es posicionado en un túnel aerodinámico en una orientación inicial determinada, y una pluralidad de sensores de orientación y de inclinación montados en dicho modelo. El modelo se desplaza a continuación en el flujo del túnel aerodinámico desde su orientación inicial hacia una segunda orientación con el fin de determinar la resistencia al avance de la aeronave en las diferentes orientaciones posibles de esta última.Document US 2007/095135 describes a method for determining the drag of an aircraft that comprises providing a model of the aircraft, which is positioned in a wind tunnel in a determined initial orientation, and a plurality of orientation and bearing sensors. inclination mounted on said model. The model is then moved in the wind tunnel flow from its initial orientation to a second orientation in order to determine the drag of the aircraft in the different possible orientations of the latter.
Sin embargo, además del hecho de que estos procedimientos son particularmente costosos teniendo en cuenta la necesidad de utilizar un túnel aerodinámico, no están adaptados para determinar la resistencia aerodinámica de un ciclista en particular.However, in addition to the fact that these procedures are particularly expensive considering the need to use a wind tunnel, they are not adapted to determine the aerodynamic drag of a particular cyclist.
Recientemente, se ha propuesto un nuevo procedimiento de medición de la resistencia aerodinámica. Este procedimiento se describe en particular en las publicaciones de T. Defraeye, B. Blocken, E. Koninckx, P. Hespel, y J. Carmeliet. "Aerodynamic study of different cyclist positions: CFD analysis and full-scale wind-tunnel tests." J Biomech, 43(7): 1262-1268, mayo de 2010 y Peter Nicholas Doval. "Aerodynamic Analysis and Drag Coefficient Evaluation of Time-Trial Bicycle Riders." PhD thesis, University of Wisconsin-Milwaukee, 2012.Recently, a new aerodynamic drag measurement procedure has been proposed. This procedure is described in particular in the publications of T. Defraeye, B. Blocken, E. Koninckx, P. Hespel, and J. Carmeliet. "Aerodynamic study of different cyclist positions: CFD analysis and full-scale wind-tunnel tests." J Biomech, 43 (7): 1262-1268, May 2010 and Peter Nicholas Doval. "Aerodynamic Analysis and Drag Coefficient Evaluation of Time-Trial Bicycle Riders." PhD thesis, University of Wisconsin-Milwaukee, 2012.
Este procedimiento consiste en acoplar un modelo 3D del ciclista y de su bicicleta obtenido gracias a un escáner con un código de cálculo numérico de dinámica de fluidos. Este procedimiento permite obtener unas resistencias al avance aerodinámicas de acuerdo con las mediciones en un túnel aerodinámico. Tiene la ventaja de estar realizado con menor coste y sin necesidad de colocarse en unas condiciones reales. Sin embargo, su principal inconveniente es que utiliza un modelo estático que no es representativo de un ciclista pedaleando.This procedure consists of coupling a 3D model of the cyclist and his bicycle obtained thanks to a scanner with a numerical calculation code of fluid dynamics. This procedure makes it possible to obtain aerodynamic drag resistances according to measurements in an aerodynamic tunnel. It has the advantage of being made at a lower cost and without the need to be placed in real conditions. However, its main drawback is that it uses a static model that is not representative of a cyclist pedaling.
Divulgación de la invenciónDisclosure of the invention
Uno de los objetivos de la invención es por lo tanto remediar estos inconvenientes proporcionando un procedimiento y un sistema de determinación de la resistencia aerodinámica de un objeto o de una persona en movimiento de diseño sencillo y poco oneroso, que permita medir sin contacto la resistencia aerodinámica de un objeto o de una persona, tal como un ciclista en movimiento, sin que dicho objeto o dicha persona se desplace con respecto al suelo. Se entiende por persona en movimiento una persona que tiene por un lado un movimiento de desplazamiento con respecto a un fluido tal como el aire por ejemplo, y por otro lado un movimiento propio, en el caso del ciclista un movimiento de pedaleo.One of the objectives of the invention is therefore to remedy these drawbacks by providing a method and a system for determining the aerodynamic resistance of an object or a person in motion with a simple and inexpensive design, which allows the aerodynamic resistance to be measured without contact. of an object or person, such as a cyclist in motion, without said object or person moving relative to the ground. A person in motion is understood as a person who has, on the one hand, a movement of movement with respect to a fluid such as air for example, and on the other hand, a proper movement, in the case of the cyclist, a pedaling movement.
Con este fin y de acuerdo con la invención, se propone un procedimiento de determinación de la resistencia aerodinámica de un objeto o de una persona que tiene un movimiento propio destacable por que comprende por lo menos las etapas siguientes:To this end and in accordance with the invention, a method for determining the aerodynamic resistance of an object or of a person that has its own remarkable movement is proposed because it comprises at least the following steps:
- adquirir entre dos instantes t¡nic y tfin unas nubes de puntos 3D que corresponden a las posiciones en el espacio de diferentes puntos posicionados en diferentes lugares del objeto o de la persona,- acquire between two instants t¡nic and tfin some clouds of 3D points that correspond to the positions in space of different points positioned in different places of the object or person,
- determinar unos modelos 3D en función del tiempo, para cada instante entre tinic y tfin, a partir de los puntos de las nubes de puntos adquiridos en cada instante entre los instantes tinic y tfin,- determine 3D models as a function of time, for each instant between tinic and tfin, from the points of the point clouds acquired at each instant between tinic and tfin,
- simular numéricamente el desplazamiento de los sólidos descritos por las nubes de puntos que corresponden al objeto o a la persona en un fluido, para cada instante entre los instantes tinic y tfin, estando dicha simulación en particular parametrizada por la dirección y la velocidad de desplazamiento del fluido con respecto al sólido representado por un vector D,- numerically simulate the displacement of the solids described by the point clouds that correspond to the object or the person in a fluid, for each instant between the instants tinic and tfin, said simulation in particular being parameterized by the direction and speed of displacement of the fluid with respect to the solid represented by a vector D,
- calcular la media de las fuerzas de resistencia aerodinámica ejercidas sobre el objeto o la persona en movimiento entre dichos instantes tinic y tfin a partir de dicha simulación numérica del desplazamiento de los sólidos.- calculating the average of the aerodynamic resistance forces exerted on the object or the person in motion between said instants tinic and tfin from said numerical simulation of the displacement of the solids.
Preferentemente, dicha etapa de adquisición de nubes de puntos consiste por lo menos en las etapas siguientes:Preferably, said point cloud acquisition stage consists of at least the following stages:
- proyectar un patrón, tal como un Moiré por ejemplo, sobre el objeto o la persona,- project a pattern, such as a Moiré for example, onto the object or person,
- adquirir varias imágenes del patrón proyectado entre los instantes tinic y tfin, - acquire several images of the pattern projected between the instants tinic and tfin,
- calcular por triangulación y registrar la posición de puntos que pertenecen al patrón, entre los instantes t¡níc y tfin. Por otro lado, la etapa de determinación de modelos 3D en función del tiempo a partir de los puntos de las nubes de puntos adquiridos en cada instante entre los instantes tinic y tfin consiste por lo menos en las etapas siguientes: - aplicar por lo menos un preprocesamiento de los datos relativos a los conjuntos de nubes de puntos registrados,- calculate by triangulation and record the position of points that belong to the pattern, between the instants t¡níc and tfin. On the other hand, the stage of determining 3D models as a function of time from the points of the point clouds acquired at each instant between the instants tinic and tfin consists of at least the following stages: - apply at least one pre-processing of data relating to the registered point cloud sets,
- determinar un enmallado a partir de los datos relativos a los conjuntos de nubes de puntos registrados y procesados.- determine a mesh from the data related to the collected and processed point cloud sets.
Dicha etapa de preprocesamiento de los datos relativos a los conjuntos de puntos registrados consiste por lo menos en una etapa de filtrado de los puntos atípicos y/o de filtrado del ruido y/o de suavizado espaciotemporal. Además, la etapa de simulación del desplazamiento de los sólidos descritos por las nubes de puntos que corresponden al objeto o a la persona en un fluido, para cada instante entre los instantes tinic y tfin, consiste por lo menos en las etapas siguientes:Said stage of pre-processing the data relating to the sets of recorded points consists of at least one stage of filtering the outliers and / or filtering the noise and / or space-time smoothing. In addition, the stage of simulating the displacement of the solids described by the point clouds corresponding to the object or the person in a fluid, for each instant between the instants tinic and tfin, consists of at least the following stages:
- descomponer, en una base ortogonal (i,j,k), el vector D que corresponde al desplazamiento relativo del sólido en el fluido como la diferencia entre el vector V que corresponde al vector característico del movimiento del fluido y el vector C que corresponde al vector característico del movimiento del objeto o de la persona,- decompose, on an orthogonal basis (i, j, k), the vector D that corresponds to the relative displacement of the solid in the fluid as the difference between the vector V that corresponds to the characteristic vector of the movement of the fluid and the vector C that corresponds to the characteristic vector of the movement of the object or person,
- determinar el vector D en función de C, Vx, Vy y Vz, siendo C la velocidad del objeto o de la persona, considerando que el movimiento del objeto o de la persona se realiza según el eje i, Vx la velocidad del viento según el eje i, Vy la velocidad del viento según el eje j y Vz velocidad del viento según el eje k, - determinar la ley de probabilidad uniforme de las variables C, Vx, Vy y Vz, siendo dichas variables consideradas como variables aleatorias no independientes, que corresponden a la densidad P(C,Vx,Vy,Vz). Dicha media de las fuerzas de resistencia aerodinámica se determina calculando la integral cuádruple del producto de la densidad P(C,Vx,Vy,Vz) y de las fuerzas de resistencia aerodinámica ejercidas sobre el objeto o la persona F(D(C,Vx,Vy,Vz)).- determine the vector D as a function of C, Vx, Vy and Vz, where C is the speed of the object or person, considering that the movement of the object or person is carried out along the i axis, Vx the wind speed according to the i axis, Vy the wind speed along the j axis and Vz the wind speed along the k axis, - determine the law of uniform probability of the variables C, Vx, Vy and Vz, these variables being considered as non-independent random variables, which correspond to the density P (C, Vx, Vy, Vz). Said mean of the aerodynamic drag forces is determined by calculating the quadruple integral of the product of the density P (C, Vx, Vy, Vz) and of the aerodynamic drag forces exerted on the object or person F (D (C, Vx , Vy, Vz)).
Ventajosamente, la media de las fuerzas de resistencia aerodinámica se puede calcular considerando que: - el flujo del fluido es paralelo al plano horizontal, es decir que la componente Vx es nula, y/oAdvantageously, the mean of the aerodynamic resistance forces can be calculated considering that: - the flow of the fluid is parallel to the horizontal plane, that is to say that the component Vx is zero, and / or
- la fuerza del viento V no depende de su dirección, y/o- the force of the wind V does not depend on its direction, and / or
- la velocidad del objeto o de la persona que tiene un movimiento propio no depende de la fuerza V del viento ni de su dirección, y/o- the speed of the object or person that has its own motion does not depend on the force V of the wind or its direction, and / or
- el movimiento propio del objeto o de la persona es simétrico con respecto al eje definido por i, y/o - las nubes de puntos son registradas en un intervalo de tiempo igual al período del movimiento propio del objeto o de la persona.- the proper motion of the object or person is symmetric with respect to the axis defined by i, and / or - the point clouds are recorded in a time interval equal to the period of the proper motion of the object or person.
Otro objeto de la invención se refiere a un sistema de determinación de la resistencia aerodinámica de un objeto o de una persona que tiene un movimiento propio destacable por que comprende por lo menos:Another object of the invention refers to a system for determining the aerodynamic resistance of an object or of a person that has its own remarkable movement because it comprises at least:
- unos medios de adquisición entre dos instantes tinic y tfin de nubes de puntos 3D que corresponden a la posición en el espacio de diferentes puntos posicionados en diferentes lugares del objeto o de la persona, - unos medios de determinación de modelos 3D en función del tiempo, para cada instante entre tinic y tfin, a partir de los puntos de las nubes de puntos adquiridos en cada instante entre tinic y tfin,- means of acquisition between two instants tinic and tfin of 3D point clouds corresponding to the position in space of different points positioned in different places of the object or of the person, - means of determining 3D models as a function of time , for each instant between tinic and tfin, from the points of the point clouds acquired at each instant between tinic and tfin,
- unos medios de simulación numérica del desplazamiento de los sólidos descritos por las nubes de puntos que corresponden al objeto o a la persona en un fluido, para cada instante entre los instantes tinic y tfin, estando dicha simulación parametrizada en particular por la dirección y la velocidad de desplazamiento del fluido con respecto al sólido representadas por un vector D,- means of numerical simulation of the displacement of the solids described by the point clouds that correspond to the object or the person in a fluid, for each instant between the instants tinic and tfin, said simulation being parameterized in particular by the direction and the speed displacement of the fluid with respect to the solid represented by a vector D,
- unos medios de cálculo de la media de las fuerzas de resistencia aerodinámica ejercidas sobre el objeto o la persona en movimiento entre dichos instantes tinic y tfin a partir de dicha simulación numérica del desplazamiento de los sólidos. - means for calculating the average of the aerodynamic resistance forces exerted on the object or the person in motion between said instants tinic and tfin from said numerical simulation of the displacement of the solids.
Preferentemente, dichos medios de adquisición de nubes de puntos comprenden por lo menos:Preferably, said point cloud acquisition means comprise at least:
- unos medios de proyección de un patrón, tal como un Moiré por ejemplo, sobre el objeto o la persona, - unos medios de adquisición de varias imágenes del patrón proyectado entre los instantes t¡nic y tfin, - unos medios de cálculo por triangulación y de registro de la posición de puntos que pertenecen al patrón, entre los instantes tinic y tfin.- means for projecting a pattern, such as a Moiré for example, onto the object or person, - means for acquiring several images of the pattern projected between the instants t¡nic and tfin, - means for calculating by triangulation and registration of the position of points that belong to the pattern, between the instants tinic and tfin.
Dichos medios de determinación de modelos 3D en función del tiempo a partir de los puntos de las nubes de puntos adquiridos en cada instante entre los instantes tinic y tfin comprenden por lo menos:Said means of determining 3D models as a function of time from the points of the point clouds acquired at each instant between the instants tinic and tfin comprise at least:
- unos medios de preprocesamiento de los datos relativos a los conjuntos de nubes de puntos registrados, - unos medios de determinación de un enmallado a partir de los datos relativos a los conjuntos de nubes de puntos registrados y procesados.- means for pre-processing the data relating to the registered and processed point cloud sets, - means for determining a mesh from the data relating to the registered and processed point cloud sets.
Por otro lado, los medios de preprocesamiento de los datos relativos a los conjuntos de puntos registrados consisten por lo menos en un filtro de los puntos atípicos y/o un filtro de ruido y/o unos medios de suavizado espaciotemporal.On the other hand, the means for preprocessing the data relating to the sets of recorded points consist at least of a filter of the outliers and / or a noise filter and / or a means of spatio-temporal smoothing.
Dichos medios de simulación del desplazamiento de los sólidos formados por las nubes de puntos que corresponden al objeto o a la persona en un fluido, para cada instante entre los instantes tinic y tfin, comprenden por lo menos:Said means of simulating the displacement of the solids formed by the point clouds that correspond to the object or the person in a fluid, for each instant between the instants tinic and tfin, comprise at least:
- unos medios de descomposición, en una base ortogonal (i,j,k), del vector D que corresponde al desplazamiento relativo del sólido en el fluido como la diferencia entre el vector V que corresponde al vector característico del movimiento del fluido y el vector C que corresponde al vector característico del movimiento del objeto o de la persona,- means of decomposition, on an orthogonal basis (i, j, k), of the vector D that corresponds to the relative displacement of the solid in the fluid as the difference between the vector V that corresponds to the characteristic vector of the movement of the fluid and the vector C corresponding to the characteristic vector of the movement of the object or person,
- unos medios de determinación del vector D en función de C, Vx, Vy y Vz, siendo C la velocidad del objeto o de la persona, considerando que el movimiento del objeto o de la persona se realiza según el eje i, Vx la velocidad del viento según el eje i, Vy la velocidad del viento según el eje j y Vz la velocidad del viento según el eje k,- means for determining vector D as a function of C, Vx, Vy and Vz, where C is the speed of the object or person, considering that the movement of the object or person is carried out along the axis i, Vx the speed of the wind along the i axis, Vy the wind speed along the j axis and Vz the wind speed along the k axis,
- unos medios de determinación de la ley de probabilidad uniforme de las variables C, Vx, Vy y Vz, siendo dichas variables consideradas como unas variables aleatorias no independientes, que corresponden a la densidad P (C,Vx,Vy,Vz).- means for determining the law of uniform probability of the variables C, Vx, Vy and Vz, said variables being considered as non-independent random variables, corresponding to the density P (C, Vx, Vy, Vz).
Dicha media de las fuerzas de resistencia aerodinámica se determina calculando la integral cuádruple del producto de la densidad P(C,Vx,Vy,Vz) y de las fuerzas de resistencia aerodinámica ejercidas sobre el objeto o la persona F(D(C,Vx,Vy,Vz)).Said mean of the aerodynamic drag forces is determined by calculating the quadruple integral of the product of the density P (C, Vx, Vy, Vz) and of the aerodynamic drag forces exerted on the object or person F (D (C, Vx , Vy, Vz)).
La función F se obtiene mediante simulación numérica. Esta simulación numérica consiste en simular un flujo de aire alrededor del objeto en movimiento, cuya dirección está definida por el parámetro D. Esta simulación permite medir así virtualmente las fuerzas de resistencia aerodinámicas ejercidas sobre el objeto en movimiento. Dicha simulación numérica se basa en las ecuaciones de Navier-Stokes. Para mejores resultados, dicha simulación puede utilizar las ecuaciones de Reynolds-Averaged Navier-Stokes (RANS) y más particularmente el modelo SST (Shear-Stress Transport), presentado en Menter, F. R. (agosto 1994), "Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications", AIAA Journal 32 (8): 1598-1605.The function F is obtained by numerical simulation. This numerical simulation consists of simulating an air flow around the moving object, the direction of which is defined by parameter D. This simulation thus makes it possible to virtually measure the aerodynamic drag forces exerted on the moving object. This numerical simulation is based on the Navier-Stokes equations. For better results, such simulation can use the Reynolds-Averaged Navier-Stokes (RANS) equations and more particularly the SST (Shear-Stress Transport) model, presented in Menter, FR (August 1994), "Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications ", AIAA Journal 32 (8): 1598-1605.
Por otro lado, la media de las fuerzas de resistencia aerodinámica se puede calcular considerando que:On the other hand, the mean of the aerodynamic drag forces can be calculated considering that:
- el flujo del fluido es paralelo al plano horizontal, es decir que la componente Vx es nula, y/o- the fluid flow is parallel to the horizontal plane, that is, the component Vx is zero, and / or
- la fuerza del viento V no depende de su dirección, y/o- the force of the wind V does not depend on its direction, and / or
- la velocidad del objeto o de la persona que tiene un movimiento propio no depende de la fuerza V del viento ni de su dirección, y/o- the speed of the object or person that has its own motion does not depend on the force V of the wind or its direction, and / or
- el movimiento propio del objeto o de la persona es simétrico con respecto al eje definido por i, y/o - las nubes de puntos se registran en un intervalo de tiempo igual al período del movimiento propio del objeto o de la persona. - the proper motion of the object or person is symmetrical with respect to the axis defined by i, and / or - the point clouds are recorded in a time interval equal to the period of the proper motion of the object or person.
Breve descripción de los dibujosBrief description of the drawings
Otras ventajas y características se desprenderán mejor de la descripción siguiente, de varias variantes de realización, dadas a título de ejemplos no limitativos, del procedimiento de determinación de la resistencia aerodinámica de un objeto o de una persona en movimiento según la invención, a partir de los dibujos adjuntos, en los que:Other advantages and characteristics will emerge better from the following description, of various embodiment variants, given by way of non-limiting examples, of the method for determining the aerodynamic resistance of an object or a person in motion according to the invention, based on the accompanying drawings, in which:
- la figura 1 es una representación esquemática de las distintas etapas del procedimiento según la invención, - figure 1 is a schematic representation of the different stages of the process according to the invention,
- la figura 2 es una representación esquemática de un sistema de adquisición de nubes de puntos 3D que corresponden a la posición de puntos en la superficie del sólido en movimiento para la realización del procedimiento según la invención.Figure 2 is a schematic representation of a 3D point cloud acquisition system corresponding to the position of points on the surface of the moving solid for carrying out the method according to the invention.
Modo de realización de la invenciónMode of carrying out the invention
Se describirá a continuación un procedimiento y un sistema de determinación de la resistencia aerodinámica de un ciclista que pedalea en una bicicleta, teniendo dicho ciclista un movimiento propio, es decir su movimiento de pedaleo; sin embargo, resulta evidente que el procedimiento y el sistema se podrán adaptar a cualquier objeto o a cualquier persona que tiene un movimiento propio sin apartarse por ello del marco de la invención.A method and a system for determining the aerodynamic resistance of a cyclist pedaling a bicycle will now be described, said cyclist having his own movement, that is, his pedaling movement; however, it is clear that the method and the system can be adapted to any object or to any person that has its own movement without departing from the framework of the invention.
Con referencia a la figura 1, el sistema según la invención que permite medir sin contacto y sin desplazamiento la resistencia aerodinámica de un sólido en movimiento comprende varios subsistemas que interactúan. Más particularmente, comprende un sistema que permite estudiar el movimiento del sólido sin que este se desplace con respecto al suelo. Así, según la invención, la medición del aerodinamismo no necesita un movimiento real del sólido estudiado con respecto al suelo. Por otro lado, esta medición tiene en cuenta el aspecto dinámico del movimiento del sólido. Por lo tanto, es necesario colocar el sólido estudiado en unas condiciones que permitan un movimiento sin tensión, pero sin desplazamiento con respecto al suelo. En la aplicación particular de la determinación de la resistencia aerodinámica de un ciclista pedaleando en una bicicleta, dicho sistema consiste por ejemplo en un home trainer de rodillos o similar.With reference to FIG. 1, the system according to the invention that makes it possible to measure without contact and without displacement the aerodynamic resistance of a solid in motion comprises several interacting subsystems. More particularly, it comprises a system that makes it possible to study the movement of the solid without it moving with respect to the ground. Thus, according to the invention, the measurement of aerodynamics does not require a real movement of the studied solid with respect to the ground. On the other hand, this measurement takes into account the dynamic aspect of the movement of the solid. Therefore, it is necessary to place the solid studied in conditions that allow a movement without tension, but without displacement with respect to the ground. In the particular application of determining the aerodynamic resistance of a cyclist pedaling a bicycle, said system consists, for example, of a home trainer with rollers or the like.
Por otro lado, la invención comprende asimismo un sistema que permite medir en el tiempo unas nubes de puntos 3D que corresponden a la posición de los puntos en la superficie del sólido en movimiento, un sistema que permite obtener un modelo 3D+t correcto del sólido en movimiento en el que se interesa y un sistema que permite simular el desplazamiento del sólido en movimiento en un fluido y medir las fuerzas ejercidas sobre este sólido.On the other hand, the invention also comprises a system that makes it possible to measure in time 3D point clouds that correspond to the position of the points on the surface of the moving solid, a system that makes it possible to obtain a correct 3D + t model of the solid in motion in which you are interested and a system that allows you to simulate the displacement of the solid in motion in a fluid and measure the forces exerted on this solid.
Dicho sistema que permite medir en el tiempo unas nubes de puntos 3D que corresponden a la posición de puntos sobre la superficie del sólido en movimiento está constituido por un conjunto de sensores 3D sin contacto habituales tales como unas cámaras de tiempo de vuelo, unas cámaras estereoscópicas, unas cámaras de proyección de patrones, o similares, colocadas convenientemente de manera que se obtenga un conjunto de puntos 3D que describen totalmente el sólido, en este caso el ciclista, estudiado en un tiempo t. Este sistema es capaz de adquirir unas nubes de puntos 3D a una cadencia compatible con la dinámica del movimiento estudiado. Se obtiene así un conjunto de nubes de puntos 3D que representan el sólido en diferentes tiempos t.Said system that allows 3D point clouds that correspond to the position of points on the surface of the moving solid to be measured in time is made up of a set of common non-contact 3D sensors such as time-of-flight cameras, stereoscopic cameras , pattern projection cameras, or the like, conveniently placed so as to obtain a set of 3D points that fully describe the solid, in this case the cyclist, studied at time t. This system is capable of acquiring 3D point clouds at a rate compatible with the dynamics of the movement studied. Thus, a set of 3D point clouds is obtained that represent the solid at different times t.
Con referencia a la figura 2, dicho sistema que permite medir en el tiempo unas nubes de puntos 3D que corresponden a la posición de puntos sobre la superficie del sólido en movimiento está constituido, por ejemplo, por cámaras 1 que proyectan unos patrones sobre un ciclista 2 que pedalea en una bicicleta 3 posicionada en un home trainer de rodillos 4. Dichas cámaras 1 comprenden asimismo unos sensores CCD aptos para registrar los patrones proyectados por las cámaras y reflejados por el ciclista 2. Estas cámaras 1 están conectadas a una unidad de procesamiento 5, que consiste por ejemplo en un ordenador de sobremesa denominado PC o similar, en la que los datos que corresponden a las nubes de puntos son registrados y después procesados por un programa informático que ejecuta las diferentes etapas del sistema de acuerdo con la invención.With reference to Figure 2, said system that allows 3D point clouds that correspond to the position of points on the surface of the moving solid to be measured in time is constituted, for example, by cameras 1 that project patterns onto a cyclist 2 pedaling on a bicycle 3 positioned on a roller home trainer 4. Said cameras 1 also comprise CCD sensors capable of recording the patterns projected by the cameras and reflected by the cyclist 2. These cameras 1 are connected to a processing unit 5, consisting for example of a desktop computer called a PC or the like, in which the data corresponding to the point clouds are recorded and then processed by a computer program that executes the different stages of the system according to the invention.
Dicho sistema que permite obtener un modelo 3D+t correcto del sólido en movimiento permite procesar el conjunto de nubes de puntos 3D con el fin de obtener un modelo 3D+t del sólido en movimiento. Para ello, se aplica un conjunto de preprocesamientos a estas nubes de puntos: filtrado de los puntos atípicos, filtrado del ruido, suavizado espaciotemporal, etc. Por último, se aplica una operación de enmallado con el fin de obtener un modelo 3D+t explotable por el sistema de simulación de mecánica de fluidos.Said system that allows obtaining a correct 3D + t model of the moving solid allows processing the set of 3D point clouds in order to obtain a 3D + t model of the moving solid. To do this, a set of preprocesses are applied to these point clouds: filtering outliers, filtering noise, space-time smoothing, etc. Finally, a meshing operation is applied in order to obtain a 3D + t model that can be exploited by the fluid mechanics simulation system.
El sistema que permite simular el desplazamiento del sólido en movimiento en un fluido y medir las fuerzas ejercidas sobre el ciclista utiliza unos procedimientos de mecánica de fluidos clásicos numérica con el fin de simular el desplazamiento del sólido estudiado en un fluido (aire en el caso del ciclista). Esta simulación está parametrizada en particular por la dirección y la velocidad de desplazamiento del fluido con respecto al sólido, representadas por el vector D. Permite in fine calcular la media F(D) de las fuerzas de resistencia aerodinámica ejercidas sobre el sólido durante toda la secuencia de movimiento registrada por los sensores.The system that allows simulating the displacement of the solid in motion in a fluid and measuring the forces exerted on the cyclist uses classical numerical fluid mechanics procedures in order to simulate the displacement of the solid studied in a fluid (air in the case of the cyclist). This simulation is parametrized in particular by the direction and the speed of displacement of the fluid with respect to the solid, represented by the vector D. It allows in fine to calculate the average F (D) of the aerodynamic resistance forces exerted on the solid during all the motion sequence recorded by sensors.
En la medida en que se busca estudiar el aerodinamismo en unas condiciones reales (como por ejemplo las de una práctica del ciclismo en exterior), es importante tener en cuenta el hecho de que la velocidad relativa del sólido con respecto al fluido puede cambiar en el tiempo. Con este fin, el procedimiento según la invención propone modelizar las diferentes variables del movimiento mediante unas leyes de probabilidad que permitan describir un entorno determinado (por ejemplo un recorrido ciclista particular).To the extent that it is sought to study aerodynamics in real conditions (such as those of cycling outdoors), it is important to take into account the fact that the relative velocity of the solid with respect to the fluid can change over time. To this end, the method according to the invention proposes to model the different variables of the movement by means of probability laws that make it possible to describe a specific environment (for example a particular cycling route).
Para ello, se proporciona el espacio de una base ortonormal (i; j; k) tal que el sólido estudiado se desplaza según el eje definido por i y que el suelo es paralelo al plano definido por i y j, por ejemplo.For this, the space of an orthonormal base (i; j; k) is provided such that the studied solid moves along the axis defined by i and that the ground is parallel to the plane defined by i and j, for example.
Se puede descomponer entonces el vector D, desplazamiento relativo del sólido en el fluido de la manera siguiente:The vector D, relative displacement of the solid in the fluid, can then be decomposed as follows:
en la que V designa el vector característico del viento y C designa el vector característico del movimiento del ciclista.in which V designates the characteristic vector of the wind and C designates the characteristic vector of the movement of the cyclist.
Se puede escribir a continuación:It can be written below:
ú or - r 7 ' - - r 7 '- v j v j ^ i ^ i , i1 - c í , i 1 - c í
í> í> - ¡r, r r j - - r, r r j - \\í\\ í
Se puede escribir por lo tanto D como una función de C, Vx, Vy y Vz, en la que C designa la velocidad del sólido, Vx designa la velocidad del viento según el eje i, Vy designa la velocidad del viento según el eje j, y Vz designa la velocidad del viento según el eje k.Therefore, D can be written as a function of C, Vx, Vy and Vz, in which C designates the speed of the solid, Vx designates the wind speed along the i axis, Vy designates the wind speed along the j axis , and Vz designates the wind speed along the k axis.
Se consideran a continuación las variables C, Vx, Vy y Vz como unas variables aleatorias no independientes cuya ley de probabilidad uniforme viene dada por la densidad P(C;Vx;Vy;Vz), permitiendo esta ley de probabilidad describir estadísticamente el entorno particular estudiado.The variables C, Vx, Vy and Vz are considered below as non-independent random variables whose uniform probability law is given by the density P (C; Vx; Vy; Vz), allowing this probability law to statistically describe the particular environment studied.
Considerando por un lado que la cinética del movimiento estudiado es muy rápida con respecto a los cambios de condiciones (es decir a una evolución de las variables C, Vx, Vy y Vz) y por otro lado, que la secuencia de movimiento no está influida por las variables C, Vx, Vy y Vz, se puede calcular entonces una media Fm de las fuerzas de resistencia aerodinámica F ejercidas sobre el sólido en el entorno estudiado de la manera siguiente:Considering on the one hand that the kinetics of the movement studied is very fast with respect to changes in conditions (that is, to an evolution of the variables C, Vx, Vy and Vz) and on the other hand, that the sequence of movement is not influenced Using the variables C, Vx, Vy and Vz, an average Fm of the aerodynamic drag forces F exerted on the solid in the studied environment can then be calculated as follows:
F se obtiene por simulación numérica. Esta simulación numérica consiste en simular un flujo de aire alrededor del objeto en movimiento, cuya dirección está definida por el parámetro D. Esta simulación permite medir así virtualmente las fuerzas de resistencia aerodinámica ejercidas sobre el objeto en movimiento. Dicha simulación numérica se basa en las ecuaciones de Navier-Stokes. Para obtener mejores resultados, dicha simulación puede utilizar las ecuaciones de Reynolds-Averaged Navier-Stokes (RANS) y más particularmente el modelo SST (Shear-Stress Transport), presentado en Menter, F. R. (agosto 1994), "Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications", AIAA Journal 32 (8): 1598-1605.F is obtained by numerical simulation. This numerical simulation consists of simulating an air flow around the moving object, the direction of which is defined by parameter D. This simulation thus makes it possible to virtually measure the aerodynamic drag forces exerted on the moving object. This numerical simulation is based on the Navier-Stokes equations. For better results, such simulation can use the Reynolds-Averaged Navier-Stokes (RANS) equations and more particularly the SST (Shear-Stress Transport) model, presented in Menter, FR (August 1994), "Two-Equation Eddy- Viscosity Turbulence Models for Engineering Applications ", AIAA Journal 32 (8): 1598-1605.
Resulta muy evidente que la medición virtual de las fuerzas de resistencia aerodinámica se podrá realizar mediante cualquier otro procedimiento bien conocido por el experto en la materia, en particular a partir de un cálculo directo en función de las mediciones que permiten la determinación del modelo 3D del objeto en movimiento, sin apartarse por ello del marco de la invención.It is very evident that the virtual measurement of the aerodynamic resistance forces can be carried out by any other method well known to the person skilled in the art, in particular from a direct calculation based on the measurements that allow the determination of the 3D model of the object in motion, without thereby departing from the scope of the invention.
El modelo es simulado utilizando un procedimiento numérico (procedimiento de los elementos finitos, procedimiento de los volúmenes finitos, procedimientos de las diferencias finitas, procedimiento espectral, etc.) utilizado en un ordenador.The model is simulated using a numerical procedure (finite element procedure, finite volume procedure, finite difference procedure, spectral procedure, etc.) used in a computer.
Con el fin de hacer que el modelo descrito anteriormente se pueda utilizar más fácilmente en un ordenador, ventajosamente este último se discretiza y se simplifica.In order to make the model described above more easily usable on a computer, the latter is advantageously discretized and simplified.
En el caso particular del estudio de la aerodinámica de un ciclista, se pueden aportar las simplificaciones siguientes: Viento paralelo a la carreteraIn the particular case of studying the aerodynamics of a cyclist, the following simplifications can be provided: Wind parallel to the road
Se puede considerar que el viento es paralelo a la carretera, o sea:It can be considered that the wind is parallel to the road, that is:
Vz = 0 Vz = 0
Dirección y fuerza del viento independienteIndependent wind direction and force
Se puede expresar Vx y Vy en función de la fuerza V del viento y del ángulo del viento con respecto al eje definido por i:Vx and Vy can be expressed as a function of the force V of the wind and the angle of the wind with respect to the axis defined by i:
VX= V cosaVX = V thing
Vy= V sinaVy = V sina
Se puede considerar que la fuerza V del viento no depende de su dirección a y a la inversa.It can be considered that the force V of the wind does not depend on its direction a and vice versa.
Velocidad del ciclista y viento independientesIndependent rider and wind speed
Se puede considerar que la velocidad C del ciclista no depende de la fuerza V del viento ni de su dirección a. De esta manera, utilizando las tres simplificaciones anteriores, se puede reescribir la ley de probabilidad uniforme en la forma simplificada siguiente:It can be considered that the speed C of the cyclist does not depend on the force V of the wind or its direction a. In this way, using the three previous simplifications, the law of uniform probability can be rewritten in the following simplified form:
La independencia de las variables aleatorias permite expresar la distribución uniforme como el producto de las leyes marginales asociadas a cada una de ellas. Estas leyes se pueden definir de manera mucho más simple que la ley uniforme. Por ejemplo, PV se define habitualmente como la densidad de probabilidad asociada a una ley de Weibull; Pa se define como la densidad de probabilidad asociada a una ley uniforme; y PC tiene un perfil que depende específicamente del tipo de prueba ciclista estudiada (sprint, carrera por etapas, etc.).The independence of the random variables allows the uniform distribution to be expressed as the product of the marginal laws associated with each of them. These laws can be defined much simpler than the uniform law. For example, PV is commonly defined as the probability density associated with a Weibull law; Pa is defined as the probability density associated with a uniform law; and PC has a profile that depends specifically on the type of cycling event studied (sprint, stage race, etc.).
Por lo tanto, es posible reescribir la función F, que corresponde a las fuerzas de resistencia aerodinámica ejercidas sobre el sólido en el entorno estudiado, de la forma siguiente:Therefore, it is possible to rewrite the function F, which corresponds to the aerodynamic resistance forces exerted on the solid in the studied environment, as follows:
Y teniendo en cuenta las simplificaciones, se puede escribir:And taking into account the simplifications, you can write:
Con el fin de limitar el número de simulaciones y discretizar el modelo, se puede considerar que el ciclista y su movimiento son simétricos con respecto al eje definido por i. Se tiene entonces:In order to limit the number of simulations and to discretize the model, the cyclist and his movement can be considered to be symmetric with respect to the axis defined by i. It then has:
F\Dic . r . o ); = Fí Di( r . - m ;• íF \ Dec. r. o); = Fí Di (r. - m ; • í
Puede bastar así con variar a de 0 a n cuando tienen lugar las simulaciones.It may thus be sufficient to vary a from 0 to n when the simulations take place.
Por otro lado, se puede considerar asimismo que el movimiento del ciclista es periódico (período es el tiempo de una vuelta de pedal). Por lo tanto, puede bastar con registrar una secuencia de movimiento de duración igual a este período.On the other hand, it can also be considered that the cyclist's movement is periodic (period is the time of one pedal turn). Therefore, it may be sufficient to record a movement sequence of duration equal to this period.
Para discretizar el modelo, se elige en un primer momento por un lado Cmax y Vmax y por otro lado, SC, SV y Sa tales que SC divide Cmax, SV divide Vmax y Sa = n/k, siendo k un entero.To discretize the model, we first choose Cmax and Vmax on the one hand and, on the other hand, SC, SV and Sa such that SC divides Cmax, SV divides Vmax and Sa = n / k, where k is an integer.
Se tiene entonces:It then has:
Se puede entonces aproximar Fm de la manera siguiente:We can then approximate Fm as follows:
con with
De esta manera, el modelo necesita solamente n = ncnvna simulaciones. Thus, the model needs only n = ncnvna simulations.
Claims (15)
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR1556815A FR3038980B1 (en) | 2015-07-17 | 2015-07-17 | METHOD AND SYSTEM FOR DETERMINING THE AERODYNAMIC RESISTANCE OF A CYCLIST |
| PCT/EP2016/066531 WO2017012923A1 (en) | 2015-07-17 | 2016-07-12 | Method and system for determining the aerodynamic resistance of a cyclist |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| ES2874575T3 true ES2874575T3 (en) | 2021-11-05 |
Family
ID=54015094
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| ES16741578T Active ES2874575T3 (en) | 2015-07-17 | 2016-07-12 | Procedure and system for determining the aerodynamic resistance of a cyclist |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US20180188134A1 (en) |
| EP (1) | EP3325936B1 (en) |
| AU (1) | AU2016294674B2 (en) |
| ES (1) | ES2874575T3 (en) |
| FR (1) | FR3038980B1 (en) |
| WO (1) | WO2017012923A1 (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113505544A (en) * | 2021-06-18 | 2021-10-15 | 清华大学 | Bicycle motion virtual numerical value wind tunnel system based on finite volume method |
| CN117893706B (en) * | 2024-01-04 | 2024-06-14 | 南京科技职业学院 | Roadway wind resistance calculation method based on point cloud scanning |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7254998B2 (en) * | 2005-10-31 | 2007-08-14 | The Boeing Company | Method for determining drag characteristics of aircraft and system for performing the method |
| US7997130B1 (en) * | 2009-03-27 | 2011-08-16 | The Boeing Company | System and method for measuring deformation of an object in a fluid tunnel |
| JP5054746B2 (en) * | 2009-10-01 | 2012-10-24 | 川崎重工業株式会社 | Wind tunnel relative distance measurement system and wind tunnel relative distance measurement method |
-
2015
- 2015-07-17 FR FR1556815A patent/FR3038980B1/en not_active Expired - Fee Related
-
2016
- 2016-07-12 AU AU2016294674A patent/AU2016294674B2/en active Active
- 2016-07-12 EP EP16741578.5A patent/EP3325936B1/en active Active
- 2016-07-12 ES ES16741578T patent/ES2874575T3/en active Active
- 2016-07-12 WO PCT/EP2016/066531 patent/WO2017012923A1/en active Application Filing
- 2016-07-12 US US15/743,224 patent/US20180188134A1/en not_active Abandoned
Also Published As
| Publication number | Publication date |
|---|---|
| WO2017012923A1 (en) | 2017-01-26 |
| EP3325936A1 (en) | 2018-05-30 |
| FR3038980B1 (en) | 2017-08-25 |
| AU2016294674B2 (en) | 2022-02-24 |
| US20180188134A1 (en) | 2018-07-05 |
| EP3325936B1 (en) | 2021-03-10 |
| FR3038980A1 (en) | 2017-01-20 |
| AU2016294674A1 (en) | 2018-02-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| TWI724168B (en) | Method and apparatus for monitoring fluid dynamic drag | |
| CN106394561B (en) | A kind of method of estimation and device of longitudinal speed of vehicle | |
| CN103632062B (en) | Method and device for determining uphill and downhill running states of vehicle | |
| De Roo et al. | A semi-empirical model of the energy balance closure in the surface layer | |
| CN110187143B (en) | Chromatography PIV reconstruction method and device based on deep neural network | |
| CN104598753A (en) | Bridge moving vehicle load recognition method based on Brakhage V method | |
| CN104008229A (en) | Method for building street pollutant concentration distribution model | |
| CN105421200B (en) | A kind of measurement method of planeness and detection means based on vibration detection | |
| ES2874575T3 (en) | Procedure and system for determining the aerodynamic resistance of a cyclist | |
| CN110114643A (en) | Axle load measuring device and axle load measuring method | |
| CN108839657A (en) | A method of online recognition road roughness information is responded based on automobile vibration | |
| Zou et al. | Response surface methodology and improved interval analysis method—For analyzing uncertainty in accident reconstruction | |
| Tengattini et al. | Validation of an outdoor coast-down test to measure bicycle resistance parameters | |
| ITTO20110255A1 (en) | METHOD OF DETERMINING THE STRESSES THAT MUST BE APPLIED TO A TIRE DURING A TEST OF COVERED DURATION CARRIED OUT IN A TEST BENCH | |
| Li et al. | Efficient method based on recurrent neural networks for pavement evenness detection | |
| Roa et al. | Onboard wind speed and road grade measurements for determination of model coefficients on cycling | |
| ES2564543T3 (en) | Method of determining the stress that should be applied to a tire during a high efficiency indoor strength test | |
| Kraft et al. | Parameter identification of multi-body railway vehicle models–Application of the adjoint state approach | |
| CN103968832B (en) | The analogy method of dynamic hangover star image centrode and device | |
| Lemaitre et al. | Experimental evaluation of a computer-vision based method to assess the aerodynamic drag of cyclists | |
| Dembski et al. | Development of a terrain severity measurement system utilizing optical lasers | |
| Enriquez | Subgrid-scale turbulence modeling for improved large-eddy simulation of the atmospheric boundary layer | |
| Alcázar-Vargas | Study of the fundamental parameters affecting tire models for use in active safety systems in vehicles | |
| Chalupa et al. | Simulation of track vehicle passability | |
| Andersen et al. | Introducing functional data analysis to coast down modeling for rolling resistance estimation |