Transcript
Control Engineering Laboratory
PORT-BASED MODELING OF DYNAMIC SYSTEMS: fundamental concepts and bond graphs Peter Breedveld Control Engineering Laboratory, Faculty of Electrical Engineering, Mathematics and Computer Science University of Twente, Netherlands
[email protected] 7/18/2005
Summerschool Bertinoro 17-23 July 2005
© P.C. Breedveld 1
Some remarks • Role of concepts (‘intellect’, multiple views) • Ports from a historical, energy point of view (energy states) • Focus on linear motion and fixed axis rotation, simple configuration constraints: no need for abstract geometrical concepts (like in the rest of the course) • Separation between configuration state and energy state: structures influence of kinematic constraints
Summerschool Bertinoro, 17-23 July 2005 (2)
© Peter Breedveld
General Outline • Background and history • Introduction to ports, bonds and physical domains • Introduction to bond graph modeling: causality, domain-independence, etc.
Summerschool Bertinoro, 17-23 July 2005 (3)
© Peter Breedveld
Assumed prerequisites • Some experience in: – – – –
Block diagram models Numerical simulation (methods) Linear and nonlinear analysis Basic physics (electrical circuits, simple rigid body mechanisms, principles of thermodynamics, etc.
Summerschool Bertinoro, 17-23 July 2005 (4)
© Peter Breedveld
Introduction – Modular modeling (prevailing trend in modeling and simulation of complex physical systems):
• represented as a network interconnection of ideal elements
– Advantages
• flexibility • re-usability of model parts • physical insight • support for automated modeling • fundamental to systems and control theory
Summerschool Bertinoro, 17-23 July 2005 (5)
© Peter Breedveld
Introduction – Disadvantages
• equations of motion obtained directly from network modeling are often – complicated – without apparent structure – will easily contain algebraic constraints arising from the interconnection of the sub-systems
• no problem for digital simulation, but • may not be very well-suited to analysis and control purposes, (nonlinear systems!)
Summerschool Bertinoro, 17-23 July 2005 (6)
© Peter Breedveld
Course objectives – particular type of network modeling, viz. port-based modeling:
• sub-systems are interacting through power exchange represented by pairs of conjugated variables • unified way to deal with systems from different physical domains (e.g., mechanical and electrical) • resulting systems of equations possess an underlying generalized Hamiltonian structure, leading to the geometric notion of a port-Hamiltonian system
Summerschool Bertinoro, 17-23 July 2005 (7)
© Peter Breedveld
Course objectives • exploiting the underlying Hamiltonian structure offers many possibilities for – – –
analysis simulation control
of complex physical systems
• illustrated in variety of application areas • possible extension to classes of distributed-parameter systems
Summerschool Bertinoro, 17-23 July 2005 (8)
© Peter Breedveld
Appetizer: Newton cradle example • Experiment
• 20-sim demo • Animation
Summerschool Bertinoro, 17-23 July 2005 (9)
© Peter Breedveld
Dual role example: Falling and bouncing object • generally considered as conceptually complex:
– switching models reinitialization, timing – energy bookkeeping – prone to sign errors
• port-based approach with explicit structure, allows conceptually simple solution Summerschool Bertinoro, 17-23 July 2005 (10)
© Peter Breedveld
Aim • System theoretic approach to physical system dynamics based on – classification of phenomena in terms of energy – fundamental principles of thermodynamics
• shown how, as a result: – variables and relations describing physical systems may be classified – models may be organized as (‘port-based approach’) • multiport elements • interconnected in an interconnection structure corresponding to a generalized network
– multiport elements describe basic behaviors with respect to energy and entropy
• special attention: – role of analogies – analogue behavior Summerschool Bertinoro, 17-23 July 2005 (11)
© Peter Breedveld
Fundamental issues • introduction & some philosophy: • Why physical systems modeling? • What is physical systems modeling? • Context of explanation (focus of engineering physics) versus justification (focus of math) • elements and components, theory building • synthesis between classical approaches
• choice of variables: – mechanical versus thermodynamic framework Summerschool Bertinoro, 17-23 July 2005 (12)
© Peter Breedveld
‘Why?:’ Three-world meta-model • Real world (assuming existence of ‘objective’ environment’) • Conceptual world (in our brain) • ‘Paper’ world, including electronically stored data – Only in/via the paper world: • • •
Communication Support Systematization
– exchangeable abstractions/concepts >>> – importance of symbols & notation (representation) Summerschool Bertinoro, 17-23 July 2005 (13)
© Peter Breedveld
‘Why?’: Are physical concepts ‘real’ ? • For example:
– energy, time, momentum, causality,…
• In a context of justification: NO! • In a context of discovery/explanation: YES, but rather ‘really useful’ than just ‘real’
Summerschool Bertinoro, 17-23 July 2005 (14)
© Peter Breedveld
‘Why?:’ The useless ‘quest for truth’ • A model is necessarily incomplete: ‘all models are wrong’, but ‘truth’ is not what counts • The issue is whether a model is competent (to solve a problem in a given problem context in the most generic sense of the word)
Summerschool Bertinoro, 17-23 July 2005 (15)
© Peter Breedveld
What is a model? Some sort of abstraction (in the ‘paper’ world) that enables • • • • • • • •
insight in the real world counterpart communication about the real world counterpart observation of the real world counterpart troubleshooting of the real world counterpart design of new aspects related to the real world counterpart modifications of the real world counterpart ‘explanation’ of functionality of the real world counterpart measurement of the real world counterpart
but, most importantly, that is :
• competent to solve a given problem and make decisions related to the real world
Summerschool Bertinoro, 17-23 July 2005 (16)
© Peter Breedveld
Modeling • Given a specific problem context, the decision process to obtain a competent model to solve this problem • Approaches between two extremes: – a priori knowledge – ‘black box’ (also an axiomatic concept with a priori
assumptions!!!) • concept ‘input’ implicitly contains model of imposing some action with negligible back effect (‘high input impedance’) • concept ‘output’ implicitly contains model of measurement with negligible effect on the system being observed (‘low output impedance’) • after a competent input-output relation is found, it is not open for modifications or physical interpretation
Summerschool Bertinoro, 17-23 July 2005 (17)
© Peter Breedveld
‘Special’ states Position/displacement • has a dual nature:
– energy state (related to a conservation or symmetry principle like all other states) – configuration state
• does not transform as a tensor Matter • convects all matter-bound properties (not volume!) • conjugate intensity depends on other intensities • boundary criterion Volume • boundary criterion Entropy • can be ‘locally’ produced (only ‘locally’ conserved)
Summerschool Bertinoro, 17-23 July 2005 (18)
© Peter Breedveld
Model representation • model representation:
– symbols used to represent • the concepts being used • (the structure of) their relations (interconnection)
• model manipulation (as opposed to modeling!) :
– transformation to different representations (including the ‘solution’) to • increase insight • draw conclusions, etc.
Summerschool Bertinoro, 17-23 July 2005 (19)
© Peter Breedveld
Role of representations Sequential:
x1
a
T
A
x2
b c
T x3
d
A T x4
A
AT x5
e
AT t
Simultaneous: T
Process time versus processing time Summerschool Bertinoro, 17-23 July 2005 (20)
A
A
x1
T
x2
x
x3
x5 x4
© Peter Breedveld
What are we trying to describe? • (dynamic) behavior • engineering systems
• herein constrained to (for sake of ‘simplicity’):
– deterministic models of systems that obey basic principles of macroscopic physics: • energy conservation (‘first law’) • positive entropy production (‘second law’)
• that describe the temporal trajectory of the common physical properties:
– mechanical (incl. hydraulic and pneumatic), electrical, magnetic, chemical, material, thermal, etc. >>> multidomain – configuration
Summerschool Bertinoro, 17-23 July 2005 (21)
© Peter Breedveld
Physical systems modeling • all concepts used in the model are, or have a direct relation to, physically relevant concepts (use of a priori knowledge!) • physical relations are maintained as much as possible • herein constrained to (for sake of ‘simplicity’): – deterministic mathematical models of macroscopic systems that • obey basic principles of macroscopic physics: – energy conservation – positive entropy production
• describe the behavior in time of the common physical properties: – – – – –
mechanical (incl. hydraulic and pneumatic) electrical magnetic chemical material
Summerschool Bertinoro, 17-23 July 2005 (22)
© Peter Breedveld
time & uncertainty • time:
– derived measure for ‘regularity’ based on counting ‘ticks’ of a ‘time-base’ based on repetitive behavior requiring state and change (in order to be able to count) – within the smallest unit used: necessarily
uncertainty
• cf. Heisenberg u.r. for displacement (= elastic state, kinetic state of change) and momentum (= elastic state of change, kinetic state)
– dialectic concepts!
Summerschool Bertinoro, 17-23 July 2005 (23)
© Peter Breedveld
concepts of ‘state’ & ‘change’ • dialectic fundamental concepts • basis for any dynamic model (more than time!) • within a context of discovery a shift to space-time is understandable and useful, like the shift from positionmomentum in Hamiltonian mechanics to position-velocity in Lagrangian mechanics Summerschool Bertinoro, 17-23 July 2005 (24)
© Peter Breedveld
Contents of the sequel • • • • •
Modeling pitfalls Port-based modeling Basic Concepts (ports, bonds) Dynamic conjugation (effort, flow) Multidomain modeling and the role of energy • (Computational) Causality
Summerschool Bertinoro, 17-23 July 2005 (25)
© Peter Breedveld
Modeling pitfalls • ‘Every model is wrong’ • Model depends on problem context • Competent models • Analogies are not identities • Avoid implicit assumptions • Avoid model extrapolation • Confusion of components with elements
Summerschool Bertinoro, 17-23 July 2005 (26)
© Peter Breedveld
Port-based modeling • multidomain approach:
– ‘mechatronics’ (and beyond)
• multiple view approach: other graphical representations: iconic diagrams, (linear graphs), block diagrams, bond graphs, etc. and equations • domain independent notation using ports: – bond graphs (& some other benefits…)
• port-based approach:
– underlying structure of 20-sim, ideal tool for demonstration
• what are ports and what are bond graphs? Summerschool Bertinoro, 17-23 July 2005 (27)
© Peter Breedveld
Physical components versus ideal elements
Summerschool Bertinoro, 17-23 July 2005 (28)
© Peter Breedveld
Physical component versus ideal elements
Summerschool Bertinoro, 17-23 July 2005 (29)
© Peter Breedveld
Physical component versus ideal elements
Summerschool Bertinoro, 17-23 July 2005 (30)
© Peter Breedveld
Component mounter
Summerschool Bertinoro, 17-23 July 2005 (31)
© Peter Breedveld
Physical components versus conceptual elements physical component: piece of rubber hose dominant behavior: •when falling: ideal mass •when pulling load: ideal spring •for vibration isolation: ideal resistor •etc. Summerschool Bertinoro, 17-23 July 2005 (32)
© Peter Breedveld
‘Parasitic’ elements • Next to dominant behavior:
Summerschool Bertinoro, 17-23 July 2005 (33)
© Peter Breedveld
In engineering models: • Avoid implicit assumptions, e.g. about – problem context – reference – orientation – coordinates – metric – ‘negligible’ phenomena, etc.
• Avoid model extrapolation – danger of ignoring earlier assumptions
• Focus at competence, not ‘truth’ Summerschool Bertinoro, 17-23 July 2005 (34)
© Peter Breedveld
Intuitive introduction of the ‘port’ concept
J
ideal current source
ideal motor
Dominant behavior
ideal inertia
ideal transmission
P
potentiometer
(not necessarily competent in each context) Summerschool Bertinoro, 17-23 July 2005 (35)
© Peter Breedveld
Simple model
Summerschool Bertinoro, 17-23 July 2005 (36)
© Peter Breedveld
Addition of relevant parasitic behavior Depending on problem context:
Polymorphic modeling Summerschool Bertinoro, 17-23 July 2005 (37)
© Peter Breedveld
What are ports? ∆ u2
∆u1
i
P =∆u⋅ i
(power)
i
∆ u3
i
i ∆ u4
∆T1
ω
electrical ports
∆T2
∆T3 J
ω
ω
MECH
P =∆T⋅ω (power)
mechanical ports
Summerschool Bertinoro, 17-23 July 2005 (38)
© Peter Breedveld
Comparison with familiar model views Iconic diagrams (‘ideal physical models’): R
L
i
Usource
v C
F
R
m
Fext K=1/C
electrical
Summerschool Bertinoro, 17-23 July 2005 (39)
mechanical
© Peter Breedveld
Iconic diagram symbols
Port -based, but Port-based, domain dependent
Summerschool Bertinoro, 17-23 July 2005 (40)
© Peter Breedveld
Iconic diagram symbols C-type storage I -type storage (M)R (dissipation, irreversible transduction) Se (effort source) Sf (flow source) TF (transformer) Summerschool Bertinoro, 17-23 July 2005 (41)
© Peter Breedveld
Iconic diagrams • Icons do not (always) represent physical structure:
20-sim icon: 1 i
TF (transformer) Summerschool Bertinoro, 17-23 July 2005 (42)
© Peter Breedveld
Common block diagram models R
L
Usource
i
v C
F
R
m
Fext K=1/C
Summerschool Bertinoro, 17-23 July 2005 (43)
© Peter Breedveld
Common block diagram models
Note: not all signals are physically meaningful variables Summerschool Bertinoro, 17-23 July 2005 (44)
© Peter Breedveld
Ports in iconic diagrams
R
L
Usource
i
v C
F
R
m
Fext K=1/C
Summerschool Bertinoro, 17-23 July 2005 (45)
© Peter Breedveld
Ports in iconic diagrams
R
m
L
F
i Usource
R Fext C
v K=1/C
Summerschool Bertinoro, 17-23 July 2005 (46)
© Peter Breedveld
Ports in block diagrams
Note: each signal is a physically meaningful variable
Summerschool Bertinoro, 17-23 July 2005 (47)
© Peter Breedveld
Structure as multiport
Summerschool Bertinoro, 17-23 July 2005 (48)
© Peter Breedveld
Different forms
Summerschool Bertinoro, 17-23 July 2005 (49)
© Peter Breedveld
Different forms
Summerschool Bertinoro, 17-23 July 2005 (50)
© Peter Breedveld
Compact & domain independent Domain dependent
Coil
Mass
Source
Spring
Source
Capacitor
Resistor
Summerschool Bertinoro, 17-23 July 2005 (51)
Damper
© Peter Breedveld
Compact & domain independent Domain independent
I
I
Se
1
Structure explicitly
R
C
Se
1
C
R
represented as multiport: JUNCTION!
Summerschool Bertinoro, 17-23 July 2005 (52)
© Peter Breedveld
Ports in bond graph view I
Se
1
C
R
Summerschool Bertinoro, 17-23 July 2005 (53)
© Peter Breedveld
Bond graphs
Inventor (MIT, 1959): Prof. Henry M. Paynter (1923-2002) Summerschool Bertinoro, 17-23 July 2005 (54)
© Peter Breedveld
Analog simulation with OP-AMPS
Summerschool Bertinoro, 17-23 July 2005 (55)
© Peter Breedveld
Bonds Storage of kinetic energy
Storage of elastic energy
Compact notation:
C
I Energy exchange = power
•Terminology and notation induced by chemical bonds:
Summerschool Bertinoro, 17-23 July 2005 (56)
© Peter Breedveld
Ports and power bonds • (power) bond connecting two elements via (power) ports (Harold Wheeler, 1949)
Summerschool Bertinoro, 17-23 July 2005 (57)
© Peter Breedveld
mass-spring system
Summerschool Bertinoro, 17-23 July 2005 (58)
© Peter Breedveld
Dynamic conjugation • Between signals of bilateral signal flow of relation: – rate of change: ‘flow’ (zero in equilibrium) • e.g. molar rate during diffusion
– equilibrium determining variable: ‘effort’ • e.g. concentration
Summerschool Bertinoro, 17-23 July 2005 (59)
© Peter Breedveld
Power conjugation = special case of dynamic conjugation : – –
‘effort’ and ‘flow’ relate to power functional relation is commonly a product – sum in case of scattering variables
dE ∂E dqi ∂E dp j P= =∑ +∑ = ∑ ei fi + ∑ e j f j dt d d t t i ∂qi j ∂p j i j dp j ∂E 'effort' ei = ej = dt ∂qi fi =
dqi dt
Summerschool Bertinoro, 17-23 July 2005 (60)
fj =
∂E ∂p j
'flow'
© Peter Breedveld
Conjugate variables & corresponding states
Summerschool Bertinoro, 17-23 July 2005 (61)
© Peter Breedveld
Basic dynamic behaviors & mnemonic codes • Storage (reversible) – C, I (q-type and p-type storage) • Irreversible transformation (‘dissipation’): – (M)R(S) • Distribution – 0-junction, 1-junction • Supply and demand: – (M)Se, (M)Sf • Reversible transformation – (M)TF, (M)GY • >> 9 basic elements Summerschool Bertinoro, 17-23 July 2005 (62)
© Peter Breedveld
Orientation conventions • 1-ports: into most elements (R, C, I), except sources (Se, Sf) • 2-port transducers: 1 in, 1 out • junction structure elements: arbitrary • multiport generalizations: same as simple form • MOST IMPORTANT: obeying grammar rules minimizes sign errors!
Summerschool Bertinoro, 17-23 July 2005 (63)
© Peter Breedveld
Convention on effort and flow positions
Summerschool Bertinoro, 17-23 July 2005 (64)
© Peter Breedveld
Basic one-port elements
Summerschool Bertinoro, 17-23 July 2005 (65)
© Peter Breedveld
Basic two- and multiport elements
Summerschool Bertinoro, 17-23 July 2005 (66)
© Peter Breedveld
Constitutive relations not necessarily linear! e.g. zener diode:
dominant behavior: irreversible transduction (resistor)
with nonlinear constitutive relation
R Summerschool Bertinoro, 17-23 July 2005 (67)
© Peter Breedveld
Other form of non-linearity: Modulation e.g. modulation of a transducer:
iconic diagram (IPM)
capstan MR Summerschool Bertinoro, 17-23 July 2005 (68)
© Peter Breedveld
Bilateral signal flow (computational causality) 2 possibilities:
Summerschool Bertinoro, 17-23 July 2005 (69)
© Peter Breedveld
Notation Causal stroke:
Summerschool Bertinoro, 17-23 July 2005 (70)
© Peter Breedveld
Back to the example I
Se
1
C
R
Summerschool Bertinoro, 17-23 July 2005 (71)
© Peter Breedveld
Back to the example I
Se
1
C
R
20-sim demo Summerschool Bertinoro, 17-23 July 2005 (72)
© Peter Breedveld
Causal port properties • Fixed causality
• Preferred (integral) causality
Summerschool Bertinoro, 17-23 July 2005 (73)
© Peter Breedveld
Causal constraints
Summerschool Bertinoro, 17-23 July 2005 (74)
© Peter Breedveld
Arbitrary causality
Summerschool Bertinoro, 17-23 July 2005 (75)
© Peter Breedveld
Causality assignment Algorithmic (SCAP):
– 1) fixed causal ports with propagation via constraints causal ‘conflict’: ‘ill-posedness’
– 2) preferred causal ports with propagation via constraints
causal ‘conflict’: dependent state(s)
– 3) choice of arbitrary causal port with propagation via constraints means existence of algebraic loop(s)
Summerschool Bertinoro, 17-23 July 2005 (76)
© Peter Breedveld
Example u
u
u T
i
i
u
u
i
J ω
i
i
I
1
I
GY
ω •Important modeling feedback
2nd order loop
Se
T
•Automatic in 20-sim •Visible in bond graph causality •‘Hidden’ in case of iconic diagrams
1
1st order loop
1st order loop
R Summerschool Bertinoro, 17-23 July 2005 (77)
R © Peter Breedveld
Positive orientation orientation NOT THE SAME AS direction!
Summerschool Bertinoro, 17-23 July 2005 (78)
© Peter Breedveld
Causality does NOT affect orientation!
Summerschool Bertinoro, 17-23 July 2005 (79)
© Peter Breedveld
Basic form of physical models interfaces and structure have to satisfy basic laws in degenerate way storage = 0, production = 0
Basic ‘laws’
boundary addition
structure constraint relations distribution, reversible transformation transport
first ‘law’ conservation ‘laws’ storage
second ‘law’ irreversible transformation
constitutive relations determine particular occurrence (instantiation) Summerschool Bertinoro, 17-23 July 2005 (80)
© Peter Breedveld
Thermodynamic approach System boundary Energy balance
GJS GJS
for the total system
Irreversible thermodynamics of ‘forces’ and ‘fluxes’ (irreversible part)
Gibbs’ relation (Energy ‘balance’ of the reversible part)
Power continuous junction structure
• Difference is used to find the entropy production: Structure is conceptual ! • Power continuous structure not made explicit, except for instantaneous Carnot engines (1-junction for entropy flow) Summerschool Bertinoro, 17-23 July 2005 (81)
© Peter Breedveld
Mechanical versus Thermodynamic framework Mechanics:
Thermodynamics:
Two types of storage One type of storage Oscillatory behavior Only relaxation behavior: C-R (damped): C-I(-R) Split domains (therm.) and couple by SGY (mech.): C-SGY-C
Summerschool Bertinoro, 17-23 July 2005 (82)
© Peter Breedveld
Symplectic gyrator • Unit gyrator:
C
0
GY
0
C
0
C
r=1
• SGY: • Multiport:
C
C11 C12 C12 C22
0
C
SGY
0
SGY
0 -1 +1 0
20-sim demo
Summerschool Bertinoro, 17-23 July 2005 (83)
© Peter Breedveld
Symplectic GYrator • Mechanical (x,p)
ekin =
∂E ∂p
C
ekin = v
0 f kin =
dp dt
SGY F =−
e pot =
e pot = F
dp dt
0 f pot =
f pot = v
∂E ∂x dx dt
C
dp
– Only in inertial frames: F = − (Newton's 2nd law) dt
• Electrical network (q,λ): emag =
C f mag =
∂E ∂λ dλ dt
∂E eelec = u eelec = ∂q
emag = i
0 −
dλ =u dt
SGY
– Only quasi-stationary (non-radiating): Summerschool Bertinoro, 17-23 July 2005 (84)
−
0 f elec = i
dq felec = dt
dλ =u dt
emag = i
C
© Peter Breedveld
Unit gyrator (SGY) as ‘dualizer’ Original:
Port equivalent:
Original:
Port equivalent:
0
SGY
I
0
C
Se
SGY
1
Sf
1
1
SGY
C
1
I
Sf
SGY
0
Se
0
1
SGY
R
1
R 1
SGY
1 1
SGY
TF
SGY
0
1
TF
SGY
0
1
0
0
0 SGY 0
SGY
1
GY
SGY
1
1
GY
1 0
1
SGY
TF
1
1
GY
1 1
1
SGY
GY
0
1
TF
SGY
0
SGY
0
1
1
0 SGY
1
0 1
Summerschool Bertinoro, 17-23 July 2005 (85)
© Peter Breedveld
0
Basic model structure (generalized bond graph): Generalized Anti Junction -reciprocal Structure (power part (generalized Boundary conditions continuous)gyrator) (sources)
Irreversible Energy storage Junction Weighted transformation Structure (reciprocal) (entropy production) Summerschool Bertinoro, 17-23 July 2005 (86)
Weighted part (generalized transformer)
Simple Junction Structure (generalized Kirchhoff laws) © Peter Breedveld
Mechanical framework of variables
Summerschool Bertinoro, 17-23 July 2005 (87)
© Peter Breedveld
Generalized thermodynamic framework of variables electric
magnetic
thermal
chemical
f flow
e effort
i current
u voltage
u voltage
i current
T temperat ure
fS entropy flow
µ
fN molar flow
chemical potential
Summerschool Bertinoro, 17-23 July 2005 (88)
∫
q = f dt
generalized state
q = ∫ idt
charge
λ = ∫ udt
magnetic flux linkage
S=
∫ fS dt
entropy
N=
∫ fN dt
number of moles
© Peter Breedveld
Generalized thermodynamic framework of variables f flow
e effort
elastic/potential translation
v velocity
F force
kinetic translation
F force
v velocity
elastic/potential rotation
ω angular velocity
T torque
kinetic rotation
T torque
ω angular velocity
ϕ volume flow
p pressur e
p pressure
ϕ volume flow
elastic hydraulic
kinetic hydraulic
Summerschool Bertinoro, 17-23 July 2005 (89)
∫
q = f dt
generalized state x = ∫ v dt
displacement p = ∫ Fdt
momentum θ = ∫ ωdt
angular displacement b = ∫ Tdt
angular momentum V = ∫ ϕdt
volume Γ=
∫ pdt
momentum of a flow tube
© Peter Breedveld
Energy as starting point • Energy
( )
– E = E q = E ( q1 ,..., qi ,...qn ) – Homogeneous function of set of state variables
• In (‘generalized’) mechanics : Hamiltonian
( )
– E = H q , p = H ( q1 ,..., qi ,...qk , p1 ,..., pi ,... pk )
• In thermodynamics of ‘simple’ systems: U = U ( q ) = U ( S , q ,..., q ,...q ) internal energy
( )
2
S = S (U , q2 ,..., qi ,...qn ) – E = U q = U (V , S , N ) – more species: E = U q = U (V , S , N ) =
i
n
( )
= U (V , S , N1 ,..., N i ,..., N m ) = = U (V , S , N1 ,..., N i ,..., N m −1 , N ) Summerschool Bertinoro, 17-23 July 2005 (90)
© Peter Breedveld
Configuration influence on energetic structure • can become an additional energy state: geometric parameter in storage relation results in a force
F ( q, x ) =
∂E ( q, x ) ∂x
– examples: LVDT, electret microphone, relay, (ideal) gas etc. (MP C) – cycles allow transduction – changes of causality correspond to Legendre transforms!! (relation to dissipation)
• can modulate an energy relation
– examples: crank-slider mechanism, etc. (MTF)
• can switch a contact (behavior) Summerschool Bertinoro, 17-23 July 2005 (91)
© Peter Breedveld
Modulated capacitor?
Summerschool Bertinoro, 17-23 July 2005 (92)
© Peter Breedveld
System (boundary) definition • Open (thermodynamic) systems • dN = 0: ‘Lagrangian coordinates’: N not a state • dV = 0: ‘Eulerian coordinates’: V not a state • Almost always mixed boundaries: N and V remain states! • Network: globally Eulerian, locally Lagrangian (e.g. network of tubes) Summerschool Bertinoro, 17-23 July 2005 (93)
© Peter Breedveld
‘Classical network case’ • Both N and V :
– constant – NOT CONSIDERED STATES, but parameters
• As a result:
– energy function not necessarily first order homogenous – no dependency between storage port efforts (no generalized Gibbs-Duhem equation) – no clear distinction between intensive and extensive description
Summerschool Bertinoro, 17-23 July 2005 (94)
© Peter Breedveld
If not a network: • Globally Lagrangian and locally Eulerian possible! • Example: balloon
Summerschool Bertinoro, 17-23 July 2005 (95)
© Peter Breedveld
Part of a tube network daxialV=0 (V not a convected property) daxial N≠0
dradialN=0 dradial V≠0 (flexible tube, local changes)
Summerschool Bertinoro, 17-23 July 2005 (96)
© Peter Breedveld
System boundary definition • Mechanical system of masses, springs, dampers, etc. (network): – globally Eulerian coordinates
• because the network (topology) of spatial elements ‘cells’ is constant
– locally Lagrangian coordinates
• e.g. the constant mass with its motion restricted to the element cell: ‘one element passing a neighboring element destroys the topology’
Summerschool Bertinoro, 17-23 July 2005 (97)
© Peter Breedveld
Summary •
System theoretic approach to physical system dynamics based on – classification of phenomena in terms of energy – fundamental principles of thermodynamics
•
as a result: – variables and relations describing physical systems may be classified – models may be organized as • multiport elements • interconnected in an interconnection structure corresponding to a generalized network
– multiport elements describe basic behaviors with respect to energy and entropy
•
Explicit conceptual separation between two roles of displacement variable: – energy state (stored potential or elastic energy) – configuration state (‘position in space’)
•
strongly reduces conceptual complexity even though numerical complexity may become temporarily higher
Summerschool Bertinoro, 17-23 July 2005 (98)
© Peter Breedveld