Repeat Example 30.3, but for
To calculate: The temperature distribution of a long, thin aluminum rod with the given properties (refer Example 30.3, in the text book) when the rod is initially at
Answer to Problem 5P
Solution: The temperature distribution is,
Explanation of Solution
Given Information:
Consider, time is represented by variable t.
Length of the rod is 10 cm.
Here,
At,
Constant
Formula used:
At boundary node
From Crank Nicolson method, the difference equation for the nodes (except first and the last node) is expressed as,
Here,
Calculation:
For the first step at
Difference equation at node
Substitute
Since,
Therefore, the difference equations for node
Further, the system of equations (2), (3),(4), (5), (6), (7), (8), and (9) can be written in matrix form as,
Thus,
Execute the following code in MATLAB to evaluate the results in the above matrix equation
The output values thus obtained are:
Therefore,
Furthermore, for the second step to get the values at
Thus, solve as,
Execute the following code in MATLAB to evaluate the results in the above matrix equation
The output values thus obtained are:
Therefore,
Want to see more full solutions like this?
Chapter 30 Solutions
EBK NUMERICAL METHODS FOR ENGINEERS
Additional Engineering Textbook Solutions
Basic Technical Mathematics
Advanced Engineering Mathematics
Fundamentals of Differential Equations (9th Edition)
Calculus and Its Applications (11th Edition)
- Please solve this numerical problem.arrow_forwardSuppose a aluminum rod of 20 ft and solid square cross section of 0.5 inches, isfixed at a wall on the left end. On the right end it is unsupported. The rod will bendbased on it’s own weight given that it has uniform density.a.) Find the Youngs modulus, Moment of Inertia, and Weight per unit necessary tosolve this problem. Then put them together to find the coefficient on the x4term.b.) Solve for coefficients c1 through c4 and report the equation governing the shapethis rod will take.c.) How far down from level does the end of the aluminum rod hang in inches?arrow_forwardFor the figure shown below (not to scale)... (a) Write equations for the internal forces at J (axial (F), shear (V), and bending (M)) in terms of only P and a. (b) Graph your equations as functions of P, using a =45°, from 0arrow_forwardCalculate the stiffness matrix for the Q8 element shown below. Make use of the shape functions in natural coordinates and the Jacobian calculated in previous problems. Use 2x2 Gaussian numerical integration. Use plane strain conditions (for unit thickness) with E=200E9 (MPa) and nu=0.25. I have attached my matlab code below. Please look thru it and correct it. Thanks. u=sym('u') %%ztav=sym('v') %%eta x1=1x2=2x3=3.5x4=3.6x5=4.0x6=2.3x7=1.2x8=0.8 y1=1y2=1y3=1.3y4=2.3y5=3.3y6=3.2y7=3.0y8=2.0X=[x1;x2;x3;x4;x5;x6;x7;x8]Y=[y1;y2;y3;y4;y5;y6;y7;y8]%% diff of ztaDN1Du=1/4*(1-v)*(2*u+v);DN2Du=-u*(1-v);DN3Du=1/4*(1-v)*(2*u-v);DN4Du=1/2*(1-v^2);DN5Du=1/4*(1+v)*(2*u+v);DN6Du=-u*(1+v);DN7Du=1/4*(1+v)*(2*u-v);DN8Du=-1/2*(1-v^2); %% diff of etaDN1Dv=1/4*(1-u)*(2*v+u);DN2Dv=-1/2*(1-u^2);DN3Dv=1/4*(1+u)*(2*v-u);DN4Dv=-v*(1+u);DN5Dv=1/4*(1+u)*(2*v+u);DN6Dv=1/2*(1-u^2);DN7Dv=1/4*(1-u)*(2*v+u);DN8Dv=-v*(1-u);%% find dx/du dx/dv, dy/du and dy/dva=[DN1Du, DN2Du, DN3Du, DN4Du, DN5Du, DN6Du,DN7Du,…arrow_forwardPlease solve and include units in each step. Your work and solution will be appraciated much. Thank You!arrow_forwardA rectangular tube of outer width w = 7.5 m, outer height h = 5 m and thickness t = 0.17 m experiences bending through the x axis. What is the moment of inertia of the cross-sectional area In? y -- X MATLAB input variables: format shortEng W = 7.5; h = 5; t = 0.17; moment of inertia I, = m4arrow_forwardPlease help with the following problem. Thank you!arrow_forwardIn the Fig. 2 below, let Ki = K2 = K and ti = t=t. %3D T -T X Fig. 2 (a) Let T= 0 °C and T= 200 °C. Solve for T: and unknown rates of heat flow in term of k and t. MEC_AMO_TEM_035_02 Page 2 of 11 Finite Element Analysis (MECH 0016.1) – Spring - 2021 -Assignment 2-QP (b) Let T- 400 °C and let fs have the prescribed value f. What are the unknowns? Solve for them in term of K, t, and f.arrow_forwardCalculate the stiffness matrix for the Q8 element shown below. Make use of the shape functions in natural coordinates and the Jacobian calculated in previous problems. Use 2x2 Gaussian numerical integration. Use plane strain conditions (for unit thickness) with E=200E9 (MPa) and nu=0.25. The final answers is F1= 0.18452e12. I have attached my matlab code below. Please look thru it and correct it. Thanks. u=sym('u') %%ztav=sym('v') %%eta x1=1x2=2x3=3.5x4=3.6x5=4.0x6=2.3x7=1.2x8=0.8 y1=1y2=1y3=1.3y4=2.3y5=3.3y6=3.2y7=3.0y8=2.0X=[x1;x2;x3;x4;x5;x6;x7;x8]Y=[y1;y2;y3;y4;y5;y6;y7;y8]%% diff of ztaDN1Du=1/4*(1-v)*(2*u+v);DN2Du=-u*(1-v);DN3Du=1/4*(1-v)*(2*u-v);DN4Du=1/2*(1-v^2);DN5Du=1/4*(1+v)*(2*u+v);DN6Du=-u*(1+v);DN7Du=1/4*(1+v)*(2*u-v);DN8Du=-1/2*(1-v^2); %% diff of etaDN1Dv=1/4*(1-u)*(2*v+u);DN2Dv=-1/2*(1-u^2);DN3Dv=1/4*(1+u)*(2*v-u);DN4Dv=-v*(1+u);DN5Dv=1/4*(1+u)*(2*v+u);DN6Dv=1/2*(1-u^2);DN7Dv=1/4*(1-u)*(2*v+u);DN8Dv=-v*(1-u);%% find dx/du dx/dv, dy/du and dy/dva=[DN1Du, DN2Du, DN3Du,…arrow_forwardConsider the following linear equations,arrow_forwardThe Philippine ten-peso coin (P10) is the second largest denomination coin of the Philippine peso. The current version, issued since 2018, has a diameter of d, = 27.0000 mm and a thickness of họ = 2.0500 mm. When the coin is subjected to a temperature change of 100 C°, the diameter of the coin increases by 0.11 percent. Suppose the coefficient of linear expansion of the coin is unknown, which of the following expressions should you use to compute for it? *arrow_forwardThe natural frequency (f) of a piezoelectric cantilever beam is described by the following equation:, where E = Modulus of elasticity, I = Area moment of inertia, L = Length of the cantilever beam, and m = Mass of proof mass (or end mass). The area moment of inertia for an arbitrary shape R with respect to an arbitrary axis (for example, x-axis) is related by: (a) What is the dimension of an area moment of inertia? (b) If the thickness of a piezoelectric cantilever beam increases, does the natural frequency increase or decrease? (c) Both the thickness and width of a piezoelectric cantilever beam can change the natural frequency. Which parameter (thickness or width) will lead to the natural frequency change more significantly? Briefly explain your answer.arrow_forwardarrow_back_iosSEE MORE QUESTIONSarrow_forward_ios
- 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