AC/DC CONVERTER WITH THREE TO SINGLE PHASE MATRIX CONVERTER, FULL-BRIDGE AC/DC CONVERTER AND HF TRANSFORMER
20170229972 · 2017-08-10
Inventors
- Diogo André CERQUEIRA PINTO BEZERRA VARAJÃO (Ponte de Lima, PT)
- Luís Miguel FARIA MIRANDA (Moure BCL, PT)
- Rui Manuel ESTEVES ARAÚJO (Porto, PT)
Cpc classification
H02J3/38
ELECTRICITY
G01R25/005
PHYSICS
H02J3/322
ELECTRICITY
Y02E40/30
GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
H02J3/18
ELECTRICITY
International classification
H02M1/42
ELECTRICITY
G01R25/00
PHYSICS
H02M1/14
ELECTRICITY
H02M1/12
ELECTRICITY
Abstract
The present application presents a solution that intends to solve the problem of providing a single-stage bidirectional power conversion system (PCS) with a controllable power factor and the capability to regulate the current in the DC side. Disclosed is a single-stage, bidirectional and high-frequency isolated PCS, comprising a high-frequency transformer (HFT), a three-phase-to-single-phase matrix converter (MC), a full-bridge (FB) AC to DC converter, and a control system, where the control system outputs are connected to the switches of the MC and the FB converter. Moreover, the PCS output can also form a DC network for energy supply of several devices. This system converts three-phase AC power input from the network into DC power output that can be used for example to charge an energy storage device or supply a direct current distribution system. It is also possible to convert DC power input into AC power output to supply the network.
Claims
1. A single-stage power conversion system, comprising: a high-frequency transformer; a three-phase-to-single-phase matrix converter; a full-bridge AC to DC converter; and a control system, where the control system outputs are connected to the switches of the matrix converter and the full-bridge converter.
2. The power conversion system according to the previous claim, wherein the full-bridge converter comprises a parallel association of two legs, where each leg is composed by two switches in series.
3. The power conversion system according to any of the previous claims, wherein the control system comprises: a peak value input of the three-phase input voltages V.sub.a,b,c; a set-point input φ.sub.i; and a phase-shift input φ.
4. The power conversion system according to any of the previous claims, further comprising: a phase detector module, configured to calculate the angle α.sub.1 of the input voltage space vector of the three-phase AC side; an adder module, configured to calculate a desired phase angle β.sub.i by subtracting angle α.sub.1 with the set-point input φ.sub.i; a sector detector module, configured to calculate a sector K that comprises an input phase angle β.sub.i; a switching space vectors selection module, configured to select a switching space vector SSV.sub.s based on the sector K; an angle normalization module, configured to normalize the phase angle, β.sub.i, as an angle within the sector, θ.sub.i; a duty cycle calculation module, configured to calculate duty cycles d.sub.1 and d.sub.2 for the two general switching states based on angle θ.sub.i; and a command signals generation module, configured to output ON/OFF signals to the matrix converter and the full-bridge converter switches based on the SSV.sub.s, d.sub.1, d.sub.2 and the phase-shift input φ.
5. The power conversion system according to any of the previous claims, further comprising a current controller configured to regulate an output current i.sub.dc imposed by a phase-shift input φ.
6. The power conversion system according to the previous claim, wherein the current controller further comprises an outer control loop configured to set the reference for the output current i.sub.dc based on a comparison between a reference V.sub.dc,ref and a DC side voltage V.sub.dc.
7. The power conversion system according to claim 5, wherein the control system further comprises a control loop configured to regulate the power factor in the AC side based on the power calculation and the comparison between a reference cos(φ.sub.i).sub.,ref and a cos(φ.sub.i).
8. The power conversion system according to any of the claims 1 to 5, further comprising a PQ controller configured to regulate the active and reactive power level based on a power calculation and the comparison with respective reference values.
9. Use of the power conversion system described in any of the previous claims as an inverter for energy harvesting from DC sources.
10. Use of the power conversion system described in any of the claims 1 to 8 as a battery-based stationary energy storage systems.
11. Use of the power conversion system described in any of the claims 1 to 8 as an interface between AC networks and DC networks.
12. Use of the power conversion system described in any of the claims 1 to 8 as an electric vehicle battery charger.
13. Use of the power conversion system described in the previous claim, wherein the charger is: on-board; or off-board.
14. A distributed energy storage system, wherein the distributed energy storage system comprises the power conversion system described in any of the claims 1 to 8.
15. A Vehicle-to-Grid system, wherein the vehicle-to-grid system comprises the power conversion system described in any of the claims 1 to 8.
Description
BRIEF DESCRIPTION OF DRAWINGS
[0050] Without intent to limit the disclosure herein, this application presents attached drawings of illustrated embodiments for an easier understanding.
[0051]
[0052]
[0053]
[0054]
[0055]
[0056]
[0057]
[0058]
[0059]
[0060]
[0061]
[0062]
[0063]
[0064]
[0065]
[0066]
[0067]
MODE(S) FOR CARRYING OUT EMBODIMENTS
[0068] Referring to the drawings, herein are described optional embodiments in more detail, which however are not intended to limit the scope of the present application.
[0069] In this document, the topology for the PCS as well the proposed modulation and control strategies are described.
Topology
[0070] The main circuit comprises an input filter, a three-phase-to-single-phase matrix converter (MC), a high-frequency transformer (HFT), a full-bridge (FB) converter and an output filter, as represented in
[0071] The front-end converter for the proposed PCS consists in a three-phase-to-single-phase MC. This converter is composed by six bidirectional switches that can be built with the configurations presented in
[0072] In the topology of
[0073] The MC works as a current source rectifier (CSR) with the purpose of generating sinusoidal input currents and high-frequency output voltages to be applied in the transformer. Due to the high harmonic content present in the input currents of the MC itself, an input filter is necessary to allow the connection with the grid in compliance with power quality standards. A generalized structure for the input filter is represented in
[0074] In the following analysis is considered that symmetrical, three-phase rotary voltages are supplied to the MC. The dynamic of the input filter is considered negligible. Input line to neutral (g point in
[0075] The transformer is used to provide galvanic isolation between the grid and the battery pack and its parameters are used to optimize the switching losses. Moreover, voltage adaptation is also an important feature when there are a significant difference between the amplitude of the input and the output voltages. Power density increases as the operating frequency of the transformer raise, resulting in a reduction of volume and weight. In order to simplify the explanation, in this analysis the transformer is modeled by primary-referred leakage inductance L, and an ideal transformer with turns ratio n. Naturally, other elements must then be considered for modeling the transformer, such as the magnetizing inductance, the primary and secondary windings resistance, inter-winding capacitance, and the core losses.
[0076] The FB establishes the connection between the secondary of the HFT and the output filter. With the proposed control system, the FB works as a current source in its output, represented by current i.sub.o. DC voltage v.sub.o is held almost constant through a small filter capacitor.
Modulation
[0077] The MC is fed by a voltage source, thus to avoid a short circuit of the input voltages, only one input phase may be connected to each P or N bar in
TABLE-US-00001 TABLE I Feasible switching states of the three-phase-to-single-phase MC. State S.sub.aP S.sub.bP S.sub.cP S.sub.aN S.sub.bN S.sub.cN v.sub.P i.sub.a i.sub.b i.sub.c SSV 1 1 0 0 0 1 0 V.sub.ab i.sub.L −i.sub.L 0 I.sub.1 2 1 0 0 0 0 1 −V.sub.ca i.sub.L 0 −i.sub.L I.sub.2 3 0 1 0 0 0 1 V.sub.bc 0 i.sub.L −i.sub.L I.sub.3 4 0 1 0 1 0 0 −V.sub.ab −i.sub.L i.sub.L 0 I.sub.4 5 0 0 1 1 0 0 V.sub.ca −i.sub.L 0 i.sub.L I.sub.5 6 0 0 1 0 1 0 −V.sub.bc 0 −i.sub.L i.sub.L I.sub.6 7 1 0 0 1 0 0 0 0 0 0 I.sub.7 8 0 1 0 0 1 0 0 0 0 0 I.sub.8 9 0 0 1 0 0 1 0 0 0 0 I.sub.9
[0078] The space vector approach is based on the instantaneous space vector representation of voltages and currents. The current SSV have fixed directions and are π/3 radians evenly spaced in the αβ reference frame. Sector numbers 1 to 6, given by variable K, are assigned for the vector spaces between two adjacent SSV. For example, sector 1 is limited by I.sub.1 and I.sub.2, and corresponds to the space between −π/6 and π/6 radians measured from the real axis ().
[0079] At any given sampling instant, the input current displacement angle φ.sub.i is a set-point of the modulation module. The input line-to-neutral voltage vector
β.sub.i=α.sub.i−φ.sub.i Equation 4
[0080] At any instant of time, the reference current vector will be in a particular sector of the αβ reference frame. Using the same transformation previously described, the reference current space vector is given by equation 5, where the average current of each input phase a, b, c can be defined as equation 6, where I.sub.i is the peak of the average input currents of the MC, and β.sub.i is the angle of the average input current space vector.
[0081] The space vector modulation is based on the selection of space vectors in order to approximate the current reference ī.sub.ref. It will be assumed that the input current space vector is lying in sector 1, but this approach can be extended to the other sectors without loss of generality. The two adjacent nonzero SSV (I.sub.K, K=1, . . . , 6) are selected as represented in
[0082] Vector ī.sub.i,ref is decomposed into components ī.sub.γ and ī.sub.δ along the two adjacent current SSV. For this example, space vectors I.sub.1 and I.sub.2 are used to synthesize ī.sub.γ and ī.sub.δ respectively. This can be expressed by the relations in equation 7.
ī.sub.i,ref=ī.sub.γ+ī.sub.δ
ī.sub.γ=d.sub.γI.sub.1
ī.sub.δ=d.sub.δI.sub.2 Equation 7
[0083] Input current reference vector is synthesized by considering the proper duty cycles d.sub.γ and d.sub.δ. The angle θ.sub.i represents the angle of ī.sub.i,ref within the actual sector, as show is
[0084] The duty cycles for the two general switching states are calculated using the relationships in equation 10, where m.sub.γδ is the basic modulation index that is set by the user or by an external system to adjust the amplitude of ī.sub.i,ref. The zero space vectors (I.sub.z=I.sub.7, I.sub.8 or I.sub.9) are used to complete the modulation period, whose duty cycle is calculated by equation 11.
[0085] Since d.sub.0γδ≧0, d.sub.γ≧0, and d.sub.δ≧0, the maximum basic modulation index is limited to 1, so 0≦m.sub.γδ≦1.
[0086]
[0087] The MC changes the low-frequency input voltage into a high-frequency voltage, v.sub.p. This voltage is applied to the primary of the HFT which provides galvanic isolation to the PCS, and makes it possible to store energy using the transformer leakage inductance, L. The MC applies the current SSV of table I in a square wave fashion. For this analysis, it is also considered that the amplitude of voltages v.sub.p and v.sub.s is constant during one period of modulation, T.sub.s. In practice, this can be approximately achieved if the switching frequency is sufficiently higher compared to ω.sub.1. Table II shows the applied SSV and the corresponding voltage v.sub.p according to the sector K. The resulting voltage v.sub.p that is applied in the HFT is the temporal reference of the modulation. The modulation period of v.sub.p can be divided into five different intervals as represented in
TABLE-US-00002 TABLE II SSV and corresponding voltage v.sub.p applied by the MC according to the sector. Duty Positive pulse Negative pulse Zero Sector cycle SSV v.sub.P SSV v.sub.P SSV 1 d1 I.sub.1 V.sub.ab I.sub.4 −V.sub.ab I.sub.7 d2 I.sub.2 −V.sub.ca I.sub.5 V.sub.ca 2 d1 I.sub.2 −V.sub.ca I.sub.5 V.sub.ca I.sub.9 d2 I.sub.3 V.sub.bc I.sub.6 −V.sub.bc 3 d1 I.sub.3 V.sub.bc I.sub.6 −V.sub.bc I.sub.8 d2 I.sub.4 −V.sub.ab I.sub.1 V.sub.ab 4 d1 I.sub.4 −V.sub.ab I.sub.1 V.sub.ab I.sub.7 d2 I.sub.5 V.sub.ca I.sub.2 −V.sub.ca 5 d1 I.sub.5 V.sub.ca I.sub.2 −V.sub.ca I.sub.9 d2 I.sub.6 −V.sub.bc I.sub.3 V.sub.bc 6 d1 I.sub.6 −V.sub.bc I.sub.3 V.sub.bc I.sub.8 d2 I.sub.1 V.sub.ab I.sub.4 −V.sub.ab
[0088] The FB applies voltage v.sub.s to the secondary of the HFT. The difference between voltages v.sub.s and v.sub.p result in voltage v.sub.L which is applied to the leakage inductance L. Thus, the power transfer between the two voltage sources (MC and FB) can be controlled by the amplitude and the difference of phase between voltages v.sub.s and v.sub.p. The amplitude of v.sub.p is imposed by the input three-phase voltage source, and the amplitude of v.sub.s is dependent of the energy storage device. Thus, power flow is controlled by introducing a phase-shift, φ, between voltages v.sub.s and v.sub.p, changing in this way the amplitude of I.sub.i and consequently the current transferred to the energy storage device. The phase-shift is limited by equation 12, and can be normalized by making
resulting in the limit defined in equation 13.
[0089]
TABLE-US-00003 TABLE III Feasible switching states of the FB converter. State S.sub.1 S.sub.2 S.sub.3 S.sub.4 v.sub.s i.sub.o 11 1 1 0 0 V.sub.o n .Math. i.sub.L 12 0 0 1 1 −V.sub.o −n .Math. i.sub.L 13 1 0 1 0 0 0 14 0 1 0 1 0 0
[0090]
TABLE-US-00004 TABLE IV Intermediate signals routing. Switch Sector S.sub.aP S.sub.bP S.sub.cP S.sub.aN S.sub.bN S.sub.cN 1 X.sub.1 Y.sub.1 Z.sub.1 X.sub.2 Y.sub.2 Z.sub.2 2 Y.sub.2 Z.sub.2 X.sub.2 Y.sub.1 Z.sub.1 X.sub.1 3 Z.sub.1 X.sub.1 Y.sub.1 Z.sub.2 X.sub.2 Y.sub.2 4 X.sub.2 Y.sub.2 Z.sub.2 X.sub.1 Y.sub.1 Z.sub.1 5 Y.sub.1 Z.sub.1 X.sub.1 Y.sub.2 Z.sub.2 X.sub.2 6 Z.sub.2 X.sub.2 Y.sub.2 Z.sub.1 X.sub.1 Y.sub.1
[0091] Using a positive phase-shift the power is transferred from the grid to the energy storage device or dc network. If a negative phase-shift is used then the power is transferred from the energy storage device to the grid.
[0092]
[0093] Considering voltages v.sub.p, v.sub.s and the resulting voltage v.sub.L applied in the inductance L, it is possible to establish a set of equations to analytically represent current i.sub.L. By table I is possible to know how the current i.sub.L is reflected in the MC input. Thus, the average currents in the input of the MC for the sector 1 over one modulation period can be given by equation 16, where ω.sub.S is the modulation frequency in radians per second and is given by ω.sub.S=2π/T.sub.S. In order to modulate the reference input current vector of equation 5, it is necessary to choose duty cycles d.sub.1 and d.sub.2 as equation 17, where m is the modulation index. The modulation period is completed through the duty cycle d.sub.0 that is calculated by equation 18.
[0094] The zero voltage is applied by the MC using the zero SSV (I.sub.z=I.sub.7, I.sub.8 or I.sub.9), and by the FB through states 13 and 14 of table III. Since d.sub.0≧0, d.sub.1≧0, and d.sub.2≧0, the modulation index is limited to 1/√{square root over (2)}, so 0≦m≦1/√{square root over (2)}. Since the amplitude of current I.sub.i can be controlled by the phase-shift φ, the modulation index can be set to its maximum value to not limit the transferred power. In the following description the modulation index will be considered m=1/√{square root over (2)} without loss of generality.
[0095] By averaging the three-phase currents in the input of the MC with these duty cycles over one modulation period, and considering the space vector transformation, it is possible to obtain the average value of I.sub.i as defined in equation 19.
[0096] This result is also valid for the other sectors in the αβ reference frame. Considering equation 2, equation 5, and equation 19, the active and reactive power is given by equation 20, using the well-known definitions for active power (P) and reactive power (Q).
[0097] Through table III is possible to know how current i.sub.L is reflected in the FB output. Assuming that the power losses in the conversion system are zero, the average current in the output of the converter over one modulation period can be expressed as defined in equation 21.
[0098] With this result is possible to determine the average output power of the PCS as defined in equation 22.
[0099] As expected, the input active power (equation 20) is equal to the output power (equation 22) since it was considered that the power losses during the conversion process are zero. Moreover, the modulation executed by the control system allows to control simultaneously the power factor in the AC side interface through the variable φ.sub.i (equation 20), and the current regulation in the DC side through phase-shift φ (equation 21).
[0100] With these control variables there are flexibility to implement several control strategies for the PCS. One possibility consist in having a current controller that imposes the proper phase-shift in order to regulate the output current i.sub.dc, as shown in
[0101] If the objective is to regulate the output voltage v.sub.dc, an outer control loop can be added to the previous strategy, as shown in
[0102]
[0103] The combination of the developed control system with the proposed main circuit results in a single-stage PCS applied, for example, to energy storage devices. The PCS, when compared with conventional solutions, has as advantages a higher power density, longer service life, higher power quality in the grid interconnection, power factor control capability, safe operation provided by galvanic isolation, wide DC voltage range of energy storage devices, and DC current regulation with low current ripple.
Examples
[0104] The disclosed technology can be applied in diverse PCSs.
[0105] As unidirectional power converter, this technology can be used in both—AC to DC or DC to AC conversion. Major DC to AC converter applications are based on grid-tied or standalone inverters for energy harvesting from renewable or other direct current sources, requiring galvanic isolation. In AC to DC conversion, the proposed technology can be used for conventional battery chargers and power supplies with galvanic isolation. Solar or wind power inverters, battery chargers for industrial equipment, high power DC sources for industrial processes, and energy supply in DC for electronic equipment, datacenters, variable frequency drives, etc., are examples of possible applications.
[0106] As bidirectional power converter it can be used for battery-based stationary energy storage systems, electric vehicle battery chargers with reverse power flow functionality, energy supply of DC loads with regenerative capabilities, interface between low voltage AC and low voltage DC networks.
[0107] A particular application is in distributed energy storage systems under the so called “Community Energy Storage” concept. The need to integrate renewable based energy sources and heavy loads such as electric vehicles, in the low voltage grid is demanding new levels of flexibility and reliability. As a relatively new concept, the community energy storage can be a key solution to mitigate the technical impacts of this large integration and has been a solution that leads to high expectations, particularly from distribution system operators. Some advanced functionalities can be provided by these systems, such as voltage support, frequency regulation, peak shaving and load leveling, buffer and smooth renewable energy variability, grid upgrade deferral or even power factor correction and harmonic filtering.
[0108] Another particular example of application can be the electric vehicles market exploring the so called “Vehicle-to-Grid” (V2G) concept. The electric vehicles are typically parked 95% of the time. Since they have energy stored in its batteries, some of this energy can be used to support the electric grid by supplying power at some demanding periods of the day. A bidirectional charger is necessary to supply the batteries during charging operation, and to discharge the batteries during V2G operation. Due to space limitation in the electric vehicle, the PCS disclosed has a great potential for V2G applications since can result in a compact converter. Currently, the available chargers in the market don't have this functionality, but is expected to be a common solution in the future, particularly in the smart home paradigm.
[0109] Additional application examples are the DC supply networks, which will become a tendency for the near future. This type of solution has been demonstrated in large computational systems, such as datacenters, and in modern networks with combined production of renewable energy with energy storage. In all these cases there have been significant gains in energy efficiency. Being power electronic converters the key enabler of this new type of networks, the proposed technology can be applied in the implementation of a bidirectional power interface between the existing AC low voltage network and the emerging DC networks.
[0110] In the case of the variable frequency drives, a major industrial market for power electronics, this solution can be applied to supply the DC bus of large groups of devices in low to mid power range. This way, the regenerative braking will become economically feasible in this power range, since the drive units itself are simpler and cheaper, and the regenerative braking functionality lies in a central bidirectional converter.
[0111] Naturally, the present embodiments are not in any way limited to the embodiments described in this document and a person with average knowledge in the field will be able to predict many possible changes to it without deviating from the main idea, as described in the claims.
REFERENCES
[0112] [1] S. Vazquez, S. M. Lukic, E. Galvan, L. G. Franquelo, and J. M. Carrasco, “Energy Storage Systems for Transport and Grid Applications,” Industrial Electronics, IEEE Transactions on, vol. 57, pp. 3881-3895, 2010. [0113] [2] Y. Gao, M. Li, P. Kuang, W. Wang, and J. Chen, “Battery charge and discharge device,” China Patent CN101604853, 2009. [0114] [3] M. N. Kheraluwala, R. W. Gascoigne, D. M. Divan, and E. D. Baumann, “Performance characterization of a high-power dual active bridge,” Industry Applications, IEEE Transactions on, vol. 28, pp. 1294-1301, 1992. [0115] [4] C. Ga-Gang, J. Doo-Yong, C. Sung-Chop, W. Chung-Yuen, J. Yong-Chae, and Y. Jang-Hyoun, “10 kW rapid-charger for electric vehicle considering vehicle to grid(V2G),” in Power Electronics and Motion Control Conference (IPEMC), 2012 7th International, 2012, pp. 2668-2672. [0116] [5] “Power Conversion Device and Charging device,” Japan Patent JP2013021867, 2013. [0117] [6] Y. Zhaoyang, J. Minli, Z. Chunjiang, and W. Weiyang, “An Integration SPWM Strategy for High-Frequency Link Matrix Converter With Adaptive Commutation in One Step Based on De-Re-Coupling Idea,” Industrial Electronics, IEEE Transactions on, vol. 59, pp. 116-128, 2012. [0118] [7] S. Ratanapanachote, C. Han Ju, and P. N. Enjeti, “A digitally controlled switch mode power supply based on matrix converter,” Power Electronics, IEEE Transactions on, vol. 21, pp. 124-130, 2006. [0119] [8] Y. Wei, W. Wu, C. Yan, D. Zhang, W. Zhang, Z. Zhang, et al., “Synchronous and frequency-doubling jagged carrier generating circuit of high-frequency link matrix type inverter,” China Patent CN101702588, 2010. [0120] [9] D. Sha, Z. Qin, D. Wu, and X. Liao, “A Digitally Controlled Three-Phase Cycloconverter Type High Frequency AC Link Inverter Using Space Vector Modulation,” Journal of Power Electronics (KIPE), vol. 11, pp. 28-36, 2011. [0121] [10] D. Sha, Z. Qin, and X. Liao, “Energy storage system and control method thereof,” China Patent CN101944745, 2011. [0122] [11] D. Sha, Z. Qin, D. Wu, and X. Liao, “High-frequency isolated three-phase cycloconverter type two-way converter and control method thereof,” China Patent CN102075109, 2011. [0123] [12] N. D. Weise, K. Basu, and N. Mohan, “Advanced modulation strategy for a three-phase AC-DC dual active bridge for V2G,” in Vehicle Power and Propulsion Conference (VPPC), 2011 IEEE, 2011, pp. 1-6. [0124] [13] R. K. Gupta, K. Mohapatra, N. Mohan, G. Castelino, K. Basu, and W. Nathan, “Soft Switching Power Electronic Transformer,” United States of America Patent US2011007534, 2011. [0125] [14] H. S. Krishnamoorthy, P. Garg, and P. N. Enjeti, “A matrix converter-based topology for high power electric vehicle battery charging and V2G application,” in IECON 2012-38th Annual Conference on IEEE Industrial Electronics Society, 2012, pp. 2866-2871. [0126] [15] P. W. Wheeler, J. Rodriguez, J. C. Clare, L. Empringham, and A. Weinstein, “Matrix converters: a technology review,” Industrial Electronics, IEEE Transactions on, vol. 49, pp. 276-288, 2002. [0127] [16] A. Trentin, P. Zanchetta, J. Clare, and P. Wheeler, “Automated Optimal Design of Input Filters for Direct AC/AC Matrix Converters,” Industrial Electronics, IEEE Transactions on, vol. 59, pp. 2811-2823, 2012. [0128] [17] L. Huber and D. Borojevic, “Space vector modulated three-phase to three-phase matrix converter with input power factor correction,” Industry Applications, IEEE Transactions on, vol. 31, pp. 1234-1246, 1995.