From 91ca6068b591e2b6d022c2f459c2fa78b637cf2a Mon Sep 17 00:00:00 2001 From: Thomas Dehaeze Date: Tue, 17 Dec 2024 16:38:55 +0100 Subject: [PATCH] Update Content - 2024-12-17 --- ...inroy02_model_desig_flexur_joint_stewar.md | 8 +-- content/article/mcinroy99_dynam.md | 4 +- .../book/du19_multi_actuat_system_contr.md | 4 +- content/book/ewins00_modal.md | 72 +++++++++---------- content/book/hatch00_vibrat_matlab_ansys.md | 6 +- content/book/morrison16_groun_shiel.md | 6 +- ..._vibrat_contr_activ_struc_fourt_edition.md | 22 +++--- .../book/skogestad07_multiv_feedb_contr.md | 18 ++--- .../li01_simul_fault_vibrat_isolat_point.md | 2 +- content/phdthesis/rankers98_machin.md | 26 +++---- content/posts/test.md | 4 +- content/posts/test2.md | 2 +- content/zettels/test.md | 4 +- 13 files changed, 89 insertions(+), 89 deletions(-) diff --git a/content/article/mcinroy02_model_desig_flexur_joint_stewar.md b/content/article/mcinroy02_model_desig_flexur_joint_stewar.md index 740d531..2df12ff 100644 --- a/content/article/mcinroy02_model_desig_flexur_joint_stewar.md +++ b/content/article/mcinroy02_model_desig_flexur_joint_stewar.md @@ -115,7 +115,7 @@ where: - \\(\mathcal{F}\_e\\) represents a vector of exogenous generalized forces applied at the center of mass - \\(g\\) is the gravity vector -By combining , and , a single equation describing the dynamics of a flexure jointed hexapod can be found: +By combining \eqref{eq:strut\_dynamics\_vec}, \eqref{eq:payload\_dynamics} and \eqref{eq:generalized\_force}, a single equation describing the dynamics of a flexure jointed hexapod can be found: \begin{equation} {}^UJ^T [ f\_m - M\_s \ddot{l} - B \dot{l} - K(l - l\_r) - M\_s \ddot{q}\_u - M\_s g\_u + M\_s v\_2] + \mathcal{F}\_e - \begin{bmatrix} mg \\\ 0\_{3\times 1} \end{bmatrix} = M\_x \ddot{\mathcal{X}} + c(\omega) \label{eq:eom\_fjh} @@ -200,7 +200,7 @@ In order to get dominance at low frequencies, the hexapod must be designed so th This puts a limit on the rotational stiffness of the flexure joint and shows that as the strut is made softer (by decreasing \\(k\\)), the spherical flexure joint must be made proportionately softer. -By satisfying , \\(f\_p\\) can be aligned with the strut for frequencies much below the spherical joint's resonance mode: +By satisfying \eqref{eq:cond\_stiff}, \\(f\_p\\) can be aligned with the strut for frequencies much below the spherical joint's resonance mode: \\[ \omega \ll \sqrt{\frac{k\_r}{m\_s l^2}} \rightarrow x\_{\text{gain}\_\omega} \approx \frac{k}{k\_r/l^2} \gg 1 \\] At frequencies much above the strut's resonance mode, \\(f\_p\\) is not dominated by its \\(x\\) component: \\[ \omega \gg \sqrt{\frac{k}{m\_s}} \rightarrow x\_{\text{gain}\_\omega} \approx 1 \\] @@ -225,14 +225,14 @@ In this case, it is reasonable to use:
-By designing the flexure jointed hexapod and its controller so that both and are met, the dynamics of the hexapod can be greatly reduced in complexity. +By designing the flexure jointed hexapod and its controller so that both \eqref{eq:cond\_stiff} and \eqref{eq:cond\_bandwidth} are met, the dynamics of the hexapod can be greatly reduced in complexity.
## Relationships between joint and cartesian space {#relationships-between-joint-and-cartesian-space} -Equation is not suitable for control analysis and design because \\(\ddot{\mathcal{X}}\\) is implicitly a function of \\(\ddot{q}\_u\\). +Equation \eqref{eq:eom\_fjh} is not suitable for control analysis and design because \\(\ddot{\mathcal{X}}\\) is implicitly a function of \\(\ddot{q}\_u\\). This section will derive this implicit relationship. Let denote: diff --git a/content/article/mcinroy99_dynam.md b/content/article/mcinroy99_dynam.md index 124845f..6b2cdb4 100644 --- a/content/article/mcinroy99_dynam.md +++ b/content/article/mcinroy99_dynam.md @@ -53,7 +53,7 @@ In order to provide low frequency passive vibration isolation, the hard actuator
Table 1: - Definition of quantities on 2 + Definition of quantities on 2
| **Symbol** | **Meaning** | @@ -142,7 +142,7 @@ where: - \\(\mathcal{F}\_e\\) represents a vector of exogenous generalized forces applied at the center of mass - \\(g\\) is the gravity vector -By combining , and , a single equation describing the dynamics of a flexure jointed hexapod can be found: +By combining \eqref{eq:strut\_dynamics\_vec}, \eqref{eq:payload\_dynamics} and \eqref{eq:generalized\_force}, a single equation describing the dynamics of a flexure jointed hexapod can be found: \begin{aligned} & {}^UJ^T [ f\_m - M\_s \ddot{l} - B \dot{l} - K(l - l\_r) - M\_s \ddot{q}\_u\\\\ diff --git a/content/book/du19_multi_actuat_system_contr.md b/content/book/du19_multi_actuat_system_contr.md index 46215d8..faebcca 100644 --- a/content/book/du19_multi_actuat_system_contr.md +++ b/content/book/du19_multi_actuat_system_contr.md @@ -206,7 +206,7 @@ is satisfied, where \\(T\_{zw}\\) is the transfer function from \\(w\\) to \\(z\ {{< figure src="/ox-hugo/du19_h_inf_diagram.png" caption="Figure 6: Block diagram for \\(\mathcal{H}\_\infty\\) loop shaping method to design the controller \\(C(s)\\) with the weighting function \\(W(s)\\)" >}} -Equation [ 1](#orgc94b8e5) means that \\(S(s)\\) can be shaped similarly to the inverse of the chosen weighting function \\(W(s)\\). +Equation [ 1](#orgcf76ccd) means that \\(S(s)\\) can be shaped similarly to the inverse of the chosen weighting function \\(W(s)\\). One form of \\(W(s)\\) is taken as \begin{equation} @@ -339,7 +339,7 @@ A decoupled control structure can be used for the three-stage actuation system ( The overall sensitivity function is \\[ S(z) = \approx S\_v(z) S\_p(z) S\_m(z) \\] -with \\(S\_v(z)\\) and \\(S\_p(z)\\) are defined in equation [ 1](#org5626095) and +with \\(S\_v(z)\\) and \\(S\_p(z)\\) are defined in equation [ 1](#org40d0f02) and \\[ S\_m(z) = \frac{1}{1 + P\_m(z) C\_m(z)} \\] Denote the dual-stage open-loop transfer function as \\(G\_d\\) diff --git a/content/book/ewins00_modal.md b/content/book/ewins00_modal.md index 01ff6f4..633d711 100644 --- a/content/book/ewins00_modal.md +++ b/content/book/ewins00_modal.md @@ -221,7 +221,7 @@ This process itself falls into two stages: Most of the effort goes into this second stage, which is widely referred to as "modal parameter extraction", or simply as "modal analysis". We have seen that we can predict the form of the FRF plots for a multi degree-of-freedom system, and that these are directly related to the modal properties of that system. -The great majority of the modal analysis effort involves **curve-fitting** an expression such as equation to the measured FRF and thereby finding the appropriate modal parameters. +The great majority of the modal analysis effort involves **curve-fitting** an expression such as equation \eqref{eq:frf\_modal} to the measured FRF and thereby finding the appropriate modal parameters. A completely general curve-fitting approach is possible but generally inefficient. Mathematically, we can take an equation of the form @@ -477,11 +477,11 @@ where \\(\eta\\) is the **structural damping loss factor** and replaces the crit #### Alternative Forms of FRF {#alternative-forms-of-frf} -So far we have defined our receptance frequency response function \\(\alpha(\omega)\\) as the ratio between a harmonic displacement response and the harmonic force . +So far we have defined our receptance frequency response function \\(\alpha(\omega)\\) as the ratio between a harmonic displacement response and the harmonic force \eqref{eq:receptance}. This ratio is complex: we can look at its **amplitude** ratio \\(|\alpha(\omega)|\\) and its **phase** angle \\(\theta\_\alpha(\omega)\\). -We could have selected the response velocity \\(v(t)\\) as the output quantity and defined an alternative frequency response function . -Similarly we could use the acceleration parameter so we could define a third FRF parameter . +We could have selected the response velocity \\(v(t)\\) as the output quantity and defined an alternative frequency response function \eqref{eq:mobility}. +Similarly we could use the acceleration parameter so we could define a third FRF parameter \eqref{eq:inertance}.
@@ -588,7 +588,7 @@ This type of display is not widely used as we cannot use logarithmic axes (as we ##### Real part and Imaginary part of reciprocal FRF {#real-part-and-imaginary-part-of-reciprocal-frf} -It can be seen from the expression of the inverse receptance that the Real part depends entirely on the mass and stiffness properties while the Imaginary part is a only function of the damping. +It can be seen from the expression of the inverse receptance \eqref{eq:dynamic\_stiffness} that the Real part depends entirely on the mass and stiffness properties while the Imaginary part is a only function of the damping. [ 5](#org-target--fig-inverse-frf-mixed) shows an example of a plot of a system with a combination of both viscous and structural damping. The imaginary part is a straight line whose slope is given by the viscous damping rate \\(c\\) and whose intercept at \\(\omega = 0\\) is provided by the structural damping coefficient \\(d\\). @@ -629,7 +629,7 @@ This makes the Nyquist plot very effective for modal testing applications. #### Free Vibration Solution - The modal Properties {#free-vibration-solution-the-modal-properties} -For an undamped MDOF system, with \\(N\\) degrees of freedom, the governing equations of motion can be written in matrix form . +For an undamped MDOF system, with \\(N\\) degrees of freedom, the governing equations of motion can be written in matrix form \eqref{eq:undamped\_mdof}.
@@ -645,7 +645,7 @@ where \\([M]\\) and \\([K]\\) are \\(N\times N\\) mass and stiffness matrices, a We shall consider first the free vibration solution by taking \\(f(t) = 0\\). In this case, we assume that a solution exists of the form \\(\\{x(t)\\} = \\{X\\} e^{i \omega t}\\) where \\(\\{X\\}\\) is an \\(N \times 1\\) vector of time-independent amplitudes. -Substitution of this condition into leads to +Substitution of this condition into \eqref{eq:undamped\_mdof} leads to \begin{equation} \label{eq:free\_eom\_mdof} \left( [K] - \omega^2 [M] \right) \\{X\\} e^{i\omega t} = \\{0\\} @@ -655,7 +655,7 @@ for which the non trivial solutions are those which satisfy \\[ \det \left| [K] - \omega^2 [M] \right| = 0 \\] from which we can find \\(N\\) values of \\(\omega^2\\) corresponding to the undamped system's **natural frequencies**. -Substituting any of these back into yields a corresponding set of relative values for \\(\\{X\\}\\): \\(\\{\psi\\}\_r\\) the so-called **mode shape** corresponding to that natural frequency. +Substituting any of these back into \eqref{eq:free\_eom\_mdof} yields a corresponding set of relative values for \\(\\{X\\}\\): \\(\\{\psi\\}\_r\\) the so-called **mode shape** corresponding to that natural frequency.
@@ -792,7 +792,7 @@ An alternative means of deriving the FRF parameters is used which makes use of t \\[ [K] - \omega^2 [M] = [\alpha(\omega)]^{-1} \\] Pre-multiply both sides by \\([\Phi]^T\\) and post-multiply both sides by \\([\Phi]\\) to obtain \\[ [\Phi]^T ([K] - \omega^2 [M]) [\Phi] = [\Phi]^T [\alpha(\omega)]^{-1} [\Phi] \\] -which leads to . +which leads to \eqref{eq:receptance\_modal}.
@@ -802,7 +802,7 @@ Receptance FRF matrix - Modal Properties; [\alpha(\omega)] = [\Phi] \left[ \bar{\omega}\_r^2 - \omega^2 \right]^{-1} [\Phi]^T \label{eq:receptance\_modal} \end{equation} -Equation permits us to compute any individual FRF parameters \\(\alpha\_{jk}(\omega)\\) using the following formula +Equation \eqref{eq:receptance\_modal} permits us to compute any individual FRF parameters \\(\alpha\_{jk}(\omega)\\) using the following formula \begin{align} \alpha\_{jk}(\omega) &= \sum\_{r=1}^N \frac{\phi\_{jr} \phi\_{kr}}{\bar{\omega}\_r^2 - \omega^2}\\\\ @@ -816,7 +816,7 @@ where \\({}\_rA\_{jk}\\) is called the **modal constant**.
-It is clear from equation that the receptance matrix \\([\alpha(\omega)]\\) is **symmetric** and this will be recognized as the **principle of reciprocity**. +It is clear from equation \eqref{eq:receptance\_modal} that the receptance matrix \\([\alpha(\omega)]\\) is **symmetric** and this will be recognized as the **principle of reciprocity**. This principle of reciprocity applies to many structural characteristics. @@ -954,7 +954,7 @@ From this full matrix equation, we have: Having derived an expression for the general term in the frequency response function matrix \\(\alpha\_{jk}(\omega)\\), it is appropriate to consider next the analysis of a situation where the system is **excited simultaneously at several points**. -The general behavior for this case is governed by equation with solution . +The general behavior for this case is governed by equation \eqref{eq:force\_response\_eom} with solution \eqref{eq:force\_response\_eom\_solution}. However, a more explicit form of the solution is \begin{equation} \label{eq:ods} @@ -977,7 +977,7 @@ The properties of the normal modes of the undamped system are of interest becaus
-We are seeking an excitation vector \\(\\{F\\}\\) such that the **response** \\(\\{X\\}\\) **consists of a single modal component** so that all terms in but one is zero. +We are seeking an excitation vector \\(\\{F\\}\\) such that the **response** \\(\\{X\\}\\) **consists of a single modal component** so that all terms in \eqref{eq:ods} but one is zero. This can be attained if \\(\\{F\\}\\) is chosen such that \\[ \\{\phi\_r\\}^T \\{F\\}\_s = 0, \ r \neq s \\] @@ -1059,7 +1059,7 @@ where \\(\omega\_r\\) is the **natural frequency** and \\(\xi\_r\\) is the **cri When the modes \\(r\\) and \\(q\\) are a complex conjugate pair: \\[ s\_r = \omega\_r \left( -\xi\_r - i\sqrt{1 - \xi\_r^2} \right); \quad \\{\psi\\}\_q = \\{\psi\\}\_r^\* \\] -From equations , we can obtain +From equations \eqref{eq:viscous\_damping\_orthogonality}, we can obtain \begin{align} 2 \omega\_r \xi\_r &= \frac{\\{\psi\\}\_r^H [C] \\{\psi\\}\_r}{\\{\psi\\}\_r^H [M] \\{\psi\\}\_r} = \frac{c\_r}{m\_r} \\\\ @@ -1352,7 +1352,7 @@ One these two series are available, the FRF can be defined at the same set of fr ##### Analysis via Fourier transform {#analysis-via-fourier-transform} -For most transient cases, the input function \\(f(t)\\) will satisfy the **Dirichlet condition** and so its Fourier Transform \\(F(\omega)\\) can be computed from . +For most transient cases, the input function \\(f(t)\\) will satisfy the **Dirichlet condition** and so its Fourier Transform \\(F(\omega)\\) can be computed from \eqref{eq:fourier\_transform}. \begin{equation} \label{eq:fourier\_transform} F(\omega) = \frac{1}{2 \pi} \int\_{-\infty}^\infty f(t) e^{i\omega t} dt @@ -1500,10 +1500,10 @@ However, the same equation can be transform to the frequency domain \tcmbox{ S\_{xx}(\omega) = \left| H(\omega) \right|^2 S\_{ff}(\omega) } \end{equation} -Although very convenient, equation does not provide a complete description of the random vibration conditions. +Although very convenient, equation \eqref{eq:psd\_input\_output} does not provide a complete description of the random vibration conditions. Further, it is clear that **is could not be used to determine the FRF** from measurement of excitation and response because it **contains only the modulus** of \\(H(\omega)\\), the phase information begin omitted from this formula. -A second equation is required and this may be obtain by a similar analysis, two alternative formulas can be obtained . +A second equation is required and this may be obtain by a similar analysis, two alternative formulas can be obtained \eqref{eq:cross\_relation\_alternatives}.
@@ -1517,8 +1517,8 @@ A second equation is required and this may be obtain by a similar analysis, two ##### To derive FRF from random vibration signals {#to-derive-frf-from-random-vibration-signals} -The pair of equations provides the basic of determining a system's FRF properties from the measurements and analysis of a random vibration test. -Using either of them, we have a simple formula for determining the FRF from estimates of the relevant spectral densities . +The pair of equations \eqref{eq:cross\_relation\_alternatives} provides the basic of determining a system's FRF properties from the measurements and analysis of a random vibration test. +Using either of them, we have a simple formula for determining the FRF from estimates of the relevant spectral densities \eqref{eq:H1} \eqref{eq:H2}.
@@ -1555,7 +1555,7 @@ Then in [ 13](#org-target--fig-frf-feedback-model) is given a more detailed and | width=\linewidth | width=\linewidth | In this configuration, it can be seen that there are two feedback mechanisms which apply. -We then introduce an alternative formula which is available for the determination of the system FRF from measurements of the input and output quantities . +We then introduce an alternative formula which is available for the determination of the system FRF from measurements of the input and output quantities \eqref{eq:H3}.
@@ -1694,7 +1694,7 @@ First, if we have a **modal incompleteness** (\\(m to re-construct the system mass and stiffness matrices. +In both reduced-model cases, it is not possible to use equation \eqref{eq:spatial\_model\_from\_modal} to re-construct the system mass and stiffness matrices. First of all because the eigen matrices are generally singular and even if it is not, the obtained mass and stiffness matrices produced have no physical significance and should not be used. @@ -2699,7 +2699,7 @@ Then, the full \\(6 \times 6\\) mobility matrix can be measured, however this pr Other methods for measuring rotational effects include specially developed rotational accelerometers and shakers. However, there is a major problem that is encountered when measuring rotational FRF: the translational components of the structure's movement tends to overshadow those due to the rotational motions. -For example, the magnitude of the difference in equation is often of the order of \\(\SI{1}{\\%}\\) of the two individual values which is similar to the transverse sensitivity of the accelerometers: potential errors in rotations are thus enormous. +For example, the magnitude of the difference in equation \eqref{eq:rotational\_diff} is often of the order of \\(\SI{1}{\\%}\\) of the two individual values which is similar to the transverse sensitivity of the accelerometers: potential errors in rotations are thus enormous. ### Multi-point excitation methods {#multi-point-excitation-methods} @@ -2738,7 +2738,7 @@ The two vectors are related by the system's FRF properties as: \\{X\\}\_{n\times 1} = [H(\omega)]\_{n\times p} \\{F\\}\_{p\times 1} \end{equation} -However, it is not possible to derive the FRF matrix from the single equation , because there will be insufficient data in the two vectors (one of length \\(p\\), the other of length \\(n\\)) to define completely the \\(n\times p\\) FRF matrix. +However, it is not possible to derive the FRF matrix from the single equation \eqref{eq:mpss\_equation}, because there will be insufficient data in the two vectors (one of length \\(p\\), the other of length \\(n\\)) to define completely the \\(n\times p\\) FRF matrix. What is required is to make a series of \\(p^\prime\\) measurements of the same basic type using different excitation vectors \\(\\{F\\}\_i\\) that should be chosen such that the forcing matrix \\([F]\_{p\times p^\prime} = [\\{F\\}\_1, \dots, \\{F\\}\_p]\\) is non-singular. This can be assured if: @@ -3226,13 +3226,13 @@ For any frequency \\(\omega\\), we have the following relationship: \end{aligned} \end{equation} -From , we obtain: +From \eqref{eq:modal\_circle\_tan}, we obtain: \begin{equation} \label{eq:modal\_circle\_omega} \omega^2 = \omega\_r^2 \left(1 - \eta\_r \tan\left(\frac{\theta}{2}\right) \right) \end{equation} -If we differentiate with respect to \\(\theta\\), we obtain: +If we differentiate \eqref{eq:modal\_circle\_omega} with respect to \\(\theta\\), we obtain: \begin{equation} \frac{d\omega^2}{d\theta} = \frac{-\omega\_r^2 \eta\_r}{2} \frac{\left(1 - (\omega/\omega\_r)^2\right)^2}{\eta\_r^2} @@ -3317,10 +3317,10 @@ The sequence is: Then we obtain the **center** and **radius** of the circle and the **quality factor** is the mean square deviation of the chosen points from the circle. 3. **Locate natural frequency, obtain damping estimate**. The rate of sweep through the region is estimated numerically and the frequency at which it reaches the maximum is deduced. - At the same time, an estimate of the damping is derived using . + At the same time, an estimate of the damping is derived using \eqref{eq:estimate\_damping\_sweep\_rate}. A typical example is shown on [Figure 29](#figure--fig:circle-fit-natural-frequency). 4. **Calculate multiple damping estimates, and scatter**. - A set of damping estimates using all possible combination of the selected data points are computed using . + A set of damping estimates using all possible combination of the selected data points are computed using \eqref{eq:estimate\_damping}. Then, we can choose the damping estimate to be the mean value. We also look at the distribution of the obtained damping estimates as is permits a useful diagnostic of the quality of the entire analysis: - Good measured data should lead to a smooth plot of these damping estimates, any roughness of the surface can be explained in terms of noise in the original data. @@ -3426,7 +3426,7 @@ we now have sufficient information to extract estimates for the four parameters 3. Plot graphs of \\(m\_R(\Omega)\\) vs \\(\Omega^2\\) and of \\(m\_I(\Omega)\\) vs \\(\Omega^2\\) using the results from step 1., each time using a different measurement points as the fixing frequency \\(\Omega\_j\\) 4. Determine the slopes of the best fit straight lines through these two plots, \\(n\_R\\) and \\(n\_I\\), and their intercepts with the vertical axis \\(d\_R\\) and \\(d\_I\\) -5. Use these four quantities, and equation , to determine the **four modal parameters** required for that mode +5. Use these four quantities, and equation \eqref{eq:modal\_parameters\_formula}, to determine the **four modal parameters** required for that mode This procedure which places more weight to points slightly away from the resonance region is likely to be less sensitive to measurement difficulties of measuring the resonance region. @@ -3491,7 +3491,7 @@ From the sketch, it may be seen that within the frequency range of interest: - the first term tends to approximate to a **mass-like behavior** - the third term approximates to a **stiffness effect** -Thus, we have a basis for the residual terms and shall rewrite equation : +Thus, we have a basis for the residual terms and shall rewrite equation \eqref{eq:sum\_modes:} \begin{equation} H\_{jk}(\omega) \simeq -\frac{1}{\omega^2 M\_{jk}^R} + \sum\_{r=m\_1}^{m\_2} \left( \frac{{}\_rA\_{jk}}{\omega\_r^2 - \omega^2 + i \eta\_r \omega\_r^2} \right) + \frac{1}{K\_{jk}^R} @@ -3551,7 +3551,7 @@ We can write the receptance in the frequency range of interest as: In the previous methods, the second term was assumed to be a constant in the curve-fit procedure for mode \\(r\\). However, if we have good **estimates** for the coefficients which constitutes the second term, for example by having already completed an SDOF analysis, we may remove the restriction on the analysis. -Indeed, suppose we take a set of measured data points around the resonance at \\(\omega\_r\\), and that we can compute the magnitude of the second term in , we then subtract this from the measurement and we obtain adjusted data points that are conform to a true SDOF behavior and we can use the same technique as before to obtain **improved estimated** to the modal parameters of more \\(r\\). +Indeed, suppose we take a set of measured data points around the resonance at \\(\omega\_r\\), and that we can compute the magnitude of the second term in \eqref{eq:second\_term\_refinement}, we then subtract this from the measurement and we obtain adjusted data points that are conform to a true SDOF behavior and we can use the same technique as before to obtain **improved estimated** to the modal parameters of more \\(r\\). This procedure can be repeated iteratively for all the modes in the range of interest and it can significantly enhance the quality of found modal parameters for system with **strong coupling**. @@ -3610,7 +3610,7 @@ If we further increase the generality by attaching a **weighting factor** \\(w\_ is minimized. -This is achieved by differentiating with respect to each unknown in turn, thus generating a set of as many equations as there are unknown: +This is achieved by differentiating \eqref{eq:error\_weighted} with respect to each unknown in turn, thus generating a set of as many equations as there are unknown: \begin{equation} \frac{d E}{d q} = 0; \quad q = {}\_1A\_{jk}, {}\_2A\_{jk}, \dots @@ -3661,7 +3661,7 @@ leading to the modified, but more convenient version actually used in the analys \end{equation} In these expressions, only \\(m\\) modes are included in the theoretical FRF formula: the true number of modes, \\(N\\), is actually one of the **unknowns** to be determined during the analysis. -Equation can be rewritten as follows: +Equation \eqref{eq:rpf\_error} can be rewritten as follows: \begin{equation} \begin{aligned} @@ -3715,7 +3715,7 @@ where \\([X], [Y], [Z], \\{G\\}\\) and \\(\\{F\\}\\) are known measured quantiti \end{equation} Once the solution has been obtained for the coefficients \\(a\_k, \dots , b\_k, \dots\\) then the second stage of the modal analysis can be performed in which the required **modal parameters are derived**. -This is usually done by solving the two polynomial expressions which form the numerator and denominator of equations and : +This is usually done by solving the two polynomial expressions which form the numerator and denominator of equations \eqref{eq:frf\_clasic} and \eqref{eq:frf\_rational:} - the denominator is used to obtain the natural frequencies \\(\omega\_r\\) and damping factors \\(\xi\_r\\) - the numerator is used to determine the complex modal constants \\(A\_r\\) @@ -3943,7 +3943,7 @@ First, we note that from a single FRF curve, \\(H\_{jk}(\omega)\\), it is possib Now, although this gives us the natural frequency and damping properties directly, it does not explicitly yield the mode shape: only a modal constant \\({}\_rA\_{jk}\\) which is formed from the mode shape data. In order to extract the individual elements \\(\phi\_{jr}\\) of the mode shape matrix \\([\Phi]\\), it is necessary to make a series of measurements of specific FRFs including, especially, the point FRF at the excitation position. -If we measure \\(H\_{kk}\\), then by using , we also obtain the specific elements in the mode shape matrix corresponding to the excitation point: +If we measure \\(H\_{kk}\\), then by using \eqref{eq:modal\_model\_from\_frf}, we also obtain the specific elements in the mode shape matrix corresponding to the excitation point: \begin{equation} H\_{kk}(\omega) \longrightarrow \omega\_r, \eta\_r, {}\_rA\_{jk} \longrightarrow \phi\_{kr}; \quad r=1, m @@ -4417,7 +4417,7 @@ In this respect, the demands of the response model are more stringent that those #### Synthesis of FRF curves {#synthesis-of-frf-curves} -One of the implications of equation is that **it is possible to synthesize the FRF curves which were not measured**. +One of the implications of equation \eqref{eq:regenerate\_full\_frf\_matrix} is that **it is possible to synthesize the FRF curves which were not measured**. This arises because if we measured three individual FRF such as \\(H\_{ik}(\omega)\\), \\(H\_{jk}(\omega)\\) and \\(K\_{kk}(\omega)\\), then modal analysis of these yields the modal parameters from which it is possible to generate the FRF \\(H\_{ij}(\omega)\\), \\(H\_{jj}(\omega)\\), etc. However, it must be noted that there is an important **limitation to this procedure** which is highlighted in the example below. @@ -4541,7 +4541,7 @@ from which is would appear that we can write \end{aligned} \end{equation} -However, equation is **only applicable when we have available the complete \\(N \times N\\) modal model**. +However, equation \eqref{eq:m\_k\_from\_modes} is **only applicable when we have available the complete \\(N \times N\\) modal model**. It is much more usual to have an incomplete model in which the eigenvector matrix is rectangle and, as such, is non-invertible. One step which can be made using the incomplete data is the construction of "pseudo" flexibility and inverse-mass matrices. diff --git a/content/book/hatch00_vibrat_matlab_ansys.md b/content/book/hatch00_vibrat_matlab_ansys.md index 1b16ed8..762dd3b 100644 --- a/content/book/hatch00_vibrat_matlab_ansys.md +++ b/content/book/hatch00_vibrat_matlab_ansys.md @@ -252,7 +252,7 @@ where: #### Eigenvalues / Characteristic Equation {#eigenvalues-characteristic-equation} -Re-injecting normal modes into the equation of motion gives the eigenvalue problem: +Re-injecting normal modes \eqref{eq:principal\_mode} into the equation of motion \eqref{eq:tdof\_eom} gives the eigenvalue problem: \begin{equation} (\bm{k} - \omega\_i^2 \bm{m}) \bm{z}\_{mi} = 0 @@ -310,7 +310,7 @@ Virtual interpretation of the eigenvectors are shown in [Figure 5](#figure--fig: #### Modal Matrix {#modal-matrix} -The modal matrix is an \\(n \times m\\) matrix with columns corresponding to the \\(m\\) system eigenvectors as shown in Eq. +The modal matrix is an \\(n \times m\\) matrix with columns corresponding to the \\(m\\) system eigenvectors as shown in Eq. \eqref{eq:modal\_matrix} \begin{equation} \bm{z}\_m = \begin{bmatrix} @@ -873,7 +873,7 @@ If modes have some damping: \frac{z\_j}{F\_k} = \sum\_{i = 1}^m \frac{z\_{nji} z\_{nki}}{s^2 + 2 \xi\_i \omega\_i s + \omega\_i^2} \label{eq:general\_add\_tf\_damp} \end{equation} -Equations and shows that in general every transfer function is made up of **additive combinations of single degree of freedom systems**, with each system having its DC gain determined by the appropriate eigenvector entry product divided by the square of the eigenvalue, \\(z\_{nji} z\_{nki}/\omega\_i^2\\), and with resonant frequency defined by the eigenvalue \\(\omega\_i\\). +Equations \eqref{eq:general\_add\_tf} and \eqref{eq:general\_add\_tf\_damp} shows that in general every transfer function is made up of **additive combinations of single degree of freedom systems**, with each system having its DC gain determined by the appropriate eigenvector entry product divided by the square of the eigenvalue, \\(z\_{nji} z\_{nki}/\omega\_i^2\\), and with resonant frequency defined by the eigenvalue \\(\omega\_i\\).
diff --git a/content/book/morrison16_groun_shiel.md b/content/book/morrison16_groun_shiel.md index 33ce41e..9227a06 100644 --- a/content/book/morrison16_groun_shiel.md +++ b/content/book/morrison16_groun_shiel.md @@ -182,7 +182,7 @@ And we see that \\(H\\) has units of amperes per meter. The magnetic field of a solenoid is shown in [Figure 7](#figure--fig:morrison16-solenoid). The field intensity inside the solenoid is nearly constant, while outside its intensity falls of rapidly. -Using Ampere's law : +Using Ampere's law \eqref{eq:ampere\_law:} \begin{equation} \oint H dl \approx n I l @@ -244,12 +244,12 @@ V = n^2 A k \mu\_0 \frac{dI}{dt} = L \frac{dI}{dt} where \\(k\\) relates to the geometry of the coil. -Equation states that if \\(V\\) is one volt, then for an inductance of one henry, the current will rise at the rate of one ampere per second. +Equation \eqref{eq:inductance\_coil} states that if \\(V\\) is one volt, then for an inductance of one henry, the current will rise at the rate of one ampere per second. ### The energy stored in an inductance {#the-energy-stored-in-an-inductance} -One way to calculate the work stored in a magnetic field is to use Eq. . +One way to calculate the work stored in a magnetic field is to use Eq. \eqref{eq:inductance\_coil}. The voltage \\(V\\) applied to a coil results in a linearly increasing current. At any time \\(t\\), the power \\(P\\) supplied is equal to \\(VI\\). Power is the rate of change of energy or \\(P = d\bm{E}/dt\\) where \\(\bm{E}\\) is the stored energy in the inductance. diff --git a/content/book/preumont18_vibrat_contr_activ_struc_fourt_edition.md b/content/book/preumont18_vibrat_contr_activ_struc_fourt_edition.md index 1f2022a..295955d 100644 --- a/content/book/preumont18_vibrat_contr_activ_struc_fourt_edition.md +++ b/content/book/preumont18_vibrat_contr_activ_struc_fourt_edition.md @@ -331,9 +331,9 @@ If we left multiply the equation by \\(\Phi^T\\) and we use the orthogonalily re If \\(\Phi^T C \Phi\\) is diagonal, the **damping is said classical or normal**. In this case: \\[ \Phi^T C \Phi = diag(2 \xi\_i \mu\_i \omega\_i) \\] -One can verify that the Rayleigh damping complies with this condition with modal damping ratios \\(\xi\_i = \frac{1}{2} ( \frac{\alpha}{\omega\_i} + \beta\omega\_i )\\). +One can verify that the Rayleigh damping \eqref{eq:rayleigh\_damping} complies with this condition with modal damping ratios \\(\xi\_i = \frac{1}{2} ( \frac{\alpha}{\omega\_i} + \beta\omega\_i )\\). -And we obtain decoupled modal equations . +And we obtain decoupled modal equations \eqref{eq:modal\_eom}.
@@ -366,15 +366,15 @@ Typical values of the modal damping ratio are summarized on table and lies in the change of coordinates. +If one accepts the assumption of classical damping, the only difference between equation \eqref{eq:general\_eom} and \eqref{eq:modal\_eom} lies in the change of coordinates. However, in physical coordinates, the number of degrees of freedom is usually very large. -If a structure is excited in by a band limited excitation, its response is dominated by the modes whose natural frequencies are inside the bandwidth of the excitation and the equation can often be restricted to theses modes. +If a structure is excited in by a band limited excitation, its response is dominated by the modes whose natural frequencies are inside the bandwidth of the excitation and the equation \eqref{eq:modal\_eom} can often be restricted to theses modes. Therefore, the number of degrees of freedom contribution effectively to the response is **reduced drastically** in modal coordinates. #### Dynamic Flexibility Matrix {#dynamic-flexibility-matrix} -If we consider the steady-state response of equation to harmonic excitation \\(f=F e^{j\omega t}\\), the response is also harmonic \\(x = Xe^{j\omega t}\\). The amplitude of \\(F\\) and \\(X\\) is related by: +If we consider the steady-state response of equation \eqref{eq:general\_eom} to harmonic excitation \\(f=F e^{j\omega t}\\), the response is also harmonic \\(x = Xe^{j\omega t}\\). The amplitude of \\(F\\) and \\(X\\) is related by: \\[ X = G(\omega) F \\] Where \\(G(\omega)\\) is called the **Dynamic flexibility Matrix**: @@ -634,7 +634,7 @@ With: - \\(T\_{me}\\) is the transduction coefficient representing the force acting on the mechanical terminals to balance the electromagnetic force induced per unit current input (in \\(\si{\newton\per\ampere}\\)) - \\(Z\_m\\) is the mechanical impedance measured when \\(i=0\\) -Equation shows that the voltage across the electrical terminals of any electromechanical transducer is the sum of a contribution proportional to the current applied and a contribution proportional to the velocity of the mechanical terminals. +Equation \eqref{eq:gen\_trans\_e} shows that the voltage across the electrical terminals of any electromechanical transducer is the sum of a contribution proportional to the current applied and a contribution proportional to the velocity of the mechanical terminals. Thus, if \\(Z\_ei\\) can be measured and substracted from \\(e\\), a signal proportional to the velocity is obtained. To do so, the bridge circuit as shown on [Figure 15](#figure--fig:bridge-circuit) can be used. @@ -706,7 +706,7 @@ With: #### Constitutive Relations of a Discrete Transducer {#constitutive-relations-of-a-discrete-transducer} -The set of equations can be written in a matrix form: +The set of equations \eqref{eq:piezo\_eq} can be written in a matrix form: \begin{equation} \begin{bmatrix}D\\\S\end{bmatrix} @@ -748,7 +748,7 @@ It measures the efficiency of the conversion of the mechanical energy into elect
-If one assumes that all the electrical and mechanical quantities are uniformly distributed in a linear transducer formed by a **stack** (see [Figure 17](#figure--fig:piezo-stack)) of \\(n\\) disks of thickness \\(t\\) and cross section \\(A\\), the global constitutive equations of the transducer are obtained by integrating over the volume of the transducer: +If one assumes that all the electrical and mechanical quantities are uniformly distributed in a linear transducer formed by a **stack** (see [Figure 17](#figure--fig:piezo-stack)) of \\(n\\) disks of thickness \\(t\\) and cross section \\(A\\), the global constitutive equations of the transducer are obtained by integrating \eqref{eq:piezo\_eq\_matrix\_bis} over the volume of the transducer: \begin{equation} \begin{bmatrix}Q\\\\Delta\end{bmatrix} @@ -773,7 +773,7 @@ where {{< figure src="/ox-hugo/preumont18_piezo_stack.png" caption="Figure 17: Piezoelectric linear transducer" >}} -Equation can be inverted to obtain +Equation \eqref{eq:piezo\_stack\_eq} can be inverted to obtain \begin{equation} \begin{bmatrix}V\\\f\end{bmatrix} @@ -801,7 +801,7 @@ The total power delivered to the transducer is the sum of electric power \\(V i\ {{< figure src="/ox-hugo/preumont18_piezo_discrete.png" caption="Figure 18: Discrete Piezoelectric Transducer" >}} -By integrating equation and using the constitutive equations , we obtain the analytical expression of the stored electromechanical energy for the discrete transducer: +By integrating equation \eqref{eq:piezo\_work} and using the constitutive equations \eqref{eq:piezo\_stack\_eq\_inv}, we obtain the analytical expression of the stored electromechanical energy for the discrete transducer: \begin{equation} W\_e(\Delta, Q) = \frac{Q^2}{2 C (1 - k^2)} - \frac{n d\_{33} K\_a}{C(1-k^2)} Q\Delta + \frac{K\_a}{1-k^2}\frac{\Delta^2}{2} @@ -853,7 +853,7 @@ And one can see that \frac{z^2 - p^2}{z^2} = k^2 \end{equation} -Equation constitutes a practical way to determine the electromechanical coupling factor from the poles and zeros of the admittance measurement ([Figure 20](#figure--fig:piezo-admittance-curve)). +Equation \eqref{eq:distance\_p\_z} constitutes a practical way to determine the electromechanical coupling factor from the poles and zeros of the admittance measurement ([Figure 20](#figure--fig:piezo-admittance-curve)). diff --git a/content/book/skogestad07_multiv_feedb_contr.md b/content/book/skogestad07_multiv_feedb_contr.md index 595dadf..eb47342 100644 --- a/content/book/skogestad07_multiv_feedb_contr.md +++ b/content/book/skogestad07_multiv_feedb_contr.md @@ -4447,7 +4447,7 @@ The peak value is close to 6, meaning that even with 6 times less uncertainty, t We here consider the relationship between \\(\mu\\) for RP and the condition number of the plant or of the controller. We consider unstructured multiplicative uncertainty (i.e. \\(\Delta\_I\\) is a full matrix) and performance is measured in terms of the weighted sensitivity. -With \\(N\\) given by , we have: +With \\(N\\) given by \eqref{eq:n\_delta\_structure\_clasic}, we have: \begin{equation\*} \overbrace{\mu\_{\tilde{\Delta}}(N)}^{\text{RP}} \le [ \overbrace{\overline{\sigma}(w\_I T\_I)}^{\text{RS}} + \overbrace{\overline{\sigma}(w\_P S)}^{\text{NP}} ] (1 + \sqrt{k}) @@ -5439,8 +5439,8 @@ for a specified \\(\gamma > \gamma\_\text{min}\\), is given by L &= (1-\gamma^2) I + XZ \end{align} -The Matlab function `coprimeunc` can be used to generate the controller in . -It is important to emphasize that since we can compute \\(\gamma\_\text{min}\\) from we get an explicit solution by solving just two Riccati equations and avoid the \\(\gamma\text{-iteration}\\) needed to solve the general \\(\mathcal{H}\_\infty\\) problem. +The Matlab function `coprimeunc` can be used to generate the controller in \eqref{eq:control\_coprime\_factor}. +It is important to emphasize that since we can compute \\(\gamma\_\text{min}\\) from \eqref{eq:gamma\_min\_coprime} we get an explicit solution by solving just two Riccati equations and avoid the \\(\gamma\text{-iteration}\\) needed to solve the general \\(\mathcal{H}\_\infty\\) problem. #### A Systematic \\(\hinf\\) Loop-Shaping Design Procedure {#a-systematic-hinf-loop-shaping-design-procedure} @@ -6273,7 +6273,7 @@ The selection of \\(u\_2\\) and \\(y\_2\\) for use in the lower-layer control sy Consider the conventional cascade control system in [ 7](#org-target--fig-cascade-extra-meas) where we have additional "secondary" measurements \\(y\_2\\) with no associated control objective, and the objective is to improve the control of \\(y\_1\\) by locally controlling \\(y\_2\\). The idea is that this should reduce the effect of disturbances and uncertainty on \\(y\_1\\). -From , it follows that we should select \\(y\_2\\) and \\(u\_2\\) such that \\(\\|P\_d\\|\\) is small and at least smaller than \\(\\|G\_{d1}\\|\\). +From \eqref{eq:partial\_control}, it follows that we should select \\(y\_2\\) and \\(u\_2\\) such that \\(\\|P\_d\\|\\) is small and at least smaller than \\(\\|G\_{d1}\\|\\). These arguments particularly apply at high frequencies. More precisely, we want the input-output controllability of \\([P\_u\ P\_r]\\) with disturbance model \\(P\_d\\) to be better that of the plant \\([G\_{11}\ G\_{12}]\\) with disturbance model \\(G\_{d1}\\). @@ -6289,7 +6289,7 @@ A set of outputs \\(y\_1\\) may be left uncontrolled only if the effects of all
-To evaluate the feasibility of partial control, one must for each choice of \\(y\_2\\) and \\(u\_2\\), rearrange the system as in and , and compute \\(P\_d\\) using . +To evaluate the feasibility of partial control, one must for each choice of \\(y\_2\\) and \\(u\_2\\), rearrange the system as in \eqref{eq:partial\_control\_partitioning} and \eqref{eq:partial\_control}, and compute \\(P\_d\\) using \eqref{eq:tight\_control\_y2}. #### Measurement Selection for Indirect Control {#measurement-selection-for-indirect-control} @@ -6459,7 +6459,7 @@ We then derive **necessary conditions for stability** which may be used to elimi For decentralized diagonal control, it is desirable that the system can be tuned and operated one loop at a time. Assume therefore that \\(G\\) is stable and each individual loop is stable by itself (\\(\tilde{S}\\) and \\(\tilde{T}\\) are stable). -Using the **spectral radius condition** on the factorized \\(S\\) in , we have that the overall system is stable (\\(S\\) is stable) if +Using the **spectral radius condition** on the factorized \\(S\\) in \eqref{eq:S\_factorization}, we have that the overall system is stable (\\(S\\) is stable) if \begin{equation} \rho(E\tilde{T}(j\omega)) < 1, \forall\omega @@ -6684,7 +6684,7 @@ When a reasonable choice of pairings have been made, one should rearrange \\(G\\ |1 + L\_i| > \max\_{k,j}\\{|\tilde{g}\_{dik}|, |\gamma\_{ij}|\\} \end{equation} - To achieve stability of the individual loops, one must analyze \\(g\_{ii}(s)\\) to ensure that the bandwidth required by is achievable. + To achieve stability of the individual loops, one must analyze \\(g\_{ii}(s)\\) to ensure that the bandwidth required by \eqref{eq:decent\_contr\_one\_loop} is achievable. Note that RHP-zeros in the diagonal elements may limit achievable decentralized control, whereas they may not pose any problems for a multivariable controller. Since with decentralized control, we usually want to use simple controllers, the achievable bandwidth in each loop will be limited by the frequency where \\(\angle g\_{ii}\\) is \\(\SI{-180}{\degree}\\) 6. Check for constraints by considering the elements of \\(G^{-1} G\_d\\) and make sure that they do not exceed one in magnitude within the frequency range where control is needed. @@ -6700,7 +6700,7 @@ When a reasonable choice of pairings have been made, one should rearrange \\(G\\ If the plant is not controllable, then one may consider another choice of pairing and go back to Step 4. If one still cannot find any pairing which are controllable, then one should consider multivariable control. -7. If the chosen pairing is controllable, then tells us how large \\(|L\_i| = |g\_{ii} k\_i|\\) must be. +7. If the chosen pairing is controllable, then \eqref{eq:decent\_contr\_one\_loop} tells us how large \\(|L\_i| = |g\_{ii} k\_i|\\) must be. This can be used as a basis for designing the controller \\(k\_i(s)\\) for loop \\(i\\) @@ -6720,7 +6720,7 @@ Thus sequential design may involve many iterations. #### Conclusion on Decentralized Control {#conclusion-on-decentralized-control} -A number of **conditions for the stability**, e.g. and , and **performance**, e.g. and , of decentralized control systems have been derived. +A number of **conditions for the stability**, e.g. \eqref{eq:decent\_contr\_cond\_stability} and \eqref{eq:decent\_contr\_necessary\_cond\_stability}, and **performance**, e.g. \eqref{eq:decent\_contr\_cond\_perf\_dist} and \eqref{eq:decent\_contr\_cond\_perf\_ref}, of decentralized control systems have been derived. The conditions may be useful in **determining appropriate pairings of inputs and outputs** and the **sequence in which the decentralized controllers should be designed**. diff --git a/content/phdthesis/li01_simul_fault_vibrat_isolat_point.md b/content/phdthesis/li01_simul_fault_vibrat_isolat_point.md index 2d339da..aadf403 100644 --- a/content/phdthesis/li01_simul_fault_vibrat_isolat_point.md +++ b/content/phdthesis/li01_simul_fault_vibrat_isolat_point.md @@ -131,7 +131,7 @@ Define a new input and a new output: u\_1 = J^T f\_m, \quad y = J^{-1} (l - l\_r) \end{equation} -Equation can be rewritten as: +Equation \eqref{eq:hexapod\_eq\_motion} can be rewritten as: \begin{equation} \label{eq:hexapod\_eq\_motion\_decoup\_1} \begin{split} diff --git a/content/phdthesis/rankers98_machin.md b/content/phdthesis/rankers98_machin.md index d2ddb03..3be630b 100644 --- a/content/phdthesis/rankers98_machin.md +++ b/content/phdthesis/rankers98_machin.md @@ -335,17 +335,17 @@ These eigenvectors have the following orthogonality properties, or can always be \phi\_i^T M \phi\_j = 0 \quad (i \neq j) \end{equation} -For \\(i=j\\) the result of the multiplication according to equation yields a non-zero result, which is normally indicated as **modal mass** \\(\mathit{m}\_i\\): +For \\(i=j\\) the result of the multiplication according to equation \eqref{eq:eigenvector\_orthogonality\_mass} yields a non-zero result, which is normally indicated as **modal mass** \\(\mathit{m}\_i\\): \begin{equation} \label{eq:modal\_mass} \phi\_i^T M \phi\_i = \mathit{m}\_i \end{equation} -Because only the direction but not the length of an eigenvector is defined, several scaling methods are used, all based on equation : +Because only the direction but not the length of an eigenvector is defined, several scaling methods are used, all based on equation \eqref{eq:modal\_mass:} -- \\(|\phi\_i| = 1\\). Each eigenvector \\(\phi\_i\\) is scaled such that its length is equal to \\(1\\). The modal mass are then calculated from equation . -- \\(\max(\phi\_i) = 1\\). Each eigenvector \\(\phi\_i\\) is scaled such that its largest element is equation to \\(1\\). The modal mass is then calculated from equation . -- \\(m\_i = 1\\). The modal mass \\(\mathit{m}\_i\\) is set to \\(1\\). The scaling of the mode vector \\(\phi\_i\\) follows from equation . +- \\(|\phi\_i| = 1\\). Each eigenvector \\(\phi\_i\\) is scaled such that its length is equal to \\(1\\). The modal mass are then calculated from equation \eqref{eq:modal\_mass}. +- \\(\max(\phi\_i) = 1\\). Each eigenvector \\(\phi\_i\\) is scaled such that its largest element is equation to \\(1\\). The modal mass is then calculated from equation \eqref{eq:modal\_mass}. +- \\(m\_i = 1\\). The modal mass \\(\mathit{m}\_i\\) is set to \\(1\\). The scaling of the mode vector \\(\phi\_i\\) follows from equation \eqref{eq:modal\_mass}. The orthogonality properties also apply to the stiffness matrix \\(K\\): @@ -450,7 +450,7 @@ In such an analysis one is typically interested in the transfer function between Applying the principle of modal decomposition, any transfer function can be derived by first calculating the behavior of the individual modes, and then **summing all modal contributions**. The contribution of one single mode \\(i\\) to the transfer function \\(x\_l/f\_k\\) can be derived by first considering the response of the modal DoF \\(q\_i\\) to a force vector \\(f\\) with only one non-zero component \\(f\_k\\). -In that case, equation is reduced to: +In that case, equation \eqref{eq:eoq\_modal\_i} is reduced to: \begin{equation} m\_i \ddot{q}\_i(t) + k\_i q\_i(t) = \phi\_{ik} f\_k(t) @@ -481,7 +481,7 @@ The overall transfer function can be found by summation of the individual modal ### Graphical Representation {#graphical-representation} -Due to the equivalence with the differential equations of a single mass spring system, equation is often represented by a single mass spring system on which a force \\(f^\prime = \phi\_i^T f\\) acts. +Due to the equivalence with the differential equations of a single mass spring system, equation \eqref{eq:eoq\_modal\_i} is often represented by a single mass spring system on which a force \\(f^\prime = \phi\_i^T f\\) acts. However, this representation implies an important loss of information because it neglects all information about the mode-shape vector. Consider the system in [Figure 8](#figure--fig:rankers98-mode-trad-representation) for which the three mode shapes are depicted in the traditional graphical representation. @@ -516,7 +516,7 @@ The resulting moment of inertia \\(J\_i\\) of the i-th modal lever then is: J\_i = \sum\_{k=1}^n m\_k \phi\_{ik}^2 \end{equation} -This result is identical to the modal mass \\(m\_i\\) found with Equation , because the mass matrix \\(M\\) is a diagonal matrix of physical masses \\(m\_k\\), and consequently the expression for the modal mass \\(m\_i\\) yields: +This result is identical to the modal mass \\(m\_i\\) found with Equation \eqref{eq:modal\_mass}, because the mass matrix \\(M\\) is a diagonal matrix of physical masses \\(m\_k\\), and consequently the expression for the modal mass \\(m\_i\\) yields: \begin{equation} m\_i = \phi\_j^T M \phi\_j = \sum\_{k=1}^n m\_k \phi\_{ik}^2 @@ -637,7 +637,7 @@ whereas the modal stiffnesses follow from \\(k\_i = \omega\_i^2 m\_i\\). {{< figure src="/ox-hugo/rankers98_example_2dof_modal.png" caption="Figure 15: Graphical representation of modes and modal parameters of the two mass spring system" >}} -From these results, the effective modal parameters for each mode, and for each individual DoF can be defined using equations and . +From these results, the effective modal parameters for each mode, and for each individual DoF can be defined using equations \eqref{eq:m\_modal\_eff} and \eqref{eq:k\_modal\_eff}. The results are summarized in [Table 2](#table--tab:2dof-example-modal-params-eff). @@ -716,7 +716,7 @@ f\_{\text{new},i}(\Delta k) &= \frac{1}{2\pi}\sqrt{\frac{k\_{\text{eff},i} + \De Let's use the two mass spring system in [Figure 14](#figure--fig:rankers98-example-2dof) as an example. In order to analyze the effect of an extra mass at \\(x\_2\\), the effective modal mass at that DoF needs to be known for both modes (see [Table 2](#table--tab:2dof-example-modal-params-eff)). -Then using equation , one can estimate the effect of an extra mass \\(\Delta m = 1 kg\\) added to \\(m\_2\\). +Then using equation \eqref{eq:sensitivity\_add\_m}, one can estimate the effect of an extra mass \\(\Delta m = 1 kg\\) added to \\(m\_2\\). To estimate the influence of extra stiffness between the two DoF, one needs to calculate the effective modal stiffness that corresponds to the relative motion between \\(x\_2\\) and \\(x\_1\\). This can be graphically done as shown in [Figure 19](#figure--fig:rankers98-example-sensitivity-2dof): @@ -726,7 +726,7 @@ k\_{\text{eff},1,(2-1)} &= 0.46 \cdot 10^7 / 0.07^2 = 93.9 \cdot 10^7 N/m \\\\ k\_{\text{eff},2,(2-1)} &= 1.23 \cdot 10^7 / 1.1^2 = 1.0 \cdot 10^7 N/m \end{align} -And using equation , the effect of additional stiffness on the frequency of the two modes can be computed. +And using equation \eqref{eq:sensitivity\_add\_m}, the effect of additional stiffness on the frequency of the two modes can be computed. The results are summarized in [Table 3](#table--tab:example-sensitivity-2dof-results). @@ -860,13 +860,13 @@ Let's introduce a variable \\(\alpha\\), which **relates the high-frequency cont \alpha = \frac{\frac{\phi\_{i,\text{servo}} \phi\_{i,\text{force}}}{m\_i}}{\frac{1}{m}} \end{equation} -which simplifies equation to: +which simplifies equation \eqref{eq:effect\_one\_mode} to: \begin{equation} \label{eq:effect\_one\_mode\_simplified} \boxed{\frac{x\_{\text{servo}}}{F\_{\text{servo}}} = \frac{1}{ms^2} + \frac{\alpha}{m s^2 + m \omega\_i^2}} \end{equation} -Equation will be the basis for the discussion of the various patterns that can be observe in the frequency response functions and the effect of resonances on servo stability. +Equation \eqref{eq:effect\_one\_mode\_simplified} will be the basis for the discussion of the various patterns that can be observe in the frequency response functions and the effect of resonances on servo stability. Three different types of intersection pattern can be found in the amplitude plot as shown in [Figure 25](#figure--fig:rankers98-frf-effect-alpha). Depending on the absolute value of \\(\alpha\\) one can observe: diff --git a/content/posts/test.md b/content/posts/test.md index 826e32a..4d92f3e 100644 --- a/content/posts/test.md +++ b/content/posts/test.md @@ -87,13 +87,13 @@ Here is some inline mathematics: \\(z = 2\\). Unumbered equation: \\[ F(x) = \int\_0^x f(t) dt \\] -Using the `equation` environment in Eq. . +Using the `equation` environment in Eq. \eqref{eq:numbered}. \begin{equation} \label{eq:numbered} F(s) = \int\_0^\infty f(t) e^{-st} dt \end{equation} -Using the `align` environment Equations and . +Using the `align` environment Equations \eqref{eq:align\_1} and \eqref{eq:align\_2}. \begin{align} \mathcal{F}(a) &= \frac{1}{2\pi i}\oint\_\gamma \frac{f(z)}{z - a}\\,dz \label{eq:align\_1} \\\\ diff --git a/content/posts/test2.md b/content/posts/test2.md index 41bed74..c4bdc71 100644 --- a/content/posts/test2.md +++ b/content/posts/test2.md @@ -2,7 +2,7 @@ title = "Second Blog Post" author = ["Dehaeze Thomas"] date = 2021-05-01T00:00:00+02:00 -lastmod = 2024-12-17T15:36:59+01:00 +lastmod = 2024-12-17T16:38:28+01:00 tags = ["hugo", "org"] categories = ["emacs"] draft = false diff --git a/content/zettels/test.md b/content/zettels/test.md index e885b24..b523be2 100644 --- a/content/zettels/test.md +++ b/content/zettels/test.md @@ -465,12 +465,12 @@ Numbering can be continued by using `+n` option as shown below. ```
Code Snippet 4: - Tikz code that is used to generate 4 + Tikz code that is used to generate 4
-{{< figure src="/ox-hugo/test_general_control_names.png" caption="Figure 3: General Control Configuration. With some mathematics in the caption: \\(z = 2\\)" >}} +{{< figure src="/ox-hugo/test_general_control_names.png" caption="Figure 3: General Control Configuration. With some mathematics in the caption: \\(z^2, \sqrt{y}\\)" >}} ```md #+name: fig:general_control_names