Hamiltonian system for the inhomogeneous plane elasticity of dodecagonal quasicrystal plates and its analytical solutions

A Hamiltonian system is derived for the plane elasticity problem of two-dimensional dodecagonal quasicrystals by introducing the simple state function. By using symplectic elasticity approach, the analytic solutions of the phonon and phason displacements are obtained further for the quasicrystal plates. In addition, the effectiveness of the approach is verified by comparison with the data of the finite integral transformation method.

Analytical solution for the effective elastic properties of rocks with the tilted penny-shaped cracks in the transversely isotropic background

Seismic prediction of cracks is of great significance in many disciplines,for which the rock physics model is indispensable.However,up to now,multitudinous analytical models focus primarily on the cracked rock with the isotropic background,while the explicit model for the cracked rock with the anisotropic background is rarely investigated in spite of such case being often encountered in the earth.Hence,we first studied dependences of the crack opening displacement tensors on the crack dip angle in the coordinate systems formed by symmetry planes of the crack and the background anisotropy,respectively,by forty groups of numerical experiments.Based on the conclusion from the experiments,the analytical solution was derived for the effective elastic properties of the rock with the inclined penny-shaped cracks in the transversely isotropic background.Further,we comprehensively analyzed,according to the developed model,effects of the crack dip angle,background anisotropy,filling fluid and crack density on the effective elastic properties of the cracked rock.The analysis results indicate that the dip angle and background anisotropy can significantly either enhance or weaken the anisotropy degrees of the P-and SH-wave velocities,whereas they have relatively small effects on the SV-wave velocity anisotropy.Moreover,the filling fluid can increase the stiffness coefficients related to the compressional modulus by reducing crack compliance parameters,while its effects on shear coefficients depend on the crack dip angle.The increasing crack density reduces velocities of the dry rock,and decreasing rates of the velocities are affected by the crack dip angle.By comparing with exact numerical results and experimental data,it was demonstrated that the proposed model can achieve high-precision estimations of stiffness coefficients.Moreover,the assumption of the weakly anisotropic background results in the consistency between the proposed model and Hudson's published theory for the orthorhombic rock.

Analytical solutions to the precession relaxation of magnetization with uniaxial anisotropy

Based on the Landau-Lifshitz-Gilbert(LLG)equation,the precession relaxation of magnetization is studied when the external field H is parallel to the uniaxial anisotropic field H_(k).The evolution of three-component magnetization is solved analytically under the condition of H=nH_(k)(n=3,1 and 0).It is found that with an increase of H or a decrease of the initial polar angle of magnetization,the relaxation time decreases and the angular frequency of magnetization increases.For comparison,the analytical solution for H_(k)=0 is also given.When the magnetization becomes stable,the angular frequency is proportional to the total effective field acting on the magnetization.The analytical solutions are not only conducive to the understanding of the precession relaxation of magnetization,but also can be used as a standard model to test the numerical calculation of LLG equation.

An improved analytical solution for solving the shield tunnel uplift problem

The problem of shield tunnel uplift is a common issue in tunnel construction.Due to the decrease in shear stiffness at the joints between the rings,uplift is typically observed as bending and dislocation deformation at these joints.Existing modeling methods typically rely on the Euler-Bernoulli beam theory,only considering the bending effect while disregarding shear deformation.Furthermore,the constraints on the shield tail are often neglected in existing models.In this study,an improved theoretical model of tunnel floating is proposed.The constraint effect of the shield machine shell on the tunnel structure is considered using the structural forms of two finite long beams and one semi-infinite long beam.Furthermore,the Timoshenko beam theory is adopted,providing a more accurate description of tunnel deformation,including both the bending effect and shear deformation,than existing models.Meanwhile,the buoyancy force and stratum resistance are calculated in a nonlinear manner.A reliable method for calculating the shear stiffness correction factor is proposed to better determination of the calculation parameters.The proposed theoretical model is validated through five cases using sitemonitored data.Its applicability and effectiveness are demonstrated.Furthermore,the influences of soil type,buried depth,and buoyancy force on the three key indicators of tunnel floating(i.e.the maximum uplift magnitude,the ring position with the fastest uplift race,and the ring position with the maximum uplift magnitude)are analyzed.The results indicate that the proposed model can provide a better understanding of the floating characteristics of the tunnel structure during construction.

Bifurcations, Analytical and Non-Analytical Traveling Wave Solutions of (2 + 1)-Dimensional Nonlinear Dispersive Boussinesq Equation

For the (2 + 1)-dimensional nonlinear dispersive Boussinesq equation, by using the bifurcation theory of planar dynamical systems to study its corresponding traveling wave system, the bifurcations and phase portraits of the regular system are obtained. Under different parametric conditions, various sufficient conditions to guarantee the existence of analytical and non-analytical solutions of the singular system are given by using singular traveling wave theory. For certain special cases, some explicit and exact parametric representations of traveling wave solutions are derived such as analytical periodic waves and non-analytical periodic cusp waves. Further, two-dimensional wave plots of analytical periodic solutions and non-analytical periodic cusp wave solutions are drawn to visualize the dynamics of the equation.

Analytical solutions fractional order partial differential equations arising in fluid dynamics

This article describes the solution procedure of the fractional Pade-Ⅱ equation and generalized Zakharov equation(GSEs)using the sine-cosine method.Pade-Ⅱ is an important nonlinear wave equation modeling unidirectional propagation of long-wave in dispersive media and GSEs are used to model the interaction between one-dimensional high,and low-frequency waves.Classes of trigonometric and hyperbolic function solutions in fractional calculus are discussed.Graphical simulations of the numerical solutions are flaunted by MATLAB.

Analytical Solutions for Testing of Baroclinic and Wind-Driven Three-Dimensional Hydrodynamic Models

It is always a challenge for a model developer to verify a three-dimensional hydrodynamic model, especially for the baroclinic term over variable topography, due to a lack of observational data sets or suitable analytical solutions. In this paper, exact solutions for the periodic forcing by surface heat flux and wind stress are given by solving the linearized equations of motion neglecting the rotation, advection and horizontal diffusion terms. The temperature at the bottom is set to a prescribed periodic value and a slip condition on flow is enforced at the bottom. The geometry of the quarter annulus, which has been extensively studied for two- and three-dimensional analytical solutions of unstratified water bodies, is used with a general power law variation of the bottom slope in the radial direction and is constant in the azimuthal direction. The analytical solutions are derived in a cylindrical coordinate system, which describes the three-dimensional fluid field in a Cartesian coordinate system. The results presented in this paper should provide a foundation for studying and verifying the baroclinic term over a varied topography in a three-dimensional numerical model.

Three-dimensional analytical solution of acoustic emission source location for cuboid monitoring network without pre-measured wave velocity

To find analytical solutions of nonlinear systems for locating the acoustic emission/microseismic（AE/MS） source without knowing the wave velocity of structures, the sensor location coordinates were simplified as a cuboid monitoring network. Different locations of sensors on upper and lower surfaces were considered and used to establish nonlinear equations. Based on the proposed functions of time difference of arrivals, the analytical solutions were obtained using five sensors under three networks. The proposed analytical solutions were validated using authentic data of numerical tests and experiments. The results show that located results are consistent with authentic data, and the outstanding characteristics of the new solution are that the solved process is not influenced by the wave velocity knowledge and iterated algorithms.

Three-dimensional analytical solution of acoustic emission or microseismic source location under cube monitoring network

To find the analytical solution of the acoustic emission/microseismic(AE/MS) source location coordinates, the sensor location coordinates were optimized and simplified. A cube monitoring network of sensor location was selected, and the AE/MS source localization equations were established. A location method with P-wave velocity by analytical solutions (P-VAS) was obtained with these equations. The virtual location tests show that the relocation results of analytical method are fully consistent with the actual coordinates for events both inside and outside the monitoring network; whereas the location error of traditional time difference method is between 0.01 and 0.03 m for events inside the sensor array, and the location errors are larger, which is up to 1080986 m for events outside the sensor array. The broken pencil location tests were carried out in the cross section of 100 mm×98 mm, 350 mm-length granite rock specimen using five AE sensors. Five AE sources were relocated with the conventional method and the P-VAS method. For the four events outside monitoring network, the positioning accuracy by P-VAS method is higher than that by the traditional method, and the location accuracy of the larger one can be increased by 17.61 mm. The results of both virtual and broken pencil location tests show that the proposed analytical solution is effective to improve the positioning accuracy. It can locate the coordinates of AE/MS source only using simple four arithmetic operations, without determining the fitting initial value and iterative calculation, which can be solved by a conventional calculator or Microsoft Excel.

ANALYTICAL SOLUTION OF AXIAL SYMMETRICAL WALL TEMPERATURE DISTRIBUTION FOR HOLLOW CYLINDER

Using the variable transformation method,the formulae of the axial symmetrical wall temperature distribution during steady heat conduction of a hollow cylinder are derived in this paper.The wall temperature distribution and the wall heat flux distribution in both axial and radial direction can be calculated by the temperature distribution of the liquid medium both inside and outside the cylinder with temperature changing in axial direction.The calculation results are almost consistent with the experience results.The applicative condition of the formulae in this paper consists with most of practice.They can be applied to the engineering calculation of the steady heat conduction.The calculation is simple and accurate.

Scattering of plane P waves by circular-arc layered alluvial valleys: An analytical solution

An analytical solution for scattering of plane P waves by circular-arc layered alluvial valleys was derived by Fourier-Bessel series expansion technique, and the solution was utilized to analyze the effects of alluvial sequence and their relative stiffness on the scattering of incident waves.

Analytic Solutions for Hindered Rotation Quantum Problem

Normalizable analytic solutions of the quantum rotor problem with divergent potential are presented here as solution of the Schrödinger equation. These solutions, unknown to the literature, represent a mathematical advance in the description of physical phenomena described by the second derivative operator associated with a divergent interaction potential and, being analytical, guarantee the optimal interpretation of such phenomena.

Analytical solutions for characteristic radii of circular roadway surrounding rock plastic zone and their application

In the non-uniform stress field, the surrounding rock plastic zone of the circular roadway shows different shapes under the different confining pressure conditions. Based on the boundary shape characteristics of the plastic zone, the characteristic radii of the plastic zone were proposed, namely the horizontal,longitudinal and medial axis radii, which could reflect the plastic zone shapes characteristics and classify the sizes of the key parts. On the theoretical basis of elastic-plastic mechanics, analytical solutions for the characteristic radii were obtained by theoretical deduction, and the relationships between the characteristic radii and key influencing factors were analyzed. Finally, the evaluation criterion of the circular roadway surrounding rock plastic zone shapes, evaluation criterion of the location of potential hazards caused by the roadway surrounding rock and evaluation critical points of roadway dynamic disasters based on characteristic radii were proposed. This work could provide a theoretical basis for stability analysis of the surrounding rock, support design, and guide the prevention and control of dynamic roadway disasters.

Unified analytical solution for deep circular tunnel with consideration of seepage pressure,grouting and lining

A new unified analytical solution is presented for predicting the range of plastic zone and stress distributions around a deep circular tunnel in a homogeneous isotropic continuous medium. The rock mass, grouting zone and lining are assumed as elastic-perfectly plastic and governed by the unified strength theory(UST). This new solution has made it possible to consider the interaction between seepage pressure, lining, grouting and rock mass, and the intermediate principal stress effect together. Moreover, parametric analysis is carried out to identify the influence of the related parameters on the plastic zone radius. Under the given conditions, the results show that the plastic zone radius decreases with an increasing cohesion, internal friction angle and hydraulic conductivity of lining and unified failure criterion parameter, respectively; whereas the plastic zone radius increases with the growth of elasticity modulus of lining. Comparison of results from the new solution and the other published one shows well agreement and provides confidence in the new solution proposed.

Scattering of plane P waves by a semi-cylindrical hill: analytical solution

An analytical solution for scattering of plane P waves by a semi-cylindrical hill was derived by using the wave function expansion method, and convergence of the solution and accuracy of truncation were verified. The effect of incident frequency and incident angle on the surface motion of the hill was discussed, and it was shown that a hill greatly amplifies incident plane P waves, and maximum horizontal displacement amplitudes appear mostly at the inclined incidence of waves, which are located at the half-space; and maximum vertical displacement amplitudes emerge mostly at the vertical incidence of waves, which are situated at the hill.

Symplectic system based analytical solution for bending of rectangular orthotropic plates on Winkler elastic foundation

This paper analyses the bending of rectangular orthotropic plates on a Winkler elastic foundation.Appropriate definition of symplectic inner product and symplectic space formed by generalized displacements establish dual variables and dual equations in the symplectic space.The operator matrix of the equation set is proven to be a Hamilton operator matrix.Separation of variables and eigenfunction expansion creates a basis for analyzing the bending of rectangular orthotropic plates on Winkler elastic foundation and obtaining solutions for plates having any boundary condition.There is discussion of symplectic eigenvalue problems of orthotropic plates under two typical boundary conditions,with opposite sides simply supported and opposite sides clamped.Transcendental equations of eigenvalues and symplectic eigenvectors in analytical form given.Analytical solutions using two examples are presented to show the use of the new methods described in this paper.To verify the accuracy and convergence,a fully simply supported plate that is fully and simply supported under uniformly distributed load is used to compare the classical Navier method,the Levy method and the new method.Results show that the new technique has good accuracy and better convergence speed than other methods,especially in relation to internal forces.A fully clamped rectangular plate on Winkler foundation is solved to validate application of the new methods,with solutions compared to those produced by the Galerkin method.

Two dimensional analytical solution for a partially vegetated compound channel flow

The theory of an eddy viscosity model is applied to the study of the flow in a compound channel which is partially vegetated. The governing equation is constituted by analyzing the longitudinal forces acting on the unit volume where the effect of the vegetation on the flow is considered as a drag force item, The compound channel is divided into 3 sub-regions in the transverse direction, and the coefficients in every region＇s differential equations were solved simultaneously. Thus, the analytical solution of the transverse distribution of the depth-averaged velocity for uniform flow in a partially vegetated compound channel was obtained. The results can be used to predict the transverse distribution of bed shear stress, which has an important effect on the transportation of sediment. By comparing the analytical results with the measured data, the analytical solution in this paper is shown to be sufficiently accurate to predict most hydraulic features for engineering design purposes.

Analytical solution to one-dimensional consolidation in unsaturated soils

This paper presents an analytical solution of the one-dimensional consolidation in unsaturated soil with a finite thickness under vertical loading and confinements in the lateral directions. The boundary contains the top surface permeable to water and air and the bottom impermeable to water and air. The analytical solution is for Fredlund＇s one-dimensional consolidation equation in unsaturated soils. The transfer relationship between the state vectors at top surface and any depth is obtained by using the Laplace transform and Cayley-Hamilton mathematical methods to the governing equations of water and air, Darcy＇s law and Fick＇s law. Excess pore-air pressure, excess pore-water pressure and settlement in the Laplace-transformed domain are obtained by using the Laplace transform with the initial conditions and boundary conditions. By performing inverse Laplace transforms, the analytical solutions are obtained in the time domain. A typical example illustrates the consolidation characteristics of unsaturated soil from an- alytical results. Finally, comparisons between the analytical solutions and results of the finite difference method indicate that the analytical solution is correct.

An Analytical Solution for Advance of Solute Front in Soils

Based on the assumption that solute transport in a semi-infinite soil columnor in a field soil profile can be described by the boundary-layer method, an analytical solution ispresented for the advance of a solute front with time. The traditional convection-dispersionequation (CDE) subjected to two boundary conditions: 1) at the soil surface (or inlet boundary) and2) at the solute front, was solved using a Laplace transformation. A comparison of residentconcentrations using a boundary-layer method and an exact solution (in a semi-infinite-domain)showed that both were in good agreement within the range between the two boundaries. This led to anew method for estimating solute transport parameters in soils, requiring only observation ofadvance of the solute front with time. This may be corroborated visually using a tracer solutionwith marking-dye or measured utilizing time domain reflectometry (TDR). This method is applicable toboth laboratory soil columns and field soils. Thus, it could be a step forward for modeling solutetransport in field soils and for better understanding of the transport processes in soils.

ANALYTICAL SOLUTIONS TO STRESS CONCENTRATION PROBLEM IN PLATES CONTAINING RECTANGULAR HOLE UNDER BIAXIAL TENSIONS

The stress concentration problem in structures with a circular or elliptic hole can be investigated by analytical methods. For the problem with a rectangular hole, only approximate results are derived. This paper deduces the analytical solutions to the stress concentration problem in plates with a rectangular hole under biaxial tensions. By using the U-transformation technique and the finite element method, the analytical displacement solutions of the finite element equations are derived in the series form. Therefore, the stress concentration can then be discussed easily and conveniently. For plate problem the bilinear rectangular element with four nodes is taken as an example to demonstrate the applicability of the proposed method. The stress concentration factors for various ratios of height to width of the hole are obtained.