Analysis of static and dynamic behaviour of machine-tool structure: Development of numerical method and computer program

Summary A number of machine-tool structures and subassemblies may be approximated by a set of parallel and orthogonal beam elements suitably interconnected. The recognized advantages offered by such a family of models for solving the static and dynamic analysis problems stimulated the development of the computer program described further in the work. This work is made up by three main parts, namely, the first dealing with theoretical foundations of the equilibrium method, the second with the development of numerical solution method of the fundamental matrix equation, and eventually the third in which program structure, input and output are described and a simple example of application is given.

---

Sommario Tenuto conto che un gran numero di macchine utensili e/o di componenti di queste possono schematizzarsi mediante una serie di elementi trabeiformi fra loro paralleli od ortogonali e rigidamente collegati e che il calcolo automatico del comportamento statico e dinamico di tali modelli non comporta costi troppo elevati assicurando risultati soddisfacenti, è stato sviluppato il programma di calcolo descritto nel lavoro in oggetto. Questo lavoro è costituito da tre parti fondamentali, la primq tratta dei fondamenti teorici del metodo di “equilibrio”, la seconda sviluppa una soluzione per via numerica della equazione matriciale fondamentale, infine la terza analizza la struttura del programma, descrivendone insieme l'input e l'output e ne riporta un esempio applicativo.

     

A computational model for the limit analysis of three-dimensional masonry structures

This paper extends previous work on the limit analysis of ductile frames and plane masonry arches to the limit analysis of three-dimensional masonry structures. A lower-bound approach is developed which can handle three-dimensional collapse mechanisms involving any combination of sliding, twisting and hingeing at the block interfaces. A computer program for determining the collapse load of such structures is used to study (a) the equilibrium limits of a block with four contact points resting on an inclined plane and (b) the collapse of a semicircular arch of four blocks. The paper also describes experimental and computational work on a radially symmetric model dome of 380 blocks subject to foundation settlement.

---

Sommario Il presentre contributo estende al campo delle structture tridimensionali in muratura un precedente lavoro sull'analisi limite di telai duttili ed archi in muratura piani. Si e' sviluppato un approccio statico che analizza meccanismi di collasso tridimensionale ottenuti per combinazione dei meccanismi semplici di scorrimento e rotazione nel piano e fuori dal piano delle superfici di interfaccia tra i blocchi. Si descrivono (a) i limiti di equilibrio di un blocco con 4 punti di contatto su base inclinata, (b) le condizioni di collasso di un arco semicircolare costituito da quattro blocchi, applicando un programma di calcolo redatto per l'analisi e la definizione del carico di collasso di tali strutture. La terza parte dell'articolo presenta il lavoro sperimentale e di calcolo sviluppato su un modello di cupola a simmetria radiale costituita da 380 blocchi soggetta a cedimenti fondali.

     

A method for approximate solutions of stationary creep problems

Summary The paper presents a matrix method for the stationary creep analysis of structures.The method implies discretization of the system, piecewise linearization of the (generalized) stress-strain rate laws, some linear elastic calculations which are useful also for other purposes, and the solution of a quadratic program. The procedure adapts to various degrees of accuracy, and is highly suitable for computers.

---

Sommario Si propone un metodo matriciale per lo studio di strutture in fase di scorrimento viscoso stazionario. Il metodo comporta discretizzazione della struttura e linearizzazione a tratti delle leggi di deformabilità degli elementi, richiede calcoli elastici lineari, utili anche ad altri scopi, e la risoluzione di un programma quadratico. Il procedimento presenta notevole adattabilità a vari gradi di approssimazione e si presta particolarmente all'uso di mezzi automatici di calcolo.

---

---

The results presented in the paper form part of a series of studies supported by the National (Italian) Research Council (C.N.R.).     

Synthesis of Stephenson's six-bar linkages by an interactive technique

Summary An interactive approach to the computer aided synthesis of the planar Stephenson's six-bar linkages is reported. As known, in the applications, six-bar linkages may meet more sophisticated requirements than simpler mechanisms could allow.The design method consists in carefully selecting a suitable coupler point of a planar four-bar linkage: the coupler curve must fit the characteristics required by the law of motion of the follower. Knowledge of this curve supplies the needed directions to complete the six-bar linkage which is then automatically analysed.The analysis method and the computer program are fully explained. The program was implemented on a mini-computer: it needs essentially a video-display unit with grafic input-output capabilities.

---

Sommario E' riportata una tecnica per il progetto degli esalateri articolati piani di Stephenson, basata sull'impiego interattivo dell'elaboratore elettronico. E'noto che gli esalateri articolati, rispetto a meccanismi più semplici, consentono di soddisfare a richieste più sofisticate in certe applicazioni.Il metodo di progetto proposto consiste nella ricerca di un opportuno punto di biella di un quadrilatero articolato piano, la cui traiettoria soddisfi a precise caratteristiche imposte dal movimento del cedente. Il meccanismo finale viene costruito a partire dal quadrilatero così ottenuto e viene analizzato automaticamente.Vengono illustrati il metodo di analisi e il programma di calcolo. Quest'ultimo è stato implementato per un minicomputer e richiede essenzialmente la disponibilità di un video grafico interattivo.

     

Finite element incremental analysis of elastoplastic plate bending

Summary The paper discusses the numerical solution of the elastoplastic bending of plates with any polygonal boundary subject to a generic load and to distributed dislocations.For the discretization of the plate Herrmann's and Hellan's constant moment triangular equilibrium model has been chosen. Compatibility is imposed through a minimum principle, proposed by Capurso and Maier, which can be considered a generalization of the minimum complementary energy principle for elastoplastic workhardening continua allowing for distributed dislocations.By means of a completely automatic computation program, the cases of simply supported and clamped square plate under uniformly distributed load and under central concentrated load are solved. Present results are partly compared with existing solutions.

---

Sommario Viene studiato il problema della flessione di una piastra a contorno poligonale qualsiasi costituita da materiale elastico-incrudente e sottoposta ad un carico generico e a distorsioni diffuse. Per la discretizzazione della piastra viene scelto il modello triangolare equilibrato a momento costante di Herrmann e Hellan. La congruenza viene imposta attraverso un principio di minimo proposto da Capurso e Maier che può essere riguardato come la generalizzazione del principio della minima energia complementare per materiali elastico-incrudenti in presenza di distorsioni diffuse.Viene sviluppato un programma di calcolo automatico, che stende sulla piastra un reticolo di densità prefissata, numera e collega fra di loro nodi, lati e triangoli. In tal modo il programma genera le condizioni di equilibrio nodale e costruisce la funzione dell'energia da minimizzare. Vengono risolti, per il confronto con le soluzioni esistenti i casi di piastra quadrata semplicemente appoggiata e incastrata al contorno sottoposte a carico uniforme e ad una forza concentrata nel centro.

---

---

This study forms part of a research program supported by the National (Italian) Research Council (C.N.R. Gruppo Problemi avanzati di ingegneria strutturale).     

Complementary energy,neutral equilibrium and buckling

Summary A general analytical approach of the non-linear problem of a linear elastic and isotropic straight bar of variable stiffness is indicated. The linear equivalence method, introduced by one of the authors, is applied to two fundamental cases for the isostatic straight bar, i.e. the cantilever bar (a Cauchy type problem) and the simply supported bar (a bilocal problem). Some numerical examples concerning moderate deformations and rotations are presented.

---

Sommario Si propone un approccio analitico generate per la soluzione del problema non lineare di un'asta rettilinea di rigidezza non uniforme di materiale elastico-lineare isotropo.Il metodo dell'equivalenza lineare introdotto da uno degli autori è applicato a due casi fondamentali della trave isostatica rettilinea, cioè la mensola (problema alla Cauchy) e la trave appoggiata (problema bilocale). Vengono presentati alcuni esempi numerici concernenti deformazioni e rotazioni moderatamente grandi.

     

On the motion of a variable-mass rotor in a particular transient

Summary This paper considers a variable mass and diameter rotor having constant peripheral speed. Such rotor sketches a bobbin on which tape or wire are wound. Two methods of solution are being compared: the first one performs the integration of the movement equations following the Runge-Kutta method; the second, using the Bessel functions, reduces the integration itself to the solution of four differential equations of the first order independent from each other. The latter method offers greater numerical stability and requires less computing time. A numerical example is given.

---

Sommario Si studia un rotore a massa e diametro variabile in cui viene mantenuta costante la velocità periferica. Tale rotore schematizza un rocchetto sul quale si avvolge filo o nastro. Si confrontano due metodi di calcolo: uno che, per la integrazione delle equazioni del moto, fa uso del metodo di Runge-Kutta, l'altro che, con l'impiego delle funzioni di Bessel, riconduce l'integrazione stessa alla integrazione di quattro equazioni differenziali del primo ordine del tutto indipendenti fra loro. Il secondo metodo presenta maggiore stabilità e richiede minore tempo di calcolo. Un esempio conclude la Nota.

     

A two-step godunov-type scheme for the euler equations

A second-order Godunov-type scheme for the Euler equations in conservation form is derived. The method is based on the ENO formulation proposed by Harten et al. The fundamental difference lies in the use of a two-step scheme to compute the time evolution. The scheme is TVD in the linear scalar case, and gives oscillation-free solutions when dealing with nonlinear hyperbolic systems. The admissible time step is twice that of classical Godunovtype schemes. This feature makes it computationally cheaper than one-step schemes, while requiring the same computer storage.

---

Sommario Viene data una nuova estensione al secondo ordine del metodo di Godunov per la soluzione delle equazioni di Eulero in forma conservativa. Il metodo é basato sulla formulazione ENO proposta da Harten et al. La differenza fondamentale consiste nel calcolo dell'evoluzione temporale, ottenuta mediante uno schema a due passi. Questo consente l'uso di un passo di integrazione nel tempo doppio rispetto agli altri schemi alla Godunov ad un solo passo. Il metodo proposto risulta quindi piú efficiente e puó inoltre essere implementato senza alcun aumento dell'occupazione di memoria. Viene dimostrato che lo schema é TVD nel caso lineare, e che fornisce soluzioni prive di oscillazioni spurie nel caso di sistemi non-lineari.

     

An algorithm for the shakedown analysis of plane frames

Summary A quite straightforward procedure for the shakedown analysis of plane frames is presented. Plastic deformability is assumed as due to the bending moment only and lumped in a finite number of sections of the structure. This procedure is also suitable for the limit analysis problems of frames.

---

Sommario Si propone un procedimento speditivo per il calcolo a shakedown di strutture intelaiate piane, sotto l'ipotesi di deformabilità plastica concentrata in sezioni isolate e dipendente dal solo momento flettente. Il procedimento è anche applicabile alla risoluzione di problemi di calcolo a rottura di telai.

     

Elastic-plastic behaviour with work-hardening: An appropriate model for structural software

Summary A method is proposed for the stress analysis of structures subject to elastic and plastic deformations. The method can be adapted to exploit existing computing programs for elastic stress analysis and is sufficiently subtle to describe a variety of types of work-hardening, on the basis of a limited number of essential material parameters which can be read off usually available test graphs.

---

Sommario Si propone un metodo per l'analisi di strutture soggette a deformazioni elasto-plastiche. Il metodo può essere adattato in modo da sfruttare programmi di calcolo, di cui già si disponga, per l'analisi elastica di strutture ed è sufficientemente duttile da consentire la rappresentazione di tipi diversi di incrudimento, con l'uso di un numero limitato di parametri materiali che possono dedursi da grafici sforzo-deformazione generalmente disponibili.

     

Lyapunov Characteristic Exponents for smooth dynamical systems and for hamiltonian systems; a method for computing all of them. Part 1: Theory

Summary Since several years Lyapunov Characteristic Exponents are of interest in the study of dynamical systems in order to characterize quantitatively their stochasticity properties, related essentially to the exponential divergence of nearby orbits. One has thus the problem of the explicit computation of such exponents, which has been solved only for the maximal of them. Here we give a method for computing all of them, based on the computation of the exponents of order greater than one, which are related to the increase of volumes. To this end a theorem is given relating the exponents of order one to those of greater order. The numerical method and some applications will be given in a forthcoming paper.

---

Sommario Da diversi anni gli esponenti caratteristici di Lyapunov sono divenuti di notevole interesse nello studio dei sistemi dinamici al fine di caratterizzare quantitativamente le proprietà di stocasticità, legate essenzialmente alla divergenza esponenziale di orbite vicine. Si presenta dunque il problema del calcolo esplicito di tali esponenti, già risolto solo per il massimo di essi. Nel presente lavoro si dà un metodo per il calcolo di tutti tali esponenti, basato sul calcolo degli esponenti di ordine maggiore di uno, legati alla crescita di volumi. A tal fine si dà un teorema che mette in relazione gli esponenti di ordine uno con quelli di ordine superiore. Il metodo numerico e alcune applicazioni saranno date in un sucessivo articolo.

     

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.

     

Universal solutions and energy minimization in the saint-venant cylinder

In this paper we study the role played, from an energy view-point, by the explicit solution of Saint-Venant for a right circular cylinder. By taking into account a harmonic solution with constant divergence, we consider Maisonneuve's solution and investigate the assumptions under which the explicit solutions form a particular case.

---

Sommario In questo lavoro si esamina, da un punto di vista energetico, il ruolo che le soluzioni esplicite di Saint-Venant hanno all'interno dell'intera classe di soluzioni equipollenti per un cilindro circolare retto. A tale scopo facendo riferimento alla soluzione ottenuta da Maisonneuve, che minimizza l'energia di deformazione nell'intera classe delle soluzioni equipollenti, si vuole studiare se esiste una intersezione non nulla tra la classe delle soluzioni di Maisonneuve e la classe delle soluzioni armoniche a divergenza costante.

     

A quadratic programming approach for certain classes of non linear structural problems

Summary The search for the deviation from the linear elastic, of the response of a nonlinear discrete structure to given loads, is formulated for various classes of cases. All the problems considered are shown to be reducible to a single mathematical model: the quadratic program with sign constraints only. The Kubn-Tucker theorem for quadratic programming directly supplies the appropriate extremum principles, for which some physical interpretations are proposed. Other results of quadratic programming theory prove to be useful in the discussion of the structural problems considered.

---

Sommario Si considerano certe categorie di strutture costituite da un numero discreto di parti a comportamento non lineare e si formula il problema della ricerca dello scostamento del loro stato di sforzo e deformazione per dati carichi, da quello che si avrebbe in regime elastico lineare. Nei vari casi considerati il problema si riduce ad un programma quadratico con variabili soggette solo a restrizioni di segno. Il teorema di Kuhn e Tucker per programmazione quadratica fornisce gli opportuni principi di estremo, per i quali si propongono interpretazioni fisiche. Altri teoremi di programmazione quadratica si rivelano utili nel discutere la soluzione dei problemi formulati.

---

---

This research was conducted during the author's stay at the Engineering Department of Cambridge University, England, with the financial support of a N.A.T.O. fellowship 1965–66. The author wishes to thank Professor Sir John Baker, Head of the Department, for the kind hospitality, and the C.N.R. (National Research Council) for the award of the fellowship.     

Multibody dynamics. A rapidly developing field of applied mechanics

Summary Using the principle of virtual power nonlinear equations of motion are developed for arbitrary systems of rigid bodies interconnected by kinematical joints, springs, dampers and active force elements. An essential system parameter is the so-called path matrix which describes the system topology. A computer program generates the equations in symbolical form for subsequent numerical simulations of dynamical system behaviour. The program finds applications in various branches of industry as well as in scientific research.

---

Sommario Usando il principio delle potenze virtuali si sviluppano equazioni nonlineari di moto per sistemi arbitrari di corpi rigidi interconnessi con giunti cinematici, molle, ammortizzatori ed elementi che esercitano forze attive. Un parametro essenziale del sistema è la cosiddetta matrice del cammino che descrive la topologia del sistema. Un programma di calcolo genera le equazioni in forma simbolica per successive simulazioni numeriche del comportamento del sistema dinamico. Il programma trova applicazioni sia in vari settori dell'industria sia nella ricerca scientifica.

---

---

Paper presented at the 8th National Congress of AIMETA Turin, September 29–October 3, 1986 (invited lecture).     

Finite element equilibrium analysis of creep using the mean value of the equivalent shear modulus

Summary The paper deals with the nonlinear creep analysis using equilibrium finite elements. In adopting a concept of the mean value of the equivalent shear modulus, a direct iterative procedure is proposed with the minimum calculation in the modification of the global stiffness matrix. The difficulties due to the incompressibility of the creep flow are overcome by a numerical artifice where the hydrostatic stress component is considered as a Lagrange parameter. Some numerical results are shown and compared with the analytical method.

---

Sommario La memoria tratta l'analisi dello scorrimento viscoso non lineare usando l'equilibrio di elementi finiti. Adottando il concetto di valore medio di un modulo di taglio equivalente, si propone un procedimento iterattivo diretto con il minimo di calcolo nella modificazione della matrice di rigidezza globale. Le difficoltà dovute alla incomprimibilità dello scorrimento viscoso sono superate con un artificio numerico in cui la componente di sollecitazione idro-statica è considerata come un parametro di Lagrange. Vengono mostrati alcuni risultati numerici, e confrontati con il metodo analitico.

     

Propagation of periodic long waves

Summary A method is developed for the computation of the steady solution of the shallow water equations with quasi periodic boundary conditions. Because of dissipation the influence of the initial conditions becomes negligible with increasing time and the solution finally depends on the boundary conditions. The unknowns variables (velocity and surface elevation) and the boundary conditions are developed in power series of a small perturbation parameter. The problem is then transformed in a sequence of linear problems which have the same associated homogeneous problem. By separating the time and space variables in the homogeneous problem we obtain an homogeneous elliptic problem of which we compute the first eigenvalues and eigensolutions. These are related to the characteristic oscillations of the water in the basin. The solution of each linear problem is then obtained as an eigensolution expansion with time dependent coefficients. These coefficients are solutions of ordinary differential equations which can be solved directly without proceeding step by step in time. In this way we are reduced to a stationary problem i.e. the determination of the eigenvalues and eigensolutions of the elliptic problem and to the computation of several integrals needed for the determination of the time dependent coefficients. A first test of the method has been carried out for a one-dimensional problem i.e. the tidal wave in a canal of finite length and constant depth. In this case the various steps of the procedure outlined above can be performed analitically. The results have been compared with those obtained by a step by step numerical integration of the shallow water equations. The agreement between these sets of results is good for the range of values of the parameters currently used in the applications.

---

Sommario Viene presentato un metodo per il calcolo della soluzione a regime delle equazioni delle onde lunghe, dipendente dalle sole condizioni al contorno quasi-periodiche, dopo aver mostrato che l'influenza delle condizioni iniziali diventa col tempo trascurabile a causa del termine di resistenza. Il metodo si basa sullo sviluppo in serie di potenze di un piccolo parametro sia delle incognite (velocità ed elevazioni del pelo libero) sia delle condizioni al contorno al fine di trasformare il problema non lineare in una serie di problemi lineari aventi lo stesso problema omogeneo associato. Con la separazione delle variabili spazio, tempo quest'ultimo problema viene ricondotto ad un problema ellittico omogeneo di cui si calcolano i primi autovalori ed autosoluzioni. Da ultimo la soluzione di ciascun problema lineare è ottenuta come sviluppo in serie di autosoluzioni a coefficienti dipendenti dal tempo: questi si ricavano risolvendo analiticamente delle equazioni differenziali ordinarie. Si elimina cosi la necessità di procedere passo passo nel tempo analogamente ai classici metodi armonici di soluzione di sistemi lineari. Riassumendo, l'applicazione del metodo riconduce alla soluzione di un problema stazionario (determinazione di autovalori ed autosoluzione del problema ellittico) e quindi al calcolo dei vari integrali necessari per la determinazione dei coefficienti temporali. Il metodo è stato provato nel caso semplice della propagazione di onde di marea in un canale di lunghezza finita e sezione costante. Per questo esempio i vari passi di calcolo possono essere svolti analiticamente. I risultati sono stati confrontati con quelli ottenuti dalla integrazione numerica delle equazioni col metodo delle caratteristiche, ottenendo un buon accordo.

---

---

Work sponsored by the CNR (National Research Council) in the framework of Project ?Conservazione del suolo? Subproject ?Dinamica dei Litorali? Publ. n° 46.     

On plastic dynamics of discrete structural models

Summary This paper discusses the dynamic loading of rigid-perfectly plastic structures by adopting a structural model discretized with constant stress finite elements. This assumption together with the hypothesis of small displacements and of a piecewise linear yield surface leads to the formulation of a problem in linear inequalities, which is, implicitly, time discretized. It is shown how the non linearity of the associated plastic flow laws leads to an equivalent extremal formulation. A numerical algorithm for solving the problem is directly derived from the mechanical statements. Numerical examples illustrate this approach. In the appendix an approximate technique is developed, which takes into account the influence of the strain-hardening and the strain rate sensitivity of material.

---

Sommario Il problema del caricamento dinamico delle strutture rigido-plastiche viene affrontato attraverso un modello discretizzato per elementi finiti a stato di sollecitazione costante. Tale assunzione insieme alle ipotesi di superficie di plasticizzazione lineare a tratti e di piccoli spostamenti, conduce alla formulazione di un problema retto da diseguaglianze lineari che si presenta discretizzato, implicitamente, nel tempo. Si mostra che la non linearità delle leggi del flusso plastico associato conduce ad una equivalente formulazione di estremo.È esposto un algoritmo risolvente affatto aderente alla formulazione meccanica. Alcuni esempi illustrano il procedimento. Una proposta di adattamento al caso di materiali incrudenti e sensibili alle velocità di deformazione è presentata in appendice.

---

---

Research supported by the National Research Council (C.N.R. - P.A.d.I.S.) of Italy.     

Finite element analysis of axisymmetric systems without ignoring shear strain effects

Summary The finite element method is used to calculate the elastic displacements and internal forces in shells with axisymmetric geometry and load; the shell is assumed to be strainable by shear as well as by bending and normal forces. It is assumed that the finite elements have a simple curvature whose displacement function is analytically expressed by power polynomials of the meridian coordinate; the chosen parameters are the displacement components of the nodes; their number may be smaller than that of the polynomial terms, whose dependence is therefore uniquely defined by the minimum total potential energy condition. By this means it is possible to obtain an accurate numeric solution even with a few elements as long as they are sufficient to represent the shell geometry and follow the load distribution. The exposition is completed by the block diagram of the numeric computing program and the results in some examples.

---

Sommario Col metodo degli elementi finiti viene impostato il calcolo degli spostamenti elastici e delle azioni interne in lastre assialsimmetriche per geometria e per carico; si suppone la lastra deformabile oltre che alle azioni flettenti e normali anche all'azione di taglio.Si assumono elementi finiti a semplice curvatura con funzione spostamento analiticamente espressa da polinomi di potenze della coordinata meridiana; si scelgono quali parametri le componenti di spostamento dei nodi; il loro numero può essere inferiore a quello dei coefficienti polinomiali considerati la cui dipendenza viene pertanto univocamente definita con la condizione di minima energia potenziale totale.Seguendo questa via è possibile ottenere una soluzione numerica accurata anche con pochi elementi finiti purchè sufficienti a rappresentare convenientemente la geometria della lastra e a seguire la distribuzione dei carichi.L'esposizione è completata dallo schema a blocchi del programma di calcolo numerico per l'elaborazione elettronica e dai risultati ottenuti in alcuni casi esemplificativi.

     

Lifting Surface in Subsonic Unsteady Regime

The paper refers to the solution of the integral equation for the acceleration (or pressure) potential for the study of subsonic linearized unsteady flow in view of aeroelastic applications. The case considered is relevant to a trapezoidal wing infinitely thin surface without discontinuities. As is well known [1, 2], the kernel of the integral equation exhibits three singularities, two of which are integrable in elementary form, whereas, for the third one integration in principal part according to Hadamard's rule is necessary. The kernel is therefore reworked in such a way that all the singularities are separated from the regular part, and eventually the discretization is performed in such a way that only the regular part is to be recalculated for each new value of the reduced frequency. Convergence tests, comparison with other methods of solution, and time saving associated with the technique of separation are also shown.Sommario. II lavoro tratta la risoluzione del problema relativo alla equazione integrale nel potenziale di accelerazione (o di pressione) per lo studio di una corrente subsonica linearizzata nonstazionaria, in vista di applicazioni aeroelastiche. Il caso considerato è quello di una superficie alare a pianta trapezoidale in assenza di discontinuità di spessore infinitesimo. Come è noto [1, 2], il nucleo della equazione integrale in parola presenta tre singolarità, due sole delle quali sono integrabili in forma elementare, (o riconducibili ad essa), mentre per la terza è necessario far ricorso alla integrazione in parte principale alla Hadamard. Il nucleo stesso viene quindi rielaborato in modo da isolare tutte le singolarità dalla componente regolare del nucleo; si procede così alla discretizzazione dell' equazione integrale, e, per ogni valore della frequenza ridotta, va ricalcolata solo la parte regolare della matrice risolvente. Vengono poi effettuati tests di convergenza, confronti con altri metodi di soluzione, analisi sui tempi di calcolo e risparmio di tempo di calcolo dovuto alla tecnica di separazione.