An advanced implicit meshless approach for fractional partial differential equation in computational mechanics

Recently, the numerical modelling and simulation for fractional partial differential equations (FPDE), which have been found with widely applications in modern engineering and sciences, are attracting increased attentions. The current dominant numerical method for modelling of FPDE is the explicit Finite Difference Method (FDM), which is based on a pre-defined grid leading to inherited issues or shortcomings. This paper aims to develop an implicit meshless approach based on the radial basis functions (RBF) for numerical simulation of time fractional diffusion equations. The discrete system of equations is obtained by using the RBF meshless shape functions and the strong-forms. The stability and convergence of this meshless approach are then discussed and theoretically proven. Several numerical examples with different problem domains are used to validate and investigate accuracy and efficiency of the newly developed meshless formulation. The results obtained by the meshless formations are also compared with those obtained by FDM in terms of their accuracy and efficiency. It is concluded that the present meshless formulation is very effective for the modelling and simulation for FPDE.

A sub‒domain smoothed Galerkin method for solid mechanics problems

A sub‒domain smoothed Galerkin method is proposed to integrate the advantages of mesh‒free Galerkin method and FEM. Arbitrarily shaped sub‒domains are predefined in problems domain with mesh‒free nodes. In each sub‒domain, based on mesh‒free Galerkin weak formulation, the local discrete equation can be obtained by using the moving Kriging interpolation, which is similar to the discretization of the high‒order finite elements. Strain smoothing technique is subsequently applied to the nodal integration of sub‒domain by dividing the sub‒domain into several smoothing cells. Moreover, condensation of DOF can also be introduced into the local discrete equations to improve the computational efficiency. The global governing equations of present method are obtained on the basis of the scheme of FEM by assembling all local discrete equations of the sub‒domains. The mesh‒free properties of Galerkin method are retained in each sub‒domain. Several 2D elastic problems have been solved on the basis of this newly proposed method to validate its computational performance. These numerical examples proved that the newly proposed sub‒domain smoothed Galerkin method is a robust technique to solve solid mechanics problems based on its characteristics of high computational efficiency, good accuracy, and convergence.

Mechanics of microfilaments networks : a cross-scales study

The mechanical properties of microfilament networks are systematically summarized at different special scales in this paper. We have presented the mechanical models of single microfilaments and microfilament networks at microscale. By adopting a coarse-grained simulation strategy, the mechanical stability of microfilaments related cellular structures are analysed. Structural analysis is conducted to microfilament networks to understand the stress relaxation under compression. The nanoscale molecular mechanisms of the microfilaments deformation is also summarized from the viewpoint of molecular dynamics simulation. This paper provides the fundaments of multiscale modelling framework for the mechanical behaviours simulation of hierarchical microfilament networks.

Uncertainty quantification applied to the performance analysis of a conical diffuser

This paper offers an uncertainty quantification (UQ) study applied to the performance analysis of the ERCOFTAC conical diffuser. A deterministic CFD solver is coupled with a non-statistical generalised Polynomial Chaos(gPC)representation based on a pseudo-spectral projection method. Such approach has the advantage to not require any modification of the CFD code for the propagation of random disturbances in the aerodynamic field. The stochactic results highlihgt the importance of the inlet velocity uncertainties on the pressure recovery both alone and when coupled with a second uncertain variable. From a theoretical point of view, we investigate the possibility to build our gPC representation on arbitray grid, thus increasing the flexibility of the stochastic framework.

Prediction of masonry compressive behaviour using a damage mechanics inspired modelling method

Masonry under compression is affected by the properties of its constituents and their interfaces. In spite of extensive investigations of the behaviour of masonry under compression, the information in the literature cannot be regarded as comprehensive due to ongoing inventions of new generation products – for example, polymer modified thin layer mortared masonry and drystack masonry. As comprehensive experimental studies are very expensive, an analytical model inspired by damage mechanics is developed and applied to the prediction of the compressive behaviour of masonry in this paper. The model incorporates a parabolic progressively softening stress-strain curve for the units and a progressively stiffening stress-strain curve until a threshold strain for the combined mortar and the unit-mortar interfaces is reached. The model simulates the mutual constraints imposed by each of these constituents through their respective tensile and compressive behaviour and volumetric changes. The advantage of the model is that it requires only the properties of the constituents and considers masonry as a continuum and computes the average properties of the composite masonry prisms/wallettes; it does not require discretisation of prism or wallette similar to the finite element methods. The capability of the model in capturing the phenomenological behaviour of masonry with appropriate elastic response, stiffness degradation and post peak softening is presented through numerical examples. The fitting of the experimental data to the model parameters is demonstrated through calibration of some selected test data on units and mortar from the literature; the calibrated model is shown to predict the responses of the experimentally determined masonry built using the corresponding units and mortar quite well. Through a series of sensitivity studies, the model is also shown to predict the masonry strength appropriately for changes to the properties of the units and mortar, the mortar joint thickness and the ratio of the height of unit to mortar joint thickness. The unit strength is shown to affect the masonry strength significantly. Although the mortar strength has only a marginal effect, reduction in mortar joint thickness is shown to have a profound effect on the masonry strength. The results obtained from the model are compared with the various provisions in the Australian Masonry Structures Standard AS3700 (2011) and Eurocode 6.

Effect of interface friction on the mechanics of indentation of a finite layer resting on a rigid substrate

Copper strips of 2.5 mm thickness resting on stainless steel anvils were normally indented by wedges under nominal plane strain conditions. Inflections in the hardness-penetration characteristics were identified. Inflections separate stages where each stage has typical mechanics of deformation. These are arrived at by studying the distortion of 0.125 mm spaced grids inscribed on the deformation plane of the strip. The sensitivity of hardness and deformation mechanics to wedge angle and the interfacial friction between strip and anvil were investigated within the framework of existing slip line field models of indentation of semi-infinite and finite blocks.

Spectral Finite Element Formulation for Nanorods via Nonlocal Continuum Mechanics

In this article, the Eringen's nonlocal elasticity theory has been incorporated into classical/local Bernoulli-Euler rod model to capture unique properties of the nanorods under the umbrella of continuum mechanics theory. The spectral finite element (SFE) formulation of nanorods is performed. SFE formulation is carried out and the exact shape functions (frequency dependent) and dynamic stiffness matrix are obtained as function of nonlocal scale parameter. It has been found that the small scale affects the exact shape functions and the elements of the dynamic stiffness matrix. The results presented in this paper can provide useful guidance for the study and design of the next generation of nanodevices that make use of the wave dispersion properties of carbon nanotubes.

Down-to-earth temperatures: The mechanics of the thermal environment

This lecture describes some recent attempts at unravelling the mechanics of the temperature distribution near ground, especially during calm, clear nights. In particular, a resolution is offered of the so-called Ramdas paradox, connected with observations of a temperature minimum some decimetres above bare soil on calm clear nights, in apparent defiance of the Rayleigh criterion for instability due to thermal convection. The dynamics of the associated temperature distribution is governed by radiative and convective transport and by thermal conduction, and is characterised by two time constants, involving respectively quick radiative adjustments and slow diffusive relaxation. The theory underlying the work described here suggests that surface parameters like ground emissivity and soil thermal conductivity can exert appreciable influence on the development of nocturnal inversions.

An adaptive singular ES-FEM for mechanics problems with singular field of arbitrary order

This paper presents a singular edge-based smoothed finite element method (sES-FEM) for mechanics problems with singular stress fields of arbitrary order. The sES-FEM uses a basic mesh of three-noded linear triangular (T3) elements and a special layer of five-noded singular triangular elements (sT5) connected to the singular-point of the stress field. The sT5 element has an additional node on each of the two edges connected to the singular-point. It allows us to represent simple and efficient enrichment with desired terms for the displacement field near the singular-point with the satisfaction of partition-of-unity property. The stiffness matrix of the discretized system is then obtained using the assumed displacement values (not the derivatives) over smoothing domains associated with the edges of elements. An adaptive procedure for the sES-FEM is proposed to enhance the quality of the solution with minimized number of nodes. Several numerical examples are provided to validate the reliability of the present sES-FEM method. (C) 2012 Elsevier B.V. All rights reserved.

Estimation of mechanical properties of single wall carbon nanotubes using molecular mechanics approach

Molecular mechanics based finite element analysis is adopted in the current work to evaluate the mechanical properties of Zigzag, Armchair and Chiral Single wall Carbon Nanotubes (SWCNT) of different diameters and chiralities. Three different types of atomic bonds, that is Carbon Carbon covalent bond and two types of Carbon Carbon van der Waals bonds are considered in the carbon nanotube system. The stiffness values of these bonds are calculated using the molecular potentials, namely Morse potential function and Lennard-Jones interaction potential function respectively and these stiffness's are assigned to spring elements in the finite element model of the CNT. The geometry of CNT is built using a macro that is developed for the finite element analysis software. The finite element model of the CNT is constructed, appropriate boundary conditions are applied and the behavior of mechanical properties of CNT is studied.

A thermodynamic correlation between damage and fracture as applied to concrete fatigue

Fatigue damage in concrete is characterized by the simultaneous presence of micro and macrocracks. The theory of fracture mechanics conveniently handles the propagation of macrocracks, whereas damage mechanics precisely describes the state of microcracking. This paper provides a platform to correlate fracture mechanics and damage mechanics theories through an energy equivalence within a thermodynamic framework by equating the energy dissipated according to each theory. Through this correlation, damage corresponding to a given crack length could be obtained, and alternatively a discrete crack could be transformed into an equivalent damage zone. The results are validated using available experimental data on concrete fatigue including stiffness degradation and acoustic emission. (C) 2015 Elsevier Ltd. All rights reserved.

A thermodynamic correlation between damage and fracture as applied to concrete fatigue

Fatigue damage in concrete is characterized by the simultaneous presence of micro and macrocracics. The theory of fracture mechanics conveniently handles the propagation of macrocracks, whereas damage mechanics precisely describes the state of microcracking. This paper provides a platform to correlate fracture mechanics and damage mechanics theories through an energy equivalence within a thermodynamic framework by equating the energy dissipated according to each theory. Through this correlation, damage corresponding to a given crack length could be obtained, and alternatively a discrete crack could be transformed into an equivalent damage zone. The results are validated using available experimental data on concrete fatigue including stiffness degradation and acoustic emission. (C) 2015 Elsevier Ltd. All rights reserved.

Fracture mechanics of piezoelectric materials

This paper presents an analysis of crack problems in homogeneous piezoelectrics or on the interfaces between two dissimilar piezoelectric materials based on the continuity of normal electric displacement and electric potential across the crack faces. The explicit analytic solutions are obtained for a single crack in an infinite piezoelectric or on the interface of piezoelectric bimaterials. For homogeneous materials it is found that the normal electric displacement D-2, induced by the crack, is constant along the crack faces which depends only on the remote applied stress fields. Within the crack slit, the perturbed electric fields induced by the crack are also constant and not affected by the applied electric displacement fields. For bimaterials, generally speaking, an interface crack exhibits oscillatory behavior and the normal electric displacement D-2 is a complex function along the crack faces. However, for bimaterials, having certain symmetry, in which an interface crack displays no oscillatory behavior, it is observed that the normal electric displacement D-2 is also constant along the crack faces and the electric field E-2 has the singularity ahead of the crack tip and has a jump across the interface. Energy release rates are established for homogeneous materials and bimaterials having certain symmetry. Both the crack front parallel to the poling axis and perpendicular to the poling axis are discussed. It is revealed that the energy release rates are always positive for stable materials and the applied electric displacements have no contribution to the energy release rates.

Environmental Mechanics Research in China

By recalling mankind's path during past 50 years in the present article, we mainly highlight the significance of environmental issues today. In particular, two major factors leading to environment deterioration in China such as water resources and coal burning are stressed on. Present-day environmental issues are obviously interdisciplinary, of multiple scales and multi-composition in nature. Therefore, a process-based approach for environment research is absolutely necessarily. A series of sub-processes, either physical, chemical or biological, are subsequently analyzed in order to established reasonable parameterization scheme and credible comprehensive model. And we are now in a position to answer questions still open to us, improve existing somewhat empirical engineering approaches and enhance quantitative accuracy in prediction. To illustrate this process-based research approach, three typical examples associated with the Yangtze River Estuary, Loess Plateau and Tenggeli Desert environments have been dealt with respectively. A theoretical model of vertical flow field accounting for runoff and tide interaction has been established to delineate salinity and sediment motion which are responsible for the formation of mouth bar at the outlet and the ecological evolution there. A kinematic wave theory combined with the revised Green-Ampt infiltration formula is applied to the prediction of runoff generation and erosion in three types of erosion region on the Loess Plateau. Three approaches describing water motion in SPAC system in arid areas at different levels have been improved by introducing vegetation sub-models. However, we have found that the formation of a dry sandy layer and biological crust skin are additional primary causes leading to deterioration of water supply and succession of ecological system.