Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Bond Graph

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 31

Escuela Técnica Superior de Ingeniería y Diseño Industrial

Universidad Politécnica de Madrid

INTRODUCCIÓN A LA TÉCNICA DE BOND-GRAPH

Prof. Dr. José A. Lozano Ruiz


Universidad Politécnica de Madrid
E.T.S. de Ingeniería y Diseño Industrial
Departamento de Ingeniería Mecánica, Química y Diseño Industrial
E-mail: joseantonio.lozano@upm.es

Palabras clave:
Bond-Graph, simulación, mecánica, dinámica, vehículos, automóviles.

Objetivos:
El objetivo de este curso o seminario es dar una visión general del uso de la técnica de
simulación de sistemas dinámicos Bond-Graph, aplicada al caso particular de la simulación de
mecanismos y vehículos automóviles.
Este objetivo general, por razones docentes se particulariza en los siguientes objetivos
concretos:
a) Introducción teórica a los conceptos básicos utilizados en la técnica de simulación de
sistemas dinámicos Bond-Graph. Se centrará la atención en los conceptos utilizados para
la simulación de la dinámica vehicular.
b) Aprendizaje aplicado a la programación de modelos de sistemas dinámicos desarrollados
mediante la técnica de Bond-Graph.
c) Obtención de las ecuaciones de comportamiento de los sistemas, a partir del diagrama de
Bond-Graph.
d) Realización de casos prácticos de aplicación de los conocimientos adquiridos,
consistentes en la simulación por ordenador de algunos casos de dinámica vehicular.
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

Contenidos
1) Introducción a la técnica de Bond-Graph.
2) Modelado de sistemas mecánicos aplicando la técnica de Bond-Graph. Ejercicios.
3) Programación y simulación de dinámica vehicular utilizando Bond-Graph.
4) Ejercicios. Conclusiones.

Metodología:
Se realizarán exposiciones orales teóricas de los conceptos básicos y técnicas utilizadas.
Después de cada exposición teórica se realizan ejercicios prácticos utilizando programas
informáticos de apoyo.
Será necesaria la utilización de equipos informáticos, tanto por parte del profesor como por
parte de los alumnos. Estos equipos informáticos son necesarios para la presentación de las
exposiciones orales, para el manejo de los programas informáticos utilizados de apoyo a las
clases, así como para la realización de los ejercicios prácticos.

A quien va dirigido:
Estudiantes, técnicos e ingenieros del sector del diseño industrial e ingeniería mecánica.

Duración y programación
El curso se programa en clases de 1,5 horas aproximadamente, con una duración total de 9
horas, estructuradas de la siguiente manera:
Clase 1: (1,5 horas). Presentación. Introducción a la técnica de Bond-Graph.
Clase 2: (1,5 horas) Aplicación de la técnica de Bond-Graph a la simulación de sistemas
mecánicos.
Clase 3: (1,5 horas) Ejercicio 1, dinámica vertical de un vehículo automóvil.
Clase 4: (1,5 horas) Ejercicio 2, dinámica longitudinal de un vehículo. Simulación de una
maniobra de frenada.
Clase 5: (1,5 horas) Ejercicio 3, simulación de un sistema ABS.

Pag. T-2 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

Bibliografía:
[1] Félez, J. (1989). “Un método para el análisis unificado de la cinemática y dinámica de sistemas
vehiculares basado en la técnica de los grafos de unión”. Tesis doctoral de la E.T. Superior de
Ingenieros Industriales de Zaragoza.
[2] Félez, J., Vera, C., San José, I., Cacho, R., (1990). “BONDYN: A Bond Graph Based
Simulation Program for Multibody Systems”. Transactions of de ASME, Journal of Dynamic
Systems, Measurement and Control, Vol. 112, 717-727.
[3] Karnopp, D., Margolis, D., Y Rosenberg, R., (1990). “System Dynamics: A Unified Approach”.
Wiley-Interscience, New York, NY, 2ª edición.
[4] Lozano, J.A. (1996). “Modelo Teórico para la reconstrucción de accidentes de tráfico por
ordenador”. Tesis Doctoral. Universidad Politécnica de Madrid.
[5] Margolis, D., (1986). “Introduction into Bond Graphs”. Third Seminar on Advanced Vehicle
Dynamics, Amalfi, mayo 1986.
[6] Mathews, J. H., (1992). “Numerical Methods for Mathematics, Science, and Engineering”.
Prentice-Hall International Editions, 2ª edición.
[7] Pacejka, H.B., (1986). “Bond Graphs in Vehicle Dynamics”. Third Seminar on Advanced
Vehicle Dynamics, Amalfi, mayo 1986.
[8] Petzold, L., (1982). “A Description of DASSL: A Differential/Algebraic System Solver”. Proc.
10th IMACS Congress, Montreal, vol 1, 430-432.
[9] Rosenberg, R. C., Wang, Y. Y., (1992). “Some structuring issues in modeling”. DSC-Vol. 41,
Automated Modeling, ASME.
[10] Van Dixhoorn, J. J., (1977). “Simulation of Bond Graphs on Minicomputers”. Journal of
Dynamic Systems, Measurement, and Control, marzo 1977.
[11] Vera, C., Aparicio, F., Félez, J., (1993). “Simulación de Sistemas Dinámicos mediante la
técnica del Bond Graph”. E.T. Superior de Ingenieros Industriales de Madrid.
[12] Vera, C., Aparicio, F., Félez, J., San José, I., (1992). “Bond Graph simulation with the
BONDYN program. Safety of Transportation”. Delft University of Technology.
[13] Vera, C., Aparicio, F., Félez, J., San José, I., (1993). “Bond Graph simulation with the
BONDYN program. Safety of Transportation”. Western Simulation Multiconference, La Jolla
(California).
[14] Vera, C., Aparicio, F., Peso, S., (1992). Theoretic model for the computer analysis of vehicle
collisions. Safety of Transportation. Delft University of Technology.
[15] Vera, C., Lozano, J.A., Aparicio, F., Félez, J., (1994). “SINRAT IV. A Program for Traffic
Accident Reconstruction”. International Conference on Bond-Graph Modeling ICBGM’94.
Las Vegas, U.S.A.
[16] Vera, C., Lozano, J.A., Félez, J., (1996). “Sistema informático de reconstrucción de
accidentes de tráfico”. II Symposium sobre Ingeniería de los Transportes. Madrid, 1996.
[17] Vera, C.; Félez, J.; San José, I.; Cacho, R., (1987). “Introducción del método de los grafos de
unión a la dinámica vehicular”. Anales de Ingeniería Mecánica, año 5, nº. 1.

Pag. T-3 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

INTRODUCCIÓN A LA TÉCNICA DE BOND-GRAPH

Indice

1.- INTRODUCCIÓN A LA TÉCNICA DE BOND-GRAPH ........................................................................................... 5


1.1.- Introducción ........................................................................................................................................................... 5
1.2.- Elementos básicos del Bond-Graph ....................................................................................................................... 6
1.2.1.- Elementos FUENTE ........................................................................................................................................ 7
1.2.2.- Elementos que almacenan energía. ................................................................................................................. 8
1.2.3.- Elementos de disipan energía. ....................................................................................................................... 10
1.2.4.- Nodos de UNIÓN .......................................................................................................................................... 11
1.2.5.- Elementos TRANSFORMER........................................................................................................................ 13
1.2.6.- Elementos GYRATOR .................................................................................................................................. 15
1.3.- La causalidad........................................................................................................................................................ 16
1.4.- Planteamiento de las ecuaciones del sistema........................................................................................................ 17
1.4.1.- Introducción .................................................................................................................................................. 17
1.4.2.- Desarrollo de las ecuaciones de estado del sistema en el diagrama de Bond-Graph. .................................... 19
1.4.3.- Ecuaciones de comportamiento del sistema a partir del diagrama de Bond-Graph. ..................................... 28

Pag. T-4 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

1.- INTRODUCCIÓN A LA TÉCNICA DE BOND-GRAPH

1.1.- Introducción

Bond-Graph es una técnica de simulación de sistemas dinámicos en cualquier dominio de la


física y de la energía. Esta técnica está basada en la representación del flujo de energía que se
produce en el seno de los sistemas e incluso entre ellos, mediante líneas o flechas, también llamadas
BOND, (figura 1.1.). La punta de la flecha del BOND, indica el sentido de transmisión de la
energía.
Esfuerzo e(t)
Regla nemo f de flujo , e de esfuerzo
Flujo f(t)
Figura 1.1. Bond.
Cada BOND lleva asociados dos valores instantáneos de esfuerzo e(t) y flujo f(t), cuyo
producto es la energía instantánea o potencia W(t):
W(t) = e(t) f(t) ec.1.1.
En el BOND-GRAPH a las variables e(t) y f(t), se les denomina variables del sistema y sus
valores definen el BOND. El significado físico de estas variables dependerá del dominio físico en
que se encuadre el sistema en estudio. Por ejemplo, en el caso de la mecánica, e(t) es la fuerza y f(t)
es la velocidad. En la tabla 1.1, se muestra el significado de estas variables para diferentes dominios
de la Física.
Además de las ya mencionadas, se utilizan dos variables más que se les denomina variables
energéticas o dinámicas. Estas dos variables son en momento P(t), y el desplazamiento q(t), que se
definen como la integral en el tiempo de la variable esfuerzo y la variable flujo, respectivamente:
τ
P( t ) = ∫ e( t )dt La suma de los esfuerzos ( integral ) de fuerza por
el tiempo es igual al momento lineal
ec. 1.2.
0

τ
q ( t ) = ∫ f ( t )dtEsto viene de que la velocidad es igual al espacio entre el tiempo, luego la
ec. 1.3.
integral de la velocidad con respecto al tiempo da el desplazamiento
0

En el dominio de la mecánica, P(t) es la cantidad de movimiento y q(t) el desplazamiento del


sistema. En la tabla 1.1 también se muestra el significado de las variables P(t) y q(t) para diferentes
Cantidad movimiento es igual a masa por velocidad
dominios de la Física.
Por otra parte, la energía transmitida E(t) por el BOND, en un periodo de tiempo t, es la
integral en el tiempo de la potencia W(t), por lo que:
τ τ
E( t ) = ∫ W ( t )dt = ∫ e( t )f ( t )dt ec. 1.4.
0 0

Energía = Potencia * tiempo Entonces si sumas el producto de la potencia por el tiempo obtienes la energía,
porque la integral representa un sumatorio

Sentido de las flechas: representa el sentido de la energía. Es positiva la


que entra al sistema y la que sale de los elementos que almacenan

Pag. T-5 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

Tabla 1.1. Significado de las variables del Bond-Graph para diferentes dominios de la Física.

Sistema Esfuerzo Flujo Momento Desplazamiento


Físico e(t) f(t) P(t) q(t)
Fuerza por tiempo
Mecánica Fuerza Velocidad Momento lineal Desplazamiento
Traslación F (N) V (m/s) P (N.s) x (m)

Mecánica Par Velocidad Momento Ángulo girado


Rotación M (N.m) angular w angular H θ (rad)
(rad/s) (N.m.s)

Presión Caudal Momento Volumen


Hidráulica 2 3 2 3
p (N/m ) Q (m /s) P (N.s/m ) V (m )

Voltaje Intensidad Flujo Carga


Electricidad
e (V) I (A) F (V.s) q (C)

Sistemas Temperatura Flujo de calor Calor


Térmicos T (K) q (w) Q (cal.)

Presión Velocidad Momento Volumen


Acústica 2 3 2 3
p (N/m ) v (m /s) P (N.s/m ) V (m )

El momento siempre se calcula como el esfuerzo por el


tiempo

1.2.- Elementos básicos del Bond-Graph

A continuación se presentan un conjunto de elementos básicos que podrán ser usados para
modelizar los sistemas físicos mediante la técnica de Bond-Graph. Estos elementos aparecerán
como componentes elementales del sistema físico, representando su idealización matemática.
La técnica de Bond-Graph permite que, con la clasificación de variables presentada en el
apartado anterior sólo se necesite una pequeña cantidad de elementos para poder representar
modelos de cualquier dominio físico.
Además, como se verá posteriormente, la notación de Bond-Graph permite visualizar el
comportamiento del sistema más fácilmente que con el conjunto de las ecuaciones de estado.

Pag. T-6 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

1.2.1.- Elementos FUENTE

Los elementos fuente son aquellos que suministran potencia al sistema. Esta potencia puede
ser suministrada a partir de un flujo conocido o mediante un esfuerzo conocido. Así pues,
respectivamente, habrá dos tipos de fuentes:
a) Fuentes de Flujo, con flujo conocido.
b) Fuentes de Esfuerzo, con esfuerzo conocido.

FUENTES DE ESFUERZO.
Son aquellas que provocan una entrada de esfuerzo al sistema. Se representan tal como se
muestra en la figura 1.2.
S proviene de source, fuente
SE
El esfuerzo tiene la barra delante ( regla nemo ED )
Figura 1.2. Fuente de esfuerzo.

Ejemplos físicos de fuentes de esfuerzo son las fuerzas de excitación, pesos, etc., en sistemas
mecánicos, o las pilas, generadores de tensión, fuentes de tensión, etc., en circuitos eléctricos.
Como regla nemo para recordar que las fuentes de esfuerzo son pesos recordar EP, que es muy parecido a las
siglas de España
FUENTES DE FLUJO.
Son aquellas que provocan una entrada de flujo al sistema. Se representan tal como se muestra
en la figura 1.3.

SF S de source y F de flujo

Figura 1.3. Fuente de flujo.

Ejemplos físicos de fuentes de flujo son las fuentes de intensidad en circuitos eléctricos, las
bombas hidráulicas que, realmente, son fuentes de caudal, o la velocidad de excitación de un
sistema mecánico.

El suelo se considera como una fuente de flujo caracterizada por una


velocidad

Pag. T-7 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

1.2.2.- Elementos que almacenan energía.

En este caso consideraremos aquellos elementos que son susceptibles de almacenar energía, y
de devolverla íntegramente al sistema sin ninguna pérdida.
En Bond-Graph, se consideran dos tipos básicos de elementos que almacenan energía: las
Puertas Resorte y las Puertas Inercia. Sus características se describen a continuación.
Lo que se pone a la derecha de los dos puntos es
PUERTA RIGIDEZ porque está caracterizado ese elemento
O puertas resorte
e(t)
Para recordar la nomenclatura de la puerta resorte recordar C :k
que como es un muelle está comprimido luego la barra va f(t)
detrás
Figura 1.4. Puerta Rigidez.

Estas puertas se dibujan en Bond-Graph tal como muestra la figura 1.4., y representan
elementos donde el esfuerzo que aparece se puede expresar como la integral a lo largo del tiempo
del producto del flujo que se produce por un parámetro k(t):
τ
e( t ) = ∫ k ( t )f ( t )dt El esfuerzo es igual a la integral del flujo ec. 1.5.
0

En sistemas dinámicos lineales, el parámetro k suele ser un parámetro constante no


dependiente del tiempo, escribiéndose:
τ Al integrar la velocidad con respecto al tiempo obtengo el desplazamiento
e( t ) = k ∫ f ( t )dt ec. 1.6.
0

Dado que la integral del flujo, es el desplazamiento q(t), se puede escribir:


e( t ) = k ⋅ q ( t ) Recordar que la ecuación de un muelle es que la fuerza esec.
igual1.7.
a la constante
por el desplazamiento
El nombre de este tipo de puertas, se debe a que el comportamiento expresado por la ecuación
1.7., es el que presenta, en el dominio de la Mecánica, la rigidez de un resorte lineal. En este caso,
el esfuerzo es la fuerza “F” aplicada al resorte y el desplazamiento es la deformación total “∆L”
experimentada por el resorte entre sus extremos, (Figura 1.5.). Pero se pueden encontrar otros tipos
de elementos con este tipo de comportamiento en otros dominios de la Física. Por ejemplo, un
condensador en el dominio de la Electricidad, o un acumulador hidráulico.

F F
Un condensador almacena energía
en forma de campo eléctrico y una
bobina en forma de campo
magnético L ∆L

Figura 1.5. Puerta Rigidez.


El desplazamiento es igual a la posición final menos la posición inicial

Pag. T-8 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

PUERTA INERCIA
La notacion de la puerta de inercia es al contrario que la puerta resorte

e(t) La puerta de inercia está caracterizada por la


I :m masa
f(t)

Figura 1.6. Puerta de Inercia.

Estas puertas se dibujan en Bond-Graph tal como muestra la figura 1.6., y representan
elementos donde el esfuerzo que aparece se puede expresar como la derivada temporal del producto
del flujo que se produce por un parámetro m(t):
d[m( t )f ( t )]
e( t ) = ec. 1.8.
dt
En sistemas dinámicos lineales, el parámetro m suele ser un parámetro constante no
dependiente del tiempo, escribiéndose:
d[f ( t )]
e( t ) = m F = m*a ec. 1.9.
dt
El nombre de este tipo de puertas, se debe a que el comportamiento expresado por la ecuación
1.9., es el que presenta, en el dominio de la Mecánica Newtoniana, una masa inercial, (figura 1.7).
En este caso, el esfuerzo es la fuerza “F” aplicada sobre una masa “m” y la derivada del flujo es la
aceleración que la masa experimenta en la dirección y sentido de aplicación de la fuerza.
Se pueden encontrar otros tipos de elementos con este tipo de comportamiento en otros
dominios de la Física. Por ejemplo, una bobina en el dominio de la Electricidad.

f≡ V

e≡F
m

Figura 1.7. Puerta Inercia.

Pag. T-9 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

1.2.3.- Elementos de disipan energía.

PUERTA RESISTENCIA.
Este elemento va a representar a aquellos dispositivos o situaciones donde se produce una
disipación de energía.
Un ejemplo claro es el de la resistencia eléctrica, que al paso de una intensidad de corriente
produce una caída de tensión. Otro ejemplo claro lo constituye el amortiguador hidráulico, que
opone una resistencia, o fuerza opuesta a su movimiento, proporcional a la velocidad a la que se
quiere desplazar. Otros ejemplos son la pérdida de energía por rozamiento mecánico de un sólido
que se desplaza apoyado en otro sólido, o la pérdida de carga que se produce en una tubería debido
al paso de un fluido.
Este elemento se representa en Bond-Graph mediante el grafo mostrado en la figura 1.4., y en
el que el esfuerzo y el flujo están relacionados por la expresión:
e(t) = R f(t) Esto sería la ley de Ohm, E= I * R ec.1.10.
Donde R podría ser un parámetro “p” lineal o no lineal.

La resistencia está caracterizada por un


parámetro p que puede ser lineal o no lineal
La resistencia no lleva barra
R:p
Figura 1.8. Puerta resistencia

Pag. T-10 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

1.2.4.- Nodos de UNIÓN

En aquellos puntos del sistema donde existen una o varias entradas de energía y una o varias
salidas de ésta, se forma un nodo cuya característica principal es que en él se cumple el principio de
conservación de la energía.

e4 f4

e3 e1
f3 • f1
e2
f2
Figura 1.9.- Nodos de Unión.

Como ya es sabido que la energía instantánea o potencia es el producto del esfuerzo “e(t)”,
por el flujo “f(t)” asociado a cada grafo. Debido al principio de conservación de la energía se
cumplirá que la energía que entra al nudo es igual a la energía que sale
n m

∑ e i f i = ∑ e jf j
i =1 j=1
ec. 1.11.

Donde “n” es el número de grafos entrantes al nodo y “m” es el número de grafos salientes
del nodo.
Existen dos tipos de uniones que se denominan tipo 0, y tipo 1. En las primeras todos los
grafos que entran y salen del nodo llevan asociado el mismo valor de esfuerzo. Por el contrario, en
las segundas todos los grafos conectados al nodo llevan asociado el mismo valor del flujo.
Como regla nemo puedo utilizar que el 1 se parece a la f, luego los nudos que tienen un 1 tienen el mismo flujo

NODO DE UNIÓN 0

e4 f4
Todos los grafos que concurren en
el nudo tienen el mismo valor de
esfuerzo e3 e1
f3 0 f1
e2
f2

Figura 1.10.- Nodos de Unión 0.

Pag. T-11 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

Tal como se ha mencionado, en los nodos de unión tipo 0 se cumple:


e1 = e2 = e3 = e4 ec.1.12.
Todos estos esfuerzos son iguales en valor y signo, independientemente del sentido que
presenten los grafos correspondientes respecto de la unión.
Y debido al principio de conservación de la energía expresado en la ecuación 1.6.:
n m
ei ∑ f i = e j ∑ f j
i =1 j=1 Como los esfuerzos en el nudo son iguales, saco factor común y
cancelo, por lo tanto me queda que la suma de flujos entrantes al
n m nudo es igual a la suma de flujos salientes del nudo.

∑ fi = ∑ f j
i =1 j=1
ec. 1.13.

NODO DE UNIÓN 1
Está caracterizado porque todos los grafos
que concurren en el nudo tienen el mismo
valor de flujo. Como regla nemo recordar que e4 f4
el 1 es similar a la f de flujo.

e3 e1
f3 1 f1
e2
f2
Figura 1.11.- Nodo de Unión 1.

Tal como se ha mencionado, en los nodos de unión tipo 0 se cumple:


f1 = f2 = f3 = f4 ec.1.14.
Todos estos flujos son iguales en valor y signo, independientemente del sentido que presenten
los grafos correspondientes respecto de la unión.
Y debido al principio de conservación de la energía expresado en la ecuación 1.6.:
n m
f i ∑ ei = f j ∑ e j
i =1 j=1

n m

∑ ei = ∑ e j
i =1 j=1
ec. 1.15.

Como los flujos en el nudo son iguales, saco factor común y cancelo, por lo tanto me queda
que la suma de esfuerzos entrantes es igual a la suma de esfuerzos salientes.

Pag. T-12 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

1.2.5.- Elementos TRANSFORMER

Dada la palanca de la figura 1.12, vamos a estudiar la relación existente entre las velocidades
y esfuerzos de sus extremos.

FA FB
r VB
A B

VA O
a b

Figura 1.12.- Palanca.

Supongamos que la palanca gira alrededor de su punto de apoyo “O”, con una velocidad
angular “r”. Se tendrá:
VA = r ⋅ a ec. 1.16.
Porque la velocidad angular es igual a la velocidad entre el radio
VB = r ⋅ b ec. 1.17.
Por otra parte, si la palanca gira a velocidad angular constante, hay equilibrio de esfuerzos y
entonces:
Se cumplirá que sumatoria de momentos cero
FA ⋅ a = FB ⋅ b ec. 1.18.
De todas estas ecuaciones se pueden deducir las relaciones existentes entre las velocidades y
esfuerzos de los extremos de la palanca:
VA / a = r = VB / b de donde VB = VA ⋅ (b/a) ec. 1.19.
FA ⋅ a = FB ⋅ b de donde FB = FA ⋅ (a/b) ec. 1.20.
Como se deduce de estas ecuaciones, al pasar del punto A al B de la palanca, se han
transformado tanto velocidades como esfuerzos, pero no ha habido pérdida de potencia y en
consecuencia:
FB VB = FA ⋅ (a/b) VA ⋅ (b/a) = FA VA

Estas transformaciones de esfuerzos y flujos sin pérdida de potencia se representan en Bond-


Graph mediante el elemento TRANSFORMER, de la forma indicada en la figura 1.13.

Pag. T-13 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

p
e1 e2
TF
f1 f2
Figura 1.13.- Elemento TRANSFORMER.

El elemento transformer viene definido por un grafo de entrada y otro de salida, así como por
un valor o razón del transformer “p”. En estos elementos se cumple que el flujo de salida es igual al
flujo de entrada multiplicado por la razón del transformer y, a la inversa, el esfuerzo de salida es el
esfuerzo de entrada dividido por la razón:
Sean:
e1 , f1 = esfuerzo y flujo de entrada.
e2 , f2 = esfuerzo y flujo de salida.
p = b/a = razón del transformer
f2 = p f1 = (b/a) f1 ec. 1.21.
e2 = (1/p) e1 = (a/b) e1 ec. 1.22.
Como no hay pérdidas, siempre se cumple que:
e2 ⋅ f2 = (1/p) ⋅ e1 ⋅ p ⋅ f1 = e1 ⋅ f1 La potencia a la entrada es igual a la de la salida

Partiendo de las expresiones ec. 1.21 y ec. 1.22, conociendo el esfuerzo y flujo de salida,
también se pueden calcular el esfuerzo y flujo de entrada. Con lo cual podríamos invertir el sentido
de actuación del elemento transformer.
En el ejemplo tratado, la razón del transformer “p” se considera constante. Si no es así,
aunque dicha razón “p” varíe con el tiempo, el elemento Transformer “TF” sigue funcionando igual.
Sin embargo, en tal caso se le denomina METATRANSFORMER “MTF”.

Otros ejemplos de aplicación del elemento TRANSFORMER son, en la Mecánica Rotacional


los engranajes rueda-piñón; en la Hidráulica, los actuadores hidráulicos donde una fuerza F aplicada
a una velocidad V se transforma en una presión P y un caudal Q, etc.

Pag. T-14 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

1.2.6.- Elementos GYRATOR

Estos elementos transforman flujos en esfuerzos y esfuerzos en flujos. Están formados por un
grafo de entrada, otro de salida y la razón del gyrator “p”. Como en el caso de los transformer, no
hay pérdida de potencia.

p
e1 e2
GY
f1 f2
Figura 1.14.- Elemento GYRATOR.

Por definición, en los elementos gyrator se cumplirá que:


e2 = p f1 ec. 1.23.
e1 = p f2 ec. 1.24.
Siempre se cumple que:
e2 ⋅ f2 = p ⋅ f1 ⋅ (1/p) ⋅ e1 = e1 ⋅ f1
Y, por tanto, no hay pérdida de potencia

Al igual que en el caso del elemento transformer, si la razón “p” del elemento Gyrator “GY”
no es constante en el tiempo, el funcionamiento sigue siendo el mismo. Pero, en tal caso, se le
denomina elemento METAGYRATOR “MGY”.

Ejemplos de aplicación son, la transformación de las velocidades en fuerzas de inercia en el


dominio de la Mecánica General cuando se utilizan sistemas de referencia no inerciales. Otro
ejemplo muy claro del elemento gyrator, en el dominio de la Electromecánica, es el motor eléctrico,
ya que el paso de una intensidad produce un par motriz, y se cumple que el par es igual a un
parámetro del motor por la intensidad.

Pag. T-15 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

1.3.- La causalidad

Se entiende por hallar la causalidad de los grafos del diagrama de Bond-Graph, el determinar,
de las dos variables que llevan asociadas todos los grafos, cual de ellas es función de la otra. Nunca
ambas serán independientes y siempre una dependerá de la otra.
Por ejemplo, en Mecánica cuando la CAUSA es esfuerzo el EFECTO es una velocidad , a la
inversa, cuando la CAUSA es velocidad, el EFECTO es un esfuerzo. En el primer caso, el esfuerzo
sería la variable independiente y la velocidad la variable dependiente, diciendo que el grafo
asociado al fenómeno tiene CAUSALIDAD ESFUERZO. En el segundo caso, la variable
independiente es la velocidad y el esfuerzo será la variable dependiente, diciendo que el grafo
asociado al fenómeno tiene CAUSALIDAD FLUJO.
En la tabla 1.2. se indican algunas reglas generales, que aplicadas adecuadamente permiten
obtener la CAUSALIDAD DE UN DIAGRAMA BOND-GRAPH de forma sencilla y sistemática.
La representación de la causalidad en Bond-Graph se realiza mediante un trazo transversal en
el extremo del grafo en el sentido de la causalidad ESFUERZO. Este criterio se entenderá mejor
durante el desarrollo de los ejemplos que aparecen en los siguientes apartados.

Tabla 1.2. Reglas de aplicación de la causalidad en Bond-Graph.

Elemento de Bond-Graph CAUSALIDAD

Fuente de flujo FLUJO

Fuente de esfuerzo ESFUERZO

AMBAS
Puerta de Inercia
Causalidad flujo en general

AMBAS
Puerta rigidez
Causalidad esfuerzo en general

Puerta resistencia AMBAS

Transformer AMBAS

Gyrator AMBAS

En cuanto exista un grafo con causalidad esfuerzo entrante


Unión 0 al nudo, el resto de grafos tendrá causalidad esfuerzo en
sentido saliente del nudo.
Estos son en los que todos tienen el mismo esfuerzo

En cuanto exista un grafo con causalidad flujo entrante al


Unión 1 nudo, el resto de grafos tendrá causalidad flujo en sentido
Estos son aquellos en los que todos tienen la misma flujo
saliente del nudo.

Pag. T-16 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

1.4.- Planteamiento de las ecuaciones del sistema

1.4.1.- Introducción

Dado un sistema físico cualquiera, el planteamiento de sus ecuaciones de estado puede


hacerse de diferentes formas.
Un sistema de N grados de libertad puede ser representado por:
a) Una ecuación de orden 2N en términos de una variable desconocida.
b) Varias combinaciones de ecuaciones y variables desconocidas con su orden
correspondiente, (no necesariamente el mismo).
c) 2N ecuaciones acopladas de primer orden en términos de N variables desconocidas.
Estas ecuaciones suelen ser tipo diferencial.

Por ejemplo, para el sistema de 2 grados de libertad de la figura 1.15., en el que K1 , K2 son
las rigideces de los resortes, m1 , m2 las masas y x1 , x2 son los desplazamientos verticales de las
citadas masas, respectivamente.

K1 x1

m1 x2

x3
K2

m2

Figura 1.15. Sistema de 2 grados de libertad.

Tomando como variable x1, para expresar el sistema en términos de este desplazamiento, la
ecuación de estado tiene la forma:

d 4 x1  K1  1 1  d 2 x 1  K 1 ⋅ K 2   1 1 
− + + K 2  +  ⋅ 2 +   ⋅ x 1 − K 2  +  ⋅ g = 0
 m1 ⋅ m 2
4
dt  m1  m1 m 2  dt   m1 m 2 
Siendo g la aceleración de la gravedad.
Como se observa, en este caso se obtiene una sola ecuación de cuarto orden.

Pag. T-17 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

Una segunda forma de plantear las ecuaciones, consiste en elegir como variables
desconocidas x1 y x3 como también se muestra en la figura 1.15.. En tal caso, el modelo viene
representado por dos ecuaciones de segundo grado acopladas:
d2x3  1 1  K
2
+ K 2  +  x 3 − 1 x 1 = 0
dt  m1 m 2  m1
d 2 x1 K 2 K
2
− x 3 + 1 x1 = 0
dt m1 m1

Finalmente, si se asocia la energía del sistema a la cantidad de movimiento de las masas y al


desplazamiento producido entre los extremos de los resortes, entonces se tiene un sistema de cuatro
ecuaciones de primer orden acopladas.
Denominando P1 y P2, las cantidades de movimiento asociadas a las masas m1 y m2 ,
respectivamente, el sistema de ecuaciones que se obtiene es el siguiente:
dP1
= K 2 ⋅ x 3 − K 1 ⋅ x 1 + m1 ⋅ g
dt
dP2
= −K 2 ⋅ x 3 + m 2 ⋅ g
dt
dx 1 P1
=
dt m1
dx 3 P P
= 2 − 1
dt m 2 m1
En todo caso, es posible transformar una cualquiera de las formas en otra, tanto si el sistema
es lineal, como si no lo es. Sin embargo, para sistemas no lineales la transformación es más difícil.
La técnica de Bond-Graph, permite aplicar esta última formulación de una forma muy sencilla
y sistemática, si se seleccionan adecuadamente las variables del sistema. Este será el caso de
formulación que se tratará a continuación en los siguientes apartados. Como se verá, se tomarán
como variables las cantidades de movimiento de las puertas de inercia y los desplazamientos en las
puertas rigidez.

Pag. T-18 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

1.4.2.- Desarrollo de las ecuaciones de estado del sistema en el diagrama de Bond-Graph.

Para tratar el tema de obtención de las ecuaciones de estado de un sistema modelizado en Bond-
Graph vamos a utilizar el ejemplo que aparece en la figura 1.16.. En esta figura se muestran dos masas
suspendidas, en el caso de la superior mediante un muelle-amortiguador en paralelo, y en la inferior
solamente con un muelle.
Se expondrá la obtención del diagrama de Bond-Graph del sistema, que se muestra en la figura
1.17, considerando incluso la causalidad.
A partir del diagrama del Bond-Graph obtenido, se irán obteniendo las ecuaciones de estado del
sistema, que son las expresiones de los flujos (velocidades) y de los esfuerzos (fuerzas) que se
desarrollan en los diferentes puntos o elementos que constituyen el sistema.
En el apartado siguiente se obtendrán las ecuaciones de comportamiento del sistema, que son
un sistema de ecuaciones diferenciales respecto del tiempo de las variables independientes del sistema,
que integradas permiten calcular o simular el comportamiento del sistema partiendo de determinadas
condiciones iniciales.

V2
m2

K2 C2
V1
m1

K1

Figura 1.16. Modelo de dos masas suspendidas.

8 E 7
12 11
Se : m 2g 1 I:m
2
6
10 8
D C :K
2
C 0
7
1
9 5

6 R :C 2
3 B 4
5 4
Se : m 1 g 1 I:m
1

3
2
2
A 0 C :K
1

1 Sf : 0

Figura 1.17. Bond Graph del modelo de dos masas suspendidas, con causalidad.

Pag. T-19 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

Como se observa en la figura 1.16., y según ya se ha citado, el sistema está constituido por dos
masas suspendidas, en el caso de la superior mediante un muelle-amortiguador en paralelo, y en la
inferior solamente con un muelle, que está apoyado en un suelo. Se quiere estudiar el movimiento de
las masas en sentido vertical solamente.
El desarrollo que se expone a continuación para explicar la obtención del diagrama de Bond-
Graph, así como las ecuaciones de estado del sistema, se refiere continuamente a las figuras 1.16. y
1.17, aunque no se haga referencia explícita a ellas. Como se observa en estas figuras, todos los
elementos están numerados o nombrados, para su mejor seguimiento.
La velocidad del extremo inferior del resorte inferior “K1” vendrá definida por la velocidad del
suelo. Este aspecto se modeliza mediante una fuente de flujo, SF : 0 , (1). Las fuentes de flujo siempre
tienen causalidad flujo hacia el sistema. En este caso el flujo f1 es la variable independiente, causa
conocida:
f1 = V0 = velocidad vertical del suelo ec. 1.25.
Los resortes K1 y K2 se representan mediante puertas Compliance o Rigidez, (2) y (6),
respectivamente. Estos resortes están trabajando sometidos en sus extremos a la diferencia de
velocidad que imponen el suelo y la masa m1, en el caso de K1; y la diferencia de velocidad entre las
masas m1 y m2, en el caso de K2. Por tanto, la causalidad desde el sistema a estas puertas es flujo. Y
como estos resortes responden al sistema con un esfuerzo, su causalidad hacia el sistema es de tipo
esfuerzo. Efectivamente, en la figura 1.18, puede observarse como el sistema informa a la puerta 2 de
la diferencia de velocidades a la que trabaja el resorte, y como éste devuelve al sistema el esfuerzo que
se genera en él debido al desplazamiento entre sus extremos.

Causalidad
flu jo, V1 (t)
4
4
1 I:m
1

Causalidad Causalidad flujo, V0-V1 (t)


flujo, V1 (t) 3
2
2
A 0 C :K
1
Causalidad
flujo, V0 (t) 1 Causalidad esfuerzo, F 1(t)

1 Sf : 0

Figura 1.18. Causalidad hacia la puerta 2 y de ella hacia el sistema.

Pag. T-20 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

Tomando como variable independiente la deformación relativa en los extremos de los resortes,
(que será la integral temporal de la velocidad relativa entre los extremos), las ecuaciones de
comportamiento de éstos son:
e2 = esfuerzo de la puerta 2 = K1 x1 ec. 1.26.
e8 = esfuerzo de la puerta 6 = K2 x2 ec. 1.27.
Siendo x1 y x2, las deformaciones de los resortes 1 y 2, respectivamente.
Es de destacar el hecho de que sobre estas puertas entran flujos y esta devuelve esfuerzo. Todas
las puertas de Bond Graph reciben un tipo de causalidad, (flujo o esfuerzo), y devuelven la otra al
sistema.

El amortiguador C2 se representa mediante una puerta Resistencia, (5). Este amortiguador está
trabajando sometido en sus extremos a la diferencia de velocidad que existe entre las masas m1 y m2.
Por tanto, la causalidad desde el sistema a esta puerta es flujo. Y en contra, este amortiguador responde
al sistema con causalidad de tipo esfuerzo, definiendo el esfuerzo:
e9 = esfuerzo de la puerta 5 = C2 f9 ec. 1.28.

Las masas m1 y m2 se representan mediante puertas Inercia, (4) y (6), respectivamente. Ambas
masas están sometidas a los esfuerzos de los resortes, amortiguador y los pesos correspondientes. Por
tanto, la causalidad del sistema sobre las puertas de Inercia es de tipo esfuerzo. Por otra parte, las masas
responderán al sistema con las velocidades correspondientes, presentando causalidad flujo hacia el
sistema. (Ver la figura 1.19).
Tomando como variable independiente el momento Pi = mi vi , o Cantidad de Movimiento en las
masas, el flujo que definen estas hacia el sistema es:
f4 = flujo en la puerta de inercia 4 = P1 / m1 ec. 1.29.
f11 = flujo en la puerta de inercia 7 = P2 / m2 ec. 1.30.

Las masas están sometidas a su peso propio, lo cual constituye una entrada de esfuerzo externo
al sistema. Este hecho se representa mediante sendas Fuentes de Esfuerzo. Las fuentes de esfuerzo
poseen causalidad ESFUERZO hacia el sistema, siendo los esfuerzos e5 y e8 los valores del peso de las
masas, causas conocidas:
e5 = peso de la masa 1 = m1 g ec. 1.31.
e12 = peso de la masa 2 = m2 g ec. 1.32.

Pag. T-21 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

Causalidad esfuerzo F2 (t )

Causalidad esfuerzo F2 (t )+Famo rt 6


8
7
D C : K2
C 0 1
5
9
Causalidad
R : C2
esfuerzo
F2 (t )+Famo rt Causalidad esfuerzo Famo rt.

6
Causalidad esfuerzo
3 B 4 Σ F(t )
5 4
Se : m 1 g 1 I : m1

Causalidad esfuerzo m 1 .g 3 Causalidad flujo v 1 (t )

Causalidad
esfuerzo F 1 (t )

2
2
A 0 C : K1

Causalidad esfuerzo F1 (t )

Figura 1.19. Causalidad hacia la puerta 4 y de ella al sistema.

Una vez descritas las diferentes puertas que forman parte del diagrama de Bond-Graph,
seguidamente vamos a analizar las conexiones entre ellas.
La conexión del resorte K1 con el resto del sistema se hace mediante la Unión 0, (A), pues esta
unión se produce a esfuerzo constante definido por el resorte (2). (Ver figura 1.20):
e1 = e3 = e2 ec. 1.33.
Sobre esta unión entra el grafo de la fuente de flujo y sale, además del de la Compliance, el que
se dirige hacia la unión 1 de la Inercia (4). Al ser de tipo 0 se plantea que la suma de flujos de entrada
es igual a la suma de flujos de salida:
f1 = f2 + f3 ec. 1.34.
El flujo del grafo inferior viene definido por la fuente de flujo 1, y en cuanto se refiere al
superior proviene de una unión tipo 1 que como es conocido significa igualdad de flujos. El flujo de
todos los grafos que salen de esta unión vendrá definido por el del grafo que entra en la unión con flujo
conocido. (Ver figura 1.19). En este caso entra en la unión con flujo conocido el procedente de la
puerta 4. El flujo que se dirige hacia la puerta 2 será el del grafo inferior, menos el del superior que es
igual al flujo del grafo 4 asociado a la puerta de Inercia (4):

Pag. T-22 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

f2 = f1 - f3 ==> f2 = f1 - f4 = Puerta 1 - Puerta 4


En la figura 1.21, se muestra en esquema un resumen de la relación entre los esfuerzos y los
flujos en la Unión A de tipo 0.

Entra con flujo conocido en


la unión 0. El flujo de este
grafo es el de la puerta 4.
Sale de la unión con flujo
3
conocido
2
2
A 0 C :K
1

1
Entra con flujo conocido en
1 Sf : 0
la unión 0. El flujo de este
grafo es el de la puerta 1.
Figura 1.20. Planteamiento de causalidad sobre la unión A.

6
3 4
5 4
Se : m 1 g 1 I:m
1

3
2
2
A 0 C :K
1

En el nudo A ==> f2 + f3 = f1
1
f1 = Flujo del grafo 1
1 Sf : 0
f2 = Flujo del grafo 2
f3 = f4 = Flujo de definido por la puerta 4
f2 = f1 - f3

e1 = e3 = e2 = Esfuerzo definido
por la puerta 2

Nota : fi hace referencia al flujo del grafo “i”, no de la puerta “i”.


Figura 1.21. Esquema de flujos y esfuerzos en la unión A.

Continuando en sentido hacia arriba sobre el sistema, (figuras 1.16 y 1.17), el siguiente elemento
que nos encontramos es la masa m1. La conexión de la masa m1 con el resto del sistema se hace
mediante la Unión 1, (B), pues esta unión se produce a flujo constante definido por masa:
f3 = f5 = f6 = f4 = Flujo definido por la masa m1 ec. 1.35.
Sobre esta unión entra el grafo de la fuente de esfuerzo de la acción del peso y el grafo de la
acción del resorte K1, y salen el grafo de la Inercia y el que se dirige hacia el resorte y amortiguador

Pag. T-23 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

superiores. La Unión 1 (B), al ser de tipo 1 se plantea que la suma de esfuerzos de entrada es igual a la
suma de esfuerzos de salida de salida:
e4 + e6 = e3 + e5 ec. 1.36.
En cuanto se refiere a la puerta 4, la causalidad de entrada en esta puerta es “esfuerzo conocido”
y proviene de una unión tipo 1, en donde se plantea que la suma de esfuerzos de entrada es igual a la
suma de esfuerzos de salida. Como puede observarse en la figura 1.20, los grafos 3, 5 y 6, entran en la
unión con esfuerzo conocido, mientras que el grafo 4, sale de ella con esfuerzo conocido. Por tanto:
e4 = e3 + e5 - e6 ec. 1.37.
El grafo 3, sale de la unión A tipo 0 analizada anteriormente. En las uniones 0, todos los grafos
menos uno salen de la unión con esfuerzo conocido, y el que queda es el que entra con causalidad
esfuerzo y por lo tanto define el esfuerzo de la unión 0.
En la unión A, el grafo que define el esfuerzo de esta unión es el que entra de la puerta 2, y por
lo tanto el esfuerzo del grafo 3 viene definido por la puerta 2.
El esfuerzo del grafo 5, viene de la puerta 3, fuente de esfuerzo, el cual se define directamente
como el peso de la masa m1. Y, por último el del grafo 6 viene de la unión C superior que, al tratarse de
una unión 0, el planteamiento es idéntico al realizado anteriormente con la unión A. Se comienza
averiguando el grafo que entra en la unión con causalidad esfuerzo, que en este caso es el grafo 7 y
proviene de una unión C tipo 1 cuyo esfuerzo se define por la suma de los de las puertas 5 y 6.
En la figura 1.22, se resume la relación de esfuerzos y los flujos en la Unión B de tipo 1.

6
8
D C :K
2
C 0
7
1
9 5

6 R :C 2
3 B 4
5 4
Se : m 1 g 1 I:m
1

3 2
2
A 0 C :K
1

f3 = f5 = f6 = f4 = Flujo definido por la puerta 4

En el nudo B se cumple: e4 + e6 = e3 + e5 (1)


En el nudo A: e3 = e2
En el nudo C: e6 = e7
En el nudo D: e7 = e8+e9
Luego, de (1) : e4 = e3+e5-e6 = e2 + e5 - e8 - e9
Además:
e2 = Esfuerzo puerta 2
e5 = Esfuerzo puerta 3
e8 = Esfuerzo puerta 6
e9 = Esfuerzo puerta 5
Por tanto:
Figura 1.22. Esquema de flujos y esfuerzos en la unión B tipo 1 de la puerta de Inercia 4.

A continuación estudiamos las uniones C y D.

Pag. T-24 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

La unión D es de tipo 1, (flujo constante), pues estamos conectando dos elementos mecánicos,
la puerta Resistencia nº 5 y la puerta Rigidez nº 6, que trabajan en paralelo sometidos a la misma
velocidad relativa entre sus extremos, que es la velocidad relativa entre las masas m1 y m2. Y, por
tanto, con causalidad de entrada a ellas “Flujo conocido”. Este flujo es el que proviene de la unión D
de tipo 1, y el grafo que entra en esta unión con flujo conocido proviene de la unión C de tipo 0, pues
es necesario plantear la ecuación de suma de flujos definidos por las puertas de Inercia nº 4 y nº 7. En
la figura 1.23, se van indicando los pasos a seguir para determinar los flujos de los grafos que entran en
la puerta 5 y en la puerta 6.
En dicha figura 1.23 también se muestra el cálculo de esfuerzos en los grafos que confluyen en
la unión D. Como se observa, puesto que el sistema actúa sobre las puertas 5 y 6 con causalidad flujo,
estas puertas reaccionan contra el sistema con causalidad esfuerzo. Estos esfuerzos definidos por
ambas puertas llegan, a través de los grafos 8 y 9, a la unión D de tipo 1. En esta unión se plantea la
suma de esfuerzos, quedando definido el esfuerzo del grafo 7:
e7 = e8 + e9 = Puerta 5 + Puerta 6 ec.1.38.
Dicho grafo 7 entra en la unión C de tipo 0, con esfuerzo conocido y, por tanto, queda definido
el esfuerzo de esta unión C. Este esfuerzo de la unión C, definido por el grafo 7, como acabamos de
ver, es la suma de los esfuerzos en las puertas Resistencia nº 5 y Rigidez nº 6, que es la fuerza que se
transmite a las masas m1 y m2 como reacción del resorte y amortiguador a la acción del movimiento
relativo entre ambas masas.

E 7
11
1 I:m
2
6
10 8
D C :K
2
C 0
7
1
9 5

6 R :C 2
B 4
4
1 I:m
1

En el nudo D, por ser de tipo 1, se cumple: f8 = f9 = f7


e7 = e8 + e 9
e8 = Esfuerzo definido por la puerta 6
e9 = Esfuerzo definido por la puerta 5
En el nudo C, de tipo 0, se tiene: e6 = e7 = e10
f7 + f10 = f6
Por el nudo B, de tipo 1: f6 = f4
Por el nudo E, tipo 1: f10 = f11
Además, en las puertas:
f4 = flujo de la puerta 4
f11 = flujo de la puerta 7
Luego el flujo entrante a las puertas 5 y 6, queda:
f8 = f9 = f6 - f10 = Puerta 4 - Puerta 7
Figura 1.23. Esquema de flujos y esfuerzos en las uniones C y D, y flujos que se dirigen a
las puertas 5 y 6.

Pag. T-25 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

Por último, queda la unión con sistema de la masa m2 , (puerta de Inercia 7), a través de la unión
E de tipo 1, pues se hace con flujo definido por la masa hacia el sistema. Por tanto, la causalidad de
entrada hacia ella es “esfuerzo conocido”. Este esfuerzo proviene de la unión E de tipo 1, en donde
será necesario plantear la suma de esfuerzos. Sobre esta unión entran dos grafos con causalidad
esfuerzo, el que viene de la puerta 8, fuente de esfuerzo definido por el peso de la masa m2, y el que
procede de la unión C inferior de tipo 0, que como ya hemos estudiado, es el esfuerzo del grafo 7, que
también es la suma de los esfuerzos de las puertas resistencia nº 5 y Rigidez nº 6.
En la figura 1.24, se plantean los flujos y esfuerzos en los grafos asociados a la unión E, hasta
llegar a la puerta nº 7.

8 E 7
12 11
Se : m2g 1 I:m
2
6
10 8
D C :K
2
C 0
7
1
9 5

R :C 2

En el nudo E, se cumple: f10 = f12 = f11


e11 = e10 + e12
Por el nudo C: e10 = e7
Y, por el nudo D: e7 = e8 + e9
Además, se tiene:
e12 = Esfuerzo puerta 8
e8 = Esfuerzo puerta 6
e9 = Esfuerzo puerta 5
Por tanto, el esfuerzo entrante a la puerta 7, es:
e11 = Puerta 8 + Puerta 5 + Puerta 6

Figura 1.24. Esfuerzos y flujos en la unión E, y esfuerzos que se dirigen a la puerta 7.

Pag. T-26 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

Como conclusión de este apartado, se muestra el listado con las ecuaciones de estado del sistema
ya obtenidas anteriormente:
Valores de los flujos
f1 = V0 Por ser fuente de flujo conocida.
f2 = f1 - f3 Por ser grafos de la unión A tipo 0
f3 = f5 = f6 = f4 Por ser grafos de la unión B tipo 1.
f4 = P1 / m1 Por ser puerta de inercia, tomando P1 como variable independiente.
f7 = f6 - f10 Por ser grafos de la unión C tipo 0.
f8 = f9 = f7 Por ser grafos de la unión D tipo 1.
f10 = f12 = f11 Por ser grafos de la unión E tipo 1.
f11 = P2 / m2 Por ser puerta de inercia, tomando P2 como variable independiente.

Valores de los esfuerzos


e1 = e3 = e2 Por ser grafos de la unión A tipo 0.
e2 = K1 x1 Por ser puerta Rigidez, con variable independiente la deformación x1.
e4 = e3 + e5 - e6 Por ser grafos de la unión B tipo 1.
e5 = -m1 g Por ser fuente de esfuerzo.
e6 = e7 = e10 Por ser grafos de la unión C tipo 0.
e7 = e8 + e9 Por ser grafos de la unión D tipo 1.
e8 = K2 x2 Por ser puerta Rigidez, con variable independiente la deformación x2.
e9 = C2 f9 Por ser puerta resistencia, con flujo f9 conocido.
e11 = e10 + e12 Por ser grafos de la unión E tipo 1.
e12 = -m2 g Por ser fuente de esfuerzo.

Pag. T-27 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

1.4.3.- Ecuaciones de comportamiento del sistema a partir del diagrama de Bond-Graph.

Partiendo de las ecuaciones de estado obtenidas en el apartado anterior, en este apartado


vamos a tratar la manera de obtener las ecuaciones de comportamiento del sistema. Ya se ha
mencionado en el apartado anterior que las ecuaciones de comportamiento del sistema son un
sistema de ecuaciones diferenciales respecto del tiempo de las variables independientes del sistema,
que integradas permiten calcular o simular el comportamiento del sistema partiendo de determinadas
condiciones iniciales.
En Bond-Graph, se toman como variables independientes del sistema:
A) Los momentos (P = m f), en las puertas de Inercia, (que significa la Cantidad de Movimiento
de las masas en el dominio de la Mecánica).
B) La deformación x en las puertas Rigidez, (que es la deformación de los resortes en el dominio
de la Mecánica).
Por tanto, ahora tomaremos las ecuaciones de estado obtenidas en el apartado anterior, en las
cuales aparecen las variables independientes mencionadas:
f4 = P1 / m1
f11 = P2 / m2
e2 = K1 x1
e8 = K2 x2 ecs. 1.39.

El significado físico de los flujos y esfuerzos que aparecen en estas ecuaciones es el siguiente:
f4 = es la velocidad vertical V1 de la masa m1
f11 = es la velocidad vertical V2 de la masa m2
e2 = es la fuerza F1 desarrollada por el resorte de rigidez K1
e8 = es la fuerza F2 desarrollada por el resorte de rigidez K2

Tomamos las expresiones ecs. 1.39, y despejamos las variables independientes consideradas:
P1 = m1 f4 = m1 V1
P2 = m2 f11 = m2 V2
x1 = e2 / K1 = F1 / K1
x2 = e8 / K2 = F2 / K2 ecs. 1.40.

Y, ahora, aplicamos derivadas respecto del tiempo, obteniendo:


dP1 dV
= m1 1 = m1 ⋅ a 1
dt dt
dP2 dV2
= m2 = m2 ⋅ a2
dt dt

Pag. T-28 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

dx 1 1 dF1
=
dt K 1 dt
dx 2 1 dF2
= ecs. 1.41.
dt K 2 dt

Analicemos el sentido físico de los términos de las ecuaciones ecs. 1.41. que acabamos de
obtener:
a) La primera derivada del momento o cantidad de movimiento de las masas del sistema,
dPi/dt, es igual al producto de la masa por la aceleración correspondiente, ai, que, por la
Segunda Ley de Newton, (F=ma), es la fuerza o suma de fuerzas que llegan a la masa.
Mirando en el diagrama de Bond-Graph obtenido en el apartado anterior, (figura 1.17), la
fuerza o suma de fuerzas que llegan a cada masa es el esfuerzo del grafo asociado a la puerta
de Inercia que representa el comportamiento de cada una de las masas, es decir:
e4 en el caso de la masa m1.
e11 en el caso de la masa m2.
b) La derivada de los desplazamientos o deformaciones de los resortes, dxi/dt, es la velocidad
relativa entre los extremos de los mismos.
Mirando en el diagrama de Bond-Graph obtenido en el apartado anterior, (figura 1.17), la
velocidad relativa entre los extremos de los resortes es el flujo de los grafos asociado a las
puertas Rigidez que representan el comportamiento de cada una de los resortes, es decir:
f2 en el caso del resorte de rigidez K1.
f8 en el caso del resorte de rigidez K2.

Por tanto, las ecuaciones 1.41, quedan:


dP1
= e4
dt
dP2
= e11
dt
dx 1
= f2
dt
dx 2
= f8 ecs. 1.42.
dt

Los flujos y esfuerzos, e4 e11 f2 f8 , los tenemos calculados en la ecuaciones de estado en el


apartado anterior:
e4 = e3 + e5 - e6 Por ser grafos de la unión B tipo 1.
e3 = e2 Por ser grafos de la unión A tipo 0.
e2 = K1 x1 ecs. 1.39.
e5 = -m1 g Por ser fuente de esfuerzo.
e6 = e7 Por ser grafos de la unión C tipo 0.

Pag. T-29 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

e7 = e8 + e9 Por ser grafos de la unión D tipo 1.


e8 = K2 x2 ecs. 1.39.
e9 = C2 f9 Por ser puerta resistencia, con flujo f9 conocido.
f9 = f8 , por ser grafos de la unión D tipo 1.
(f8 se calcula más adelante).

e11 = e10 + e12 Por ser grafos de la unión E tipo 1.


e10 = e7 Por ser grafos de la unión C tipo 0.
e7 = e8 + e9 Por ser grafos de la unión D tipo 1.
e8 = K2 x2 ecs. 1.39.
e9 = C2 f9 Por ser puerta resistencia, con flujo f9 conocido.
f9 = f8 , por ser grafos de la unión D tipo 1.
(f8 se calcula más adelante).
e12 = -m2 g Por ser fuente de esfuerzo.

f2 = f1 - f3 Por ser grafos de la unión A tipo 0


f1 = V0 Por ser fuente de flujo conocida.
f3 = f4 Por ser grafos de la unión B tipo 1.
f4 = P1 / m1 ecs. 1.39.

f8 = f9 = f7 Por ser grafos de la unión D tipo 1.


f7 = f6 - f10 Por ser grafos de la unión C tipo 0.
f6 = f4 Por ser grafos de la unión B tipo 1.
f4 = P1 / m1 ecs. 1.39.
f10 = f11 Por ser grafos de la unión E tipo 1.
f11 = P2 / m2 ecs. 1.39.

Sustituyendo sucesivamente los valores de las distintas expresiones anteriores, se obtiene:


e4 = K1 x1 - m1 g – (K2 x2 + C2 f8)
e11 = (K2 x2 + C2 f8) - m2 g
f2 = V0 - P1 / m1
f8 = P1 / m1 - P2 / m2 ecs. 1.43.

Sustituyendo, nuevamente, estas últimas ecuaciones en las ecuaciones ecs. 1.42., queda:
dP1   P P 
= K 1 x 1 - m1 g - K 2 x 2 + C 2  1 - 2 
dt   m 1 m 2 
dP2   P P 
= K 2 x 2 + C 2  1 - 2  - m 2 g
dt   m1 m 2 
dx 1 P
= V0 − 1
dt m1
dx 2 P P
= 1 − 2 ecs. 1.44.
dt m1 m 2

Pag. T-30 de 31
J.A. Lozano – ETSIDI - UPM Introducción a técnica de Bond-Graph

Las ecuaciones ecs. 1.44, son las ecuaciones de comportamiento del sistema de dos grados de
libertad mostrado en la figura 1.16, que en forma matricial, finalmente quedan:
dP1
dt − C 2 m1 C2 m2 K1 - K 2 P1 - m1 g
dP2
dt = C 2 m1 − C2 m2 0 K 2 P2 - m2 g
+ ec. 1.45.
dx 1 - 1 m1 0 0 0 x1 V0
dt 1 m1 -1 m 2 0 0 x2 0
dx 2
dt

Partiendo de unas condiciones iniciales dadas, aplicando algoritmos numéricos por ordenador,
resulta muy sencillo resolver o integrar en el tiempo este sistema de ecuaciones diferenciales. De esta
manera se puede simular el comportamiento del sistema y obtener su respuesta temporal en función de
unas condiciones iniciales.

Pag. T-31 de 31

También podría gustarte