next up previous
Next: 3 Analyzing Model Discontinuities Up: 2 The Hybrid Modeling Previous: 2.1 Controlled Junctions

2.2 Mode Switching in the Hybrid Model

 

Discontinuous effects in physical system models occur when energy or power variables cross a certain threshold value. For example, a diode is modeled to come on when the voltage drop across it exceeds 0.6V. These discontinuous effects establish or break energetic connections in the model, and this may cause related signals to change discontinuously. The tank system, Fig. 2, has two pressure valves that open when either pressure tex2html_wrap_inline686 or tex2html_wrap_inline688 crosses a pre-set threshold value. If the specified threshold value for the valve controlled by tex2html_wrap_inline686 is tex2html_wrap_inline692 the threshold value of the second valve, the opening of the first valve will cause the second to open. If the resistance, capacity, and flow inertia of the connecting pipe are abstracted out of the model, these two discrete changes occur instantaneously and the mode where only the valve controlled by tex2html_wrap_inline686 is open occurs for a brief instant in time (according to the modeling assumptions it is instantaneous). Therefore, discontinuous changes may cause other signal value changes which result in a sequence of discontinuous changes. Since modeling assumptions require that discontinuous changes occur instantaneously, these transitions are called mythical.

   figure96
Figure 2: Tank with two pressure controlled valves.

In this paper we establish that state variables in the system do not change during mythical changes. The system does undergo discrete state changes, determined by the finite state automata associated with the controlled junctions. Eventually, a sequence of discrete switches terminates in a real mode (system behavior again evolves as a function of time), and the continuous state vector for this new mode has to be derived. This is illustrated in Fig. 3. Mythical modes are depicted in a white background and real modes in a dark background. In real mode tex2html_wrap_inline696 a signal value crosses a threshold at time tex2html_wrap_inline698 , which causes a discontinuous change to model configuration tex2html_wrap_inline700 , represented by the discrete state vector tex2html_wrap_inline702 . The power variable values tex2html_wrap_inline704 in this new configuration are calculated from the original energy distribution tex2html_wrap_inline706 values. If the new values cause another instantaneous mode change, a new mode tex2html_wrap_inline708 is reached, where the new power variables values, ( tex2html_wrap_inline710 , tex2html_wrap_inline712 ), are calculated from the original energy distribution tex2html_wrap_inline706 . Further mythical mode changes may occur till a real mode, tex2html_wrap_inline716 , is reached. The final step involves mapping the energy distribution, or continuous state variable values, of the departed real mode to the new real mode. This issue is non-trivial, and discussed in detail in a later section. Real time continuous simulation resumes at tex2html_wrap_inline718 so system behavior in real time implies mode tex2html_wrap_inline716 follows tex2html_wrap_inline696 . The formal Mythical Mode Algorithm (MMA) is outlined below.

  1. Calculate the energy values ( tex2html_wrap_inline724 ) and signal values ( tex2html_wrap_inline726 ) for bond graph model tex2html_wrap_inline696 using ( tex2html_wrap_inline730 , tex2html_wrap_inline732 ), values at the previous simulation step as initial values.
  2. Use CSPEC to infer a possible new mode given ( tex2html_wrap_inline726 ).
  3. If one or more controlled junctions switch states then:
    1. Derive the bond graph for this configuration.
    2. Assign causal strokes using the SCAP algorithm[12].
    3. Calculate the signals ( tex2html_wrap_inline736 ) for the new mode, tex2html_wrap_inline738 , based on the initial values ( tex2html_wrap_inline724 ).
    4. Use CSPEC again to infer a possible new mode based on ( tex2html_wrap_inline736 ) for the new mode, tex2html_wrap_inline738 .
    5. Repeat step 3 till no more mode changes occur.
  4. Establish the final mode, tex2html_wrap_inline716 , as the new system configuration.
  5. Map tex2html_wrap_inline748 to the energy distribution for tex2html_wrap_inline716 , tex2html_wrap_inline752 .

   figure125
Figure 3: State Evolution: Mythical + Real Modes.

A complete simulation system that incorporates continuous simulation and the MMA algorithm been implemented and tested on a number of physical system examples[9].


next up previous
Next: 3 Analyzing Model Discontinuities Up: 2 The Hybrid Modeling Previous: 2.1 Controlled Junctions

Pieter J. Mosterman
Mon Jul 21 19:58:19 CDT 1997