a large class of physical and engineering problems can be described in terms of lagrangian and hamiltonian systems , such as controlled mechanical systems , externally forced and dissipative systems , collisions , electric circuits , stochastic systems , and field theories such as electromagnetism and elasticity .
such systems are usually categorized into lagrangian or hamiltonian systems with constraints , whether holonomic or nonholonomic ; however , one may be faced with serious problems on ever larger scales becoming more and more heterogeneous , involving a mixture of mechanical and electrical components , with flexible and rigid parts , and elastic or magnetic couplings ( see , for instance , @xcite and @xcite ) .
dirac structures are known as a powerful tool for providing a natural geometric framework for describing the structure of such diverse lagrangian and hamiltonian systems on the product of configuration spaces by introducing nontrivial constraints that model interactions between subsystems , which is called an _
interconnection_. the notion of an interconnection was originally developed by @xcite in his book of `` diakoptics '' for modeling complicated physical systems .
the word `` diakoptics '' denotes a procedure of _ tearing _ a dynamical system into disconnected subsystems as well as _ interconnecting _ them to reconstruct the original system .
the original system may be regarded as an interconnected system of subsystems , with constraints induced from the interconnection .
interconnection has been studied through the lens of power conservation motivated by the bond graph theory established by @xcite . in electric circuit theory , the interconnection network has been widely employed through nonenergic multiports ( see @xcite ) , where kirchhoff s current law gives a constraint distribution on a configuration charge space such that tellegen s theorem holds . in mechanics ,
it was shown by @xcite that kinematical constraints due to mechanical joints , coordinate transformations as well as force equilibrium conditions in dalembert principle can be represented by the interconnection modeled by nonenergic multiports using the bond graph theory in the context of lagrangian mechanics .
it was demonstrated by @xcite that the interconnection can be represented by dirac structures induced from poisson structures and also that nonholonomic systems and l - c circuits can be represented by _
implicit hamiltonian systems_. on the lagrangian side , it was shown by @xcite that nonholonomic mechanical systems and l - c circuits ( as degenerate lagrangian systems ) can be formulated by _
implicit lagrangian systems _ associated with the dirac structures on the cotangent bundles induced from given constraint distributions . from the view point of variational structures ,
it was shown by @xcite that the standard implicit lagrangian systems , namely , _
implicit euler - lagrange equations _ for unconstrained systems can be derived from the _ hamilton - pontryagin principle _ and also that the lagrange - dirac dynamical systems can be formulated in the context of _ lagrange - dalembert - pontryagin principle_.    from the viewpoint of control theory , a notion of _ implicit port - controlled hamiltonian ( ipch ) systems _ ( an implicit hamiltonian system with external control inputs ) was developed by @xcite ( see also @xcite , @xcite and @xcite ) and much effort has been devoted to the passivity based control for the interconnected mechanical systems ( @xcite ) . for the case of regular lagrangian systems , the equivalence between _ controlled lagrangian ( cl ) systems _ and _ controlled hamiltonian ( ch ) _ systems was shown by @xcite and the reduction theory for the cl systems was shown by @xcite . for the general case
in which a given lagrangian is degenerate , an implicit lagrangian analogue of ipch systems , namely , a notion of _ implicit port - controlled lagrangian ( ipcl ) systems _ for electrical circuits was constructed by @xcite and @xcite , where it was shown that l - c transmission lines can be represented in the context of the iplc system by employing induced dirac structures with lagrange multipliers .
recently , it was shown by @xcite that any power - conserving interconnection of ipch systems is another port - hamiltonian system , with the dirac structure being the _ composition _ of the dirac structures of its constituent parts , hamiltonian the sum of the hamiltonians , and resistive relations determined by the individual resistive relations .
as a control application , the feedback interconnection of a plant port - hamiltonian system with a controller port - hamiltonian system can be represented by the composition of a plant dirac structure with a controller dirac structure .
[ [ goals - and - organizations - of - the - paper . ] ] goals and organizations of the paper .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    one of the main goals in this paper is to develop a simple expression for the interconnection of distinct dirac structures and their associated mechanical systems ; namely , we will show how to obtain an interconnection of distinct dirac structures @xmath4 on distinct manifolds @xmath5 as another dirac structure @xmath6 on @xmath7 by `` interconnecting '' @xmath4 . to do this ,
we will introduce two mathematical ingredients , namely , the _ dirac sum _ @xmath8 and the _ bowtie product _ @xmath2 of dirac structures and we will show how two distinct dirac structures can be interconnected to yield another dirac structure via an _ interconnection _ dirac structure ( usually labeled @xmath1 in this paper ) under an assumption of transversal intersections . in particular , the interconnection of dirac structures is @xmath9 .
the interconnection dirac structure @xmath1 provides an interface to interconnect `` distinct '' dirac structures @xmath4 to yield another dirac structure @xmath6 .
we will also show how the associated interconnection of lagrange - dirac dynamical systems can be developed in the context of variational structures .
we will demonstrate our theory by some examples of circuits and mechanical systems .
this paper consists of the following parts :    1 .
we briefly review dirac structures and we show how associated lagrange - dirac dynamical systems can be developed by the _ lagrange - dalembert - pontryagin principle _ ( see @xcite ) .
we show how distinct dirac structures can be interconnected to yield another dirac structure on the product manifold through an interconnection dirac structure .
to do this , we develop a map called the dirac sum @xmath8 and a tensor product called the bowtie product @xmath2 between dirac structures .
we explore the variational structure of interconnection of lagrange - dirac dynamical systems ; namely , the interconnected lagrange - dirac dynamical system can be formulated in the context of the lagrange - dalembert -pontryagin principle .
we demonstrate some illustrative examples of circuits and nonholonomic systems , which can be generalized to more complex systems .
[ [ acknowledgements . ] ] acknowledgements .
+ + + + + + + + + + + + + + + + +    we are also very grateful to henrique bursztyn , mathieu desbrun , melvin leok , tomoki ohsawa , arjan van der schaft and alan weinstein for useful remarks and suggestions .
the research of h. j. is supported by caltech and jst - crest in part .
the research of h.  y. is partially supported by jsps grant - in - aid 20560229 , jst - crest and waseda university grant for sr 2010a-606 .
[ [ linear - dirac - structures . ] ] linear dirac structures .
+ + + + + + + + + + + + + + + + + + + + + + + +    first , let us recall the definition of a linear dirac structure , namely , a _ dirac structure on a vector space _ @xmath10 ; we assume that @xmath10 is finite dimensional for simplicity ( see , @xcite ) . let @xmath11 be the dual space of @xmath10 , and @xmath12 be the natural paring between @xmath11 and @xmath10 .
define the symmetric paring @xmath13 on @xmath14 by    @xmath15    for @xmath16 .
a _ dirac structure _ on @xmath10 is a subspace @xmath17 such that @xmath18 , where @xmath19 is the orthogonal of @xmath6 relative to the pairing @xmath20 .
[ [ dirac - structures - on - manifolds . ] ] dirac structures on manifolds .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    let @xmath21 be a smooth manifold and let @xmath22 denote the whitney sum bundle over @xmath21 , namely , the bundle over the base @xmath21 and with fiber over the point @xmath23 equal to @xmath24 .
a subbundle @xmath25 is called an almost dirac structure on @xmath21 , when @xmath26 is a dirac structure on the vector space @xmath27 at each point @xmath28 .
a given two - form @xmath29 on @xmath21 together with a distribution @xmath30 on @xmath21 determines an _ ( almost ) dirac structure _ on @xmath21 as follows : for each @xmath28 , @xmath31 where @xmath32 is the restriction of @xmath29 to @xmath30 .
we can also restate the dirac structure as , for each @xmath28 , @xmath33 where @xmath34 is an annihilator of @xmath30 .
we call @xmath6 an _ integrable dirac structure _ if the following integrability condition @xmath35 is satisfied for all pairs of vector fields and one - forms @xmath36 , @xmath37 , @xmath38 that take values in @xmath6 , where @xmath39 denotes the lie derivative along the vector field @xmath40 on @xmath21 .    [ [ remark . ] ] remark
. + + + + + + +    let @xmath41 be a space of local sections of @xmath22 , which is endowed with the skew - symmetric bracket @xmath42 : \gamma(tm \oplus t^{\ast}m ) \times   \gamma(tm \oplus t^{\ast}m ) \to   \gamma(tm \oplus t^{\ast}m)$ ] defined by @xmath43 & : = \left ( \left[x_{1 } , x_{2 } \right ] ,   \pounds_{x_1 } \alpha_{2}- \pounds_{x_2 } \alpha_{1 } + \mathbf{d}\left < \alpha_{2 } , x_{1 } \right>\right)\\ & \;=\left ( \left[x_{1 } , x_{2 } \right ] ,   \mathbf{i}_{x_1 } \mathbf{d}\alpha_{2}- \mathbf{i}_{x_2 } \mathbf{d}\alpha_{1 } + \mathbf{d } \left < \alpha_{2 } , x_{1 } \right>\right ) . \end{split}\ ] ] this bracket is the one that @xcite originally developed and it does not necessarily satisfy the jacobi identity .
it was shown by @xcite that the integrability condition of the dirac structure @xmath44 given in equation can be expressed as @xmath45 \subset \gamma(d),\ ] ] which is the closedness condition of the courant bracket ( see also @xcite and @xcite ) .
the following lemma is useful and the proof may be found in @xcite and @xcite .
[ lem : isotropic ] a subbundle @xmath46 is isotropic with respect to @xmath20 if and only if @xmath47 .    [ [ induced - dirac - structures . ] ] induced dirac structures .
+ + + + + + + + + + + + + + + + + + + + + + + + +    one of the most important and interesting dirac structures in mechanics is an _ induced dirac structure _ from kinematic constraints , whether holonomic or nonholonomic .
such constraints are generally given by a distribution on a configuration manifold ( as to the details , refer to @xcite ) .
let @xmath48 be a configuration manifold .
let @xmath49 and @xmath50 be the tangent and cotangent bundles of @xmath48 .
let @xmath51 be a regular distribution on @xmath48 and define a lifted distribution on @xmath50 by @xmath52 where @xmath53 is the cotangent projection so that its tangent is a map @xmath54 . let @xmath29 be the canonical two - form on @xmath50 . the _ induced dirac structure _
@xmath6 on @xmath50 is the subbundle of @xmath55 , whose fiber is given for each @xmath56 as @xmath57 this is , of course , a special instance of the construction .
notice that the induced dirac structure can be restated by using the bundle map @xmath58 as follows : @xmath59 where @xmath60 is the annihilator of @xmath61 .
[ [ remark.-1 ] ] remark .
+ + + + + + +    if there exists no constraint , namely , the case in which @xmath62 , the dirac structure @xmath6 corresponds to the canonical dirac structure , which is given by a graph of the bundle map @xmath63 associated to @xmath29 .
[ [ local - expressions . ] ] local expressions .
+ + + + + + + + + + + + + + + + + +    let @xmath10 be a model space for @xmath48 and @xmath64 be the range in @xmath10 of a chart on @xmath48 .
then , the range of the induced chart on @xmath49 is represented by @xmath65 , while @xmath50 is locally given by @xmath66 . in this representation ,
@xmath67 are local coordinates of @xmath49 and @xmath68 are local coordinates of @xmath50 .
further , @xmath69 is represented by @xmath70 , while @xmath71 is locally given by @xmath72 . in this local representation ,
@xmath73 are the corresponding coordinates of @xmath69 , while @xmath74 are the local coordinates of @xmath71 .    using @xmath75
locally denoted by @xmath76 and @xmath77 , it follows that @xmath78 and the annihilator of @xmath61 is locally represented as @xmath79 since we have the local formula @xmath80 , the condition @xmath81 reads @xmath82 and @xmath83 thus , the induced dirac structure is locally represented by @xmath84 where @xmath85 is the polar of @xmath86 .
in this section , let us recall the definition of implicit lagrangian systems or lagrange - dirac dynamical systems by following @xcite .
[ [ lagrangedirac - dynamical - systems . ] ] lagrange  dirac dynamical systems .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    let @xmath87 be a lagrangian , possibly degenerate .
let @xmath6 be a dirac structure on @xmath50 induced from a given distribution @xmath51 .
recall that a _ partial vector field _
@xmath88 is defined as a mapping such that @xmath89 , where @xmath90 and @xmath91 is the tangent projection . as a mapping @xmath92 such that @xmath93 the partial vector field @xmath94 is a map that assigns to each point @xmath95 , a vector in @xmath69 at the point @xmath96 ;
we write @xmath40 as @xmath97 where we note that @xmath98 and @xmath99 are understood to be functions of @xmath100 .
let @xmath101 be the _
generalized energy _ given by @xmath102    an implicit lagrangian system or a lagrange
dirac dynamical system is the triple @xmath103 that satisfies , for each @xmath95 and with @xmath104 , namely , @xmath105 , @xmath106 where @xmath107 is the legendre transformation . in the above , @xmath108 is a linear function on @xmath109 at each point @xmath95 and hence the restriction @xmath110 may be regarded as a linear function on @xmath111 such that @xmath112 .
[ [ local - expressions.-1 ] ] local expressions .
+ + + + + + + + + + + + + + + + + +    it follows from that the lagrange  dirac dynamical system @xmath113 is locally given by @xmath114    [ [ remark.-2 ] ] remark .
+ + + + + + +    if a partial vector field @xmath115 satisfies the condition for a lagrange  dirac dynamical system @xmath113 , then the legendre transformation @xmath112 is consistent with the equality of base points and the condition itself gives the second - order condition @xmath116 .
in other words , the partial vector field of the lagrange  dirac dynamical system is uniquely given on the graph of the legendre transformation .
for the case in which no kinematic constraint is imposed , i.e. , @xmath62 , we can develop the standard implicit lagrangian system , which is expressed in local coordinates as @xmath117 which we shall call implicit euler
lagrange equations .
note that the implicit euler
lagrange equations include the euler  lagrange equations @xmath118 , the legendre transformation @xmath119 as well as the second - order condition @xmath116 .
[ [ example - harmonic - oscillators . ] ] example : harmonic oscillators .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    we illustrate the implicit euler
lagrange equations by the example of a linear harmonic oscillator . in this case
, the configuration space is @xmath120 representing the position of a particle on the line .
the lagrangian is given by @xmath121 .
we find the generalized energy to be @xmath122 .
the canonical dirac structure @xmath123 can be easily constructed by using the canonical symplectic structure @xmath29 on @xmath50 .
the lagrange - dirac system is given by @xmath103 , which satisfies @xmath124 and it follows @xmath125 .
we see @xmath126 and @xmath127 .
so , the dynamics is equivalent to the following equations : @xmath128 together with the legendre transform @xmath129 .
[ [ lagrange - dirac - dynamical - systems - with - external - forces . ] ] lagrange - dirac dynamical systems with external forces .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    here , we consider the case of a system with external force fields .
let @xmath130 be a natural projection .
an external force field @xmath131 induces a horizontal one - form @xmath132 on @xmath133 such that , for @xmath95 , @xmath134 where @xmath135 and @xmath136 .    given a lagrangian @xmath87 ( possibly degenerate ) , a lagrange - dirac dynamical system with an external force field is defined by a _ quadruple _
@xmath137 , which satisfies , in coordinates @xmath95 , @xmath138 it follows that the local expression of a lagrange - dirac system in equation may be given by @xmath139 and with @xmath140 the curve @xmath141 , @xmath142 in @xmath133 that satisfies the condition is a solution curve of @xmath143 .
[ [ power - balance - law . ] ] power balance law .
+ + + + + + + + + + + + + + + + + +    the time derivative of @xmath144 restricted to the solution curve @xmath141 reads @xmath145 where we have employed @xmath146 and @xmath147 .    [ [ example - the - harmonic - oscillator - with - damping . ] ] example : the harmonic oscillator with damping .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    as before , consider the harmonic oscillator in the setting of @xmath120 , @xmath148 , @xmath149 and @xmath150 and add a damping force @xmath151 defined by @xmath152 , where @xmath153 is a positive damping constant .
then , @xmath154 . the formulas in equation
give us the equations : @xmath155 with the legendre transformation @xmath129 .
[ [ the - hamilton - pontryagin - principle . ] ] the hamilton - pontryagin principle .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    lagrange - dirac dynamical systems can describe general scenarios such as non - holonomic degenerate lagrangian systems because of a firm basis in variational structures .
given a lagrangian @xmath87 , unconstrained equations of motion are described by hamilton s principle : @xmath156 from the viewpoint of dirac geometry , we start rather from the hamilton - pontryagin principle , which is given by the stationary condition of the action integral on the space of curves @xmath157 $ ] in @xmath133 as @xmath158 where @xmath159 .
the above variational principle can be restated by using the generalized energy @xmath160 as @xmath161 keeping the endpoints @xmath162 and @xmath163 of @xmath164 fixed , we can obtain the _ implicit euler  lagrange equations_.    [ [ the - lagrange - dalembert - pontryagin - principle . ] ] the lagrange - dalembert - pontryagin principle .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    consider a mechanical system with kinematic constraints that are given by a constraint distribution @xmath165 on @xmath48 .
then , the motion of the mechanical system @xmath166 \to q$ ] is said to be constrained if @xmath167 for all @xmath168 .
assume that the distribution @xmath169 is not involutive ; that is , @xmath170 \notin \delta_{q}(q)$ ] for any two vector fields @xmath171 on @xmath48 with values in @xmath169 .
further , let @xmath172 be a ( possibly degenerate ) lagrangian on @xmath49 and let @xmath131 be an external force field .
the lagrange - dalembert - pontryagin principle for a curve @xmath141 , @xmath173 in @xmath133 is represented by @xmath174 for a given variation @xmath175 and with the constraint @xmath176 . keeping the endpoints of @xmath164 fixed
, it follows @xmath177 for a chosen variation @xmath175 , for all @xmath178 and @xmath179 , and with @xmath180 .
the lagrange - dalembert - pontryagin principle gives the local expression of equations of motion for nonholonomic mechanical systems with external forces such that @xmath181    from equation , it reads @xmath182 which is satisfied for a given variation @xmath175 , for all @xmath178 and @xmath179 , and with the constraint @xmath180 .
thus , we obtain equation .    [ [ coordinate - representation . ] ] coordinate representation .
+ + + + + + + + + + + + + + + + + + + + + + + + + +    let @xmath48 be a finite dimensional smooth manifold whose dimension is @xmath183 .
choose local coordinates @xmath184 on @xmath48 so that @xmath48 is locally represented by an open set @xmath185 .
then , @xmath49 is locally given by @xmath186 .
similarly @xmath50 may be locally denoted by @xmath187 .
the constraint set @xmath169 defines a subspace on each fiber of @xmath49 , which we denote by @xmath188 at each point @xmath189 .
if the dimension of @xmath190 is @xmath191 , then we can choose a basis @xmath192 of @xmath193 .
recall that the constraint sets can be also represented by the annihilator @xmath194 , which is spanned by @xmath195 one - forms @xmath196 .
it follows that equation can be represented , in coordinates , by employing the lagrange multipliers @xmath197 , as follows : @xmath198 where we employ the local expression @xmath199 .
there exist various types of the power conserving interconnections between distinct dynamical systems such as two rigid bodies connected by ball - socket joints , mechanical systems connected by transmission gears , transducers between electric circuits and mechanical systems such as d - c motors , fluid and rigid body interactions , etc .
all the examples listed above may be represented by dirac structures . in this section
, we will show how distinct dirac structures can be interconnected via an _ interconnection _ dirac structure to yield another dirac structure of a fully interconnected dynamical system .
[ [ the - direct - sum - of - dirac - structures . ] ] the direct sum of dirac structures .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    before going into details on the case of @xmath183 distinct dirac structures , we will begin with the case of @xmath200 to introduce the direct sum of dirac structures on manifolds .
let @xmath201 be the set of ( almost ) dirac structures over a manifold @xmath21 .
is called `` daleth '' , which is the fourth letter of the semitic alphabet corresponding to `` d '' . ]
let @xmath202 and @xmath203 be distinct smooth manifolds .
to transition from systems on @xmath204 and @xmath205 to one on @xmath206 , we shall define a _
direct sum _ associated to distinct dirac structures .
let @xmath207 and @xmath208 be given distinct dirac structures on @xmath202 and @xmath203 respectively .
define the dirac sum of @xmath209 and @xmath210 by @xmath211    we now show that @xmath8 is a mapping from @xmath212 .    if @xmath213 , then @xmath214 .
let @xmath215 and let @xmath216 .
notice that @xmath217 for @xmath218 . by definition , @xmath219 and @xmath220 for some @xmath221 and @xmath222 such that @xmath223 and @xmath224 .
then , one has , for each @xmath225 , @xmath226 since @xmath223 and @xmath224 .
thus , noting that @xmath195 is arbitrary , one can prove the theorem by lemma [ lem : isotropic ] .
[ [ the - dirac - sum - of - induced - dirac - structures . ] ] the dirac sum of induced dirac structures .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    let @xmath227 and @xmath228 be distinct configuration spaces .
let @xmath229 and @xmath230 be constraint distributions , and we can define the induced dirac structures @xmath209 and @xmath210 , where we assume that @xmath231 and @xmath232 .
let @xmath233 be an extended configuration manifold and @xmath234 be a distribution on @xmath48 .
as before , define the induced distribution @xmath61 on @xmath50 by @xmath235 where @xmath54 is the tangent map of the cotangent bundle projection @xmath236 , while the annihilator of @xmath61 is defined by , for each @xmath96 , @xmath237    let @xmath238 be the canonical symplectic structures on @xmath239 and @xmath240 be the associated bundle maps .
let @xmath241 be the natural projections with @xmath242 . ] .
denote by @xmath243 the dirac structures on @xmath239 induced from @xmath244 .
then , the _ dirac sum _
@xmath245 can be represented by the induced dirac structure on @xmath50 from @xmath246 , which is given by , for each @xmath96 , @xmath247 where @xmath248 is the associated bundle map with @xmath249 .
let @xmath250 for each @xmath96 .
let @xmath251 .
now , we may decompose @xmath252 as @xmath253 and @xmath254 such that @xmath255 and @xmath256 .
define the distributions @xmath257 and @xmath258 .
then , it follows that @xmath259 and @xmath260 , where @xmath261 and @xmath262 . noting that @xmath263 and @xmath264 , one has @xmath265
. additionally , @xmath266 .
it follows on each fiber @xmath267 that @xmath268 and it is obvious that @xmath269 .
putting these all together gives us @xmath270 and @xmath271 , for each @xmath96 . since @xmath272 is the canonical symplectic form of @xmath50 we find @xmath245 is the induced dirac structure from @xmath165 .
[ [ an - interconnection - dirac - structure . ] ] an interconnection dirac structure .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    in order to intertwine distinct dirac structures , we will introduce a special dirac structure called _ interconnecting dirac structure_. a simple case may be given by a distribution and its annihilator .
suppose there exists a constraint distribution @xmath273 due to the interconnection , such that @xmath274 where @xmath275 is a one - form on @xmath48 associated with the interconnection .
the annihilator @xmath276 is defined by @xmath277 we can define the lifted distribution @xmath278 on @xmath50 by @xmath279 where @xmath54 is the tangent map of the cotangent bundle projection @xmath236 , while the annihilator of @xmath280 is given by , for each @xmath96 , @xmath281 then , the interconnection dirac structure associated to the distribution @xmath280 on @xmath50 can be defined by , for each @xmath96 , @xmath282    [ [ example - two - particles - moving - with - contact . ]
] example : two particles moving with contact .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    consider two masses moving with contact in a vector space @xmath283 , whose velocities are given by @xmath284 .
since the two particles are in contact and their velocities are common , it follows @xmath285 where @xmath286 is a constraint subspace of @xmath10 .
this constraint is enforced through the associated constraint forces @xmath287 at the contact point , where @xmath288 is the force associated with the @xmath289-th particle .
this implies newton s third law and it follows @xmath290 where @xmath291 is the annihilator of @xmath278 .
note that we can define a dirac structure on @xmath10 by @xmath292 .
then , it follows that @xmath293 .
[ [ remark.-3 ] ] remark .
+ + + + + + +    we can consider a more general interconnection dirac structure than stated in the above , such as those given by gyrators , transformers and so on , which are expressed by two - forms or poisson structures . in this paper , we will mainly focus on the simple case of the interconnection dirac structure @xmath1 induced from a given distribution @xmath278    [ [ strategy - for - the - interconnection - of - distinct - dirac - structures . ] ] strategy for the interconnection of distinct dirac structures .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    before going into details , let us outline our strategy for interconnecting dirac structures .
let @xmath209 and @xmath210 be two distinct dirac structures on distinct manifolds @xmath202 and @xmath203 .
we will interconnect @xmath209 and @xmath210 by introducing an _ interconnection _ dirac structure @xmath1 on @xmath294 , by which constraints due to the interconnection are imposed .
to do this , we first make the direct sum of dirac structures , namely @xmath245 , which is a single dirac structure on @xmath21 .
then , by introducing the _ bowtie product _ @xmath295 of dirac structures on @xmath21 , we define an interconnected dirac structure as @xmath296 . the bowtie product will be shown to make @xmath297 a commutative category with identity element @xmath298 . under the condition
that the fiber dimension of @xmath299 is constant , it will be shown that @xmath300 is indeed a dirac structure on @xmath21 .
[ [ the - bowtie - product . ] ] the bowtie product .
+ + + + + + + + + + + + + + + + + + +    let @xmath301 and let @xmath302 be a submanifold given as @xmath303 .
given dirac structures @xmath304 , @xmath305 on @xmath21 , define the direct product of dirac structures @xmath304 and @xmath305 by @xmath306 on @xmath307 .
recall from @xcite that one may pull back the dirac structure @xmath308 on @xmath307 to @xmath309 as @xmath310 which is a lagrangian submanifold of @xmath311 .
assume that @xmath312 is a smooth sub - bundle by either of the following conditions :    * @xmath313 has constant dimension , * @xmath314 has constant dimension .
the subbundle @xmath312 is a dirac structure on @xmath309 .
let @xmath315 be the map @xmath316 ( this is a diffeomorphism between @xmath21 and @xmath309 ) .
then @xmath317 given by @xmath318 is called the _ diagonal embedding_. given dirac structures @xmath319 , we may define the dirac structure @xmath320 .
noting @xmath321 , it follows that @xmath322 \in \daleth(m).\ ] ]    let @xmath304 and @xmath323 be dirac structures over @xmath21 .
then , we define a tensor product called the bowtie product of @xmath304 and @xmath323 by @xmath324 or alternatively given by : @xmath325    under the assumption that @xmath326 has constant dimension , @xmath327 is a dirac structure on @xmath21 .
[ [ remark.-4 ] ] remark .
+ + + + + + +    in this paper , we focus on how to use the bowtie product @xmath2 to make a dirac structure for an interconnected lagrange
dirac dynamical system .
while the bowtie product was initially proposed in @xcite for this purpose using formula , it was revealed that an equivalent _ tensor product _ ( denoted by @xmath328 ) was also developed in the context of generalized complex geometry ( see @xcite and @xcite ) using formula .
we will make explicit the equivalence of @xmath2 and @xmath328 ( i.e. and ) in the appendix .
[ [ remark.-5 ] ] remark .
+ + + + + + +    for standard lagrange - dirac mechanics , we usually employ @xmath329 , @xmath304 is the canonical dirac structure on @xmath21 , and @xmath330 for some distribution @xmath331 is an interconnection dirac structure . though we will not explore this in detail ,
one of the particularly interesting interconnection dirac structures is the case in which @xmath305 can be developed by lagrangians that are linear in velocities of the form @xmath332 , where @xmath333 and @xmath334 is a one - form on @xmath48 causing magnetic terms to appear .    in order to show interconnection of dirac structures
is associative , we will use a special restricted two - form @xmath32 induced by a dirac structure @xmath6 on @xmath21 with a distribution @xmath335 , where @xmath336 .
this two - form will be employed in the proofs for theorems [ thm : bowtie ] and [ associative ] .
[ lem : two - form ] let @xmath337 and set @xmath335 . on each fiber of @xmath338 at @xmath28 there exists a bilinear anti - symmetric map @xmath339 defined as @xmath340    this two - form was initially introduced by @xcite for linear dirac structures ( see also @xcite ) . we can easily generalize it to the case of general manifolds since @xmath32 is defined fiber - wise , and hence we have that @xmath26 is a linear dirac structures on each fiber @xmath341 for all @xmath28 . here , we present a proof for completeness .
we desire to show that equation defines a unique bi - linear anti - symmetric map @xmath32 .
first we prove linearity .
let @xmath342 and let @xmath343 . since @xmath26 is a linear space at each @xmath28 we have @xmath344 . by definition of @xmath32 , it follows , for all @xmath345 , @xmath346 thus , @xmath347 is linear in the first argument .
second , in order to prove that @xmath32 is anti - symmetric , take any @xmath348 .
since @xmath6 is isotropic , we find @xmath349 and we compute @xmath350 thus , we have shown anti - symmetry from which bi - linearity follows .
given a dirac structure @xmath337 , we call the two - form @xmath32 of lemma [ lem : two - form ] the dirac two - form .
recall from equation that , for each @xmath28 , @xmath26 may be given by @xmath351 which may be also stated by @xmath352 where @xmath353 is the skew - symmetric bundle map that is a natural extension of the skew - symmetric map @xmath354 , which is defined by @xmath355 on @xmath356 .
[ thm : bowtie ]    let @xmath304 and @xmath357 .
let @xmath358 and let @xmath359 and @xmath360 be the dirac two - forms for @xmath304 and @xmath305 respectively . if @xmath361 has constant rank , then @xmath362 is a dirac structure with distribution @xmath363 and with the dirac two - form @xmath364 .
[ pf : bowtie ] let @xmath365 for each @xmath28 .
then , there exists @xmath366 such that @xmath367 for each @xmath28 .
this is equivalent to @xmath368 where @xmath369 and @xmath370 .
this means @xmath371 and @xmath372 . but @xmath373 .
therefore , upon setting @xmath374 and @xmath375 , we can write @xmath376 and @xmath377 ; namely , @xmath378 , where @xmath379 is a dirac structure with @xmath278 and @xmath380 .
then , it follows that @xmath381 .
equality follows from noting that both @xmath382 and @xmath383 are subspaces of @xmath384 with the same dimension .
[ cor : bowtie2 ] if @xmath385 , then it follows that @xmath386 and also that @xmath387 is induced by the constraint @xmath388 and the two - form @xmath359 .    [ associative ]
let @xmath389 with distributions @xmath390 , @xmath391 , and @xmath392 .
assume that @xmath361 , @xmath393 and @xmath394 have constant ranks .
then the bowtie product @xmath2 is associative and commutative ; namely we have @xmath395 and @xmath396    first we prove commutativity .
recall that any dirac structure may be constructed by its associated constraint distribution @xmath335 and the dirac two - form @xmath347 .
let @xmath397,and @xmath398 be the dirac two - forms corresponding to @xmath399 , and @xmath400 respectively .
then we find by theorem [ thm : bowtie ] that @xmath401 is defined by the distribution @xmath402 and the dirac two - form @xmath403 . by commutativity of @xmath404 and @xmath405
, we find the same distribution and the two - form for @xmath406 , we have @xmath407 .
now we prove associativity .
let @xmath408 and @xmath409 and it follows @xmath410 if @xmath411 and @xmath412 are respectively the dirac two - forms for @xmath413 and @xmath414 , we find @xmath415|_{\delta_{(ab)c } } \\          & = ( \omega_{\delta_{a } } + \omega_{\delta_{b } } + \omega_{\delta_{c } } ) |_{\delta_{(ab)c } } \\          & = ( \omega_{\delta_{a } } + \omega_{\delta_{b } } + \omega_{\delta_{c } }   ) |_{\delta_{a(bc)}}\\          & = \omega_{\delta_{a(bc)}}.\end{aligned}\ ] ] thus , we obtain @xmath416    [ [ remark.-6 ] ] remark . + + + + + + +    we have shown that @xmath2 acts on pairs of dirac structures with clean intersections to give a new dirac structure and also that that it is an associative and commutative product .
it is easy to verify from equation that the dirac structure @xmath417 satisfies the property of the identity element as @xmath418 for every @xmath337 .
thus we find that @xmath297 is a commutative category ( it lacks totality and an inverse element , so it is neither a monoid nor a groupoid ) .
[ [ interconnection - of - two - distinct - induced - dirac - structures . ] ] interconnection of two distinct induced dirac structures .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    let @xmath419 and @xmath420 be dirac structures on distinct manifolds @xmath227 and @xmath228 induced from constraint distributions @xmath421 and @xmath422 as before .
let @xmath278 be a given distribution on @xmath423 due to the interconnection of @xmath209 and @xmath210 .
we have mathematical ingredients for _ interconnecting _ the induced dirac structures @xmath209 and @xmath210 through the interconnection dirac structure @xmath1 on @xmath48 .
recall that an interconnection dirac structure is given by @xmath424 , where @xmath425 is the constraint distribution associated to the interconnection .
recall that @xmath426 and @xmath427 .
it follows from theorem [ thm : bowtie ] that @xmath428 may be constructed from the constraint @xmath429 and @xmath430 , where @xmath431 and @xmath432 are respectively the canonical two - forms on @xmath433 and @xmath434 .
[ prop : velocity_constraint_bowtie ] assuming @xmath30 has constant rank , the interconnection of two distinct induced dirac structures @xmath209 and @xmath210 through @xmath435 , namely , @xmath436 is represented by , @xmath437    this may be viewed as a corollary to theorem [ thm : bowtie ] and proposition [ induced_ds_dirac_sum ] .    [
[ interconnection - of - n - distinct - dirac - structures . ] ] interconnection of @xmath183 distinct dirac structures .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    let us consider the interconnection of @xmath183 distinct dirac structures @xmath438 on distinct manifolds @xmath439 , in which @xmath440 .
recall that the dirac sum @xmath8 is associative as @xmath441 so that we may define the iterated direct sum by @xmath442 by choosing the appropriate interconnection dirac structure @xmath443 and assuming that @xmath444 is constant , we can simply develop the interconnected dirac structure for @xmath445 distinct dirac structures as @xmath446 the condition that @xmath447 is constant is equivalent to the assumption that @xmath448 is constant , which can be shown by resorting to local coordinate expressions .    [ [ a - link - between - composition - and - interconnection - of - dirac - structures . ] ] a link between composition and interconnection of dirac structures .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    the notion of _ composition _ of dirac structures was introduced in @xcite in the context of port - hamiltonian systems .
let @xmath449 and @xmath450 be vector spaces .
let @xmath419 be a linear dirac structure on @xmath451 and @xmath420 be a linear dirac stucture on @xmath452 .
the _ composition _ of @xmath419 and @xmath420 is given by @xmath453 where @xmath454 and @xmath455 denote the dual space of @xmath456 and @xmath457 .
it was also shown that the set @xmath458 is itself a dirac structure on @xmath459 . in this section
, we shall see the link between the composition and the interconnection of dirac structures .
the composition developed in @xcite is constructed on vector - spaces .
so , we also focus on the vector space case in this paper .
let @xmath460 and @xmath461 . before going into details
, we establish an important fact on the projection from @xmath10 to @xmath462 .
set @xmath463 to be the projection @xmath464 .
the mapping @xmath465 dual to @xmath466 is given by @xmath467    next , we recall the push - forward map associated to a dirac structure .
let @xmath64 and @xmath468 be vector spaces .
let @xmath469 be a linear map and let @xmath470 be a linear dirac structure on @xmath64 .
the push - forward of @xmath470 to @xmath468 by @xmath471 is the set @xmath472 it is worth noting that the push - forward of a dirac structure @xmath473 is itself a dirac structure @xmath474 ; namely , we write it as @xmath475 see , for instance , @xcite and @xcite .
we are now ready to provide the link between the bowtie product and the composition of dirac structures .
[ thm : composition ] let @xmath476 .
set @xmath477 .
set @xmath463 to be the projection @xmath464 . for linear dirac structures
@xmath419 on @xmath478 and @xmath420 on @xmath479 , set @xmath480 and set also @xmath481 then , one has @xmath482
first , we observe that @xmath483 .
we also observe @xmath484 utilizing the fact that @xmath485 if and only if @xmath486 and @xmath487 , we may restate the above as @xmath488 finally , we have @xmath489 but this is nothing but the composition of dirac structures , namely , @xmath490 .
given distinct mechanical systems , one can develop dynamics of each mechanical system in the context of both the lagrange - dalembert - pontryagin variational structures and the induced dirac structures over separate configuration manifolds . in this section , we will show how dynamics of an interconnected systems can be formulated in the context of the induced dirac structures as well as in the context of the variational structures .
[ [ variational - structures - for - interconnected - mechanical - systems . ] ] variational structures for interconnected mechanical systems .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    let us introduce an interconnection of the lagrange - dalembert - pontryagin principles for distinct mechanical systems .
let @xmath491 be distinct configuration manifolds corresponding to distinct mechanical systems @xmath492 .
let @xmath493 be the generalized energies associated to given lagrangians ( possibly degenerate ) @xmath494 on @xmath495 and @xmath496 be external force fields .
let @xmath497 be a distribution on @xmath491 for @xmath492 .
the lagrange - dalembert - pontryagin principles for curves @xmath498 in @xmath499 is denoted by @xmath500 for @xmath501 and with @xmath502 .
the interconnection of distinct lagrange - dalembert - pontryagin principles may be given by imposing some distribution @xmath278 on the extended configuration manifold @xmath503 such that @xmath504    it follows from equations and that , by using variations @xmath505 vanishing at the endpoints , one can obtain @xmath506    by direct computations and integration by parts , it follows that @xmath507 for @xmath508 and with @xmath509 . imposing the interconnection condition ,
we obtain implicit lagrange - dalembert equations for the interconnected system .
[ [ an - alternative - lagrange - dalembert - pontryagin - principle . ] ] an alternative lagrange - dalembert - pontryagin principle .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    here , we show that there exists an alternative lagrange - dalembert - pontryagin principle .
let @xmath503 and let @xmath510 , @xmath511 and @xmath512 .
define the generalized energy for the interconnected system by @xmath513 , where @xmath514 is the lagrangian for the interconnected system .
set the constraint distribution @xmath515 .
the lagrange - dalembert - pontryagin principle for a curve @xmath516 in @xmath133 may be given by @xmath517 for the chosen variation @xmath518 , for all @xmath178 and @xmath179 , and with @xmath519 .
+ keeping the endpoints of @xmath164 fixed , it follows that one can simply derive the lagrange - dalembert - pontryagin equations as follows : @xmath520    the lagrange - dalembert - pontryagin principles in for curves @xmath521 @xmath522 , @xmath523 in @xmath499 together with the interconnection condition in is equivalent with the lagrange - dalembert - pontryagin principle in .
it follows from the lagrange - dalembert - pontryagin principles that one can obtain equation as @xmath524 imposing the interconnection condition @xmath525 and @xmath526 , one has @xmath527 and @xmath528 where we note that the final distribution is given by @xmath529 and hence that its annihilator is denoted by @xmath530    the following statements are equivalent :    1 .
a curve @xmath531 in @xmath133 satisfies the lagrange
pontryagin principle : @xmath532 with respect to variations @xmath533 and @xmath534 arbitrary and with fixed endpoints of @xmath164 .
curves @xmath535 satisfy the lagrange  dalembert  pontryagin principles @xmath536 under variations @xmath537 , arbitrary @xmath538 and @xmath539 , together with @xmath502 and fixed endpoints of @xmath540 , such that the force vector @xmath541 and the velocity vector @xmath542 satisfy @xmath543 where @xmath544 denote force fields .
[ [ interconnected - lagrange - dirac - dynamical - systems . ] ] interconnected lagrange - dirac dynamical systems . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    next , consider the interconnection of distinct lagrange - dirac dynamical systems corresponding to the above variational structures . to do this
, we first consider the interconnection of distinct dirac structures @xmath243 on @xmath545 , each of which is induced from given distinct distributions @xmath546 .
let @xmath547 be @xmath183 distinct lagrange - dirac systems , where @xmath548 denote the generalized energies on @xmath499 associated to given lagrangians @xmath494 on @xmath495 , @xmath549 are the partial vector fields , and the external force fields @xmath550 .
recall that @xmath551 induce the horizontal one - forms @xmath552 on @xmath499 such that @xmath553 for all @xmath554 , where @xmath555 are natural projections .
as before , for an interconnected system , let @xmath503 and @xmath278 be a given distribution on @xmath48 associated to the interconnection .
let @xmath515 .
then , define an interconnection dirac structure by @xmath556 as in , where @xmath557 . set a generalized energy @xmath513 associated to the lagrangian @xmath558 .
furthermore , define a partial vector field on @xmath133 by @xmath559 define also an external force field by @xmath560 which induces the horizontal one - form @xmath561 on @xmath133 such that @xmath562 for all @xmath563 , where @xmath564 is the natural projection . recall from that the interconnection of distinct dirac structures @xmath565 through @xmath1
is given by @xmath566    the interconnected lagrange - dirac dynamical system is a quadruple @xmath103 that satisfies @xmath567 where @xmath568 .
using local coordinates @xmath510 , @xmath569 , and @xmath570 , one has @xmath571    [ prop : ldds_splitting ] the following statements are equivalent :    1 .
a curve @xmath572 satisfies @xmath573 where @xmath574 .
2 .   curves @xmath575 satisfy @xmath576 where @xmath577 , such that @xmath578 where @xmath579 .
recall that @xmath580 , the condition @xmath581 is equivalent to the existence of some covector @xmath582 such that @xmath583 and @xmath584 it follows from the condition that @xmath585 and also from condition it follows that @xmath586 and @xmath587 .
furthermore , noting @xmath588 , it follows that the curves @xmath589 satisfy the conditions @xmath590 together with @xmath591 .
one can easily check the converse .
thus , we obtain the following theorem :    [ thm : ldds_splitting ] the following statements are equivalent :    1 .
a curve @xmath592 in @xmath593 satisfies @xmath594 where @xmath595 .
2 .   curves @xmath596 satisfy @xmath597 where @xmath598 and @xmath599 is the horizontal one - forms on @xmath600 of force fields @xmath601 , such that the condition of interconnection holds @xmath602 3 .
a curve @xmath592 satisfies the hamilton
pontryagin principle : @xmath603 with respect to variations @xmath604 and @xmath534 arbitrary , with @xmath605 and with fixed endpoints of @xmath164 .
curves @xmath596 satisfy the lagrange  dalembert  pontryagin principles : @xmath606 where @xmath501 , and @xmath607 are arbitrary , @xmath502 and with fixed endpoints of @xmath164 , such that the interconnection condition for the force vector @xmath541 and the velocity vector @xmath542 holds as @xmath608
in this section we provide several examples of how interconnection can be carried out with the bowtie product .
we have chosen rather simple examples , however the advantage of using the bowtie product is based on the procedure of tearing and interconnecting systems , which enables us to treat more complicated systems in a hierarchical structure . in future , we will seek more complex uses of the bowtie product .    *
( i ) example : mass - spring mechanical systems . *
+ let us consider an illustrative example of a mass - spring system in the context of the interconnection of dirac structures and associated lagrange - dirac systems .
in this example , there exists _ no external force_. let @xmath609 and @xmath610 be the @xmath289-th mass and spring @xmath611 .        [ [ tearing - and - interconnection . ] ] tearing and interconnection .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + +    inspired by the concept of _ tearing and interconnection _ originally developed by @xcite , the mechanical system can be torn apart into two subsystems as in fig [ tornapart ] , each of which can be regarded as a modular unit of the whole system .
let us call each disconnected system a `` primitive system '' by following @xcite .
note that the procedure of tearing inevitably yields _ interactive boundaries _ associated to the two disconnected subsystems , where it induces the following condition : @xmath612 without the above condition , there exists no energetic interaction between the disconnected systems .
in other words , _ the original mechanical system can be reconstructed by interconnecting the subsystems through the continuity conditions_.    the continuity conditions in imply the continuity of power flow ; namely , the _ power invariance _ holds as @xmath613 where @xmath614 and @xmath615 needless to say , equation can be understood as the defining condition for an interconnection dirac structure .    [ [ primitive - lagrangian - systems . ] ] primitive lagrangian systems .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + +    let us consider how dynamics of the disconnected subsystems can be formulated in the context of lagrange - dirac dynamical systems .
the configuration space of the subsystem 1 may be given by @xmath616 with local coordinates @xmath617 , while the configuration space of the subsystem 2 is @xmath618 with local coordinates @xmath619 .
we can naturally define the canonical dirac structures @xmath620 and @xmath621 in this example . for primitive system 1 ,
the lagrangian @xmath622 is given by @xmath623 while the lagrangian @xmath624 for primitive system 2 is given by @xmath625 then , we can define the generalized energy @xmath626 on @xmath627 by @xmath628 and also define the generalized energy @xmath629 on @xmath630 by @xmath631 .
though the original system has no external force , each disconnected primitive system has an interconnection constraint force at the interactive boundary . when viewing each system separately , the constraint force acts as an external force .
this is because _ tearing always yields constraint forces at the boundaries associated with the disconnected primitive systems _ , as shown in fig .
[ tornapart ]    the constraint force fields at the boundaries @xmath632 and @xmath633 are given by @xmath634 an @xmath635 .
further , let @xmath636 be the partial vector field , which is defined at each point @xmath637 as @xmath638 similarly , let @xmath639 be the partial vector field , which is defined at each point @xmath640 as @xmath641 together with the _ consistency condition _ @xmath642 .
primitive system 1 : we can formulate dynamics of primitive system 1 in the context of the lagrange - dirac dynamical system @xmath643 as @xmath644 which may be given in coordinates by @xmath645 and with @xmath646 and @xmath647 .    primitive system 2 : similarly , we can also formulate dynamics of primitive system 2 in the context of the lagrange - dirac dynamical system @xmath648 as @xmath649 which may be given in coordinates by @xmath650 together with @xmath651 and @xmath652 as well as @xmath642 .    recall that each primitive system is physically independent of the other , which implies that there exists no energetic interaction between them .
next , let us see how the disconnected primitive systems can be interconnected through a dirac structure .    [ [ interconnection - of - distinct - dirac - structures . ] ] interconnection of distinct dirac structures .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    let @xmath653 be an extended configuration space with local coordinates @xmath654 .
recall that the dirac sum of the dirac structures is given by @xmath245 on @xmath50 .
now , the constraint distribution due to the interconnection is given by @xmath655 where @xmath656 is a one - form on @xmath48 .
on the other hand , the annihilator @xmath657 is defined by @xmath658 it follows from this codistribution that @xmath659 , @xmath660 and @xmath661 .
hence , we obtain the conditions for the interconnection as in equation ; namely , @xmath662 and @xmath663 .
recall that the distribution @xmath280 on @xmath50 is obtained by @xmath664 and hence the dirac structure @xmath665 can be defined as in .
further , recall that the interconnected dirac structure @xmath6 on @xmath50 is given by @xmath666    [ [ interconnected - lagrange - dirac - systems . ] ] interconnected lagrange - dirac systems .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    now , let us see how two primitive lagrange - dirac systems , namely , @xmath667 and @xmath668 can be interconnected to be a lagrange - dirac dynamical system .
define the lagrangian @xmath669 for the interconnected system by @xmath670 and hence the generalized energy is given by @xmath671 .
. set a partial vector field by @xmath673 , which is defined at points @xmath674 , where the consistency condition @xmath642 holds .
finally , the interconnected lagrange - dirac system is given by @xmath675 , which satisfies , for each @xmath676 , @xmath677 and with @xmath678 .
[ [ the - lagrange - dalembert - pontryagin - principle.-1 ] ] the lagrange - dalembert - pontryagin principle .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    alternatively , we can develop the interconnected lagrange - dirac system by the lagrange - dalembert - pontryagin principle : @xmath679 for all @xmath680 , for all @xmath681 and @xmath682 , and with @xmath683 .
thus , we can obtain the dynamics of the interconnected lagrange - dirac system . *
( ii ) example : electric circuits * + consider the electric circuit depicted in figure [ vrlc_circuit_pic ] , where @xmath684 denotes a resistor , @xmath172 an inductor , and @xmath685 a capacitor .        as in figure [ torn_circuit_pic ] , we decompose the circuit into the disconnected two circuits , each of which we shall call `` primitive circuit 1 '' and `` primitive circuit 2 '' , in which @xmath686 and @xmath687 denote external ports . in order to reconstruct the original circuit in figure [ vrlc_circuit_pic ] ,
the external ports may be connected by equating currents across each .
[ [ primitive - circuit-1 . ] ] primitive circuit 1 .
+ + + + + + + + + + + + + + + + + + + +    the configuration manifold for circuit 1 is denoted by @xmath688 with local coordinates @xmath689 , where @xmath690 and @xmath691 are the charge fluxes through the resistor @xmath684 , inductor @xmath172 and port @xmath686 respectively .
there exists the kcl constraint distribution @xmath229 in circuit 1 , which is given by for each @xmath692 , @xmath693 where @xmath694 denotes the current vector at each @xmath695 , while the kvl constraint is given by its annihilator @xmath696 .
then , we can naturally define the induced dirac structure @xmath419 on @xmath697 from @xmath698 as before .
the lagrangian for circuit 1 , namely , @xmath699 on @xmath700 is given by @xmath701 which is degenerate .
define the generalized energy @xmath702 on @xmath627 by @xmath703 on @xmath627 .
circuit 1 also has the external force field due to the resistor and the external port , @xmath704 as @xmath705 so , we can set up the lagrange - dirac dynamical system @xmath706 .
[ [ primitive - circuit-2 . ] ] primitive circuit 2 .
+ + + + + + + + + + + + + + + + + + + +    the configuration manifold for circuit 2 is @xmath707 with local coordinates @xmath708 , where @xmath709 is the charge flux through the port @xmath710 and @xmath711 is the charge stored in the capacitor .
the kcl constraint distribution is given by , for each @xmath712 , @xmath713 and so the kvl space is given by its annihilator @xmath714 .
this gives us the dirac structure @xmath210 on @xmath715 .
set the lagrangian @xmath716 for circuit 2 as @xmath717 and set the generalized energy @xmath718 on @xmath630 .
circuit 2 has the external force field due to the resistor and the external port @xmath687 , namely , @xmath719 as @xmath720 then , we can formulate the lagrange - dirac dynamical system @xmath721 .
[ [ the - interconnection - dirac - structure ] ] the interconnection dirac structure + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    set @xmath722 . the interconnection constraint is given by @xmath723 by the tangent lift @xmath724 , we can define the interconnection dirac structure as @xmath725 which is denoted , in local coordinates , by @xmath726    [ [ the - interconnected - system . ] ] the interconnected system .
+ + + + + + + + + + + + + + + + + + + + + + + + + +    now , we can define the interconnected dirac structure by @xmath727 which may be viewed as the induced dirac structure from the constraint space @xmath728 which is given , in coordinates @xmath729 , by @xmath730    set the lagrangian for the interconnected system as @xmath731 and the external force field @xmath732 .
set also @xmath733 .
let @xmath734 be a partial vector field .
the interconnected lagrange - dirac dynamical system is given by the quadruple @xmath735 , which satisfies @xmath736 where @xmath737 .
* ( iii ) example : a rolling ball on rotating tables * + consider the mechanical system depicted in fig .
[ fig : gear_system ] , where there are rotating tables via a gear and a ball is rolling on one of the tables without slipping .
we assume that there is no loss of energy and non - slipping in the gear and also that the external torque is constant .
let @xmath738 be moments of inertia for the table @xmath739 .
let us decompose the system into distinct three systems ; ( 1 ) a rotating ( small ) table , ( 2 ) a rotating ( large ) table , and ( 3 ) a rolling ball on a table .
[ [ primitive - system-1 . ] ] primitive system 1 .
+ + + + + + + + + + + + + + + + + + +    for primitive system 1 the configuration manifold is the circle , @xmath740 .
the lagrangian is the rotational energy of the system . @xmath741 and set the generalized energy @xmath742 on @xmath627 .
we employ the canonical dirac structure on @xmath697 as @xmath743 there is a constant torque acting on primitive system 1 , given by a map @xmath744 defined as @xmath745 , where @xmath746 is constant in @xmath747 .    [ [ primitive - system-2 . ] ] primitive system 2 .
+ + + + + + + + + + + + + + + + + + +    for primitive system 2 , the configuration manifold is also the circle , @xmath748 and the lagrangian is again the rotational energy @xmath749 and set the generalized energy @xmath750 on @xmath630 .
again , we have the canonical dirac structure @xmath751    [ [ primitive - system-3 . ] ] primitive system 3 .
+ + + + + + + + + + + + + + + + + + +    primitive system 3 is a rolling sphere of uniform density and radius @xmath752.the sphere moves in space by changing its position and orientation relative to a reference configuration .
the configuration manifold is given by the special euclidean group @xmath753 , which we parametrize as @xmath754 where @xmath755 .
following @xcite , let @xmath756 be the set of points of the sphere in the reference configuration . for configuration @xmath757 , a point @xmath758
is transformed into @xmath759 by the action @xmath760 .
the lagrangian is given by the kinetic energy as @xmath761 where @xmath762 .
we use body coordinates such that the center of the sphere in the reference configuration is at the origin so that @xmath763 .
substituting these relations , the above lagrangian is to be @xmath764 setting @xmath765 and noting that @xmath766 when @xmath767 , one finally obtains @xmath768 and with the generalized energy @xmath769 on @xmath770 .
since the motion along the @xmath771-direction is constrained so that the ball does not leave the plane of table 2 , we get the holonomic constraint @xmath772 this yields the induced dirac structure @xmath773    [ [ interconnection - constraints . ] ] interconnection constraints .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + +    let @xmath774 .
in order to interconnection the three primitive systems , we need to impose the constraints due to the non slip conditions .
first , consider the mapping @xmath775 for @xmath776 .
the map @xmath777 sends the tangent fiber @xmath778 to the lie - algebra @xmath779 , where @xmath780 is the identity element on @xmath781 viewed as a group . from this point on ,
we do not explicitly invoke @xmath777 but use it to interpret @xmath782 as @xmath783 . the interconnection constraint between system 1 and system 2
is given by @xmath784 and with its annihilator @xmath785 where @xmath786 .
this constraint ensures that the gears rotate ( without slipping ) at the same speed in opposite directions .
next , we consider the interconnection constraint between systems 2 and 3 .
note that the velocity of a point located at the bottom of the sphere is given by @xmath787 note also that a point rotating on the gear of system 2 with the axle taken to be the origin has velocity @xmath788 so the non slip condition between system 2 and 3 is given by @xmath789 where @xmath790 are the basis on @xmath759 . setting the interconnection constraint distribution @xmath791
let @xmath792 and with its annihilator @xmath793 .
then , the interconnection dirac structure is given by @xmath794    [ [ the - interconnected - lagrange - dirac - system . ] ] the interconnected lagrange - dirac system .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    the dirac structure for the interconnected system is given by @xmath795 and we use the lagrangian @xmath796 .
the dynamics of the system are implied by the constraint : @xmath797 where @xmath798 is a partial vector field , and @xmath101 is the generalized energy of @xmath172 .
this gives us the interconnected lagrangian system @xmath103 . by theorem [ thm : bowtie ] , these dynamics use the constraint distribution @xmath799 and the canonical two - form on @xmath50 restricted to @xmath30 as the dirac two - form .
noting the annihilator is given by @xmath800 we finally obtain the dynamics as @xmath801
in this paper , we provided a notion of interconnection for dirac structures and their associated lagrange - dirac dynamical systems .
the principal idea lies in the fact that the _ interconnection of distinct dirac structures itself is a dirac structure_. we also clarified the interconnection of the associated lagrange - dirac dynamical systems with the variational structure by the lagrange - dalembert - pontryagin principle . to do this
, we developed the bowtie product , @xmath2 , which serves us well in the case of interconnection by constraints .
we then showed how the bowtie product can effectively carry the lagrange - dalembert - pontryagin principle to the interconnected system .
lastly , we demonstrated our theory with some illustrative examples of a mass - spring mechanical system , an interconnection in circuits , and a non - holonomic mechanical system .
we raise the following topics as arena s for further investigation :    * applications to complicated systems such as guiding central motion problems , multibody systems , fluid - rigid body interactions , passivity controlled interconnected systems , etc .
( see , for example , @xcite and @xcite ) .
* reduction and symmetry for interconnected lagrange - dirac systems ( @xcite , @xcite ) . * the use of more general interconnection dirac structures with nontrivial dirac two - forms such as the one given with symplectic gyrators ( in this paper , we restrict @xmath1 to the form @xmath802 ) . as to examples of such general interconnections in physical systems , see for instance @xcite and @xcite . *
the integrability condition for the bowtie product of dirac structures . as to the integrability condition for dirac structures , see @xcite and @xcite .
the link with _
symplectic categories _ ( see @xcite ) . * discrete version of interconnected lagrange - dirac mechanics ; namely , discretizing the hamilton - pontryagin principle one arrives at a discrete mechanics version of dirac structures ( see @xcite and @xcite ) .
we will show the equivalence between the bowtie product @xmath2 and the dirac tensor product @xmath328 developed by @xcite , where the definition of the tensor product @xmath328 is given as follows .
let @xmath803 .
let @xmath804 be the diagonal embedding in @xmath805 .
the tensor product of @xmath304 and @xmath305 is defined as @xmath806 where @xmath807 and @xmath808 is the orthogonal complement in @xmath809 .
first note that @xmath812 so that @xmath813
. then we see that @xmath814 and the equivalence class @xmath815 is @xmath816 where we chose the element @xmath817 as the representative for the equivalence class @xmath818 = \ { ( v , v,\alpha + \xi,\alpha-\xi ) , \forall \xi \in t^{\ast}p \ }   \in d_a \boxtimes d_b.\ ] ] splitting the expression @xmath306 into @xmath819 and @xmath323 , we obtain @xmath820 we see that the assignment @xmath821 is a bijection and also that @xmath822 this is just the definition of @xmath401 .
150 afshari , e. , h.  s. bhat , a. hajimiri , and j.  e. marsden [ 2006 ] , extremely wideband signal shaping using one- and two - dimensional nonuniform nonlinear transmission lines , _ journal of applied physics _ , * 99 * , 05401116 .
bloch , a.  m. and p.  e. crouch [ 1997 ] , representations of dirac structures on vector spaces and nonlinear l - c circuits .
differential geometry and control ( boulder , co , 1997 ) , _ proc .
pure math . _ * 64*. amer .
soc . , providence , ri , pages 103117 .
bloch , a.  m. [ 2003 ] , _ nonholonomic mechanics and control _ , volume  24 of _ interdisciplinary applied mathematics_. springer - verlag , new york . with the collaboration of j. baillieul , p. crouch and j. marsden , and with scientific input from p. s. krishnaprasad , r. m. murray and d. zenkov .
chang , d. , a.  m. bloch , n.  leonard , j.  e. marsden , and c.  woolsey [ 2002 ] , the equivalence of controlled lagrangian and controlled hamiltonian systems , _ control and the calculus of variations ( special issue dedicated to j.l .
lions ) _ , * 8 * , 393422 .
courant , t. and a.  weinstein [ 1988 ] , beyond poisson structures . in _
action hamiltoniennes de groupes .
troisime thorme de lie ( lyon , 1986 ) _ , volume  27 of _ travaux en cours _ ,
pages 3949 .
hermann , paris .
ortega ,  r. ,  j.a.l .
perez ,  p. j. nicklasson ,  h. j. sira - ramirez and h. sira - ramirez [ 1998 ] , _ passivity - based control of euler - lagrange systems : mechanical , electrical , and electromechanical applications_. communications and control engineering series , springer  verlag .
ortega , r. , a. j. van der schaft , b. maschke , and g.  escobar [ 2002 ] , interconnection and damping assignment passivity - based control of port - controlled hamiltonian systems , _ automatica _ , * 38*(4 ) , 585596 .
yoshimura , h. and j.  e. marsden [ 2006c ] , dirac structures and implicit lagrangian systems in electric networks , proc . of the 17th international symposium on mathematical theory of networks and systems , paper wea08.5 , 2006 , pp .
16 , july 24 - 28 , 2006 , kyoto .
yoshimura , h. and j.  e. marsden [ 2007a ] , dirac structures and the legendre transformation for implicit lagrangian and hamiltonian systems , _ lagrangian and hamiltonian methods in nonlinear control 2006 _ , volume * 366 * of _ lecture notes in control and information science series _ , pages 233247 , springer - verlag .
yoshimura , h. and j.  e. marsden [ 2008 ] , representations of dirac structures and implicit port - controlled lagrangian systems , _
proc.of international symposium on mathematical theory of networks and systems 2008 , blacksburg , virginia _ , pages 112 .
yoshimura , h. , h. jacobs and j.  e.  marsden [ 2010 ] , _ interconnection of dirac structures and lagrange - dirac dynamical systems _
, _ proceedings of mathematical theory on networks and systems in 2010 , budapest _ , july 59 . pages 16 .