Acessibilidade / Reportar erro

Processing and transmission of timing signals in synchronous networks

Abstract

In order to have accurate operation, synchronous telecommunication networks need a reliable time basis signal extracted from the line data stream in each node. When the nodes are synchronized, routing and detection can be performed, guaranteeing the correct sequence of information distribution among the several users of a transmission trunk. Consequently, an auxiliary network is created inside the main network, a sub-network, dedicated to the distribution of the clock signals. There are different solutions for the architecture of the time distribution sub-network and choosing one of them depends on cost, precision, reliability and operational security. In this work we analyze the possible time distribution networks and formulate problems related to precision and stability of the timing signals by using the qualitative theory of differential equations. Correspondences between constitutive parameters of the networks and the dynamics of the spatial phase and frequency errors are established.

bifurcation; master-slave network; phase-locked loop; synchronous network


Processing and transmission of timing signals in synchronous networks

J.R.C. Piqueira* * Sponsored by CNPq of Brazil. ; E.Y. Takada

Escola Politécnica da Universidade de São Paulo, Departamento de Engenharia de Telecomunicações e Controle, Av. Prof. Luciano Gualberto, travessa 3, n. 158, 05508-900 São Paulo, SP, Brasil, E-mails: piqueira@lac.usp.br

ABSTRACT

In order to have accurate operation, synchronous telecommunication networks need a reliable time basis signal extracted from the line data stream in each node. When the nodes are synchronized, routing and detection can be performed, guaranteeing the correct sequence of information distribution among the several users of a transmission trunk. Consequently, an auxiliary network is created inside the main network, a sub-network, dedicated to the distribution of the clock signals.

There are different solutions for the architecture of the time distribution sub-network and choosing one of them depends on cost, precision, reliability and operational security.

In this work we analyze the possible time distribution networks and formulate problems related to precision and stability of the timing signals by using the qualitative theory of differential equations. Correspondences between constitutive parameters of the networks and the dynamics of the spatial phase and frequency errors are established.

Mathematical subject classification: 70K20, 70K45.

Key words: bifurcation, master-slave network, phase-locked loop, synchronous network.

1 Introduction

The analysis of geographically separated oscillators started to become an important problem for telecommunications in the sixties with the introduction of the first digital trunks which required synchronous time basis for demodulation and regeneration of pulse code modulation (PCM) signals [11, 14].

The phase-locked loop (PLL) is a device introduced by Belescize [3] in 1932 to extract timing signals. Nowadays, it is used in integrated circuit versions with high precision and low cost [4]. This device can extract the clock from digital signals corrupted by distortion and noise in transmission media. They were initially used in regenerators and termination units of digital multiplexing equipment [19, 17].

With the development of higher hierarchy multiplexing systems, mechanisms to guarantee synchronization between digital streams of lower hierarchy were required in the terminal stations. Then two different synchronization strategies were developed: PDH (Plesiochronous Digital Hierarchy) e SDH (Synchronous Digital Hierarchy) [2, 23].

Nodes in PDH systems operate with precise and independent clocks, corrected by operators from time to time. This is an expensive network as it requires precise oscillators in all nodes. With cheaper oscillators, the operational result is unsatisfactory implying bad performance.

In SDH systems, synchronization between nodes can be achieved with a few nodes with precise clocks. The others use PLLs for extracting the clock signal from the line with good precision and low cost.

In this work we are interested only in synchronous networks. We are going to discuss the several possible solutions for the architecture of the clock distribution network taking the dynamics of the PLL as the basis of our analysis.

The idea is to show that, in spite of the problem complexity, using Dynamical Systems theory is an interesting tool in order to obtain conditions for existence and stability of synchronous states.

2 Phase-locking problem

The problem of phase-locking consists of controlling the phase of a local oscillator by the phase of an external oscillation, making them coincide or, at least, differ by a constant. From the point of view of electronic engineering, PLL is the device that accomplishes it. It is a closed loop system connecting three basic elements: a phase detector (PD), a filter (F) and a voltage-controlled oscillator (VCO) [4, 10]. A basic PLL is shown in figure 1.


The input and output signals are, respectively, given by:

vi(t) = Vi sin(w0t + qi(t)),

v0(t) = V0 cos(w0t + q0(t))

In these expressions, w0 is the central frequency here named free-running frequency of the loop, qi(t) and q0(t) are the instantaneous phases, and Vi and V0 are the amplitudes of vi(t) and v0(t).

We consider the loop in a locked or synchronous state when it reaches an equilibrium state, with constant phase error j = qi – q0 and null frequency error = - [4, 10].

As the phase detector is a signal multiplier, the PD output is given by:

where Km is the phase detector gain.

The filter is supposed to eliminate high frequency terms. So, if the double frequency term is sufficiently attenuated by the filter [10], equation (1) is reduced to:

with Kd = KmViV0, in volts per radian.

Here we consider the model of PD and its output vd(t) given by equation (2). Being simple, we take the filter F as an all-pole low-pass with zeros in infinite [20, 5] and transfer function:

where Vc(s) and Vd(s) represent the Laplace transforms of signals vc(t) and vd(t), respectively.

The combination of equations (2) and (3) yields:

The output phase of VCO q0 is controlled by vc(t) and satisfies = K0vc, where K0 is a VCO constant, in radians per volt per second [10]. Thus, equation (4) can be rewritten as:

Defining L(·)

and by taking the phase error j(t) = qi – q0 as the dynamic variable, equation (5) becomes:

The ordinary differential equation (6) describes the behaviour of a PLL that is the main component of circuits for extracting time signals.

3 Distribution of timing signals

The problem of time distribution along networks consists of controlling frequency and phase of clock signals spreading over a wide area. The idea is synchronizing the frequency and phase scales of several oscillators in a network by using the data communication capacity of the links.

This problem has several applications [16]:

  • Establishing a world wide time distribution system;

  • Synchronizing clocks located at different points in a digital communication network;

  • Distributing time signals in a network in order to apply control actions and commands at specific times;

  • Establishing a supercomputer by interconnecting several computers in a network.

These items are sufficient to justify the relevance of timing distribution in applications related to control and communication engineering.

In real problems, objective comparisons among the several possibilities are needed. Then, a precise mathematical treatment is necessary.

3.1 Problem formulation

As we have already stated, our intention is to discuss the several strategies for spreading clock signals and the synchronization of several oscillators distributed over a wide geographic area.

There are situations in which precision in synchronization is not a critical point. In these cases, independent clocks manually adjusted are used. This strategy originated the plesiochronous networks.

When synchronization results from interactions between the oscillators of the network we say that the network is synchronous.

Synchronous networks with a clock priority mechanism are called master-slave networks. When all the clocks in a network have equal relevance in determining the synchronous state, we say that the network is mutually synchronized.

Master-slave and mutually synchronized networks may include delay compensation techniques. Figure 2 shows a schematic diagram with the strategies of clock distribution [16].


In what follows, the phases of local oscillators, denoted by F, are composed by a free-running term wt, a forcing term q(t) and a perturbation P(t), i.e.,

F(t) = wt + q(t) + P(t).

3.2 Master-slave network classification

Master-slave networks are classified according to the transmission direction of the time basis in One-Way Master-Slave (OWMS) and Two-Way Master-Slave (TWMS).

In OWMS networks, the master clock has its own and independent time basis. Slave clocks have their basis depending on a unique node, the master or another slave. Besides, these networks are classified according to the topology in chain and star.

In TWMS networks, the master clock has its own time basis but the control signal sent to the slave clocks is adjusted according to the basis of other nodes. Slave clocks may have their time basis dependent on several nodes.

According to the topology, TWMS networks can be classified as chain, star or loop.

The possible strategies for master-slave networks are detailed in figure 3.


Master-slave networks are extensively adopted in public telecommunication networks due to simple implementation, good timing performance, reliability, and low cost [6]. They also have applications in parallel distributed computation [24], robotics [15], and multimedia systems [25].

3.2.1 Single-star OWMS network

The topology of a single-star OWMS network is illustrated in figure 4. Master-node, denoted by M, sends its time basis for all other nodes.


Master clock signal is independent on the other nodes. Its frequency is given by

where wM is the frequency of normal operation of the master clock, and Wt represents the deviation during the operation.

Slave-nodes are PLLs with input signal phase equal to the master-node phase, delayed by the propagation time from master to the considered slave node.

3.2.2 Single-chain OWMS network

The topology of this network is shown in figure 5. The master clock, denoted by M, sends its time basis to node-2, which sends its time basis to node-3, and so on, up to the last node.


Master-node, in this case, operates according to equation 7. Each slave clock can be considered a PLL.

As we have seen, the input signal phase in a node will be equal to the phase of the former node VCO, delayed by the propagation time.

3.2.3 Double-star TWMS network

This topology is illustrated in figure 6, with a master-node controlling the time basis of all slave nodes.


Master node M has an accurate and independent time basis. However the control signal that it sends to the slaves considers its own phase and the phase of all slaves.

3.2.4 Double-chain TWMS network

Figure 7 illustrates the topology of a double-chain TWMS network, with node i – 1 working as a master for node i. When establishing its time basis, each slave uses the phase of its master and slave nodes.


The master node M has an accurate and independent time basis. The control signal generated by the master M and sent to node 2 considers the phases of nodes M and 2.

3.2.5 Single-loop TWMS network

In this topology, node i - 1 works as a master for node i, as shown in figure 8.


The node M has an accurate and independent time basis. The control signal generated by the master M and sent to node 2, considers the phases obtained from M and from slave N.

3.2.6 Double-loop TWMS network

The topology of this network is shown in figure 9, with node i - 1 working as a master for node i.


A slave-node establishes its time basis by considering the signals from its master and slave nodes.

The master node M has an accurate and independent time basis. The control signal generated by the master M and sent to nodes 2 and N, considers the phases obtained from M and from slaves 2 and N.

3.3 Master node in TWMS networks

Figure 10 shows a scheme of master nodes in TWMS networks, indicating the mechanism for generating control signals considering the phase of the master FM and the phase of the slaves Fi.


Control signals F sent by the master to the network is submitted to a weighting process that considers all the phase of the slaves with coefficients ai,j such that

According to the network strategy the coefficients ai,j are:

double chain:
double-star:
simple-loop:
double-loop: a1,2 = 1 and a1,j = 0, Ñj ¹ 2
a1,j = 1/(N – 1), Ñj = 2,..., N
a1,N = 1 and a1,j 0, Ñj ¹ N
a1,2 = a1,N = 1/2 and a1,j = Ñj = 3,..., N – 1

3.4 Slave node in TWMS networks

In a TWMS network, PLL belonging to the i-slave node has an input signal with phase F resulting from a linear combination of phases from the several nodes, as shown in figure 11.


The linear combination follows the condition

in each ith-slave.

According to the clock distribution strategy, we have:

double chain: ai,i+1 = ai,i–1 = 1/2, Ñi = 2,...,N – 1 and
aN,1 = aN,N–1 = 1/2 double-star:
simple-loop:
double-loop: ai,1 = 1 and ai,j = Ñi, j = 2,..., N
ai,i–1 = 1
same as double chain

4 Center Manifold Theorem applied to the comparison between OWMS topologies

As explained in the last section, there are a lot of possible solutions regarding to the architecture of a time distribution network topology.

When designing the network it is important to know the robustness of the synchronous state depending on the constitutive parameters of the nodes.

Reasoning in this way, the Center Manifold Theorem from Dynamical Systems Theory is useful to study the behaviour of solutions near equilibrium states providing important hints about the synchronous state stability.

As an example, in this section we compare two OWMS topologies, single-star and single-chain, by applying the Center Manifold Theorem [18].

The nodes of the network are supposed to have second order PLLs as timing detectors. Slow instabilities and Doppler effect of master signal are neglected in our modeling.

4.1 Single-star OWMS network

The network to be analyzed consists of a master and two slaves in a single-star OWMS network.

Node 1 is the master with a free oscillator with phase F1(t). Nodes 2 and 3 are the slaves and they are second order PLLs with µ1 being the filter time constant and µ2, the loop gain.

Phases of slaves, F2(t) and F3(t), respectively, are controlled by the master phase signal. Transmission delays ti,1 from master node to ith-node, for simplicity, will be considered equal t2,1 = t3,1 = t.

Then, the network dynamics is given by:

The subscript j stands for the master, and i for the slaves. Consequently, jii is the difference between the phase of the VCO output signal and the phase of the PD input signal in the ith-slave and is called local phase error.

In order to represent the spatial phase error we use jij that is the difference between the phase of the VCO output signal of node i and the phase of the master.

By considering that the phase of the master output signal is a step at t=0, we have:

1(t) = img1(t – t) = 0 and 1(t) = 1(t – t) = 0, considering t > t.

Re-scaling the time variable T = µ1t, equation (8) becomes:

where µ = µ2/ and x' = dx/dT.

Equations (9) show that the dynamics of the interaction between nodes 1 and 2 and between nodes 1 and 3 are identical and there are two pairs of non coupled differential equations.

Without loss of generality, we are going to study the equations related only to the first pair of nodes.

By choosing the state variables as x1 = , x2 = j22 and x3 = :

This system admits a cylindrical phase surface, so equilibrium states of (10) are P1 = (0, 0, 0) e P2 = (0, –p, 0), and the eigenvalues of the Jacobian matrix associated, calculated in P1, are:

By observing the eigenvalues, we conclude that:

  • for µ > 0: dim(Es) = 3, P1 is asymptotically stable.

  • for µ < 0: dim(Es) = 2 and dim(Eu) = 1, P1 is unstable.

  • for µ = 0: dim(Es) = 2, dim(Ec) = 1, and nothing can be said about P1 stability observing only the eigenvalues of the Jacobian matrix.

Analogously, the eigenvalues of Jacobian matrix in P2 are:

Therefore:

  • for µ < 0: dim(Es) = 3, P2 is asymptotically stable.

  • for µ > 0: dim(Es) = 2 and dim(Eu) = 1, P2 is unstable.

  • for µ = 0: dim(Es) = 2 and dim(Ec) = 1, nothing can be said about P2 stability observing only the eigenvalues of the Jacobian matrix.

When µ = 0 the stability of P1 and P2 changes, i.e., there is a bifurcation.

In this case, we have to analyze how the system behavior depends on µ, restricted to its central manifold. In order to do this, we rewrite the equations (10) including the parameter in the dynamics [27].

The eigenvalues of the Jacobian matrix associated to system (11) calculated in (x1, x2, x3, µ) = (0, 0, 0, 0) are l1 = -1, l2 = -1 and l3 = 0.

By using the Taylor approximation sin x2 = x2 - + and Jordan canonical form, we have:

with

So, from Central Manifold Theorem, the stability of (0, 0, 0) near µ = 0 can be determined by analyzing the vector field restricted to a central manifold Wc(0), 0 Î 3. In our case, we can write:

Wc(0) = {(v1, v2, v3, µ) Î 4 / v1 = h1(x, µ), v2 = h2(x, µ), v3 = x,

hi(0,0) = 0, Dhi(0,0) = 0, i = 1,2}

for x and µ sufficiently small.

By applying the center manifold theorem [7, 27], and considering polynomial approximations, we have:

Replacing (13) in the third equation of (12) and by considering the forth order terms, we have the vector field, reduced to the center manifold Wc(0), giving by:

When we plot the equilibrium states of (14) we can observe from the bifurcation diagram (figure 12) that x = 0 is a stable equilibrium state for µ > 0 and unstable for µ < 0. When µ > 1, two new unstable equilibrium states given by x2 = 1/(µ – 1) are created.


An analogous reasoning can be conducted for P2, and its stability near µ = 0 can be studied by analyzing the bifurcation diagram shown in (figure 13).


4.2 Single-chain OWMS network

The dynamics of a network consisting of a master and two slaves in a single-chain topology is just the same as the considered in equation 8, with F1(t) being the phase of the master. F2(t) and F3(t) are the phases of the slaves.

Changing the time scale by T = µ1t, and considering a step with finite amplitude as the output signal of master, the dynamics of a single-chain OWMS network is given by:

With these equations, we conclude that the dynamics of a single-chain network and of a single-star network, submitted to step inputs, are the same because these networks are described by the same equations.

4.3 Engineering conclusion

By using the Center Manifold Theory we might formulate the problem of the dynamics of single-chain and single-star OWMS networks and useful engineering conclusions could be:

  • Delays are irrelevant when we consider a step with finite amplitude as the input in our problem.

  • For a given phase and frequency initial state, parameter m that determines the existence and stability [13] of the synchronous state changes the dynamics as illustrated in bifurcation diagrams shown in figures 12 and 13.

5 Stability of equilibrium states and frequency errors in a double-star TWMS network

When a TWMS strategy of clock distribution is chosen we have a more robust and accurate performance for the network. But, in this situation, due to the feedback loops between the nodes, frequency errors, even low, propagate along the whole network spoiling the performance.

In this section, we study the problem of frequency error propagation in double-star TWMS networks by using techniques from dynamical systems theory [22] obtaining conditions for existence and stability of synchronous states [13].

The slaves considered are second order PLLs with a time constant µ. The architecture is the double-star with a master M and N – 1 slaves.

The master is an oscillator with phase FM(t). Signal propagation time from the master to the ith slave is indicated by t1i, and, from the ith slave to the master, by ti1, for i = 2,..., N.

Phases of oscillators output in this network are defined as follows:

  • Master oscillator

  • ith slave-node oscillator

  • Master output phase

Modeling each ith slave, i = 2,3,4,...,N, with a PLL equation, as we have seen in section 2, the dynamics can be described as follows:

where µ is the filter cut-off frequency in all nodes, and µi is the ith slave-node PLL gain.

Defining frequency and phase spatial errors by:

Considering phase perturbations of second order with master acceleration WM and slave acceleration Wi, the substitution of equations (16), (17) and (18) in (19), taking into account equation (20), results:

The dynamics is non-linear and depends explicitly on time, so there is no equilibrium state.

Consequently, the oscillator degradation combined with the delays does not allow the system to be locked in the steady state [9, 21].

If we take the derivatives in equation (21) and consider a linear approximation by expanding the non-linear terms in Taylor series [10, 9, 21], we have:

Considering the state variables:

the system becomes:

This system admits an equilibrium state which corresponds to constant frequency spatial errors

M, i and non-limited phase spatial errors jM, i.

Acceleration spatial error

M, i tends to a zero stationary state. After the transient states, the acceleration of any slave follows the acceleration of the master.

The linear part of the new system, around the equilibrium state, can be represented by:

If we assume the simplifying but realist fact that delays are close to the gain of the PLLs, with

t1i = ti1 = t, µi = v and Wi = W, i = 2, 3, ..., N,

the equilibrium state is given by:

Developing the expressions and verifying the influence of the number of nodes, we can write the equilibrium state by using mathematical induction as:

The linear part of this new system around the equilibrium state is represented by matrix A:

Calculating the eigenvalues of A by using MAPLE V [1]:

Examining the eigenvalues we observe that the equilibrium point of the system, having constant acceleration error, is asymptotically stable for any physical possible value of the parameters.

5.1 Another engineering conclusion

As a result of long-term instability of the master oscillator in a double star TWMS network, slave nodes do not synchronize in phase with master. Phase spatial error jM,i between a slave and the master is unlimited as a consequence of equation (21).

However, in some practical situations, we can relate the propagation times with the gains of PLLs making the frequency errors controllable.

Summarizing:

  • It follows from equation (21) that the double-star TWMS network does not present synchronous solution when oscillators suffer a phase acceleration and frequency errors propagate along the whole network.

  • It follows from (25) that the frequency spatial errors do not dependent on the number of slaves but only on the acceleration of the master and the slaves.

  • Restrictions on the stability domain of equilibrium state depend on the relation between PLLs gains and signal delays.

6 Final comment

The problem of choosing a good circuit for time distribution networks has many aspects mainly related to topology and parameter design.

Modeling the several possible networks gives a nonlinear high dimension ordinary differential equation and analytical solutions are hard to be obtained.

Using Dynamical System Theory is a very useful tool for this kind of problem providing existence and stability conditions for the synchronous state of the networks, relating circuit parameters, transmission delays and deviations.

Received:05/ I/04. Accepted:15/IV/04

#592/04.

  • [1] M.L. Abell and J.P. Braselton, "Maple V by example", Academic Press, London, (1999).
  • [2] J.C. Bellamy, Digital network synchronization, IEEE Communications Magazine (1995), 70-83.
  • [3] H. de Bellescize, La reception synchrone, Onde Electrique, 11 (1932), 230-240.
  • [4] R.E. Best, "Phase-locked loops - 4th edition", McGraw Hill, New York, (1999).
  • [5] H.J. Blinchkoff, All-Pole Phase-Locked Tracking Filters, IEEE Trans. Commun., vol COM-30, pp. 2312-2318, Oct. 1982.
  • [6] C.S. Bregni, A historical perspective on telecomunnications network synchronization, IEEE Commun. Magazine, pp. 158-166, June 1998.
  • [7] J. Carr, "Applications of Centre Manifold Theory", Springer, New York, (1981).
  • [8] R.D. Cideciyan and W.C. Lindsey, Effects of Long-Term Clock Instability on Master-Slave Networks, IEEE Trans. Commun., vol COM-35, pp. 950-955, Sep. 1987.
  • [9] P.A. Garcia, Redes simples de malhas de sincronismo de fase: uma análise via teoria de sistemas dinâmicos, Dissertade Mestrado, Universidade Mackenzie, Săo Paulo, 2000.
  • [10] F.M. Gardner, "Phase-lock Techniques", John Wiley & Sons, New York, (1979).
  • [11] A. Gersho and B.J. Karafin, Mutual Synchronization of Geographically Separeted Oscillators, The Bell System Technical Journal (1966), 1689-1704.
  • [12] J. Guckenheimer and P. Holmes, "Nonlinear Oscillations, Dynamical Systems and Bifurcation of Vector Fields 5th edition", Springer, New York, (1997).
  • [13] T. Kailath, "Linear Systems", Prentice-Hall, New Jersey, (1980).
  • [14] M. Karnaugh, A Model for Organic Synchronization of Communications Systems, The Bell System Technical Journal (1966), 1705-1735
  • [15] W.C. Lee, J. Lee, D. Choi, M. Kim and C. Lee, The Distributted Controller Architecture for a Masterarm and its Application to Teleoperation with Force Feedback, IEEE International Conference on Robotics and Automation, pp. 213-218, (1999).
  • [16] W.C. Lindsey, F. Ghazvinian, W.C. Hagmann and K. Dessouky, Network Synchronization, Proceedings of the IEEE, vol 73, 1445-1467, Oct. 1985.
  • [17] W.C. Lindsey and M.K. Simon, "Telecommunication Systems Engineering", Dover, New York, (1973).
  • [18] C.N. Marmo, V.F. de Faria, L.H.A. Monteiro e J.R.C. Piqueira, Sincronismo em Redes Mestre-Escravo: comparaçăo de topologias, XIV Congresso Brasileiro de Automática, Natal-RN-Brasil, 2002, pp 691-696.
  • [19] M. Schwartz, "Transmissăo de Informaçăo, Modulaçăo e Ruído", Guanabara Dois, Rio de Janeiro, 1979 (translated from the edition: McGraw Hill, New York, 1970).
  • [20] K. Ogata, "Modern Control Engineering", Prentice Hall, New Jersey, (1997).
  • [21] J.R.C. Piqueira, Uma contribuiçăo ao estudo das redes com malhas de sincronismo de fase, Tese de Livre Docęncia, EPUSP, Săo Paulo, (1997).
  • [22] S.A. Castillo-Vargas, Propagaçăo de erros de freqüęncia em redes mestre-escravo, Dissertaçăo de Mestrado, EPUSP, Săo Paulo, (2002).
  • [23] M. Sexton and A. Reid, "Transmission Networking: SONET and Synchronous Digital Hierarchy", Norwood, MA, Artech House, (1992).
  • [24] G. Shao, F. Berman and R. Wolski, Master/slave Computing on the Grid, Proceedings of the IEEE COMPSAC, 2000.
  • [25] S. Sohail and G. Raj, Replication of Multimedia Data using Master-Slave architecture, Proceedings of the IEEE 21st COMPSAC, (1997).
  • [26] A. Weinberg and B. Liu, Discrete Time Analysis of Non-uniform Sampling: First and Second Order Digital Phase-Locked Loops, IEEE Trans. Commun., vol COM-22, pp. 123-137, Feb. 1974.
  • [27] S. Wiggins, "Introduction to Applied Nonlinear Dynamical Systems and Chaos", Springer, New York, (1990).
  • *
    Sponsored by CNPq of Brazil.
  • Publication Dates

    • Publication in this collection
      07 May 2009
    • Date of issue
      Apr 2005

    History

    • Accepted
      15 Apr 2004
    • Received
      05 Jan 2004
    Sociedade Brasileira de Matemática Aplicada e Computacional Sociedade Brasileira de Matemática Aplicada e Computacional - SBMAC, Rua Maestro João Seppe, nº. 900 , 16º. andar - Sala 163, 13561-120 São Carlos - SP Brasil, Tel./Fax: 55 16 3412-9752 - São Carlos - SP - Brazil
    E-mail: sbmac@sbmac.org.br