Free vibration of antisymmetric,angle-ply laminated plates including transverse shear deformation by the finite element method

A finite element formulation of the equations governing the laminated anisotropic plate theory of Yang, Norris and Stavsky, is presented. The theory is a generalization of Mindlin's theory for isotropic plates to laminated anisotropic plates and includes shear deformation and rotary inertia effects. Finite element solutions are presented for rectangular plates of antisymmetric angle-ply laminates whose material properties are typical of a highly anisotropic composite material. Two sets of material properties that are typical of high modulus fiber-reinforced composites are used to show the parametric effects of plate aspect ratio, length-to-thickness ratio, number of layers and lamination angle. The numerical results are compared with the closed form results of Bert and Chen. As a special case, numerical results are presented for thick isotropic plates, and are compared with those for 3-D linear elasticity theory and Mindlin's thick plate theory.     

ON THE ROLES OF COMPLEMENTARY AND ADMISSIBLE BOUNDARY CONSTRAINTS IN FOURIER SOLUTIONS TO THE BOUNDARY VALUE PROBLEMS OF COMPLETELY COUPLED r TH ORDER PDES

A heretofore unavailable double Fourier series based approach, for obtaining non-separable solution to a system of completely coupled linear r th order partial differential equations with constant coefficients and subjected to general (completely coupled) boundary conditions, has been presented. The method has been successfully implemented to solve a class of hitherto unsolved boundary-value problems, pertaining to free and forced vibrations of arbitrarily laminated anisotropic doubly curved thin panels of rectangular planform, with arbitrarily prescribed (both symmetric and asymmetric with respect to the panel centerlines) admissible boundary conditions and subjected to general transverse loading.Existing solutions such as those due to Navier or Levy are based on the well-known method of separation of variables. Such solutions represent particular solutions whenever the method of separation of variables work, and when these particular solution functions fortuitously satisfy the boundary conditions. For derivation of the complementary solution, the complementary boundary constraints are introduced through boundary discontinuities of some of the particular solution functions and their partial derivatives. Such discontinuities form sets of measure zero.Various cases of lamination, geometry and dynamic response (forced and free vibrations) of a class of thin anisotropic laminated shells (curved panels) have been shown to follow from the above. Six sets of boundary conditions are used to illustrate the present method for the derivation of complementary solutions. Navier-type solutions whenever available form special cases of the present general solution.     

PHYSICAL AND NUMERICAL MODELLING OF THE DYNAMIC BEHAVIOR OF A FLY LINE

The planar equations of motion for a tapered fly line subjected to tension, bending, aerodynamic drag, and weight are derived. The resulting theory describes the large non-linear deformation of the line as it forms a propagating loop during fly casting. A cast is initiated by the motion of the tip of the fly rod that represents the boundary condition at one end of the fly line. At the opposite end, the boundary condition describes the equations of motion of a small attached fly (point mass with air drag). An efficient numerical algorithm is reviewed that captures the initiation and propagation of a non-linear wave that describes the loop. The algorithm is composed of three major steps. First, the non-linear initial-boundary-value problem is transformed into a two-point boundary-value problem, using finite differencing in time. The resulting non-linear boundary-value problem is linearized and then transformed into an initial-value problem in space. Example results are provided that illustrate how an overhead cast develops from initial conditions describing a perfectly laid out back cast. The numerical solutions are used to explore the influence of two sample effects in fly casting, namely, the drag created by the attached fly and the shape of the rod tip path.     

STRUCTURAL RESPONSE OF LAMINATED COMPOSITE SHELLS SUBJECTED TO BLAST LOADING: COMPARISON OF EXPERIMENTAL AND THEORETICAL METHODS

The results from a theoretical and experimental investigation of the dynamic response of cylindrically curved laminated composite shells subjected to normal blast loading are presented. The dynamic equations of motion for cylindrical laminated shells are derived using the assumptions of Love's theory of thin elastic shells. Kinematically admissible displacement functions are chosen to represent the motion of the clamped cylindrical shell and the governing equations are obtained in the time domain using the Galerkin method. The time-dependent equations of the cylindrically curved laminated shell are then solved by the Runge-Kutta-Verner method. Finite element modelling and analysis for the blast-loaded cylindrical shell are also presented. Experimental results for cylindrically curved laminated composite shells with clamped edges and subjected to blast loading are presented. The blast pressure and strain measurements are performed on the shell panels. The strain response frequencies of the clamped cylindrical shells subjected to blast load are obtained using the fast Fourier transformation technique. In addition, the effects of material properties on the dynamic behaviour are examined. The strain-time history curves show agreement between the experimental and analysis results in the longitudinal direction of the cylindrical panels. However, there is a discrepancy between the experimental and analysis results in the circumferential direction of the cylindrical panels. A good prediction is obtained for the response frequency of the cylindrical shell panels.     

静电场边值问题的矩量法解

本文采用修正格林函数积分方程的矩量法数值解技术求解静电场边值问题。这种方法具有较边界积分法积分域小的特点,不仅有高的计算效率和精度,尤其适用于有无限延伸导体边界的边值问题。本文着重讨论与传输线问题等效的边值问题,这是典型的二维静电场问题,对高频和微波技术有重要的实用价值。本文方法也适用于三维静电场问题和交变电磁场问题。     

非线性声学谐波方程的特解及其在边值问题中的应用

在流体或固体介质中,微扰法求解非线性声波的反射和折射问题时,谐波场通常满足非齐次波动方程。应用分离变量法及拉格朗日变动参数法求它的特解,出现了待定的分离常数,给这类非线性声学边值问题带来困难。本文结果表明,为不使定解问题出现佯僇,其独立的特解只有两类,一类是沿边界法线方向有积累效应的解,另一类则是沿平行边界面方向上有积累效应的解,由定解条件来决定究竟选用哪一类解。应用这个结果研究了平面边界的反射和折射谐波,该理论对非线性声学中的平面边值问题有普遍的应用意义。     

Vibration of bimodular sandwich beams with thick facings: A new theory and experimental results

This study deals with both analytical and experimental investigations of three-layer beams with cores of polyurethane foam and facings of unidirectional cord-rubber. Both of these materials are bimodular (i.e., having different behavior in compression as compared to tension). The new theory presented is a shear-flexible laminate version of the well-known Timoshenko beam theory, which, due to the bending-stretching coupling present in the bimodular case, results in a coupled sixth-order system of differential equations. In this theory, a separate derivation is presented for the shear correction factor. Due to the discontinuities in the normal stress distribution and the bimodularity, the shear correction factor is much different than the classical homogeneous material value of $\text{5}\text{6}$. Theoretical and experimental results are presented for the frequencies of the first three modes of vibration for a pin-ended beam without axial restraint. This work is believed to be the first devoted to vibration of bimodular materials in a sandwich configuration.     

The contribution of a lateral wave in simulating low-frequency sound fields in an irregular waveguide with a liquid bottom

A further development of a previously proposed approach to calculating the sound field in an arbitrarily irregular ocean is presented. The approach is based on solving the first-order causal mode equations, which are equivalent to the boundary-value problem for acoustic wave equations in terms of the cross-section method. For the mode functions depending on the horizontal coordinate, additional terms are introduced in the cross-section equations to allow for the multilayer structure of the medium. A numerical solution to the causal equations is sought using the fundamental matrix equation. For the modes of the discrete spectrum and two fixed low frequencies, calculations are performed for an irregular two-layer waveguide model with fluid sediments, which is close to the actual conditions of low-frequency sound propagation in the coastal zone of the oceanic shelf. The calculated propagation loss curves are used as an example for comparison with results that can be obtained for the given waveguide model with the use of adiabatic and one-way propagation approximations.     

Simulation of low-frequency sound propagation in an irregular shallow-water waveguide with a fluid bottom

A further development of a previously proposed approach to calculating the sound field in an arbitrarily irregular ocean is presented. The approach is based on solving the first-order causal mode equations, which are equivalent to the boundary-value problem for acoustic wave equations in terms of the cross-section method. For the mode functions depending on the horizontal coordinate, additional terms are introduced in the cross-section equations to allow for the multilayer structure of the medium. A numerical solution to the causal equations is sought using the fundamental matrix equation. For the modes of the discrete spectrum and two fixed low frequencies, calculations are performed for an irregular two-layer waveguide model with fluid sediments, which is close to the actual conditions of low-frequency sound propagation in the coastal zone of the oceanic shelf. The calculated propagation loss curves are used as an example for comparison with results that can be obtained for the given waveguide model with the use of adiabatic and one-way propagation approximations.     

Investigation of velocity field of gas flow moving behind a piston in rectangular branch pipe

The two-dimensional boundary-value problem of the unsteady flow of an incompressible viscous gas moving behind the piston in a “long” rectangular branch pipe is solved. An analytic solution is constructed for two velocity components with a refining polynomial, which reduces to a system of nonlinear algebraic equations after the substitution into the governing system of equations. By virtue of the solution uniqueness of the boundary-value problem under study, the only solution is found from obtained values of the refining polynomial constants for each point of the branch pipe internal space for the velocity components in analytic form.     

微波强耦合论

在本文中,作者尝试从耦合波的观点来研究电磁波通过两边有导电屏的长槽(即槽耦合波导系统)的衍射问题,提出一个实际求解这个边值问题的、建立在明皙物理概念上的理论。应用这个理论的观点,研究耦合波导问题,就和研究单一波导问题一样,能够采取相同的数学途径,这个途径就是正交函数的展开理论。这样,从来都是用不同方法来处理的两类波导传输问题,现在就能用统一的、联系的观点来分析。为了就明这个理论的具体应用,作者分析了用来完成矩形波导主波和低衰减圆电波之间功率转换的长槽定向耦合器,得出了一系列的计算基本参数的原始公式。文末指出了本文所提理论和方法的若干进一步应用。     

DYNAMICS BEHAVIOR OF AXISYMMETRIC AND BEAM-LIKE ANISOTROPIC CYLINDRICAL SHELLS CONVEYING FLUID

The flow-induced vibration characteristics of anisotropic laminated cylindrical shells partially or completely filled with liquid or subjected to a flowing fluid are studied in this work for two cases of circumferential wave number, the axisymmetric, where n=0 and the beam-like, where n=1. The shear deformation effects are taken into account in this theory; therefore, the equations of motion are determined with displacements and transverse shear as independent variables. The present method is a combination of finite element analysis and refined shell theory in which the displacement functions are derived from the exact solution of refined shell equations based on orthogonal curvilinear co-ordinates. Mass and stiffness matrices are determined by precise analytical integration. A finite element is defined for the liquid in cases of potential flow that yields three forces (inertial, centrifugal and Coriolis) of moving fluid. The mass, stiffness and damping matrices due to the fluid effect are obtained by an analytical integration of the fluid pressure over the liquid element. The available solution based on Sanders' theory can also be obtained from the present theory in the limiting case of infinite stiffness in transverse shear. The natural frequencies of isotropic and anisotropic cylindrical shells that are empty, partially or completely filled with liquid as well as subjected to a flowing fluid, are given. When these results are compared with corresponding results obtained using existing theories, very good agreement is obtained.     

A comparison of some shell theories used for the dynamic analysis of cross-ply laminated circular cylindrical panels

The free vibration problem of thin elastic cross-ply laminated circular cylindrical panels is considered. For this problem, a theoretical unification as well as a numerical comparison of the thin shell theories most commonly used (in engineering applications) is presented. In more detail, the problem is formulated in such a way that by using some tracers, which have the form of Kronecker's deltas, the stress-strain relations, constitutive equations and equations of motion obtained produce, as special cases, the corresponding relations and equations of Donnell's, Love's, Sanders' and Flugge's theories. By using a closed form solution, obtained for simply supported panels, a comparison of corresponding numerical results obtained on the basis of all of the aforementioned shell theories is attempted.     

Diffraction of elastic waves by a sub-surface crack (anti-plane motion)

A rigorous theory of the diffraction of SH-waves by a stress-free crack embedded in a semi-infinite elastic medium is presented. The incident time-harmonic SH-wave is taken to be either a uniform plane wave or a cylindrical wave originating from a surface line-source. The resulting boundary-value problem for the unknown jump in the particle displacement across the crack is solved by employing an integral equation approach. The unknown quantity is expanded in a complete sequence of Chebyshev polynomials. By writing the Green function as a Fourier integral, an infinite system of linear, algebraic equations for the expansion coefficients is obtained. Numerical results are presented for the particle displacement at the surface of the half-space, the far field radiation characteristic, the scattering cross-section of the crack and the dynamic stress intensity factor at the crack tips, for a range of geometrical parameters.     

Riemann-Hilbert problems for the Ernst equation and fibre bundles

Riemann-Hilbert techniques are used in the theory of completely integrable differential equations to generate solutions that contain a free function which can be used at least in principle to solve initial or boundary-value problems. The solution of a boundary-value problem is thus reduced to the identification of the jump data of the Riemann-Hilbert problem from the boundary data. But even if this can be achieved, it is very difficult to get explicit solutions since the matrix Riemann-Hilbert problem is equivalent to an integral equation. In the case of the Ernst equation (the stationary axisymmetric Einstein equations in vacuum), it was shown in a previous work that the matrix problem is gauge equivalent to a scalar problem on a Riemann surface. If the jump data of the original problem are rational functions, this surface will be compact which makes it possible to give explicit solutions in terms of hyperelliptic theta functions. In the present work, we discuss Riemann-Hilbert problems on Riemann surfaces in the framework of fibre bundles. This makes it possible to treat the compact and the non-compact case in the same setting and to apply general existence theorems.     

FREE VIBRATIONS OF MULTILAYERED PLATES BASED ON A MIXED VARIATIONAL APPROACH IN CONJUNCTION WITH GLOBAL PIECEWISE-SMOOTH FUNCTIONS

This paper presents a two-dimensional model for the analysis of freely vibrating laminated plates. The governing differential equations, associated boundary conditions and constitutive equations are derived from Reissner's mixed variational theorem. Both the governing differential equations and the related boundary conditions are presented in terms of resultant stresses and displacements. The model is able to provide the results for the corresponding three-dimensional theory. Such a performance is guaranteed from an appropriate expansion of relevant kinetic and stress quantities through the thickness of the multilayered plate. The expansion is realized by using a novel selection of global piecewise-smooth functions (GPSFs). The number of GPSFs can be arbitrarily increased to achieve a two-dimensional plate theory which is, at least, as accurate as that of a full layerwise theory. It is also shown that GPSFs permit to deal a multilayered plate as if it was virtually made of a single layer. Indeed, the theory need not explicitly introduce continuity conditions for both displacements and relevant stresses. The performance of the present two-dimensional model in conjunction with the global piecewise-smooth functions is tested and discussed by comparing its resulting eigen-parameters, for a class of simply supported plates, with those of other two-dimensional models and with those existing of the exact three-dimensional theory.     

Order parameter oscillations in a bounded solid solution and their bifurcations upon cooling

A nonlinear boundary-value problem for the Landau-Khalatnikov equation, which models the time evolution of the order parameter in a binary solution, is considered. A parameter characterizing the closeness of the alloy temperature to the critical temperature is introduced. This parameter depends on the dimensionless diffusivity. A change in this parameter results in successive bifurcations of solutions; so the order parameter corresponds to an oscillating steady-state structure to which almost all nonsteady-state solutions of the boundary-value problem are drawn after a long period of time.     

ELECTROSTATIC POTENTIAL OF STRONGLY NONLINEAR COMPOSITES: HOMOTOPY CONTINUATION APPROACH

The homotopy continuation method is used to solve the electrostatic boundary-value problems of strongly nonlinear composite media, which obey a current-field relation of J=σ E+χ|E|²E. With the mode expansion, the approximate analytical solutions of electric potential in host and inclusion regions are obtained by solving a set of nonlinear ordinary different equations, which are derived from the original equations with homotopy method. As an example in dimension two, we apply the method to deal with a nonlinear cylindrical inclusion embedded in a host. Comparing the approximate analytical solution of the potential obtained by homotopy method with that of numerical method, we can obverse that the homotopy method is valid for solving boundary-value problems of weakly and strongly nonlinear media.     

用边界元法分析非均匀介质中的传输线

 从静电场边值问题的积分解出发,推导出用边界元法求解分区均匀介质填充传输线问题的矩阵表达式,给出传输线电容参数的计算公式,介绍用边界元法求解分区均匀介质填充传输线问题的基本原理和求解过程。对两类传输线的计算结果表明：用边界元法求解分区均匀介质填充传输线问题,不仅具有较高的计算精度,而且可以很方便地应用于各类复杂截面分区均匀介质填充传输线问题的工程设计与计算,边界元法是求解分区均匀介质填充传输线问题的一种有效方法。     

AN EXACT METHOD FOR THE FREE VIBRATION ANALYSIS OF TIMOSHENKO-KELVIN BEAMS WITH OSCILLATORS

In this paper, a modern exact method is proposed for solving the problem of free vibrations of a Timoshenko-type viscoelastic beam with discrete rigid bodies, connected to the beam by means of viscoelastic constraints. The phenomenon of free vibrations of this discrete-continuous system is described by a set of three partial and two subsystem ordinary differential equations with generalized boundary conditions and initial conditions. Vector notation of the equations allows one to identify the self-adjoint linear operators of inertia, stiffness and damping. In this case, these operators are not homothetic hence a separation of variables in this set of equations is possible only in a complex Hilbert space. Such separation of variables leads to ordinary differential equations of motion with respect to time as well as to a set of three ordinary differential equations with respect to a spatial variable and two subsystem algebraical equations. Solution of the boundary-value problem was carried out in the classical way, but its results are complex conjugated. Using these results and the fundamental principle, describing the orthogonality property of complex eigenvectors, the problem of free vibrations of the system with arbitrary initial conditions has been finally solved exactly.