Analysis and Algorithm
|
|
- Jody Greene
- 6 years ago
- Views:
Transcription
1 Analysis and Algorith Developers and sponsors assue no responsibility for the use of MIDAS Faily Progra (MIDAS/GeoX, MIDAS/Civil, MIDAS/FX+, MIDAS/Abutent, MIDAS/Pier, MIDAS/Deck, MIDAS/GS, MIDAS/Gen, MIDAS/ADS, MIDAS/SDS, MIDAS/Set ; hereinafter referred to as MIDAS package ) or for the accuracy or validity of any results obtain fro the MIDAS package. Developers and sponsors shall not be liable for loss of profit, loss of business, or financial loss which ay be caused directly or indirectly by the MIDAS package, when used for any purpose or use, due to any defect or deficiency therein.
2 Developers and distributors assue no responsibility for the use of MIDAS Faily Progra (idas Civil, idas FEA, idas FX+, idas Gen, idas Drawing, idas SDS, idas GS, SoilWorks, idas NFX ; hereinafter referred to as MIDAS package ) or for the accuracy or validity of any results obtained fro the MIDAS package. Developers and distributors shall not be liable for loss of profit, loss of business, or financial loss which ay be caused directly or indirectly by the MIDAS package, when used for any purpose or use, due to any defect or deficiency therein. Accordingly, the user is encouraged to fully understand the bases of the progra and becoe failiar with the users anuals. he user shall also independently verify the results produced by the progra.
3 CONENS. Structural Eleents. Overview russ Eleent Bea Eleent Shell Eleent Plane Stress Eleent Plane Strain Eleent Axisyetric Eleent Solid Eleent Spring Rigid Link Reinforceent Eleent. Introduction Reinforceent ypes Linear Analysis Reinforceent in Plane Strain Eleent Reinforceent in Axisyetric Eleent Reinforceent in Plane Stress Eleent Reinforceent in Solid Eleent Reinforceent in Plate Eleent Prestress of Reinforceent i
4 CONENS 3. Interface Eleent 3. Introduction Coordinate Syste and Relative Displaceent Point Interface Eleent Line Interface Eleent Surface Interface Eleent Finite Eleent Forulation Interface Eleent Output Linear Analysis 4. Linear Static Linear Static with Nonlinear Eleents Modal Analysis 5. Introduction Lanczos Iteration Method Subspace Iteration Method Optional Paraeters ie History Analysis 6. ie History Analysis Modal Superposition Method Direct Integration Method Optional Paraeters Cautionary Notes ii
5 CONENS 7. Response Spectru 7. Response Spectru Analysis Spectru Function Linear Buckling Analysis 8. Introduction Geoetric Stiffness Critical Load Factor Extraction Optial Paraeter Load and Boundary 9. Constraint Condition of Degrees of Freedo Skewed (Inclined) Support Condition Constraint Equation Nodal Load Eleent Pressure Load Body Force Prescribed Displaceent Construction Stage Analysis 0. Introduction Coposition of Construction Stages ie Dependent Material Properties iii
6 CONENS. Heat of Hydration. Introduction Heat ransfer Analysis heral Stress Analysis Hear of Hydration Analysis Considering Construction Stages Results of Hear of Hydration Analysis Material Models. Introduction Yield Criteria otal Strain Crack 3. Introduction Basic Properties Loading and Unloading Crack Strain ransforation Stiffness Matrix Copression Models ension Models Shear Models Lateral Influence Interface Nonlinear Behaviors Analysis 4. Introduction Discrete Cracking Crack Dilatancy Bond Slip iv
7 CONENS 4.5 Coulob Friction Cobined Cracking-Shearing-Crushing Geoetric Nonlinear 5. Introduction russ Eleent Plane Stress Eleent Plate Eleent Plane Strain Eleent Axisyetric Eleent Solid Eleent Iteration Method 6. Introduction Initial Stiffness Method Newton-Raphson Method Arc-Length Method Displaceent Control Convergence Criteria Auto-Switching Equation Solver 7. Introduction Direct Method Iterative Method Solver Characteristic v
8 CONENS 8. Contact 8. Introduction Contact Search Function Fatigue 9. Introduction Load Cycles Mean Stress Effects Modifying Factors Rainflow-Counting Algorith Fatigue Analysis Procedure CFD Analysis (Coputational Fluid Dynaic) 0. Introduction RANS equation and turbulence odel Spatial discretization Steady flow Unsteady flow Nuerical stability Coputational fluid analysis results vi
9 Eigenvalue analysis for analyzing the dynaic behavior of structures is also referred to as free vibration analysis. For a syste without daping and with no excitation, the otion equation is the nd order linear differential equation (5.). Mu&& ( t) + Ku( t) = 0 (5.) where, K : Stiffness atrix of structure M : Mass atrix of structure u(t) : Displaceent vector of structure ü(t) : Acceleration vector of structure If the displaceent vector u is assued to be a linear cobination of ode-shape-vectors, defined by the ode shape atrixφ, and the cobination factors for the selected odes are defined by a vector of tie-functions Y ( t), we can replace the displaceent vector u= ΦY ( t). By substituting the expression for u in equation (5.) the following equation is obtained: idasfea MΦY+KΦY && = 0 (5.) he tie function, Y( t) is defined as. { n } Y( t) = y ( t) L y ( t) L y ( t) (5.3) where n is the total nuber of degrees of freedo in the syste. When we assue that the cobination factors y (t) are haronic functions in tie y ( ) t : cos( wt+ b)
10 Analysis and Algorith Manual t he second derivative y (t) with tie of the haronic function can be written as inverse ultiplication with a constant factor λ = ω of the original function y (t) : t y ( t) =- y ( t) l (5.3.) Subsequently we can use this assuption to transfor the equation (5.3) into, (- MΦΛ+ KΦY ) = 0 (5.4) where the atrices Λ and Φ are forulated as below. idasfea Λ él ê ê ê ù ú ú ú, = l ê ú ê ê ë O O ú l ú nû [ f f f ] n l = w (5.5) Φ= L L (5.6) he equation (5.4) ust be valid for every haronic function, that eans for every y ( t ), it can be transfored into. Kf - l M f = 0 (5.7) Equation (5.7) is an eigenvalue proble, which ust satisfy the condition of equation (5.8) and fro this condition the free vibration odes can be calculated. l K- M = 0 (5.8) Equation (5.8) has n nuber of solutions (eigenvalues) l, l,, l that satisfy the condition n (5.8) and equal to the nuber of degrees of freedo of the finite eleent odel. Usually the eigenvalues are ordered such that λ is the sallest eigenvalue. For each eigenvalue, l a corresponding eigen ode, f exists. Since the ass atrix, M and the stiffness atrix, K in
11 equation (5.7) are syetrical, the eigenvalue, l and the eigen ode, f are real nubers. Because the ass atrix is positive definite, and the stiffness atrix is positive sei-definite l ³ 0. As such the circular frequency, w under the condition of undaped free vibration is a real nuber. A structure vibrates in the eigen ode shapes of f without external excitation, and the velocities are circular frequencies, w ( radian / tie ). he velocities of a structure are expressed in ters of natural frequencies, f ( cycle / tie ) or natural periods, ( tie / cycle ). he relationship betweenw, =, f f and is given below. w f = (5.9) p Generally the eigenvalue, l represents the ratio of the strain energy to the kinetic energy for the -th ode shape, and the odes are referred to as st ode, nd ode,, n-th ode fro the sallest ratio onward. Fig. 5-() shows the fr ee vibration odes of a cantilever fro the st to 3 rd ode. idasfea u u u t t t
12 Analysis and Algorith Manual MIDAS calculates indices of the dynaic properties of a structure such as odal participation factors, effective odal asses and odal direction factors. he values of the directional odal participation factors are calculation by equation (5.0), which are used for response spectru analysis or tie history analysis of a structure subjected to seisic loads. idasfea where, fmx fmy fmz G X =, G =, G = Y Z f Mf f Mf f Mf fmrx fmry fmrz G RX =, G =, G = RY RZ f Mf f Mf f Mf (5.0) : Modal participation factors in the GCS X, Y and Z GX, GY, GZ translations for the -th ode : Modal participation factors in the GCS X, Y and Z GRX, GRY, GRZ rotations for the -th ode : Directional vectors, which retain unit values only for the X, Y, Z translations in the X, Y and Z directions : Directional vectors, which retain unit values only for the RX, RY, RZ rotations in the respective X, Y and Z directions f : th ode shape he directional values of the odal participation asses for each ode are calculated according to equation (5.). Since the calculation includes the signs (positive or negative) of the odes, the values can be 0 depending on the ode shapes. he su of the directional participation asses for all the odes is equal to the total ass of the structure in each corresponding direction. General seisic design codes require that the su of the odal participation asses included in the analysis of a structure in each direction be at least 90% of the total ass. his is intended to include ost of the ajor odes, which influence the analysis results.
13 éf ù éf ù éf ù * ë MXû * *, ë MY = = û, = ë MZû X Y Z f Mf f Mf f Mf M M M éf ù éf ù éf ù * ë MRXû * *, ë MRY = = û, = ë MRZû RX RY RZ f Mf f Mf f Mf M M M (5.) where, M, M, M * * * X Y Z M, M, M * * * RX RY RZ : Modal participation asses in the GCS X, Y and Z translations for the -th ode : Modal participation asses in the GCS X, Y and Z rotations for the -th ode he directional coefficient for each ode represents the ratio of the directional participation ass to the total participation ass for the corresponding ode. MIDAS uses the Lanczos and subspace iteration ethods which are suitable for analyzing large systes of eigenvalue probles such as equation (5.7). idasfea
14 Analysis and Algorith Manual he Lanczos iteration ethod finds approxiate eigenvalues using a tri-diagonal atrix,, which can be found by defining the Krylov subspace span{ V, V,..., V k }. In order to effectively apply the Lanczos iteration ethod to the free vibration eigenvalue proble as expressed in equation (5.7), the eigenvalue, l ust be replaced by l = s+ / q. his is referred to as a shift-invert technique. s is the expected first eigenvalue. he Lanczos iteration calculation process applying the shift-invert technique is described in below. k idasfea - Assue an initial value V for the block vector in case of first iterative calculation. - Multiply the ass atrix Uk = MV k - Solve linear siultaneous equations ( K- sm) Wk = U k - Orthogonalize W k - Calculate atrix C k - Orthogonalize * W k - Noralize the block vector W = W - V B * k k k- k- C = V MW * k k k W = W - V C ** * k k k k W = V B ** k k+ k MIDAS uses the block vector, V for effectively calculating eigenvalues. his is referred to as k the Block Lanczos ethod. he block tri-diagonal atrix,, which takes place in the process of iterative calculations above, is as follows: k éc B ù ê ú êb C O ú k = ê O O O ú ê ú ê O Ck- Bk- ú ê ú ë Bk- Ck û (5.)
15 If k is used to solve the eigenvalue proble, l = s + / q. * * y = q y, * * * k * l can be obtained using * l is an approxiate value of the original eigenvalue proble of Eq. (5.7). If N b is defined to be the block size of V, the size of increases to as uch as k k Nb and converges to l, as the nuber of iterations increases in the process of iterative calculations. * he approxiate value, f of the eigenode, f can be obtained by the equation below, which * converges with l. * l [ ] f = V V V V (5.3) * * 3... k y he convergence criterion for eigenvalues and eigenodes calculated fro the Lanczos iteration ethod is noted below. Kf -l Mf * * * K e Where,. eans -nor. MIDAS uses e -6 =. 0. (5.4) idasfea
16 Analysis and Algorith Manual he Subspace ethod searches for a set of N s orthogonal vectors X k that expand over the subspace Ek that is defined by the eigen odes é ëf f.. f N ù sû For large systes the nuber of N s is uch saller than the nuber of degrees of freedo in the finite eleent odel. he iterative process of calculating the set of vectors X is explained below. k idasfea - Assue a first set of N s vectors X and call this X k - Calculate vector Y k with length N s by solving linear set of equations KYk = MX k - Calculate the projection of stiffness in E K k - Calculate the projection of ass E M k = Y K Y k+ k k = Y MY k+ k k - Solve projected eigenvalue proble of size N s Kk+ Qk+ = Mk+ Qk+ Λ k+ by using classical eigenvalue solution techniques * * * é k+ y y y ù N Q = ë... s û Λ k+ * él ù ê * ú l = ê ú ê O ú ê ú * êë ln ú sû Calculating X k+ by Xk+ = Yk Qk+ and repeating these steps until individually converge to N s eigenvalues and N s eigen odes of the syste * ln and X k, * l l, X k éf f... f ù ë Ns û (5.5) If Xk is coposed of N s vectors, the sizes of Kk and the nuber of iterations in the process of iterative calculations. M k are Ns Ns calculated by the subspace iteration ethod is evaluated by the change ofl., irrespective of he convergence of eigenvalues * l -l ( k+ ) ( k ) ( k+ ) l e (5.6)
17 wo types of ass atrices can be used in an eigenvalue analysis, which are the consistent ass atrix and the luped ass atrix. When the consistent ass atrix is applied, this generally results in higher eigenvalues copared to the theoretical values. Whereas using the luped ass atrix results in saller eigenvalues. It is known that using the consistent ass atrix results in better convergence than using the luped ass atrix depending on the nuber of eleents. he drawback however is that the consistent ass atrix requires ore calculations and eory (see figure 5-() ). idasfea When calculating eigen frequencies using the Lanczos iteration ethod in MIDAS two
18 Analysis and Algorith Manual frequencies, f and f can be specified by the user for defining the range of interest of eigen frequencies. he range of interest of frequencies is used in the calculation to define the expected eigenvalues to s = ( p f ) by the shift-invert technique. Fig. 5-(3) shows the range and sequence of calculating eigenvalues for cobinations of f and f. If N f is defined as the nuber of eigenvalues to be calculated, the nuber of each calculated eigenvalues is deterined based on f and f as, - If f= f : he N f frequencies closest to f are calculated. - If f< f : he N f frequencies closest to [ f, f] are calculated. - If f> f : he N f frequencies closest to f within the range, [ f, f] are calculated. idasfea 4 f = f 3 f f f 3 4 f f f 4 3 f he default values of f and f, in the range of interest are f= f= 0. When eigenvalues are calculated using an iteration ethod, there is a possibility that higher eigenvalues are found before the lower eigenvalues. MIDAS provides functionality, which
19 prevents that frequencies are skipped. When applying the Lanczos iteration ethod the vector W k is calculated by solving the following set of equations: ( K- sm) Wk = U k (5.7) In order to solve the equations, the atrix is decoposed by the atrices nuber of negative values in the diagonal ters in the atrix, D is calculated. n LDL and the he nuber represents the nuber of eigenvalues, l that are saller thans. For exaple, if s = l + d and the atrix is decoposed, we can find eigenvalues that are saller than nuber of eigenvalues saller than ln f ln f N f + d. If the + d is greater than N f, the process of the Lanczos iteration continues until all the oitted eigenvalues are converged. At this point, d should satisfy the relationship below. l < l + d < l + (5.8) N N N f f f his condition to δ is called the Stur-sequence check. When MIDAS perfors a Stur sequence check, the size of the Lanczos block vector, V is set to N = 7 to speed up the calculations. k For liiting the coputational effort MIDAS also provides an option to skip the Stur sequence check. b idasfea When a free-vibration analysis of a structure with insufficient boundary conditions is perfored, K will be a singular atrix and the calculations cannot be progressed withs = 0. In such a case, the shift-invert technique is applied with a negative value fors. In this case the calculated eigenvalues will include a 0 for every rigid ode. MIDAS predicts the expected first ode frequency, l using the diagonal ters of the atrices K and M and setss =- l.
20 Analysis and Algorith Manual he subspace iteration ethod in MIDAS allows the user to set the size, k N s of the subspace, X, the nuber of axiu iterations, N and the convergence criterion, e in Eq. (5.6). he actual size of subspace, s s _ 0 f f I N s used for the calculations is as follows: N = ax{ N,in( N, N + 8)} (5.9) Where, N s _ 0 is the user input N s value. he default value for the axiu nuber of iterations is N I = 30. he default value for the convergence criterion is e -6 =.0 0. idasfea
21 When a structure is subjected to a dynaic load, the behavior is described by the dynaic otion equation given by (6.). where, Mu&& ( t) + Cu& ( t) + Ku( t) = F( t) (6.) In case there is no external excitation F( t ) = 0 and the daping atrix C = 0, the equation (6.) will reduce to the free vibration equation (5.). However, if F( t) is an external force (or external displaceent, velocity, acceleration, etc.), which varies with tie, the equation (6.) becoes a forced vibration analysis proble. For this type of probles MIDAS provides a odal superposition ethod as well as a direct integration ethod.. M : Mass atrix C : Daping atrix K : Stiffness atrix F ( t) : Dynaic load vector u ( t), u& ( t) & u&& ( t) : Displaceent, Velocity & Acceleration vector idasfea
22 Analysis and Algorith Manual he odal superposition ethod akes use of the orthogonality of eigen odes and decoposes the equation (6.) into independent ode equations. his ethod assues that the daping atrix is coposed of a linear cobination of the ass atrix and the stiffness atrix. he displaceent vector is obtained by a linear cobination of a selected nuber n of orthogonal eigen odes as expressed in equation (6.). n u( t) FY( t) fy (6.) = =å i= j By substituting equation (6.) into equation (6.), we get i idasfea MΦY&& ( t) + CΦY& ( t) + KΦY( t) = F( t) (6.3) Multiplying equation (6.3) byf (-th ode shape), the equation (6.4) is obtained. f MΦY&& ( t) + f CΦY& ( t) + f KΦY( t) = f F( t) (6.4) he ass and stiffness atrices can be orthogonalized to eigen odes as, f Mf = d i j ij ij f Kf = k d i j ij ij (6.5) where, d : Krönecker Delta ( i= j ; d =, i¹ j ; d = 0 ) ij ij Using the above orthogonality, the dynaic otion equation (6.3) can be decoposed into n independent equations. f Mf Y&& ( t) + f Cf Y& ( t) + f Kf Y( t) = f F( t) ( =,,...,n) (6.6) ij he above n independent decoposed equations thus becoe dynaic otion equations for a single degree of freedo in a general coordinate syste.
23 Subsequently we rearrange equation (6.6) where, ( t) Y && (t ) + xw Y & f F (t ) + wy (t ) = (6.7) f Mf f Cf x w = f Mf w f Kf = f Mf x : th ode daping ratio w : th ode natural frequency q ( t), q& ( t) & q&& ( t) : th ode general displaceent, velocity and acceleration Displaceent response in the general coordinate syste is obtained by equation (6.8). -x (0) (0) wté xwq + q ù q( t) = e êq(0)coswdt+ sinwdtú ë wd û t -xw ( t-t ) + P ( t ) e sin wd( t t ) dt w ò - 0 D (6.8) idasfea where, w = w - x D he displaceent response of a structure is obtained by substituting the general displaceent of each ode, which is obtained by the single degree of freedo equation (6.8), into the syste equation (6.). he accuracy of the displaceent response in the odal superposition ethod depends on the selected odes used for analysis. he odal superposition ethod is ost widely used by structural analysis progras, and it is effective in linear dynaic analyses of large structures. However, this ethod cannot be used for nonlinear dynaic analysis; nor can it be used in case daping devices are included and their properties cannot be assued on the basis of a linear cobination of stiffness and ass.
24 Analysis and Algorith Manual Using the direct integration ethod is one in which the total analysis tie range is sub-divided into a nuber of finite tie steps, and nuerical integration of the dynaic equilibriu equation is perfored at each tie step. his ethod can be applied to systes reflecting nonlinearity of stiffness and/or daping. Since the direct integration ethod evaluates the dynaic otion equation at every tie step, the analysis tie increases with the nuber of tie steps. idasfea A variety of ethods can be applied for nuerical integration. MIDAS uses the average acceleration ethod of the Newark-β ethod in which the acceleration u&& ( t) in the tie range t t t + i < < i is assued to be constant at the average of i u&& and u&& i+ in Eq. (6.9). u&& i + u&& i+ u&& ( t) = = const. (6.9) Consequently, the velocity and displaceent at t = t i+ are expressed as, u&& + u&& u& i i+ i+ = u& i + Dt (6.0) u&& + u&& ui+ = ui + uid t+ Dt 4 & i i+ (6.) Expressing the equations (6.0) and (6.) with the integration variables of the Newark-β ethod, ( g) u& = u& + - D tu&& + gdtu&& (6.) i+ i i i+ i i t æ ö u + = u +D u& i + ç -b D t u&& i + bdt u&& i+ (6.3) è ø
25 where, b =0.5, g =0.5. Rearranging the equation (6.3) gives the following expression for the acceleration at the end of the tie-step, ì æ ö ü u&& i = i i t i b t + íu - u -D u& - - D i b t + ç u&& ý (6.4) D î è ø þ Substituting equation (6.4) into equation (6.) and rearranging it, gives the following expression for the velocity at the end of the tie-step, g g æ g ö æ g ö u& t i + = u b t i + - u + b t i ç - u& + - D b i ç u&& (6.5) D D è ø è bø i he equations (6.4) and (6.5) are substituted into the dynaic otion equation, which is rearranged for the displaceent response ui+ at the end of the increent as follows: idasfea æ g ö ç M+ C+K u èbdt bdt ø i+ ìï æ ö üï ìï g æg ö æ g ö üï = F+ Mí u i + u& i + ç - u&& iý+ Cí ui + ç - u& i + ç - Dtu&& iý ïî bdt bdt è b ø ïþ ïî bdt èb ø è b ø ïþ (6.6) Substituting the displaceent u i+ at the tie t i + as defined in equation (6.6) into equations (6.4) and (6.5), the velocity and acceleration at the end of the increent as function of status at the beginning of the increent only can be obtained. Rayleigh daping expressed in equation (6.7) is used for daping in the direct integration
26 Analysis and Algorith Manual ethod. Rayleigh Daping: C = a0m+ ak (6.7) where, a 0 & a : Proportional constants for ass and stiffness for daping calculation By substituting the equation (6.7) into equation (6.6), the dynaic otion equation is expressed as, idasfea ìï æ ag ö æ ag ö üï íç + M+ ç + Kýu ïî èbdt bdtø èbdt ø ïþ ìï æ ö üï =F+ Mí u i + u& i + ç - u&& i + adý+ akd ïî bdt bdt è b ø ïþ g æg ö æ g ö where, D = ui + ç - u& i + ç - Dtu&& i bdt èb ø è b ø i+ (6.8) he solutions for the dynaic otion equation that are found through tie history analysis are the relative displaceent u ( t), the relative velocity u& ( t) and the relative acceleration u&& ( t). When a structure is subjected to dynaic loads such as ground acceleration, the absolute response of the structure is obtained by adding the relative response and the ground response as expressed in equation (6.9). u&& + u&& : Absolute acceleration g, i+ i+ u& + u& : Absolute velocity (6.9) g, i+ i+ u + u : Absolute displaceent g, i+ i+
27 where, u&& g, i+, u& g, i+ & u g, i+ are ground acceleration, velocity & displaceent respectively. In MIDAS, ground velocity and displaceent due to ground acceleration are calculated by using the linear acceleration ethod (equation 6.0), which are applied to calculating the absolute response in the odal superposition and direct integration ethods. u&& - u&& u& g, i+ = u& g, i +D tu&& g, i + Dt Dt u&& u u u u 6 g, i+ g, i - u&& g, i+ g, i 3 g, i+ = g, i +D t& g, i + D t && g, i + Dt Dt (6.0) idasfea
28 Analysis and Algorith Manual In MIDAS, the following daping ethods are used depending on the ethod of dynaic analysis. Selection of daping for tie history analysis by response spectru and odal superposition: - Modal - Mass & Stiffness Proportional (Rayleigh daping) idasfea Selection of daping for tie history analysis by direct integration: - Mass & Stiffness Proportional (Rayleigh daping) he daping atrix in Rayleigh Daping is coposed of the linear cobination of the ass atrix and the stiffness atrix of a structure as shown in Fig. (6-b). If the daping integer and the natural frequency w r for the r th ode and the daping integerxs and the natural frequency w s for the s th ode are given, the Rayleigh daping atrix is expressed below. Notice that the r th and s th odes represent the two ain odes of the structure. x r C= a M+ a K (6.) 0 æ a ö 0 xi = ç + a wi èwi ø (6.) where, ( ) w w x w - x w a =, r s r s s r 0 ( ws -wr ) a = ( x w - x w ) s s r r ( ws -wr )
29 x i Mass Proportional C= a 0 M a0 xi = w i C= a K a wi x i = Stiffness Proportional w w w3 w4 Natural frequencies w i idasfea For odal daping, the user directly defines the daping ratio for each ode. Modal daping can be used for tie history analysis, by response spectru analysis and the odal superposition ethod. When the response spectru analysis and odal superposition ethod are used, the kineatic equation of the structure is decoposed by odes to which the odal daping ratios defined by the user are applied.
30 Analysis and Algorith Manual he accuracy of analysis depends greatly on the tie interval used in analysis. Iproper tie interval ay result in inaccurate solutions. In particular, the size of the tie interval is closely related to the axiu frequency of oveent of the structure. In general, one-tenth of the highest odal period under consideration is a reasonable value for the tie interval. In addition, the tie interval should be saller than that of the applied load. A dynaic load needs to sufficiently depict the change in the total loads. MIDAS linearly interpolates the excitation loads. idasfea
31 Response spectru analyses are generally carried out for seisic designs using the design spectra defined in design standards. Response spectru analysis assues the response of a ulti-degreeof-freedo (MDOF) syste as a cobination of ultiple single-degree-of-freedo (SDOF) systes. A response spectru defines the peak values of responses corresponding to and varying with natural periods (or frequencies) of vibration that have been prepared through a nuerical integration process. Displaceents, velocities and accelerations for the basis of a spectru. o predict the peak design response values, the axiu response for each ode is obtained first and then cobined by an appropriate ethod. Equation (7.) shows the dynaic otion equation for a structure subjected to a ground otion used in a response spectru analysis. M[ u&& ( t) + ru&& ( t)] + Cu& ( t) + Ku( t) = 0 g Mu&& ( t) + Cu& ( t) + Ku( t) =-Mru&& ( t) g (7.) idasfea where, M C K : Mass atrix : Daping atrix : Stiffness atrix r : Directional vector of ground acceleration u&& g ( t ) : ie history of ground acceleration u( t ), u& ( t ) & u&& ( t ) : Relative displaceent, velocity and acceleration If the displaceent, u( t) is expressed in ters of a cobination of odal displaceents using the eigen ode shapesφ, obtained fro undaped, free vibration analysis, equation (7.) is
32 Analysis and Algorith Manual applicable. u( t) = Φy( t ) (7.) Now we substitute equation. (7.) into equation (7.) and ultiply both sides byφ, resulting in Φ MΦy &&( t) + Φ CΦy& ( t) + Φ KΦy( t) =-Φ Mr u&& ( t ) (7.3) g he relationship below is found due to orthogonality of the eigen odes. i ( or or ) 0 ( i j) f M C K f = ¹ (7.4) j Accordingly, if the diensionless eigen ode shape ( Φ MΦ = ) for ass is applied to equation (7.3), we obtain the independent siultaneous differential equations for each ode. idasfea é ù éxw ù éw ù ê ú ê ú ê ú ê O ú ê O ú ê O ú ê ú && y t + ê ú y& t + ê ú y t =-Φ Mru&& t ê ú ê ú ê ú ê O ú ê O ú ê O ú ê ú ê x w ú ê w ú ë û ë n nû ë nû ( ) x w ( ) w ( ) g ( ) (7.5) And rearranging the expression for the -th ode in equation (7.5) results in. && y t x w y& t w y t u&& t ( ) + ( ) + ( ) =-G g ( ) G = f Mr (7.6) he odal participation factor, G in equation (7.6) is defined by the ultiplication of the diensionless ode shape for ass φ, the ass atrix M and the directional vector of ground acceleration r. he directional vector of ground acceleration contains unit value only for the degree of freedo in the direction of ground acceleration and all other coponents are equal to zero. he solution to the dynaic equilibriu equation of a structure under ground acceleration
33 action is obtained by solving the n equations of (7.6) and then cobining the in the sae way as equation (7.7). u( t) = Φy( t), u& ( t) = Φy& ( t), u&& ( t) = Φy && ( t ) (7.7) he concept of the dynaic equilibriu equation subjected to a ground acceleration action in response spectru analysis follows the equations (7.) to (7.7). In response spectru analysis, spectru functions are used to obtain the results for each ode rather than defining the seisic acceleration as a specific function. Generally, the spectru function eans the axiu value for each period in equation (7.5). he spectru function is defined in various design codes and specifications, which consider the probabilities of seisic accelerations, characteristics of regions and iportance of structures. he solutions to a single d.o.f syste corresponding to the -th ode like equation (7.6) are found by ultiplying the displaceent, velocity and acceleration obtained fro the spectru function by the corresponding odal participation factors as equation (7.8). Sa Sa Sd =, S = v w w y =G S, y& =G S, && y =G S d v a (7.8) idasfea he results for each ode are calculated by ultiplying the results of equation (7.7) by the ode shapes as in equation (7.9). u = f G S, u& = f G S, u&& = f G S (7.9) d v a where, S d : Spectral displaceent of th ode S v : Spectral velocity for of th ode S a : Spectral acceleration of th ode Since the analysis results for each ode pertain to only the axiu values, it is not possible to perfor linear cobinations as done in tie history analysis. herefore, the final results of a
34 Analysis and Algorith Manual response spectru analysis are obtained by a odal cobination of the analysis results for each ode in equation (7.9). MIDAS enables us to analyze response spectru in any direction on the X-Y plane and in the Z- direction in the global coordinate syste. he user can choose one of the 3 ethods for odal cobination, which are ABS (Absolute Su), SRSS (Square Root of the Su of the Squares) and CQC (Coplete Quadratic Cobination) ethods. ABS (ABsolute Su) Rax = R + R + + Rn (7.0) idasfea SRSS (Square Root of the Su of the Squares) Rax = é ër + R + + R ù n û (7.) CQC (Coplete Quadratic Cobination) é ù ax =êåå N N R Ririj R jú (7.) ë i= j= û 3 i j i + rijx rij 8 xx ( x ) rij = (- r ) + 4 r ( + r ) + 4( + ) r ij xx i j ij ij xi x j ij 8 x ( + r ) r r = ( x = x = x ) 3 ij ij ij i j (- rij ) + 4 x rij ( + rij ) wi r = w > w w ij j i j 0 rij rij = ( i= j )
35 where, R ax : Peak response R i r ij w i,w j : Peak response (spectru value) of i th ode : Eigenvalue ratio of j th ode to i th ode : Eigenvalues of i th and j th odes, x x i j : Daping ratios of i th and j th odes In equation (7.), when i = j, then ρ ij = regardless of the daping ratio ( x i, x j ). If the daping ratio becoes zero (0), both CQC and SRSS ethods produce identical results. he ABS ethod produces the largest cobination values aong the three ethods. he SRSS ethod has been widely used in the past, but it tends to overestiate or underestiate the cobination results in the cases where the values of natural frequencies are close to one another. As a result, the use of the CQC ethod is increasing recently as it accounts for probabilistic inter-relations between the odes. Exaple idasfea If we now copare the natural frequencies and displaceents for each ode for a structure having 3 DOF s with a daping ratio of 0.05, the results fro the applications of SRSS and CQC are as follows: Natural frequencies w = 0.46, w = 0.5, w =.4 3 Response spectru value for each ode: for j-th ode) D ij (displaceent coponent of i-th degree of freedo é ù é ê ú ëd ijù= û ê - - ú êë úû
36 Analysis and Algorith Manual SRSS ethod results é0.04ù Rax = ér R R ù ê ú ë + + û = ê ú êë 0.05úû CQC ethod results r = r= r3 = r3= r3= r3 = idasfea é0.046ù Rax = ér R R3 rr R r3r R3 r3rr ù ê ú ë û = ê ú êë 0.053úû Coparing the two sets of displaceents for each degree of freedo, we note that the SRSS ethod underestiates the agnitude for the first degree of freedo but overestiates the value for the second degree of freedo relative to those obtained by CQC. hus, the SRSS ethod should be used with care when natural frequencies are close to one another.
37 idasfea S7- S6 Sx = ( x - 6) + S In response spectru analysis, results for each ode are calculated by using spectru functions. Generally, the spectru functions are coposed of the axiu values of tie history analysis results obtained fro equation (7.6). Once the daping ratiox and the seisic acceleration varying with tie u&& g ( t) are defined, solutions to equation (7.6) can be found according to the natural periods of the structure. Fig. 7-() shows how the spectru functions are deterined using a plot of the results of displaceent, velocity and acceleration vertically against the periods of the
38 Analysis and Algorith Manual structure horizontally. Spectru functions used in response spectru analysis are generally provided by various design codes. MIDAS generates spectru functions used in seisic analysis according to a selected design codes by siply entering the dynaic factor, foundation factor, zoning factor, iportance factor, response odification factor, etc. Since linear or log scale interpolation is used to deterine spectru values corresponding to the natural periods of the structure, it is recoended that the data in the region of rapid changes be closely defined. And the range of the spectru function ust cover the range of the axiu and iniu periods calculated fro the eigenvalue analysis. MIDAS uses the axiu or iniu value of the spectru function if the periods of eigenvalue analysis exceed the range of the inputted range of the spectru function. idasfea
Seismic Analysis of Structures by TK Dutta, Civil Department, IIT Delhi, New Delhi.
Seisic Analysis of Structures by K Dutta, Civil Departent, II Delhi, New Delhi. Module 5: Response Spectru Method of Analysis Exercise Probles : 5.8. or the stick odel of a building shear frae shown in
More informationEVALUATION OF A SIMPLIFIED METHOD FOR THE DETERMINATION OF THE NON LINEAR SEISMIC RESPONSE OF RC FRAMES
EVALUATIO OF A SIMPLIFIED METHOD FOR THE DETERMIATIO OF THE O LIEAR SEISMIC RESPOSE OF RC FRAMES 9 Misael REQUEA And A. Gustavo AYALA SUMMARY In this paper a siplified ethod is developed for the evaluation
More informationCOMPONENT MODE SYNTHESIS, FIXED-INTERFACE MODEL Revision A
COMPONEN MODE SYNHESS, FXED-NERFACE MODEL Revision A By o rvine Eail: toirvine@aol.co February, ntroduction Coponent ode synthesis is a ethod for analyzing the dynaic behavior of a syste consisting of
More informationFeature Extraction Techniques
Feature Extraction Techniques Unsupervised Learning II Feature Extraction Unsupervised ethods can also be used to find features which can be useful for categorization. There are unsupervised ethods that
More informationMonitoring and system identification of suspension bridges: An alternative approach
Monitoring and syste identification of suspension bridges: An alternative approach Erdal Şafak Boğaziçi University, Kandilli Observatory and Earthquake Reseach Institute, Istanbul, Turkey Abstract This
More informationSEISMIC FRAGILITY ANALYSIS
9 th ASCE Specialty Conference on Probabilistic Mechanics and Structural Reliability PMC24 SEISMIC FRAGILITY ANALYSIS C. Kafali, Student M. ASCE Cornell University, Ithaca, NY 483 ck22@cornell.edu M. Grigoriu,
More informationLecture #8-3 Oscillations, Simple Harmonic Motion
Lecture #8-3 Oscillations Siple Haronic Motion So far we have considered two basic types of otion: translation and rotation. But these are not the only two types of otion we can observe in every day life.
More information821. Study on analysis method for deepwater TTR coupled vibration of parameter vibration and vortex-induced vibration
81. Study on analysis ethod for deepwater TTR coupled vibration of paraeter vibration and vortex-induced vibration Wu Xue-Min 1, Huang Wei-Ping Shandong Key aboratory of Ocean Engineering, Ocean University
More informationma x = -bv x + F rod.
Notes on Dynaical Systes Dynaics is the study of change. The priary ingredients of a dynaical syste are its state and its rule of change (also soeties called the dynaic). Dynaical systes can be continuous
More informationUsing a De-Convolution Window for Operating Modal Analysis
Using a De-Convolution Window for Operating Modal Analysis Brian Schwarz Vibrant Technology, Inc. Scotts Valley, CA Mark Richardson Vibrant Technology, Inc. Scotts Valley, CA Abstract Operating Modal Analysis
More informationACTIVE VIBRATION CONTROL FOR STRUCTURE HAVING NON- LINEAR BEHAVIOR UNDER EARTHQUAKE EXCITATION
International onference on Earthquae Engineering and Disaster itigation, Jaarta, April 14-15, 8 ATIVE VIBRATION ONTROL FOR TRUTURE HAVING NON- LINEAR BEHAVIOR UNDER EARTHQUAE EXITATION Herlien D. etio
More informationLeast Squares Fitting of Data
Least Squares Fitting of Data David Eberly, Geoetric Tools, Redond WA 98052 https://www.geoetrictools.co/ This work is licensed under the Creative Coons Attribution 4.0 International License. To view a
More informationBasic concept of dynamics 3 (Dynamics of a rigid body)
Vehicle Dynaics (Lecture 3-3) Basic concept of dynaics 3 (Dynaics of a rigid body) Oct. 1, 2015 김성수 Vehicle Dynaics Model q How to describe vehicle otion? Need Reference fraes and Coordinate systes 2 Equations
More informationChapter 6 1-D Continuous Groups
Chapter 6 1-D Continuous Groups Continuous groups consist of group eleents labelled by one or ore continuous variables, say a 1, a 2,, a r, where each variable has a well- defined range. This chapter explores:
More informationCh 12: Variations on Backpropagation
Ch 2: Variations on Backpropagation The basic backpropagation algorith is too slow for ost practical applications. It ay take days or weeks of coputer tie. We deonstrate why the backpropagation algorith
More informationChapter 1: Basics of Vibrations for Simple Mechanical Systems
Chapter 1: Basics of Vibrations for Siple Mechanical Systes Introduction: The fundaentals of Sound and Vibrations are part of the broader field of echanics, with strong connections to classical echanics,
More informationPHYSICS 110A : CLASSICAL MECHANICS MIDTERM EXAM #2
PHYSICS 110A : CLASSICAL MECHANICS MIDTERM EXAM #2 [1] Two blocks connected by a spring of spring constant k are free to slide frictionlessly along a horizontal surface, as shown in Fig. 1. The unstretched
More informationDETECTION OF NONLINEARITY IN VIBRATIONAL SYSTEMS USING THE SECOND TIME DERIVATIVE OF ABSOLUTE ACCELERATION
DETECTION OF NONLINEARITY IN VIBRATIONAL SYSTEMS USING THE SECOND TIME DERIVATIVE OF ABSOLUTE ACCELERATION Masaki WAKUI 1 and Jun IYAMA and Tsuyoshi KOYAMA 3 ABSTRACT This paper shows a criteria to detect
More informationANALYSIS ON RESPONSE OF DYNAMIC SYSTEMS TO PULSE SEQUENCES EXCITATION
The 4 th World Conference on Earthquake Engineering October -7, 8, Beijing, China ANALYSIS ON RESPONSE OF DYNAMIC SYSTEMS TO PULSE SEQUENCES EXCITATION S. Li C.H. Zhai L.L. Xie Ph. D. Student, School of
More informationAnalysis of Impulsive Natural Phenomena through Finite Difference Methods A MATLAB Computational Project-Based Learning
Analysis of Ipulsive Natural Phenoena through Finite Difference Methods A MATLAB Coputational Project-Based Learning Nicholas Kuia, Christopher Chariah, Mechatronics Engineering, Vaughn College of Aeronautics
More informationEnvelope frequency Response Function Analysis of Mechanical Structures with Uncertain Modal Damping Characteristics
Copyright c 2007 Tech Science Press CMES, vol.22, no.2, pp.129-149, 2007 Envelope frequency Response Function Analysis of Mechanical Structures with Uncertain Modal Daping Characteristics D. Moens 1, M.
More informationSupplementary Information for Design of Bending Multi-Layer Electroactive Polymer Actuators
Suppleentary Inforation for Design of Bending Multi-Layer Electroactive Polyer Actuators Bavani Balakrisnan, Alek Nacev, and Elisabeth Sela University of Maryland, College Park, Maryland 074 1 Analytical
More informationNUMERICAL MODELLING OF THE TYRE/ROAD CONTACT
NUMERICAL MODELLING OF THE TYRE/ROAD CONTACT PACS REFERENCE: 43.5.LJ Krister Larsson Departent of Applied Acoustics Chalers University of Technology SE-412 96 Sweden Tel: +46 ()31 772 22 Fax: +46 ()31
More informationModel Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon
Model Fitting CURM Background Material, Fall 014 Dr. Doreen De Leon 1 Introduction Given a set of data points, we often want to fit a selected odel or type to the data (e.g., we suspect an exponential
More informationChaotic Coupled Map Lattices
Chaotic Coupled Map Lattices Author: Dustin Keys Advisors: Dr. Robert Indik, Dr. Kevin Lin 1 Introduction When a syste of chaotic aps is coupled in a way that allows the to share inforation about each
More informationDynamic Response of SDOF System: A Comparative Study Between Newmark s Method and IS1893 Response Spectra
Dynaic Response of SDOF Syste: A Coparative Study Between Newark s Method and IS893 Response Spectra [] P. Manoj Reddy [2] Dr. A. Viala [] Post Graduate Student, [2] Associate Professor [][2] Chaitanya
More informationFrame with 6 DOFs. v m. determining stiffness, k k = F / water tower deflected water tower dynamic response model
CE 533, Fall 2014 Undaped SDOF Oscillator 1 / 6 What is a Single Degree of Freedo Oscillator? The siplest representation of the dynaic response of a civil engineering structure is the single degree of
More informationPhysics 139B Solutions to Homework Set 3 Fall 2009
Physics 139B Solutions to Hoework Set 3 Fall 009 1. Consider a particle of ass attached to a rigid assless rod of fixed length R whose other end is fixed at the origin. The rod is free to rotate about
More information2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all
Lecture 6 Introduction to kinetic theory of plasa waves Introduction to kinetic theory So far we have been odeling plasa dynaics using fluid equations. The assuption has been that the pressure can be either
More informationCOS 424: Interacting with Data. Written Exercises
COS 424: Interacting with Data Hoework #4 Spring 2007 Regression Due: Wednesday, April 18 Written Exercises See the course website for iportant inforation about collaboration and late policies, as well
More informationPhysics 2210 Fall smartphysics 20 Conservation of Angular Momentum 21 Simple Harmonic Motion 11/23/2015
Physics 2210 Fall 2015 sartphysics 20 Conservation of Angular Moentu 21 Siple Haronic Motion 11/23/2015 Exa 4: sartphysics units 14-20 Midter Exa 2: Day: Fri Dec. 04, 2015 Tie: regular class tie Section
More informationAutomated Frequency Domain Decomposition for Operational Modal Analysis
Autoated Frequency Doain Decoposition for Operational Modal Analysis Rune Brincker Departent of Civil Engineering, University of Aalborg, Sohngaardsholsvej 57, DK-9000 Aalborg, Denark Palle Andersen Structural
More informationHarmonic Standing-Wave Excitations of Simply-Supported Isotropic Solid Elastic Circular Cylinders: Exact 3D Linear Elastodynamic Response.
Haronic Standing-Wave Excitations of Siply-Supported Isotropic Solid Elastic Circular Cylinders: Exact 3D inear Elastodynaic Response Jaal Sakhr and Blaine A. Chronik Departent of Physics and Astronoy,
More informationPh 20.3 Numerical Solution of Ordinary Differential Equations
Ph 20.3 Nuerical Solution of Ordinary Differential Equations Due: Week 5 -v20170314- This Assignent So far, your assignents have tried to failiarize you with the hardware and software in the Physics Coputing
More informationAnalysis of ground vibration transmission in high precision equipment by Frequency Based Substructuring
Analysis of ground vibration transission in high precision equipent by Frequency Based Substructuring G. van Schothorst 1, M.A. Boogaard 2, G.W. van der Poel 1, D.J. Rixen 2 1 Philips Innovation Services,
More informationPhysics 215 Winter The Density Matrix
Physics 215 Winter 2018 The Density Matrix The quantu space of states is a Hilbert space H. Any state vector ψ H is a pure state. Since any linear cobination of eleents of H are also an eleent of H, it
More informationEFFECT OF MATERIAL PROPERTIES ON VIBRATIONS OF NONSYMMETRICAL AXIALLY LOADED THIN-WALLED EULER-BERNOULLI BEAMS
Matheatical and Coputational Applications, Vol. 5, No., pp. 96-07, 00. Association for Scientific Research EFFECT OF MATERIAL PROPERTIES ON VIBRATIONS OF NONSYMMETRICAL AXIALLY LOADED THIN-WALLED EULER-BERNOULLI
More informationData-Driven Imaging in Anisotropic Media
18 th World Conference on Non destructive Testing, 16- April 1, Durban, South Africa Data-Driven Iaging in Anisotropic Media Arno VOLKER 1 and Alan HUNTER 1 TNO Stieltjesweg 1, 6 AD, Delft, The Netherlands
More informationQ5 We know that a mass at the end of a spring when displaced will perform simple m harmonic oscillations with a period given by T = 2!
Chapter 4.1 Q1 n oscillation is any otion in which the displaceent of a particle fro a fixed point keeps changing direction and there is a periodicity in the otion i.e. the otion repeats in soe way. In
More informationBlock designs and statistics
Bloc designs and statistics Notes for Math 447 May 3, 2011 The ain paraeters of a bloc design are nuber of varieties v, bloc size, nuber of blocs b. A design is built on a set of v eleents. Each eleent
More informationi ij j ( ) sin cos x y z x x x interchangeably.)
Tensor Operators Michael Fowler,2/3/12 Introduction: Cartesian Vectors and Tensors Physics is full of vectors: x, L, S and so on Classically, a (three-diensional) vector is defined by its properties under
More information2nd Workshop on Joints Modelling Dartington April 2009 Identification of Nonlinear Bolted Lap Joint Parameters using Force State Mapping
Identification of Nonlinear Bolted Lap Joint Paraeters using Force State Mapping International Journal of Solids and Structures, 44 (007) 8087 808 Hassan Jalali, Haed Ahadian and John E Mottershead _ Γ
More information2. A crack which is oblique (Swedish sned ) with respect to the xy coordinate system is to be analysed. TMHL
(Del I, teori; 1 p.) 1. In fracture echanics, the concept of energy release rate is iportant. Fro the fundaental energy balance of a case with possible crack growth, one usually derives the equation where
More information13 Harmonic oscillator revisited: Dirac s approach and introduction to Second Quantization
3 Haronic oscillator revisited: Dirac s approach and introduction to Second Quantization. Dirac cae up with a ore elegant way to solve the haronic oscillator proble. We will now study this approach. The
More informationANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER
IEPC 003-0034 ANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER A. Bober, M. Guelan Asher Space Research Institute, Technion-Israel Institute of Technology, 3000 Haifa, Israel
More informationDESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS *
IJST, Transactions of Mechanical Engineering, Vol. 39, No. M1, pp 89-100 Printed in The Islaic Republic of Iran, 2015 Shira University DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS
More informationDetermining a Function for the Damping Coefficient of a laminated Stack
DOI: 10.435/UB.OVGU-017-093 TECHNISCHE MECHANIK, 37, -5, (017), 161 170 subitted: June 9, 017 Deterining a Function for the Daping Coefficient of a lainated Stack C. Zahalka, K. Ellerann The design of
More informationAn Inverse Interpolation Method Utilizing In-Flight Strain Measurements for Determining Loads and Structural Response of Aerospace Vehicles
An Inverse Interpolation Method Utilizing In-Flight Strain Measureents for Deterining Loads and Structural Response of Aerospace Vehicles S. Shkarayev and R. Krashantisa University of Arizona, Tucson,
More informationIntelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley osig 1 Winter Seester 2018 Lesson 6 27 February 2018 Outline Perceptrons and Support Vector achines Notation...2 Linear odels...3 Lines, Planes
More informationF = 0. x o F = -k x o v = 0 F = 0. F = k x o v = 0 F = 0. x = 0 F = 0. F = -k x 1. PHYSICS 151 Notes for Online Lecture 2.4.
PHYSICS 151 Notes for Online Lecture.4 Springs, Strings, Pulleys, and Connected Objects Hook s Law F = 0 F = -k x 1 x = 0 x = x 1 Let s start with a horizontal spring, resting on a frictionless table.
More informationIn this chapter we will start the discussion on wave phenomena. We will study the following topics:
Chapter 16 Waves I In this chapter we will start the discussion on wave phenoena. We will study the following topics: Types of waves Aplitude, phase, frequency, period, propagation speed of a wave Mechanical
More informationChapter 2: Introduction to Damping in Free and Forced Vibrations
Chapter 2: Introduction to Daping in Free and Forced Vibrations This chapter ainly deals with the effect of daping in two conditions like free and forced excitation of echanical systes. Daping plays an
More informationMulti-Scale/Multi-Resolution: Wavelet Transform
Multi-Scale/Multi-Resolution: Wavelet Transfor Proble with Fourier Fourier analysis -- breaks down a signal into constituent sinusoids of different frequencies. A serious drawback in transforing to the
More informationDynamic analysis of frames with viscoelastic dampers: a comparison of damper models
Structural Engineering and Mechanics, Vol. 41, No. 1 (2012) 113-137 113 Dynaic analysis of fraes with viscoelastic dapers: a coparison of daper odels R. Lewandowski*, A. Bartkowiak a and H. Maciejewski
More informationAn Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period
An Approxiate Model for the Theoretical Prediction of the Velocity... 77 Central European Journal of Energetic Materials, 205, 2(), 77-88 ISSN 2353-843 An Approxiate Model for the Theoretical Prediction
More informationRECOVERY OF A DENSITY FROM THE EIGENVALUES OF A NONHOMOGENEOUS MEMBRANE
Proceedings of ICIPE rd International Conference on Inverse Probles in Engineering: Theory and Practice June -8, 999, Port Ludlow, Washington, USA : RECOVERY OF A DENSITY FROM THE EIGENVALUES OF A NONHOMOGENEOUS
More informationHydro-Elastic Criterion for Practical Design
Hydro-Elastic Criterion for Practical Design Hannes Bogaert ), Mirek Kainski ) ) MARIN, Hydro-Structural Services, Wageningen, Netherlands & Delft University of Technology, Ship Structures Laboratory,
More informationPhysics 140 D100 Midterm Exam 2 Solutions 2017 Nov 10
There are 10 ultiple choice questions. Select the correct answer for each one and ark it on the bubble for on the cover sheet. Each question has only one correct answer. (2 arks each) 1. An inertial reference
More informationTopic 5a Introduction to Curve Fitting & Linear Regression
/7/08 Course Instructor Dr. Rayond C. Rup Oice: A 337 Phone: (95) 747 6958 E ail: rcrup@utep.edu opic 5a Introduction to Curve Fitting & Linear Regression EE 4386/530 Coputational ethods in EE Outline
More informationSimple Schemes of Multi anchored Flexible Walls Dynamic Behavior
6 th International Conference on Earthquake Geotechnical Engineering -4 Noveber 05 Christchurch, New Zealand Siple Schees of Multi anchored Flexible Walls Dynaic Behavior A. D. Garini ABSTRACT Siple schees
More informationTutorial Exercises: Incorporating constraints
Tutorial Exercises: Incorporating constraints 1. A siple pendulu of length l ass is suspended fro a pivot of ass M that is free to slide on a frictionless wire frae in the shape of a parabola y = ax. The
More information2.003 Engineering Dynamics Problem Set 2 Solutions
.003 Engineering Dynaics Proble Set Solutions This proble set is priarily eant to give the student practice in describing otion. This is the subject of kineatics. It is strongly recoended that you study
More informationANALYTICAL INVESTIGATION AND PARAMETRIC STUDY OF LATERAL IMPACT BEHAVIOR OF PRESSURIZED PIPELINES AND INFLUENCE OF INTERNAL PRESSURE
DRAFT Proceedings of the ASME 014 International Mechanical Engineering Congress & Exposition IMECE014 Noveber 14-0, 014, Montreal, Quebec, Canada IMECE014-36371 ANALYTICAL INVESTIGATION AND PARAMETRIC
More informationChapter 11 Simple Harmonic Motion
Chapter 11 Siple Haronic Motion "We are to adit no ore causes of natural things than such as are both true and sufficient to explain their appearances." Isaac Newton 11.1 Introduction to Periodic Motion
More informationFEM-Design. Verification Examples. version Motto: ,,There is singularity between linear and nonlinear world. (Dr.
FEM-Design version.3 8 Motto:,,There is singularity between linear and nonlinear world. (Dr. Ire Bojtár) StruSoft AB Visit the StruSoft website for copany and FEM-Design inforation at www.strusoft.co Copyright
More informationEE5900 Spring Lecture 4 IC interconnect modeling methods Zhuo Feng
EE59 Spring Parallel LSI AD Algoriths Lecture I interconnect odeling ethods Zhuo Feng. Z. Feng MTU EE59 So far we ve considered only tie doain analyses We ll soon see that it is soeties preferable to odel
More informationAssessment of wind-induced structural fatigue based on joint probability density function of wind speed and direction
The 1 World Congress on Advances in Civil, Environental, and Materials Research (ACEM 1) eoul, Korea, August 6-3, 1 Assessent of wind-induced structural fatigue based on oint probability density function
More informationPolygonal Designs: Existence and Construction
Polygonal Designs: Existence and Construction John Hegean Departent of Matheatics, Stanford University, Stanford, CA 9405 Jeff Langford Departent of Matheatics, Drake University, Des Moines, IA 5011 G
More informationImproved Direct Displacement-Based Design Procedure for Performance-Based Seismic Design of Structures
Iproved Direct Displaceent-Based Design Procedure for Perforance-Based Seisic Design of Structures Rakesh K. Goel, California Polytechnic State University, San Luis Obispo, and Anil K. Chopra, University
More informationHee = ~ dxdy\jj+ (x) 'IJ+ (y) u (x- y) \jj (y) \jj (x), V, = ~ dx 'IJ+ (x) \jj (x) V (x), Hii = Z 2 ~ dx dy cp+ (x) cp+ (y) u (x- y) cp (y) cp (x),
SOVIET PHYSICS JETP VOLUME 14, NUMBER 4 APRIL, 1962 SHIFT OF ATOMIC ENERGY LEVELS IN A PLASMA L. E. PARGAMANIK Khar'kov State University Subitted to JETP editor February 16, 1961; resubitted June 19, 1961
More informationAbout the definition of parameters and regimes of active two-port networks with variable loads on the basis of projective geometry
About the definition of paraeters and regies of active two-port networks with variable loads on the basis of projective geoetry PENN ALEXANDR nstitute of Electronic Engineering and Nanotechnologies "D
More informationSupplementary Material for Fast and Provable Algorithms for Spectrally Sparse Signal Reconstruction via Low-Rank Hankel Matrix Completion
Suppleentary Material for Fast and Provable Algoriths for Spectrally Sparse Signal Reconstruction via Low-Ran Hanel Matrix Copletion Jian-Feng Cai Tianing Wang Ke Wei March 1, 017 Abstract We establish
More informationWaves & Normal Modes. Matt Jarvis
Waves & Noral Modes Matt Jarvis January 19, 016 Contents 1 Oscillations 1.0.1 Siple Haronic Motion - revision................... Noral Modes 5.1 The coupled pendulu.............................. 6.1.1
More informationDistributed Subgradient Methods for Multi-agent Optimization
1 Distributed Subgradient Methods for Multi-agent Optiization Angelia Nedić and Asuan Ozdaglar October 29, 2007 Abstract We study a distributed coputation odel for optiizing a su of convex objective functions
More informationEFFECTIVE MODAL MASS & MODAL PARTICIPATION FACTORS Revision I
EFFECTIVE MODA MASS & MODA PARTICIPATION FACTORS Revision I B To Irvine Eail: to@vibrationdata.co Deceber, 5 Introduction The effective odal ass provides a ethod for judging the significance of a vibration
More informationPhysics 221B: Solution to HW # 6. 1) Born-Oppenheimer for Coupled Harmonic Oscillators
Physics B: Solution to HW # 6 ) Born-Oppenheier for Coupled Haronic Oscillators This proble is eant to convince you of the validity of the Born-Oppenheier BO) Approxiation through a toy odel of coupled
More informationIII.H Zeroth Order Hydrodynamics
III.H Zeroth Order Hydrodynaics As a first approxiation, we shall assue that in local equilibriu, the density f 1 at each point in space can be represented as in eq.iii.56, i.e. f 0 1 p, q, t = n q, t
More informationConstrained Consensus and Optimization in Multi-Agent Networks arxiv: v2 [math.oc] 17 Dec 2008
LIDS Report 2779 1 Constrained Consensus and Optiization in Multi-Agent Networks arxiv:0802.3922v2 [ath.oc] 17 Dec 2008 Angelia Nedić, Asuan Ozdaglar, and Pablo A. Parrilo February 15, 2013 Abstract We
More informationProc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES
Proc. of the IEEE/OES Seventh Working Conference on Current Measureent Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Belinda Lipa Codar Ocean Sensors 15 La Sandra Way, Portola Valley, CA 98 blipa@pogo.co
More informationMechanics Physics 151
Mechanics Physics 5 Lecture Oscillations (Chapter 6) What We Did Last Tie Analyzed the otion of a heavy top Reduced into -diensional proble of θ Qualitative behavior Precession + nutation Initial condition
More informationExperimental Based Substructuring Using a Craig-Bampton Transmission Simulator Model
Experiental Based Substructuring Using a raig-bapton ransission Siulator Model Mathew S. llen, Daniel. Kaer Departent of Engineering Physics University of Wisconsin Madison, WI 376 kaer@engr.wisc.edu,
More informationmidas Civil Dynamic Analysis
Edgar De Los Santos Midas IT August 23 rd 2017 Contents: Introduction Eigen Value Analysis Response Spectrum Analysis Pushover Analysis Time History Analysis Seismic Analysis Seismic Analysis The seismic
More informationME Machine Design I. FINAL EXAM. OPEN BOOK AND CLOSED NOTES. Friday, May 8th, 2009
ME 5 - Machine Design I Spring Seester 009 Nae Lab. Div. FINAL EXAM. OPEN BOOK AND LOSED NOTES. Friday, May 8th, 009 Please use the blank paper for your solutions. Write on one side of the paper only.
More informationReducing Vibration and Providing Robustness with Multi-Input Shapers
29 Aerican Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June -2, 29 WeA6.4 Reducing Vibration and Providing Robustness with Multi-Input Shapers Joshua Vaughan and Willia Singhose Abstract
More informationGeneralized Rayleigh Wave Dispersion in a Covered Half-space Made of Viscoelastic Materials
Copyright 7 Tech Science Press CMC vol.53 no.4 pp.37-34 7 Generalized Rayleigh Wave Dispersion in a Covered Half-space Made of Viscoelastic Materials S.D. Akbarov and M. Negin 3 Abstract: Dispersion of
More informationA Note on Scheduling Tall/Small Multiprocessor Tasks with Unit Processing Time to Minimize Maximum Tardiness
A Note on Scheduling Tall/Sall Multiprocessor Tasks with Unit Processing Tie to Miniize Maxiu Tardiness Philippe Baptiste and Baruch Schieber IBM T.J. Watson Research Center P.O. Box 218, Yorktown Heights,
More informationOptical Properties of Plasmas of High-Z Elements
Forschungszentru Karlsruhe Techni und Uwelt Wissenschaftlishe Berichte FZK Optical Properties of Plasas of High-Z Eleents V.Tolach 1, G.Miloshevsy 1, H.Würz Project Kernfusion 1 Heat and Mass Transfer
More informationKernel Methods and Support Vector Machines
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley ENSIAG 2 / osig 1 Second Seester 2012/2013 Lesson 20 2 ay 2013 Kernel ethods and Support Vector achines Contents Kernel Functions...2 Quadratic
More informationSpine Fin Efficiency A Three Sided Pyramidal Fin of Equilateral Triangular Cross-Sectional Area
Proceedings of the 006 WSEAS/IASME International Conference on Heat and Mass Transfer, Miai, Florida, USA, January 18-0, 006 (pp13-18) Spine Fin Efficiency A Three Sided Pyraidal Fin of Equilateral Triangular
More informationPeriodic Motion is everywhere
Lecture 19 Goals: Chapter 14 Interrelate the physics and atheatics of oscillations. Draw and interpret oscillatory graphs. Learn the concepts of phase and phase constant. Understand and use energy conservation
More informationPhysics 201, Lecture 15
Physics 0, Lecture 5 Today s Topics q More on Linear Moentu And Collisions Elastic and Perfect Inelastic Collision (D) Two Diensional Elastic Collisions Exercise: Billiards Board Explosion q Multi-Particle
More informationVibration Characteristics of Cardboard Inserts in Shells
2003-01-1489 Vibration Characteristics of Cardboard Inserts in Shells Martin G. Foulkes and Jaes P. De Clerck General Motors Corporation Raendra Singh he Ohio State University Copyright 2003 SAE International
More informationSuccessive Model-Updating of the dynamic behaviour of casing bodies on a practical example of an axial piston pump
Successive Model-Updating of the dynaic behaviour of casing bodies on a practical exaple of an axial piston pup Ulrich Bittner Bosch Rexroth, Horb, Gerany Suary: These days, generally all new products
More information1 (40) Gravitational Systems Two heavy spherical (radius 0.05R) objects are located at fixed positions along
(40) Gravitational Systes Two heavy spherical (radius 0.05) objects are located at fixed positions along 2M 2M 0 an axis in space. The first ass is centered at r = 0 and has a ass of 2M. The second ass
More informationFlipping Physics Lecture Notes: Free Response Question #1 - AP Physics Exam Solutions
2015 FRQ #1 Free Response Question #1 - AP Physics 1-2015 Exa Solutions (a) First off, we know both blocks have a force of gravity acting downward on the. et s label the F & F. We also know there is a
More informationINNER CONSTRAINTS FOR A 3-D SURVEY NETWORK
eospatial Science INNER CONSRAINS FOR A 3-D SURVEY NEWORK hese notes follow closely the developent of inner constraint equations by Dr Willie an, Departent of Building, School of Design and Environent,
More informationA new type of lower bound for the largest eigenvalue of a symmetric matrix
Linear Algebra and its Applications 47 7 9 9 www.elsevier.co/locate/laa A new type of lower bound for the largest eigenvalue of a syetric atrix Piet Van Mieghe Delft University of Technology, P.O. Box
More informationAccuracy of the Scaling Law for Experimental Natural Frequencies of Rectangular Thin Plates
The 9th Conference of Mechanical Engineering Network of Thailand 9- October 005, Phuket, Thailand Accuracy of the caling Law for Experiental Natural Frequencies of Rectangular Thin Plates Anawat Na songkhla
More informationPage 1 Lab 1 Elementary Matrix and Linear Algebra Spring 2011
Page Lab Eleentary Matri and Linear Algebra Spring 0 Nae Due /03/0 Score /5 Probles through 4 are each worth 4 points.. Go to the Linear Algebra oolkit site ransforing a atri to reduced row echelon for
More informationNumerical Modeling of Self-Compacting Mortar Flow Using Discrete Element Method
Nuerical Modeling of Self-Copacting Flow Using Discrete Eleent Method - Technical Paper - Miansong HUANG *1, Xuehui AN *, Takayuki OBARA *3 and Masahiro OUCHI *4 ABSTRACT A nuerical odeling of Self-Copacting
More information