-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathextendModalAnalysis.html
151 lines (112 loc) · 7.98 KB
/
extendModalAnalysis.html
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
<!DOCTYPE html>
<html>
<head>
<script src="javascripts/flot/jquery.js"></script>
<script src="javascripts/jqmath/jqmath-etc-0.2.0.min.js"></script>
<link rel="stylesheet" href="javascripts/jqmath/jqmath-0.2.0.css">
<link rel="stylesheet" type="text/css" media="screen" href="stylesheets/style.css">
<title>Generalized Modal Analysis</title>
</head>
<body>
<div id="main">
<div id="menu">
<a href="https://github.com/llarsen71/GMA1D"> Git Repository </a>
</div><!--menu-->
<div id="title">
<h1>Extending Modal Analysis to Nonlinear Systems</h1>
</div><!--title-->
<p><a href="ModalAnalysis.html" style="text-decoration: none;">← Previous Topic: Linear Modal Analysis</a>. </p>
<h2> Generalized Modal Coordinates</h2>
<p>So, is there a way to extend modal analysis to nonlinear systems? There is, and addressing that question is the purpose of this
website. The method discussed on this site applies to linear and nonlinear systems and reduces directly to classical modal
analysis for linear systems. I refer to this approach as <b class="emphasis">Generalized Modal Analysis</b> or GMA.</p>
<div class="box-right" style="margin-top: 0em; width: 520px;">
<h1>Definition</h1>
<b>Generalized Modal Coordinate</b> - A function $q(x_1,x_2,...,x_n)$ that when used as a coordinate for a system of ODEs results in a
decoupled ODE of the form ${dq}/{dt} = h(q)$ where $h(q)$ is any function of $q$ alone.
</div><!--box-right-->
<p>For nonlinear systems, the definition of modal analysis has to be generalized. First we drop the restriction that
the coordinate function be a linear function - any function which decouples the system of ODEs will be
allowed. Second, we relax the restriction on the form of the decoupled ODE - any decoupled ODE will be allowed. The key to calculating
modes for a nonlinear system is to express a generalized mode in terms of the chain rule rather than using eigenvectors for the
transformation.</p>
<h3> First Order Form </h3>
<p>In order to apply the chain rule to a system of ODEs, it is useful to write the system in first order form, where the term on the
left hand side of the equation is a lone single derivative, and the term on the right contains only state variables:</p>
<p>${dx_1}/{dt} = f_1(x_1, x_2, ..., x_n)$<br>
${dx_2}/{dt} = f_2(x_1, x_2, ..., x_n)$<br>
...<br>
${dx_n}/{dt} = f_n(x_1, x_2, ..., x_n)$</p>
<h4> Example System </h4>
<p>As an example, the same linear system will be used as on the previous page so that we can verify that we get the same linear modes.
However, we will also see that there are other generalized modes associated with a linear system. The linear system in first order
form is: </p>
<p>${dx_1}/{dt} = -2 x_1 + x_2$<br>
${dx_2}/{dt} = x_1 - 2 x_2$</p>
<h3><a name="vector">Vector Field Form</a></h3>
<p>A closely related form that is useful for expressing the chain rule is vector field form, where $∂/{∂ x_i}$ is used to
represent a coordinate vector for coordinate $x_i$. This might seem like a strange notation for a basis vector if you have done
multivariate calculus. However, there is a good reason for this notation. This is a case where the notation helps you know how
to apply a transformation of coordinates correctly. The system of ODEs represented in vector field form is:</p>
<p>(1) $V = {dx_1}/{dt} ∂/{∂ x_1} + {dx_2}/{dt} ∂/{∂ x_2} + ... + {dx_n}/{dt} ∂/{∂ x_n}$<br>
(2) $V = f_1 ∂/{∂ x_1} + f_2 ∂/{∂ x_2} + ... + f_n ∂/{∂ x_n}$</p>
<p>Note that when the ODE vector field is expressed in form (1), it kind of resembles the chain rule.</p>
<h4>Example Vector Field</h4>
<p>The example system in vector field form is:</p>
<p>$V = (-2 x_1 + x_2) {∂}/{∂ x_1} + (x_1 - 2 x_2) {∂}/{∂ x_2}$</p>
<h3>Coordinate Transformation (The Chain Rule)</h3>
<p>If a new coordinate function $q(x_1, x_2, ..., x_n)$ is defined, ${dq}/{dt}$ is calculated by:</p>
<p>(3) ${dq}/{dt} = V(q) = {dx_1}/{dt} {∂ q}/{∂ x_1} + {dx_2}/{dt} {∂ q}/{∂ x_2} + ... + {dx_n}/{dt} {∂ q}/{∂ x_n} = f_1 {∂ q}/{∂ x_1} + f_2 {∂ q}/{∂ x_2} + ... + f_n {∂ q}/{∂ x_n}$</p>
<p> Supposed we defined a complete set of new coordinates $q_1, q_2, ... , q_n$, where these coordinate are expressed as functions of
$x_1, x_2, ..., x_n$, and would like to represent the vector field in terms of the new coordinate system. The new vector field is calculated
as follows:</p>
<p>(3) $W = V(q_1) {∂}/{∂ q_1} + V(q_2) {∂}/{∂ q_2} + ... + V(q_n) {∂}/{∂ q_n} = {dq_1}/{dt} {∂}/{∂ q_1} + {dq_2}/{dt} {∂}/{∂ q_2} + ... + {dq_n}/{dt} {∂}/{∂ q_n}$</p>
<p> Combining the definition of $V(q_i)$ from (3) with (4), you can see that this is the chain rule, and the the vector field notation
make the transformation of a system of ODEs from one coordinate system more obvious. For this reason, it is a preferred notation for
generalizing modal analysis to nonlinear systems.</p>
<h4>Example Transformations</h4>
<p>Lets try a couple of random examples to show how you calculate ${dq}/{dt}$ for a function $q$. Note that the
coordinate examples chosen below are NOT modal coordinates.</p>
<hr>
<table>
<tr><td><div class="box"><b>Example 1</b></div></td><td>$q_1=x_1^2 + x_2^3$<br>
${dq_1}/{dt} = V(q_1) = (-2 x_1 + x_2) {∂ (x_1^2 + x_2^3)}/{∂ x_1} + (x_1 - 2 x_2) {∂(x_1^2 + x_2^3)}/{∂ x_2}$<br>
${dq_1}/{dt} = (-2 x_1 + x_2) 2 x_1 + (x_1 - 2 x_2) 3 x_2^2$</td></tr>
</table>
<hr>
<table>
<tr><td><div class="box"><b>Example 2</b></div></td><td>$q_2=x_1 x_2$<br>
${dq_2}/{dt} = V(q_2) = (-2 x_1 + x_2) {∂ (x_1 x_2)}/{∂ x_1} + (x_1 - 2 x_2) {∂(x_1 x_2)}/{∂ x_2} = (-2 x_1 + x_2) x_2 + (x_1 - 2 x_2) x_1$</td></tr>
</table>
<hr>
<h3>Generalized Modal Analysis</h3>
<p>For generalized modal analysis the coordinate transformation from (3) is combined with the decoupling condition
${dq}/{dt} = h(q)$. A coordinate function satifying this equation is a generalized modal coordinate and results
in a decoupled ODE which can be solved independent of any other choice of coordinate for the system of ODEs.</p>
<p>${dq}/{dt} = V(q) = f_1 {∂ q}/{∂ x_1} + f_2 {∂ q}/{∂ x_2} + ... + f_n {∂ q}/{∂ x_n} = h(q)$<br>
$f_1 {∂ q}/{∂ x_1} + f_2 {∂ q}/{∂ x_2} + ... + f_n {∂ q}/{∂ x_n} = h(q)$
</p>
<h4>Verify the Linear Modes</h4>
<p>On the <a href="ModalAnalysis.html#modes">Linear Modal Analysis</a> page, the linear modes of the example system were calculated.
Applying the chain rule gives:</p>
<hr>
<table>
<tr><td><div class="box"><b>Mode 1</b></div> <div class="pad"></td><td>$q_1 = s x_1 + s x_2$<br>
${dq_1}/{dt} = V(q_1) = (-2 x_1 + x_2) {∂ (s x_1 + s x_2)}/{∂ x_1} + (x_1 - 2 x_2) {∂(s x_1 + s x_2)}/{∂ x_2}$<br>
${dq_1}/{dt} = (-2 x_1 + x_2) s + (x_1 - 2 x_2) s = - s x_1 - s x_2$<br>
${dq_1}/{dt} = -q_1$<br><br>
This matches the <a href="ModalAnalysis.html#eqn">modal equation</a> for mode 1 calculated using the standard method.</td></tr>
</table>
<hr>
<table>
<tr><td><div class="box"><b>Mode 2</b></div> <div class="pad"></td><td>$q_2 = s x_1 - s x_2$<br>
${dq_2}/{dt} = V(q_2) = (-2 x_1 + x_2) {∂ (s x_1 - s x_2)}/{∂ x_1} + (x_1 - 2 x_2) {∂(s x_1 - s x_2)}/{∂ x_2}$<br>
${dq_2}/{dt} = (-2 x_1 + x_2) s - (x_1 - 2 x_2) s = - 3 (s x_1 - s x_2)$<br>
${dq_2}/{dt} = - 3 q_2$<br><br>
This matches the <a href="ModalAnalysis.html#eqn">modal equation</a> for mode 2 calculated using the standard method.</td></tr>
</table>
<hr>
<p>
</div><!--main-->
</body>
</html>