THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS

Size: px
Start display at page:

Download "THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS"

Transcription

1 ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June 04 THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS M. Wysocki a,b*, M. Szpieg a, P. Heström a and F. Ohsson c a Swerea SICOMP AB, Box 04, SE-43 Mönda, Sweden b Chamers University of Technoogy, Dept. of Appied Mechanics, SE-4 96 Göteborg, Sweden c Oxeon AB, Företagsgatan 4, SE Borås, Sweden *maciej.wysocki@swerea.se Keywords: Preform mechanics, textie composites, Hertzian contact, experimenta characterisation. Abstract In this paper the constitutive compressive behaviour of neary parae spread-tow textie reinforcement is studied. The striking resut of our anaysis is that the spread-tow type of reinforcement shoud obey inear reation between force and deformation. This is in contrast to standard textie reinforcements that obey a power-aw type of behaviour. To support the theoretica investigation we have deveoped an test rig who's chief purpose is to achieve compression between neary perfecty parae surfaces. This is achieved using a mechanica arrangement consisting of a ba-joint.. Introduction Presenty, the aerospace industry is ooking into spread-tow textie reinforcement for use in primary structura composites. This is mainy due to their reativey ow price and improved mechanica performance, hand abiity and manufacturabiity. As introduction of new structura composite materias into aero structures requires a timey and costy certification process, there is presenty a arge effort towards characterisation, understanding and modeing various aspects of the spread-tow based composite materias. Since forming of reinforcements pays a key roe in terms of the subsequent composite manufacturing and the performance of the fina product, we focus in this contribution on the eastic compression of such a preform. Since, the mechanica properties of a fibre mass are important in many fieds of engineering, incuding composite manufacturing, a significant effort has been spent to describe and mode these properties. In particuar, Van Wyk [] pioneered the mechanistic anaysis of the compressibiity of 3D random fibre masses. Van Wyk regarded the fibre mass as a system of bending units consisting of fibre beam eements between adjacent fibre contacts and he ignored twisting, sip, and extension of the fibres. His key assumptions were that the mean contact spacing is proportiona to the reciproca of the fibre voume fraction and that the segments deform as in bending of straight sender beams. The resut of his anaysis is a simpe power aw for the compressive stress as P k E, () 3 3 0

2 ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June 04 where k is a structure dependent constant; E is the Young s moduus of the fibres; and 0 is the imiting fibre voume fraction beow which P=0. Most of the work foowing van Wyk has accepted the form Eq. and focused on finding an appropriate expression for k by extending the theories describing the structure deveopment during deformation. It was not unti 998 that To [] suggested a more genera expression for the packing of fibre masses n n P k E 0, () where the exponent, n, was shown anayticay to take the vaue of 3 for the 3D and 5 for the D random cases. The author aso succeeded to fit experimenta data for materias where the contacts between fibres are ines rather than points, e.g., fibre bundes, by adjusting the exponent in Eq.. Nevertheess, since the assumptions used in the anayses in [, ] proceeds from an assumption of point contact and bending of fibre segment, the fitting of Eq. to data for highy parae fibre architectures in [] is to be considered as painy phenomenoogica. Figure. Geometry of the contact between eiptic parobooids with principa axes of body (x, y, z ) and the principa axes of body (x, y, z ). The ange is the ange between the x, z and x z panes, i.e. the principa radii of curvatures R and R. A striking quaity of the spread tow technoogy is the possibiity to obtain extremey parae fibre assembies. This is in contrast to a typica reinforcement where the fibre architecture are usuay somehow irreguar. Therefore in the context of mechanica behaviour, the compressive response of the highy oriented spread-tow fabrics is expected to be significanty different from the standard reinforcement. The main difference between the two types of fabrics is the constitutive assumption of oad transfer mechanism between the fibres []. Namey, in standard fabric the oad transfer is assumed via bending of fibre segments, whie in spread tow fabric the main mechanism is assumed to be Hertzian contact between adjacent fibres. Therefore, in this paper we study the constitutive behaviour of highy parae fibre network governed by Hertzian contact. The theoretica study is aso supported by experimenta characterisation of the transverse compressibiity of the highy anisotropic and packed CF spread-tow fibre assembies. The resuting force-defection curve is anaysed against the Hertzian contact modeing effort to give effective transverse moduus.

3 ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June 04. Theoretica background The probem of contact between eastic bodies has ong been of considerabe interest. Assume that two eastic soids are brought into contact at a point 0, as shown in Fig.. If coinear forces are appied so as to press the two soids together, deformation occurs, and we expect a sma contact area to repace the point of the unoaded state. If we woud abe to determine the size and shape of this contact area and the distribution of norma pressure, then the interva stresses and deformation coud be cacuated. Figure. Geometry of deformed bodies. Dashed ines show the surface as they woud be in the absence of deformation. Continuous ines show the surfaces of the deformed bodies. In the anaysis of the contact probem by Hertz, as extended by Timoshenko and Goodier [3], the foowing is assumed: () the contacting surfaces are perfecty smooth so that the actua shape can be described by a second degree equation of the form z Dx Ey Fxy where D, E and F are arbitrary constants.; the eastic imit of the materias are not exceeded during the contact; (3) ony norma forces between the contacting surfaces are considered; and (4) the contacting surfaces are sma in comparison to the entire surfaces. Based on the above assumptions and by appying potentia theory it can be showed that: () the contact area is bounded by an eipse whose semi-axes can be cacuated from the geometric parameters of the contacting bodies; and () the norma pressure distribution over this area is x a x P p0 b, where p 0 is maximum pressure at centre, a is the major axis of eipse of contact and b is the minor axis of eipse. If the two bodies are pressed together by appied norma forces (cf. Figure 3), then a deformation occurs near the origina point of contact aong the Z-axis. Here again, we consider ony forces acting parae to the z-axis where the distance from the z-axis is sma. The dispacements at a point are w and w where w is the deformation of point P of body and w is the deformation of point P for body, pane C is the origina pane of tangency; z is the distance from P to the undeformed state, and z is the distance from P to the undeformed state. For points inside the contact area, we have z w z, where P w is the deformation function we seek. To sove the probem, consider the contact surfaces and define the foowing geometrica reations [4] A B R R R R (3) 3

4 ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June 04 and B A cos. (4) R R R R R R R R The constants A and B are expressions for the combinations of the principa curvatures of the surfaces and the ange between the panes of curvature. These are used to cacuate via the reation Ax By w w. (5) The probem is now to find a pressures distribution and potentia functions to satisfy the Eq. 5. The generic soution based on the use of eiptic integras is to found in numerous textbooks onto contact mechanics, see e.g. [4]. Therefore we wi here focus on the contact mechanics of fibres in contact with fibres and pane, i.e. the specia case of ine contact. To sove this probem we sha make use of the expressions aready deveoped [3] to sove for the "pressure distribution" and size of the area of contact by aowing one axis of the eipse of contact to become infinite. To determine the deformation, the contact area wi be taken as being a finite rectange with one side very much greater than the other. The derivation wi be for the case of a pair of cyinders with their axes parae and is based on the work presented in [3]. The soution for a cyinder to pane contact can easiy be obtained by aowing the radius of one of the cyinders to become infinite. Figure 3. Contact geometry between two parae cyinders. Line contact occurs when two cyinders rest on each other with their axes parae, cf. Fig. 3, and when a cyinder rests on a pane. As the two cyinders are pressed together aong their axes, the resuting pressure area is a narrow "rectange " of width b and ength L (assuming no taper in the cyinders). In other words, the area of contact is an eongated eipse with the major axis of the eipse equa to L and the eccentricity approaching unity. Without going into detais, due to imited space, the tota deformation of a pair of cyinders with their axes parae (after some engthy agebra) can be obtained as P n 4 n L n b, (6) L 4

5 ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June 04 where b is given by b 4 P RR L R R (7) and where i i. (8) E i The Eq. 6 (together with Eq. 7) is simpified to the foowing expression for a cyinder on a fat surface 3 P L n n 4. (9) L 4 PR Machine attachment Ba joint Testing area Machine attachment Figure 4. Sef-adjusting compaction rig. The main resut of the above derivations is the observation that the reation between the deformation and the appied force is neary inear for moderate oadings such as the expected for typica manufacturing processes of composite materia. This is evident when representative materia data for carbon fibres and stee pate are used in Eqs. 6 and 9: CFs 9 transverse stiffness of E 50 Pa and Poisons ratio of 0. 35, CFs diameter of T T 6 R fib 30 m, together with data for stee. This eads to a reation between deformation and force in terms of 7 8 P n P or P (0a) for contact between two fibres, and as 5

6 ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June P n P or P (0b) for contact between a fibre and stee surface. Of course, the above derivations are for a singe fibre ony and needs to be further revised for mutipe fibres and fibre assembies. However, one can sti concude that compressibiity of parae fibre assembies is (neary) inear with respect to the appied force (foowing from a simpe assumption of P in Eq. 0a,b being a sum of individua fibre forces). This is in contrast to the commony observed behaviour in standard reinforcement where the behaviour is expected to be power aw, cf. Eqs. and.. Experimenta method and resuts Since compressive testing of extremey thin fabric (in the case of the studied TeXtreme UD spread-tow band approximatey 0.05 mm) requires an terrific paraeism of the testing equipment, we decided to deveop an sef-adjusting test rig where the adjustment is performed using a ba-joint, cf. Fig. 4. The test rig was machined from a singe piece of stee and consist of machine attachment, testing area and a ba-joint used for the paraeism adjustment. Machine compiance was subtracted from the experimenta force-defection data by assuming a spring mode. The tested materia was TeXtreme UD spread-tow bands from Oxeon consisting of TR50S 5k CF fibres. A typica force-defection data for the TeXtreme UD spread-tow, after correcting for machine compiance, resuts in an amost inear reation between Force and Defection given as P Concuding remarks In this contribution we have studied the constitutive behaviour of neary parae fibre assembies in the form of TeXtreme UD spread-tow bands from Oxeon. The resuts indicate that the compressive behaviour of the neary parae fibres shoud obey a inear reation, see Eqs. 0a and b, which is in contrast to the standard textie reinforcements that obey power aw type of behaviour. This theoretica investigation is further supported by our initia experimenta data which shows amost inear behaviour. Compared to the theoretica investigation, the experimenta data shows however somehow higher compiance. This coud be due to various phenomena such as, amongst other; imited fibre paraeism resuting in partia or eiptica-type of contact, non-homogeneous fibre distribution and theoretica deveopment performed for a singe fibre ony. This needs to be further investigated into in the near future. References [] C.M. van Wyk, Note on the Compressibiity of Woo, Journa of the Textie Institute, 85-9, 946 [] S. To, Packing mechanics of fiber reinforcements, Poymer Engineering and Science, , 998 [3] S.P. Timoshenko and. L.N. Goodier, Axisymmetric stress and deformation in a soid of revoution, in: Theory of Easticity, 3rd ed., pp , McGraw Hi, NY, 95 [4] K.L. Johnson, Contact Mechanics, Oxford University Press, Oxford, 985 6

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s

More information

3.10 Implications of Redundancy

3.10 Implications of Redundancy 118 IB Structures 2008-9 3.10 Impications of Redundancy An important aspect of redundant structures is that it is possibe to have interna forces within the structure, with no externa oading being appied.

More information

Strain Energy in Linear Elastic Solids

Strain Energy in Linear Elastic Solids Strain Energ in Linear Eastic Soids CEE L. Uncertaint, Design, and Optimiation Department of Civi and Environmenta Engineering Duke Universit Henri P. Gavin Spring, 5 Consider a force, F i, appied gradua

More information

Separation of Variables and a Spherical Shell with Surface Charge

Separation of Variables and a Spherical Shell with Surface Charge Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation

More information

In-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017

In-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017 In-pane shear stiffness of bare stee deck through she finite eement modes G. Bian, B.W. Schafer June 7 COLD-FORMED STEEL RESEARCH CONSORTIUM REPORT SERIES CFSRC R-7- SDII Stee Diaphragm Innovation Initiative

More information

Lecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling

Lecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling Lecture 9 Stabiity of Eastic Structures Lecture 1 Advanced Topic in Coumn Bucking robem 9-1: A camped-free coumn is oaded at its tip by a oad. The issue here is to find the itica bucking oad. a) Suggest

More information

UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE

UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE Juan Huang, Ronghui Wang and Tao Tang Coege of Traffic and Communications, South China University of Technoogy, Guangzhou, Guangdong 51641,

More information

MECHANICAL ENGINEERING

MECHANICAL ENGINEERING 1 SSC-JE SFF SELECION COMMISSION MECHNICL ENGINEERING SUDY MERIL Cassroom Posta Correspondence est-series16 Rights Reserved www.sscje.com C O N E N 1. SIMPLE SRESSES ND SRINS 3-3. PRINCIPL SRESS ND SRIN

More information

Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg

Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Deft University of Technoogy Marijn Drienburg October 2017 Contents 1 Introduction 2 1.1 Hand Cacuation....................................

More information

Nonlinear Analysis of Spatial Trusses

Nonlinear Analysis of Spatial Trusses Noninear Anaysis of Spatia Trusses João Barrigó October 14 Abstract The present work addresses the noninear behavior of space trusses A formuation for geometrica noninear anaysis is presented, which incudes

More information

APPENDIX C FLEXING OF LENGTH BARS

APPENDIX C FLEXING OF LENGTH BARS Fexing of ength bars 83 APPENDIX C FLEXING OF LENGTH BARS C.1 FLEXING OF A LENGTH BAR DUE TO ITS OWN WEIGHT Any object ying in a horizonta pane wi sag under its own weight uness it is infinitey stiff or

More information

SIMULATION OF TEXTILE COMPOSITE REINFORCEMENT USING ROTATION FREE SHELL FINITE ELEMENT

SIMULATION OF TEXTILE COMPOSITE REINFORCEMENT USING ROTATION FREE SHELL FINITE ELEMENT 8 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS SIMULATION OF TEXTILE COMPOSITE REINFORCEMENT USING ROTATION FREE SHELL FINITE ELEMENT P. Wang, N. Hamia *, P. Boisse Universite de Lyon, INSA-Lyon,

More information

On a geometrical approach in contact mechanics

On a geometrical approach in contact mechanics Institut für Mechanik On a geometrica approach in contact mechanics Aexander Konyukhov, Kar Schweizerhof Universität Karsruhe, Institut für Mechanik Institut für Mechanik Kaiserstr. 12, Geb. 20.30 76128

More information

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed

More information

DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE

DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE 3 th Word Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 4 Paper No. 38 DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE Bo JIN SUMMARY The dynamic responses

More information

c 2007 Society for Industrial and Applied Mathematics

c 2007 Society for Industrial and Applied Mathematics SIAM REVIEW Vo. 49,No. 1,pp. 111 1 c 7 Society for Industria and Appied Mathematics Domino Waves C. J. Efthimiou M. D. Johnson Abstract. Motivated by a proposa of Daykin [Probem 71-19*, SIAM Rev., 13 (1971),

More information

2.1. Cantilever The Hooke's law

2.1. Cantilever The Hooke's law .1. Cantiever.1.1 The Hooke's aw The cantiever is the most common sensor of the force interaction in atomic force microscopy. The atomic force microscope acquires any information about a surface because

More information

ELASTICITY PREVIOUS EAMCET QUESTIONS ENGINEERING

ELASTICITY PREVIOUS EAMCET QUESTIONS ENGINEERING ELASTICITY PREVIOUS EAMCET QUESTIONS ENGINEERING. If the ratio of engths, radii and young s modui of stee and brass wires shown in the figure are a, b and c respectivey, the ratio between the increase

More information

Section 6: Magnetostatics

Section 6: Magnetostatics agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. The

More information

Lecture 6: Moderately Large Deflection Theory of Beams

Lecture 6: Moderately Large Deflection Theory of Beams Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey

More information

Technical Data for Profiles. Groove position, external dimensions and modular dimensions

Technical Data for Profiles. Groove position, external dimensions and modular dimensions Technica Data for Profies Extruded Profie Symbo A Mg Si 0.5 F 25 Materia number.206.72 Status: artificiay aged Mechanica vaues (appy ony in pressing direction) Tensie strength Rm min. 245 N/mm 2 Yied point

More information

1D Heat Propagation Problems

1D Heat Propagation Problems Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2

More information

VTU-NPTEL-NMEICT Project

VTU-NPTEL-NMEICT Project MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 VTU-NPTE-NMEICT Project Progress Report The Project on Deveopment of Remaining Three Quadrants to NPTE Phase-I under grant in aid

More information

Problem Set 6: Solutions

Problem Set 6: Solutions University of Aabama Department of Physics and Astronomy PH 102 / LeCair Summer II 2010 Probem Set 6: Soutions 1. A conducting rectanguar oop of mass M, resistance R, and dimensions w by fas from rest

More information

Module 22: Simple Harmonic Oscillation and Torque

Module 22: Simple Harmonic Oscillation and Torque Modue : Simpe Harmonic Osciation and Torque.1 Introduction We have aready used Newton s Second Law or Conservation of Energy to anayze systems ike the boc-spring system that osciate. We sha now use torque

More information

International Journal of Advance Engineering and Research Development

International Journal of Advance Engineering and Research Development Scientific Journa of Impact Factor (SJIF): 4.4 Internationa Journa of Advance Engineering and Research Deveopment Voume 3, Issue 3, March -206 e-issn (O): 2348-4470 p-issn (P): 2348-6406 Study and comparison

More information

СРАВНИТЕЛЕН АНАЛИЗ НА МОДЕЛИ НА ГРЕДИ НА ЕЛАСТИЧНА ОСНОВА COMPARATIVE ANALYSIS OF ELASTIC FOUNDATION MODELS FOR BEAMS

СРАВНИТЕЛЕН АНАЛИЗ НА МОДЕЛИ НА ГРЕДИ НА ЕЛАСТИЧНА ОСНОВА COMPARATIVE ANALYSIS OF ELASTIC FOUNDATION MODELS FOR BEAMS СРАВНИТЕЛЕН АНАЛИЗ НА МОДЕЛИ НА ГРЕДИ НА ЕЛАСТИЧНА ОСНОВА Милко Стоянов Милошев 1, Константин Савков Казаков 2 Висше Строително Училище Л. Каравелов - София COMPARATIVE ANALYSIS OF ELASTIC FOUNDATION MODELS

More information

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization

More information

Easticity. The strain produced in the stretched spring is ) Voume Strain ) Shearing Strain 3) Tensie Strain 4) None of the above. A body subjected to strain a number of times does not obey Hooke's aw due

More information

Lobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z

Lobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z Chapter W3 Mechanica Systems II Introduction This companion website chapter anayzes the foowing topics in connection to the printed-book Chapter 3: Lumped-parameter inertia fractions of basic compiant

More information

Post-buckling behaviour of a slender beam in a circular tube, under axial load

Post-buckling behaviour of a slender beam in a circular tube, under axial load Computationa Metho and Experimenta Measurements XIII 547 Post-bucking behaviour of a sender beam in a circuar tube, under axia oad M. Gh. Munteanu & A. Barraco Transivania University of Brasov, Romania

More information

Physics 235 Chapter 8. Chapter 8 Central-Force Motion

Physics 235 Chapter 8. Chapter 8 Central-Force Motion Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion

More information

FOURIER SERIES ON ANY INTERVAL

FOURIER SERIES ON ANY INTERVAL FOURIER SERIES ON ANY INTERVAL Overview We have spent considerabe time earning how to compute Fourier series for functions that have a period of 2p on the interva (-p,p). We have aso seen how Fourier series

More information

Transverse Anisotropy in Softwoods

Transverse Anisotropy in Softwoods Transverse Anisotropy in Softwoods Modeing and Experiments CARL MODÉN Licenciate Thesis Stockhom, Sweden 2006 TRITA-AVE 2006:30 ISSN 1651-7660 ISBN 91-7178-385-7 KTH Engineering Sciences Department of

More information

$, (2.1) n="# #. (2.2)

$, (2.1) n=# #. (2.2) Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

Volume 13, MAIN ARTICLES

Volume 13, MAIN ARTICLES Voume 13, 2009 1 MAIN ARTICLES THE BASIC BVPs OF THE THEORY OF ELASTIC BINARY MIXTURES FOR A HALF-PLANE WITH CURVILINEAR CUTS Bitsadze L. I. Vekua Institute of Appied Mathematics of Iv. Javakhishvii Tbiisi

More information

Jackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

Jackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell Jackson 4.10 Homework Probem Soution Dr. Christopher S. Baird University of Massachusetts Lowe PROBLEM: Two concentric conducting spheres of inner and outer radii a and b, respectivey, carry charges ±.

More information

Problem set 6 The Perron Frobenius theorem.

Problem set 6 The Perron Frobenius theorem. Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator

More information

PhysicsAndMathsTutor.com

PhysicsAndMathsTutor.com . Two points A and B ie on a smooth horizonta tabe with AB = a. One end of a ight eastic spring, of natura ength a and moduus of easticity mg, is attached to A. The other end of the spring is attached

More information

FRIEZE GROUPS IN R 2

FRIEZE GROUPS IN R 2 FRIEZE GROUPS IN R 2 MAXWELL STOLARSKI Abstract. Focusing on the Eucidean pane under the Pythagorean Metric, our goa is to cassify the frieze groups, discrete subgroups of the set of isometries of the

More information

ANALYTICAL AND EXPERIMENTAL STUDY OF FRP-STRENGTHENED RC BEAM-COLUMN JOINTS. Abstract

ANALYTICAL AND EXPERIMENTAL STUDY OF FRP-STRENGTHENED RC BEAM-COLUMN JOINTS. Abstract ANALYTICAL AND EXPERIMENTAL STUDY OF FRP-STRENGTHENED RC BEAM-COLUMN JOINTS Dr. Costas P. Antonopouos, University of Patras, Greece Assoc. Prof. Thanasis C. Triantafiou, University of Patras, Greece Abstract

More information

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische

More information

4 1-D Boundary Value Problems Heat Equation

4 1-D Boundary Value Problems Heat Equation 4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x

More information

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is

More information

Multiple Beam Interference

Multiple Beam Interference MutipeBeamInterference.nb James C. Wyant 1 Mutipe Beam Interference 1. Airy's Formua We wi first derive Airy's formua for the case of no absorption. ü 1.1 Basic refectance and transmittance Refected ight

More information

Proceedings of Meetings on Acoustics

Proceedings of Meetings on Acoustics Proceedings of Meetings on Acoustics Voume 9, 23 http://acousticasociety.org/ ICA 23 Montrea Montrea, Canada 2-7 June 23 Architectura Acoustics Session 4pAAa: Room Acoustics Computer Simuation II 4pAAa9.

More information

CE601-Structura Anaysis I UNIT-IV SOPE-DEFECTION METHOD 1. What are the assumptions made in sope-defection method? (i) Between each pair of the supports the beam section is constant. (ii) The joint in

More information

4 Separation of Variables

4 Separation of Variables 4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE

More information

Related Topics Maxwell s equations, electrical eddy field, magnetic field of coils, coil, magnetic flux, induced voltage

Related Topics Maxwell s equations, electrical eddy field, magnetic field of coils, coil, magnetic flux, induced voltage Magnetic induction TEP Reated Topics Maxwe s equations, eectrica eddy fied, magnetic fied of cois, coi, magnetic fux, induced votage Principe A magnetic fied of variabe frequency and varying strength is

More information

Unit 48: Structural Behaviour and Detailing for Construction. Deflection of Beams

Unit 48: Structural Behaviour and Detailing for Construction. Deflection of Beams Unit 48: Structura Behaviour and Detaiing for Construction 4.1 Introduction Defection of Beams This topic investigates the deformation of beams as the direct effect of that bending tendency, which affects

More information

Laboratory Exercise 1: Pendulum Acceleration Measurement and Prediction Laboratory Handout AME 20213: Fundamentals of Measurements and Data Analysis

Laboratory Exercise 1: Pendulum Acceleration Measurement and Prediction Laboratory Handout AME 20213: Fundamentals of Measurements and Data Analysis Laboratory Exercise 1: Penduum Acceeration Measurement and Prediction Laboratory Handout AME 20213: Fundamentas of Measurements and Data Anaysis Prepared by: Danie Van Ness Date exercises to be performed:

More information

(Refer Slide Time: 2:34) L C V

(Refer Slide Time: 2:34) L C V Microwave Integrated Circuits Professor Jayanta Mukherjee Department of Eectrica Engineering Indian Intitute of Technoogy Bombay Modue 1 Lecture No 2 Refection Coefficient, SWR, Smith Chart. Heo wecome

More information

RELUCTANCE The resistance of a material to the flow of charge (current) is determined for electric circuits by the equation

RELUCTANCE The resistance of a material to the flow of charge (current) is determined for electric circuits by the equation INTRODUCTION Magnetism pays an integra part in amost every eectrica device used today in industry, research, or the home. Generators, motors, transformers, circuit breakers, teevisions, computers, tape

More information

Analysis of Cylindrical Tanks with Flat Bases by Moment Distribution Methods

Analysis of Cylindrical Tanks with Flat Bases by Moment Distribution Methods May, 958 65 Anaysis of Cyindrica Tanks with Fat Bases by Moment Distribution Methods T. SyopSis by AminGhai, HE moment distribution method is used to find the moments and the ring tension in the was and

More information

Experimental Investigation and Numerical Analysis of New Multi-Ribbed Slab Structure

Experimental Investigation and Numerical Analysis of New Multi-Ribbed Slab Structure Experimenta Investigation and Numerica Anaysis of New Muti-Ribbed Sab Structure Jie TIAN Xi an University of Technoogy, China Wei HUANG Xi an University of Architecture & Technoogy, China Junong LU Xi

More information

A Brief Introduction to Markov Chains and Hidden Markov Models

A Brief Introduction to Markov Chains and Hidden Markov Models A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,

More information

Bending Analysis of Continuous Castellated Beams

Bending Analysis of Continuous Castellated Beams Bending Anaysis of Continuous Casteated Beams * Sahar Eaiwi 1), Boksun Kim ) and Long-yuan Li 3) 1), ), 3) Schoo of Engineering, Pymouth University, Drake Circus, Pymouth, UK PL4 8AA 1) sahar.eaiwi@pymouth.ac.uk

More information

Vibrations of beams with a variable cross-section fixed on rotational rigid disks

Vibrations of beams with a variable cross-section fixed on rotational rigid disks 1(13) 39 57 Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks Abstract The work is focused on the probem of vibrating beams with a variabe cross-section fixed on a rotationa

More information

Applied Nuclear Physics (Fall 2006) Lecture 7 (10/2/06) Overview of Cross Section Calculation

Applied Nuclear Physics (Fall 2006) Lecture 7 (10/2/06) Overview of Cross Section Calculation 22.101 Appied Nucear Physics (Fa 2006) Lecture 7 (10/2/06) Overview of Cross Section Cacuation References P. Roman, Advanced Quantum Theory (Addison-Wesey, Reading, 1965), Chap 3. A. Foderaro, The Eements

More information

Identification of macro and micro parameters in solidification model

Identification of macro and micro parameters in solidification model BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vo. 55, No. 1, 27 Identification of macro and micro parameters in soidification mode B. MOCHNACKI 1 and E. MAJCHRZAK 2,1 1 Czestochowa University

More information

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method. Version 2 CE IIT, Kharagpur

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method. Version 2 CE IIT, Kharagpur odue 2 naysis of Staticay ndeterminate Structures by the atri Force ethod Version 2 E T, Kharagpur esson 12 The Three-oment Equations- Version 2 E T, Kharagpur nstructiona Objectives fter reading this

More information

UNCOMPLICATED TORSION AND BENDING THEORIES FOR MICROPOLAR ELASTIC BEAMS

UNCOMPLICATED TORSION AND BENDING THEORIES FOR MICROPOLAR ELASTIC BEAMS 11th Word Congress on Computationa Mechanics WCCM XI 5th European Conference on Computationa Mechanics ECCM V 6th European Conference on Computationa Fuid Dynamics ECFD VI E. Oñate J. Oiver and. Huerta

More information

Combining reaction kinetics to the multi-phase Gibbs energy calculation

Combining reaction kinetics to the multi-phase Gibbs energy calculation 7 th European Symposium on Computer Aided Process Engineering ESCAPE7 V. Pesu and P.S. Agachi (Editors) 2007 Esevier B.V. A rights reserved. Combining reaction inetics to the muti-phase Gibbs energy cacuation

More information

SEA Subsystem Attribute Correction Based on Stiffness Multipliers

SEA Subsystem Attribute Correction Based on Stiffness Multipliers Internationa Conference on Eectrica, Mechanica and Industria Engineering (ICEMIE 016) SEA Subsystem Attribute Correction Based on Stiffness Mutipiers Jie Gao1,* and Yong Yang 1 Coege of hysics and Information

More information

Modal analysis of a multi-blade system undergoing rotational motion

Modal analysis of a multi-blade system undergoing rotational motion Journa of Mechanica Science and Technoogy 3 (9) 5~58 Journa of Mechanica Science and Technoogy www.springerin.com/content/738-494x DOI.7/s6-9-43-3 Moda anaysis of a muti-bade system undergoing rotationa

More information

A. Distribution of the test statistic

A. Distribution of the test statistic A. Distribution of the test statistic In the sequentia test, we first compute the test statistic from a mini-batch of size m. If a decision cannot be made with this statistic, we keep increasing the mini-batch

More information

Interim Exam 1 5AIB0 Sensing, Computing, Actuating , Location AUD 11

Interim Exam 1 5AIB0 Sensing, Computing, Actuating , Location AUD 11 Interim Exam 1 5AIB0 Sensing, Computing, Actuating 3-5-2015, 14.00-15.00 Location AUD 11 Name: ID: This interim exam consists of 1 question for which you can score at most 30 points. The fina grade for

More information

Physics 506 Winter 2006 Homework Assignment #6 Solutions

Physics 506 Winter 2006 Homework Assignment #6 Solutions Physics 506 Winter 006 Homework Assignment #6 Soutions Textbook probems: Ch. 10: 10., 10.3, 10.7, 10.10 10. Eectromagnetic radiation with eiptic poarization, described (in the notation of Section 7. by

More information

More Scattering: the Partial Wave Expansion

More Scattering: the Partial Wave Expansion More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction

More information

Work and energy method. Exercise 1 : Beam with a couple. Exercise 1 : Non-linear loaddisplacement. Exercise 2 : Horizontally loaded frame

Work and energy method. Exercise 1 : Beam with a couple. Exercise 1 : Non-linear loaddisplacement. Exercise 2 : Horizontally loaded frame Work and energy method EI EI T x-axis Exercise 1 : Beam with a coupe Determine the rotation at the right support of the construction dispayed on the right, caused by the coupe T using Castigiano s nd theorem.

More information

Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.

Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n. Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

BP neural network-based sports performance prediction model applied research

BP neural network-based sports performance prediction model applied research Avaiabe onine www.jocpr.com Journa of Chemica and Pharmaceutica Research, 204, 6(7:93-936 Research Artice ISSN : 0975-7384 CODEN(USA : JCPRC5 BP neura networ-based sports performance prediction mode appied

More information

Introduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled.

Introduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled. imuation of the acoustic fied produced by cavities using the Boundary Eement Rayeigh Integra Method () and its appication to a horn oudspeaer. tephen Kirup East Lancashire Institute, Due treet, Bacburn,

More information

High Efficiency Development of a Reciprocating Compressor by Clarification of Loss Generation in Bearings

High Efficiency Development of a Reciprocating Compressor by Clarification of Loss Generation in Bearings Purdue University Purdue e-pubs Internationa Compressor Engineering Conference Schoo of Mechanica Engineering 2010 High Efficiency Deveopment of a Reciprocating Compressor by Carification of Loss Generation

More information

Instructional Objectives:

Instructional Objectives: Instructiona Objectives: At te end of tis esson, te students soud be abe to understand: Ways in wic eccentric oads appear in a weded joint. Genera procedure of designing a weded joint for eccentric oading.

More information

14 Separation of Variables Method

14 Separation of Variables Method 14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt

More information

CHAPTER XIII FLOW PAST FINITE BODIES

CHAPTER XIII FLOW PAST FINITE BODIES HAPTER XIII LOW PAST INITE BODIES. The formation of shock waves in supersonic fow past bodies Simpe arguments show that, in supersonic fow past an arbitrar bod, a shock wave must be formed in front of

More information

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA ON THE SYMMETRY OF THE POWER INE CHANNE T.C. Banwe, S. Gai {bct, sgai}@research.tecordia.com Tecordia Technoogies, Inc., 445 South Street, Morristown, NJ 07960, USA Abstract The indoor power ine network

More information

Lecture 17 - The Secrets we have Swept Under the Rug

Lecture 17 - The Secrets we have Swept Under the Rug Lecture 17 - The Secrets we have Swept Under the Rug Today s ectures examines some of the uirky features of eectrostatics that we have negected up unti this point A Puzze... Let s go back to the basics

More information

Research on liquid sloshing performance in vane type tank under microgravity

Research on liquid sloshing performance in vane type tank under microgravity IOP Conference Series: Materias Science and Engineering PAPER OPEN ACCESS Research on iquid soshing performance in vane type tan under microgravity Reated content - Numerica simuation of fuid fow in the

More information

Nonperturbative Shell Correction to the Bethe Bloch Formula for the Energy Losses of Fast Charged Particles

Nonperturbative Shell Correction to the Bethe Bloch Formula for the Energy Losses of Fast Charged Particles ISSN 002-3640, JETP Letters, 20, Vo. 94, No., pp. 5. Peiades Pubishing, Inc., 20. Origina Russian Text V.I. Matveev, D.N. Makarov, 20, pubished in Pis ma v Zhurna Eksperimenta noi i Teoreticheskoi Fiziki,

More information

SECTION A. Question 1

SECTION A. Question 1 SECTION A Question 1 (a) In the usua notation derive the governing differentia equation of motion in free vibration for the singe degree of freedom system shown in Figure Q1(a) by using Newton's second

More information

Introduction to LMTO method

Introduction to LMTO method 1 Introduction to MTO method 24 February 2011; V172 P.Ravindran, FME-course on Ab initio Modeing of soar ce Materias 24 February 2011 Introduction to MTO method Ab initio Eectronic Structure Cacuations

More information

Assignment 7 Due Tuessday, March 29, 2016

Assignment 7 Due Tuessday, March 29, 2016 Math 45 / AMCS 55 Dr. DeTurck Assignment 7 Due Tuessday, March 9, 6 Topics for this week Convergence of Fourier series; Lapace s equation and harmonic functions: basic properties, compuations on rectanges

More information

> 2 CHAPTER 3 SLAB 3.1 INTRODUCTION 3.2 TYPES OF SLAB

> 2 CHAPTER 3 SLAB 3.1 INTRODUCTION 3.2 TYPES OF SLAB CHAPTER 3 SLAB 3. INTRODUCTION Reinforced concrete sabs are one of the most widey used structura eements. In many structures, in addition to providing a versatie and economica method of supporting gravity

More information

Physics 111. Thursday, Dec. 9, 3-5pm and 7-9pm. Announcements. Tuesday, December 7, Stress Strain. For the rest of the semester

Physics 111. Thursday, Dec. 9, 3-5pm and 7-9pm. Announcements. Tuesday, December 7, Stress Strain. For the rest of the semester ics day, ember 7, 004 Ch 17: Kinetic Theory Stress Strain Ch 18: 1st Law of Thermodynamics nd Law of Thermodynamics or the rest of the semester Thursday,. 9, 3-5pm and 7-9pm Monday,. 13, 004 10:30 am 1:30

More information

Lecture 8 February 18, 2010

Lecture 8 February 18, 2010 Sources of Eectromagnetic Fieds Lecture 8 February 18, 2010 We now start to discuss radiation in free space. We wi reorder the materia of Chapter 9, bringing sections 6 7 up front. We wi aso cover some

More information

1 Equivalent SDOF Approach. Sri Tudjono 1,*, and Patria Kusumaningrum 2

1 Equivalent SDOF Approach. Sri Tudjono 1,*, and Patria Kusumaningrum 2 MATEC Web of Conferences 159, 01005 (018) IJCAET & ISAMPE 017 https://doi.org/10.1051/matecconf/01815901005 Dynamic Response of RC Cantiever Beam by Equivaent Singe Degree of Freedom Method on Eastic Anaysis

More information

Forces of Friction. through a viscous medium, there will be a resistance to the motion. and its environment

Forces of Friction. through a viscous medium, there will be a resistance to the motion. and its environment Forces of Friction When an object is in motion on a surface or through a viscous medium, there wi be a resistance to the motion This is due to the interactions between the object and its environment This

More information

Bohr s atomic model. 1 Ze 2 = mv2. n 2 Z

Bohr s atomic model. 1 Ze 2 = mv2. n 2 Z Bohr s atomic mode Another interesting success of the so-caed od quantum theory is expaining atomic spectra of hydrogen and hydrogen-ike atoms. The eectromagnetic radiation emitted by free atoms is concentrated

More information

Candidate Number. General Certificate of Education Advanced Level Examination January 2012

Candidate Number. General Certificate of Education Advanced Level Examination January 2012 entre Number andidate Number Surname Other Names andidate Signature Genera ertificate of Education dvanced Leve Examination January 212 Physics PHY4/1 Unit 4 Fieds and Further Mechanics Section Tuesday

More information

Numerical Simulation for Optimizing Temperature Gradients during Single Crystal Casting Process

Numerical Simulation for Optimizing Temperature Gradients during Single Crystal Casting Process ISIJ Internationa Vo 54 (2014) No 2 pp 254 258 Numerica Simuation for Optimizing Temperature Gradients during Singe Crysta Casting Process Aeksandr Aeksandrovich INOZEMTSEV 1) Aeksandra Sergeevna DUBROVSKAYA

More information

The Bending of Rectangular Deep Beams with Fixed at Both Ends under Uniform Load

The Bending of Rectangular Deep Beams with Fixed at Both Ends under Uniform Load Engineering,,, 8-9 doi:.6/eng..7 Pubised Onine December (ttp://.scirp.org/journa/eng) Te Bending of Rectanguar Deep Beams it Fied at Bot Ends under Uniform Load Abstract Ying-Jie Cen, Bao-Lian Fu, Gang

More information

Estimation of Lamina Stiffness and Strength of Quadriaxial Non-Crimp Fabric Composites Based on Semi-Laminar Considerations

Estimation of Lamina Stiffness and Strength of Quadriaxial Non-Crimp Fabric Composites Based on Semi-Laminar Considerations appied sciences Artice Estimation of Lamina Stiffness and Strength of Quadriaxia Non-Crimp Fabric Composites Based on Semi-Laminar Considerations Yong Cao, Yunwen Feng *, Wenzhi Wang *, Danqing Wu and

More information

Chemical Kinetics Part 2

Chemical Kinetics Part 2 Integrated Rate Laws Chemica Kinetics Part 2 The rate aw we have discussed thus far is the differentia rate aw. Let us consider the very simpe reaction: a A à products The differentia rate reates the rate

More information

Research Article Solution of Point Reactor Neutron Kinetics Equations with Temperature Feedback by Singularly Perturbed Method

Research Article Solution of Point Reactor Neutron Kinetics Equations with Temperature Feedback by Singularly Perturbed Method Science and Technoogy of Nucear Instaations Voume 213, Artice ID 261327, 6 pages http://dx.doi.org/1.1155/213/261327 Research Artice Soution of Point Reactor Neutron Kinetics Equations with Temperature

More information

XSAT of linear CNF formulas

XSAT of linear CNF formulas XSAT of inear CN formuas Bernd R. Schuh Dr. Bernd Schuh, D-50968 Kön, Germany; bernd.schuh@netcoogne.de eywords: compexity, XSAT, exact inear formua, -reguarity, -uniformity, NPcompeteness Abstract. Open

More information

An approximate method for solving the inverse scattering problem with fixed-energy data

An approximate method for solving the inverse scattering problem with fixed-energy data J. Inv. I-Posed Probems, Vo. 7, No. 6, pp. 561 571 (1999) c VSP 1999 An approximate method for soving the inverse scattering probem with fixed-energy data A. G. Ramm and W. Scheid Received May 12, 1999

More information

THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE

THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE KATIE L. MAY AND MELISSA A. MITCHELL Abstract. We show how to identify the minima path network connecting three fixed points on

More information

Faculty of Machine Building. Technical University of Cluj Napoca

Faculty of Machine Building. Technical University of Cluj Napoca Facuty of Machine Buiding Technica University of Cuj Napoca CONTRIBUTIONS TO THE CALCULATION AND ANALYSIS OF DYNAMIC ABSORBERS, WITH APPLICATIONS ON BALANCING MECHANICAL SYSTEMS PhD THESIS 11 Scientific

More information