SOME INVERSE PROBLEMS IN PENETRATION MECHANICS
|
|
- Moris Weaver
- 5 years ago
- Views:
Transcription
1 Mechanics Based Design of Structures and Machines, 38(4): , 21 Copyright Taylor & Francis Group, LLC ISSN: print/ online DOI: 1.18/ SOME INVERSE PROBLEMS IN PENETRATION MECHANICS G. Ben-Dor, A. Dubinsky, and T. Elperin Department of Mechanical Engineering, Pearlstone Center for Aeronautical Engineering Studies, Ben-Gurion University of the Negev, Beer-Sheva, Israel Many empirical and semi-empirical models that are widely used in impact dynamics employ the dependencies between the impact (initial) velocity and the depth of penetration (DOP). These dependencies are often obtained by statistical processing of experimental results and are not based on the theoretical models of penetration. In this study we develop an approach that allows using such phenomenological correlations for determining the dependencies of forces acting on a penetrator versus the instantaneous DOP and the instantaneous velocity. The developed approach can be also useful in other fields of mechanics. Keywords: Concrete; Granular media; High speed; Impact; Penetration. INTRODUCTION Inverse problems often occur in various fields of science and are a subject of numerous studies (see e.g., Denisov, 1994; Ramm, 25 and references therein). Solution of the inverse problems for ordinary differential equations (ODEs) implies determining either the right-hand side of ODEs or their coefficients (including the case when they are functions of the independent variable), using the known solution. However, inverse problems of different type also occur in applied sciences, and particular in penetration mechanics (Ambroso et al., 25; Eisler et al., 1998; Katsuragi and Durian, 27; Kennedy, 1976; Pilyugin, 24). Comprehensive theoretical foundations for investigating these problems have not been developed as yet. It is convenient to describe the essence of these problems without dwelling on a particular application. Let y = x y be a solution of the first-order ODE with the initial condition y = y for x =. Let us consider for y =ỹ the set of x 1 and y, which satisfy the equation ỹ = x 1 y. It is assumed that the empirical dependence x 1 y x 1 y ỹ = is known. Consider the problem of recovering the ODE, which includes the unknown functions such that the solution of this ODE implies the dependence, x 1 y =. In the case when ODE describes penetration of a striker Received April 16, 21; Accepted June 4, 21 Communicated by N.Banichuk. Correspondence: Tov Elperin, Department of Mechanical Engineering, Pearlstone Center for Aeronautical Engineering Studies, Ben-Gurion University of the Negev, P.O. Box 653, Beer-Sheva 8415, Israel; elperin@.bgu.ac.il 468
2 469 BEN-DOR ET AL. into a semi-infinite shield, x and y are the instantaneous depth of penetration (DOP) and the instantaneous velocity of the striker, respectively, and formula x 1 y = relates the initial velocity of the striker, y, and the maximum penetration depth, x 1, whereby the velocity of the striker vanishes (ỹ =. The goal is determining the dependence of the drag force acting at the striker as a function of the instantaneous velocity and the instantaneous DOP. In the present study we develop an approach for investigating inverse problems of this type, which was suggested by Ben-Dor et al. (29) for modeling penetration of strikers into various media. FORMULATION OF THE PROBLEM Consider normal penetration of a rigid projectile into a semi-infinite shield along the axis H, where the coordinate H determines the instantaneous location of the penetration and is defined as the distance between some reference point at the impactor and the front surface of the shield; the effects associated with the stage of the incomplete immersion of the nose of the impactor into the shield are neglected. Assuming that the total force acting at the impactor, R, depends on the velocity of the impactor, V, and the instantaneous coordinate of the impactor, H, Newton s second law describing the motion of the impactor with mass m can be written as follows: mv dv/dh = R H V (1) If V = V imp H is the solution of Eq. (1) with the initial condition V = V imp (V imp is the impact velocity) then the DOP, H max (a depth whereby the velocity of the penetrator vanishes), is determined by the following formula: V imp H max =. Assume that the empirical dependence, H max = F V imp, is known and F = F V imp > F V imp > V imp > (2) The problem is to determine a function R that implies this empirical dependence. Using the dimensionless variables Eq. (1) can be rewritten as follows: where h = H H v = V V w = v 2 = V 2 dw/dh = D h w (3) V 2 D h w = 2H R H mv 2 h V w (4) V is some reference velocity and H is a characteristic length. The empirical dependence between the squared dimensionless impact velocity, w imp, and the dimensionless DOP, h max, can be expressed in the following two forms: h max = F w w imp F w = or w imp = F h h max F h = (5)
3 SOME INVERSE PROBLEMS IN PENETRATION MECHANICS 47 where F h and F w are the known, mutually inverse functions, w imp = V imp /V 2 h max = H max /H (6) We consider expression for D h w that allows solving the equation of motion by quadratures. EQUATION WITH SEPARABLE VARIABLES General Solution Consider the case when the function D can be written as follows: D h w = D h h D w w (7) where D h and D w are some functions. Separating the variables in Eq. (3) and using initial conditions, w = w imp for h =, we obtain the equation of motion of the impactor. w w imp d w h D w w = D h h d h (8) Equation (8) can be written for the end point of the trajectory substituting w = and h = h max. wimp d w D w w = hmax D h h d h (9) Considering Eq. (9) as a definition of the functional dependence h max versus w imp and differentiating both sides of Eq. (8) with respect to w imp yields D w w imp D h h max dh max /dw imp = 1 (1) Substituting h max = F w w imp into Eq. (1) reduces Eq. (1) to the identity, D w w imp D h F w w imp F w w imp = 1, which must be satisfied for arbitrary w imp.itis convenient to rewrite this identity using the variable w instead of w imp. D w w D h F w w F w w = 1 (11) If function D h h in Eq. (7) is known then Eq. (11) implies the following expression for the unknown function, D w w : D w w = 1/ D h F w w F w w (12) Using the dependence between w imp and h max in the second form given by Eq.(5), Eq.(12) can be rewritten after the change of variables, F w x = F 1 h x, asd w F h h D h h F h h =. If D w w is known then D h h = F h h /D w F h h (13)
4 471 BEN-DOR ET AL. Equations (12) and (13) imply that the function D h w cannot be determined uniquely on the basis of the known correlation between the impact velocity and the DOP. Consequently, one can select either the function D w w or the function D h h. Substituting D w w from Eq. (12) or D h h from Eq. (13) into Eq. (7) we obtain two expressions for D. D h w = D h h D h F w w F w w D h w = D w w F h h (14) D w F h h In order to write the equation of motion of the impactor using only one unknown function, we substitute D w w from Eq. (12) into Eq. (8). wimp w h D h F w w F w w d w = D h h d h (15) After change of variable in the left-hand side, F w w = h, Eq. (15) can be rewritten as follows: Fw w imp F w w D h h d h= h D h h d h (16) Equation (16) implies that F w w imp F w w = h (17) where z = z D h h d h (18) Functions D h z and D w z can be replaced by the function z. Clearly (see Eq. (18)), and Eqs. (12) and (19) imply that D h z = z (19) D w z = 1/ F w z F w z (2) Therefore the instantaneous DOP as a function of the instantaneous velocity of the impactor as well as the instantaneous drag force can be represented using function in the following form: h w = 1( F w w imp F w w ) (21) h D h w = (22) F w w F w w where is an arbitrary function that satisfies the following conditions: = z > z > z > (23)
5 SOME INVERSE PROBLEMS IN PENETRATION MECHANICS 472 Using Eq. (21), D could be represented as a function of single argument, h or w, but corresponding expressions are cumbersome and are not presented here. Some Classes of Models Model z = z p If z = z p, where p>1 is a constant, then 1 z = z 1/p and Eqs. (21) and (22) can be written as follows: h w = F w w imp p F w w p 1/p (24) D h w = (25) F w w p 1 F w w Substituting h from Eq. (24) into Eq. (25) we obtain expression for D as a function of w. h p 1 D w = F w w imp p F w w p p 1 /p F w w p 1 F w w (26) Equations (24) and (25) remain valid for p = 1 after setting p = 1, h p 1 = F w h p 1 = 1. Formula (19) implies that q = 1 corresponds to the case D h w = 1 in Eq. (7). then Model z = F h z q F w w F w w z =F 1 w If z = F h z q = Fw 1 z q, where q>1is a constant, z = d F w w dw z =F 1 w z (27) = dwq dw = qwq 1 1 z = F w z 1/q and Eqs. (21) and (22) can be written as follows: Equation (27) implies that and the expression for D as a function of h reads ( h w = F w w q imp wq 1/q) (28) D h w = F h h q 1 F h h w1 q (29) w h = w q imp F h h q 1/q (3) D h = dw/dh = w q imp F h h q 1 q /q F h h q 1 F h h w1 q (31) Equations (28) and (31) remain valid for q = 1 after putting q = 1 and w 1 q = F h h q 1 = 1. Equations (14) and (27) imply that q = 1 corresponds to the case D w w = 1 in Eq. (7).
6 473 BEN-DOR ET AL. Particular Case Assume that the dependence between w imp and h max is given in an implicit form: w imp = G h max (32) where and G are positive, increasing functions for positive arguments and = G = (33) As previously we consider Eq. (9) as a definition of the functional dependence h max versus w imp. Equations (9) and (32) imply that wimp d w D w w = w imp hmax D h h d h = G h max (34) and, consequently, Using Eqs. (18), (32), and (35) we find that D w w = 1/ w D h h = G h (35) D h w = G h / w (36) z = G z F w w = G 1( w ) (37) Combining Eqs. (21) and (37) we arrive at the following relation: h w = G 1( w imp w ) (38) Note that Eq. (34) yields only one of the all possible models of penetration. The complete set of models can be obtained using the general method outlined in Section 3.1 and applying it to the dependence, F w w = G 1 w. High Speed Penetration into Semi-Infinite Concrete Shield NDRC Formula Consider widely known modified National Defence Research Committee (NDRC) formula (Kennedy, 1976), associated with penetration of a rigid projectile into a semi-infinite concrete shield which can be represented in the form given by Eq. (32) with w imp = wimp n n = 9 (39) { h G h max = max /2 2 if h max 2 (4) h max 1 if h max > 2 where the characteristic length H = d, d is the diameter of a projectile, the reference velocity, V, depends on d, the shape of the impactor s nose and mechanical properties of the shield.
7 SOME INVERSE PROBLEMS IN PENETRATION MECHANICS 474 Equation (32) allows us to write the following expression for the functions F w and F h : F w z = G 1 z = { 2z n/2 if z 1 z n + 1 if z>1 (41) { F h z = Fw 1 z z/2 2/n if z 2 z 1 1/n if z>2 (42) Model z = z p In the case of the model z = z p Eqs. (24) and (22) after substituting F w from Eq. (41) yield 2 w pn/2 imp wpn/2 1/p if w w imp for w imp < 1 {[ h w = w n imp + 1 p w n + 1 p]1/p if 1 w w imp (43) [ for w imp 1 w n imp + 1 p 2w n/2 p]1/p if w 1 h/2 p 1 w 2 pn /2/n if w 1 D h w = h p 1 w 1 n (44) if w>1 n w n + 1 p 1 where 1 p<2/n. Model z = F h z q If z = F h z q Eqs.(24) and (25) after substituting F w from Eq. (41) yield the following relations: 2 w { q imp wq n/ 2q if w w imp for w imp < 1 h w = 2 w q imp wq n/ 2q if w q imp w w 1 1/q imp (45) w q imp wq n/q + 1 if w w q imp for w imp 1 1 1/q { h/2 2q n /n w 1 q /n if h 2 D h w = (46) h 1 q n /n w 1 q /n if h>2 where n/2 <q 1. In the particular case q = n the model recovers the model suggested by Kennedy (1976) who derived it using the approach, which is considered in the Section Results of Calculations Figures 1(a) and (b) and 2(a) and (b) illustrate dependences w h and D h for models described in Eqs.(43) and(44) (solidcurves) and Eqs.(45) and(46) (dashcurves);correlationsd h are obtained by excluding w from Eqs. (43) (44) and (45) (46). Inspection of these figures shows that a large variety of different dependencies can be obtained by selecting the function.
8 475 BEN-DOR ET AL. Figure 1 (a) and (b). Instantaneous squared velocity of impactor versus instantaneous depth of penetration(dimensionless variables, penetration into concrete) for models described by Eqs.(43) - (46) LINEAR EQUATION General Solution Assume that the expression for the resistance force D can be written as follows: D h w = D 1 h w + D h (47) where D h and D 1 h are some functions. The solution of Eq.(3) with the initial condition w = w imp reads w h = ( w imp h [ D h / h ] ) d h h ( h = exp h ) D 1 h d h (48)
9 SOME INVERSE PROBLEMS IN PENETRATION MECHANICS 476 Figure 2 (a) and (b). Instantaneous resistance force versus instantaneous depth of penetration of impactor (dimensionless variables, penetration into concrete) for models described by Eqs.(43)-(46) Setting w = and h = h max in Eq. (48) we conclude that this model implies the following relationship between the DOP, h max, and the impact velocity: w imp = hmax D h / h d h (49) Considering w imp as a known function of h max, w imp = F h h max, differentiating Eq. (49) with respect to h max and replacing h max by h yields Equation (48) implies that D h = h F h h (5) D 1 h = h / h (51)
10 477 BEN-DOR ET AL. Therefore functions D h and D 1 h can be determined if one selects the function h, where h > = 1 (52) Substituting D h and D 1 h from Eqs. (5) and (51) into Eqs. (47) and (48) we obtain D h w = h / h w + h F h h (53) w h = w imp F h h h (54) Substituting w from Eq. (54) into Eq. (53) we obtain D h = h w imp F h h + h F h h = dw/dh (55) Low-Speed Penetration into Granular Media Consider a low-speed normal penetration of a rigid ball into a granular medium, and assume that the empirical dependence between the impact velocity and the DOP reads H max d b = k ( ) V 2 imp + H max (56) 2gd b d b where d b is the diameter of the ball and g is the acceleration of gravity. Parameters k and are determined by the choice of a particular model. In the model suggested by Uehara et al. (23), k = 14/ b / g, = 1/3, where g and b are densities of the grains and of the impacting ball, respectively, = tan is the grain grain friction coefficient, is the repose angle of the grains. In the model considered by Seguin et al. (28) which takes into account the influence of the confinement on penetration, k = A A b / g, =, where parameters A > and < < 1 correspond to the unbounded case, <A < 1 and < < 1 are the correction for the confinement and >1 is a parameter. We consider the case when the expression for the total resistance force, R, reads R H V = R 1 H V 2 + R H (57) where R and R 1 are some functions. Semi-empirical models of penetration into granular media are usually selected among the models of this type. Note that similar models are also widely used in high-speed penetration dynamics (Ben-Dor et al., 26). Then Newton s second law, Eq. (1), can be written in the form given by Eqs. (3) and (47) where D h = R H h D mg 1 h = 2H m R 1 H h (58)
11 SOME INVERSE PROBLEMS IN PENETRATION MECHANICS 478 Table 1 Three models that imply Eq. (59) for = 3. Model h Value of the parameter 1 exp c 1 h c 1 = / 1 + c 2 h 3/2 c 2 = c 3 h 2 c 3 = 3 H = d b k 1/ 1 is the initial collapse depth of the ball prior to the beginning of its motion from the state of rest (V imp =, variables h, v, w, w imp, h max are determined by Eqs. (4) and (6) and V = 2gH. Using the dimensionless variables Eq. (56) can be rewritten as follows: w imp = F h h max F h z = z z = 1/ (59) and Eqs. (53) and (55) imply the following dependencies for the drag force and squared instantaneous velocity versus the instantaneous DOP: D h = w imp h + h h + h h 1 1 (6) w h = ( w imp h + h ) h (61) Since the considered problem has a non-unique solution, it is feasible, as a rule, to generate different models having the same accuracy by selecting different h. To illustrate this, several versions of h are selected (see Table 1, columns 1 and 2) including Model 1 by Ambroso et al. (25) that can be readily obtained in the framework of the suggested approach. In comparison with the experiment, we use the approach suggested in Ambroso et al. (25) that is based on the following relationship that can be Figure 3 Comparison of models with experimental data for penetration into granular media.
12 479 BEN-DOR ET AL. obtained from Eqs. (59) and (61): ( ) V 2 ( = V imp 1 H 1 H 1 1 ) H/ (62) where H = h/h max = H/H max 1, = 1/h max = H /H max 1. Following Ambroso et al. (25), we determine parameters of the Models 2 and 3 using the relationship (62) on the basis of the same experimental data (495 points) for the same value = 3 and average value = 36. The best one-parameter fits to Eq. (62) for the Models 2 and 3 are shown in Fig. 3, as well as, for comparison, the curve corresponding to the Model 1 suggested by Ambroso et al. (25); the values of the parameters are presented in Table 1 (column 3). Inspection of this figure demonstrates that the experimental data are faithfully approximated by each of the models and the models have the same accuracy. CONCLUDING REMARKS We suggested an approach for solving some problems of penetration dynamics but the scope of application of the proposed method is considerably wider: this method can be used in modeling different phenomena, which are described by ODEs. The developed approach allows extension to other types of ODEs and other kinds of relationships that are assumed to be known. Non-uniqueness in the choice of the model, which implies the known empirical relationship, affords a wide range of possibilities for constructing the appropriate model of the considered phenomenon. We developed a general procedure for constructing a variety of models for selecting a particular model. The choice of the particular model is another problem, which is outside the scope of this study. It must be emphasized that the suggested method on its own does not guarantee the adequacy of the obtained equations. REFERENCES Ambroso, M. A., Kamien, R. D., Durian, D. J. (25). Dynamics of shallow impact cratering. Physical Review E 72:4135. Ben-Dor, G., Dubinsky, A., Elperin, T. (26). Applied High-Speed Plate Penetration Dynamics. Dordrecht: Springer. Ben-Dor, G., Dubinsky, A., Elperin, T. (29). Modeling of penetration by rigid impactors. Mechanics Research Communications 36(5): Denisov, A. M. (1994). Introduction to Theory of Inverse Problems [Vvedenije V Teoriju Obratnykh Zadach]. Moscow: Moscow State University. Eisler, R. D., Chatterjee, A. K., Burghart, G. H., Loan, P. (1998). Simulates the Tissue Damage from Small Arms Projectiles and Fragments Penetrating the Musculoskeletal System. Fountain Valley CA: Mission Res. Corp., Final Report. Katsuragi, H., Durian, D. J. (27). Unified force law for granular impact cratering. Nature Physics 3(6): Kennedy, R. P. (1976). A review of procedures for the analysis and design of concrete structures to resist missile impact effects. Nuclear Engineering and Design 37(2):
13 SOME INVERSE PROBLEMS IN PENETRATION MECHANICS 48 Pilyugin, N. N. (24). A simulation of the shape of a crater in an organic-glass target under high-velocity impact. High Temperature 42(3): Ramm, A. G. (25). Inverse Problems: Mathematical and Analytical Techniques with Applications to Engineering. New York: Springer. Seguin, A., Bertho, Y., Gondret, P. (28). Influence of confinement on granular penetration by impact. Physical Review E 78:131. Uehara, J., Ambroso, M., Ojha, R., Durian, D. (23). Low-speed impact craters in loose granular media. Physical Review Letters 9:19431.
Shape Optimization of Impactors Against a Finite Width Shield Using a Modified Method of Local Variations #
Mechanics Based Design of Structures and Machines, 35: 113 125, 27 Shape Optimization of Impactors Against a Finite Width Shield Using a Modified Method of Local Variations # G. Ben-Dor, A. Dubinsky, and
More informationOptimization of two-component composite armor against ballistic impact
Composite Structures 69 (005) 89 94 www.elsevier.com/locate/compstruct Optimization of two-component composite armor against ballistic impact G. Ben-Dor,. Dubinsky, T. Elperin * The Pearlstone Center for
More informationPeculiarities of first-order phase transitions in the presence of an electric field
Peculiarities of first-order phase transitions in the presence of an electric field Yu. Dolinsky* and T. Elperin The Pearlstone Center for Aeronautical Engineering Studies Department of Mechanical Engineering
More informationUnified force law for granular impact cratering
1 Unified force law for granular impact cratering H. Katsuragi and D. J. Durian Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, Pennsylvania 19104-6396, USA Experiments on
More informationShape optimization of high-speed penetrators: a review. Author copy
Cent. Eur. J. Eng. 2(4) 212 473-482 DOI: 1.2478/s13531-12-22-4 Central European Journal of Engineering Shape optimization of high-speed penetrators: a review Review Article Gabi Ben-Dor, Anatoly Dubinsky,
More informationEffects of abrasion on the penetration of ogival-nosed projectiles into concrete targets
7(2010) 413 422 Effects of abrasion on the penetration of ogival-nosed projectiles into concrete targets Abstract This paper investigates the effects of abrasion on the penetration of an ogival-nosed projectile
More informationQUANTUM BOUND STATES WITH ZERO BINDING ENERGY
ψ hepth@xxx/940554 LA-UR-94-374 QUANTUM BOUND STATES WITH ZERO BINDING ENERGY Jamil Daboul Physics Department, Ben Gurion University of the Negev Beer Sheva, Israel and Michael Martin Nieto 2 Theoretical
More informationEquilibrium orientation of an ellipsoidal particle inside a dielectric medium with a finite electric conductivity in the external electric field
PHYSICAL REVIEW E 71, 056611 005 Equilibrium orientation of an ellipsoidal particle inside a dielectric medium with a finite electric conductivity in the external electric field Yu. Dolinsky* and T. Elperin
More informationFundamentals Physics. Chapter 10 Rotation
Fundamentals Physics Tenth Edition Halliday Chapter 10 Rotation 10-1 Rotational Variables (1 of 15) Learning Objectives 10.01 Identify that if all parts of a body rotate around a fixed axis locked together,
More informationBallistic Impact: Recent Advances in Analytical Modeling of Plate Penetration Dynamics A Review
Anatoly Dubinsky's Home page: http://www.bgu.ac.il/~dubin Gabi Ben-Dor e-mail: bendorg@bgu.ac.il Anatoly Dubinsky e-mail: dubin@bgu.ac.il Tov Elperin e-mail: elperin@bgu.ac.il The Pearlstone Center for
More informationDevelopment of Empirical formula Prediction on Critical Impact Energy for Perforation Phenomena on Concrete Structures
Development of Empirical formula Prediction on Critical Impact Energy for Perforation Phenomena on Concrete Structures Qadir Bux alias Imran Latif Faculty of Civil and Environmental Engineering, Universiti
More informationAP PHYSICS 1 Content Outline arranged TOPICALLY
AP PHYSICS 1 Content Outline arranged TOPICALLY with o Big Ideas o Enduring Understandings o Essential Knowledges o Learning Objectives o Science Practices o Correlation to Common Textbook Chapters Much
More informationMixing at the External Boundary of a Submerged Turbulent Jet
Mixing at the External Boundary of a Submerged Turbulent Jet A. Eidelman, T. Elperin, N. Kleeorin, I. Rogachevskii, I. Sapir-Katiraie The Ben-Gurion University of the Negev, Beer-Sheva, Israel G. Hazak
More informationOrdinary Differential Equations
Chapter 10 Ordinary Differential Equations 10.1 Introduction Relationship between rate of change of variables rather than variables themselves gives rise to differential equations. Mathematical formulation
More informationLecture 10 Penetration Mechanics I (references are in lecture 11 notes) 1. Rigid Penetration of Targets
Lecture 10 Penetration Mechanics I (references are in lecture 11 notes) 1. Rigid Penetration of Targets Projectiles may be rigid or deforming. Rigid projectiles include nails, hard bullets, and low velocity
More informationChapter 8: Particle Systems and Linear Momentum
Chapter 8: Particle Systems and Linear Momentum Up to this point in our study of classical mechanics, we have studied primarily the motion of a single particle or body. To further our comprehension of
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics. Physics 8.01 Fall Problem Set 2: Applications of Newton s Second Law Solutions
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics Physics 8.01 Fall 2012 Problem 1 Problem Set 2: Applications of Newton s Second Law Solutions (a) The static friction force f s can have a magnitude
More informationModelling of a Tennis Ball Server
Mechanical Analysis and Design ME 2104 Lecture 6 Modelling of a Tennis Ball Server Prof Ahmed Kovacevic Department of Mechanical Engineering and Aeronautics Room CG25, Phone: 8780, E-Mail: a.kovacevic@city.ac.uk
More informationChapter 6. Circular Motion and Other Applications of Newton s Laws
Chapter 6 Circular Motion and Other Applications of Newton s Laws Circular Motion Two analysis models using Newton s Laws of Motion have been developed. The models have been applied to linear motion. Newton
More informationarxiv: v1 [cond-mat.soft] 10 Dec 2010
epl draft The role of interstitial gas in determining the impact response of granular beds arxiv:112.2399v1 [cond-mat.soft] 1 Dec 21 John R. Royer 1 (a), Bryan Conyers 1, Eric I. Corwin 1 (b), Peter J.
More informationVersion PREVIEW Semester 1 Review Slade (22222) 1
Version PREVIEW Semester 1 Review Slade () 1 This print-out should have 48 questions. Multiple-choice questions may continue on the next column or page find all choices before answering. Holt SF 0Rev 10A
More informationUNIT I: MECHANICS Chapter 5: Projectile Motion
IMPORTANT TERMS: Component Projectile Resolution Resultant Satellite Scalar quantity Vector Vector quantity UNIT I: MECHANICS Chapter 5: Projectile Motion I. Vector and Scalar Quantities (5-1) A. Vector
More informationPhysics 6010, Fall Relevant Sections in Text: Introduction
Physics 6010, Fall 2016 Introduction. Configuration space. Equations of Motion. Velocity Phase Space. Relevant Sections in Text: 1.1 1.4 Introduction This course principally deals with the variational
More informationAction Principles in Mechanics, and the Transition to Quantum Mechanics
Physics 5K Lecture 2 - Friday April 13, 2012 Action Principles in Mechanics, and the Transition to Quantum Mechanics Joel Primack Physics Department UCSC This lecture is about how the laws of classical
More informationBLUE PRINT FOR MODEL QUESTION PAPER 3 BLUE PRINT
Unit Chapter Number Number of teaching Hours Weightage of marks Mark Marks 3 Marks 5 Marks (Theory) 5 Marks (Numerica l Problem) BLUE PRINT FOR MODEL QUESTION PAPER 3 Class : I PUC Subject : PHYSICS (33)
More informationWork and Energy continued
Chapter 6 Work and Energy continued 6.2 The Work-Energy Theorem and Kinetic Energy Chapters 1 5 Motion equations were been developed, that relate the concepts of velocity, speed, displacement, time, and
More informationy(t) = y 0 t! 1 2 gt 2. With y(t final ) = 0, we can solve this for v 0 : v 0 A ĵ. With A! ĵ =!2 and A! = (2) 2 + (!
1. The angle between the vector! A = 3î! 2 ĵ! 5 ˆk and the positive y axis, in degrees, is closest to: A) 19 B) 71 C) 90 D) 109 E) 161 The dot product between the vector! A = 3î! 2 ĵ! 5 ˆk and the unit
More informationFirst Order ODEs (cont). Modeling with First Order ODEs
First Order ODEs (cont). Modeling with First Order ODEs September 11 15, 2017 Bernoulli s ODEs Yuliya Gorb Definition A first order ODE is called a Bernoulli s equation iff it is written in the form y
More informationCEE 271: Applied Mechanics II, Dynamics Lecture 17: Ch.15, Sec.2 4
1 / 38 CEE 271: Applied Mechanics II, Dynamics Lecture 17: Ch.15, Sec.2 4 Prof. Albert S. Kim Civil and Environmental Engineering, University of Hawaii at Manoa Tuesday, October 16, 2012 2 / 38 PRINCIPLE
More informationMultiple choice questions [60 points]
Multiple choice questions [60 points] Answer all of the following questions. Read each question carefully. Fill the correct bubble on your scantron sheet. Each correct answer is worth 4 points. Each question
More informationThe Effect of Air Resistance Harvey Gould Physics 127 2/06/07
1 The Effect of Air Resistance Harvey Gould Physics 127 2/06/07 I. INTRODUCTION I simulated the motion of a falling body near the earth s surface in the presence of gravity and air resistance. I used the
More informationAnalytical reconsideration of the von Neumann paradox in the reflection of a shock wave over a wedge
PHYSICS OF FLUIDS 20, 046101 2008 Analytical reconsideration of the von Neumann paradox in the reflection of a shock wave over a wedge Eugene I. Vasilev, 1 Tov Elperin, 2 and Gabi Ben-Dor 2 1 Department
More information3-Fold Decomposition EFB Closure for Convective Turbulence and Organized Structures
3-Fold Decomposition EFB Closure for Convective Turbulence and Organized Structures Igor ROGACHEVSKII and Nathan KLEEORIN Ben-Gurion University of the Negev, Beer-Sheva, Israel N.I. Lobachevsky State University
More informationCenter of Gravity. The location of the center of gravity is defined by: n mgx. APSC 111 Review Page 7
Center of Gravity We have said that for rigid bodies, all of the forces act at the centre of mass. This is a normally a very good approximation, but strictly speaking, the forces act at the centre of gravity,
More informationBasic Theory of Differential Equations
page 104 104 CHAPTER 1 First-Order Differential Equations 16. The following initial-value problem arises in the analysis of a cable suspended between two fixed points y = 1 a 1 + (y ) 2, y(0) = a, y (0)
More informationEmden-Fowler Equation and Inverse Analysis of Simple Flows through Variable Permeability Porous Layers
4 Int. J. Open Problems Compt. Math., Vol. 0, No., June 07 ISSN 998-66; Copyright ICSRS Publication, 07 www.i-csrs.org Emden-Fowler Equation and Inverse Analysis of Simple Flows through Variable Permeability
More informationEnergy balance in self-powered MR damper-based vibration reduction system
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES, Vol. 59, No. 1, 2011 DOI: 10.2478/v10175-011-0011-4 Varia Energy balance in self-powered MR damper-based vibration reduction system J. SNAMINA
More informationMechanical Engineering Research Journal PENETRATION AND PERFORATION OF COMPOSITE STRUCTURES
Dept. of Mech. Eng. CUET Published Online March 2015 (http://www.cuet.ac.bd/merj/index.html) Mechanical Engineering Research Journal Vol. 9, pp. 37 42, 2013 M E R J ISSN: 1990-5491 PENETRATION AND PERFORATION
More informationUniversity of Alabama Department of Physics and Astronomy. PH 105 LeClair Summer Problem Set 3 Solutions
University of Alabama Department of Physics and Astronomy PH 105 LeClair Summer 2012 Instructions: Problem Set 3 Solutions 1. Answer all questions below. All questions have equal weight. 2. Show your work
More informationEvaporation and Condensation of Large Droplets in the Presence of Inert Admixtures Containing Soluble Gas
MARCH 2007 E L P E R I N E T A L. 983 Evaporation and Condensation of Large Droplets in the Presence of Inert Admixtures Containing Soluble Gas T. ELPERIN, A. FOMINYKH, AND B. KRASOVITOV Department of
More informationPHYS 100 Mid-Term #1
D.W. Poppy Secondary School Physics 12 PHYS 100 Mid-Term #1 Name: Directions: Fill in the scantron form with the following information: 1. ID number (student number) 2. Name at top of form 3. Name bubbled
More informationDaba Meshesha Gusu and O.Chandra Sekhara Reddy 1
International Journal of Basic and Applied Sciences Vol. 4. No. 1 2015. Pp.22-27 Copyright by CRDEEP. All Rights Reserved. Full Length Research Paper Solutions of Non Linear Ordinary Differential Equations
More information0J2 - Mechanics Lecture Notes 2
0J2 - Mechanics Lecture Notes 2 Work, Power, Energy Work If a force is applied to a body, which then moves, we say the force does work. In 1D, if the force is constant with magnitude F, and the body moves
More informationNon-textbook problem #I: The kinetic energy of a body depends on its mass and speed as. K = 1 2 mv2 (1) m 1 v 2 1 = 1 2 m 2v 2 2 (2)
PHY 309 K. Solutions for Problem set # 7. Non-textbook problem #I: The kinetic energy of a body depends on its mass and speed as K = 1 mv (1) Therefore, two bodies of respective masses m 1 and m and speeds
More informationEnergy present in a variety of forms. Energy can be transformed form one form to another Energy is conserved (isolated system) ENERGY
ENERGY Energy present in a variety of forms Mechanical energy Chemical energy Nuclear energy Electromagnetic energy Energy can be transformed form one form to another Energy is conserved (isolated system)
More informationPRINCIPLE OF LINEAR IMPULSE AND MOMENTUM (Section 15.1)
PRINCIPLE OF LINEAR IMPULSE AND MOMENTUM (Section 15.1) Linear momentum: L = mv vector mv is called the linear momentum denoted as L (P in 1120) vector has the same direction as v. units of (kg m)/s or
More informationENGI 3424 First Order ODEs Page 1-01
ENGI 344 First Order ODEs Page 1-01 1. Ordinary Differential Equations Equations involving only one independent variable and one or more dependent variables, together with their derivatives with respect
More informationPhysics 20 Homework 2 SIMS 2016
Physics 20 Homework 2 SIMS 2016 Due: Saturday, August 20 th 1. In class, we ignored air resistance in our discussion of projectile motion. Now, let s derive the relevant equation of motion in the case
More informationChapter 3 Acceleration
Chapter 3 Acceleration Slide 3-1 Chapter 3: Acceleration Chapter Goal: To extend the description of motion in one dimension to include changes in velocity. This type of motion is called acceleration. Slide
More informationModule 17: Systems, Conservation of Momentum and Center of Mass
Module 17: Systems, Conservation of Momentum and Center of Mass 17.1 External and Internal Forces and the Change in Momentum of a System So far we have restricted ourselves to considering how the momentum
More information! 1.1 Definitions and Terminology
! 1.1 Definitions and Terminology 1. Introduction: At times, mathematics aims to describe a physical phenomenon (take the population of bacteria in a petri dish for example). We want to find a function
More informationChapter 4 Making Sense of the Universe: Understanding Motion, Energy, and Gravity. Copyright 2009 Pearson Education, Inc.
Chapter 4 Making Sense of the Universe: Understanding Motion, Energy, and Gravity How do we describe motion? Precise definitions to describe motion: Speed: Rate at which object moves speed = distance time
More informationwhere x 0 is arbitrary.
The forces internal to a system are of two types. Conservative forces, such as gravity; and dissipative forces such as friction. Internal forces arise from the natural dynamics of the system in contract
More informationVectors and Coordinate Systems
Vectors and Coordinate Systems In Newtonian mechanics, we want to understand how material bodies interact with each other and how this affects their motion through space. In order to be able to make quantitative
More information27 the electric field
27 the electric field With every point in space near the earth we can associate a gravitational field vector g (see Eq. 16-12). This is the gravitational acceleration that a test body, placed at that point
More informationMATH 320, WEEK 2: Slope Fields, Uniqueness of Solutions, Initial Value Problems, Separable Equations
MATH 320, WEEK 2: Slope Fields, Uniqueness of Solutions, Initial Value Problems, Separable Equations 1 Slope Fields We have seen what a differential equation is (relationship with a function and its derivatives),
More informationDynamics. Dynamics of mechanical particle and particle systems (many body systems)
Dynamics Dynamics of mechanical particle and particle systems (many body systems) Newton`s first law: If no net force acts on a body, it will move on a straight line at constant velocity or will stay at
More informationA Correlation of Pearson Physics 2014
A Correlation of Pearson Physics Ohio s New Learning Standards for Science, 2011 Physics, High School Science Inquiry and Application Course Content A Correlation of, INTRODUCTION This document demonstrates
More informationNumerical Analysis: Interpolation Part 1
Numerical Analysis: Interpolation Part 1 Computer Science, Ben-Gurion University (slides based mostly on Prof. Ben-Shahar s notes) 2018/2019, Fall Semester BGU CS Interpolation (ver. 1.00) AY 2018/2019,
More informationPhysics 0174(CHS) Exam #1 Academic Year NAME
. Physics 0174(CHS) Exam #1 Academic Year 2015-2016 NAME This exam consists of 6 pages in addition to this page; please check to see that you have all of them. Be sure to show clearly how you arrive at
More informationUNIVERSITY OF MANITOBA
PAGE NO.: 1 of 6 + Formula Sheet Equal marks for all questions. No marks are subtracted for wrong answers. Record all answers on the computer score sheet provided. USE PENCIL ONLY! Black pen will look
More informationDistance travelled time taken and if the particle is a distance s(t) along the x-axis, then its instantaneous speed is:
Chapter 1 Kinematics 1.1 Basic ideas r(t) is the position of a particle; r = r is the distance to the origin. If r = x i + y j + z k = (x, y, z), then r = r = x 2 + y 2 + z 2. v(t) is the velocity; v =
More informationSalmon: Lectures on partial differential equations
4 Burger s equation In Lecture 2 we remarked that if the coefficients in u x, y,! "! "x + v x,y,! "! "y = 0 depend not only on x,y but also on!, then the characteristics may cross and the solutions become
More information8.01T Problem Set 2 Solutions Fall 2004
Problem 1. Measurement of g. a) The ball moves in the vertical direction under the influence of the constant force of gravity 1. Hence in our approximation the ball undergoes onedimensional motion with
More informationUnit 1: Mechanical Equilibrium
Unit 1: Mechanical Equilibrium Chapter: Two Mechanical Equilibrium Big Idea / Key Concepts Student Outcomes 2.1: Force 2.2: Mechanical Equilibrium 2.3: Support Force 2.4: Equilibrium for Moving Objects
More informationImpacts of Kinematic Links with a Granular Material. Seunghun Lee
Impacts of Kinematic Links with a Granular Material Except where reference is made to the work of others, the work described in this dissertation is my own or was done in collaboration with my advisory
More informationPHYSICS 221, FALL 2010 EXAM #1 Solutions WEDNESDAY, SEPTEMBER 29, 2010
PHYSICS 1, FALL 010 EXAM 1 Solutions WEDNESDAY, SEPTEMBER 9, 010 Note: The unit vectors in the +x, +y, and +z directions of a right-handed Cartesian coordinate system are î, ĵ, and ˆk, respectively. In
More informationCombined Effect of Magnetic field and Internal Heat Generation on the Onset of Marangoni Convection
International Journal of Fluid Mechanics & Thermal Sciences 17; 3(4): 41-45 http://www.sciencepublishinggroup.com/j/ijfmts doi: 1.11648/j.ijfmts.1734.1 ISSN: 469-815 (Print); ISSN: 469-8113 (Online) ombined
More informationChapter 8. Potential Energy & Conservation of Energy
Chapter 8 Potential Energy & Conservation of Energy 8.1 Potential Energy Technically, potential energy is energy that can be associated with the configuration (arrangement) of a system of objects that
More information1.The pressure drop per unit length that develops due to friction cannot generally be solved analytically. A. True B. False
CHAPTER 07 1.The pressure drop per unit length that develops due to friction cannot generally be solved analytically. 2.A qualitative description of physical quantities can be given in terms of. YOUR ANSWER:
More informationAP Physics C Textbook Problems
AP Physics C Textbook Problems Chapter 13 Pages 412 416 HW-16: 03. A 200-kg object and a 500-kg object are separated by 0.400 m. Find the net gravitational force exerted by these objects on a 50.0-kg object
More informationStars Above, Earth Below By Tyler Nordgren Laboratory Exercise for Chapter 7
Name Section Partners By Tyler Nordgren Laboratory Exercise for Chapter 7 Equipment: Ruler Sand box Meter stick Log-log paper Small balls such as those included in the table at the end of the lab THE FORMATION
More informationPHYS 111 K SECOND HOUR EXAM 2015
PHYS 111 K SECOND HOUR EXAM 2015 This is a closed book closed note exam. Do all your writing in your blue book (s) and be sure to put your name on each blue book you use. You will not need nor are permitted
More information8 Rotational motion of solid objects
8 Rotational motion of solid objects Kinematics of rotations PHY166 Fall 005 In this Lecture we call solid objects such extended objects that are rigid (nondeformable) and thus retain their shape. In contrast
More informationProblem 1 Problem 2 Problem 3 Problem 4 Total
Name Section THE PENNSYLVANIA STATE UNIVERSITY Department of Engineering Science and Mechanics Engineering Mechanics 12 Final Exam May 5, 2003 8:00 9:50 am (110 minutes) Problem 1 Problem 2 Problem 3 Problem
More informationLecture D10 - Angular Impulse and Momentum
J. Peraire 6.07 Dynamics Fall 2004 Version.2 Lecture D0 - Angular Impulse and Momentum In addition to the equations of linear impulse and momentum considered in the previous lecture, there is a parallel
More informationEN40: Dynamics and Vibrations. Midterm Examination Thursday April
EN40: Dynamics and Vibrations Mierm Examination Thursday April 015 School of Engineering Brown University NAME: General Instructions No collaboration of any kind is permitted on this examination. You may
More informationIntroduction to Engineering Mechanics
Introduction to Engineering Mechanics Statics October 2009 () Introduction 10/09 1 / 19 Engineering mechanics Engineering mechanics is the physical science that deals with the behavior of bodies under
More information16. Friction Mechanics of Manipulation
16. Friction Mechanics of Manipulation Matt Mason matt.mason@cs.cmu.edu http://www.cs.cmu.edu/~mason Carnegie Mellon Lecture 16. Mechanics of Manipulation p.1 Lecture 16. Friction. Chapter 1 Manipulation
More informationTALLINN UNIVERSITY OF TECHNOLOGY, DIVISION OF PHYSICS 13. STOKES METHOD
13. STOKES METHOD 1. Objective To determine the coefficient of viscosity of a known fluid using Stokes method.. Equipment needed A glass vessel with glycerine, micrometer calliper, stopwatch, ruler. 3.
More informationMATHEMATICAL MODELLING, MECHANICS AND MOD- ELLING MTHA4004Y
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 2017 18 MATHEMATICAL MODELLING, MECHANICS AND MOD- ELLING MTHA4004Y Time allowed: 2 Hours Attempt QUESTIONS 1 and 2, and ONE other
More informationENGI 2422 First Order ODEs - Separable Page 3-01
ENGI 4 First Order ODEs - Separable Page 3-0 3. Ordinary Differential Equations Equations involving only one independent variable and one or more dependent variables, together with their derivatives with
More informationBoundary Conditions for the Child Langmuir Sheath Model
IEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. 28, NO. 6, DECEMBER 2000 2207 Boundary Conditions for the Child Langmuir Sheath Model Mikhail S. Benilov Abstract A collision-free space-charge sheath formed by
More informationUniversity Physics 226N/231N Old Dominion University. Exam Review, Friction and Work
University Physics 226N/231N Old Dominion University Exam Review, Friction and Work Dr. Todd Satogata (ODU/Jefferson Lab) satogata@jlab.org http://www.toddsatogata.net/2016-odu Wednesday, September 28,
More informationarxiv: v1 [cond-mat.stat-mech] 7 Mar 2019
Langevin thermostat for robust configurational and kinetic sampling Oded Farago, Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB EW, United Kingdom Department of Biomedical
More informationStudy on Estimation of Hydraulic Conductivity of Porous Media Using Drag Force Model Jashandeep Kaur, M. A. Alam
26 IJSRSET Volume 2 Issue 3 Print ISSN : 2395-99 Online ISSN : 2394-499 Themed Section: Engineering and Technology Study on Estimation of Hydraulic Conductivity of Porous Media Using Drag Force Model Jashandeep
More informationPhysics 142 Energy in Mechanics Page 1. Energy in Mechanics
Physics 4 Energy in Mechanics Page Energy in Mechanics This set of notes contains a brief review of the laws and theorems of Newtonian mechanics, and a longer section on energy, because of its central
More informationTurbulentlike Quantitative Analysis on Energy Dissipation in Vibrated Granular Media
Copyright 011 Tech Science Press CMES, vol.71, no., pp.149-155, 011 Turbulentlike Quantitative Analysis on Energy Dissipation in Vibrated Granular Media Zhi Yuan Cui 1, Jiu Hui Wu 1 and Di Chen Li 1 Abstract:
More informationProjectile Motion. Chin- Sung Lin STEM GARAGE SCIENCE PHYSICS
Projectile Motion Chin- Sung Lin Introduction to Projectile Motion q What is Projectile Motion? q Trajectory of a Projectile q Calculation of Projectile Motion Introduction to Projectile Motion q What
More informationFlow and Natural Convection Heat Transfer in a Power Law Fluid Past a Vertical Plate with Heat Generation
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.7(2009) No.1,pp.50-56 Flow and Natural Convection Heat Transfer in a Power Law Fluid Past a Vertical Plate with
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department 8.01 Physics I Fall Term 2009 Review Module on Solving N equations in N unknowns
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department 8.01 Physics I Fall Term 2009 Review Module on Solving N equations in N unknowns Most students first exposure to solving N linear equations in N
More informationChapter 12 Study Guide
Chapter 12 Study Guide Key Concepts 12.1 12.2 12.3 12.4 A force can cause a resting object to move, or it can accelerate a moving object by changing the object s speed or direction. When the forces on
More informationSection 8.0 Introduction to Boundary Value Problems. How do initial value problems (IVPs) and boundary value problems (BVPs) differ?
Section 8.0 Introduction to Boundary Value Problems Key terms/ideas How do initial value problems (IVPs) and boundary value problems (BVPs) differ? What are boundary conditions? In what type of problems
More information(c) McHenry Software. McHenry. Accident Reconstruction. by Raymond R. McHenry Brian G. McHenry. McHenry Training Seminar 2008
McHenry Training eminar 008 McHenry Accident Reconstruction 008 by Raymond R. McHenry Brian G. McHenry McHenry oftware PO Box 1716 Cary, NC 751 UA (919)-468-966 email: mchenry@mchenrysoftware.com www:
More informationAnalytical engineering models of high speed normal impact by hard projectiles on metal shields. Author's copy
Cent. Eur. J. Eng. 3(3) 2013 349-373 DOI: 10.2478/s13531-013-0108-7 Central European Journal of Engineering Analytical engineering models of high speed normal impact by hard projectiles on metal shields
More informationThe effect of disturbances on the flows under a sluice gate and past an inclined plate
J. Fluid Mech. (7), vol. 576, pp. 475 49. c 7 Cambridge University Press doi:.7/s7486 Printed in the United Kingdom 475 The effect of disturbances on the flows under a sluice gate and past an inclined
More informationRigid bodies - general theory
Rigid bodies - general theory Kinetic Energy: based on FW-26 Consider a system on N particles with all their relative separations fixed: it has 3 translational and 3 rotational degrees of freedom. Motion
More informationResistance is Futile
Resistance is Futile Joseph Hays 16 May 2014 Page 1 of 24 Abstract When introductory physics books consider projectile motion, they often suggest that the reader assume air resistance is negligible; an
More informationPHYSICS 231 INTRODUCTORY PHYSICS I
PHYSICS 231 INTRODUCTORY PHYSICS I Lecture 4 Main points of last lecture Scalars vs. Vectors Vectors A: (A x, A y ) or A & θ Addition/Subtraction Projectile Motion X-direction: a x = 0 (v x = constant)
More informationInternational ejournals
Available online at www.internationalejournals.com ISSN 0976 1411 International ejournals International ejournal of Mathematics and Engineering 30 (013) 4 55 RADIATION EFFECTS ON UNSTEADY FLOW PAST AN
More information