{VERSION 6 1 "Windows NT (unknown)" "6.1" } {USTYLETAB {PSTYLE "Left Justified Maple Output" -1 12 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }1 3 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Warning" -1 7 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 0 0 0 2 2 1 0 0 0 1 }1 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "D ash Item" -1 16 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }1 1 0 0 3 3 2 0 2 0 2 2 -1 3 }{PSTYLE "Heading 4" -1 20 1 {CSTYLE "" -1 -1 "Times" 1 10 0 0 0 1 1 1 2 2 2 2 0 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Heading 3" -1 5 1 {CSTYLE "" -1 -1 "Times " 1 12 0 0 0 1 1 1 2 2 2 2 0 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 } {PSTYLE "Error" -1 8 1 {CSTYLE "" -1 -1 "Courier" 1 10 255 0 255 1 0 0 0 2 2 1 0 0 0 1 }1 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "Author" -1 19 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }3 1 0 0 8 8 2 0 2 0 2 2 -1 1 }{PSTYLE "Heading 2" -1 4 1 {CSTYLE "" -1 -1 "Times" 1 14 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }1 1 0 0 8 2 2 0 2 0 2 2 -1 1 }{PSTYLE "Help" -1 10 1 {CSTYLE "" -1 -1 "Courier" 1 9 0 0 255 1 0 0 0 2 2 1 0 0 0 1 }1 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "Text Ou tput" -1 2 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 0 0 0 2 2 1 0 0 0 1 }1 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "Heading 1" -1 3 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }1 1 0 0 8 4 2 0 2 0 2 2 -1 1 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }0 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "M aple Plot" -1 13 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }3 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Line Printed Output" -1 6 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 0 0 0 2 2 1 0 0 0 1 } 1 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "Title" -1 18 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 2 1 1 2 2 2 0 0 0 1 }3 1 0 0 12 12 2 0 2 0 2 2 -1 1 }{PSTYLE "Diagnostic" -1 9 1 {CSTYLE "" -1 -1 "Courier" 1 10 64 128 64 1 0 0 0 2 2 1 0 0 0 1 }1 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 } {PSTYLE "Maple Output" -1 11 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }3 3 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "List Item" -1 14 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }1 1 0 0 3 3 2 0 2 0 2 2 -1 5 }{PSTYLE "Bullet Item" -1 15 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }1 1 0 0 3 3 2 0 2 0 2 2 -1 2 }{PSTYLE "Fixed Width" -1 17 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }3 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{CSTYLE "La TeX" -1 32 "" 0 1 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }{CSTYLE "Maple Comment" -1 21 "Courier" 0 1 0 0 0 1 0 0 0 2 2 2 0 0 0 1 }{CSTYLE "2D Math Bol d" -1 5 "Times" 0 1 0 0 0 0 0 1 2 2 2 2 0 0 0 1 }{CSTYLE "Help Underli ned" -1 44 "Times" 1 12 0 0 0 0 0 0 1 2 2 2 0 0 0 1 }{CSTYLE "Output L abels" -1 29 "Times" 1 8 0 0 0 1 0 0 0 2 2 2 0 0 0 1 }{CSTYLE "Maple I nput" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 2 1 2 0 0 0 1 }{CSTYLE "Page \+ Number" -1 33 "Times" 0 1 0 0 0 0 0 0 2 2 2 2 0 0 0 1 }{CSTYLE "Maple \+ Input Placeholder" -1 200 "Courier" 1 12 200 0 200 1 0 1 0 2 1 2 0 0 0 1 }{CSTYLE "Text" -1 201 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 } {CSTYLE "2D Math Italic Small" -1 202 "Times" 0 1 0 0 0 0 1 0 2 2 2 2 0 0 0 1 }{CSTYLE "Help Heading" -1 26 "" 1 14 0 0 0 0 0 1 0 2 2 2 0 0 0 1 }{CSTYLE "Help Italic Bold" -1 40 "Times" 1 12 0 0 0 0 1 1 0 2 2 2 0 0 0 1 }{CSTYLE "Help Normal" -1 30 "Times" 1 12 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }{CSTYLE "Help Nonterminal" -1 24 "Courier" 0 1 0 0 0 1 0 1 0 2 2 2 0 0 0 1 }{CSTYLE "2D Math Small" -1 7 "Times" 0 1 0 0 0 0 0 0 2 2 2 2 0 0 0 1 }{CSTYLE "Help Underlined Italic" -1 43 "Times" 1 12 0 0 0 0 1 0 1 2 2 2 0 0 0 1 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 2 2 2 0 0 0 1 }{CSTYLE "Copyright" -1 34 "Times" 1 10 0 0 0 0 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Hyperlink" -1 17 "" 0 1 0 128 128 1 2 0 1 2 2 2 0 0 0 1 }{CSTYLE "Plot Title" -1 27 "" 1 10 0 0 0 0 0 1 0 2 2 2 0 0 0 1 }{CSTYLE "2D Input" -1 19 "Times" 0 1 255 0 0 1 0 0 2 2 1 2 0 0 0 1 }{CSTYLE "Help Maple Name" -1 35 "" 0 1 104 64 92 1 0 1 0 2 2 2 0 0 0 1 }{CSTYLE "Help Italic" -1 42 "Times" 1 12 0 0 0 0 1 0 0 2 2 2 0 0 0 1 }{CSTYLE "Help Underlined Bold" -1 41 "Times" 1 12 0 0 0 0 0 1 1 2 2 2 0 0 0 1 }{CSTYLE "Default" -1 38 "" 0 1 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }{CSTYLE "Help Fixed" -1 23 "Courier" 1 10 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }{CSTYLE "Dictionary Hyperlink" -1 45 "" 0 1 147 0 15 1 2 0 1 2 2 2 0 0 0 1 }{CSTYLE "2D Output" -1 20 "Times" 0 1 0 0 255 1 0 0 2 2 2 2 0 0 0 1 }{CSTYLE "2D Math Italic" -1 3 "Times" 0 1 0 0 0 0 1 0 2 2 2 2 0 0 0 1 }{CSTYLE "Popup" -1 31 "" 0 1 0 128 128 1 1 0 1 2 2 2 0 0 0 1 }{CSTYLE "2D Math Symbol 2" -1 16 "Times" 0 1 0 0 0 0 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Plot Text" -1 28 "" 1 8 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }{CSTYLE "Help Notes" -1 37 "" 0 1 0 0 0 1 0 1 0 2 2 2 0 0 0 1 }{CSTYLE "Help Bold" -1 39 "Times" 1 12 0 0 0 0 0 1 0 2 2 2 0 0 0 1 }{CSTYLE "Help Menus" -1 36 "" 0 1 0 0 0 1 0 1 0 2 2 2 0 0 0 1 } {CSTYLE "2D Comment" -1 18 "Times" 0 1 0 0 0 0 0 0 2 2 2 2 0 0 0 1 } {CSTYLE "Help Variable" -1 25 "Courier" 0 1 0 0 0 1 2 2 0 2 2 2 0 0 0 1 }{CSTYLE "2D Math Bold Small" -1 10 "Times" 0 1 0 0 0 0 0 1 2 2 2 2 0 0 0 1 }{CSTYLE "Help Emphasized" -1 203 "" 0 1 0 0 0 0 1 2 0 2 2 2 0 0 0 1 }{CSTYLE "Prompt" -1 1 "Courier" 0 1 0 0 0 1 0 0 0 2 2 2 0 0 0 1 }{PSTYLE "_pstyle1" -1 200 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }0 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{CSTYLE "_cstyle1" -1 204 "Times" 0 1 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }{PSTYLE "_pstyle2" -1 201 1 {CSTYLE "" -1 -1 "Courier" 0 1 255 0 0 1 0 1 0 2 1 2 0 0 0 1 }0 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "_pstyle3" -1 202 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }0 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }} {SECT 0 {EXCHG {PARA 200 "" 0 "" {TEXT 204 9 "Math 2250" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 31 "Maple Lab 8: Earthquake project" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 11 "Spring 2008" }} {PARA 200 "" 0 "" {TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 66 "Nam e _____________________________________ Class Time __________" } {TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 61 "Project 8. Solve problems L8.1 to L8.5. The problem hea ders:" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 57 "_______ PROBLEM L8.1. EARHQUAKE MODEL FOR A \+ BUILDING." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 67 "_______ \+ PROBLEM L8.2. TABLE OF NATURAL FREQUENCIES AND PERIODS." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 68 "_______ PROBLEM L8.3. U NDETERMINED COEFFICIENTS STEADY-STATE SOL" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 46 "_______ PROBLEM L8.4. PRACTICAL RESONANC E." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 44 "_______ PROB LEM L8.5. EARTHQUAKE DAMAGE." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 16 "SIX FLOOR Model." } {TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 63 "Refer to the textbook of Edwards-Penney, section 7.4, page 437." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 63 "Consider a building with six floors each wei ghing 50 tons. Each" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 64 "floor corresponds to a restoring Hooke's force with constant k=5" } {TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 64 "tons/foot. Assume tha t ground vibrations from the earthquake are" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 45 "modeled by (1/4)cos(wt) with period T=2*Pi/w ." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 200 "" 0 "" {TEXT 204 47 "PROBLEM L8.1. BUILDING MODEL FOR AN EARTHQUAKE." }{TEXT 204 0 "" } }{PARA 200 "" 0 "" {TEXT 204 36 "Model the 6-floor problem in Maple. " }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 68 "Define the 6 by 6 m ass matrix M and Hooke's matrix K for this system" }{TEXT 204 0 "" }} {PARA 200 "" 0 "" {TEXT 204 64 "and convert Mx''=Kx into the system x' '=Ax where A is defined by" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 32 "textbook equation (1), page 437." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 70 "Sanity check: Mass m=3125, and the 6x6 matrix contains fraction 16/5." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 67 "Then find the eigenvalues of the matrix A to six digits, using the " }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 29 "Maple command \"e igenvals(A).\"" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 53 "Sani ty check: All six eigenvalues should be negative." }{TEXT 204 0 "" }} {PARA 200 "" 0 "" {TEXT 204 0 "" }}{PARA 201 "> " 0 "" {MPLTEXT 1 0 46 "# Sample Maple code for a model with 4 floors." }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 52 "# Use maple help to le arn about evalf and eigenvals." }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 72 "# A:=matrix([ [-20,10,0,0], [10,-20,10, 0], [0,10,-20,10],[0,0,10,-10]]);" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 36 "# with(linalg): evalf(eigenvals(A));" } {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 14 "# Problem L8.1" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 35 "# Define k, m and the 6x6 matrix \+ A." }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 36 " # with(linalg): evalf(eigenvals(A));" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 200 "" 0 "" {TEXT 204 55 "PROBLEM L8.2. TABLE OF NATURAL FREQUENCIES AND PERIODS." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 33 "Refer to figure 7.4 .17, page 437." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 0 "" }} {PARA 200 "" 0 "" {TEXT 204 65 "Find the natural angular frequencies \+ omega=sqrt(-lambda) for the" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 53 "six story building and also the corresponding periods" } {TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 68 "2PI/omega, accurate t o six digits. Display the answers in a table ." }{TEXT 204 0 "" }} {PARA 200 "" 0 "" {TEXT 204 70 "Compare with answers in Figure 7.4.17, page 437, for the 7-story case." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 0 "" }}{PARA 201 "> " 0 "" {MPLTEXT 1 0 48 "# Sample code fo r a 4x3 table, 4-story building." }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 48 "# Use maple help to learn about nops an d printf." }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 64 "# ev:=[-10,-1.206147582,-35.32088886,-23.47296354]: n:=nops(ev ):" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 33 " # Omega:=lambda -> sqrt(-lambda):" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 38 "# format:=\"%10.6f %10.6f %10.6f\\n\" : " }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 74 " # seq(printf(format,ev[i],Omega(ev[i]),2*evalf(Pi)/Omega(ev[i])), i=1. .n);" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 14 "# Problem L8.2" }{MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 34 "# ev:=[fill this in]: n:=n ops(ev):" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 69 "# Omega:=lambda -> sqrt(-lambda): format:=\"%10.6f %10.6f %10. 6f\\n\":" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 73 "# seq(printf(format,ev[i],Omega(ev[i]),2*evalf(Pi)/Omega(ev[i])) ,i=1..n);" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 200 "" 0 "" {TEXT 204 41 "PROBLEM L8.3. UNDE TERMINED COEFFICIENTS " }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 45 " STEADY-STATE PERIODIC SOLUTION." }{TEXT 204 0 "" }} {PARA 200 "" 0 "" {TEXT 204 69 "Consider the forced equation x'=Ax+co s(wt)b where b is a constant" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 65 "vector. The earthquake's ground vibration is accounted f or by the" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 47 "extra ter m cos(wt)b, which has period T=2Pi/w." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 68 "The solution x(t) is the 6-vector of excursion s from equilibrium" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 30 " of the corresponding 6 floors." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 65 "Sought here is not the general solution, which certainly contains" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 70 "transient terms, but rather the steady-state periodic solution, which " }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 64 "is known from the theory to have the form x(t)=cos(wt)c for some" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 43 "vector c that depends only on A and b." } {TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 65 "Define b:=0.25*w*w*vector([1,1,1,1,1,1]): in Maple and find the" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 69 "vector c in the undetermined coefficients solution x(t)=cos(wt)c. " }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 64 "Vector c depends on w. As outlined in the textbook, vector c " }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 68 "can be found by solving the linear algebra problem \+ -w^2 c = Ac + b;" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 66 "s ee page 433. Don't print c, as it is too complex; instead, print" } {TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 24 "c[1] as an illustrati on." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 0 "" }}{PARA 201 "> " 0 "" {MPLTEXT 1 0 56 "#Sample code for defining b and A, then solv ing for c " }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 21 "#in the 4-floor case." }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 52 "# See maple help to learn about v ector and linsolve." }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 47 "# w:='w': u:=w*w: b:=0.25*u*vector([1,1,1,1]):" } {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 73 "# A: =matrix([ [-20,10,0,0], [10,-20,10,0], [0,10,-20,10],[0,0,10,-10]]);" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 31 "# Au :=evalm(A+u*diag(1,1,1,1));" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 23 "# c:=linsolve(Au,-b): " }{MPLTEXT 1 0 0 "" } }}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 17 "# evalf(c[1],2);" } {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 14 "# PROBLEM L8.3" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 28 "# Define w, u, b, A, Au, c " } {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 16 "# eva lf(c[1],2);" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 200 "" 0 "" {TEXT 204 34 "PROBLEM L8 .4. PRACTICAL RESONANCE." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 64 "Consider the forced equation x'=Ax+cos(wt)b of L8.3 above wi th" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 34 "b:=0.25*w*w*vect or([1,1,1,1,1,1])." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 64 " Practical resonance can occur if a component of x(t) has large" } {TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 70 "amplitude compared to the vector norm of b. For example, an earthquake" }{TEXT 204 0 "" }} {PARA 200 "" 0 "" {TEXT 204 61 "might cause a small 3-inch excursion o n level ground, but the" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 66 "building's floors might have 50-inch excursions, enough to des troy" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 13 "the building." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 0 "" }}{PARA 200 "" 0 " " {TEXT 204 68 "Let Max(c) denote the maximum modulus of the component s of vector c." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 70 "Plot g(T)=Max(c(w)) with w=(2*Pi)/T for periods T=0 to T=6, ordinates" } {TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 66 "Max=0 to Max=10, the \+ vector c(w) being the answer produced in L8.3" }{TEXT 204 0 "" }} {PARA 200 "" 0 "" {TEXT 204 68 "above. Compare your figure to the tex tbook Figure 7.4.18, page 438." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 0 "" }}{PARA 201 "> " 0 "" {MPLTEXT 1 0 68 "# Sample maple c ode to define the function Max(c), 4-floor building." }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 64 "# Use maple help to \+ learn about norm, vector, subs and linsolve." }{MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 16 "# with(linalg):" } {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 47 "# w: ='w': Max:= c -> norm(c,infinity); u:=w*w:" }{MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 33 "# b:=0.25*w*w*vector([1,1 ,1,1]):" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 73 "# A:=matrix([ [-20,10,0,0], [10,-20,10,0], [0,10,-20,10], [0,0,1 0,-10]]);" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 32 "# Au:=evalm(A+u*diag(1,1,1,1));" }{MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 39 "# C:=ww -> subs(w=ww,lins olve(Au,-b)):" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 51 "# plot(Max(C(2*Pi/r)),r=0..6,0..10,numpoints=150);" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 48 "# PROBLEM L8.4. WARNING: S ave your file often!!!" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 49 "# w:='w': Max:= c -> norm(c,infinity): u:=w*w: " }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 10 "# D efine b" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 10 "# Define A" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 11 "# Define Au" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 " > " 0 "" {MPLTEXT 1 0 12 "# Define C " }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 50 "# plot(Max(C(2*Pi/r)),r=0..6,0..1 0,numpoints=150);" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 200 "" 0 "" {TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 32 "PROBLEM L8.5. EARTHQUAKE DAMAGE." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 61 "The maximum amplitude plot of L8.4 can be used to de tect the " }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 32 "of earthq uake damage for a given" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 63 "ground vibration of period T. A ground vibration (1/4)cos(wt) ," }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 38 "T=2*Pi/w, will be assumed, as in L8.4." }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 81 "(a) Replot the amplitudes in L8 .4 for periods 1.5 to 5.5 and amplitudes 5 to 10. " }{TEXT 204 0 "" }} {PARA 200 "" 0 "" {TEXT 204 27 "There will be five spikes. " }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 61 "(b) Create five zoom-in plo ts, one for each spike, choosing a" }{TEXT 204 0 "" }}{PARA 200 "" 0 " " {TEXT 204 38 "T-interval that shows the full spike. " }{TEXT 204 0 " " }}{PARA 200 "" 0 "" {TEXT 204 67 "(c) Determine from the five zoom-i n plots approximate intervals for" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 64 "the period T such that some floor in the building wil l undergo" }{TEXT 204 0 "" }}{PARA 200 "" 0 "" {TEXT 204 48 "excursion s from equilibrium in excess of 5 feet." }{TEXT 204 0 "" }}{PARA 200 " " 0 "" {TEXT 204 0 "" }}{PARA 201 "> " 0 "" {MPLTEXT 1 0 63 "# Example : Zoom-in on a spike for amplitudes 5 feet to 10 feet," }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 22 "#periods 1.97 to 2 .01." }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 59 "#with(linalg): w:='w': Max:= c -> norm(c,infinity); u:=w*w:" } {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 80 "#Au:= matrix([ [-20+u,10,0,0], [10,-20+u,10,0], [0,10,-20+u,10],[0,0,10,-10+ u]]);" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 31 "#b:=0.25*w*w*vector([1,1,1,1]):" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 37 "#C:=ww -> subs(w=ww,linsolve(Au,- b)):" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 55 "#plot(Max(C(2*Pi/r)),r=1.97..2,01,5..10,numpoints=150);" } {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 47 "# PROBLEM L8.5. WARNING: S ave your file often!!" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 " " {MPLTEXT 1 0 29 "#(a) Re-plot the five spikes." }{MPLTEXT 1 0 0 "" } }}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 54 "# plot(Max(C(2*Pi/r)),r= 1.5..5.5,5..10,numpoints=150);" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 30 "#(b) Plot five zoom-in graphs." } {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 65 "# one :=1.79..1.83:plot(Max(C(2*Pi/r)),r=one,5..10,numpoints=150);" } {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 58 "# two :=???:plot(Max(C(2*Pi/r)),r=two,5..10,numpoints=150);" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 62 "# three:=???:plot( Max(C(2*Pi/r)),r=three,5..10,numpoints=150);" }{MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 60 "# four:=???:plot(Max(C(2*P i/r)),r=four,5..10,numpoints=150);" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 60 "# five:=???:plot(Max(C(2*Pi/r)),r =five,5..10,numpoints=150);" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 25 "#(c) Print period ranges." }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 42 "# PeriodRanges:=[one,t wo,three,four,five];" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 202 "" 0 "" {TEXT -1 0 "" }}}{MARK "0 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }