Transcription

1 JOURNAL OF THEORETICAL AND APPLIED MECHANICS 51, 4, pp , Warsaw 2013 ENERGY CONSUMPTION IN MECHANICAL SYSTEMS USING A CERTAIN NONLINEAR DEGENERATE MODEL Krzysztof Jamroziak Gen. Tadeusz Kosciuszko Military Academy of Land Forces, Wrocław, Poland Miroslaw Bocian, Maciej Kulisiewicz Wroclaw University of Technology, The Institute of Materials Science and Applied Mechanics, Wrocław, Poland In modern engineering materials used for creating effective ballistic shields, the issues of evaluation of their energy consumption are extremely important. The paper presents a new way of solving this problem using a certain degenerate model with dry friction. This method involves the use of specially derived identification equations which describe the decrease in potential energy of the system during its vibratory motion induced by a single pulse load. Analytical considerations have been verified using a computer simulation technique for selected examples. Key words: composite ballistic shield, modeling, degenerated model, dissipation, impact energy 1. Introduction The amount of the energy that is absorbed in mechanical vibrating systems is usually described byasingleparameterthatisrelatedtotheadoptedmodelinformofpureviscousfriction(the friction that is proportional to the rate of deformation of the material) or so-called dry friction. In the case of the ballistic impact on shields which are made of modern(lightweight) materials, this issue becomes very complicated, because of complex and time-related strains that occur during penetration of the projectile into the shield. The pierced material is subjected to varying degrees of shearing, tension and compression, which depend primarily on the impact velocity aswellasontheshapeandmechanicalpropertiesoftheshieldandtheprojectile.thisissue was comprehensively discussed in Bourke s(2007) studies. The problem of dissipation of kinetic energy of the projectile in multilayer materials, such as laminates and composites, is currently analyzed by many researchers, such as Abrate(1998, 2010), Sanchez-Galvez et al.(2005), García- Castilloetal.(2012),Sidneyetal.(2011),Houetal.(2010),Katzetal.(2008).Themathematical approachtodissipationoftheimpactenergybasedontheseworksisanoriginalapproachtothe problem. Tabiei and Nilakantan(2008) presented the in-depth description of this phenomena through the analysis of the literature on the subject matter. The authors have synthetically presented the previous areas of the research conducted by worldwide scientists. In the papers by Jamroziak and Bocian(2008), Kulisiewicz et al.(2008), Jamroziak et al.(2010), the issue of dissipation of the impact energy was presented using the degenerate models. ItisgenerallyassumedthattheworkAdonebytheprojectileduringtheprocessofpiercing maybedescribedbydropofitskineticenergye,startingfromthezeroposition(theimpact velocity)uptothemomentitleavestheshieldorstopsintheshield.therelationsthatdescribe piercing are derived on the basis of the a priori assumed models of the constitutive relations (stress strain relations) which are highly complex for this type of materials, as it can be seen inthepapersbyjachetal.(2004),rusińskietal.(2005),buchmayretal.(2008).someofthe

2 828 K. Jamroziak et al. assumptions which are adopted for the description of these models were also included in the paper by Iluk(2012). Buchacz and his team have been conducting a long-time research aiming to develop a mathematical algorithm of analysis and synthesis of simple and complex mechanical and mechatronic systems. To realize these tasks, different categories of graphs and structural numbers were proposed by Białas(2008), Buchacz(1995), Buchacz and Wojnarowski(1995), Buchacz and Płaczek(2009). The studies included also computer-aided methods of realizing these tasks, Buchacz(2005). Vibrating mechatronic systems with piezoelectric transducers used to damp or induce vibrations were modelled and analysed in the papers by Baier and Lubczyński(2009),BuchaczandWróbel(2010).Theaimoftheseworkswastoidentifytheoptimal (according to the adopted criteria) mathematical model of the analysed systems as well as to develop mathematical tools useful to analyse these systems using approximate methods, see Białas(2010, 2012), Buchacz and Płaczek(2010), Wróbel(2012), Żółkiewski(2010, 2011) and the paper by Kulisiewicz et al.(2001) presenting the balance methods. Most of the hypotheses concerning the problem of dissipation of the impact energy take as the starting point the well known law of conservation of energy. This analysis was presented in the papers by Włodarczyk (2006), Włodarczyk and Jackowski(2008), Carlucci and Jacobson(2008). 2. Formulation of the problem The basic assumption made by the authors is that the lightweight shield acts on the piercing mass m with the resisting force S, whoose functional form is based on the mathematical analysis ofthedynamicmodel.thispaperassumesthemodelpresentedontheschemeshowninfig.1b. The standard model consists of the Maxwell element in parallel configuration with a purely elasticelementcandtheelemenththatdescribesdryfriction.itmaybenotedthattheadoption oftheconstant inthissystemresultsinobtainingthewidelyuseddynamicmodelwhich describes the vibrations of one-degree-of-freedom mechanical systems with dry friction(fig. 1a). Similarly,if c=0,theobtainedmodeltakesformofthepurelymaxwellelementinparallel configurationwiththeelementh.inthissense,thesystempresentedinfig.1bistheuniversal model, which should accurately describe the mechanical properties of many modern construction materials.theintroductionoftheelementhofdryfrictioninbothmodelshasbeenbasedon the results of the previous research of the authors, which concerned the impact process. This research was presented in the papers by Bocian et al.(2009), Jamroziak et al.(2009), Jamroziak and Bocian(2010). Fig. 1. The scheme of the analyzed dynamic models:(a) typical model with dry friction;(b) standard model(zener model) with dry friction

5 Energy consumption in mechanical systems Inaddition,for,equation(4.3)isidenticaltoequation(3.4) 1.Becauseofthat,the energy losses described by equation(4.3) are more complete and, therefore, this equation should be used in practice. Some exemplary results of computer simulations are described below. 5. Thesimulations Looking for solutions to differential equations(3.2) and(4.1), which describe vibrations of the analyzed models using the Mathematica software, simulation studies were performed. The followingvalueshavebeensubstitutedintoequations(3.3) 2,3 and(3.3) 5 ofmodel1aandinto equations(4.2) 1,2 and(4.2) 3 ofmodel1b: Model1a:k=480kg/s,c=30000kg/s 2,m=40kg,h=5 Model1b:k=480kg/s,c=30000kg/s 2, =20000kg/s 2,m=40kg,k 0 =406kg/s, h=5. Eachcaseincludedsimulationofthepulseloadintheformof: fort<0.1s wheretheforcep(t)hasbeenmodeledbyfunctionp(t)=asin(10πt). Responses of the models for pulse loads and different parameters are illustrated in the following figures(figs.4 7).Examplesoftheappliedexcitingforcep(t)areshowninFig.4.Thepulses wereone-sidedandtheirassumedtimeofdurationwasequalto t 0 =0.1s.Examplesofthe obtainedresponsesareshowninfig.5formodel(a)andinfig.6formodel(b).thephase trajectories for both models are presented in Fig. 7. Fig.4.Pulseloadsp(t):(a)formodel1a,(b)formodel1b Fig.5.Responsestothepulseloadsformodel1a,(a)velocityv=v(t),(b)displacementx=x(t)

6 832 K. Jamroziak et al. Fig.6.Responsestothepulseloadsformodel1b,(a)velocityv=v(t),(b)displacementx=x(t) Fig.7.Trajectoriesofthephaseanalysisforsimulationinthetimeinterval t:(a)formodel1a,(b)for model 1b The loop fields have been determined on the basis of the results obtained during simulations. Then, the values of coefficients from the identification equation(table 1) were generated using the linear regression. As it can be seen, these equations are generally satisfied, although in one caseofmodel(b)theestimatedvalueofthecoefficientcontainsabigmistake.inothercases, the mistakes were no greater than 10%. Table 1. The parameters assumed and derived from the linear regression Model(a) Model(b) Assumed Derived Assumed Derived 2k 2k 2k 0 m c = c = = k 0 m = c c 2k 0 ( +c) = k 0 ( +c) = c c 2h c = h c = h c = h c = Summary The presented fragment of the research work concerns the analysis of the energy consumption of dynamic rheological models in the process of ballistic impact. Two dynamic models have been adopted for this analysis. The first model describes vibrations of mechanical systems with one

7 Energy consumption in mechanical systems degree of freedom, and the second one describes vibrations of mechanical systems with one and ahalfdegreeoffreedom.thisisamodelfromagroupofthedegeneratemodels.theanalysis of the models have been conducted based on the mathematical relations which determine the dissipation of the impact energy. The models have also been subjected to computer simulation in order to verify the theoretical assumptions. The simulations resulted in obtaining time responses for the pulse force p(t). It hasbeenassumedthatatthetimepoints and t 2,thedisplacement x(t)shouldreachthe extremevaluesforthepointsoftimeinwhichthevelocity v(t)musthavethevalueofzero. The simulation has confirmed the expected objectives. Indeed, for the velocity v(t)(fig. 5), at moments its reaches zero, the displacement in these moments x(t) reaches its maximum(fig. 6). Itcanbealsonotedthatinthecaseofthedegeneratemodel,thedisplacementx(t)takesmuch lower values than in model 1a. These differences occur because the responses of this model are suppressed more effectively. This also shows that the process of energy consumption is described more accurately by the degenerate model. The obtained shapes of the phase trajectory of the analyzed systems in the time interval t shows also some significant differences. Drawing conclusions requires still some additional simulations at this stage. The presented diagrams provide rather the qualitative description of these phenomena. To reach a quantitative description, a number of research studies still have to be carried out, and their results will be presented in the future papers. To sum up, the hypothesis assumed by the authors that the degenerate models can quite accurately describe the mechanical properties of modern structural materials has been confirmed and justifies the direction of the undertaken work. References 1. Abrate S., 1998, Impact on Composite Structures, Cambrige University Press, Cambrigde 2. Abrate S., 2010, Ballistic Impacts on Composite and Sandwich Structures, Major Accomplishments in Composite Materials and Sandwich Structures, DOI: / Baier A., Lubczyński T., 2009, Virtual simulation of mechatronics laboratory, Solid State Phenomena, , Białas K., 2008, Polar graphs and structural numbers in synthesis of active and passive mechanical systems, Journal of Achievements in Materials and Manufacturing Engineering, 30, 1, Białas K., 2010, Passive and active elements in reduction of vibrations of torsional systems, Solid State Phenomena, 164, Białas K., 2012, Mechanical and electrical elements in reduction of vibrations, Journal of Vibroengineering, 14, 1, Bocian M., Jamroziak K., Kulisiewicz M., Piesiak S., 2009, Analysis of absorbsing energy process in light ballistic shield(in Polish), Zeszyty Naukowe Politechniki Świętokrzyskiej Budowa i Eksploatacja Maszyn, 12, 23-24, full text CD-ROM 8. Bourke P., 2007, Ballistic Impact on Composite Armour, Cambrige University Press, Shrivenharn Wiltshire 9. Buchacz A., 1995, Modeling, synthesis and analysis of bar systems characterized by a cascade structure represented by graphs, Mechanism and Machine Theory, 30, 7, Buchacz A., 2005, The expansion of the synthesized structures of mechanical discrete systems represented by polar graphs, Journal of Materials Processing Technology, 164, Buchacz A., Płaczek M., 2009, Damping of mechanical vibrations using piezoelements, including influence of connection layer s properties on the dynamic characteristic, Solid State Phenomena, 147/149,

APPLIED MATHEMATICS ADVANCED LEVEL INTRODUCTION This syllabus serves to examine candidates knowledge and skills in introductory mathematical and statistical methods, and their applications. For applications

CHAPTER 0 Introduction to Engineering System Dynamics 0.1 INTRODUCTION The objective of an engineering analysis of a dynamic system is prediction of its behaviour or performance. Real dynamic systems are

Physics 53 Oscillations You've got to be very careful if you don't know where you're going, because you might not get there. Yogi Berra Overview Many natural phenomena exhibit motion in which particles

The Big Idea: In most realistic situations forces and accelerations are not fixed quantities but vary with time or displacement. In these situations algebraic formulas cannot do better than approximate

Aim : To study how the time period of a simple pendulum changes when its amplitude is changed. Teacher s Signature Name: Suvrat Raju Class: XIID Board Roll No.: Table of Contents Aim..................................................1

3 GENERATORS An important application area for is in the conversion of mechanical energy into electrical energy, and this chapter describes the conditions under which should be used to convert the maximum

Quite often the measurements of pressures has to be conducted in unsteady conditions. Typical cases are those of -the measurement of time-varying pressure (with periodic oscillations or step changes) -the

SAFE A HEAD Structural analysis and Finite Element simulation of an innovative ski helmet Prof. Petrone Nicola Eng. Cherubina Enrico Goal Development of an innovative ski helmet on the basis of analyses

An Introduction to Piezoelectric Transducer Crystals Piezoelectric Materials and their Properties Certain single crystal materials exhibit the following phenomenon: when the crystal is mechanically strained,

OPTIMAL ANALYSIS, DESIGN AND TESTING OF AN ELEC- TROMAGNETIC DAMPER WITH RESONANT SHUNT CIRCUIT FOR VIBRATION CONTROL OF A CIVIL STRUCTURE Wai Kei Ao, Paul Reynolds University of Exeter, Vibration Engineering

Chapter 4: Summary and Conclusions 4.1 Summary Three different models are presented and analyzed in this research for the purpose of studying the potential of using post-buckled or pre-bent elastic struts

1. Which of the following statements about a spring-block oscillator in simple harmonic motion about its equilibrium point is false? (A) The displacement is directly related to the acceleration. (B) The

.1.1 Measure the motion of objects to understand.1.1 Develop graphical, the relationships among distance, velocity and mathematical, and pictorial acceleration. Develop deeper understanding through representations

Fall 08 Prepared by: Keivan Anbarani Abstract In this experiment four eccentric masses are used in conjunction with four springs and one damper to simulate the vibration. Masses are aligned in different

The Ship Power Supplier Guidelines to engine dynamics and vibration by Hannu Tienhaara Technology Wärtsilä Corporation In recent years shipyards and shipowners have increasingly tuned in to the negative

Section 8. 8. Elastic Strain Energy The strain energy stored in an elastic material upon deformation is calculated below for a number of different geometries and loading conditions. These expressions for

Journal of KONES Powertrain and Transport, Vol. 18, No. 4 2011 COMPARATIVE ANALYSIS OF THE VIBRATIONS OF A DIFFERENT KIND OF ENGINE MOUNTED IN THE SAME NEW MOTOR VEHICLES Lech Sitnik, Monika Magdziak Tok

1.1 Using Figure 1.6, verify that equation (1.1) satisfies the initial velocity condition. Solution: Following the lead given in Example 1.1., write down the general expression of the velocity by differentiating

11 FORCED OSCILLATIONS AND RESONANCE POINTER INSTRUMENTS Analogue ammeter and voltmeters, have CRITICAL DAMPING so as to allow the needle pointer to reach its correct position on the scale after a single

Assessment Chapter Test B Chapter: Measurements and Calculations PART I In the space provided, write the letter of the term or phrase that best completes each statement or best answers each question. 1.

Updated 2013 (Mathematica Version) M1.1 Introduction. Lab M1: The Simple Pendulum The simple pendulum is a favorite introductory exercise because Galileo's experiments on pendulums in the early 1600s are

dspace DSP DS-1104 based State Observer Design for Position Control of DC Servo Motor Jaswandi Sawant, Divyesh Ginoya Department of Instrumentation and control, College of Engineering, Pune. ABSTRACT This

SOLID MECHANICS TUTORIAL MECHANISMS KINEMATICS - VELOCITY AND ACCELERATION DIAGRAMS This work covers elements of the syllabus for the Engineering Council exams C105 Mechanical and Structural Engineering

The First Law of Thermodynamics: Closed Systems The first law of thermodynamics can be simply stated as follows: during an interaction between a system and its surroundings, the amount of energy gained

Development of New Inkjet Head Applying MEMS Technology and Thin Film Actuator Kenji MAWATARI, Koich SAMESHIMA, Mitsuyoshi MIYAI, Shinya MATSUDA Abstract We developed a new inkjet head by applying MEMS

w Technical Product Information Precision Miniature Load Cell with Overload Protection 1. Introduction The load cells in the model 8431 and 8432 series are primarily designed for the measurement of force

Advanced Higher Physics: MECHANICS Simple Harmonic Motion At the end of this section, you should be able to: Describe examples of simple harmonic motion (SHM). State that in SHM the unbalanced force is

Journal of Engineering Science and Technology Review 2 (1) (2009) 76-81 Lecture Note JOURNAL OF Engineering Science and Technology Review www.jestr.org Time of flight and range of the motion of a projectile

Dynamic Process Modeling Process Dynamics and Control 1 Description of process dynamics Classes of models What do we need for control? Modeling for control Mechanical Systems Modeling Electrical circuits

ERROR PROPAGATION For sums, differences, products, and quotients, propagation of errors is done as follows. (These formulas can easily be calculated using calculus, using the differential as the associated