Elements of Coordinate System Transformations

Size: px
Start display at page:

Download "Elements of Coordinate System Transformations"

Transcription

1 B Elements of Coordinate System Transformations Coordinate system transformation is a powerfl tool for solving many geometrical and kinematic problems that pertain to the design of gear ctting tools and the kinematics of gear machining processes. Conseqent coordinate system transformations can easily be described analytically with the implementation of matrices. The se of matrices for coordinate system transformation can be traced back to the late 1940s [41], early 1950s [11, 42 and others]. The implementation of coordinate system transformations is necessary for representation in a common coordinate system a) of the gear ctting tool, and b) of its motion relative to the tooth flank of the work-gear. At every instant of time, the configration position and orientation) of the gear ctting tool relative to the work-gear can be described analytically with the help of a homogeneos transformation matrix corresponding to the displacement of the ctting tool from its crrent location to a certain consective location. B.1 Coordinate System Transformation In this text, coordinate system transformations are discssed briefly from the standpoint of their implementation for the prposes of gear ctting tool design. The interested reader may wish to go to [66, 70, 73, 74, 80] and other advanced sorces for details. B.1.1 Introdction Homogenos coordinates tilize a mathematical trick to embed three-dimensional coordinates and transformations into a for-dimensional matrix format. As a reslt, inversions or combinations of linear transformations are simplified to inversion or mltiplication of the corresponding matrices. Geometry of Srfaces: A Practical Gide for Mechanical Engineers, First Edition. Stephen P. Radzevich John Wiley & Sons, Ltd. Pblished 2013 by John Wiley & Sons, Ltd.

2 212 Appendix B: Elements of Coordinate System Transformations Homogenos coordinate vectors. Instead of representing each point rx, y, z) in threedimensional space with a single three-dimensional vector x r = y z B.1) homogenos coordinates allow each point rx, y, z) to be represented by any of an infinite nmber of for-dimensional vectors T x r = T y T z T B.2) The three-dimensional vector corresponding to any for-dimensional vector can be compted by dividing the first three elements by the forth, and a for-dimensional vector corresponding to any three-dimensional vector can be created by simply adding a forth element and setting it eqal to one. Homogenos coordinate transformation matrices of the dimension 4 4. Homogenos coordinate transformation matrices operate on for-dimensional homogenos vector representations of traditional three-dimensional coordinate locations. Any three-dimensional linear transformation translation, rotation, etc.) can be represented by a 4 4 homogenos coordinate transformation matrix. In fact, becase of the redndant representation of three-space in a homogenos coordinate system, an infinite nmber of different 4 4 homogenos coordinate transformation matrices is available to perform any given linear transformation. This redndancy can be eliminated to provide a niqe representation by dividing all elements of a 4 4 homogenos transformation matrix by the last element which will become eqal to one). This means that a 4 4 homogenos transformation matrix can incorporate as many as 15 independent parameters. The generic format representation of a homogenos transformation eqation for mapping the three-dimensional coordinate x 1, y 1, z 1 ) to the three-dimensional coordinate x 2, y 2, z 2 )is T x 2 T a T b T c T d T x 2 T y 2 T z 2 = T e T f T g T h T i T j T k T m T y 2 T z 2 T T n T p T q T T B.3) If any two matrices or vectors of this eqation are known, the third matrix or vector) can be compted and then the redndant T element in the soltion can be eliminated by dividing all elements of the matrix by the last element. Varios transformation models can be sed to constrain the form of the matrix to transformations with fewer degrees of freedom.

3 Appendix B: Elements of Coordinate System Transformations 213 a x a 2 a z X1 X2 a) b) c) Figre B.1 Analytical description of the operators of translation Tr a x, X), Tr a y, Y ), Tr a z, Z) along the coordinate axes. B.1.2 Translations Translation of a coordinate system is one of the major linear transformations sed for the prposes of gear ctting tool design. Translations of the coordinate system along the axes of the coordinate system are illstrated in Fig. B.1. The translations can be described analytically by the homogenos transformation matrices of dimension 4 4. For the analytical description of translation along coordinate axes, the operators of translation Tr a x, X), Tr a y, Y ) and Tr a z, Z) are sed. The operators yield matrix representation in the form a x Tr a x, X) = B.4) Tr a y, Y ) = a y B.5) Tr a z, Z) = a z B.6) Here, a x, a y, a z are signed vales that denote distances of translations along corresponding axes. Consider the two coordinate systems and shifted along the axis on a x [Fig. B.1a)]. Let a point m in the coordinate system be given by the position vector r 2 m). In the coordinate system that same point m can be specified by the position vector r 1 m). Then, the position vector r 1 m) can be expressed in terms of the position vector

4 214 Appendix B: Elements of Coordinate System Transformations ϕ x ϕ y ϕ z a) b) c) Figre B.2 Analytical description of the operators of rotation Rt ϕ x, X), Rt ϕ y, Y ), Rt ϕ z, Z) abot the coordinate axes. r 2 m) by the eqation r 1 m) = Tr a x, X) r 2 m). Eqations similar to that above are valid for other operators Tr a y, Y ) and Tr a z, Z) of the coordinate system transformation [Fig. B.1b, c)]. Any coordinate system transformation that doesn t change the orientation of a geometrical object is an orientation-preserving transformation, or a direct transformation. Therefore, the transformation of translation is an example of a direct transformation. B.1.3 Rotation abot a coordinate axis Rotation of a coordinate system abot a coordinate axis is another kind of major linear transformation sed for the prpose of gear ctting tool design. Rotation of the coordinate system abot the axis of the coordinate system is illstrated in Fig. B.2. For an analytical description of rotation abot coordinate axes, the operators of rotation Rt ϕ x, X), Rt ϕ y, Y ) and Rt ϕ z, Z) are sed. The operators yield representation in the form of homogenos matrices Rt ϕ x, X) = 0 cos ϕ x sin ϕ x 0 0 sin ϕ x cos ϕ x 0 cos ϕ y 0 sin ϕ y 0 Rt ϕ y, Y ) = sin ϕ y 0 cos ϕ y 0 cos ϕ z sin ϕ z 0 0 Rt ϕ z, Z) = sin ϕ z cos ϕ z B.7) B.8) B.9)

5 Appendix B: Elements of Coordinate System Transformations 215 Here, ϕ x, ϕ y, ϕ z are signed vales that denote angles of rotation abot the corresponding axis: ϕ x is rotation arond the X-axis pitch); ϕ y is rotation arond the Y -axis roll), and ϕ z is rotation arond the Z-axis yaw). Consider two coordinate systems and trned abot the -axis throgh the angle ϕ x [Fig. B.2a)]. In the coordinate system, a certain point m is given by the position vector r 2 m). In the coordinate system, that same point m can be specified by the position vector r 1 m). Then, the position vector r 1 m) can be expressed in terms of the position vector r 2 m) by the eqation r 1 m) = Rt ϕ x, X) r 2 m). Eqations similar to that above are valid for other operators Rt ϕ y, Y ) and Rt ϕ z, Z) of the coordinate system transformation [Fig. B.2b, c)]. B.1.4 Resltant coordinate system transformation The operators of translation Tr a x, X), Tr a y, Y ), Tr a z, Z) together with the operators of rotation Rt ϕ x, X), Rt ϕ y, Y ), Rt ϕ z, Z) are sed to compose the operator Rs 1 2) of the resltant coordinate system transformation. The operator Rs 1 2) of the resltant coordinate system transformation describes analytically the transition from the initial coordinate system to a certain coordinate system. Consider three conseqent translations along the coordinate axes, and. Sppose that a point P on a rigid body goes throgh a translation describing a straight path from P 1 to P 2 with a change of coordinates a x, a y, a z ). This motion can be described with an operator of the resltant coordinate system transformation Rs 1 2), which can be expressed in terms of the operators Tr a x, X), Tr a y, Y ), Tr a z, Z) of elementary coordinate system transformations. The operator Rs 1 2) is eqal to a x Rs 1 2) = Tr a z, Z) Tr a y, Y ) Tr a x, X) = a y a z B.10) In this particlar case, the operator of the resltant coordinate system transformation Rs 1 2) can be interpreted as the operator Tr a, A) of translation along an axis A [Rs 1 2) = Tr a, A)]. Evidently, the axis A is always an axis throgh the origin. Similarly, three conseqent rotations abot coordinate axes can be described with another operator of the resltant coordinate system transformation Rs 1 2) Rs 1 2) = Rt ϕ z, Z) Rt ϕ y, Y ) Rt ϕ x, X) B.11) In this particlar case, the operator of the resltant coordinate system transformation Rs 1 2) can be interpreted as the operator Rt ϕ,a) of rotation abot an axis A [Rs 1 2) = Rt ϕ,a)]. Evidently, the axis A is always an axis throgh the origin. In practice, it is often necessary to perform coordinate system transformations comprising translations along and rotations abot coordinate axes. For example, the expression Rs 1 5) = Tr a x, X) Rt ϕ z, Z) Rt ϕ x, X) Tr a y, Y ) B.12)

6 216 Appendix B: Elements of Coordinate System Transformations a x1 Z 4 Z 5 a y4 Z 3 X1, X2 Y 5 Y 4 Z 3 X 5 X 4 Z 4 Y 5 a x1 r 1 M ) M a y4 Z 5 r 5 M ) X 5 ϕ y3 Y 3 ϕ z2 Y 3 X 3 X 3 Y 4 X 4 Figre B.3 Example of the resltant coordinate system transformation. indicates that the transition from the coordinate system to the coordinate system X 5 Y 5 Z 5 Fig. B.3) is performed in the following for steps: a) translation Tr a y, Y ), followed by b) rotation Rt ϕ x, X), followed by c) second rotation Rt ϕ z, Z),followedbyd) translation Tr a x, X). Ultimately, the eqality r 1 M) = Rt 5 1) r 5 M) is observed. When the operator Rs 1 t) of a resltant coordinate system transformation is known, then the transition in the opposite direction can be performed by means of the operator Rs t 1) of the inverse coordinate system transformation. The eqality Rs t 1) = Rs 1 1 t) B.13) is valid for the operator Rs t 1) of the inverse coordinate system transformation. B.1.5 Screw motion abot a coordinate axis Operators for analytical description of screw motions abot an axis of the Cartesian coordinate system are a particlar case of operators of the resltant coordinate system transformation. By definition Fig. B.4), the operator Sc x ϕ x, p x ) of a screw motion abot the X-axis of the Cartesian coordinate system XYZ is eqal to Sc x ϕ x, p x ) = Rt ϕ x, X) Tr a x, X) B.14)

7 Appendix B: Elements of Coordinate System Transformations 217 ax = px ϕx X1 X2 ϕ x Figre B.4 Analytical description of the operator of screw motion Sc x ϕ x, p x ). After sbstitting the operator of translation Tr a x, X) [Eq. B.4)] and the operator of rotation Rt ϕ x, X) [Eq. B.7)], Eq. B.14) casts into the expression p x ϕ x Sc x ϕ x, p x ) = 0 cos ϕ x sin ϕ x 0 0 sin ϕ x cos ϕ x 0 B.15) for comptation of the operator of the screw motion Sc x ϕ x, p x ) abot the X-axis. The operators of screw motion Sc y ϕ y, p y ) and Sc z ϕ z, p z ) abot the Y - and Z-axis correspondingly are defined in a way similar to that above. For the operator of screw motion Sc x ϕ x, p x ): Sc y ϕ y, p y ) = Rt ϕ y, Y ) Tr a y, Y ) Sc z ϕ z, p z ) = Rt ϕ z, Z) Tr a z, Z) B.16) B.17) Using Eqs B.5) and B.6) together with Eqs B.8) and B.9), one can come p with the expressions cos ϕ y 0 sin ϕ y 0 Sc y ϕ y, p y ) = p y ϕ y sin ϕ y 0 cos ϕ y 0 cos ϕ z sin ϕ z 0 0 Sc z ϕ z, p z ) = sin ϕ z cos ϕ z p z ϕ z B.18) B.19) for comptation of the operators of screw motion Sc y ϕ y, p y ) and Sc z ϕ z, p z ) abot the Y - and Z-axis.

8 218 Appendix B: Elements of Coordinate System Transformations ax = Rw ϕy R w V ϕ y ω Figre B.5 Illstration of the transformation of rolling Rl x ϕ y, Y ) of a coordinate system. Screw motions abot a coordinate axis, as well as screw srfaces, are common in the design of gear ctting tools. This makes practical se of the operators of screw motion Sc x ϕ x, p x ), Sc y ϕ y, p y ) and Sc z ϕ z, p z ) when designing gear ctting tools. If necessary, an operator of the screw motion abot an arbitrary axis either throgh the origin of the coordinate system or not throgh the origin of the coordinate system can be derived in a similar way as the operators Sc x ϕ x, p x ), Sc y ϕ y, p y ) and Sc z ϕ z, p z ). B.1.6 Rolling motion of a coordinate system One more practical combination of a rotation and a translation is often sed when designing a gear ctting tool. Consider a Cartesian coordinate system Fig. B.5). The coordinate system is traveling in the direction of the -axis. The speed of translation is denoted V. The coordinate system is rotating abot the -axis simltaneosly with the translation. The speed of rotation is denoted ω. Assme that the ratio V / ω is constant. Under sch a scenario, the resltant motion of the reference system to its arbitrary position allows interpretation in the form of rolling with no sliding of a cylinder of radis R w over the plane. The plane is parallel to the coordinate -plane, and remote from it at a distance R w. For comptation of the radis of the rolling cylinder, the expression R w = V / ω can be sed. Owing to the rolling of the cylinder of radis R w over the plane being performed with no sliding, a certain correspondence between the translation and rotation of the coordinate system is established. When the coordinate system trns throgh a certain angle ϕ y, the translation of the origin of the coordinate system along the -axis is eqal to a x = ϕ r R w. The transition from the coordinate system to the coordinate system can be described analytically by the operator of the resltant coordinate system transformation Rs 1 2), where Rs 1 2) = Rt ϕ y, ) Tr a x, ) B.20) Here, Tr a x, ) designates the operator of translation along the -axis, and Rt ϕ y, ) is the operator of rotation abot the -axis. Operators of the reslting coordinate system transformation of this kind are referred to as operators of rolling motion over a plane.

9 Appendix B: Elements of Coordinate System Transformations 219 When the translation is performed along the -axis, and the rotation is performed abot the -axis, the operator of rolling is denoted Rl x ϕ y, Y ). In this particlar case, the eqality Rl x ϕ y, Y ) = Rs 1 2) [see Eq. B.20)] is valid. Based on this eqality, the operator of rolling over a plane Rl x ϕ y, Y ) can be compted from the eqation cos ϕ y 0 sin ϕ y a x cos ϕ y Rl x ϕ y, Y ) = sin ϕ y 0 cos ϕ y a x sin ϕ y B.21) While rotation remains abot the -axis, translation can be performed not along the -axis bt along the -axis instead. For rolling of this kind, the operator of rolling is eqal to cos ϕ y 0 sin ϕ y a z sin ϕ y Rl z ϕ y, Y ) = sin ϕ y 0 cos ϕ y a z cos ϕ y B.22) For cases when the rotation is performed abot the -axis, the corresponding operators of rolling are Rl y ϕ x, X) = 0 cos ϕ x sin ϕ x a y cos ϕ x 0 sin ϕ x cos ϕ x a y sin ϕ x B.23) for the case of rolling along the -axis and Rl z ϕ x, X) = 0 cos ϕ x sin ϕ x a z sin ϕ x 0 sin ϕ x cos ϕ x a z cos ϕ x for the case of rolling along the -axis. Similar expressions can be derived for the case of rotation abot the -axis cos ϕ z sin ϕ z 0 a x cos ϕ z Rl x ϕ z, Z) = sin ϕ z cos ϕ z 0 a x sin ϕ z cos ϕ z sin ϕ z 0 a y sin ϕ z Rl y ϕ z, Z) = sin ϕ z cos ϕ z 0 a y cos ϕ z B.24) B.25) B.26)

10 220 Appendix B: Elements of Coordinate System Transformations ϕ 1 * * ϕ 2 R 1 * R 2 * O 1 P O 2 ϕ 2 ω 1 ω 2 C Figre B.6 Derivation of the operator of rolling Rr ϕ 1, ) of two coordinate systems. Use of the operators of rolling [Eqs B.21) throgh B.26)] significantly simplifies analytical description of the coordinate system transformations. B.1.7 Rolling of two coordinate systems When designing a gear ctting tool, combinations of two rotations abot parallel axes are of particlar interest. As an example, consider two Cartesian coordinate systems and as shown in Fig. B.6. The coordinate systems and are rotated abot their axes and. The axes of rotation are parallel to each other ). Rotations ω 1 and ω 2 of the coordinate systems can be interpreted so that a circle of a certain radis R 1 associated with the coordinate system is rolling with no sliding over a circle of corresponding radis R 2 associated with the coordinate system. When the centerdistance C is known, the radii R 1 and R 2 of the circles can be expressed in terms of the center-distance C and the given rotations ω 1 and ω 2. For the comptations, the formlae 1 R 1 = C 1 + R 2 = C 1 + B.27) B.28) can be sed. Here, the ratio ω 1 /ω 2 is denoted. In the initial configration, the - and -axes align with each other. The - and - axes are parallel to each other. In Fig. B.6, the initial configration of the coordinate systems and is labeled X 1 Y 1 Z 1 and X 2 Y 2 Z 2. When the coordinate system trns throgh a certain angle ϕ 1, the coordinate system trns throgh the corresponding angle ϕ 2. When the angle ϕ 1 is known, the corresponding angle ϕ 2 is eqal to ϕ 2 = ϕ 1 /.

11 Appendix B: Elements of Coordinate System Transformations 221 The transition from the coordinate system to the coordinate system can be described analytically by the operator of the resltant coordinate system transformation Rs 1 2). In the case nder consideration, the operator Rs 1 2) can be expressed in terms of the operators of the elementary coordinate system transformations Rs 1 2) = Rt ϕ 1, ) Rt ϕ 1 /, ) Tr C, ) B.29) Other eqivalent combinations of the operators of elementary coordinate system transformations can reslt in the same operator Rs 1 2) of the resltant coordinate system transformation. The interested reader may wish to derive the eqivalent expressions for the operator Rs 1 2). The operator of the resltant coordinate system transformation [see Eq. B.29)] is referred to as an operator of rolling motion over a cylinder. When rotations are performed arond the - and -axes, the operator of rolling motion over a cylinder is denoted Rr ϕ 1, ). In this particlar case, the eqality Rr ϕ 1, ) = Rs 1 2) [see Eq. B.29)] is valid. Based on this eqality, the operator of rolling Rr ϕ 1, ) over a cylinder can be compted from the eqation Rr ϕ 1, ) = sin cos ϕ ) ϕ ) sin ϕ ) cos ϕ ) 0 C B.30) For the inverse transformation, the inverse operator of rolling of two coordinate systems Rr ϕ 2, ) can be sed. It is eqal to Rr ϕ 2, ) = Rr 1 ϕ 1, ). In terms of operators of the elementary coordinate system transformations, the operator Rr ϕ 2, ) can be expressed as Rr ϕ 2, ) = Rt ϕ 1 /, ) Rt ϕ 1, ) Tr C, ) B.31) Other eqivalent combinations of the operators of elementary coordinate system transformations can reslt in the same operator Rr ϕ 2, ) of the resltant coordinate system transformation. The interested reader may wish to derive the eqivalent expressions for the operator Rr ϕ 2, ). For the comptation of the operator of rolling of two coordinate systems Rr ϕ 2, ), the eqation can be sed. cos Rr ϕ 2, ) = ϕ sin ϕ ) ) sin ϕ ) cos ϕ ) 0 C B.32)

12 222 Appendix B: Elements of Coordinate System Transformations Similar to the expression [see Eq. B.30)] derived for the comptation of the operator of rolling Rr ϕ 1, ) arond the - and -axes, corresponding formlae can be derived for the comptation of the operators of rolling Rr ϕ 1, ) and Rr ϕ 1, ) abot parallel axes and, as well as abot parallel axes and. Use of the operators of rolling abot two axes Rr ϕ 1, ), Rr ϕ 1, ) and Rr ϕ 1, ) sbstantially simplifies analytical description of the coordinate system transformations. B.2 Conversion of the Coordinate System Orientation Application of the matrix method of coordinate system transformation presmes that both coordinate systems i and i ± 1) are of the same hand. This means that it is assmed from the very beginning that both of them are either right-hand oriented or left-hand oriented Cartesian reference systems. In the event the coordinate systems i and i ± 1) are of opposite hand, one of the coordinate systems needs to be converted into an oppositely oriented Cartesian coordinate system. For conversion of a left-hand oriented Cartesian coordinate system into a right-hand oriented coordinate system and vice versa, operators of reflection are sed. In order to change the direction of the X i -axis of the initial coordinate system i to the opposite direction in this case, in the new coordinate system i ± 1) the eqalities X i±1 = X i, Y i±1 Y i and Z i±1 Z i are observed), the operator of reflection Rf x Y i Z i ) can be applied. The operator of reflection yields representation in matrix form as Rf x Y i Z i ) = B.33) Similarly, implementation of the operators of reflection Rf y X i Z i ) and Rf z X i Y i ) reslts in the directions of the Y i - and Z i -axes being reversed. The operators of reflection Rf y X i Z i ) and Rf z X i Y i ) in this case can be expressed analytically in matrix form as Rf y X i Z i ) = Rf z X i Y i ) = B.34) B.35) A linear transformation that reverses the direction of the coordinate axis is an opposite transformation. The transformation of reflection is an example of an orientation-reversing transformation.

13 Appendix B: Elements of Coordinate System Transformations 223 B.3 Transformation of Srface Fndamental Forms Every coordinate system transformation reslts in a corresponding change to the eqation of the gear tooth srface G and/or the generating srface T of the gear ctting tool. Becase of this, it is necessary to recalclate the coefficients of the first 1.g and of the second 2.g fndamental of the srfaces G as many times as the coordinate system transformation is performed. This rotine and time-consming operation can be eliminated if the operators of coordinate system transformations are sed directly on the fndamental forms 1.g and 2.g. After being compted in an initial coordinate system, the fndamental magnitdes E g, F g, G g, L g, M g and N g of the forms 1.g and 2.g can be determined in any new coordinate system sing for this prpose the operators of translation, rotation and reslting coordinate system transformation. Sch transformations of the fndamental magnitdes 1.g and 2.g become possible de to implementation of the formlae below. Consider a gear tooth srface G that is given by the eqation r g = r g U g, V g ), where U g, V g ) G. For convenience, the first fndamental form 1.g of the gear tooth srface G is represented in matrix form [66, 74] as E g F g 0 0 du g [ ] [ 1.g = dug dv g 0 0 ] F g G g dv g 0 0 B.36) Similarly, the eqation of the second fndamental form 2.g of the srface G can be given by the eqation L g M g 0 0 du g [ ] [ 2.g = dug dv g 0 0 ] M g N g dv g 0 0 B.37) The coordinate system transformation with the operator of the resltant linear transformation Rs 1 2) transfers the eqation r g = r g U g, V g ) of the gear tooth srface G, initially given in, to the eqation r g = r g U g, V g ) of that same srface P in a new coordinate system. It is clear that the position vector of a point of the tooth flank point G in the first reference system differs from the position vector of that same point in the second reference system i.e. r g r g ). The operator of the reslting linear transformation Rs 1 2) of the srface P having first 1.g and second 2.g fndamental forms from the initial coordinate system to the new coordinate system reslts in new fndamental forms expressed in the form [66, 73, 74] [ 1.g ] = Rs T 1 2) [ 1.g ] Rs 1 2) B.38) [ 2.g ] = Rs T 1 2) [ 2.g ] Rs 1 2) B.39) Eqations B.38) and B.39) reveal that after the coordinate system transformation is complete, the first 1.g and second 2.g fndamental forms of the srface G in the coordinate system

14 224 Appendix B: Elements of Coordinate System Transformations are expressed in terms of the first 1.g and second 2.g fndamental forms initially represented in the coordinate system. In order to convert the fndamental forms 1.g and 2.g to the new coordinate system, either the corresponding fndamental form 1.g or 2.g needs to be pre-mltiplied by Rs 1 2) and after that post-mltiplied by Rs T 1 2). Implementation of Eqs B.38) and B.39) significantly simplifies the formlae of transformations.eqations [ 1.c ] = Rs T 1 2) [ 1.c ] Rs 1 2) B.40) [ 2.c ] = Rs T 1 2) [ 2.c ] Rs 1 2) B.41) similar to those above are valid for the generating srface T of the gear ctting tool. Implementation of the elements of screw calcls for the prpose of coordinate system transformation is a possible way to enhance the approach. In this case, screw operators of jst one kind can be applied for all kinds of coordinate system transformations.

Advanced topics in Finite Element Method 3D truss structures. Jerzy Podgórski

Advanced topics in Finite Element Method 3D truss structures. Jerzy Podgórski Advanced topics in Finite Element Method 3D trss strctres Jerzy Podgórski Introdction Althogh 3D trss strctres have been arond for a long time, they have been sed very rarely ntil now. They are difficlt

More information

FEA Solution Procedure

FEA Solution Procedure EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations

More information

Section 7.4: Integration of Rational Functions by Partial Fractions

Section 7.4: Integration of Rational Functions by Partial Fractions Section 7.4: Integration of Rational Fnctions by Partial Fractions This is abot as complicated as it gets. The Method of Partial Fractions Ecept for a few very special cases, crrently we have no way to

More information

Geometric Image Manipulation. Lecture #4 Wednesday, January 24, 2018

Geometric Image Manipulation. Lecture #4 Wednesday, January 24, 2018 Geometric Image Maniplation Lectre 4 Wednesda, Janar 4, 08 Programming Assignment Image Maniplation: Contet To start with the obvios, an image is a D arra of piels Piel locations represent points on the

More information

Mean Value Formulae for Laplace and Heat Equation

Mean Value Formulae for Laplace and Heat Equation Mean Vale Formlae for Laplace and Heat Eqation Abhinav Parihar December 7, 03 Abstract Here I discss a method to constrct the mean vale theorem for the heat eqation. To constrct sch a formla ab initio,

More information

The Linear Quadratic Regulator

The Linear Quadratic Regulator 10 The Linear Qadratic Reglator 10.1 Problem formlation This chapter concerns optimal control of dynamical systems. Most of this development concerns linear models with a particlarly simple notion of optimality.

More information

Simplified Identification Scheme for Structures on a Flexible Base

Simplified Identification Scheme for Structures on a Flexible Base Simplified Identification Scheme for Strctres on a Flexible Base L.M. Star California State University, Long Beach G. Mylonais University of Patras, Greece J.P. Stewart University of California, Los Angeles

More information

Technical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty

Technical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty Technical Note EN-FY160 Revision November 30, 016 ODiSI-B Sensor Strain Gage Factor Uncertainty Abstract Lna has pdated or strain sensor calibration tool to spport NIST-traceable measrements, to compte

More information

Nonlinear parametric optimization using cylindrical algebraic decomposition

Nonlinear parametric optimization using cylindrical algebraic decomposition Proceedings of the 44th IEEE Conference on Decision and Control, and the Eropean Control Conference 2005 Seville, Spain, December 12-15, 2005 TC08.5 Nonlinear parametric optimization sing cylindrical algebraic

More information

Reduction of over-determined systems of differential equations

Reduction of over-determined systems of differential equations Redction of oer-determined systems of differential eqations Maim Zaytse 1) 1, ) and Vyachesla Akkerman 1) Nclear Safety Institte, Rssian Academy of Sciences, Moscow, 115191 Rssia ) Department of Mechanical

More information

Classify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled.

Classify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled. Jnction elements in network models. Classify by nmber of ports and examine the possible strctres that reslt. Using only one-port elements, no more than two elements can be assembled. Combining two two-ports

More information

An Investigation into Estimating Type B Degrees of Freedom

An Investigation into Estimating Type B Degrees of Freedom An Investigation into Estimating Type B Degrees of H. Castrp President, Integrated Sciences Grop Jne, 00 Backgrond The degrees of freedom associated with an ncertainty estimate qantifies the amont of information

More information

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University 9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)

More information

The Cross Product of Two Vectors in Space DEFINITION. Cross Product. u * v = s ƒ u ƒƒv ƒ sin ud n

The Cross Product of Two Vectors in Space DEFINITION. Cross Product. u * v = s ƒ u ƒƒv ƒ sin ud n 12.4 The Cross Prodct 873 12.4 The Cross Prodct In stdying lines in the plane, when we needed to describe how a line was tilting, we sed the notions of slope and angle of inclination. In space, we want

More information

Lecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2

Lecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 BIJU PATNAIK UNIVERSITY OF TECHNOLOGY, ODISHA Lectre Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 Prepared by, Dr. Sbhend Kmar Rath, BPUT, Odisha. Tring Machine- Miscellany UNIT 2 TURING MACHINE

More information

4 Primitive Equations

4 Primitive Equations 4 Primitive Eqations 4.1 Spherical coordinates 4.1.1 Usefl identities We now introdce the special case of spherical coordinates: (,, r) (longitde, latitde, radial distance from Earth s center), with 0

More information

Chapter 2 Introduction to the Stiffness (Displacement) Method. The Stiffness (Displacement) Method

Chapter 2 Introduction to the Stiffness (Displacement) Method. The Stiffness (Displacement) Method CIVL 7/87 Chater - The Stiffness Method / Chater Introdction to the Stiffness (Dislacement) Method Learning Objectives To define the stiffness matrix To derive the stiffness matrix for a sring element

More information

FEA Solution Procedure

FEA Solution Procedure EA Soltion Procedre (demonstrated with a -D bar element problem) MAE 5 - inite Element Analysis Several slides from this set are adapted from B.S. Altan, Michigan Technological University EA Procedre for

More information

Solving a System of Equations

Solving a System of Equations Solving a System of Eqations Objectives Understand how to solve a system of eqations with: - Gass Elimination Method - LU Decomposition Method - Gass-Seidel Method - Jacobi Method A system of linear algebraic

More information

Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introduction The transmission line equations are given by

Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introduction The transmission line equations are given by Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introdction The transmission line eqations are given by, I z, t V z t l z t I z, t V z, t c z t (1) (2) Where, c is the per-nit-length

More information

Chem 4501 Introduction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics. Fall Semester Homework Problem Set Number 10 Solutions

Chem 4501 Introduction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics. Fall Semester Homework Problem Set Number 10 Solutions Chem 4501 Introdction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics Fall Semester 2017 Homework Problem Set Nmber 10 Soltions 1. McQarrie and Simon, 10-4. Paraphrase: Apply Eler s theorem

More information

FEA Solution Procedure

FEA Solution Procedure EA Soltion rocedre (demonstrated with a -D bar element problem) MAE - inite Element Analysis Many slides from this set are originally from B.S. Altan, Michigan Technological U. EA rocedre for Static Analysis.

More information

Control Systems

Control Systems 6.5 Control Systems Last Time: Introdction Motivation Corse Overview Project Math. Descriptions of Systems ~ Review Classification of Systems Linear Systems LTI Systems The notion of state and state variables

More information

EXERCISES WAVE EQUATION. In Problems 1 and 2 solve the heat equation (1) subject to the given conditions. Assume a rod of length L.

EXERCISES WAVE EQUATION. In Problems 1 and 2 solve the heat equation (1) subject to the given conditions. Assume a rod of length L. .4 WAVE EQUATION 445 EXERCISES.3 In Problems and solve the heat eqation () sbject to the given conditions. Assme a rod of length.. (, t), (, t) (, ),, > >. (, t), (, t) (, ) ( ) 3. Find the temperatre

More information

called the potential flow, and function φ is called the velocity potential.

called the potential flow, and function φ is called the velocity potential. J. Szantr Lectre No. 3 Potential flows 1 If the flid flow is irrotational, i.e. everwhere or almost everwhere in the field of flow there is rot 0 it means that there eists a scalar fnction ϕ,, z), sch

More information

Restricted Three-Body Problem in Different Coordinate Systems

Restricted Three-Body Problem in Different Coordinate Systems Applied Mathematics 3 949-953 http://dx.doi.org/.436/am..394 Pblished Online September (http://www.scirp.org/jornal/am) Restricted Three-Body Problem in Different Coordinate Systems II-In Sidereal Spherical

More information

The Dual of the Maximum Likelihood Method

The Dual of the Maximum Likelihood Method Department of Agricltral and Resorce Economics University of California, Davis The Dal of the Maximm Likelihood Method by Qirino Paris Working Paper No. 12-002 2012 Copyright @ 2012 by Qirino Paris All

More information

4 Exact laminar boundary layer solutions

4 Exact laminar boundary layer solutions 4 Eact laminar bondary layer soltions 4.1 Bondary layer on a flat plate (Blasis 1908 In Sec. 3, we derived the bondary layer eqations for 2D incompressible flow of constant viscosity past a weakly crved

More information

Setting The K Value And Polarization Mode Of The Delta Undulator

Setting The K Value And Polarization Mode Of The Delta Undulator LCLS-TN-4- Setting The Vale And Polarization Mode Of The Delta Undlator Zachary Wolf, Heinz-Dieter Nhn SLAC September 4, 04 Abstract This note provides the details for setting the longitdinal positions

More information

1. Tractable and Intractable Computational Problems So far in the course we have seen many problems that have polynomial-time solutions; that is, on

1. Tractable and Intractable Computational Problems So far in the course we have seen many problems that have polynomial-time solutions; that is, on . Tractable and Intractable Comptational Problems So far in the corse we have seen many problems that have polynomial-time soltions; that is, on a problem instance of size n, the rnning time T (n) = O(n

More information

A Decomposition Method for Volume Flux. and Average Velocity of Thin Film Flow. of a Third Grade Fluid Down an Inclined Plane

A Decomposition Method for Volume Flux. and Average Velocity of Thin Film Flow. of a Third Grade Fluid Down an Inclined Plane Adv. Theor. Appl. Mech., Vol. 1, 8, no. 1, 9 A Decomposition Method for Volme Flx and Average Velocit of Thin Film Flow of a Third Grade Flid Down an Inclined Plane A. Sadighi, D.D. Ganji,. Sabzehmeidani

More information

1 Differential Equations for Solid Mechanics

1 Differential Equations for Solid Mechanics 1 Differential Eqations for Solid Mechanics Simple problems involving homogeneos stress states have been considered so far, wherein the stress is the same throghot the component nder std. An eception to

More information

Solutions to Math 152 Review Problems for Exam 1

Solutions to Math 152 Review Problems for Exam 1 Soltions to Math 5 Review Problems for Eam () If A() is the area of the rectangle formed when the solid is sliced at perpendiclar to the -ais, then A() = ( ), becase the height of the rectangle is and

More information

Math 263 Assignment #3 Solutions. 1. A function z = f(x, y) is called harmonic if it satisfies Laplace s equation:

Math 263 Assignment #3 Solutions. 1. A function z = f(x, y) is called harmonic if it satisfies Laplace s equation: Math 263 Assignment #3 Soltions 1. A fnction z f(x, ) is called harmonic if it satisfies Laplace s eqation: 2 + 2 z 2 0 Determine whether or not the following are harmonic. (a) z x 2 + 2. We se the one-variable

More information

APPENDIX B MATRIX NOTATION. The Definition of Matrix Notation is the Definition of Matrix Multiplication B.1 INTRODUCTION

APPENDIX B MATRIX NOTATION. The Definition of Matrix Notation is the Definition of Matrix Multiplication B.1 INTRODUCTION APPENDIX B MAIX NOAION he Deinition o Matrix Notation is the Deinition o Matrix Mltiplication B. INODUCION { XE "Matrix Mltiplication" }{ XE "Matrix Notation" }he se o matrix notations is not necessary

More information

Motion in One Dimension. A body is moving with velocity 3ms towards East. After s its velocity becomes 4ms towards North. The average acceleration of the body is a) 7ms b) 7ms c) 5ms d) ms. A boy standing

More information

WEAR PREDICTION OF A TOTAL KNEE PROSTHESIS TIBIAL TRAY

WEAR PREDICTION OF A TOTAL KNEE PROSTHESIS TIBIAL TRAY APPLIED PHYSICS MEDICAL WEAR PREDICTION OF A TOTAL KNEE PROSTHESIS TIBIAL TRAY L. CÃPITANU, A. IAROVICI, J. ONIªORU Institte of Solid Mechanics, Romanian Academy, Constantin Mille 5, Bcharest Received

More information

III. Demonstration of a seismometer response with amplitude and phase responses at:

III. Demonstration of a seismometer response with amplitude and phase responses at: GG5330, Spring semester 006 Assignment #1, Seismometry and Grond Motions De 30 Janary 006. 1. Calibration Of A Seismometer Using Java: A really nifty se of Java is now available for demonstrating the seismic

More information

PROBLEMS

PROBLEMS PROBLEMS------------------------------------------------ - 7- Thermodynamic Variables and the Eqation of State 1. Compter (a) the nmber of moles and (b) the nmber of molecles in 1.00 cm of an ideal gas

More information

3.3 Operations With Vectors, Linear Combinations

3.3 Operations With Vectors, Linear Combinations Operations With Vectors, Linear Combinations Performance Criteria: (d) Mltiply ectors by scalars and add ectors, algebraically Find linear combinations of ectors algebraically (e) Illstrate the parallelogram

More information

L = 2 λ 2 = λ (1) In other words, the wavelength of the wave in question equals to the string length,

L = 2 λ 2 = λ (1) In other words, the wavelength of the wave in question equals to the string length, PHY 309 L. Soltions for Problem set # 6. Textbook problem Q.20 at the end of chapter 5: For any standing wave on a string, the distance between neighboring nodes is λ/2, one half of the wavelength. The

More information

CONTENTS. INTRODUCTION MEQ curriculum objectives for vectors (8% of year). page 2 What is a vector? What is a scalar? page 3, 4

CONTENTS. INTRODUCTION MEQ curriculum objectives for vectors (8% of year). page 2 What is a vector? What is a scalar? page 3, 4 CONTENTS INTRODUCTION MEQ crriclm objectives for vectors (8% of year). page 2 What is a vector? What is a scalar? page 3, 4 VECTOR CONCEPTS FROM GEOMETRIC AND ALGEBRAIC PERSPECTIVES page 1 Representation

More information

Second-Order Wave Equation

Second-Order Wave Equation Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order

More information

Computational Fluid Dynamics Simulation and Wind Tunnel Testing on Microlight Model

Computational Fluid Dynamics Simulation and Wind Tunnel Testing on Microlight Model Comptational Flid Dynamics Simlation and Wind Tnnel Testing on Microlight Model Iskandar Shah Bin Ishak Department of Aeronatics and Atomotive, Universiti Teknologi Malaysia T.M. Kit Universiti Teknologi

More information

Discontinuous Fluctuation Distribution for Time-Dependent Problems

Discontinuous Fluctuation Distribution for Time-Dependent Problems Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation

More information

Evaluation of the Fiberglass-Reinforced Plastics Interfacial Behavior by using Ultrasonic Wave Propagation Method

Evaluation of the Fiberglass-Reinforced Plastics Interfacial Behavior by using Ultrasonic Wave Propagation Method 17th World Conference on Nondestrctive Testing, 5-8 Oct 008, Shanghai, China Evalation of the Fiberglass-Reinforced Plastics Interfacial Behavior by sing Ultrasonic Wave Propagation Method Jnjie CHANG

More information

Dynamic Optimization of First-Order Systems via Static Parametric Programming: Application to Electrical Discharge Machining

Dynamic Optimization of First-Order Systems via Static Parametric Programming: Application to Electrical Discharge Machining Dynamic Optimization of First-Order Systems via Static Parametric Programming: Application to Electrical Discharge Machining P. Hgenin*, B. Srinivasan*, F. Altpeter**, R. Longchamp* * Laboratoire d Atomatiqe,

More information

Regression Analysis of Octal Rings as Mechanical Force Transducers

Regression Analysis of Octal Rings as Mechanical Force Transducers Regression Analysis of Octal Rings as Mechanical Force Transdcers KHALED A. ABUHASEL* & ESSAM SOLIMAN** *Department of Mechanical Engineering College of Engineering, University of Bisha, Bisha, Kingdom

More information

1. State-Space Linear Systems 2. Block Diagrams 3. Exercises

1. State-Space Linear Systems 2. Block Diagrams 3. Exercises LECTURE 1 State-Space Linear Sstems This lectre introdces state-space linear sstems, which are the main focs of this book. Contents 1. State-Space Linear Sstems 2. Block Diagrams 3. Exercises 1.1 State-Space

More information

STEP Support Programme. STEP III Hyperbolic Functions: Solutions

STEP Support Programme. STEP III Hyperbolic Functions: Solutions STEP Spport Programme STEP III Hyperbolic Fnctions: Soltions Start by sing the sbstittion t cosh x. This gives: sinh x cosh a cosh x cosh a sinh x t sinh x dt t dt t + ln t ln t + ln cosh a ln ln cosh

More information

Efficiency Increase and Input Power Decrease of Converted Prototype Pump Performance

Efficiency Increase and Input Power Decrease of Converted Prototype Pump Performance International Jornal of Flid Machinery and Systems DOI: http://dx.doi.org/10.593/ijfms.016.9.3.05 Vol. 9, No. 3, Jly-September 016 ISSN (Online): 188-9554 Original Paper Efficiency Increase and Inpt Power

More information

LINEAR COMBINATIONS AND SUBSPACES

LINEAR COMBINATIONS AND SUBSPACES CS131 Part II, Linear Algebra and Matrices CS131 Mathematics for Compter Scientists II Note 5 LINEAR COMBINATIONS AND SUBSPACES Linear combinations. In R 2 the vector (5, 3) can be written in the form

More information

Lecture 3. (2) Last time: 3D space. The dot product. Dan Nichols January 30, 2018

Lecture 3. (2) Last time: 3D space. The dot product. Dan Nichols January 30, 2018 Lectre 3 The dot prodct Dan Nichols nichols@math.mass.ed MATH 33, Spring 018 Uniersity of Massachsetts Janary 30, 018 () Last time: 3D space Right-hand rle, the three coordinate planes 3D coordinate system:

More information

Einstein Classes, Unit No. 102, 103, Vardhman Ring Road Plaza, Vikas Puri Extn., Outer Ring Road New Delhi , Ph. : ,

Einstein Classes, Unit No. 102, 103, Vardhman Ring Road Plaza, Vikas Puri Extn., Outer Ring Road New Delhi , Ph. : , PK K I N E M A T I C S Syllabs : Frame of reference. Motion in a straight line : Position-time graph, speed and velocity. Uniform and non-niform motion, average speed and instantaneos velocity. Uniformly

More information

arxiv: v1 [physics.flu-dyn] 11 Mar 2011

arxiv: v1 [physics.flu-dyn] 11 Mar 2011 arxiv:1103.45v1 [physics.fl-dyn 11 Mar 011 Interaction of a magnetic dipole with a slowly moving electrically condcting plate Evgeny V. Votyakov Comptational Science Laboratory UCY-CompSci, Department

More information

PhysicsAndMathsTutor.com

PhysicsAndMathsTutor.com . Two smooth niform spheres S and T have eqal radii. The mass of S is 0. kg and the mass of T is 0.6 kg. The spheres are moving on a smooth horizontal plane and collide obliqely. Immediately before the

More information

OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIELD OF A POLYHEDRAL BODY WITH LINEARLY INCREASING DENSITY 1

OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIELD OF A POLYHEDRAL BODY WITH LINEARLY INCREASING DENSITY 1 OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIEL OF A POLYHERAL BOY WITH LINEARLY INCREASING ENSITY 1 V. POHÁNKA2 Abstract The formla for the comptation of the gravity field of a polyhedral body

More information

Conditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane

Conditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane Filomat 3:2 (27), 376 377 https://doi.org/.2298/fil7276a Pblished by Faclty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Conditions for Approaching

More information

Lab Manual for Engrd 202, Virtual Torsion Experiment. Aluminum module

Lab Manual for Engrd 202, Virtual Torsion Experiment. Aluminum module Lab Manal for Engrd 202, Virtal Torsion Experiment Alminm modle Introdction In this modle, o will perform data redction and analsis for circlar cross section alminm samples. B plotting the torqe vs. twist

More information

Homogeneous Liner Systems with Constant Coefficients

Homogeneous Liner Systems with Constant Coefficients Homogeneos Liner Systems with Constant Coefficients Jly, 06 The object of stdy in this section is where A is a d d constant matrix whose entries are real nmbers. As before, we will look to the exponential

More information

1 The space of linear transformations from R n to R m :

1 The space of linear transformations from R n to R m : Math 540 Spring 20 Notes #4 Higher deriaties, Taylor s theorem The space of linear transformations from R n to R m We hae discssed linear transformations mapping R n to R m We can add sch linear transformations

More information

Modelling by Differential Equations from Properties of Phenomenon to its Investigation

Modelling by Differential Equations from Properties of Phenomenon to its Investigation Modelling by Differential Eqations from Properties of Phenomenon to its Investigation V. Kleiza and O. Prvinis Kanas University of Technology, Lithania Abstract The Panevezys camps of Kanas University

More information

Momentum Equation. Necessary because body is not made up of a fixed assembly of particles Its volume is the same however Imaginary

Momentum Equation. Necessary because body is not made up of a fixed assembly of particles Its volume is the same however Imaginary Momentm Eqation Interest in the momentm eqation: Qantification of proplsion rates esign strctres for power generation esign of pipeline systems to withstand forces at bends and other places where the flow

More information

MAT389 Fall 2016, Problem Set 6

MAT389 Fall 2016, Problem Set 6 MAT389 Fall 016, Problem Set 6 Trigonometric and hperbolic fnctions 6.1 Show that e iz = cos z + i sin z for eer comple nmber z. Hint: start from the right-hand side and work or wa towards the left-hand

More information

Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers

Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers D.R. Espinoza-Trejo and D.U. Campos-Delgado Facltad de Ingeniería, CIEP, UASLP, espinoza trejo dr@aslp.mx Facltad de Ciencias,

More information

u P(t) = P(x,y) r v t=0 4/4/2006 Motion ( F.Robilliard) 1

u P(t) = P(x,y) r v t=0 4/4/2006 Motion ( F.Robilliard) 1 y g j P(t) P(,y) r t0 i 4/4/006 Motion ( F.Robilliard) 1 Motion: We stdy in detail three cases of motion: 1. Motion in one dimension with constant acceleration niform linear motion.. Motion in two dimensions

More information

A New Approach to Direct Sequential Simulation that Accounts for the Proportional Effect: Direct Lognormal Simulation

A New Approach to Direct Sequential Simulation that Accounts for the Proportional Effect: Direct Lognormal Simulation A ew Approach to Direct eqential imlation that Acconts for the Proportional ffect: Direct ognormal imlation John Manchk, Oy eangthong and Clayton Detsch Department of Civil & nvironmental ngineering University

More information

Graphs and Their. Applications (6) K.M. Koh* F.M. Dong and E.G. Tay. 17 The Number of Spanning Trees

Graphs and Their. Applications (6) K.M. Koh* F.M. Dong and E.G. Tay. 17 The Number of Spanning Trees Graphs and Their Applications (6) by K.M. Koh* Department of Mathematics National University of Singapore, Singapore 1 ~ 7543 F.M. Dong and E.G. Tay Mathematics and Mathematics EdOOation National Institte

More information

A Survey of the Implementation of Numerical Schemes for Linear Advection Equation

A Survey of the Implementation of Numerical Schemes for Linear Advection Equation Advances in Pre Mathematics, 4, 4, 467-479 Pblished Online Agst 4 in SciRes. http://www.scirp.org/jornal/apm http://dx.doi.org/.436/apm.4.485 A Srvey of the Implementation of Nmerical Schemes for Linear

More information

Study on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom

Study on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom EPJ Web of Conferences 80, 0034 (08) EFM 07 Stdy on the implsive pressre of tank oscillating by force towards mltiple degrees of freedom Shigeyki Hibi,* The ational Defense Academy, Department of Mechanical

More information

4.2 First-Order Logic

4.2 First-Order Logic 64 First-Order Logic and Type Theory The problem can be seen in the two qestionable rles In the existential introdction, the term a has not yet been introdced into the derivation and its se can therefore

More information

Replacement of Einstein s Relativity Theory with a New One: Why the Second Postulate is Superfluous?

Replacement of Einstein s Relativity Theory with a New One: Why the Second Postulate is Superfluous? International Jornal of Physics, 16, Vol. 4, No. 5, 14-145 Available online at http://pbs.sciepb.com/ijp/4/5/5 Science and Edcation Pblishing DOI:1.1691/ijp-4-5-5 Replacement of Einstein s Relativity Theory

More information

CONSIDER an array of N sensors (residing in threedimensional. Investigating Hyperhelical Array Manifold Curves Using the Complex Cartan Matrix

CONSIDER an array of N sensors (residing in threedimensional. Investigating Hyperhelical Array Manifold Curves Using the Complex Cartan Matrix IEEE JOURNAL OF SELECTED TOPICS IN SIGNAL PROCESSING (SPECIAL ISSUE), VOL.?, NO.?,???? 1 Investigating yperhelical Array Manifold Crves Using the Complex Cartan Matrix Athanassios Manikas, Senior Member,

More information

CHANNEL SELECTION WITH RAYLEIGH FADING: A MULTI-ARMED BANDIT FRAMEWORK. Wassim Jouini and Christophe Moy

CHANNEL SELECTION WITH RAYLEIGH FADING: A MULTI-ARMED BANDIT FRAMEWORK. Wassim Jouini and Christophe Moy CHANNEL SELECTION WITH RAYLEIGH FADING: A MULTI-ARMED BANDIT FRAMEWORK Wassim Joini and Christophe Moy SUPELEC, IETR, SCEE, Avene de la Bolaie, CS 47601, 5576 Cesson Sévigné, France. INSERM U96 - IFR140-

More information

IMPROVED ANALYSIS OF BOLTED SHEAR CONNECTION UNDER ECCENTRIC LOADS

IMPROVED ANALYSIS OF BOLTED SHEAR CONNECTION UNDER ECCENTRIC LOADS Jornal of Marine Science and Technology, Vol. 5, No. 4, pp. 373-38 (17) 373 DOI: 1.6119/JMST-17-3-1 IMPROVED ANALYSIS OF BOLTED SHEAR ONNETION UNDER EENTRI LOADS Dng-Mya Le 1, heng-yen Liao, hien-hien

More information

A State Space Based Implicit Integration Algorithm. for Differential Algebraic Equations of Multibody. Dynamics

A State Space Based Implicit Integration Algorithm. for Differential Algebraic Equations of Multibody. Dynamics A State Space Based Implicit Integration Algorithm for Differential Algebraic Eqations of Mltibody Dynamics E. J. Hag, D. Negrt, M. Ianc Janary 28, 1997 To Appear Mechanics of Strctres and Machines Abstract.

More information

Constructive Root Bound for k-ary Rational Input Numbers

Constructive Root Bound for k-ary Rational Input Numbers Constrctive Root Bond for k-ary Rational Inpt Nmbers Sylvain Pion, Chee Yap To cite this version: Sylvain Pion, Chee Yap. Constrctive Root Bond for k-ary Rational Inpt Nmbers. 19th Annal ACM Symposim on

More information

3 2D Elastostatic Problems in Cartesian Coordinates

3 2D Elastostatic Problems in Cartesian Coordinates D lastostatic Problems in Cartesian Coordinates Two dimensional elastostatic problems are discssed in this Chapter, that is, static problems of either plane stress or plane strain. Cartesian coordinates

More information

A Fully-Neoclassical Finite-Orbit-Width Version. of the CQL3D Fokker-Planck code

A Fully-Neoclassical Finite-Orbit-Width Version. of the CQL3D Fokker-Planck code A Flly-Neoclassical Finite-Orbit-Width Version of the CQL3 Fokker-Planck code CompX eport: CompX-6- Jly, 6 Y. V. Petrov and. W. Harvey CompX, el Mar, CA 94, USA A Flly-Neoclassical Finite-Orbit-Width Version

More information

A Contraction of the Lucas Polygon

A Contraction of the Lucas Polygon Western Washington University Western CEDAR Mathematics College of Science and Engineering 4 A Contraction of the Lcas Polygon Branko Ćrgs Western Washington University, brankocrgs@wwed Follow this and

More information

Pulses on a Struck String

Pulses on a Struck String 8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a

More information

Finite Difference Method of Modelling Groundwater Flow

Finite Difference Method of Modelling Groundwater Flow Jornal of Water Resorce and Protection, 20, 3, 92-98 doi:0.4236/warp.20.33025 Pblished Online March 20 (http://www.scirp.org/ornal/warp) Finite Difference Method of Modelling Grondwater Flow Abstract Magns.

More information

Linear System Theory (Fall 2011): Homework 1. Solutions

Linear System Theory (Fall 2011): Homework 1. Solutions Linear System Theory (Fall 20): Homework Soltions De Sep. 29, 20 Exercise (C.T. Chen: Ex.3-8). Consider a linear system with inpt and otpt y. Three experiments are performed on this system sing the inpts

More information

Incompressible Viscoelastic Flow of a Generalised Oldroyed-B Fluid through Porous Medium between Two Infinite Parallel Plates in a Rotating System

Incompressible Viscoelastic Flow of a Generalised Oldroyed-B Fluid through Porous Medium between Two Infinite Parallel Plates in a Rotating System International Jornal of Compter Applications (97 8887) Volme 79 No., October Incompressible Viscoelastic Flow of a Generalised Oldroed-B Flid throgh Poros Medim between Two Infinite Parallel Plates in

More information

EXCITATION RATE COEFFICIENTS OF MOLYBDENUM ATOM AND IONS IN ASTROPHYSICAL PLASMA AS A FUNCTION OF ELECTRON TEMPERATURE

EXCITATION RATE COEFFICIENTS OF MOLYBDENUM ATOM AND IONS IN ASTROPHYSICAL PLASMA AS A FUNCTION OF ELECTRON TEMPERATURE EXCITATION RATE COEFFICIENTS OF MOLYBDENUM ATOM AND IONS IN ASTROPHYSICAL PLASMA AS A FUNCTION OF ELECTRON TEMPERATURE A.N. Jadhav Department of Electronics, Yeshwant Mahavidyalaya, Ned. Affiliated to

More information

RIGID BODY MOTION (Section 16.1)

RIGID BODY MOTION (Section 16.1) RIGID BODY MOTION (Section 16.1) There are cases where an object cannot be treated as a particle. In these cases the size or shape of the body must be considered. Rotation of the body about its center

More information

Konyalioglu, Serpil. Konyalioglu, A.Cihan. Ipek, A.Sabri. Isik, Ahmet

Konyalioglu, Serpil. Konyalioglu, A.Cihan. Ipek, A.Sabri. Isik, Ahmet The Role of Visalization Approach on Stdent s Conceptal Learning Konyaliogl, Serpil Department of Secondary Science and Mathematics Edcation, K.K. Edcation Faclty, Atatürk University, 25240- Erzrm-Trkey;

More information

PLANETARY ORBITS. According to MATTER (Re-examined) Nainan K. Varghese,

PLANETARY ORBITS. According to MATTER (Re-examined) Nainan K. Varghese, PLANETARY ORBITS According to MATTER (Re-examined) Nainan K. Varghese, matterdoc@gmail.com http://www.matterdoc.info Abstract: It is an established fact that sn is a moving macro body in space. By simple

More information

Multivariable Ripple-Free Deadbeat Control

Multivariable Ripple-Free Deadbeat Control Scholarly Jornal of Mathematics and Compter Science Vol. (), pp. 9-9, October Available online at http:// www.scholarly-jornals.com/sjmcs ISS 76-8947 Scholarly-Jornals Fll Length Research aper Mltivariable

More information

STABILIZATIO ON OF LONGITUDINAL AIRCRAFT MOTION USING MODEL PREDICTIVE CONTROL AND EXACT LINEARIZATION

STABILIZATIO ON OF LONGITUDINAL AIRCRAFT MOTION USING MODEL PREDICTIVE CONTROL AND EXACT LINEARIZATION 8 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES STABILIZATIO ON OF LONGITUDINAL AIRCRAFT MOTION USING MODEL PREDICTIVE CONTROL AND EXACT LINEARIZATION Čeliovsý S.*, Hospodář P.** *CTU Prage, Faclty

More information

3.4-Miscellaneous Equations

3.4-Miscellaneous Equations .-Miscellaneos Eqations Factoring Higher Degree Polynomials: Many higher degree polynomials can be solved by factoring. Of particlar vale is the method of factoring by groping, however all types of factoring

More information

Sensitivity Analysis in Bayesian Networks: From Single to Multiple Parameters

Sensitivity Analysis in Bayesian Networks: From Single to Multiple Parameters Sensitivity Analysis in Bayesian Networks: From Single to Mltiple Parameters Hei Chan and Adnan Darwiche Compter Science Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwiche}@cs.cla.ed

More information

A New Method for Calculating of Electric Fields Around or Inside Any Arbitrary Shape Electrode Configuration

A New Method for Calculating of Electric Fields Around or Inside Any Arbitrary Shape Electrode Configuration Proceedings of the 5th WSEAS Int. Conf. on Power Systems and Electromagnetic Compatibility, Corf, Greece, Agst 3-5, 005 (pp43-48) A New Method for Calclating of Electric Fields Arond or Inside Any Arbitrary

More information

Introdction Finite elds play an increasingly important role in modern digital commnication systems. Typical areas of applications are cryptographic sc

Introdction Finite elds play an increasingly important role in modern digital commnication systems. Typical areas of applications are cryptographic sc A New Architectre for a Parallel Finite Field Mltiplier with Low Complexity Based on Composite Fields Christof Paar y IEEE Transactions on Compters, Jly 996, vol 45, no 7, pp 856-86 Abstract In this paper

More information

UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL

UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL 8th International DAAAM Baltic Conference "INDUSTRIAL ENGINEERING - 19-1 April 01, Tallinn, Estonia UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL Põdra, P. & Laaneots, R. Abstract: Strength analysis is a

More information

Sources of Non Stationarity in the Semivariogram

Sources of Non Stationarity in the Semivariogram Sorces of Non Stationarity in the Semivariogram Migel A. Cba and Oy Leangthong Traditional ncertainty characterization techniqes sch as Simple Kriging or Seqential Gassian Simlation rely on stationary

More information

Vectors. February 1, 2010

Vectors. February 1, 2010 Vectors Febrary 1, 2010 Motivation Location of projector from crrent position and orientation: direction of projector distance to projector (direction, distance=magnitde) vector Examples: Force Velocity

More information

Lesson 81: The Cross Product of Vectors

Lesson 81: The Cross Product of Vectors Lesson 8: The Cross Prodct of Vectors IBHL - SANTOWSKI In this lesson yo will learn how to find the cross prodct of two ectors how to find an orthogonal ector to a plane defined by two ectors how to find

More information

Approach to a Proof of the Riemann Hypothesis by the Second Mean-Value Theorem of Calculus

Approach to a Proof of the Riemann Hypothesis by the Second Mean-Value Theorem of Calculus Advances in Pre Mathematics, 6, 6, 97- http://www.scirp.org/jornal/apm ISSN Online: 6-384 ISSN Print: 6-368 Approach to a Proof of the Riemann Hypothesis by the Second Mean-Vale Theorem of Calcls Alfred

More information

A new integral transform on time scales and its applications

A new integral transform on time scales and its applications Agwa et al. Advances in Difference Eqations 202, 202:60 http://www.advancesindifferenceeqations.com/content/202//60 RESEARCH Open Access A new integral transform on time scales and its applications Hassan

More information