On Mechanical and Chaotic Problem Modeling and Numerical Simulation Using Electric Networks

: After reviewing the use of electrical circuit elements to model dynamic processes or the operation of devices or equipment, both in real laboratory implementations and through ideal circuits implemented in simulation software, a network model design protocol is proposed. This approach, following the basic rules of circuit theory, makes use of controlled generators to implement any type of nonlinearity contained in the governing equations. Such a protocol constitutes an interesting educational tool that makes it possible for nonexpert students in mathematics to design and numerically simulate complex physical processes. Three applications to mechanical and chaotic problems are presented to illustrate the versatility of the proposed protocol.


Introduction
The simulation of physical or engineering processes by means of models designed in the laboratory is a very old resource that allows predicting both the steady state solution and the time evolution of these systems in a simpler, more economical, and often reliable way.For the elaboration of these models, based on the mathematical equivalence of the physical equations of the real process and the equations of the laboratory prototype, different techniques were used, such as the use of membranes (soap film analogy), chemical solutions (electrolytic tank), graph paper (analog field plotter), construction of electrical circuits, and reduced models based on scale factors.Thus, and sticking to the subject of electrical analogy, based on the laws of electric circuits, we should mention the first attempts to successfully simulate processes governed by Laplace's equation through the construction of an electrolytic tank [1], or a century later with the use of graphite paper (analog field plotter [2]).More recently, Arvinti et al. [3] implemented a laboratory electrical model to solve the Laplace equation in the whole domain, approaching the solutions using Lagrange polynomials.During the decades between 1940 and 1970, large analog equipment consisting of resistors and capacitors were developed that allowed the simulation of heat and mass flow processes both linear, with the 'heat and mass analyzer' of Paschkis and Heisler [4], and nonlinear, by means of the 'differential analyzer' of Karplus and Soroka [5].After the development of computers, these physical models were replaced by numerical computational techniques that directly address the solution to governing equations using a variety of precise methods, such as finite elements, finite differences, and variational techniques.Table 1 lists the different models developed in historical order with their applications.sentences, which allows for implementing any type of physical problem in the models, particularly the so-called controlled generators, which allow for implementing any nonlinear or coupled term that is part of the governing equations; (iii) the computer algorithms developed in these programs being perhaps the most up-to-date, optimized, and computationally powerful, which results in the reliability of the numerical solutions and the reduction of computation time [24]; and (iv) the programming rules for preparing the text files of the models being relatively few and established on the basic theory of electric circuits, i.e., on the constitutive laws of their elements and on the theorems of uniqueness of the electric potential and conservation of electric charge (Kirchhoff's theorems) [22].In fact, the researcher only has to worry about the correct design of the network model-or equivalent circuit-which has to collect the boundary and initial conditions of the problem, forgetting about the algorithms for numerical computation.
Numerous researchers have used this analogy by applying it to complex problems in different fields of physics and engineering, as shown in Table 1.Horno et al. [22] apply it to charge transfer processes in membranes, López-García et al. [7] to the study of the electric double layer in colloids, Cánovas et al. [14] to the Bénard convection cell problem, Bég et al. [11] to magneto-hydrodynamics processes, Serna et al. [12] to lid cavity problems, and García-Ros et al. [15] to the nonlinear consolidation of soils.
The contributions of this work include (i) establishing a protocol for the design of electrical circuit models (network models) ruled by ordinary differential equations, linear or not; (ii) adapting the above design to the programming language of the Pspice software, exploring the use of controlled generators as essential elements for the implementation of second-order derivative terms and other complex terms of the equation; and (iii) to make use of the above protocol as an educational tool to make it possible for students who are not experts in mathematical and/or numerical calculation to simulate the solution of any kind of ordinary differential equations.
Section 2 illustrates, step by step, the application of the protocol of the network simulation method.The analytical solution to such problems, when it exists, is neither immediate nor simple.In this section, constitutive laws of the basic passive elements (resistor and capacitor) and the types of controlled sources or generators are described, explaining also the use of theorems for the model design.Section 3 includes three applications to dynamic problems, presenting the network models of each problem, the numerical simulation, the graphical solutions with Pspice, and other aspects of the physical behavior of the system, added or derived from these solutions.Finally, the conclusions are presented in the last section.

Design of the Network Models-The Electrical Components of the Model
The starting point is the mathematical model of the problem, i.e., the governing equation and the initial conditions.Each term of the governing equation is considered an electric current that is implemented in the model by means of a single component that is arranged in parallel-between the main node and the common or mass-with the components associated with the other terms.All components are connected between the main node and the common or ground node in such a way the balance of the equation is ensured by the conservation of the currents converging at the main node (Kirchhoff's first theorem, which is valid in quasi-stationary regime) [25].There are as many components as there are terms in the differential equation.The equilibrium imposed by the law of conservation of the free electric charges (inherent in the circuit theory software itself) forces the potential at the main node (value of the variable sought) to be univocally the one that satisfies this equilibrium, i.e., the solution to the problem.
Whatever the terms of the governing equation and its complexity, its simulation can be carried out by means of a network model (electrical circuit) consisting of only a few electrical components: (i) capacitors and resistors to implement the linear terms; and (ii) constant or controlled sources to implement the rest of the terms.Table 2 lists the symbols of all components and their constitutive equations.For resistors and capacitors, the constitutive equations or relationships between the electric current (i R ) and the potential difference at their ends (v R ) are i R = v R R and i C = C{ dv C dt }, respectively.For the correct implementation of the component in the model, the directions of the electric current and the potential difference must be consistent with those indicated in the table.In the constitutive equations of the controlled sources, the output variable (voltage or current, v out or i out ) is an arbitrary function-which is defined by programming-of the input variable (voltage or current, v in or i in ).Constant sources implement constant terms of the equation and time-dependent sources implement time-dependent terms.

Resistor
Whatever the terms of the governing equation and its complexity, its simulation can be carried out by means of a network model (electrical circuit) consisting of only a few electrical components: (i) capacitors and resistors to implement the linear terms; and (ii) constant or controlled sources to implement the rest of the terms.Table 2 lists the symbols of all components and their constitutive equations.For resistors and capacitors, the constitutive equations or relationships between the electric current (i ) and the potential difference at their ends (v ) are i = and i = C{ }, respectively.For the correct implementation of the component in the model, the directions of the electric current and the potential difference must be consistent with those indicated in the table.In the constitutive equations of the controlled sources, the output variable (voltage or current, v or i ) is an arbitrary function-which is defined by programming-of the input variable (voltage or current, v or i ).Constant sources implement constant terms of the equation and time-dependent sources implement time-dependent terms.

Component
Symbol Constitutive Equation There are four types of controlled sources, although the most commonly used are the voltage-controlled current sources, to implement terms of the equation (linear or not) that are a function of the dependent variable, and current-controlled current sources to implement the second-and higher-order derivative terms.For some applications, for example, when the governing equation contains time-dependent terms, it is necessary to implement in the model an auxiliary circuit to define the time variable.

Basic Circuits
Among the most common summands or terms of a differential equation, we first distinguish the derivative terms.The order of the largest existing derivative classifies the type of equation: first order, second order, etc.Since the solution sought is going to be the voltage v (t) at the main node of the model (node I), the implementation of a capacitor (C ) in one of the branches of the circuit of that node sets a current of value i = C (dv /dt), which is equivalent to implement the equivalent term of the first derivative (when it exists).The value of C allows us to adjust the coefficient of this term to its particular value Whatever the terms of the governing equation and its complexity, its simulation can be carried out by means of a network model (electrical circuit) consisting of only a few electrical components: (i) capacitors and resistors to implement the linear terms; and (ii) constant or controlled sources to implement the rest of the terms.Table 2 lists the symbols of all components and their constitutive equations.For resistors and capacitors, the constitutive equations or relationships between the electric current (i ) and the potential difference at their ends (v ) are i = and i = C{ }, respectively.For the correct implementation of the component in the model, the directions of the electric current and the potential difference must be consistent with those indicated in the table.In the constitutive equations of the controlled sources, the output variable (voltage or current, v or i ) is an arbitrary function-which is defined by programming-of the input variable (voltage or current, v or i ).Constant sources implement constant terms of the equation and time-dependent sources implement time-dependent terms.

Component
Symbol Constitutive Equation There are four types of controlled sources, although the most commonly used are the voltage-controlled current sources, to implement terms of the equation (linear or not) that are a function of the dependent variable, and current-controlled current sources to implement the second-and higher-order derivative terms.For some applications, for example, when the governing equation contains time-dependent terms, it is necessary to implement in the model an auxiliary circuit to define the time variable.

Basic Circuits
Among the most common summands or terms of a differential equation, we first distinguish the derivative terms.The order of the largest existing derivative classifies the type of equation: first order, second order, etc.Since the solution sought is going to be the voltage v (t) at the main node of the model (node I), the implementation of a capacitor (C ) in one of the branches of the circuit of that node sets a current of value i = C (dv /dt), which is equivalent to implement the equivalent term of the first derivative (when it exists).The value of C allows us to adjust the coefficient of this term to its particular value Whatever the terms of the governing equation and its complexity, its simulation can be carried out by means of a network model (electrical circuit) consisting of only a few electrical components: (i) capacitors and resistors to implement the linear terms; and (ii) constant or controlled sources to implement the rest of the terms.Table 2 lists the symbols of all components and their constitutive equations.For resistors and capacitors, the constitutive equations or relationships between the electric current (i ) and the potential difference at their ends (v ) are i = and i = C{ }, respectively.For the correct implementation of the component in the model, the directions of the electric current and the potential difference must be consistent with those indicated in the table.In the constitutive equations of the controlled sources, the output variable (voltage or current, v or i ) is an arbitrary function-which is defined by programming-of the input variable (voltage or current, v or i ).Constant sources implement constant terms of the equation and time-dependent sources implement time-dependent terms.

Component
Symbol Constitutive Equation There are four types of controlled sources, although the most commonly used are the voltage-controlled current sources, to implement terms of the equation (linear or not) that are a function of the dependent variable, and current-controlled current sources to implement the second-and higher-order derivative terms.For some applications, for example, when the governing equation contains time-dependent terms, it is necessary to implement in the model an auxiliary circuit to define the time variable.

Basic Circuits
Among the most common summands or terms of a differential equation, we first distinguish the derivative terms.The order of the largest existing derivative classifies the type of equation: first order, second order, etc.Since the solution sought is going to be the voltage v (t) at the main node of the model (node I), the implementation of a capacitor (C ) in one of the branches of the circuit of that node sets a current of value i = C (dv /dt), which is equivalent to implement the equivalent term of the first derivative (when it exists).The value of C allows us to adjust the coefficient of this term to its particular value Whatever the terms of the governing equation and its complexity, its simulation can be carried out by means of a network model (electrical circuit) consisting of only a few electrical components: (i) capacitors and resistors to implement the linear terms; and (ii) constant or controlled sources to implement the rest of the terms.Table 2 lists the symbols of all components and their constitutive equations.For resistors and capacitors, the constitutive equations or relationships between the electric current (i ) and the potential difference at their ends (v ) are i = and i = C{ }, respectively.For the correct implementation of the component in the model, the directions of the electric current and the potential difference must be consistent with those indicated in the table.In the constitutive equations of the controlled sources, the output variable (voltage or current, v or i ) is an arbitrary function-which is defined by programming-of the input variable (voltage or current, v or i ).Constant sources implement constant terms of the equation and time-dependent sources implement time-dependent terms.
Table 2. Electrical components of network models and their constitutive equations.

Component
Symbol Constitutive Equation There are four types of controlled sources, although the most commonly used are the voltage-controlled current sources, to implement terms of the equation (linear or not) that are a function of the dependent variable, and current-controlled current sources to implement the second-and higher-order derivative terms.For some applications, for example, when the governing equation contains time-dependent terms, it is necessary to implement in the model an auxiliary circuit to define the time variable.

Basic Circuits
Among the most common summands or terms of a differential equation, we first distinguish the derivative terms.The order of the largest existing derivative classifies the type of equation: first order, second order, etc.Since the solution sought is going to be the voltage v (t) at the main node of the model (node I), the implementation of a capacitor (C ) in one of the branches of the circuit of that node sets a current of value i = C (dv /dt), which is equivalent to implement the equivalent term of the first derivative (when it exists).The value of C allows us to adjust the coefficient of this term to its particular value Whatever the terms of the governing equation and its complexity, its simulation can be carried out by means of a network model (electrical circuit) consisting of only a few electrical components: (i) capacitors and resistors to implement the linear terms; and (ii) constant or controlled sources to implement the rest of the terms.Table 2 lists the symbols of all components and their constitutive equations.For resistors and capacitors, the constitutive equations or relationships between the electric current (i ) and the potential difference at their ends (v ) are i = and i = C{ }, respectively.For the correct implementation of the component in the model, the directions of the electric current and the potential difference must be consistent with those indicated in the table.In the constitutive equations of the controlled sources, the output variable (voltage or current, v or i ) is an arbitrary function-which is defined by programming-of the input variable (voltage or current, v or i ).Constant sources implement constant terms of the equation and time-dependent sources implement time-dependent terms.

Component
Symbol Constitutive Equation There are four types of controlled sources, although the most commonly used are the voltage-controlled current sources, to implement terms of the equation (linear or not) that are a function of the dependent variable, and current-controlled current sources to implement the second-and higher-order derivative terms.For some applications, for example, when the governing equation contains time-dependent terms, it is necessary to implement in the model an auxiliary circuit to define the time variable.

Basic Circuits
Among the most common summands or terms of a differential equation, we first distinguish the derivative terms.The order of the largest existing derivative classifies the type of equation: first order, second order, etc.Since the solution sought is going to be the voltage v (t) at the main node of the model (node I), the implementation of a capacitor (C ) in one of the branches of the circuit of that node sets a current of value i = C (dv /dt), which is equivalent to implement the equivalent term of the first derivative (when it exists).The value of C allows us to adjust the coefficient of this term to its particular value Whatever the terms of the governing equation and its complexity, its simulation can be carried out by means of a network model (electrical circuit) consisting of only a few electrical components: (i) capacitors and resistors to implement the linear terms; and (ii) constant or controlled sources to implement the rest of the terms.Table 2 lists the symbols of all components and their constitutive equations.For resistors and capacitors, the constitutive equations or relationships between the electric current (i ) and the potential difference at their ends (v ) are i = and i = C{ }, respectively.For the correct implementation of the component in the model, the directions of the electric current and the potential difference must be consistent with those indicated in the table.In the constitutive equations of the controlled sources, the output variable (voltage or current, v or i ) is an arbitrary function-which is defined by programming-of the input variable (voltage or current, v or i ).Constant sources implement constant terms of the equation and time-dependent sources implement time-dependent terms.

Component
Symbol Constitutive Equation There are four types of controlled sources, although the most commonly used are the voltage-controlled current sources, to implement terms of the equation (linear or not) that are a function of the dependent variable, and current-controlled current sources to implement the second-and higher-order derivative terms.For some applications, for example, when the governing equation contains time-dependent terms, it is necessary to implement in the model an auxiliary circuit to define the time variable.

Basic Circuits
Among the most common summands or terms of a differential equation, we first distinguish the derivative terms.The order of the largest existing derivative classifies the type of equation: first order, second order, etc.Since the solution sought is going to be the voltage v (t) at the main node of the model (node I), the implementation of a capacitor (C ) in one of the branches of the circuit of that node sets a current of value i = C (dv /dt), which is equivalent to implement the equivalent term of the first derivative (when it exists).The value of C allows us to adjust the coefficient of this term to its particular value Whatever the terms of the governing equation and its complexity, its simulation can be carried out by means of a network model (electrical circuit) consisting of only a few electrical components: (i) capacitors and resistors to implement the linear terms; and (ii) constant or controlled sources to implement the rest of the terms.Table 2 lists the symbols of all components and their constitutive equations.For resistors and capacitors, the constitutive equations or relationships between the electric current (i ) and the potential difference at their ends (v ) are i = and i = C{ }, respectively.For the correct implementation of the component in the model, the directions of the electric current and the potential difference must be consistent with those indicated in the table.In the constitutive equations of the controlled sources, the output variable (voltage or current, v or i ) is an arbitrary function-which is defined by programming-of the input variable (voltage or current, v or i ).Constant sources implement constant terms of the equation and time-dependent sources implement time-dependent terms.

Component
Symbol Constitutive Equation There are four types of controlled sources, although the most commonly used are the voltage-controlled current sources, to implement terms of the equation (linear or not) that are a function of the dependent variable, and current-controlled current sources to implement the second-and higher-order derivative terms.For some applications, for example, when the governing equation contains time-dependent terms, it is necessary to implement in the model an auxiliary circuit to define the time variable.

Basic Circuits
Among the most common summands or terms of a differential equation, we first distinguish the derivative terms.The order of the largest existing derivative classifies the type of equation: first order, second order, etc.Since the solution sought is going to be the voltage v (t) at the main node of the model (node I), the implementation of a capacitor (C ) in one of the branches of the circuit of that node sets a current of value i = C (dv /dt), which is equivalent to implement the equivalent term of the first derivative (when it exists).The value of C allows us to adjust the coefficient of this term to its particular value Whatever the terms of the governing equation and its complexity, its simulation can be carried out by means of a network model (electrical circuit) consisting of only a few electrical components: (i) capacitors and resistors to implement the linear terms; and (ii) constant or controlled sources to implement the rest of the terms.Table 2 lists the symbols of all components and their constitutive equations.For resistors and capacitors, the constitutive equations or relationships between the electric current (i ) and the potential difference at their ends (v ) are i = and i = C{ }, respectively.For the correct implementation of the component in the model, the directions of the electric current and the potential difference must be consistent with those indicated in the table.In the constitutive equations of the controlled sources, the output variable (voltage or current, v or i ) is an arbitrary function-which is defined by programming-of the input variable (voltage or current, v or i ).Constant sources implement constant terms of the equation and time-dependent sources implement time-dependent terms.

Component
Symbol Constitutive Equation There are four types of controlled sources, although the most commonly used are the voltage-controlled current sources, to implement terms of the equation (linear or not) that are a function of the dependent variable, and current-controlled current sources to implement the second-and higher-order derivative terms.For some applications, for example, when the governing equation contains time-dependent terms, it is necessary to implement in the model an auxiliary circuit to define the time variable.

Basic Circuits
Among the most common summands or terms of a differential equation, we first distinguish the derivative terms.The order of the largest existing derivative classifies the type of equation: first order, second order, etc.Since the solution sought is going to be the voltage v (t) at the main node of the model (node I), the implementation of a capacitor (C ) in one of the branches of the circuit of that node sets a current of value i = C (dv /dt), which is equivalent to implement the equivalent term of the first derivative (when it exists).The value of C allows us to adjust the coefficient of this term to its particular value There are four types of controlled sources, although the most commonly used are the voltage-controlled current sources, to implement terms of the equation (linear or not) that are a function of the dependent variable, and current-controlled current sources to implement the second-and higher-order derivative terms.For some applications, for example, when the governing equation contains time-dependent terms, it is necessary to implement in the model an auxiliary circuit to define the time variable.

Basic Circuits
Among the most common summands or terms of a differential equation, we first distinguish the derivative terms.The order of the largest existing derivative classifies the type of equation: first order, second order, etc.Since the solution sought is going to be the voltage v I (t) at the main node of the model (node I), the implementation of a capacitor (C 1 ) in one of the branches of the circuit of that node sets a current of value i C 1 = C 1 (dv I /dt), which is equivalent to implement the equivalent term of the first derivative (when it exists).The value of C 1 allows us to adjust the coefficient of this term to its particular value in the equation.Each of the following derivatives is implemented by two components: a current-controlled voltage generator (H ccvs ) and a new capacitor (C).For example, the second derivative (node II) does with the pair H ccvs,1 and C 2 .The input of H ccvs,1 , defined by programming in the source specification, is the current i C 1 , while the output voltage (connected to the ends of C 2 ) has the same numerical value as that current, i.e., i.e., the second derivative term of the equation.Again, the value of C 2 allows for the adjustment of the coefficient of such a term.The third and successive derivative terms follow the same rule for their implementation in the model, using the pair C 3 − H ccvs,2 to implement the third derivative, and so on.Figure 1a shows these components within the network model.The independent term of constant value is implemented directly by a constant-current generator (I) with a sense consistent with the algebraic sign of the term in the equation.Figure 1b implements a constant-valued term of numerical value equal to that of the gener-ator current, i G .A term of the form a o x (with x the dependent variable, v I ) is implemented, according to its constitutive equation, by a resistor (R) of value R = 1/a o , as shown in Figure 1c.The remaining terms are always implemented with a new H ccvs , and as many as there are summands in the equation.The expression of these terms, sometimes complex, may require auxiliary circuits for implementation.Figure 1d

Text Files
The model can be entered into the program using the schematics tool included in the software itself, which uses the standard circuit theory symbols, or through a text file using a minimal set of writing rules.For resistors, capacitors, and constant sources, it is sufficient to indicate their name, the nodes of the circuit to which they are connected, and their value, and for controlled sources, their name, the nodes or element from which they read the input, the nodes where their output is connected, and the control function that is specified by programming.Table 3 shows the specification of the text lines corresponding to the elements listed in Table 2 (Pspice [23]).In the Section 3, the text files of the studied models are shown and explained.

Applications and Simulation
Three applications to mechanical and chaotic dynamic processes have been selected that sufficiently illustrate the application of the protocol of the network simulation method.These processes contain in their governing equations both linear and nonlinear summands-terms in the first and second derivative, terms as a function of the dependent variable, rational exponents, etc.-and can therefore be taken as very representative of any kind of dynamic process.For each application, the design of the network model is described in detail, showing its text file to be read and executed in the software.The simulation results obtained with the graphical output environment of the software itself are also briefly discussed.

Mass Falling in Air or Viscous Fluid
This simple process, ruled by a differential equation of three addends and nonlinear because of the existence of the exponent b o in general different from unity, governs the motion of a mass falling in air or viscous fluid.It is also named the skydiver equation, as it governs the movement of a parachutist falling towards the earth by the action of gravitational force but slowed down by air friction with a force dependent on the instantaneous velocity.The balance of forces (Newton's second law) is given by the equation dt 2 .With P = mg, and expressing the frictional force in the form f friction = γv a o = γ dx dt a o , the balance equation yields with m the mass of the parachutist, v its instantaneous velocity, g the acceleration of gravity, γ and a o constants of the friction force expression, and x the dependent variable position.The initial conditions are: x t=0 = v t=0 = 0.The network model, shown in Figure 2, consists of a main circuit (node I) with three branches corresponding to the three summands of the equation and two auxiliary circuits (nodes II and III), which implement the first and second derivatives.The solution to the problem is the main node voltage, x(t) = v I (t).In the first auxiliary circuit, formed by E vcvs,1 and C 1 , with C 1 = 1, the generator output current is equal to the input voltage, i E vcvs,1 = v I , so i with m the mass of the parachutist, v its instantaneous velocity, g the acceleration of gravity, γ and a constants of the friction force expression, and x the dependent variable position.The initial conditions are: The network model, shown in Figure 2, consists of a main circuit (node I) with three branches corresponding to the three summands of the equation and two auxiliary circuits (nodes II and III), which implement the first and second derivatives.The solution to the problem is the main node voltage, x(t) = v (t).In the first auxiliary circuit, formed by E , and C , with C = 1, the generator output current is equal to the input voltage, i   The text file of the model is as follows: .Finally, the constant term g is implemented with the constant current generator I 1 towards the main node (because of the negative sign), i out,I 1 = g.Figure 2 shows the network model of the problem.The initial conditions of position and velocity, x ini and v ini , are applied as initial voltage on capacitors C 1 and C 2 , v ini,C 1 = v t=0 = 0, and v ini,C 2 = x t=0 = 0, respectively.The flow chart in Figure 3   The text file of the model is as follows: The asterisk implies that this line of text is not executed by the program.In the voltage source, time is read is v time , while the sentence '.TRANS' specifies the time interval to be simulated.The curves in Figure 4, in the Pspice graphical environment, show the simulation results-position (x), velocity (v) and acceleration (a), as a function of time-for the parameters g = 9.81, m = 1, x ini = 0, γ = 2 and 1.8, a o = 1.5 and 0.6, and v ini = 0 and 5.Note that in both cases, as expected, the location increases monotonically until lineal, and velocity increases progressively until a steady value is reached and acceleration diminished, converging to a zero value.Initial conditions determine the starting point of the curves.The asterisk implies that this line of text is not executed by the program.In the voltage source, time is read is v , while the sentence '.TRANS' specifies the time interval to be simulated.The curves in Figure 4, in the Pspice graphical environment, show the simulation results-position (x), velocity (v) and acceleration (a), as a function of time-for the parameters g = 9.81, m = 1, x = 0, γ = 2 and 1.8, a = 1.5 and 0.6, and v = 0 and 5.Note that in both cases, as expected, the location increases monotonically until lineal, and velocity increases progressively until a steady value is reached and acceleration diminished, converging to a zero value.Initial conditions determine the starting point of the curves.

Crimped Bead Sliding on a Parabolic Shaped Wire
In Figure 5, by the action of gravity, the ball attached to the wire falls (sliding without friction) following the parabolic trajectory y = a x .The balance between the gravitational force and the normal reaction of the wire (Newton's law) allows us to write the governing equation in the following form: The mathematical model is completed with the initial conditions that we will choose simply, as x = 5 and v = 0.The position y(t) is given by the trajectory equation once the solution x(t) is obtained.

Crimped Bead Sliding on a Parabolic Shaped Wire
In Figure 5, by the action of gravity, the ball attached to the wire falls (sliding without friction) following the parabolic trajectory y = a o x 2 .The balance between the gravitational force and the normal reaction of the wire (Newton's law) allows us to write the governing equation in the following form: The mathematical model is completed with the initial conditions that we will choose simply, as x t=0 = 5 and v t=0 = 0.The position y(t) is given by the trajectory equation once the solution x(t) is obtained.hand, as it is an undamped harmonic motion, there is no loss of energy due to friction, so the height reached by the ball at the ends of the motion is the same and its value is  hand, as it is an undamped harmonic motion, there is no loss of energy due to friction, so the height reached by the ball at the ends of the motion is the same and its value is Curves show a clear influence of the nonlinear terms of the equation.On the one hand, as it is an undamped harmonic motion, there is no loss of energy due to friction, so the height reached by the ball at the ends of the motion is the same and its value is a o (x t=0 ) 2 .Because of the symmetry of the parabola, the maximum and minimum horizontal positions are also the same and their absolute value is x t=0 .On the other hand, Table 4 shows the influence of the parameters a o and b o on the maximum velocity (which occurs at x = 0) and on the oscillation period (τ).While the coefficient b o associated with the elastic restoring force is a clear determinant of the maximum velocity, the coefficient a o , which affects the concavity of the parabola, does not influence this velocity.As for the period of oscillation, it decreases with increasing b o and increases with increasing a o .
oscillation, it decreases with increasing bo and increases with increasing ao.

The van der Pol Oscillator
The mathematical model of this nonlinear oscillator [25] is with the initial conditions v = v and x = x , respectively.The network model, shown in Figure 8, which retains the auxiliary circuits of the first application (to implement the current ) and of the second (to implement the voltage ), has four branches

The van der Pol Oscillator
The mathematical model of this nonlinear oscillator [26] is with the initial conditions v t=0 = v o and x t=0 = x o , respectively.The network model, shown in Figure 8, which retains the auxiliary circuits of the first application (to implement the current d 2 x dt 2 ) and of the second (to implement the voltage dx dt ), has four branches in its main circuit (node I).The first one implements the term d 2 x dt 2 through the controlled source F cccs .The second branch, which implement the term x 2 ( dx dt ) through the controlled source H ccvs,2 , is controlled by the voltages v I = x and v III = dx dt , according to the mathematical expression of the term.The summand x is implemented directly by a resistor (R 3 ) of unity value, while the independent term is implemented by the constant current generator Io, whose output is i i o = 1.
source Fcccs.The second branch, which implement the term x ( ) through the controlled source Hccvs,2, is controlled by the voltages v = x and v = , according to the mathematical expression of the term.The summand x is implemented directly by a resistor (R ) of unity value, while the independent term is implemented by the constant current generator Io, whose output is i = 1.The solutions for position, velocity, and acceleration evolution over time and the position versus velocity phase diagram are shown in Figure 9.For the initial conditions imposed, the motion of this oscillator converges over time to a point where both velocity and acceleration cancel out after a run where these variables oscillate chaotically.The solutions for position, velocity, and acceleration evolution over time and the position versus velocity phase diagram are shown in Figure 9.For the initial conditions imposed, the motion of this oscillator converges over time to a point where both velocity and acceleration cancel out after a run where these variables oscillate chaotically.

Discussions and Conclusions
The numerical solution of differential equations is not currently a problem for students or researchers thanks to the existence of specific software capable of reliably solving any type of equation, linear or nonlinear.However, the application of this software does not allow the user to deepen complementary aspects of knowledge such as the physical

Discussions and Conclusions
The numerical solution of differential equations is not currently a problem for students or researchers thanks to the existence of specific software capable of reliably solving any type of equation, linear or nonlinear.However, the application of this software does not allow the user to deepen complementary aspects of knowledge such as the physical meaning of each term of the equation and its direct effect on the solution.The network models based on the electrical analogy proposed in this work, and the necessary association between physical and electrical variables that such an analogy entails for the elaboration of the model, allows the user to delve into and better understand the phenomenological and physical aspects involved in the governing equation.For this reason, the protocol proposed to design the network model can be considered a powerful educational tool, as it allows students unfamiliar with numerical or differential calculus techniques to easily access the solution and understand the physical processes once the numerical simulation has been carried out.
The design of network models (electrical circuits) has proven to be a useful and accurate tool for simulating a wide range of dynamic processes.The design protocol presented in this work, based on the elementary rules of circuit theory and the constitutive equations of its basic components, includes the incorporation of so-called controlled current or voltage generators.These are capable of implementing in the model any type of summands of the governing equation, such as nonlinear second-order derivative terms, terms that depend on time or on the main variable, etc.Each term of the equation is assumed as an electric current that balances at a common node with the currents of the other terms.The fulfillment of such balancing imposes a unique instantaneous value on the voltage at that node, a value that is the solution to the problem.
The numerical solution is carried out by standard circuit simulation software, such as Pspice.There are two main advantages in the use of electric models: (i) few rules are needed for the design, as the elements that compose the model are very few (resistors, capacitors, and constant or controlled sources); and (ii) the use of powerful computer algorithms contained in these programs, meaning the numerical computation provides a quasi-exact solution to the network model.In addition, the graphical output environment of these programs offers the user an immediate representation of the solutions from which the temporal evolution of the variables, phase diagrams, spectral representation of the harmonic responses, etc., can be selected.
The applications of the proposed protocol to two problems of nonlinear dynamics and one of chaotic motion illustrate the advantages mentioned in the previous paragraph.In all three cases, it is immediate to infer first the dependence of the unknowns of interest, position, velocity, and acceleration on the parameters of the equation, and second, the convergence or not of these unknowns to a stationary final position.

Modelling 2024, 5 , 6 Figure 1 .
Figure 1.Implementation of (a) derivative terms, (b) positive constant term, (c) first-degree term of the dependent variable, and (d) nonlinear first-degree term in the dependent variable.

Figure 1 .
Figure 1.Implementation of (a) derivative terms, (b) positive constant term, (c) first-degree term of the dependent variable, and (d) nonlinear first-degree term in the dependent variable.
implements a term of the form b o x c o .The output current of the generator H ccvs (defined by programming) is i H ccvs = b o x c o = b o v c o I .Other examples are shown in the Section 3.
C 1 = C 1 dv I dt = dx dt .The second auxiliary circuit consists of H ccvs,1 and C 2 , with C 2 = 1.The output current of this controlled source (i H ccvs,1,out = i C 2 ) has the same value as the current i C 1 .Thus, i H ccvs,1,out = C 2 dv NodeIII dt = d 2 x dt 2 .

,.
= v , so i = C = .The second auxiliary circuit consists of H , and C , with C = 1.The output current of this controlled source (i , , = i ) has the same value as the current i .Thus, the first derivatives as currents in C and C , i = and i = , it is immediate to implement the branches of the main circuit corresponding to the summands and γ The first contains generator F , , whose input and output have the same value, i , , = i , , = i = .The second contains generator F , , whose input is i , , = i = , and with output a function of this value defined by programming, i , , = i in,F cccs,3 = γ .Finally, the constant term g is implemented with the constant current generator I towards the main node (because of the negative sign), i , = g.Figure 2 shows the network model of the problem.The initial conditions of position and velocity, x and v , are applied as initial voltage on capacitors C and C , v , = v = 0, and v , = x = 0, respectively.The flow chart in Figure 3 explains the procedure for creating the network model text file.

Figure 2 .
Figure 2. Network model of the mass falling in air or viscous fluid.

Figure 3 .
Figure 3. Flow chart for creating the network model text file.

Figure 2 ..
Figure 2. Network model of the mass falling in air or viscous fluid.Having the first derivatives as currents in C 1 and C 2 , i C 1 = dx dt and i C 2 = d 2 x dt 2 , it is immediate to implement the branches of the main circuit corresponding to the summands d 2 x dt 2 and γ m dx dt a o .The first contains generator F cccs,2 , whose input and output have the same value, i out,F cccs,2 = i in,F cccs,2 = i C 2 = d 2 x dt 2 .The second contains generator F cccs,3 , whose input is i in,F cccs,3 = i C 1 = dx dt , and with output a function of this value defined by programming, i out,F cccs,3 = γ m i in,F cccs,3 a o = γ m dx dt a o

Figure 2 .
Figure 2. Network model of the mass falling in air or viscous fluid.

Figure 3 .
Figure 3. Flow chart for creating the network model text file.

Figure 3 .
Figure 3. Flow chart for creating the network model text file.

Figure 5 .
Figure 5. Sliding bead strung on a fixed wire.

Figure 6 .Figure 7
Figure 6.Network model of the bead crimped to the wire.The text file of the model is as follows:

Figure 5 . 10 Figure 5 .
Figure 5. Sliding bead strung on a fixed wire.The network model is shown in Figure 6.The derivative factors dx dt and d 2 x dt 2 are implemented with circuits similar to those explained in the previous application, auxiliary circuits of nodes II and III.The terms d 2 x dt 2 and a o x 2 d 2 x dt 2 are implemented in the main circuit by the current generators F cccs,2 and F cccs,3 , the term a o x dx dt

Figure 6 .Figure 7
Figure 6.Network model of the bead crimped to the wire.The text file of the model is as follows:

Figure 6 . 2 + 2 }Figure 7
Figure 6.Network model of the bead crimped to the wire.The text file of the model is as follows: *Solution of ordinary differential equations *Governing equation: d 2 x dt 2 1 + a o x 2 + a o x dx dt

Figure 7 .
Figure 7. Dynamic solutions of the motion of the bead strung on a wire.

Figure 7 .
Figure 7. Dynamic solutions of the motion of the bead strung on a wire.

Figure 8 .
Figure 8. Network model of the van der Pol oscillator.

Figure 8 .
Figure 8. Network model of the van der Pol oscillator.

Table 2 .
Electrical components of network models and their constitutive equations.

Table 2 .
Electrical components of network models and their constitutive equations.

Table 2 .
Electrical components of network models and their constitutive equations.

Table 2 .
Electrical components of network models and their constitutive equations.

Table 2 .
Electrical components of network models and their constitutive equations.

Table 2 .
Electrical components of network models and their constitutive equations.

Table 2 .
Electrical components of network models and their constitutive equations.

Table 2 .
Electrical components of network models and their constitutive equations.

Table 3 .
Lines of the text file specifying the electrical components of the network model.

Table 4 .
Influence of parameters ao and bo on solutions x

Table 4 .
Influence of parameters a o and b o on solutions x max , y max , v max , and τ.