Meshless Integral Method for Analysis of Elastoplastic Geotechnical Materials

Author(s):  
Jianfeng Ma ◽  
Joshua David Summers ◽  
Paul F. Joseph

The meshless integral method based on regularized boundary equation [1][2] is extended to analyze elastoplastic geotechnical materials. In this formulation, the problem domain is clouded with a node set using automatic node generation. The sub-domain and the support domain related to each node are also generated automatically using algorithms developed for this purpose. The governing integral equation is obtained from the weak form of elastoplasticity over a local sub-domain and the moving least-squares approximation is employed for meshless function approximation. The geotechnical materials are described by pressure-sensitive multi-surface Drucker-Prager/Cap plasticity constitutive law with hardening. A generalized collocation method is used to impose the essential boundary conditions and natural boundary conditions are incorporated in the system governing equations. A comparison of the meshless results with the FEM results shows that the meshless integral method is accurate and robust enough to solve geotechnical materials.

1972 ◽  
Vol 39 (4) ◽  
pp. 1041-1046 ◽  
Author(s):  
M.-J. Yan ◽  
E. H. Dowell

For constrained-layer damping a simple differential equation for nonsymmetric sandwich plates or beams made of isotropic and homogeneous layers is deduced. The natural boundary conditions associated with this equation are also derived. Typical numerical results are presented including a preliminary comparison with experiment.


1968 ◽  
Vol 35 (4) ◽  
pp. 689-696 ◽  
Author(s):  
J. D. Achenbach ◽  
C. T. Sun ◽  
G. Herrmann

A continuum theory for a laminated medium is further developed in this paper. Constitutive equations, stress equations of motion, and natural boundary conditions are presented, and sufficient conditions for a unique solution are discussed. The governing equations and boundary conditions are employed to study the thickness-twist motion of a laminated layer. For every nodal number there is a low frequency acoustic mode and a high frequency optical mode. The frequencies of the acoustic modes are compared with the corresponding frequencies predicted by the effective modulus theory, and the relative magnitudes of the material parameters for which these frequencies are substantially at variance are indicated.


The matter of appropriate boundary conditions for open or truncated outflow regions in internal flow is still focus of discussion and research. In most practical applications, one can at best estimate mean pressure values or flow rates at such outlets. In the context of finite element methods, it is known that enforcing mean pressures through the pseudo-tractions arising from the Laplacian Navier-Stokes formulation yields accurate, physically consistent solutions. Nevertheless, when generalised Newtonian fluid models are considered, the resulting non-uniform viscosity fields render the classical Laplacian formulation inadequate. Thus, it is common practice to use the socalled stress-divergence formulation with natural boundary conditions known for causing nonphysical outflow behaviour. In order to overcome such a limitation, this work presents a novel mixed variational formulation that can be seen as a generalisation of the Laplacian Navier-Stokes form to fluids with shear-rate-dependent viscosity, as appearing in hemodynamic and polymeric flows. By appropriately manipulating the viscous terms in the variational formulation and employing a simple projection of the constitutive law, it is possible to devise a formulation with the desired natural boundary conditions and low computational complexity. Several numerical examples are presented to showcase the potential of our method, revealing improved accuracy and robustness in comparison with the state of the art.


2020 ◽  
Vol 8 (4) ◽  
pp. 994-1010
Author(s):  
Weizhi Xu

This paper investigates one kind of interpolation for scattered data by bi-cubic polynomial natural spline, in which the integral of square of partial derivative of two orders to x and to y for the interpolating function is minimal (with natural boundary conditions). Firstly, bi-cubic polynomial natural spline interpolations with four kinds of boundary conditions are studied. By the spline function methods of Hilbert space, their solutions are constructed as the sum of bi-linear polynomials and piecewise bi-cubic polynomials. Some properties of the solutions are also studied. In fact, bi-cubic natural spline interpolation on a rectangular domain is a generalization of the cubic natural spline interpolation on an interval. Secondly, based on bi-cubic polynomial natural spline interpolations of four kinds of boundary conditions, and using partition of unity technique, a Partition of Unity Interpolation Element Method (PUIEM) for fitting scattered data is proposed. Numerical experiments show that the PUIEM is adaptive and outperforms state-of-the-art competitions, such as the thin plate spline interpolation and the bi-cubic polynomial natural spline interpolations for scattered data.


Sign in / Sign up

Export Citation Format

Share Document