Semi-analytical finite element method applied for characterizing micropolar fibrous composites

A semi-analytical finite element method(SAFEM),based on the two-scale asymptotic homogenization method(AHM)and the finite element method(FEM),is implemented to obtain the effective properties of two-phase fiber-reinforced composites(FRCs).The fibers are periodically distributed and unidirectionally aligned in a homogeneous matrix.This framework addresses the static linear elastic micropolar problem through partial differential equations,subject to boundary conditions and perfect interface contact conditions.The mathematical formulation of the local problems and the effective coefficients are presented by the AHM.The local problems obtained from the AHM are solved by the FEM,which is denoted as the SAFEM.The numerical results are provided,and the accuracy of the solutions is analyzed,indicating that the formulas and results obtained with the SAFEM may serve as the reference points for validating the outcomes of experimental and numerical computations.

HAMILTONIAN SYMPLECTIC SEMI-ANALYTICAL METHOD AND ITS APPLICATION IN INHOMOGENEOUS ELECTROMAGNETIC WAVEGUIDES

Dual vectors are applied in Hamilton system of applied mechanics. Electric and magnetic field vectors are the dual vectors in electromagnetic field. The Hamilton system method is introduced into the analysis of electromagnetism waveguide with inhomogeneous materials. The transverse electric and magnetic fields are regarded as the dual. The basic equations are solved in Hamilton system and symplectic geometry. With the Hamilton variational principle, the symplectic semi_analytical equations are derived and preserve their symplectic structures. The given numerical example demonstrates the solution of LSE (Longitudinal Section Electric) mode in a dielectric waveguide. (

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.

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.

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.

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.

Revealing the Hidden Mathematical Beauties of the Cayley-Hamilton Method

The inversion of a non-singular square matrix applying a Computer Algebra System (CAS) is straightforward. The CASs make the numeric computation efficient but mock the mathematical characteristics. The algorithms conducive to the output are sealed and inaccessible. In practice, other than the CPU timing, the applied inversion method is irrelevant. This research-oriented article discusses one such process, the Cayley-Hamilton (C.H.) [1]. Pursuing the process symbolically reveals its unpublished hidden mathematical characteristics even in the original article [1]. This article expands the general vision of the original named method without altering its practical applications. We have used the famous CAS Mathematica [2]. We have briefed the theory behind the method and applied it to different-sized symbolic and numeric matrices. The results are compared to the named CAS’s sealed, packaged library commands. The codes are given, and the algorithms are unsealed.

Hamilton-Jacobi method for solving ordinary differential equations

The Hamilton-Jacobi method for solving ordinary differential equations is presented in this paper. A system of ordinary differential equations of first order or second order can be expressed as a Hamilton system under certain conditions. Then the Hamilton-Jacobi method is used in the integration of the Hamilton system and the solution of the original ordinary differential equations can be found. Finally, an example is given to illustrate the application of the result.

Stress analysis of anisotropic thick laminates in cylindrical bending using a semi-analytical approach

Semi-analytical elasticity solutions for bending of angle-ply laminates in cylindrical bending are presented using the state-space-based differential quadrature method (SSDQM). Partial differential state equation is derived from the basic equations of elasticity based on the state space concept. Then, the differential quadrature (DQ) technique is introduced to discretize the longitu- dinal domain of the plate so that a series of ordinary differential state equations are obtained at the discrete points. Meanwhile, the edge constrained conditions are handled directly using the stress and displacement components without the Saint-Venant principle. The thickness domain is solved analytically based on the state space formalism along with the continuity conditions at interfaces. The present method is validated by comparing the results to the exact solutions of Pagano’s problem. Numerical results for fully clamped thick laminates are presented, and the influences of ply angle on stress distributions are discussed.

Development and Verification of an SP3 Code Using Semi-Analytic Nodal Method for Pin-by-Pin Calculation

SP3 （simplified P3） theory is widely used in LWR （light water reactor） analyses to partly capture the transport effect, especially for pin-by-pin core analysis with pin size homogenization. In this paper, a SP3 code named STELLA is developed and verified at SNERDI （Shanghai Nuclear Engineering Research and Design Institute）. For SP3 method, neutron transport equation can be transformed into two coupled equations in the same mathematical form as diffusion equation. In this work, SANM （semi-analytic nodal method） is used to solve diffusion-like equation, due to its easy to handle multi-group problem. Whole core nodal boundary net current coupling is used to improve convergence stability in SANM, instead of solving two-node problem. CMFD （coarse-mesh finite difference） acceleration method is employed for 0-th SP3 equation, which represents the neutron balance relationship. Three benchmarks are used to verify the SP3 code, STELLA. The first one is a self-defined one dimensional problem, which demonstrates SP3 method is extremely accurate, due to no academic approximation in one dimensional for SP3. The second one is a two dimensional one-group problem cited from Larsen＇s paper, which is usually used to verify and prove the SP3 code correct and accurate. And the third one is modified from 2D C5G7-MOX benchmark, whose numerical results indicate that STELLA is accurate and efficient in pin size level, compared to diffusion model.

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.

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.

Modified domain decomposition method for Hamilton-Jacobi-Bellman equations

This paper presents a modified domain decomposition method for the numerical solution of discrete Hamilton-Jacobi-Bellman equations arising from a class of optimal control problems using diffusion models. A convergence theorem is established. Numerical results indicate the effectiveness and accuracy of the method.

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.

共反射点轨迹的Hamilton方法

本文采用Ｈａｍｉｌｔｏｎ理论方法，研究反射地震记录域中地下同一反射点对应的反射走时随炮检距变化的特性，这在叠前地震资料处理中具有重要意义．由于走时函数可以表示成炮点位置、检波点位置、地震波射线出射角和接收角的函数，通过对共反射点走时与炮点位置变化关系的分析，提出了复杂介质中共反射点轨迹可用Ｈａｍｉｌｔｏｎ正则方程描述．在线性变速介质中，走时函数可以解析给出，由此可导出该介质中共反射点轨迹的Ｈａｍｉｌｔｏｎ正则方程．文中结合常速度、横向变速和线性变速介质模型，计算了地下不同反射面上各点相应的共反射点轨迹，阐明Ｈａｍｉｌｔｏｎ方法的一些特点，并首次引入Ｈａｍｉｌｔｏｎ理论方法描述共反射点走时随炮检距的变化特性，开创了该领域研究的一条新途径．

基于Hamilton-Jacobi方程的飞行器机动动作可达集分析

为了给驾驶员完成标准机动动作提供决策支持,提出一种使用哈密尔顿–雅克比(Hamilton-Jacobi)方程求解机动动作可行状态空间的研究方法.使用关键点将机动动作划分为不同阶段,将各关键点的标准状态约束作为目标集,逆时间求解目标集对应的可达集得到各阶段的边界状态范围,目标集和可达集均由零水平集表示.使用该方法得到斤斗动作三维度运动模型下各阶段的可达集及斤斗动作的可行状态空间,为了使运动模型的控制量与驾驶员实际操纵更为接近,构建了以迎角变化率为控制量的四维度运动模型,在此基础上对斤斗动作各阶段的可达集进行了分析.

基于Hamilton体系的分离变量法

利用Hamilton体系下的分离变量法 ,对可用传统分离变量法分离变量的一类偏微分方程进行了求解 ,得到了完全相同的结果 ,从而说明了Hamilton体系下分离变量法的正确性和潜在能力 ,为Hamilton体系下分离变量法的广泛应用打下了基础 .