A SEMI-ANALYTICAL AND SEMI-NUMERICAL METHOD FOR SOLVING 2-D SOUND-STRUCTURE INTERACTION PROBLEMS

Based on the transfer matrix method and the virtual source simulation technique, this paper proposes a novel semi-analytical and semi-numerical method for solving 2-D sound- structure interaction problems under a harmonic excitation.Within any integration segment, as long as its length is small enough,along the circumferential curvilinear coordinate,the non- homogeneous matrix differential equation of an elastic ring of complex geometrical shape can be rewritten in terms of the homogeneous one by the method of extended homogeneous capacity proposed in this paper.For the exterior fluid domain,the multi-circular virtual source simulation technique is adopted.The source density distributed on each virtual circular curve may be ex- panded as the Fourier's series.Combining with the inverse fast Fourier transformation,a higher accuracy and efficiency method for solving 2-D exterior Helmholtz's problems is presented in this paper.In the aspect of solution to the coupling equations,the state vectors of elastic ring induced by the given harmonic excitation and generalized forces of coefficients of the Fourier series can be obtained respectively by using a high precision integration scheme combined with the method of extended homogeneous capacity put forward in this paper.According to the superposition princi- ple and compatibility conditions at the interface between the elastic ring and fluid,the algebraic equation of system can be directly constructed by using the least square approximation.Examples of acoustic radiation from two typical fluid-loaded elastic rings under a harmonic concentrated force are presented.Numerical results show that the method proposed is more efficient than the mixed FE-BE method in common use.

SEMI-ANALYTICAL AND SEMI-NUMERICAL METHOD FOR DYNAMIC ANALYSIS OF FOUNDATION

A semi-analytical and semi-numerical method is proposed for the dynamic analysis of foundations. The Lamb＇s solution and the approximate formulae were used to establish the relation of the contact force and deflection between the foundation and soil. Therefore, the foundation can be separated from soil and analyzed by FEM as for the static cases. The plate can be treated as that the known forces are acting on the upper surface, and the contact pressure from soil can be represented as the deflection. So that only the plate needs to be divided into elements in the analysis. By this method, a series of vibration problems, including various shapes and rigidities of foundations, different excitation frequencies, were analyzed. Furthermore, it can be used for the embedded foundation. The numerical examples show that this method has simplicity, highly accurate and versatile. It is an effective method for the dynamic analysis of foundations.

Semi-analytical solution for internal forces of tunnel lining with multiple longitudinal cracks

Longitudinal cracks on the tunnel lining significantly influence the performance of tunnels in operation.In this study,we propose a semi-analytical method that provides a simple and effective way to calculate the internal forces of tunnel linings with multiple cracks.The semi-analytical solution is obtained using structural analysis considering the flexural rigidity for the cracked longitudinal section of the tunnel lining.Then the proposed solution is verified numerically.Using the proposed method,the influences of the crack depth and the number of cracks on the bending moment and modified crack tip stress are investigated.With the increase in crack depth,the bending moment of lining scetion adjacent to the crack decreases,while the bending moment of lining scetion far away from the crack increases slightly.The more the number of cracks in a tunnel lining,the easier the new cracks initiated.

A NOVEL SEMI-ANALYTICAL METHOD FOR SOLVING ACOUSTIC RADIATION FROM LONGITUDINALLY STIFFENED INFINITE NON-CIRCULAR CYLINDRICAL SHELLS IN WATER

Based on the extended homogeneous capacity high precision integration method and the spectrum method of virtual boundary with a complex radius vector, a novel semi-analytical method, which has satisfactory computation efectiveness and precision, is presented for solving the acoustic radiation from a submerged infnite non-circular cylindrical shell stifened by longitudinal ribs by means of the Fourier integral transformation and stationary phase method. In this work, besides the normal interacting force, which is commonly adopted by some researchers, the other interacting forces and moments between the longitudinal ribs and the non-circular cylindrical shell are considered at the same time. The efects of the number and the size of the cross-section of longitudinal ribs on the characteristics of acoustic radiation are investigated. Numerical results show that the method proposed is more efcient than the existing mixed FE-BE method.

Nonlinear dynamic response and buckling of laminated cylindrical shells with axial shallow groove based on a semi-analytical method

The effect of axial shallow groove on the nonlinear dynamic response and buckling of laminated cylindrical shells subjected to radial compression loading was investigated. Based on the first-order shear deformation theory （FSDT）, the nonlinear dynamic equations involving the transverse shear deformation and initial geometric imperfections were derived with the Hamilton philosophy. The axial shallow groove of the laminated composite cylindrical shell was treated as the initial geometric imperfections in the dynamic equations. A semi-analytical method of expanding displacements and loads along the circumferential direction and employing the finite difference method along the axial direction and in the time domain is used to solve the governing equations and obtain the dynamic response of the laminated shell. The B-R criterion was employed to determine the critical loads of dynamic buckling of the shell. The effects of the parameters of the shallow groove on the dynamic response and buckling were discussed in this paper and the results show that the axial shallow grooves greatly affect the dynamic response and buckling.

Numerical Study on Motion Responses of Two Parallel Ships During Underway Replenishment by the Semi-Analytical HOTP Method

Motion responses of two ships advancing parallel in waves with hydrodynamic interactions are investigated in this paper. Within the framework of the frequency-domain potential flow theory, a semi-analytical higher-order translating-pulsating source(HOTP) method is presented to solve the problems of coupled radiation and diffraction potential. The method employs nine-node bi-quadratic curvilinear elements to discretize the boundary integral equations(BIEs) constructed over the mean wetted surface of the two ship hulls. In order to eliminate the numerical oscillation, analytical quadrature formulas are derived and adopted to evaluate the integrals related to the Froudedependent part of the Green’s function along the horizontal direction in the BIEs. Based on the method, a numerical program is originally coded. Through the calculations of hydrodynamic responses of single ships, the numerical implementation is proved successful. Then the validated program is applied in the investigations on the hydrodynamic interactions of two identical Wigley Ⅲ hulls and the underway replenishment of a frigate and a supply ship in waves with and without stagger, respectively. The comparison between the present computed results with experimental data and numerical solutions of other methods shows that the semi-analytical HOTP method is of higher accuracy than the pulsating source Green’s function method with speed correction and better stability than the traditional HOTP method based on Gauss quadrature. In addition, for two ships with obviously different dimensions,the influence of hydrodynamic interactions on the smaller ship is found to be more noticeable than that on the larger ship, which leads to the differences between the motions of frigate with and without the presence of supply ship.

SEMI-ANALYTICAL FINITE ELEMENT METHOD FOR FICTITIOUS CRACK MODEL IN FRACTURE MECHANICS OF CONCRETE

Based on the Hamiltonian governing equations of plane elasticity for sectorial domain, the variable separation and eigenfunction expansion techniques were employed to develop a novel analytical finite element for the fictitious crack model in fracture mechanics of concrete. The new analytical element can be implemented into FEM program systems to solve fictitious crack propagation problems for concrete cracked plates with arbitrary shapes and loads. Numerical results indicate that the method is more efficient and accurate than ordinary finite element method.

FORMULATION OF A SEMI-ANALYTICAL APPROACH BASEDON GURTIN VARIATIONAL PRINCIPLE FOR DYNAMICRESPONSE OF GENERAL THIN PLATES

4 semi-analytical approach for the dynamic response of general thin plates whichemployes finite element discretization in space domain and a series of representation intime demain is developed on the basis of Curtin variational principles.The formulationof time series is also investigated so that the dynamic response of plates with arbitraryshape and boundary constraints can be achieved with adequate accuracy.

Substructure Method of Three-Dimensional Semi-Analytic Boundary Element and Its Application to Seismic Scattering of Hill Topography

In this paper,a substructure method of three-dimensional semi-analytic boundary element is established.The seismic scattering by three-dimensional topography of a hill can be analyzed by the method in frequency domain.Using this method,the computational effort and storage space are reduced considerably.Finally,analytic results are given.

ELASTO-PLASTIC COUPLED ANALYSIS OF BURIED STRUCTURE AND SOIL MEDIUM BY PERTURBATIONAL SEMI-ANALYTIC METHOD

In this paper, an effective numerical method for physically nonlinear interaction analysis is studied, in which the elasto-plastic problem of coupled analysis between the structure and medium may be transformed into several linear problems by means of the perturbation technique, then, the finite strip method and finite layer method are used to analyze the underground structure and rock medium, respectively, for their corresponding linear problems, so the purpose of simplifying the calculation can be achieved. This kind of method has made use of the twice semi-analytical technique: the perturbation and semi-analytic solution function to simplify 3-D nonlinear coupled problem into 1-D linear numerical one. In addition, this method is a new advance of semi-analytical method in the application to nonlinear problems by means of combinating with the analytical perturbation method, and it is also a branch of the perturbational numerical method developed in last years.

Axisymmetric Finite Element Method for Analysis of Plate on Layered Soil

To obtain the fundamental solution of soil has become the key problem for the semi-analytical and semi-numerical (SASN) method in analyzing plate on layered soil. By applying axisymmetric finite element method (FEM),an expression relating the surface settlement and the reaction of the layered soil can be obtained. Such a reaction can be treated as load acting on the applied external load. Having the plate modelled by four-node elements,the governing equation of the plate can be formed and solved. In this case, the fundamental solution can be introduced into the global soil stiffness matrix and five-node or nine-node element soil stiffness matrix.The existing commercial FEM software can be used to solve the fundamental solution of soil, which can bypass the complicated formula derivation and boasts high computational efficiency as well.

Two-dimensional fracture analysis of piezoelectric material based on the scaled boundary node method

A scaled boundary node method （SBNM） is developed for two-dimensional fracture analysis of piezoelectric material, which allows the stress and electric displacement intensity factors to be calculated directly and accurately. As a boundary- type meshless method, the SBNM employs the moving Kriging （MK） interpolation technique to an approximate unknown field in the circumferential direction and therefore only a set of scattered nodes are required to discretize the boundary. As the shape functions satisfy Kronecker delta property, no special techniques are required to impose the essential boundary conditions. In the radial direction, the SBNM seeks analytical solutions by making use of analytical techniques available to solve ordinary differential equations. Numerical examples are investigated and satisfactory solutions are obtained, which validates the accuracy and simplicity of the proposed approach.

Thermal field in water pipe cooling concrete hydrostructures simulated with singular boundary method

The embedded water pipe system is often used as a standard cooling technique during the construction of large-scale mass concrete hydrostructures. The prediction of the temperature distribution considering the cooling effects of embedded pipes plays an essential role in the design of the structure and its cooling system. In this study, the singular boundary method, a semi-analytical meshless technique, was employed to analyze the temperature distribution. A numerical algorithm solved the transient temperature field with consideration of the effects of cooling pipe specification, isolation of heat of hydration, and ambient temperature. Numerical results are verified through comparison with those of the finite element method, demonstrating that the proposed approach is accurate in the simulation of the thermal field in concrete structures with a water cooling pipe.

Semi-analytic Solution of Steady Temperature Fields During Thin Plate Welding

Usually, it is very difficult to find out an analytical solution to thermal conduction problems during high temperature welding. Therefore, as an important numerical approach, the method of lines （MOLs） is introduced to solve the temperature field characterized by high gradients. The basic idea of the method is to semi-discretize the governing equation of the problem into a system of ordinary differential equations （ODEs） defined on discrete lines by means of the finite difference method, by which the thermal boundary condition with high gradients are directly embodied in formulation. Thus the temperature field can be obtained by solving the ODEs. As a numerical example, the variation of an axisymmetrical temperature field along the plate thickness can be obtained.

Study of the natural vibration characteristics of water motion in the moon pool by the semi-analytical method

The three-dimensional natural vibration characteristics of water inside a moon pool of an ocean structures are studied. The governing equations are derived based on the linear potential flow theory, and the boundary condition of the total opening bottom suggested by Molin is adopted. A semi-analytical method is used to solve the governing equations, and the natural frequencies and the motion modes are obtained. Two types of motions are studied： （1） the piston motion in the vertical direction, and （2） the sloshing motion of the free surface. The influences of moon pool＇s structural parameters on the natural frequencies, and the modal shapes are analyzed.

Application of semi-analytical finite element method to analyze the bearing capacity of asphalt pavements under moving loads

To facilitate long term infrastructure asset management systems, it is necessary to determine the bearing capacity of pavements. Currently it is common to conduct such measurements in a stationary manner, however the evaluation with stationary loading does not correspond to reality a tendency towards continuous and high speed measurements in recent years can be observed. The computational program SAFEM was developed with the objective of evaluating the dynamic response of asphalt under moving loads and is based on a semi-analytic element method. In this research project SAFEM is compared to commercial finite element software ABAQUS and field measurements to verify the computational accuracy. The computational accuracy of SAFEM was found to be high enough to be viable whilst boasting a computational time far shorter than ABAQUS. Thus, SAFEM appears to be a feasible approach to determine the dynamic response of pavements under dynamic loads and is a useful tool for infrastructure administrations to analyze the pavement bearing capacity.

Application of semi-analytical finite element method coupled with infinite element for analysis of asphalt pavement structural response

A specific computational program SAFEM was developed based on semi-analytical finite element （FE） method for analysis of asphalt pavement structural responses under static loads. The reliability and efficiency of this FE program was proved by comparison with the general commercial FE software ABAQUS. In order to further reduce the computational time without decrease of the accuracy, the infinite element was added to this program. The results of the finite-infinite element coupling analysis were compared with those of finite element analysis derived from the verified FE program, The study shows that finite-infinite element coupling analysis has higher reliability and efficiency.

A semi-analytical method for forced vibration analysis of cracked laminated composite beam with general boundary condition

In this paper,a semi-analytical method for the forced vibration analysis of cracked laminated composite beam(CLCB)is investigated.One computational model is formulated by Timoshenko beam theory and its dynamic solution is solved using the Jacobi-Ritz method.The boundary conditions(BCs)at both ends of the CLCB are generalized by the application of artificial elastic springs,the CLCB is separated into two elements along the crack,the flexibility coefficient of fracture theory is used to model the essential continuous condition of the connective interface.All the allowable displacement functions used to analyze dynamic characteristics of CLCB are expressed by classical Jacobi orthogonal polynomials in a more general form.The accuracy of the proposed method is verified through the compare with results of the finite element method(software ABAQUS is used in this paper).On this basis,the parametric study for dynamic analysis characteristics of CLCB is performed to provide reference datum for engineers.

Numerical study on butterfly wings around inclusion based on damage evolution and semi-analytical method

Butterfly wings are closely related to the premature failure of rolling element bearings.In this study,butterfly formation is investigated using the developed semi-analytical three-dimensional(3D)contact model incorporating inclusion and material property degradation.The 3D elastic field introduced by inhomogeneous inclusion is solved by using numerical approaches,which include the equivalent inclusion method(EIM)and the conjugate gradient method(CGM).The accumulation of fatigue damage surrounding inclusions is described using continuum damage mechanics.The coupling between the development of the damaged zone and the stress field is considered.The effects of the inclusion properties on the contact status and butterfly formation are discussed in detail.The model provides a potential method for quantifying material defects and fatigue behavior in terms of the deterioration of material properties.

Efficient semi-analytical methods for the vibration response of and acoustic radiation from a periodical orthogonally rib-stiffened plate

In order to theoretically predict and analyze the vibration response and acoustic radiation characteristics of a periodical orthogonally rib-stiffened plate,its vibro-acoustic equations of an underwater infinite model are established.The rib-stiffened plate is stimulated by a harmonic plane pressure.By using the Fourier transforms,Poisson’s summation formula and space harmonic method,the structural vibration response and acoustic radiation pressure are expressed as functions of displacement harmonic components.Efficient semi-analytical methods are proposed in this work,and then approximate solutions for finite terms of the harmonic components are obtained by employing the truncation technique.Effects of the vibration response,rib spacing and torsional moment of the ribs on the radiation pressure are examined,and the validity of the present methods is also verified.Theoretical results show that the torsional moment of the ribs affects the modal frequencies of the stiffened plate,which should not be neglected in engineering applications with high precision requirement.With attachment of the ribs to the thin plate,its far field radiation pressure can be reduced in the low frequency range by adjusting rib spacing and cross sectional size of the ribs.