Sloshing and energy dissipation in an egg: SPH simulations and experiments

Tall structures, such as towers and bridges, can oscillate at excessive magnitudes when subjected to wind and earthquake loads. Liquid sloshing absorbers can be used to suppress these excessive oscillations by tuning the frequency of the sloshing to the critical frequency of the structure. Sloshing absorbers are simple structures consisting of a partially full container of liquid with a free surface. Tuning ensures that significant amounts of harmful energy can be extracted from the structure to the sloshing liquid. However, there needs to be a rapid means of dissipating this energy to avoid its returning back to the structure (then back to the liquid periodically).A hen׳s egg seems to have evolved to efficiently dissipate energy to protect its embryo using sloshing of its liquid content. Hence, the potential to implement the egg׳s unique properties as a sloshing absorber for structural control, is the main focus of this study. Numerical simulations, using Smoothed Particle Hydrodynamics (SPH), and experimental comparisons are presented in this paper. One objective is to demonstrate the ability of SPH to simulate complex free surface behaviour in three dimensions. Such a tool is then useful to identify different dissipation modes. Effects of fill volume and viscosity on the rate of dissipation, are also investigated.     

Thermoelastic interaction with energy dissipation in a transversely isotropic thin circular disc

The distribution of stresses due to step input of temperature on the boundaries of a homogeneous transversely isotropic circular disc is investigated by applying Laplace transform technique in the context of generalized theories of thermo-elasticity. The inverse of the transformed solution is carried out by applying a method of Bellman et al. The stresses are computed numerically and presented graphically in a number of figures. A comparison of the results for different theories (CTE, CCTE, TRDTE(GL), TEWED(GN)) and the effect of anisotropy on the stresses are also presented. When the material is isotropic and outer radius of the disc tends to infinity, the corresponding result agrees with that of existing literature.     

Effect of particle polydispersity on micromechanical properties and energy dissipation in granular mixtures

A series of numerical tests was conducted to study the micromechanical properties and energy dissipation in polydisperse assemblies of spherical particles subjected to uniaxial compression. In general, distributed particle size assemblies with standard deviations ranging from 0% to 80% of the particle mean diameter were examined. The microscale analyses included the trace of the fabric tensor, magnitude and orientation of the contact forces, trace of stress, number of contacts and degree of mobilization of friction in contacts between particles. In polydisperse samples, the average coordination numbers were lower than in monodisperse assemblies, and the mobilization of friction was higher than in monodisperse assemblies due to the non-uniform spatial rearrangement of spheres in the samples and the smaller displacements of the particles. The effect of particle size heterogeneity on both the energy density and energy dissipation in systems was also investigated.     

Thermoelastic interactions without energy dissipation in an unbounded body with a spherical cavity subjected to harmonically varying temperature

The present work is concerned with the thermally induced vibration in a homogeneous and isotropic unbounded body with a spherical cavity. The Green and Nagdhi model of thermoelasticity without energy dissipation is employed. The closed form solutions for distributions of displacement, temperature and stresses are obtained. The solutions valid in the case of small frequency are deduced and the results are compared with the corresponding results obtained in other generalized thermoelasticity theories. Numerical results applicable to a copper-like material are also presented graphically and the nature of variations of the physical quantities with radial coordinate and with frequency of vibration is analyzed.     

Evaluation of the drag force by integrating the energy dissipation rate in stokes flow for 2D domains using the FEM

The total drag force on the surface of a body, which is the sum of the form drag and the skin friction drag in a 2D domain, is numerically evaluated by integrating the energy dissipation rate in the whole domain for an incompressible Stokes fluid. The finite element method is used to calculate both the energy dissipation rate in the whole domain as well as the drag on the boundary of the body. The evaluation of the drag and the energy dissipation rate are post-processing operations which are carried out after the velocity field and the pressure field for the flow over a particular profile have been obtained. The results obtained for the flow over three different but constant area profiles—a circle, an ellipse and a cross-section of a prolate spheroid—with uniform inlet velocity are presented and it is shown that the total drag force times the velocity is equal to the total energy dissipation rate in the entire finite flow domain. Hence, by calculating the energy dissipation rate in the domain with unit velocity specified at the far-field boundary enclosing the domain, the drag force on the boundary of the body can be obtained.     

Modeling rock specimens through 3D printing: Tentative experiments and prospects

Current developments in 3D printing (3DP) technology provide the opportunity to produce rock-like specimens and geotechnical models through additive manufacturing, that is, from a file viewed with a computer to a real object. This study investigated the serviceability of 3DP products as substitutes for rock specimens and rock-type materials in experimental analysis of deformation and failure in the laboratory. These experiments were performed on two types of materials as follows: (1) compressive experiments on printed sand-powder specimens in different shapes and structures, including intact cylinders, cylinders with small holes, and cuboids with pre-existing cracks, and (2) compressive and shearing experiments on printed polylactic acid cylinders and molded shearing blocks. These tentative tests for 3DP technology have exposed its advantages in producing complicated specimens with special external forms and internal structures, the mechanical similarity of its product to rock-type material in terms of deformation and failure, and its precision in mapping shapes from the original body to the trial sample (such as a natural rock joint). These experiments and analyses also successfully demonstrate the potential and prospects of 3DP technology to assist in the deformation and failure analysis of rock-type materials, as well as in the simulation of similar material modeling experiments.     

Theoretical analysis of electric, magnetic and magnetoelectric properties of nano-structured multiferroic composites

Electric, magnetic and magnetoelectric properties of the nano-structured multiferroic composites were studied by using an energy formulation with the consideration of the surface, interface, and size effect. Coupled thermodynamic evolution equations with respect to the spontaneous polarization and magnetization were established, in which the elastic fields in the matrix and inclusions were solved based on the Eshelby's equivalent inclusion concept and the Mori–Tanaka method. Physical properties of the composite, such as the spontaneous order parameters, piezoelectric/piezomagnetic properties, and the magnetoelectric coupling effect are highly dependent on the stress state and the microstructures of the nano-composites. Magnetoelectric coupling voltage coefficient was unstable in the vicinity of the critical size and disappeared below the critical size. The model is versatile enough for various composite structures.     

Interfacial energy and dissipation in martensitic phase transformations. Part I: Theory

A model of evolving martensitic microstructures is formulated that incorporates the interfacial energy and dissipation on three different scales corresponding to the grain boundaries attained by martensite plates, the interfaces between austenite and martensite plates, and the twin interfaces within martensite plates. Three different time scales are also considered in order to clarify the meaning of rate-independent dissipation related to instabilities at more refined temporal and spatial scales. On the slowest time scale, the process of deformation and martensitic phase transformation is investigated as being composed of segments of smooth quasi-static evolution separated by sudden jumps associated with creation or annihilation of interfaces. A general evolution rule is used in the form of minimization of the incremental energy supply to the whole compound thermodynamic system, including the rate-independent dissipation. Close relationship is shown between the evolution rule and the thermodynamic condition for stability of equilibrium, with no a priori assumption on convexity of the dissipation function. It is demonstrated that the extension of the minimum principle from the first-order rates to small but finite increments requires a separate symmetry restriction imposed on the state derivative of the dissipation function. Formulae for the dissipation associated with annihilation of interfaces are proposed that exhibit limited path-independence and satisfy that symmetry requirement. By exploiting the incremental energy minimization rule with the help of the transport theorems, local propagation conditions are derived for both planar and curved phase transformation fronts. The theory serves as a basis for the algorithm for calculation of the stress-induced evolution of martensitic microstructures along with their characteristic dimensions and related hysteresis loops in shape memory alloys; the examples are given in Part II of the paper.     

A preliminary study of 3D difference scheme with energy dynamic equilibrium

In this paper, a difference scheme with energy dynamic equilibrium (DS-EDE) is presented, which can be used for the simulation of long-term atmosphere and sea motion. Based on three dimensional nonlinear evolution equations for atmosphere and sea motion, a three dimensional compact upwind scheme (CUWS) is constructed, as the basis of the DS-EDE. The DS-EDE satisfies the following condition of energy dynamic equilibrium (EDE): the total work of external forces on the region boundary is equal to the sum of the total effective variation of the kinetic energy and the energy dissipation in the average flow motion and the effective variation of the potential energy per unit time within the region of interest. It really reflects the basic mechanism of the action of external forces and dissipation in atmosphere and sea movement. Therefore, the DS-EDE developed in this paper is a suitable model for simulating long-term atmosphere and sea movement with forcing and dissipation.     

三维编织复合材料渐进损伤的非线性模型及强度分析

建立了考虑周期性位移边界条件的细观体胞模型,对三维编织复合材料的渐进损伤过程进行数值模拟。采用Eshelby-Mori—Tanaka方法计算含损伤裂纹的材料的剐度矩阵,并将有限元网格尺寸和单元裂纹尺寸引入损伤演化方程,有效地降低了模拟结果对有限元网格的依赖程度。通过计算得到了材料应力应变的非线性关系和失效时的极限强度,并分析了材料的破坏机理。结果表明,大编织角材料的破坏模式主要是基体失效与纤维横向拉剪破坏,模拟计算结果与文献中的实验值吻合较好。     

Interfacial energy and dissipation in martensitic phase transformations. Part II: Size effects in pseudoelasticity

This paper is a continuation of the Part I (H. Petryk, S. Stupkiewicz, Interfacial energy and dissipation in martensitic phase transformations. Part I: Theory. J. Mech. Phys. Solids, 2010, doi:10.1016/j.jmps.2009.11.003). A fully three-dimensional model of an evolving martensitic microstructure is examined, taking into account size effects due to the interfacial energy and also dissipation related to annihilation of interfaces. The elastic micro-strain energy at microstructured interfaces is determined with the help of finite element computations and is approximated analytically. Three interface levels are examined: of grain boundaries attained by parallel martensite plates, of interfaces between austenite and twinned martensite, and of twin interfaces within the martensite phase. Minimization of the incremental energy supply, being the sum of the increments in the free energy and dissipation of the bulk and interfacial type at all levels, is used as the evolution rule, based on the theory presented in Part I. An example of the formation and evolution of a rank-three laminated microstructure of finite characteristic dimensions in a pseudoelastic CuAlNi shape memory alloy is examined quantitatively.     

Finite element solution of heat and mass transfer flow with Hall current, heat source, and viscous dissipation

The aim of the paper is to investigate the effect of heat and mass transfer on the unsteady magnetohydrodynamic free convective flow with Hall current, heat source, and viscous dissipation. The problem is governed by the system of coupled non-linear partial differential equations whose exact solution is difficult to obtain. Therefore, the problem is solved by using the Galerkin finite element method. The effects of the various parameters like Hall current, Eckert number, heat source parameter, Prandtl number, and Schmidt number on the velocity components, the temperature, and the concentration are also examined through graphs.     

Nonlinear viscoelastic multi-scale repetitive unit cell model of 3D woven composites with damage evolution

Three-dimensional (3D) textile composites have great potential applications to aircrafts and high speed vehicles because of the high strength/weight ratios and the capabilities of manufacturing complex, net-shape preforms. This paper reports the nonlinear viscoelastic responses and damage mechanisms of one kind of 3D textile composites, named as 3D orthogonal woven composite (3DOWC) under quasi-static tensile loading based on a micro/meso-scale repetitive unit cells (RUCs) model. In the RUCs model, the resin is described with a nonlinear viscoelastic material and the fibers/tows with an elastic material. The damage initiation and propagation in resin are simulated by the post-damage constitutive models with maximum principal theory failure criteria. The fibers/tows impregnated with resin are defined by elastic transverse-isotropic material model with ultimate strengths failure of ‘expanded smeared crack’ both along and perpendicular to fibers/tows axis direction. The engineering parameters and ultimate strengths of homogenized fibers/tows filled with matrix in meso-RUCs model are transferred from the numerical analysis of the micro-RUCs. The results are compared with experimental and theoretical values of RUC deformation and damage initiation and propagation under monotonic axial loading. The methodology of establishing the nonlinear visco-elastic multi-scale model of 3D textile composites without introducing the real fabric architecture in finite element analyses is explained. With the multi-scale RUCs model, the mechanical behaviors of other kinds of 3D textile composites can also be predicted.     

Damping for large-amplitude vibrations of plates and curved panels,Part 1: Modeling and experiments

Theoretical and experimental non-linear vibrations of thin rectangular plates and curved panels subjected to out-of-plane harmonic excitation are investigated. Experiments have been performed on isotropic and laminated sandwich plates and panels with supported and free boundary conditions. A sophisticated measuring technique has been developed to characterize the non-linear behavior experimentally by using a Laser Doppler Vibrometer and a stepped-sine testing procedure. The theoretical approach is based on Donnell's non-linear shell theory (since the tested plates are very thin) but retaining in-plane inertia, taking into account the effect of geometric imperfections. A unified energy approach has been utilized to obtain the discretized non-linear equations of motion by using the linear natural modes of vibration. Moreover, a pseudo arc-length continuation and collocation scheme has been used to obtain the periodic solutions and perform bifurcation analysis. Comparisons between numerical simulations and the experiments show good qualitative and quantitative agreement. It is found that, in order to simulate large-amplitude vibrations, a damping value much larger than the linear modal damping should be considered. This indicates a very large and non-linear increase of damping with the increase of the excitation and vibration amplitude for plates and curved panels with different shape, boundary conditions and materials.     

An explicit finite volume scheme on staggered grids for the Euler equations: Unstructured meshes,stability analysis,and energy conservation

We set up a numerical strategy for the simulation of the Euler equations, in the framework of finite volume staggered discretizations where numerical densities, energies, and velocities are stored on different locations. The main difficulty relies on the treatment of the total energy, which mixes quantities stored on different grids. The proposed method is strongly inspired, on the one hand, from the kinetic framework for the definition of the numerical fluxes, and, on the other hand, from the discrete duality finite volume (DDFV) framework, which has been designed for the simulation of elliptic equations on complex meshes. The time discretization is explicit and we exhibit stability conditions that guaranty the positivity of the discrete densities and internal energies. Moreover, while the scheme works on the internal energy equation, we can define a discrete total energy which satisfies a local conservation equation. We provide a set of numerical simulations to illustrate the behavior of the scheme.     

A Binary Model of textile composites: III high failure strain and work of fracture in 3D weaves

Prior experiments have revealed exceptionally high values of the work of fracture (0.4-) in carbon/epoxy 3D interlock woven composites. Detailed destructive examination of specimens suggested that much of the work of fracture arose when the specimens were strained well beyond the failure of individual tows yet still carried loads . A mechanism of lockup amongst broken tows sliding across the final tensile fracture surface was suggested as the means by which high loads could still be transferred after tow failure. In this paper, the roles of weave architecture and the distribution of flaws in the mechanics of tow lockup are investigated by Monte Carlo simulations using the so-called Binary Model. The Binary Model was introduced previously as a finite element formulation specialised to the problem of simulating relatively large, three-dimensional segments of textile composites, without any assumption of periodicity or other symmetry, while preserving the architecture and topology of the tow arrangement. The simulations succeed in reproducing all qualitative aspects of measured stress-strain curves. They reveal that lockup can indeed account for high loads being sustained beyond tow failure, provided flaws in tows have certain spatial distributions. The importance of the interlock architecture in enhancing friction by holding asperities on sliding fibre tows into firm contact is highlighted.     

Thermally Developing Forced Convection in a Porous Medium: Circular Duct with Walls at Constant Temperature,with Longitudinal Conduction and Viscous Dissipation Effects

A modified Graetz methodology is applied to investigate the thermal development of forced convection in a circular duct filled by a saturated porous medium, with walls held at constant temperature, and with the effects of longitudinal conduction and viscous dissipation included. The Brinkman model is employed. The analysis leads to expressions for the local Nusselt number, as a function of the dimensionless longitudinal coordinate and other parameters (Darcy number, Péclet number, Brinkman number).     

Ideomechanics of transitory and dissipative systems associated with length, velocity, mass and energy

Validation of Ideomechanics (IDM) is manifested by removing the inconsistency of applying open system test data to closed system theories. Instead, the available open system test data can be used rightly to determine the physical parameters of the transitional functions defined by the mean values of length (free path), velocity, mass and energy. Multiscaling and size/time effects are considered. Ambiguities are mitigated when energy takes precedent in lieu of the concept of force.Determined directly from IDM is the energy density function from the velocity that can represent the magnitude of the energy sink and source. The formulation involves grouping pairs of variables of opposing poles that can be constructed as ideograms, much like yin-yang of I-Chin. The flow of Chi implicates the arrow of time and irreversibility. Mass activation/inactivation (AIA) is assumed to be related to the expansion/contraction (EXCO) of matter. Inadvertently, physical systems are identified with inhaling and exhaling of energy corresponding, respectively, to direct-absorption and self-dissipation (DASD). They are postulated to be the basic process for determining the integrity of the system.In contrast to Newtonian/Einsteinian mechanics (NEM) that uses field equations for determining the behavior of the whole everywhere for all time, IDM considers the mean behavior at any given size/time scale, however, large and small. Uncertainties are addressed by the scale transitional functions. The new paradigm can be applied to scaling shifting and to construct equivalence relations for open systems and to the use of existing test data free of ambiguities. Classical conservation laws for closed systems are reducible from the equivalence principles of open systems. The same holds for the classical kinetic molecular theory of matter that can be modified to include dissipation.