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.

Failure mechanism of directional roof cutting and design method optimization

Directional roof cutting(DRC)is one of the key techniques in non-pillar coal mining with self-formed entries(NCMSE)mining method.Due to the inability to accurately measure the expansion coefficient of the goaf rock mass,the implementation of this technology often encounters design challenges,leading to suboptimal results and increased costs.This paper establishes a structural analysis model of the goaf working face roof,revealing the failure mechanism of DRC,and clarifies the positive role of DRC in improving the stress of the roadway surrounding rock and reducing the subsidence of the roof through numerical simulation experiments.On this basis,the paper further analyses the roadway pressure and roof settlement under different DRC design heights,and ultimately proposes an optimized design method for the DRC height.The results indicate that the implementation of DRC can significantly optimize the stress environment of the working face roadway surrounding rock.At the same time,during the application of DRC,three scenarios may arise:insufficient,reasonable,and excessive DRC height.Insufficient height will significantly reduce the effectiveness of the technology,while excessive height has little impact on the implementation effect but will greatly increase construction costs and difficulty.Engineering verification shows that the optimized DRC design method proposed in this paper reduces the peak stress of the protective coal pillar in the roadway by 27.2%and the central subsidence of the roof by 41.8%,demonstrating excellent application results.This method provides technical support for the further promotion of NCMSE mining method.

Study on Measurement Methods of Multi-Directional Waves—Wave Gauge Array Method

Measuring multi-directional waves with the wave gauge array is one of the fundamental and easily realised methods. In this paper, the wave gauge array is described and the effects of the gauge spacing, the array orientations, etc. of the three array arrangements, i. e., linear array, T-type array and pentagon array, on the resolution of the directional spreading of waves, are investigated experimentally. This study can be used as a reference in the experimental study and the field measurement of directional waves.

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.

Study of the pressure transient behavior of directional wells considering the effect of non-uniform flux distribution

During the production,the fluid in the vicinity of the directional well enters the wellbore with different rates,leading to non-uniform flux distribution along the directional well.However,in all existing studies,it is oversimplified to a uniform flux distribution,which can result in inaccurate results for field applications.Therefore,this paper proposes a semi-analytical model of a directional well based on the assumption of non-uniform flux distribution.Specifically,the direction well is discretized into a carefully chosen series of linear sources,such that the complex well trajectory can be captured and the nonuniform flux distribution along the wellbore can be considered to model the three-dimensional flow behavior.By using the finite difference method,we can obtain the numerical solutions of the transient flow within the wellbore.With the aid of Green's function method,we can obtain the analytical solutions of the transient flow from the matrix to the wellbore.The complete flow behavior of a directional well is perfectly represented by coupling the above two types of transient flow.Subsequently,on the basis of the proposed model,we conduct a comprehensive analysis of the pressure transient behavior of a directional well.The computation results show that the flux variation along the direction well has a significant effect on pressure responses.In addition,the directional well in an infinite reservoir may exhibit the following flow regimes:wellbore afterflow,transition flow,inclined radial flow,elliptical flow,horizontal linear flow,and horizontal radial flow.The horizontal linear flow can be observed only if the formation thickness is much smaller than the well length.Furthermore,a dip region that appears on the pressure derivative curve indicates the three-dimensional flow behavior near the wellbore.

A laboratory study of directional spectra with maximum likelihood method─I Developing wind wave

From measurements by a circular array consisting of 18 wave gauges in a large wave tank, directional spectra of wind-generated waves in deep water are systematically determined by using maximum likehood method.The investigations reveal that the angular spreading of the wave energy is consistent with cos2s(θ/2) proposed by Longuet-Higgins et al. (1963, Ocean Wad Spectra,11~136), if the bimodal distributions of wave energy are not taken into account. Bimodality occurring on higher frequency than peak frequency is too rare to affect our whole results. Surprisingly, a much broader directional spreading than that of the field, which is interpreted by the strongly nonlinear energy transfer because of the very young waves in laboratory, is found. The parameter s depends on frequency in the same way as observed by Mitsuyasu et al. (1975, Journal of Physical Oceanography, 5, 750~760)and Hasselmann et al. (1980, Journal of physical Oceanography, 10, 1264~1280) in the field, and the relationship between the four nondimensional parameters sm, fo, b1 and b2, determining the directional width, and (corresponding to the inverse of wave age) are given respectively. The observed distributions are found to agree well with the suggestion of Donelan et al. (1985, Philosophical Transaction of Royal Society of London, A315, 509~562) when applied to field waves.

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.

Superconvergence of Direct Discontinuous Galerkin Methods:Eigen-structure Analysis Based on Fourier Approach

This paper investigates superconvergence properties of the direct discontinuous Galerkin(DDG)method with interface corrections and the symmetric DDG method for diffusion equations.We apply the Fourier analysis technique to symbolically compute eigenvalues and eigenvectors of the amplification matrices for both DDG methods with different coefficient settings in the numerical fluxes.Based on the eigen-structure analysis,we carry out error estimates of the DDG solutions,which can be decomposed into three parts:(i)dissipation errors of the physically relevant eigenvalue,which grow linearly with the time and are of order 2k for P^(k)(k=2,3)approximations;(ii)projection error from a special projection of the exact solution,which is decreasing over the time and is related to the eigenvector corresponding to the physically relevant eigenvalue;(iii)dissipative errors of non-physically relevant eigenvalues,which decay exponentially with respect to the spatial mesh sizeΔx.We observe that the errors are sensitive to the choice of the numerical flux coefficient for even degree P^(2)approximations,but are not for odd degree P^(3)approximations.Numerical experiments are provided to verify the theoretical results.

A joint method to retrieve directional ocean wave spectra from SAR and wave spectrometer data

This paper proposes a joint method to simultaneously retrieve wave spectra at dif ferent scales from spaceborne Synthetic Aperture Radar(SAR) and wave spectrometer data. The method combines the output from the two dif ferent sensors to overcome retrieval limitations that occur in some sea states. The wave spectrometer sensitivity coeffi cient is estimated using an ef fective signifi cant wave height(SWH), which is an average of SAR-derived and wave spectrometer-derived SWH. This averaging extends the area of the sea surface sampled by the nadir beam of the wave spectrometer to improve the accuracy of the estimated sensitivity coeffi cient in inhomogeneous sea states. Wave spectra are then retrieved from SAR data using wave spectrometer-derived spectra as fi rst guess spectra to complement the short waves lost in SAR data retrieval. In addition, the problem of 180° ambiguity in retrieved spectra is overcome using SAR imaginary cross spectra. Simulated data were used to validate the joint method. The simulations demonstrated that retrieved wave parameters, including SWH, peak wave length(PWL), and peak wave direction(PWD), agree well with reference parameters. Collocated data from ENVISAT advanced SAR(ASAR), the airborne wave spectrometer STORM, the PHAROS buoy, and the European Centre for Medium-Range Weather Forecasting(ECMWF) were then used to verify the proposed method. Wave parameters retrieved from STORM and two ASAR images were compared to buoy and ECMWF wave data. Most of the retrieved parameters were comparable to reference parameters. The results of this study show that the proposed joint retrieval method could be a valuable complement to traditional methods used to retrieve directional ocean wave spectra, particularly in inhomogeneous sea states.

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.

Radiation heat transfer model for complex superalloy turbine blade in directional solidification process based on finite element method

For the sake of a more accurate shell boundary and calculation of radiation heat transfer in the Directional Solidification(DS) process, a radiation heat transfer model based on the Finite Element Method(FEM)is developed in this study. Key technologies, such as distinguishing boundaries automatically, local matrix and lumped heat capacity matrix, are also stated. In order to analyze the effect of withdrawing rate on DS process,the solidification processes of a complex superalloy turbine blade in the High Rate Solidification(HRS) process with different withdrawing rates are simulated; and by comparing the simulation results, it is found that the most suitable withdrawing rate is determined to be 5.0 mm·min^(-1). Finally, the accuracy and reliability of the radiation heat transfer model are verified, because of the accordance of simulation results with practical process.

Tip-splitting instability in directional solidification based on bias field method

Tip splitting instability of cellular interface morphology in directional solidification is analyzed based on the bias field method proposed recently by Glicksman. The physical mechanism of tip instability is explained by analyzing the interface potential, the tangential energy flux, and the normal energy flux. A rigorous criterion for tip-splitting instability is established analytically, i.e., the ratio of the cellular tip radius to the cellular width α 〉3/2/π≈ 0.3899, which is in good agreement with simulation results. This study also reveals that the cellular tip splitting instability is attributable to weak Gibbs–Thomson energy acting on the interface.

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.

Fracture Mechanism of 3D, Five-directional Braided (SiO_2)_f /SiO_2 Composites Prepared by Silicasol-infiltration-sintering Method

Three-dimensional （3D） five-directional braided （SiO2）/SiO2 composites were prepared by silicasol-infiltration-sintering （SIS） method. The flexural properties and microstructures were studied. The flexural strength and flexural elastic modulus were found to be 73 MPa and 12 GPa, respectively. The results of stress vs deflection curve and SEM examinations revealed that the fracture mechanism of 3D, five-directional braided （SiO2）/SiO2 composite was a mixture mode of ductile and brittle. The ductile mode was attributed to the weak bonding strength of fiber/matrix at low temperature. The brittle fracture might be caused by the propagation of micro defect or crack, which existed in the as-prepared composites for the ten-cycle process.

Finite element numerical simulation of 2.5D direct current method based on mesh refinement and recoarsement

To deal with the problem of low computational precision at the nodes near the source and satisfy the requirements for computational efficiency in inversion imaging and finite-element numerical simulations of the direct current method, we propose a new mesh refinement and recoarsement method for a two-dimensional point source. We introduce the mesh refinement and mesh recoarsement into the traditional structured mesh subdivision. By refining the horizontal grids, the singularity owing to the point source is minimized and the topography is simulated. By recoarsening the horizontal grids, the number of grid cells is reduced significantly and computational efficiency is improved. Model tests show that the proposed method solves the singularity problem and reduces the number of grid cells by 80% compared to the uniform grid refinement.

THREE-DIMENSIONAL TRAJECTORY OPTIMIZATION WITH DIRECT METHOD

The principle of direct method used in optimal control problem is introduced. Details of applying this method to flight trajectory generation are presented including calculation of velocity and controls histories. And capabilities of flight and propulsion systems are considered also. Combined with digital terrain map technique, the direct method is applied to the three dimensional trajectory optimization for low altitude penetration, and simplex algorithm is used to solve the parameters in optimization. For the small number of parameters, the trajectory can be optimized in real time on board.

Fast alternating direction method of multipliers for total-variation-based image restoration

A novel algorithm, i.e. the fast alternating direction method of multipliers （ADMM）, is applied to solve the classical total-variation （ TV ）-based model for image reconstruction. First, the TV-based model is reformulated as a linear equality constrained problem where the objective function is separable. Then, by introducing the augmented Lagrangian function, the two variables are alternatively minimized by the Gauss-Seidel idea. Finally, the dual variable is updated. Because the approach makes full use of the special structure of the problem and decomposes the original problem into several low-dimensional sub-problems, the per iteration computational complexity of the approach is dominated by two fast Fourier transforms. Elementary experimental results indicate that the proposed approach is more stable and efficient compared with some state-of-the-art algorithms.

Direct computing methods for turn flows in traffic assignment

Two methods based on a slight modification of the regular traffic assignmentalgorithms are proposed to directly compute turn flows instead of estimating them from link flows orobtaining them by expanding the networks. The first one is designed on the path-turn incidencerelationship, and it is similar to the computational procedure of link flows. It applies to thetraffic assignment algorithms that can provide detailed path structures. The second utilizes thelink-turn incidence relationship and the conservation of flow on links, a law deriving from thisrelationship. It is actually an improved version of Dial's logit assignment algorithm. The proposedapproaches can avoid the shortcomings both of the estimation methods, e. g. Furness's model andFrator's model, and of the network-expanding method in precision, stability and computation scale.Finally, they are validated by numerical examples.