![Numerical Methods for Engineers](https://www.bartleby.com/isbn_cover_images/9780073397924/9780073397924_largeCoverImage.gif)
The dynamics of a forced spring-mass-damper system can be represented by the following second-order ODE:
where
(a)
![Check Mark](/static/check-mark.png)
To calculate: The displacement and velocity as a function of time for linear system where
Answer to Problem 49P
Solution:
The first few solutions for displacement and velocity as a function of time for liner system is,
t | x | v |
0 | 0 | 0 |
0.125 | 0.003836 | 0.060764 |
0.25 | 0.015019 | 0.117402 |
0.375 | 0.032976 | 0.169015 |
0.5 | 0.05703 | 0.214828 |
0.625 | 0.086414 | 0.254195 |
0.75 | 0.120289 | 0.286602 |
0.875 | 0.157759 | 0.311673 |
1 | 0.197891 | 0.329164 |
1.125 | 0.239729 | 0.338967 |
1.25 | 0.282313 | 0.341104 |
1.375 | 0.324691 | 0.335717 |
1.5 | 0.365939 | 0.323069 |
1.625 | 0.405171 | 0.303527 |
1.75 | 0.441553 | 0.277555 |
1.875 | 0.474314 | 0.245705 |
2 | 0.50276 | 0.208601 |
2.125 | 0.526274 | 0.16693 |
2.25 | 0.544332 | 0.121428 |
2.375 | 0.556503 | 0.072863 |
2.5 | 0.562454 | 0.02203 |
2.625 | 0.56195 | -0.03027 |
2.75 | 0.554859 | -0.08324 |
2.875 | 0.541146 | -0.13608 |
3 | 0.520874 | -0.18805 |
3.125 | 0.494199 | -0.23842 |
3.25 | 0.461364 | -0.28651 |
3.375 | 0.422694 | -0.33168 |
3.5 | 0.378588 | -0.37339 |
3.625 | 0.329513 | -0.41112 |
3.75 | 0.275993 | -0.44444 |
3.875 | 0.218601 | -0.47301 |
4 | 0.157952 | -0.49653 |
4.125 | 0.094688 | -0.5148 |
4.25 | 0.029475 | -0.52771 |
4.375 | -0.03701 | -0.53518 |
4.5 | -0.1041 | -0.53725 |
4.625 | -0.1711 | -0.53399 |
4.75 | -0.23738 | -0.52558 |
4.875 | -0.30229 | -0.51221 |
5 | -0.36524 | -0.49416 |
Explanation of Solution
Given Information:
The dynamic of a forced spring-mass-damper system is given as,
The values,
The initial condition,
Formula used:
The fourth-order RK method for
Where,
Calculation:
Consider the dynamic of a forced spring-mass-damper system,
As
Divide both the sides of above equation by m,
Now, substitute the values
For linear, substitute
Use VBA code for RK4 method as below to solve for x and v,
Code:
Output:
To draw the graph of the above results, follow the steps in excel sheet as given below,
Step 1: Select the cell from A4 to A205 and cell B4 to B205. Then, go to the Insert and select the scatter with smooth lines from the chart.
Step 2: Select the cell from A4 to A205 and cell C4 to C205. Then, go to the Insert and select the scatter with smooth lines from the chart.
Step 3: Select one of the graphs and paste it on another graph to merge the graphs.
The graph obtained is,
And, to draw the phase plane plot follow the steps as below,
Step 4: Select the column B and column C. Then, go to the Insert and select the scatter with smooth lines from the chart.
The graph obtained is,
(b)
![Check Mark](/static/check-mark.png)
To calculate: The displacement and velocity as a function of time for non-linear system where
Answer to Problem 49P
Solution:
t | x | v |
0 | 0 | 0 |
0.125 | 0.003836 | 0.060764 |
0.25 | 0.015019 | 0.117402 |
0.375 | 0.032976 | 0.169014 |
0.5 | 0.05703 | 0.214821 |
0.625 | 0.086412 | 0.254166 |
0.75 | 0.120279 | 0.286506 |
0.875 | 0.157729 | 0.311417 |
1 | 0.19781 | 0.328581 |
1.125 | 0.239542 | 0.337784 |
1.25 | 0.28192 | 0.338923 |
1.375 | 0.323936 | 0.332004 |
1.5 | 0.36459 | 0.31716 |
1.625 | 0.402907 | 0.294658 |
1.75 | 0.437953 | 0.264919 |
1.875 | 0.468859 | 0.228524 |
2 | 0.494837 | 0.186221 |
2.125 | 0.515205 | 0.138911 |
2.25 | 0.5294 | 0.087637 |
2.375 | 0.536996 | 0.033543 |
2.5 | 0.537718 | -0.02217 |
2.625 | 0.531437 | -0.07829 |
2.75 | 0.518177 | -0.13366 |
2.875 | 0.498099 | -0.18721 |
3 | 0.47149 | -0.23801 |
3.125 | 0.438743 | -0.28531 |
3.25 | 0.400333 | -0.32852 |
3.375 | 0.3568 | -0.36723 |
3.5 | 0.308725 | -0.40117 |
3.625 | 0.256712 | -0.43023 |
3.75 | 0.201374 | -0.45436 |
3.875 | 0.143325 | -0.47361 |
4 | 0.083172 | -0.48804 |
4.125 | 0.021513 | -0.49771 |
4.25 | -0.04106 | -0.50268 |
4.375 | -0.10396 | -0.50296 |
4.5 | -0.1666 | -0.49854 |
4.625 | -0.2284 | -0.48938 |
4.75 | -0.28875 | -0.47544 |
4.875 | -0.34705 | -0.45667 |
5 | -0.40271 | -0.43306 |
Explanation of Solution
Given Information:
The dynamic of a forced spring-mass-damper system is given as,
The values,
And,
The initial condition,
Formula used:
The fourth-order RK method for
Where,
Calculation:
Consider the dynamic of a forced spring-mass-damper system,
As
Divide both the sides of above equation by m,
Now, substitute the values
For linear, substitute
Use VBA code for RK4 method as below to solve for x and v,
Code:
Output:
Few data are shown below,
To draw the graph of the above results, follow the steps in excel sheet as given below,
Step 1: Select the cell from A4 to A205 and cell B4 to B205. Then, go to the Insert and select the scatter with smooth lines from the chart.
Step 2: Select the cell from A4 to A205 and cell C4 to C205. Then, go to the Insert and select the scatter with smooth lines from the chart.
Step 3: Select one of the graphs and paste it on another graph to merge the graphs.
The graph obtained is,
And, to draw the phase plane plot follow the steps as below,
Step 4: Select the column B and column C. Then, go to the Insert and select the scatter with smooth lines from the chart.
The graph obtained is,
Want to see more full solutions like this?
Chapter 28 Solutions
Numerical Methods for Engineers
- A press that delivers 115 strokes per minute, each stroke providing a force of 7826 N throughout a distance of 18 mm. The press efficiency is 90% and is driven by a 1749-rpm motor. Determine average torque that must be provided by the motor in the units of N-m.arrow_forward·3) find the force (P) for the figures (1) and (2) 15cm 10cm 15 h=10mm h2=6mm // Call = 90 N/2 P Agate Fig (i) Ans: 1)P=112614N 2) P=1956.5 N 25cm 25 cm الفترة أو الحجم تمر بالتي عثر اكو تورشن (ک Fig (2) h₁ = 10mm 42=6mm Cmarrow_forwardI want a human solutionarrow_forward
- (Read Image)arrow_forward47 14 16 12 34 10 12 12 33arrow_forward3. A steam power plant has an average monthly net power delivery of 740 MW over the course of a year. This power delivery is accomplished by burning coal in the boiler. The coal has a heating value of 9150 Btu/lbm. The cost of the coal is $14.20/ton. The overall thermal efficiency of the plant is, nth = Wnet Qboiler = 0.26 = 26% Determine the annual cost of the coal required to deliver the given average monthly power.arrow_forward
- 47 14 16 12 34 10 12 12 33arrow_forward= The forces F₁ = 590 lb, F₂ = 380 lb, F3 = 240 lb and F 330 lb. Determine the forces in each member of the truss. Use positive values to indicate tension and negative values to indicate compression. a a a D b F₁ A 000 B. 779977 F₂V H G E F4 b BY NC SA 2013 Michael Swanbom Values for dimensions on the figure are given in the following table. Note the figure may not be to scale. Variable Value a 6 ft b 10.1 ft The force in member AB is lb. The force in member AH is lb. The force in member GH is lb. The force in member BH is lb. The force in member BC is lb. The force in member BG is lb. The force in member EG is lb. The force in member CD is lb. The force in member DE is lb. The force in member CE is lb. The force in member CG is lb.arrow_forwardMultiple Choice Circle the best answer to each statement. 1. Which type of surface deviation is controlled by a cy- lindricity tolerance but not by a circularity tolerance? A. B. C. Ovality Taper Lobing D. None of the above 2. When verifying a cylindricity tolerance, the inspec- tion method must be able to collect a set of points and determine the: A. Distance between two coaxial cylinders that con- tain the set of points B. Cylinder that circumscribes the set of points C. Cylinder that inscribes the set of points D. Distance between two coaxial circles that contain the set of points 3. Where Rule #1 applies to a cylindrical regular feature of size, the tolerance value of a cylindricity tolerance applied to the feature of size must be tolerance. A. Less than B. Equal to C. Greater than D. None of the above the size 4. Which of the following modifiers may be applied with a cylindricity tolerance? A. M B. C. ℗ D. Ø 5. Which geometric tolerance can provide an indirect cylindricity…arrow_forward
- The beam AB is attached to the wall in the xz plane by a fixed support at A. A force of F = (−129î + 69.0ĵ + 3591) N is applied to the end of the beam at B. The weight of the beam can be modeled with a uniform distributed load of intensity w = 85.0 N/m acting in the negative z direction along its entire length. Find the support reactions at A. Z с A b a B F y Cc 10 BY NC SA 2016 Eric Davishahl X Values for dimensions on the figure are given in the following. table. Note the figure may not be to scale. Variable Value a 5.60 m b 5.00 m C 3.70 m A II = MA = ( m 2.> ~.> + + k) N k) N-arrow_forwardneed help?arrow_forwardA bent pipe is attached to a wall with brackets as shown. A force of F = 180 lb is applied to the end of the tube with direction indicated by the dimensions in the figure. Determine the support reactions at the brackets B, C, and D. Model these brackets as journal bearings (only force reactions perpendicular to the axis of the tube) and neglect couple moment reactions. Assume the distance between the supports at B and C and the tube bends nearby are negligible such that the support at C is directly above the support at D and the dimension g gives the distance between supports B and C. Enter your answers in Cartesian components. 2013 Michael Swanbom cc 10 BY NC SA g h א B 8° У A C x каж Values for dimensions on the figure are given in the table below. Note the figure may not be to scale. Variable Value a 6.72 in b 11.8 in с 14.8 in d 42.0 in h 26.6 in g 28.0 in → The reaction at B is B = lb. The reaction at C is C = lb. The reaction at D is D = lb. + << + + 2. + + 557 〈んarrow_forward
- Elements Of ElectromagneticsMechanical EngineeringISBN:9780190698614Author:Sadiku, Matthew N. O.Publisher:Oxford University PressMechanics of Materials (10th Edition)Mechanical EngineeringISBN:9780134319650Author:Russell C. HibbelerPublisher:PEARSONThermodynamics: An Engineering ApproachMechanical EngineeringISBN:9781259822674Author:Yunus A. Cengel Dr., Michael A. BolesPublisher:McGraw-Hill Education
- Control Systems EngineeringMechanical EngineeringISBN:9781118170519Author:Norman S. NisePublisher:WILEYMechanics of Materials (MindTap Course List)Mechanical EngineeringISBN:9781337093347Author:Barry J. Goodno, James M. GerePublisher:Cengage LearningEngineering Mechanics: StaticsMechanical EngineeringISBN:9781118807330Author:James L. Meriam, L. G. Kraige, J. N. BoltonPublisher:WILEY
![Text book image](https://www.bartleby.com/isbn_cover_images/9780190698614/9780190698614_smallCoverImage.gif)
![Text book image](https://www.bartleby.com/isbn_cover_images/9780134319650/9780134319650_smallCoverImage.gif)
![Text book image](https://www.bartleby.com/isbn_cover_images/9781259822674/9781259822674_smallCoverImage.gif)
![Text book image](https://www.bartleby.com/isbn_cover_images/9781118170519/9781118170519_smallCoverImage.gif)
![Text book image](https://www.bartleby.com/isbn_cover_images/9781337093347/9781337093347_smallCoverImage.gif)
![Text book image](https://www.bartleby.com/isbn_cover_images/9781118807330/9781118807330_smallCoverImage.gif)