An analytical approach to the dynamics of rotating shafts

An analytical procedure, based on the dynamic stiffness method, is proposed for the study of rotor dynamics problems. On the grounds of the governing differential equations of a continuous beam, the dynamic stiffness matrix of the rotating Timoshenko beam is derived, including the effects of translational and rotational inertia, gyroscopic moments and shear deformation of the shaft. Concentrated disks and isotropic, elastic bearings are taken into account in the element formulation. The results obtained by the proposed method are compared with both classical closed form solutions and finite element analyses taken from the literature.

---

Sommario Nel presente lavoro viene esposto un procedimento analitico, basato sul metodo della rigidezza dinamica, per lo studio del comportamento dinamico dei rotori. In particolare, è ricavata la matrice di rigidezza dinamica della trave di Timoshenko rotante, che consente di tenere conto degli effetti dell'inerzia alla traslazione e alla rotazione, del momento giroscopico distribuito e della deformazione a taglio della trave. Nell'ambito della formulazione proposta, è possibile considerare, mediante inserimento diretto nella matrice di rigidezza dinamica, gli effetti dovuti alla presenza di dischi e di supporti elastici isotropi. Il metodo proposto è applicato alla determinazione delle velocità critiche di sistemi rotanti, con riferimento sia a casi semplici, per i quali è nota la soluzione analitica, sia al caso piu' complesso di un rotore già studiato in letteratura; per quest'ultimo esempio, in particolare, è presentato un confronto con i risultati delle analisi agli elementi finiti svolte da altri autori.

     

A numerical solution of the nonlinear wave resistance problem for simple shaped submerged bodies

Summary This paper presents a numerical method to compute the nonlinear wave resistance of a body, submerged in a free stream of infinite depth. The solution is given in the form of a potential of single layer and the nonlinear boundary conditions on the free surface are fulfilled by means of an iterative procedure. Results are shown for the 2-D case of an elliptical cylinder and for the 3-D case of prolate spheroids. All the presented results are discussed taking into account some experimental data and analytical solutions available in literature.

---

Sommario In questo lavoro viene sviluppato un metodo per il calcolo della resistenza ondosa di corpi immersi bidimensionali e tridimensionali, tenendo conto della non linearità delle condizioni sulla superficie libera. La soluzione cercata è posta nella forma di potenziale di strato semplice e la non linearità presente nelle condizioni al contorno viene trattata con un procedimento iterativo. Vengono presentati dei risultati numerici relativi al caso bidimensionale di un cilindro ellittico ed al caso tridimensionale di sferoidi prolati di diversa eccentricità. I risultati vengono discussi in relazione a soluzioni analitiche e dati sperimentali ottenuti dalla letteratura.

     

Rational criteria for quantifying the performance of mechanical presses

Summary In the practical mechanics a great interest have for the mechanical engineers the questions concerning a correct definition of the performance of a mechanical press and the theoretical and practical methods to determine its allowable limit in service.After a preliminary study on the parameters really significant, this paper shows an analytic method to get, in quantitative form, the solution of the problems concerning the determination of the duty conditions compatible with the regular operation of a mechanical press.As a conclusion, some general relations are proposed in order to describe quantitatively the performance of a mechanical press and to give a rational form to the set of mechanical quantities characterising it.

---

Sommario Nel campo della ingegneria meccanica, un grande interesse rivestono le questioni concernenti la corretta definizione delle prestazioni di una pressa e i metodi teorici e pratici per determinare i corrispondenti limiti ammissibili di esercizio.Dopo uno studio preliminare sui parametri realmente significativi, viene esposto un metodo analitico per ottenere la soluzione in forma quantitativa di problemi concernenti con la determinazione delle condizioni di esercizio compatibili con il regolare funzionamento della pressa.Come conclusione, si propongono alcune generali relazioni per la formulazione quantitativa delle prestazioni delle presse meccaniche e per descrivere in forma razionale le caratteristiche operative e di impiego di detta macchina.

     

Elastic stability of rectangular shells in membrane equilibrium

Summary This study furnishes an algorithm for determining, in accordance with Dirichelt's (energetic) principle, the critical load of rectangular membranes of any shape. The algorithm is formulated with reference to an extrinsic system of orthogonal Cartesian coordinates, that is the system used in statics for Pücher's method. This involved obtaining the expressions of the deformation characteristics and of the curvature variations of the shell up to the second order in the displacement components. These quantities yielded the final expression of the deformation energy and the second order work of internal stresses. In view of the applications, possible displacement terns are proposed characterizing the unstability secondary buckled lines, which are further detailed for the rectangular plan paraboloid. Finally, some numerical applications are given for the square plan shell.

---

Sommario Col presente studio viene fornito un appropriato algoritmo per la determinazione, secondo il principio energetico del Dirichlet, dei moltiplicatori critici per le membrane di forma qualsiasi su pianta rettangolare. La formulazione è fatta con riferimento ad un sistema di coordinate cartesiane ortogonali estrinseco e cioè lo stesso adottato per il metodo di Pücher nella statica. Allo scopo è stato necessario ricavare le espressioni delle caratteristiche della deformazione e delle variazioni di curvatura della volta sino al secondo ordine nelle componenti dello spostamento. Tali quantità hanno permesso di dare in forma esplicita l'espressione dell'energia di deformazione ed il lavoro del secondo ordine degli sforzi interni. In vista delle applicazioni sono proposte delle possibili terne di spostamenti caratterizzanti le deformate secondarie instabilizzanti, particolarizzandole infine per il caso del paraboloide su pianta rettangolare.Si riportano infine alcune applicazioni numeriche relative alla volta su pianta quadrata.

     

Tensorial formulation of neutron diffusion theory

Summary In this paper a tensorial formulation of monoenergetic neutron diffusion theory is presented as can be derived starting from the integral form of the Boltzmann equation. A result of this tensorial approach is that, as a consequence of a spatial variation of the macroscopic cross sections or of the finite dimensions of the body under examination, the diffusion coefficient is no longer a constant scalar quantity, but a second order symmetric tensor, whose components are functions of both the position and the properties of the material the body is made of. The components of the diffusion coefficient tensor are explicitly evaluated in the case of some homogeneous convex bodies.

---

Sommario In questa nota viene presentata una formulazione tensoriale della teoria della diffusione di neutroni monoenergetici, quale può essere dedotta a partire dalla forma integrale dell'equazione di Boltzmann. Nell'ambito di questa formulazione si precisa che il coefficiente di diffusione, in conseguenza di una variazione regolare delle sezioni d'urto macroscopiche con la posizione o delle dimensioni finite del mezzo in esame, non è più una quantità scalare costante, ma un tensore simmetrico del secondo ordine, le cui componenti sono funzioni della posizione e delle proprietà del materiale. Queste componenti vengono calcolate esplicitamente per alcuni mezzi omogenei convessi.

     

Shakedown analysis of limited-ductility structures

Summary The paper deals with the shakedown analysis problem of discrete elastoplastic structures, subjected to variable and repeated quasi-static loads, account being given to the limited structural ductility by imposing a-priori limit values to some relevant parameters of plastic deformation. The key-concept of limited-ductility load amplifier is assessed and the problem is given a unified formulation as a nonlinear mathematical programming problem, where ductility constraints — of generalized character — are obtained by using a so-called perturbation method which leads to a-priori upper bounds of the relevant plastic deformation measures. Lower bounds to the load amplifier are so obtained and three different computationals stratagies are herein outlined: the linear parametric programming one is applied to three case-study structures.

---

Sommario La memoria affronta lo studio del problema di adattamento elastoplastico di strutture discrete, sottoposte a carichi comunque variabili in modo quasi-statico, tenendo conto di un comportamento a limitata duttilità della struttura, il quale richiede che uno o più parametri della deformazione plastica risultino non superiori di assegnati valori. Del problema di analisi posto, utilizzando il concetto di moltiplicatore dei carichi per l'adattamento con limitata duttilità, viene data una formulazione unitaria come problema di programmazione matematica non lineare, i cui vincoli di duttilità sono ottenuti con tecniche perturbative di maggiorazione a-priori delle deformazioni plastiche. Il metodo conduce a limiti inferiori del moltiplicatore dei carichi dipendenti dal livello assegnato alla deformazione plastica. Tre differenti strategie computazionali vengono delineate per la soluzione del problema, una delle quali di tipo programmazione lineare parametrica viene applicata a tre esempi.

---

---

Ricerca finanziata dal CNR e dal MPI.     

Accurate solutions and stability criterion for periodic oscillations in hysteretic systems

Summary The study of the periodic response of hysteretic oscillators is reduced to that of nonlinear elastic oscillators by assuming an incremental formulation for the constitutive relationship. The harmonic balance method with many components allows for accurate periodic solution computation. The Floquet theory can be used to check stability. Developed equations are applied to the study of frequency response curves of a hysteretic oscillator that, although simple, shows both degrading and non degrading behaviour. The results reported clearly show the shortcomings of traditional methods; the influence of higher harmonics is far from negligible.

---

Sommario Nello studio della risposta periodica degli oscillatori isteretici si incontrano difficoltà dovute alla natura non olonoma del legame costitutivo. Nel presente lavoro queste difficoltà vengono superate in parte adottando una formulazione incrementale che riconduce il caso isteretico a quello elastico non lineare seppure al prezzo di un allargamento dello spazio delle variabili di stato. Le soluzioni periodiche sono determinate seguendo il metodo del bilancio armonico con molte componenti; la stabilità viene analizzata con la teoria di Floquet. Le equazioni sviluppate sono utilizzate per determinare le curve di risposta in frequenza sotto una forzante sinusoidale, di un oscillatore che, benchè semplice, presenta una notevole varietà di comportamento. I risultati ottenuti mostrano chiaramente le lacune dei metodi tradizionali; l'influenza delle armoniche superiori in molti casi è ben lungi dall'essere trascurabile.

     

On the general solution of stress equation for the membrane theory of shells of revolution

Summary A new method is presented to solve the equation of Neményi-Truesdell, in the membrane theory of thin shells of revolution. It is an ordinary linear differential equation of the second order, for which resolution methods of Truesdell and Mitrinovitch are known. The proposed method is based on the transformation of the given equation into a linear integral equation of Volterra type. This allows to write the general solution of stress equation only in terms of equation of the meridian generating the shell.

---

Sommario Si presenta un nuovo metodo per risolvere l'equazione di Neményi-Truesdell, nella teoria delle membrane dei gusci sottili di rivoluzione. Si tratta di un'equazione differenziale lineare ordinaria del secondo ordine, della quale sono noti metodi risolutivi dovuti a Truesdell ed a Mitrinovitch. Il metodo proposto è basato sulla trasformazione dell'equazione data in un'equazione integrale lineare di Volterra e consente di esprimere l'integrale generale dell'equazione degli sforzi esclusivamente in termini della equazione del meridiano generatore del guscio.

     

A differential equation for the minors of the state-transition matrix of linear systems

Summary A differential equation for the minors of the state transition matrix of linear time-invariant systems is established.An example is then presented, where the differential equation is employed in order to check the existence of the extremal control of a linear system with a performance index given by the integral of a quadratic form.

---

Sommario Si stabilisce un'equazione differenziale per i minori della matrice di transizione dei sistemi lineari stazionari.Si illustra quindi un esempio, in cui tale equazione differenziale viene impiegata per verificare l'esistenza del controllo estremale di un sistema lineare con indice di comportamento dato dall'integrale di una forma quadratica.

     

On the formal structure of continuum mechanics Part I: Deformation theory

Summary The purpose of this and of following papers is to give an operator form to the equations of statics of continuous media in the case of small deformations (linear case) and of large deformations (non linear), aiming to obtain a unitary view of many known results in linear and nonlinear mechanics. In this paper an operator version of compatibility equations for large deformations, is given introducing a nonlinear operator that will be shown to be of the potential type. Some new results follow, namely a variational formulation of compatibility conditions extending a similar result obtained by the author in the small deformation theory. This formulation in three-dimensional space is similar to that of Einstein describing the gravitational field in space-time. An operational version of Bianchi's identities and a construction of boundary compatibility conditions are also given.

---

Sommario In questa nota e in altre che seguiranno ci si propone di dare una versione operatoriale delle equazioni della statica dei continui deformabili sia a piccole deformazioni (caso lineare) che a deformazioni finite (non lineare), al fine di ottenere una visione unitaria di molti risultati noti nella meccanica lineare e non lineare. In questa nota si dà forma operatoriale alle condizioni di congruenza per deformazioni finite, introducendo un operatore nonlineare che sarà mostrato essere di tipo potenziale. Di qui seguono alcuni risultati nuovi: la formulazione variazionale delle condizioni di compatibilità che estende un risultato precedente dell'autore nella teoria delle piccole deformazioni. Questa formulazione nello spazio geometrico tridimensionale è del tutto identica a quella di Einstein che descrive il campo gravitazionale nello spazio-tempo. Si da altresì interpretazione operazionale alle identità di Bianchi e una costruzione delle condizioni di congruenza al contorno.

---

---

This work has been supported by C.N.R.     

Higher Order Alongwind Response Statistics of Linear SDOF Structures Including the Quadratic Term of the Fluctuating Wind Velocity

Abstract. A method evaluating the response correlations of a SDOF system subjected to the wind action, once that the first order statistics and the transition matrix of the system are known, is presented. The method does not require further integrations and allows to take into account the presence of the quadratic term of the velocity in the expression of the wind action. This fact makes the response a non-Gaussian process as it really is. The application of the method to some examples shows both the good accuracy of the results, if compared with those obtained by Monte Carlo simulations, and that neglecting the quadratic term may lead to relevant errors even for the correlations, as well as for the first order statistics.Sommario. In questo lavoro viene presentato un metodo che consente la valutazione delle funzioni di correlazione di qualunque ordine della risposta di un sistema strutturale soggetto all'azione del vento. Tale metodo richiede la conoscenza delle caratteristiche statistiche del primo ordine della risposta e della matrice di transizione. Il metodo non richiede ulteriori integrazioni rispetto a quelle eventualmente necessarie per determinare le statistiche del primo ordine, e consente di tenere conto della presenza del termine quadratico della parte fluttuante della velocità del vento nell'espressione dell'eccitazione sulla struttura. Ciòrende la risposta strutturale non Gaussiana, come realmente è L'applicazione del metodo ad alcuni esempi mostra, da un lato, la buona accuratezza dei risultati ottenuti, confrontandoli con quelli ottenuti tramite simulazioni alla Monte Carlo, e dall'altro che trascurare il termine quadratico nella forzante può condurre ad errori rilevanti nella valutazione delle correlazioni della risposta, coscome era stato mostrato in altri lavori per le statistiche del primo ordine.     

Hemivariational inequalities for linear and non linear elastic materials

Summary In the present paper we study hemivariational inequalities arising in linear and nonlinear elastic bodies which are subjected to boundary conditions of nonmonotone multivalued type. After the formulation of the boundary conditions and the derivation of the variational principles we study a hemivariational inequality for a plane elastic body concerning the approximation of its solution. The proposed method is illustrated by an example dealing with the adhesive contact problem of two cylinders and another one concerning an adhesively supported beam.

---

Sommario Nel presente lavoro si studiano diseguaglianze semivariazionali che nascono in relazione a corpi elastici lineari e non lineari, sottoposti a condizioni al contorno di tipo a più valori non monotono. Dopo la formulazione delle condizioni al contorno e la derivazione dei principi variazionali si studia una diseguaglianza semivariazionale per un corpo elastico piano, con particolare riguardo per l'approssimazione della sua soluzione. Il metodo proposto viene illustrato con un esempio che tratta il problema del contatto con adesione di due cilindri.

---

---

On leave from School of Technology. Aristotle University Thessaloniki, Greece. Fulbright Visiting Scholar.     

The diffraction of a plane wave by an isolated breakwater

Summary The diffraction of a plane wave incident on an isolated breakwater is here studied: the exact solution of the problem is briefly reported and a general method involving energies is used to determine comparative importance of the terms of the series which appear in the solution, which becomes of much greater facility in use in practical problems. Numerical calculations have been carried out for twelve different cases, with the wavelength of the incident wave comparable to the length of the breakwater. Results are graphically shown in a set of tables.

---

Sommario Si studia la diffrazione di un'onda piana che investe una diga foranea isolata: si riporta brevemente la soluzione esatta del problema e si utilizza un metodo energetico generale per la determinazione dei pesi relativi dei termini delle serie che compaiono nella soluzione, che viene resa di molto più semplice impiego per i problemi pratici. Si sono effettuati calcoli numerici relativi a dodici casi diversi, con lunghezze d'onda incidente paragonabili alla lunghezza della diga. I risultati sono stati tradotti graficamente in una serie di tavole.

     

Uniqueness of the concentrated-load problem in the linear theory of couple-stress elasticity

Summary It is pointed out that there exist at least two different solutions of the problem of concentrated loads in the two-dimensional, linear couple-stress theory when the formulation is based on the usual uniqueness theorem. An extension of this uniqueness theorem is proved. A set of conditions sufficient for uniqueness is found and is used in a formulation of the concentrated load problem which results in a unique solution. The significant new condition is that the order of the stress singularity is limited to O(r^–1), where r is the distance from the concentrated load.

---

Sommario Si fa notare che esistono almeno due soluzioni diverse del problema dei carichi concentrati nella teoria lineare, a due dimensioni, delle coppie di volume quando la formulazione è basata sul teorema di unicità.Si dimostra una estensione di questo teorema di unicità. Si trova un gruppo di condizioni sufficienti per l'unicità; queste condizioni vengono usate nella formulazione del problema del carico concentrato che dà luogo ad un'unica soluzione.La nuova condizione significativa è che l'ordine della singolarità dello sforzo è limitato a O(r^–1), dove r è la distanza dal carico concentrato.

---

---

This work is a result of research sponsored by the Office of Naval Research, U.S. Navy, under Contract Nonr-610(06).     

Separation points in laminar flows

This paper presents an approximate method for determining separation points in laminar flows. The method of solution of the boundary layer equations uses two integrated equations: the main characteristic of the method is that its approximations (based on a representation of the velocity profiles by means of the Taylor formula) are only in integrals. The first-order approximation gives the solution in closed form and is accurate for accelerated flows and for decelerated flows when the non-dimensional axial derivative of the pressurep _x is less than 0.5: for stronger adverse pressure gradients accurate results are obtained by means of the third approximation.The same method was applied successfully to the energy equation.

---

Sommario Questo lavoro presenta un metodo approssimato per determinare i punti di separazione in flussi laminari. Il metodo di soluzione delle equazioni dello strato limite utilizza le derivate dell'equazione della quantita' di moto e due equazioni integrate: la caratteristica fondamentale del metodo risiede nell'approssimare (mediante una rappresentazione dei profili di velocita' basata sulla formula di Taylor) solo integrali. La prima approssimazione da' la soluzione in forma chiusa ed e' accurata per flussi accelerati e per flussi decelerati quando la derivata della pressione assiale in forma adimensionalep _x e' minore di 0.5: per gradienti di pressione piu' sfavorevoli risultati accurati si hanno con la terza approssimazione. Lo stesso metodo e' stato applicato con successo all'equazione dell'energia.

     

The differential geometry of bending

Summary The deformation associated with uniform bending has been treated using the techniques of differential geometry. The Burgers circuit as well as various related tensor quantities such as torsion and dislocation density have all been determinated with respect to these various states of deformation.

---

Sommario La deformazione associata con la flessione uniforme viene trattata usando le tecniche della geometria differenziale. Il circuito di Burgers e altre quantità tensoriali e la densità di dislocazione, vengono determinate nei vari stati di deformazione.

     

Weighted residuals as a basis of a general solution method in elasticity

Summary It is shown that the general integral form of the elastic equilibrium equations obtainable through the weighted residuals agrees with the variational formulation given by the extremum conditions of the Washizu functional allowing a complete relaxation of the interelement continuity requirements.

---

Sommario Si mostra che la forma generate integrale del sistema di equazioni dell'equilibrio elastico a cui conduce il metodo dei residui con fattori ponderali concorda con la formulazione variazionale fornita dalle condizioni di estremo del funzionale di Washizu che deriva da quello di Hellinger-Reissner mediante sblocco completo dei vincoli di continuità fra regioni contigue.

     

Bias in modal parameters due to aliasing in the time domain

This paper deals with some aliasing effects in the time domain that can lead to unacceptable misestimations of modal parameters. When a frequency response function of a vibrating system is sampled and inverse Fourier transformed, the resulting impulse response is given by an infinite geometric series, the single term of which is the impulse response itself shifted in time. For this reason, some modal parameters, if estimated in the time domain, are biased; in particular, while the damping factor and the natural frequency are not influenced by the aliasing phenomenon, the magnitude and phase of the residue can be highly biased. Corrective terms are theoretically evaluated and their efficiency is shown in numerical simulations.

---

Sommario In questo articolo e' presentato un metodo per correggere gli errori che si compiono nella stima di alcuni parametri modali, quando essi vengono ricavati nel dominio del tempo. Infatti se la funzione di risposta in frequenza e' ottenuta con eccitazioni particolari—quali ad esempio lo pseudo-random, lo stepped-sine a passo costante o lo sweep in frequenza-la risposta impulsiva, ottenuta per mezzo dell'antitrasformata discreta di Fourier, risulta periodica, con periodo pari all'inverso della spaziatura tra le righe spettrali. Cio' comporta un errore nella stima dell'ampiezza e della fase del residuo, mentre nessuna conseguenza si ha sulla frequenza naturale e sul fattore di smorzamento. L'errore sulle stime risulta tanto piu' grande, quanto meno smorzata e' la risposta impulsiva all'interno dell'intervallo di osservazione.Partendo dall'inviluppo e dalla fase istantanea dei segnali complessi, ottenibili per mezzo della trasformata di Hilbert, si sono ricavati i termini correttivi sia per il modulo, che per la fase del residuo. La validita' delle correzioni e' mostrata con esempi numerici.

     

Theoretical studies on heterogeneous deflagration waves. 1. A partial differential equation formulation of the problem

Summary This is the first of a three — part theoretical study on heterogeneous deflagration wave propagation. The theory is of a general nature; but specific reference to a composite solid rocket propellant is made. The ultimate objective of this line of research is to define conditions for statically and dynamically stable deflagration propagation. In this first paper, a quite general model of the problem in terms of a partial differential equation is reviewed. A flame model is used for the gas phase. The important assumptions are quasi-steady gas phase and collapsed reacting layer in the condensed phase. The model developed is the starting point for a transformation into an ordinary differential equation formulation of the problem (second part of the study); from this, the stability features of heterogeneous deflagration waves are immediately defined (last part of the study). The theory is verified by computer and experimental work, presently under progress.

---

Sommario Questa è la prima di una serie di tre memorie dedicate all'analisi teorica della propagazione di onde eterogenee di deflagrazione. La teoria è di natura generale, ma per concretezza viene applicata al caso specifico di un propellente solido composito per endoreattori. Lo scopo finale della ricerca è di definire le condizioni necessarie per la stabilità statica e dinamica della propagazione di onde eterogenee di deflagrazione. In questa prima memoria viene rivisto un modello abbastanza generale del problema basato su di una equazione differenziale alle derivate parziali. La fase gassosa è trattata secondo un modello di fiamma. Le ipotesi più importanti riguardano la quasi-stazionarietà della fase gassosa e lo spessore trascurabile dello strato reagente in fase condensata. Il modello cosi sviluppato costituisce il punto di partenza per una formulazione del problema basata su di una equazione differenziale ordinaria (seconda memoria), che permette di definire immediatamente le proprietà di stabilità di onde eterogenee di deflagrazione (terza memoria). Sono in corso di svolgimento le verifiche numeriche e sperimentali della teoria proposta.

---

Support by CNPM (Centro di Studio per Ricerche sulla Propulsione e sull'Energetica) is gratefully aknowledged.     

Some remarks on neutron transport in homogeneous slabs and spheres

Summary In this work we are concerned with the stationary neutron transport Boltzmann equation (in its integral form) in plane and spherical symmetry; both the study of the properties of the exact and approximate solution are based on functional analysis. We have taken advantage of the fact that the kernel of the integral equation is even, to semplify the general discussion; among the results, we quote that relating to the continuity of the solution with respect to the optical thickness. Moreover we emphasized that the Ritz method seems the most convenient one for the approximate evaluation of the solution; in this connexion, some observations have been made about the convergence of the approximate solutions to the exact one.

---

Sommario In questo lavoro consideriamo l'equazione stazionaria di Boltzmann (nella forma integrale) in simmetria piana e sferica; l'uso di alcuni metodi dell'analisi funzionale ci consente di studiare sia le proprietà della soluzione esatta che di quella approssimata. La semplice osservazione che il nucleo dell'equazione integrale è pari, ci ha permesso di semplificare la discussione generale; tra i risultati segnaliamo quello riguardante la continuità della soluzione rispetto allo spessore ottico. Inoltre abbiamo sottolineato che il metodo di Ritz ci sembra il più opportuno punto di partenza per il calcolo approssimato della soluzione; in relazione a ciò, vengono fatte alcune osservazioni a proposito della convergenza delle soluzioni approssimate verso quella esatta.

---

---

Work performed under contract C.N.R.