- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Physics Research International
Volume 2012 (2012), Article ID 375892, 12 pages
Computer Modeling and Simulation of Ultrasonic Signal Processing and Measurements
Department of Electronics, Adarsha Science J.B. Arts and BIRLA Commerce, Mahavidyalaya, Dhamangaon Rly 444709, India
Received 29 August 2011; Revised 18 October 2011; Accepted 14 December 2011
Academic Editor: Anand Pathak
Copyright © 2012 Y. B. Gandole. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
The system for simulation, measurement, and processing in Graphical User Interface implementation is presented. The received signal from the simulation is compared to that of an actual measurement in the time domain. The comparison of simulated, experimental data clearly shows that acoustic wave propagation can be modeled. The feasibility has been demonstrated in an ultrasound transducer setup for material property investigations. The results of simulation are compared to experimental measurements. The simulation result has good agreement with the experimental data which confirms the validity of the model. The simulation tool therefore provides a way to predict the received signal before anything is built. Furthermore, the use of an ultrasonic simulation package allows for the development of the associated electronics to amplify and process the received ultrasonic signals. Such a virtual design and testing procedure not only can save us time and money, but also can provide better understanding on design failures and allow us to modify designs more efficiently and economically.
Numerical simulations—that is, the use of computers to solve problems by simulating theoretical models—is part of new technology that has taken place alongside pure theory and experiment during the last few decades. Numerical simulations permit one to solve problems that may be inaccessible to direct experimental study or too complex for theoretical analysis. Computer simulations can bridge the gap between analysis and experiment. Numerical simulations analysis and experiment cover mutual weakness of both experiment and theory. These simulations will remain a third dimension in ultrasonic measurements, of equal status and importance to experiment and analysis. It has taken a permanent place in all aspects of ultrasonic measurements from basic research to engineering design. The computer experiment is a new and potentially powerful tool. By combining conventional theory, experiment, and computer simulation, one can discover new and unsolved aspects of natural process. These aspects could often neither have been understood nor revealed by analysis or experiments alone.
There might be many use of ultrasound but a common one is its application to nondestructive evaluation. Pulsed ultrasonic is finding an increasing number of applications in research and industrial nondestructive testing. In such evaluation, one tries to obtain information about the inner parts of an ensemble without dismantling it. In an ultrasonic system, a transducer consists of a collection of material layers. The design and optimization of a multilayered transducer is a complicated engineering task that involves knowledge of physical acoustics, analog electronics, and the acoustical properties of the materials involved. This task is made even more difficult by the lack of available information about frequency and thermal dependencies of these materials characteristics. The optimal combination of suitable materials can be found by trial and error, but not without considerable time and cost, both of which can be minimized through the use of simulations. The aim of this paper is to present a tool which provides a simulation of the received signal prior to construction. Of the different ways to model the electroacoustic system, a total electrical simulation tool is used for the following reasons. First, the modeling of acoustic wave propagation in one dimension by electrical lines can be handled with a certain ease; second the associated electronics used to excite, receive, amplify, and process the signals can be designed to meet the application’s specifications prior to building system. This paper presents a simulation solution to ease the selection process. The electronic software simulation package used is PSPICE . The use of PSPICE provides an opportunity to simulate the complex set of excitation electronics, the ultrasonic transducer, the material under investigation, and the receiving electronics. Electrical analogies of one-dimensional acoustic phenomenon have studied over the years. Mason  modeled electromechanical transducers with a lumped equivalent circuit. Redwood  incorporated a transmission line into Mason’s model to obtain useful information about the transient response of a piezoelectric transducer. With the transmission line, one can represent the time delay necessary for a mechanical signal to travel from one side of the transducer to the other. In the case of a plate transducer, the derivation of both models includes a negative capacitor. Using SPICE and an equivalent circuit approximating the negative capacitor, Morris and Hutchens  simulated Redwood’s implementation of Mason’s model. Krimholtz et al.  presented another equivalent circuit for elementary piezoelectric transducer. Leach  used controlled current and voltage sources instead of transformers. Leach mathematically derives his model by adding terms equal to zero in one of the devices electromechanical equations to obtain the form of the telegraphist’s equation. Püttmer et al.  used a lossy transmission line in Leach’s model to account for acoustical attenuation. Benny et al.  outline a method that has been implemented to predict and measure the acoustic radiation generated by ultrasonic transducers operating into air in continuous wave mode. A comparison of experimental and simulated results for piezoelectric composite, piezoelectric polymer, and electrostatic transducers is then presented to demonstrate some quite different airborne ultrasonic beam-profile characteristics. San Emeterio et al.  present an approximate frequency domain electro-acoustic model for pulsed piezoelectric ultrasonic transmitters which, by integrating partial models of the different stages, allows the computation of the emission transfer function and output force temporal waveform. Hirsekorn et al.  perform numerical simulations of acoustic wave propagation through sonic crystals consisting of local resonators using the local interaction simulation approach (LISA). The current work applies the approach of Püttmer et al.  to liquids and piezoelectric transducer to obtain an electrical analogue of one-dimensional, acoustic wave propagation through such materials. In order to keep things at a manageable level, the following simplifications and assumptions are made. The acoustic propagation travels along one direction and consists of planner longitudinal waves, which are normal to the direction of propagation. The amplitudes are small enough to keep things in linear regions of the devices such that the principle of superposition is not violated. From available material data, such as the modulus of elasticity and Poisson’s ratio, the necessary electrical parameters are deduced. Validation of the theory is achieved by comparing experimental data obtained from different liquids at fixed frequency and temperature.
Here is the thickness, is the force, is the particle velocity, is the electrical voltage, is the electrical current, is the piezoelectric constant, and is the Laplace operator.
The equivalent circuit consists of the static capacitance (capacitance between the electrodes), a transmission line (representing the mechanical part of the piezoelectric transducer), and two controlled sources for coupling between the electrical and mechanical parts of the circuit. Suppose an ultrasonic pulse travels through a medium with a finite speed (m/s). This pulse can be pictured as a disturbance to which the medium reacts. In the case of longitudinal wave, the disturbance is a compression or rarefaction of matter, which the medium displaces to return to its equilibrium state. The compression and rarefaction within the medium is related to its density (kg/m3), and the restoring force is related to the medium’s bulk modulus (pa) . Their relationship to the speed of sound is
Similarly, in an electrical transmission line an electrical pulse can travel through it. These pulses are received at the other end of the line after a very short but finite time. The pulses travel at a certain velocity. Similar to the acoustic wave, the electrical pulses are concentrations and rarefaction of electrons within the transmission line . A distributed-parameter network with the circuit parameters distributed throughout the line can approximate a lossy transmission line. One line segment with the length can be approximated with an electric circuit as in the following (Figure 2).
This lossy transmission line model is described by four lumped parameters.(i) is the resistance in both conductors per unit length in Ω/m,(ii) is the inductance in both conductors per unit length in H/m,(iii) is the conductance of the dielectric media per unit length in S/m,(iv) is the capacitance between the conductors per unit length in F/m,
where and are zero under lossless conditions.
To derive the parameters we begin by using Kirchhoff’s voltage law on the circuit in Figure 1: which can be written as and then letting , we get
Then we have one equation containing and . To get another equation relating and , we apply Kirchhoff’s current law on the circuit and get and letting in this equation also, we get
The first order partial differential equations, (4) and (6), are called the general transmission line equations. These equations can be simplified if the voltage and the current are time-harmonic cosine functions: where is the angular frequency. Using (7), the general transmission line equations in (4) and (5) can be written as These equations are called the time-harmonic transmission line equations. Equation (8) can be used to derive the propagation constant and the characteristic impedance of the line. By differentiating (8) with respect to , we get where is the propagation constant: The real part, , of the propagation constant is called the attenuation constant in Np/m and the imaginary part, , is called the phase constant of the line in rad/m. The general solution to the differential equation (10) is The time dependence into (11) can be obtained by multiplying to (11): Equation (12) describes two traveling waves; one traveling in the positive direction with an amplitude which decays at a rate , while the other with an amplitude travels in the opposite direction with the same rate of decay. The same type of differential equations governs the propagation of an acoustical wave. In the case of harmonic waves, corresponding to (10) and (9), we have the lossy linearized acoustic plane wave equations: where (Pa) is the pressure and (m/s) is the particle velocity . Equivalent to is the complex wave number composed of an attenuation constant (Np/m) and a wave number (rad/m). The general solution of the wave equation (15) is and is identical to the solutions obtained in transmission line’s represented by (12). Equation (13) has a solution of the same form. The complex wave number is given as where is the relaxation time .
In order to unify the two theories, an impedance type analogy is chosen where mechanical force is represented by voltage and current represents particle velocity. The characteristic impedance becomes important at the boundaries because there the continuity conditions have to be satisfied. Pressure and normal particle velocity must be continuous at the boundaries as voltage and current must be continuous at connections. For the lossy transmission line , the characteristic impedance is and for the lossy acoustic medium , the characteristic acoustic impedance is where is the density of the medium.
Expanding (17) and (10), we obtain Considering small but nonnegligible losses, where , , and the second term of (23) is negligible, leaving the characteristic impedance as , Similarly from (18), the low loss acoustical characteristic impedance can be approximated as . Also, the wave number from (16) become . To correlate the two characteristic impedances, we choose an impedance type analogy (Figure 1), in which the force (and not pressure) is represented by voltage and the particle velocity is represented by current. The equivalence between the two systems is where (m2) is the cross-sectional area of the acoustic beam.
Assisted with the definition of the low loss characteristic impedances equation, the following relationships can be obtained The real part of (19) is the attenuation constant Draw a parallel with the classical theory of acoustic attenuation where is the coefficient of attenuation due to viscous losses, and is the coefficient of attenuation due to thermal conduction. From (21) and (22), we can solve for and to model the attenuations such that Because the material layers used in this paper have a low heat conductance, the loss due to thermal conduction is negligible, and we let the conductance .
2.1. Piezoelectric Transducer
An important task when designing ultrasonic transducers and complete transducer system is the simulation of possible configuration prior to construction. Ultrasonic transducer usually consists of a piezoelectric element and nonpiezoelectric layers for encapsulation and acoustic matching. Points of interest are effects of backing materials, matching layers, piezoelectric materials, layer thickness, electrical matching, and coupling on transducer characteristics like bandwidth, time response on different excitation signals, and ring down behavior. As presented in the introduction, different models have been developed over the years to simulate these transducers. Mechanically, a transmission line of length len (m) represents the acoustical layer. The length is selected to achieve the desired center frequency (Hz) of the transducer. With fixed ends, the piezoelectric plate has a fundamental resonant frequency as where is the velocity of sound through it at temperature .
Using (21) and the piezoceramic’s density , required for transmission line, and values can be calculated. The mechanical factor describes the shape of the resonance peak in the frequency domain. The relation between angular frequency , inductance , and the resistance is given as 
In the electrical section, the static capacitance is calculated as where (C2/Nm2) is the permittivity with constant strain .
The latter is related to the permittivity with constant stress (free) as where is the piezoelectric coupling constant.
The mechanical and electrical sections interact with two current controlled sources. From the mechanical side, the deformation itself is not measurable, but the current representing the rate of deformation is the difference between the velocity of each surface normal to the propagation path, represented by the currents and , which is the rate of deformation. This current controls the current source . It has a gain equal to the product of the transmitting constant and the capacitance . is the ratio of the piezoelectric stress constant in the direction of propagation and the permittivity with zero or constant strain . In the thickness mode, it is 
This source’s output is in parallel with the capacitor . The result is a potential difference across the capacitor that is proportional to the deformation. In the electrical section, the current through the capacitor controls the current source . The gain for this second current source is . Its output needs to be integrated to obtain the total charge on the electrodes which proportionally deforms the transducer. The integration is performed by the capacitor . The voltage controlled voltage source with unity gain is a one-way isolation for the integrator.
The speed of sound in a liquid is given as where is the adiabatic bulk modulus .
Starting with the Navier-Stokes equation, Kinsler et al.  derive the following formula for the coefficient of attenuation due to viscosity: where is the viscosity of liquid. One of the assumptions made in this derivation is that the bulk viscosity is negligible. Parameters that are temperature sensitive, like viscosity, need to be obtained at the desired temperature for single temperature simulation or as a function of temperature for parametric simulation over a temperature range.
3. Experimental Setup
The block diagram of pulser-receiver system is shown in Figure 3. The experimental circuit which is designed by Gandole et al.  is used in this set up. The RF pulse generator generates sharp radio frequency pulses of various frequencies in the range 1 MHz to 10 MHz, having pulse width 2 to 60 microseconds. The repetition rate of the pulses is 1 KHz. The radio frequency pulse is fed to the ultrasonic transducer (piezoelectric). The transducer is excited and sends ultrasonic pulse through the sample. The receiver transducer, which is at the other end of the sample, converts the received ultrasonic pulse into an electrical signal (RF pulse). The details of pulser circuit are given in Gandole et al. . This RF pulse is fed to an amplifier consisting of single stage which is assembled by using transistor BF195. The output of single stage amplifier is again amplified with the help of another amplifier using IC CA3028 in cascade mode. The overall characteristics of amplifier are gain = 50 dB; bandwidth = 15 MHz; input impedance = 10.5 kΩ, and low noise. The stability of IC 3028 amplifier is much higher because of small reverse feedback . After the detection of signal through LM 393 (wide band zero cross detector), the detected signal is fed to unity gain buffer, designed using high speed, low power Op-amp AD 826. The characteristics of AD826 are 50 MHz unity gain bandwidth, 350 V/μs slew rate, 70 ns settling time to 0.01%, and 2.0 mV max input offset voltage. This detected signal is then given to reset input of RS flip-flop, while the outputs of IC 74121 of RF pulse generator are used to set the flip-flop. Therefore, a single pulse is obtained at the output, whose width is the time taken by the pulse to travel through the sample. Time meter having accuracy of 0.01 μs measures the width of the pulse. For the measurement of velocity and attenuation, the sender transducer is firmly fixed at one end of the measuring cell, while receiving transducer is fixed to movable scale (Griffin and Tatlock Ltd. London, Make) having least count 0.0001 cm. The system communicates with personal computer (PC) through 12-bit DAS card. For communication, control and set of parameters of ultrasonic system Graphical User Interface (GUI) were created (Figures 5(a)–5(c)). GUI provides general data representation. In GUI it is possible to change almost all of parameters of measured signal and parameters of signal processing. For link-up with system, functions for control drivers in Dynalog system were used. Then parameters for communication with system were set up. In windows (Figure 5), it is possible to select channel number, polarity (unipolar/Bipolar), and range using DAQ configuration frame. “RUN DAQ” frame is used to open the card using “open” command button, close the card using “close” command button. “Start Scan” command button samples the data (voltage level and time) and stores in buffer memory. “Stop Scan” command button stops the sampling process. Measurement and analysis of received ultrasonic signal is achieved by “Measurement and Analysis” frame. The command button “Transient Response” opens the new window, which displays the transient response of the sampled data. Using the cursor point, it is possible to measure the time and pulse height of received signal. This data is stored in database for further analysis.
Matched pair of quartz transducers completely sealed (supplied by Electrosonic Industries, New Delhi.) was used for ultrasonic generation and detection. The frequency of measurement was 2 MHz. Diameter of transducer was 25.4 mm. Also the second set of matched pair of transducer PZT-5A (supplied by Panametrics Videoscan V3456) with center frequency 5 MHz was used. Diameter of transducer was 12.5 mm.
The liquid sample was contained in a measurement cell as shown in Figure 3. A glass bottle of suitable size was cut to form a measuring cell. The bottle was fixed with adhesive araldite in an inner space of double-walled chamber which was made up of thick galvanized sheet. Water circulation arrangement was made through thermostat. Dimensions of double-walled chamber are as height of doubled-walled chamber 6.5 inches, outer diameter 5.5 inches, inner diameter 3.25 inches, height of glass bottle 7.5 inches. The double walled chamber was provided with inlet and outlet for constant temperature water circulation. The lower surface of cell (glass bottle) and double-walled chamber were in the same plane. The double-walled chamber was kept on disc to which sender transducer was fixed. Applying silicone grease made the contact of sender transducer and measuring cell. The lower surface of measuring cell acted as acoustic window through which ultrasonic waves could enter in the measuring cell. The clamps were provided between double-walled chamber and disc to avoid movement of doubled-walled chamber and hence of measuring cell during the ultrasonic measurements.
4. Simulation Setup
The analogous circuit setup is described in Figure 6. Simulation code is given below, to describe the circuit.
Here static capacitance is .
4.1. The Pulser-Receiver
The pulser circuitry, shown in Figure 6, consists of a pulse generator using pulse source having 5 V pulsed voltage, 0 sec delay, 1 ns rise and fall time, 2 μs pulse width and 1 ms period. Another source is sinusoidal source having 0 offset voltage, 5 V peak voltage, 0 sec delay, 5 MHz frequency, 0 damping factor and phase delay followed by TTL7400, 7407, 1 KΩ resistor , 20 V dc supply , and 2nf capacitor . The pulse start at 20 V and decays to 0 V in 2 μs. The output of the capacitor is connected to the electrical part of the transducer and a parallel damping resistor . The receiving and amplifying electronics are not implemented in this schematic.
4.2. The Transducer and Liquid Samples
The matched pair of piezoelectric material PZT-5A, whose material data is obtained from Benny et al. , was chosen and given in Table 1. For simplicity, the front wear plate is omitted, and to reduce the ringing of the echoes a backing material is selected with a characteristic acoustic impedance of 15.8 MPas/m. The thickness of the backing layer should be selected such that no echoes return from it. This fact permits us to model the backing layer with a resistor ( KΩ). Another matched pair of quartz transducers supplied by Electrosonic Industries, New Delhi, was chosen and whose material data is given in Table 1. The data obtained for the lower surface of cell (glass bottle) is given in Table 1. The data in Table 1, for liquid samples was obtained from standard references [15, 16, 18].
5. Graphical User Interface for Simulation Analysis
Figure 7 shows the Graphical User Interface (“UT Measurements”) for easy access to the basic functionality of the modeled ultrasonic system using PSPICEA_D. The results of mixed-signal simulations can then be plotted in the same Probe window with little effort. The newly developed Graphical User Interface heightens the level of clarity, accelerates the speed of navigation, and improves control of information flow enabling end users to take advantage of the systems full potential.
6. Signal Processing
6.1. Time Domain
In this section, the time domain received signal and experimental signal are compared. The simulation signal is shown in Figure 8(a) and the measured signal is shown in Figure 8(b). Figure 8 shows the first 30 μs pulse received by the 5 MHz transducer with an ethanol sample at 25°C.
Figure 9 shows the first 32 μs pulse received by the 5 MHz transducer in methanol at 25°C. Figure 10 shows the first 36 μs pulse received by the 5 MHz transducer in carbon tetra chloride at 25°C. Figure 11 shows the first 30 μs pulse received by the 5 MHz transducer in acetone at 25°C. Figure 12 shows the first 27 μs pulse received by the 5 MHz transducer in benzene at 25°C. Figure 13 shows the first 25 μs pulse received by the 5 MHz transducer with a distilled water sample at 25° C. Table 2 shows the comparison between the literature values and simulation values of ultrasonic velocity and attenuation for the liquid samples. Figures 14(a)–14(e) show the first 50 μs pulse received by the 2 MHz crystal transducer with binary liquid mixtures of methanol and benzene sample at 25°C.
Figures 15(a)–15(e) show the first 50 μs pulse received by the 2 MHz crystal transducer with binary liquid mixtures of methanol and benzene sample at 25°C. Tables 3(a)-3(b) show the comparison between the experimental values and simulation values of ultrasonic velocity and attenuation for the binary liquid mixture samples.
Analogy between acoustic media and transmission lines is reviewed, and an analogous electrical model of an ultrasonic transducer using controlled sources is discussed. A simulation model of a complete ultrasonic system is presented. The received signal from the simulation is compared to that of an actual measurement in the time domain. The comparison of simulated, experimental data clearly shows that temperature and frequency dependencies of parameters of relevance to acoustic wave propagation can be modeled. The feasibility has been demonstrated in an ultrasound transducer setup for material property investigations. Comparisons were made for attenuation and velocity of sound for ethanol, methanol, carbon tetrachloride, acetone, benzene, and distilled water. For these materials, the agreement is good. The simulation tool therefore provides a way to predict the received signal before anything is built. Furthermore, the use of an ultrasonic simulation package allows for the development of the associated electronics to amplify and process.
- M. H. Rashid, SPICE for Circuits and Electronics Using PSPICE, Printice Hall, New Delhi, India, 2nd edition, 2002.
- W. P. Mason, Electromechanical Transducers and Wave Filters, Van Nostrand, New York, NY, USA, 1942.
- M. Redwood, “Transient performance of a piezoelectric transducer,” Journal of the Acoustical Society of America, vol. 33, pp. 527–536, 1961.
- S. A. Morris and C. G. Hutchens, “Implementation of Mason's model on circuit analysis programs,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 33, no. 3, pp. 295–298, 1986.
- R. Krimholtz, D. A. Leedom, and G. L. Matthei, “New equivalent circuits for elementary piezoelectric transducers,” Electronics Letters, vol. 6, no. 12, pp. 398–399, 1970.
- W. M. Leach, “Controlled-source analogous circuits and SPICE models for piezoelectric transducers,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 41, no. 1, pp. 60–66, 1994.
- A. Püttmer, P. Hauptmann, R. Lucklum, O. Krause, and B. Henning, “SPICE model for lossy piezoceramic transducers,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 44, no. 1, pp. 60–66, 1997.
- G. Benny, G. Hayward, and R. Chapman, “Beam profile measurements and simulations for ultrasonic transducers operating in air,” Journal of the Acoustical Society of America, vol. 107, no. 4, pp. 2089–2100, 2000.
- J. L. San Emeterio, A. Ramos, P. T. Sanz, A. Ruíz, and A. Azbaid, “Modeling NDT piezoelectric ultrasonic transmitters,” Ultrasonics, vol. 42, no. 1–9, pp. 277–281, 2004.
- M. Hirsekorn, P. P. Delsanto, N. K. Batra, and P. Matic, “Modelling and simulation of acoustic wave propagation in locally resonant sonic materials,” Ultrasonics, vol. 42, no. 1–9, pp. 231–235, 2004.
- L. E. Kinsler, A. R. Frey, A. B. Coppens, and J. V. Sanders, Fundamentals of Acoustics, John Wiley & Sons, New York, NY, USA, 3rd edition, 1982.
- D. K. Cheng, Field and Wave Electromagnetics, Addison-Wesley, Reading, Mass, USA, 2nd edition, 1989.
- G. S. Kino, Acoustic Waves: Devices, Imaging, and Analog Signal Processing, Prentice Hall, Englewood Cliffs, NJ, USA, 1988.
- A. D. Pierce, Acoustics. An Introduction to Its Physical Principles and Applications, Acoustic Society of America, Woodbury, NY, USA, 1989.
- R. C. Weast, Ed., Handbook of Chemistry and Physics, Chemical Rubber, Cleveland, Ohio, USA, 45th edition, 1964.
- D. E. Gray, American Institute of Physics Handbook, McGraw-Hill, New York, NY, USA, 3rd edition, 1972.
- A. A. Berdyev and B. Khemraev, “Method of investigating the acoustical properties of liquids at frequencies of 300–1000 mc,” Russian Journal of Physical Chemistry, vol. 41, article 1490, 1967.
- W. Schaaff, “Numerical data and functional relationalships in science and technology, new seris group II: atomic and molecular physics,” in Molecular Acoustics, K. H. Hellwege and A. M. Hellweg, Eds., vol. 5, Springer, Berlin, Germany, 1967.
- D. Fennell Evans, J. Thomas, J. A. Nadas, and M. A. Matesich, “The conductance of electrolytes in acetone and in 1-propanol-acetone mixtures at ,” The Journal of Physical Chemistry, vol. 75, no. 11, pp. 1714–1722, 1971.
- Y. B. Gandole, S. P. Yawale, and S. S. Yawale, “Simplified instrumentation for ultrasonic measurements,” Electronic Technical Acoustics, vol. 35, 2005.
- M. G. S. Ali, “Analysis of broadband piezoelectric transducers by discreat tiem model,” Egyptian Journal of Solids, vol. 23, pp. 287–295, 2000.