AHermitian C^(2) Differential Reproducing Kernel Interpolation Meshless Method for the 3D Microstructure-Dependent Static Flexural Analysis of Simply Supported and Functionally Graded Microplates

This work develops a Hermitian C^(2) differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredependent static flexural behavior of a functionally graded(FG)microplate subjected to mechanical loads and placed under full simple supports.In the formulation,we select the transverse stress and displacement components and their first-and second-order derivatives as primary variables.Then,we set up the differential reproducing conditions(DRCs)to obtain the shape functions of the Hermitian C^(2) differential reproducing kernel(DRK)interpolant’s derivatives without using direct differentiation.The interpolant’s shape function is combined with a primitive function that possesses Kronecker delta properties and an enrichment function that constituents DRCs.As a result,the primary variables and their first-and second-order derivatives satisfy the nodal interpolation properties.Subsequently,incorporating ourHermitianC^(2)DRKinterpolant intothe strong formof the3DCCST,we develop a DRKIM method to analyze the FG microplate’s 3D microstructure-dependent static flexural behavior.The Hermitian C^(2) DRKIM method is confirmed to be accurate and fast in its convergence rate by comparing the solutions it produces with the relevant 3D solutions available in the literature.Finally,the impact of essential factors on the transverse stresses,in-plane stresses,displacements,and couple stresses that are induced in the loaded microplate is examined.These factors include the length-to-thickness ratio,the material length-scale parameter,and the inhomogeneity index,which appear to be significant.

Approximate solution for the Klein Gordon-Schrdinger equation by the homotopy analysis method

The Homotopy analysis method is applied to obtain the approximate solution of the Klein-Gordon Schrodinger equation. The Homotopy analysis solutions of the Klein-Gordon Schrodinger equation contain an auxiliary parameter which provides a convenient way to control the convergence region and rate of the series solutions. Through errors analysis and numerical simulation, we can see the approximate solution is very close to the exact solution.

Approximate Analytical Solutions of Fractional Coupled mKdV Equation by Homotopy Analysis Method

In this paper, the approximate analytical solutions of the fractional coupled mKdV equation are obtained by homotopy analysis method (HAM). The method includes an auxiliary parameter which provides a convenient way of adjusting and controlling the convergence region of the series solution. The suitable value of auxiliary parameter is determined and the obtained results are presented graphically.

STRESS FIELD ANALYSIS OF EXTRA-HEIGHT FORGING DIE USING FINITE ELEMENT METHOD

ＴｈｅｅｘｔｒａｈｅｉｇｈｔｆｏｒｇｉｎｇｄｉｅｉｎＳｏｕｔｈｅｗｅｓｔＡｌｕｍｉｎｉｕｍＦａｂｒｉｃａｔｉｏｎＰｌａｎｔｉｓａｓｐｅｃｉａｌｄｉｅｆｏｒｍａｋｉｎｇａｒｍａｍｅｎｔａｎｄｓｐｅｃｉａｌｃｉｖｉｌｐｒｏｄｕｃｔｓｉｎＣｈｉｎａ．Ａ...

Size-dependent effect on biaxial and shear nonlinear buckling analysis of nonlocal isotropic and orthotropic micro-plate based on surface stress and modified couple stress theories using differential quadrature method

The size-dependent effect on the biaxial and shear nonlinear buckling analysis of an isotropic and orthotropic micro-plate based on the surface stress, the modified couple stress theory （MCST）, and the nonlocal elasticity theories using the differential quadrature method （DQM） is presented. Main advantages of the MCST over the classical theory （CT） are the inclusion of the asymmetric couple stress tensor and the consideration of only one material length scale parameter. Based on the nonlinear von Karman assumption, the governing equations of equilibrium for the micro-classical plate consid- ering midplane displacements are derived based on the minimum principle of potential energy. Using the DQM, the biaxial and shear critical buckling loads of the micro-plate for various boundary conditions are obtained. Accuracy of the obtained results is validated by comparing the solutions with those reported in the literature. A parametric study is conducted to show the effects of the aspect ratio, the side-to-thickness ratio, Eringen＇s nonlocal parameter, the material length scale parameter, Young＇s modulus of the surface layer, the surface residual stress, the polymer matrix coefficients, and various boundary conditions on the dimensionless uniaxial, biaxial, and shear critical buckling loads. The results indicate that the critical buckling loads are strongly sensitive to Eringen＇s nonlocal parameter, the material length scale parameter, and the surface residual stress effects, while the effect of Young＇s modulus of the surface layer on the critical buckling load is negligible. Also, considering the size dependent effect causes the increase in the stiffness of the orthotropic micro-plate. The results show that the critical biaxial buckling load increases with an increase in G12/E2 and vice versa for E1/E2. It is shown that the nonlinear biaxial buckling ratio decreases as the aspect ratio increases and vice versa for the buckling amplitude. Because of the most lightweight micro-composite materials with high strength/weight and stiffness/weight ratios, it is anticipated that the results of the present work are useful in experimental characterization of the mechanical properties of micro-composite plates in the aircraft industry and other engineering applications.

Effective stress analysis method of seismic response for high tailings dam

Based on the analysis method for tailings dam in upstream raising method presently used in metallurgy and nonferrous metals tailings depository in the world, an effective stress analysis method of seismic response for high tailings dam was developed according to the results of engineering geological exploration, static and dynamic test and stability analysis on Baizhishan tailing dam 113.5 m high. The law of generation, diffusion and dissipation of seismic pore water pressure during and after earthquake was investigated, and the results of tailings dam’s acceleration, seismic dynamic stress and pore water pressure were obtained. The results show that the seismic stability and liquefaction resistance of high tailings dam are strengthened remarkably, and the scope and depth of liquefaction area at the top of dam are reduced greatly. The interior stress is compressive stress, the stress level of every element is less than 1.0 and the safety coefficient of every element is greater than 1.0. The safety coefficient against liquefaction of every element of tailing dam is greater than 1.5 according to the effective stress analysis of seismic response by finite element method. The calculated results prove that liquefaction is the main reason of seismic failure of high tailing dams, and the effect of seismic inertia forces on high tailing dams’ stability during earthquake is secondary reason.

Thermal stress analysis method considering geometric effect of risers in sand mold casting process

Solidif ication and f luid f low analysis using computer simulation is a current common practice. There is also a high demand for thermal stress analysis in the casting process because casting engineers want to control the defects related to thermal stresses, such as large deformation and crack generation during casting. The riser system is an essential part of preventing the shrinkage defects in the casting process, and it has a great inf luence on thermal phenomena. The analysis domain is dramatically expanded by attaching the riser system to a casting product due to its large volume, and it makes FEM mesh generation diff icult. However, it is diff icult to study and solve the above proposed problem caused by riser system using traditional analysis methods which use single numerical method such as FEM or FDM. In this paper, some research information is presented on the effects of the riser system on thermal stress analysis using a FDM/FEM hybrid method in the casting process simulation. The results show the optimal conditions for stress analysis of the riser model in order to save computation time and memory resources.

ANALYSIS OF INTER-LAMINAR STRESSES FOR COMPOSITE LAMINATED PLATE WITH INTERFACIAL DAMAGE

A constitutive model for composite laminated plates with the damage effect of the intra-layers and inter-laminar interface is presented. The model is based on the general six-degrees-of-freedom plate theory, the discontinuity of displacement on the interfaces are depicted by three shape functions, which are formulated according to solutions satisfying three equilibrium equations, By using the variation principle, the three-dimensional non-linear equilibrium differential equations of the laminated plates with two different damage models are derived. Then, considering a simply supported laminated plate with damage, an analytical solution is presented using finite difference method to obtain the inter-laminar stresses.

General displacement arch-cantilever element method for stress analysis of arch dam

Based on the general displacement method and the basic hypothesis of the trial load method, a new advanced trial load method, the general displacement arch-cantilever element method, was proposed to derive the transformation relation of displacements and loads between the surface nodes and middle plane nodes. This method considers the nodes on upstream and downstream surfaces of the arch dam to be exit nodes （master nodes）, and the middle plane nodes to be slave nodes. According to the derived displacement and load transformation matrices, the equilibrium equation treating the displacement of middle plane nodes as a basic unknown variable is transformed into one that treats the displacement of upstream and downstream nodes as a basic unknown variable. Because the surface nodes have only three degrees of freedom （DOF）, this method can be directly coupled with the finite element method （FEM）, which is used for foundation simulation to analyze the stress of the arch dam with consideration of dam-foundation interaction. Moreover, using the FEM, the nodal load of the arch dam can be easily obtained. Case studies of a typical cylindrical arch dam and the Wudongde arch dam demonstrate the robustness and feasibility of the proposed method.

A NEW METHOD FOR STRESS ANALYSIS OF A CYCLICALLY SYMMETRIC STRUCTURE

In this paper, a computational method for finite element stress analysis of a cyclically symmetric structure subjected to arbitrary loads is provided. At first, using discrete Fourier transformation technique, the complete structure is analyzed by considering only one sector with appropriate complex constraints on its boundary with the adjacent sectors. Next, an imaginary structure which is composed of two identically overlapping sectors is constructed, and that the complex constraints mentioned above can be equivalently replaced by a set of real constraints on this imaginary structure is proved. Therefore, the stress analysis of a cyclically symmetric structure can be solved conveniently by most of finite element programs.

A FINITE ELEMENT METHOD FOR STRESS ANALYSIS OR ELASTOPLASTIC BODY WITH POLYGONAL LINE STRAIN——HARDENING

In this paper, the stress-strain curve of material is fitted by polygonal line composed of three lines. According to the theory of proportional loading in elastoplasticity, we simplify the complete stress-strain relations, which are given by the increment theory of elastoplasticity. Thus, the finite element equation with the solution of displacement is derived. The assemblage elastoplastic stiffness matrix can be obtained by adding something to the elastic matrix, hence it will shorten the computing time. The determination of every loading increment follows the von Mises yield criteria. The iterative method is used in computation. It omits the redecomposition of the assemblage stiffness matrix and it will step further to shorten the computing time. Illustrations are given to the high-order element application departure from proportional loading, the computation of unloading fitting to the curve and the problem of load estimation.

THE APPLICATION OF BOUNDARY ELEMENT METHOD IN STRESS ANALYSIS OF AUTOFRETTED TUBE WITH NOTCH

In this paper, the boundary element method is applied to investigate the internal state of stress of autofretted tube with notch and the calculated results are important in the practical design.

THE APPLICATION OF COMPATIBLE STRESS ITERATIVE METHOD IN DYNAMIC FINITE ELEMENT ANALYSIS OF HIGH VELOCITY IMPACT

There is a common difficulty in elastic-plastic impact codes such as EPIC[2,3] NONSAP[4], etc.. Most of these codes use the simple linear functions usually taken from static problem to represent the displacement components. In such finite element formulation, the stress components are constant in each element and they are discontinuous in any two neighboring elements. Therefore, the bases of using the virtual work principle in such elements are unreliable. In this paper, we introduce a new method, namely, the compatible stress iterative method, to eliminate the above-said difficulty. The calculated examples show that the calculation using the new method in dynamic finite element analysis of high velocity impact is valid and stable, and the element stiffness can be somewhat reduced.

Photoelastic Stress Analysis by Use of Hybrid Technique and Fringe Phase Shifting Method

Photoelastic fringe patterns for stress analysis are investigated by use of hybrid technique and fringe phase shift method. The first one is a hybrid method which combines the conformal mapping technique and measured data away from the edge of a geometric discontinuity. Photoelastic data are hybridized with complex variable/mapping techniques to calculate photoelastic stress-field around a circular hole or an elliptical hole in plates under uniaxial tensile loading. This method determines full-field stresses in perforated finite tensile plates containing either a circular or an elliptical hole. The second one is a fringe phase shift method to separate isochromatics and isoclinics from photoelastic fringes of a circular disk under diametric compression by use of phase shift method. Digitally determined isochromatics and isoclinics are agreed well with those of manual measurements.

DDM regression analysis of the in-situ stress field in a non-linear fault zone

A multivariable regression analysis of the in-situ stress field, which considers the non-linear deformation behavior of faults in practical projects, is presented based on a newly developed three-dimensional displacement discontinuity method （DDM） program. The Bar- ton-Bandis model and the Kulhaway model are adopted as the normal and the tangential deformation model of faults, respectively, where the Mohr-Coulomb failure criterion is satisfied. In practical projects, the values of the mechanical parameters of rock and faults are restricted in a bounded range for in-situ test, and the optimal mechanical parameters are obtained from this range by a loop. Comparing with the traditional finite element method （FEM）, the DDM regression results are more accurate.

Design and Analysis of Tapered Stress Joint for Top Tensioned Riser System

Stress Joint （SJ） plays a key role in the Top Tensioned Riser （TTR） system for deep water engineering. A preliminary design method of tapered SJ is proposed in the paper, which could help designers obtain accurate design data. After a further sensitive analysis is carried out, the related parameters choice and control methods are recommended in the engineering practice. By taking the extreme environment conditions into consideration, the effects of bending stress reduction and curve control are analyzed, and the 3-D FE models are established by ABQOUS for numerical evaluation to verify the correctness of design results. At last, dynamic analysis and fatigue analysis, based on actual project, are carried out with designed stress joint. The analysis results prove the feasibility and guidance of this method in the practical engineering applications.

A Method Combining Numerical Analysis and Limit Equilibrium Theory to Determine Potential Slip Surfaces in Soil Slopes

This paper describes a precise method combining numerical analysis and limit equilibrium theory to determine potential slip surfaces in soil slopes. In this method, the direction of the critical slip surface at any point in a slope is determined using the Coulomb’s strength principle and the extremum principle based on the ratio of the shear strength to the shear stress at that point. The ratio, which is considered as an analysis index, can be computed once the stress field of the soil slope is obtained. The critical slip direction at any point in the slope must be the tangential direction of a potential slip surface passing through the point. Therefore, starting from a point on the top of the slope surface or on the horizontal segment outside the slope toe, the increment with a small distance into the slope is used to choose another point and the corresponding slip direction at the point is computed. Connecting all the points used in the computation forms a potential slip surface exiting at the starting point. Then the factor of safety for any potential slip surface can be computed using limit equilibrium method like Spencer method. After factors of safety for all the potential slip surfaces are obtained, the minimum one is the factor of safety for the slope and the corresponding potential slip surface is the critical slip surface of the slope. The proposed method does not need to pre-assume the shape of potential slip surfaces. Thus it is suitable for any shape of slip surfaces. Moreover the method is very simple to be applied. Examples are presented in this paper to illustrate the feasibility of the proposed method programmed in ANSYS software by macro commands.

LOWER BOUND LIMIT ANALYSIS OF THREE-DIMENSIONAL ELASTOPLASTIC STRUCTURES BY BOUNDARY ELEMENT METHOD

Based on the lower bound theorem of limit analysis, a solution procedure for limit analysis of three_dimensional elastoplastic structures was established using conventional boundary element method (BEM). The elastic stress field for lower bound limit analysis was computed directly by three_dimensional boundary element method (3_D BEM). The self_equilibrium stress field was constructed by the linear combination of several self_equilibrium “basis vectors” which can be computed by elastic_plastic incremental iteration of 3_D BEM analysis. The lower bound limit analysis problem was finally reduced to a series of nonlinear programming sub_problems with relatively few optimal variables. The complex method was used to solve the nonlinear programming sub_problems. The numerical results show that the present solution procedure has good accuracy and high efficiency.

Solitary Solution of Discrete mKdV Equation by Homotopy Analysis Method

In this paper, we apply homotopy analysis method to solve discrete mKdV equation and successfully obtain the bell-shaped solitary solution to mKdV equation. Comparison between our solution and the exact solution shows that homotopy analysis method is effective and validity in solving hybrid nonlinear problems, including solitary solution of difference-differential equation.

Finite element analysis of slope stability by expanding the mobilized principal stress Mohr's circles-Development, encoding and validation

In recent years,finite element analysis is increasingly being proposed in slope stability problems as a competitive method to traditional limit equilibrium methods(LEMs)which are known for their inherent deficiencies.However,the application of finite element method(FEM)to slope stability as a strength reduction method(SRM)or as finite element limit analysis(FELA)is not always a success for the drawbacks that characterize both methods.To increase the performance of finite element analysis in this problem,a new approach is proposed in this paper.It consists in gradually expanding the mobilized stress Mohr’s circles until the soil failure occurs according to a prescribed non-convergence criterion.The present approach called stress deviator increasing method(SDIM)is considered rigorous for three main reasons.Firstly,it preserves the definition of the factor of safety(FOS)as the ratio of soil shear strength to the mobilized shear stress.Secondly,it maintains the progressive development of shear stress resulting from the increase in the principal stress deviator on the same plane,on which the shear strength takes place.Thirdly,by introducing the concept of equivalent stress loading,the resulting trial stresses are checked against the violation of the actual yield criterion formed with the real strength parameters rather than those reduced by a trial factor.The new numerical procedure was encoded in a Fortran computer code called S^(4)DINA and verified by several examples.Comparisons with other numerical methods such as the SRM,gravity increasing method(GIM)or even FELA by assessing both the FOS and contours of equivalent plastic strains showed promising results.