Bending of Euler–Bernoulli nanobeams based on the strain-driven and stress-driven nonlocal integral models: a numerical approach

Eringen’s nonlocal elasticity theory is extensively employed for the analysis of nanostructures because it is able to capture nanoscale effects. Previous studies have revealed that using the differential form of the strain-driven version of this theory leads to paradoxical results in some cases, such as bending analysis of cantilevers, and recourse must be made to the integral version. In this article, a novel numerical approach is developed for the bending analysis of Euler–Bernoulli nanobeams in the context of strain- and stress-driven integral nonlocal models. This numerical approach is proposed for the direct solution to bypass the difficulties related to converting the integral governing equation into a differential equation. First, the governing equation is derived based on both strain-driven and stress-driven nonlocal models by means of the minimum total potential energy. Also, in each case, the governing equation is obtained in both strong and weak forms. To solve numerically the derived equations, matrix differential and integral operators are constructed based upon the finite difference technique and trapezoidal integration rule. It is shown that the proposed numerical approach can be efficiently applied to the strain-driven nonlocal model with the aim of resolving the mentioned paradoxes. Also, it is able to solve the problem based on the strain-driven model without inconsistencies of the application of this model that are reported in the literature.     

Modelling cell wall growth using a fibre-reinforced hyperelastic–viscoplastic constitutive law

A fibre-reinforced hyperelastic–viscoplastic model using a finite strain Finite Element (FE) analysis is presented to study the expansive growth of cell walls. Based on the connections between biological concepts and plasticity theory, e.g. wall-loosening and plastic yield, wall-stiffening and plastic hardening, the modelling of cell wall growth is established within a framework of anisotropic viscoplasticity aiming to represent the corresponding biology-controlled behaviour of a cell wall. In order to model in vivo growth, special attention is paid to the differences between a living cell and an isolated wall. The proposed hyperelastic–viscoplastic theory provides a unique framework to clarify the interplay between cellulose microfibrils and cell wall matrix and how this interplay regulates sustainable growth in a particular direction while maintaining the mechanical strength of the cell walls by new material deposition. Moreover, the effect of temperature is taken into account. A numerical scheme is suggested and FE case studies are presented and compared with experimental data.     

Nonlinear vibration analysis of fractional viscoelastic Euler–Bernoulli nanobeams based on the surface stress theory

The nonlinear vibrations of viscoelastic Euler–Bernoulli nanobeams are studied using the fractional calculus and the Gurtin–Murdoch theory. Employing Hamilton's principle, the governing equation considering surface effects is derived. The fractional integro-partial differential governing equation is first converted into a fractional–ordinary differential equation in the time domain using the Galerkin scheme. Thereafter, the set of nonlinear fractional time-dependent equations expressed in a state-space form is solved using the predictor–corrector method. Finally, the effects of initial displacement, fractional derivative order, viscoelasticity coefficient, surface parameters and thickness-to-length ratio on the nonlinear time response of simply-supported and clamped-free silicon viscoelastic nanobeams are investigated.     

A modified beam model based on Gurtin–Murdoch surface elasticity theory

In this work, a modified surface-effect incorporated beam model based on Gurtin and Murdoch (GM) surface elasticity theory is established by satisfying the required balance equations on surfaces, which is often overlooked by researchers in this field. With the refinement, the proposed model is more rigorous in mathematics and mechanics compared with GM theory-based beam models in literature. To demonstrate the model, the problem for static bending of simply supported beam considering surface effects is solved by applying the general equations derived, and numerical results are obtained and discussed.

     

Kinematics of Bulkley–Herschel fluid flow with a free surface during the filling of a channel

The non-Newtonian fluid flow with a free surface occurring during the filling of a plane channel in the gravity field is modeled. The mathematical formulation of the problem using the rheological Bulkley–Herschel model is presented. A numerical finite-difference algorithm for solving this problem is developed. A parametric investigation of the main characteristics of the process as functions of the control parameters is performed. The effect of the rheological parameters of the fluid on the distribution of the quasisolid motion zones is demonstrated.     

Sliding contact conditions using the master–slave approach with application on geometrically non-linear beams

Frictionless sliding conditions between two bodies are usually defined using either the method of Lagrangian multipliers or by prescribing an artificial (penalty) stiffness which resists the penetration at the contact point. Both of these methods impose the condition that the contact force should be normal to the contact surface, with the Lagrangian multiplier or the penalty parameter serving as a measure of this force. In this work, an alternative approach is undertaken: the frictionless sliding condition is defined through a relationship between nodal parameters of the virtual displacements of a discretised principle of virtual work. This method, which does not involve additional force parameters or degrees of freedom, is known as the master–slave or the minimum-set method and is particularly convenient for displacement-based finite-element implementation. The method is analysed in detail in context of bilateral sliding constraints characteristic of prismatic and cylindrical joints in flexible beam assemblies undergoing large overall motion. Two numerical examples are presented and assessed against the results in the literature.     

Synchronization control for networked mobile robot systems based on Udwadia–Kalaba approach

This paper addresses the problem of synchronization control for networked multi-mobile robot systems from the perspective of analytical mechanics. By reformulating the task requirement as a constrained motion problem, a unified synchronization algorithm for networked multi-mobile robot systems with or without leaders is proposed in combination with algebraic graph theory and the Udwadia–Kalaba approach. With the proposed algorithm, the networked mobile robot system can achieve synchronization from arbitrary initial conditions for the leaderless case and realize accurate trajectory tracking with explicitly given reference trajectories for the leader-following case. Numerical simulations of a networked wheeled mobile robot system are performed under different network structures and various trajectory requirements to show the performance of the proposed control algorithm.

     

Waves on the free surface of a two‐phase medium

A boundaryvalue problem is posed to determine the wave motion caused by propagation of a gravity wave on the free surface of a layer of a twophase medium. The problem is solved analytically in the linear approximation. The shape of the free surface, the phase velocity, and the frequency and damping factor of the wave are determined. An example of the solution of the problem is given.     

The energy momentum tensor in the presence of body forces and the Peach–Koehler force on a dislocation

A modified energy momentum tensor, in the presence of body forces, is introduced and used to construct the nonconserved J, M, and L integrals, and to derive the energetic forces associated with a defect motion within the material. The J integral is then applied to evaluate the Peach–Koehler force on an inclined edge dislocation within a large block due to its own weight. The equilibrium position of the dislocation is determined for different boundary conditions of interest in geomechanics.     

Hybrid LES–RANS technique based on a one-equation near-wall model

In order to reduce the high computational effort of wall-resolved large-eddy simulations (LES), the present paper suggests a hybrid LES–RANS approach which splits up the simulation into a near-wall RANS part and an outer LES part. Generally, RANS is adequate for attached boundary layers requiring reasonable CPU-time and memory, where LES can also be applied but demands extremely large resources. Contrarily, RANS often fails in flows with massive separation or large-scale vortical structures. Here, LES is without a doubt the best choice. The basic concept of hybrid methods is to combine the advantages of both approaches yielding a prediction method, which, on the one hand, assures reliable results for complex turbulent flows, including large-scale flow phenomena and massive separation, but, on the other hand, consumes much fewer resources than LES, especially for high Reynolds number flows encountered in technical applications. In the present study, a non-zonal hybrid technique is considered (according to the signification retained by the authors concerning the terms zonal and non-zonal), which leads to an approach where the suitable simulation technique is chosen more or less automatically. For this purpose the hybrid approach proposed relies on a unique modeling concept. In the LES mode a subgrid-scale model based on a one-equation model for the subgrid-scale turbulent kinetic energy is applied, where the length scale is defined by the filter width. For the viscosity-affected near-wall RANS mode the one-equation model proposed by Rodi et al. (J Fluids Eng 115:196–205, 1993) is used, which is based on the wall-normal velocity fluctuations as the velocity scale and algebraic relations for the length scales. Although the idea of combined LES–RANS methods is not new, a variety of open questions still has to be answered. This includes, in particular, the demand for appropriate coupling techniques between LES and RANS, adaptive control mechanisms, and proper subgrid-scale and RANS models. Here, in addition to the study on the behavior of the suggested hybrid LES–RANS approach, special emphasis is put on the investigation of suitable interface criteria and the adjustment of the RANS model. To investigate these issues, two different test cases are considered. Besides the standard plane channel flow test case, the flow over a periodic arrangement of hills is studied in detail. This test case includes a pressure-induced flow separation and subsequent reattachment. In comparison with a wall-resolved LES prediction encouraging results are achieved.      

Displacement of a viscous fluid from an annular Hele–Shaw cell with a sink within the framework of the Brinkman model

The stability of the radial front of viscous fluid displacement from an annular Hele–Shaw cell with a sink of finite radius is analyzed. It is shown that in the absence of the surface tension and at a minimal manifestation of molecular diffusion the role of the stabilizing factor determining the displacement front structure can be played by small viscous forces in the cell plane. The viscous fingers formed in this case turn out to be wider than those in a rectangular cell.     

Effects of tension–compression asymmetry on the surface wrinkling of film–substrate systems

Many soft materials and biological tissues are featured with the tension–compression asymmetry of constitutive relations. The surface wrinkling of a stiff thin film lying on a compliant substrate is investigated through theoretical analysis and numerical simulations. It is found that the tension–compression asymmetry of the soft substrate not only affects the critical strain of buckling but, more importantly, may also influence the wrinkling pattern that occurs in the film–substrate system under specified loading conditions. Due to this mechanism, the thin film subjected to equi-biaxial compression may first buckle into a hexagonal array of dimples or bulges, instead of the checkerboard pattern, and consequently evolve into labyrinths with further loading. Under non-equi-biaxial compression, the system may buckle either into a parallel bead-chain pattern or a stripe pattern, depending on the substrate nonlinearity and the loading biaxiality. Phase diagrams are established for the wrinkling patterns in a wide range of geometric and mechanical parameters, which facilitate the design of surface patterns with desired properties and functions.     

Direct simulation of liquid–gas–solid flow with a free surface lattice Boltzmann method

Direct numerical simulation of liquid–gas–solid flows is uncommon due to the considerable computational cost. As the grid spacing is determined by the smallest involved length scale, large grid sizes become necessary – in particular, if the bubble–particle aspect ratio is on the order of 10 or larger. Hence, it arises the question of both feasibility and reasonability. In this paper, we present a fully parallel, scalable method for direct numerical simulation of bubble–particle interaction at a size ratio of 1–2 orders of magnitude that makes simulations feasible on currently available super-computing resources. With the presented approach, simulations of bubbles in suspension columns consisting of more than 100,000 fully resolved particles become possible. Furthermore, we demonstrate the significance of particle-resolved simulations by comparison to previous unresolved solutions. The results indicate that fully resolved direct numerical simulation is indeed necessary to predict the flow structure of bubble–particle interaction problems correctly.     

Minimizers of Caffarelli–Kohn–Nirenberg Inequalities with the Singularity on the Boundary

Let Ω be a bounded smooth domain in ${{\bf R}^N, N\geqq 3}Let Ω be a bounded smooth domain in R^N, N\geqq 3{{\bf R}^N, N\geqq 3}, and D_a^1,2(W){D_a^{1,2}(\Omega)} be the completion of C₀^￥(W){C_0^\infty(\Omega)} with respect to the norm:

||u||a2=òW |x|-2a|?u|2dx.||u||_a^2=\int_\Omega |x|^{-2a}|\nabla u|^2{d}x.      16. On the minimal speed of traveling waves for a nonlocal delayed reaction–diffusion equation      M. Aguerrea  G. Valenzuela 《Nonlinear Oscillations》2010,13(1):1-9 In this note, we give constructive upper and lower bounds for the minimal speed of propagation of traveling waves for a nonlocal delayed reaction–diffusion equation.      17. Tension–torsion fracture experiments – Part II: Simulations with the extended Gurson model and a ductile fracture criterion based on plastic strain      Zhenyu Xue  Jonas Faleskog  John W. Hutchinson 《International Journal of Solids and Structures》2013,50(25-26):4258-4269 An extension of the Gurson model that incorporates damage development in shear is used to simulate the tension–torsion test fracture data presented in Faleskog and Barsoum (2013) (Part I) for two steels, Weldox 420 and 960. Two parameters characterize damage in the constitutive model: the effective void volume fraction and a shear damage coefficient. For each of the steels, the initial effective void volume fraction is calibrated against data for fracture of notched round tensile bars and the shear damage coefficient is calibrated against fracture in shear. The calibrated constitutive model reproduces the full range of data in the tension–torsion tests thereby providing a convincing demonstration of the effectiveness of the extended Gurson model. The model reinforces the experiments by highlighting that for ductile alloys the effective plastic strain at fracture cannot be based solely on stress triaxiality. For nominally isotropic alloys, a ductile fracture criterion is proposed for engineering purposes that depends on stress triaxiality and a second stress invariant that discriminates between axisymmetric stressing and shear dominated stressing.      18. A new flux-limiting approach–based kinetic scheme for the Euler equations of gas dynamics      Raushan Kumar  Anoop K. Dass 《国际流体数值方法杂志》2019,90(1):22-56 This paper proposes a new kinetic-theory-based high-resolution scheme for the Euler equations of gas dynamics. The scheme uses the well-known connection that the Euler equations are suitable moments of the collisionless Boltzmann equation of kinetic theory. The collisionless Boltzmann equation is discretized using Sweby's flux-limited method and the moment of this Boltzmann level formulation gives a Euler level scheme. It is demonstrated how conventional limiters and an extremum-preserving limiter can be adapted for use in the scheme to achieve a desired effect. A simple total variation diminishing criteria relaxing parameter results in improving the resolution of the discontinuities in a significant way. A 1D scheme is formulated first and an extension to 2D on Cartesian meshes is carried out next. Accuracy analysis suggests that the scheme achieves between first- and second-order accuracy as is expected for any second-order flux-limited method. The simplicity and the explicit form of the conservative numerical fluxes add to the efficiency of the scheme. Several standard 1D and 2D test problems are solved to demonstrate the robustness and accuracy.      19. Computing flow-induced stresses of injection molding based on the Phan–Thien–Tanner model      Wei Cao  Changyu Shen  Chunjie Zhang  Lixia Wang 《Archive of Applied Mechanics (Ingenieur Archiv)》2008,78(5):363-377 A numerical approach is introduced to solve the viscoelastic flow problem of filling and post-filling in injection molding. The governing equations are in terms of compressible, non-isothermal fluid, and the constitutive equation is based on the Phan–Thien–Tanner model. By introducing some hypotheses according to the characteristics of injection molding, a quasi-Poisson type equation about pressure is derived with part integration. Besides, an analytical form of flow-induced stress is also generalized by using the undermined coefficient method. The conventional Galerkin approach is employed to solve the derived pressure equation, and the ‘upwind’ difference scheme is used to discrete the energy equation. Coupling is achieved between velocity and stress by super relax iteration method. The flow in the test mold is investigated by comparing the numerical results and photoelastic photos for polystyrene, showing flow-induced stresses are closely related to melt temperatures. The filling of a two-cavity box is also studied to investigate the viscoelastic effects on real injection molding.      20. Modelling Surface Runoff and Infiltration of Rain by an Elliptic-Parabolic Equation Coupled with a First‐Order Equation on the Boundary      Ján Filo  Stephan Luckhaus 《Archive for Rational Mechanics and Analysis》1999,146(2):157-182 . (Accepted April 9, 1998)

On the minimal speed of traveling waves for a nonlocal delayed reaction–diffusion equation

In this note, we give constructive upper and lower bounds for the minimal speed of propagation of traveling waves for a nonlocal delayed reaction–diffusion equation.     

Tension–torsion fracture experiments – Part II: Simulations with the extended Gurson model and a ductile fracture criterion based on plastic strain

An extension of the Gurson model that incorporates damage development in shear is used to simulate the tension–torsion test fracture data presented in Faleskog and Barsoum (2013) (Part I) for two steels, Weldox 420 and 960. Two parameters characterize damage in the constitutive model: the effective void volume fraction and a shear damage coefficient. For each of the steels, the initial effective void volume fraction is calibrated against data for fracture of notched round tensile bars and the shear damage coefficient is calibrated against fracture in shear. The calibrated constitutive model reproduces the full range of data in the tension–torsion tests thereby providing a convincing demonstration of the effectiveness of the extended Gurson model. The model reinforces the experiments by highlighting that for ductile alloys the effective plastic strain at fracture cannot be based solely on stress triaxiality. For nominally isotropic alloys, a ductile fracture criterion is proposed for engineering purposes that depends on stress triaxiality and a second stress invariant that discriminates between axisymmetric stressing and shear dominated stressing.     

A new flux-limiting approach–based kinetic scheme for the Euler equations of gas dynamics

This paper proposes a new kinetic-theory-based high-resolution scheme for the Euler equations of gas dynamics. The scheme uses the well-known connection that the Euler equations are suitable moments of the collisionless Boltzmann equation of kinetic theory. The collisionless Boltzmann equation is discretized using Sweby's flux-limited method and the moment of this Boltzmann level formulation gives a Euler level scheme. It is demonstrated how conventional limiters and an extremum-preserving limiter can be adapted for use in the scheme to achieve a desired effect. A simple total variation diminishing criteria relaxing parameter results in improving the resolution of the discontinuities in a significant way. A 1D scheme is formulated first and an extension to 2D on Cartesian meshes is carried out next. Accuracy analysis suggests that the scheme achieves between first- and second-order accuracy as is expected for any second-order flux-limited method. The simplicity and the explicit form of the conservative numerical fluxes add to the efficiency of the scheme. Several standard 1D and 2D test problems are solved to demonstrate the robustness and accuracy.     

Computing flow-induced stresses of injection molding based on the Phan–Thien–Tanner model

A numerical approach is introduced to solve the viscoelastic flow problem of filling and post-filling in injection molding. The governing equations are in terms of compressible, non-isothermal fluid, and the constitutive equation is based on the Phan–Thien–Tanner model. By introducing some hypotheses according to the characteristics of injection molding, a quasi-Poisson type equation about pressure is derived with part integration. Besides, an analytical form of flow-induced stress is also generalized by using the undermined coefficient method. The conventional Galerkin approach is employed to solve the derived pressure equation, and the ‘upwind’ difference scheme is used to discrete the energy equation. Coupling is achieved between velocity and stress by super relax iteration method. The flow in the test mold is investigated by comparing the numerical results and photoelastic photos for polystyrene, showing flow-induced stresses are closely related to melt temperatures. The filling of a two-cavity box is also studied to investigate the viscoelastic effects on real injection molding.     

Modelling Surface Runoff and Infiltration of Rain by an Elliptic-Parabolic Equation Coupled with a First‐Order Equation on the Boundary

. (Accepted April 9, 1998)