Nonlinear dynamics of a circular curved cantilevered pipe conveying pulsating fluid based on the geometrically exact model

Due to the novel applications of flexible pipes conveying fluid in the field of soft robotics and biomedicine,the investigations on the mechanical responses of the pipes have attracted considerable attention.The fluid-structure interaction(FSI)between the pipe with a curved shape and the time-varying internal fluid flow brings a great challenge to the revelation of the dynamical behaviors of flexible pipes,especially when the pipe is highly flexible and usually undergoes large deformations.In this work,the geometrically exact model(GEM)for a curved cantilevered pipe conveying pulsating fluid is developed based on the extended Hamilton's principle.The stability of the curved pipe with three different subtended angles is examined with the consideration of steady fluid flow.Specific attention is concentrated on the large-deformation resonance of circular pipes conveying pulsating fluid,which is often encountered in practical engineering.By constructing bifurcation diagrams,oscillating shapes,phase portraits,time traces,and Poincarémaps,the dynamic responses of the curved pipe under various system parameters are revealed.The mean flow velocity of the pulsating fluid is chosen to be either subcritical or supercritical.The numerical results show that the curved pipe conveying pulsating fluid can exhibit rich dynamical behaviors,including periodic and quasi-periodic motions.It is also found that the preferred instability type of a cantilevered curved pipe conveying steady fluid is mainly in the flutter of the second mode.For a moderate value of the mass ratio,however,a third-mode flutter may occur,which is quite different from that of a straight pipe system.

On the thermodynamics of plasticity during quasi-isentropic compression of metallic glass

Entropy production in quasi-isentropic compression (QIC) is critically important for understanding the properties of materials under extremeconditions. However, the origin and accurate quantification of entropy in this situation remain long-standing challenges. In this work, a framework is established for the quantification of entropy production and partition, and their relation to microstructural change in QIC. Cu50Zr50is taken as a model material, and its compression is simulated by molecular dynamics. On the basis of atomistic simulation-informed physicalproperties and free energy, the thermodynamic path is recovered, and the entropy production and its relation to microstructural change aresuccessfully quantified by the proposed framework. Contrary to intuition, entropy production during QIC of metallic glasses is relativelyinsensitive to the strain rate ˙γ when ˙γ ranges from 7.5 × 10^(8) to 2 × 10^(9)/s, which are values reachable in QIC experiments, with a magnitudeof the order of 10^(−2)kB/atom per GPa. However, when ˙γ is extremely high (>2 × 10^(9)/s), a notable increase in entropy production rate with˙γ is observed. The Taylor–Quinney factor is found to vary with strain but not with strain rate in the simulated regime. It is demonstrated thatentropy production is dominated by the configurational part, compared with the vibrational part. In the rate-insensitive regime, the increase inconfigurational entropy exhibits a linear relation to the Shannon-entropic quantification of microstructural change, and a stretched exponential relation to the Taylor–Quinney factor. The quantification of entropy is expected to provide thermodynamic insights into the fundamentalrelation between microstructure evolution and plastic dissipation.

Analysis and design for the comprehensive ballistic and blast resistance of polyurea-coated steel plate

Fragments and blast waves generated by explosions pose a serious threat to protective structures.In this paper,the impact resistance of polyurea-coated steel plate under complex dynamic loading is analyzed and designed for improving comprehensive ballistic and blast resistance using the newly established computational evaluating model.Firstly,according to the thickness and placement effects of the coating on the impact resistance,the steel-core sandwich plates are designed,which are proved to own outstanding comprehensive ballistic and blast resistance.Besides,the distribution diagram of ballistic and blast resistance for different polyurea-coated steel plates is given to guide the design of protective structures applying in different explosion scenarios.Furthermore,the dynamic response of designed plates under two scenarios with combined fragments and blast loading is studied.The results show that the synergistic effect of the combined loading reduces both the ballistic and blast resistance of the polyurea-coated steel plate.Besides,the acting sequence of the fragments and blast affects the structural protective performance heavily.It is found that the first loading inducing structural large deformation or damage is dominant.When fragments impact first,the excellent unit-thickness ballistic performance of the structural front part is strongly needed for improving the comprehensive ballistic and blast resistance.

Large eddy simulation of aircraft wake vortex with self-adaptive grid method

A self-adaptive-grid method is applied to numerical simulation of the evolu- tion of aircraft wake vortex with the large eddy simulation （LES）. The Idaho Falls （IDF） measurement of run 9 case is simulated numerically and compared with that of the field experimental data. The comparison shows that the method is reliable in the complex atmospheric environment with crosswind and ground effect. In addition, six cases with different ambient atmospheric turbulences and Brunt V^iis/il^i （BV） frequencies are com- puted with the LES. The main characteristics of vortex are appropriately simulated by the current method. The onset time of rapid decay and the descending of vortices are in agreement with the previous measurements and the numerical prediction. Also, sec-ondary structures such as baroclinic vorticity and helical structures are also simulated. Only approximately 6 million grid points are needed in computation with the present method, while the number can be as large as 34 million when using a uniform mesh with the same core resolution. The self-adaptive-grid method is proved to be practical in the numerical research of aircraft wake vortex.

Extension of covariant derivative(Ⅰ): From component form to objective form

This paper extends the covariant derivative un der curved coordinate systems in 3D Euclid space. Based on the axiom of the covariant form invariability, the classical covariant derivative that can only act on components is ex tended to the generalized covariant derivative that can act on any geometric quantity including base vectors, vectors and tensors. Under the axiom, the algebra structure of the gen eralized covariant derivative is proved to be covariant dif ferential ring. Based on the powerful operation capabilities and simple analytical properties of the generalized covariant derivative, the tensor analysis in curved coordinate systems is simplified to a large extent.

Extension of covariant derivative(Ⅱ): From flat space to curved space

This paper extends the classical covariant deriva tive to the generalized covariant derivative on curved sur faces. The basement for the extension is similar to the pre vious paper, i.e., the axiom of the covariant form invariabil ity. Based on the generalized covariant derivative, a covari ant differential transformation group with orthogonal duality is set up. Through such orthogonal duality, tensor analy sis on curved surfaces is simplified intensively. Under the covariant differential transformation group, the differential invariabilities and integral invariabilities are constructed on curved surfaces.

Application of the quadrilateral area coordinate method:a new element for laminated composite plate bending problems

Recently, some new quadrilateral finite elements were successfully developed by the Quadrilateral Area Coordinate （QAC） method. Compared with those traditional models using isoparametric coordinates, these new models are less sensitive to mesh distortion. In this paper, a new displacement-based, 4-node 20-DOF （5-DOF per node） quadrilateral bending element based on the first-order shear deformation theory for analysis of arbitrary laminated composite plates is presented. Its bending part is based on the element AC-MQ4, a recent-developed high-performance Mindlin-Reissner plate element formulated by QAC method and the generalized conforming condition method; and its in-plane displacement fields are interpolated by bilinear shape functions in isoparametric coordinates. Furthermore, the hybrid post-rocessing procedure, which was firstly proposed by the authors, is employed again to improve the stress solutions, especially for the transverse shear stresses. The resulting element, denoted as AC-MQ4-LC, exhibits excellent performance in all linear static and dynamic numerical examples. It demonstrates again that the QAC method, the generalized conforming condition method, and the hybrid post-processing procedure are efficient tools for developing simple, effective and reliable finite element models.

Extension of covariant derivative(Ⅲ): From classical gradient to shape gradient

This paper further extends the generalized covari ant derivative from the first covariant derivative to the sec ond one on curved surfaces. Through the linear transforma tion between the first generalized covariant derivative and the second one, the second covariant differential transformation group is set up. Under this transformation group, the sec ond class of differential invariants and integral invariants on curved surfaces is made clear. Besides, the symmetric struc ture of the tensor analysis on curved surfaces are revealed.

A Comprehensive Review of Isogeometric Topology Optimization:Methods,Applications and Prospects

Topology Optimization(TO)is a powerful numerical technique to determine the optimal material layout in a design domain,which has accepted considerable developments in recent years.The classic Finite Element Method(FEM)is applied to compute the unknown structural responses in TO.However,several numerical deficiencies of the FEM significantly influence the effectiveness and efficiency of TO.In order to eliminate the negative influence of the FEM on TO,IsoGeometric Analysis(IGA)has become a promising alternative due to its unique feature that the Computer-Aided Design(CAD)model and Computer-Aided Engineering(CAE)model can be unified into a same mathematical model.In the paper,the main intention is to provide a comprehensive overview for the developments of Isogeometric Topology Optimization(ITO)in methods and applications.Finally,some prospects for the developments of ITO in the future are also presented.

Numerical study on the turbulent mixing of planar shock-accelerated triangular heavy gases interface

The interaction of a planar shock wave with a triangle-shaped sulfur hexafluoride (SF6) cylinder surrounded by air is numerically studied using a high resolution finite volume method with minimum dispersion and controllable dissipation reconstruction.The vortex dynamics of the Richtmyer-Meshkov instability and the turbulent mixing induced by the KelvinHelmholtz instability are discussed.A modified reconstruction model is proposed to predict the circulation for the shock triangular gas-cylinder interaction flow.Several typical stages leading the shock-driven inhomogeneity flow to turbulent mixing transition are demonstrated.Both the decoupled length scales and the broadened inertial range of the turbulent kinetic energy spectrum in late time manifest the turbulent mixing transition for the present case.The analysis of variable-density energy transfer indicates that the flow structures with high wavenumbers inside the Kelvin-Helmholtz vortices can gain energy from the mean flow in total.Consequently,small scale flow structures are generated therein by means of nonlinear interactions.Furthermore,the occasional 'pairing' between a vortex and its neighboring vortex will trigger the merging process of vortices and,finally,create a large turbulent mixing zone.

From the second gradient operator and second class of integral theorems to Gaussian or spherical mapping invariants

By combining of the second gradient operator, the second class of integral theorems, the Gaussian-curvature-based integral theorems and the Gaussian （or spherical） mapping, a series of invariants or geometric conservation quantities under Gaussian （or spherical） mapping are revealed. From these mapping invariants important transformations between original curved surface and the spherical surface are derived. The potential applications of these invariants and transformations to geometry are discussed

Analysis on high-temperature oxidation and growth stress of iron-based alloy using phase field method

High-temperature oxidation is an important property to evaluate thermal protection materials. However, since oxidation is a complex process involving microstructure evolution, its quantitative analysis has always been a challenge. In this work, a phase field method （PFM） based on the thermodynamics theory is developed to simulate the oxidation behavior and oxidation induced growth stress. It involves microstructure evolution and solves the problem of quantitatively computational analysis for the oxidation behavior and growth stress. Employing this method, the diffusion process, oxidation performance, and stress evolution axe predicted for Fe-Cr-A1-Y alloys. The numerical results agree well with the experimental data. The linear relationship between the maximum growth stress and the environment oxygen concentration is found. PFM provides a powerful tool to investigate high-temperature oxidation in complex environments.

Gigacycle fatigue behaviors of two SNCM439 steels with different tensile strengthes

Gigacycle fatigue behaviors of two SNCM439 steels with different tensile strengthes were experimentally studied by rotating bending tests,to investigate the effects of the tensile strength obtained by different heat treatment processes on very high cycle fatigue failure mechanisms.The material with higher tensile strength of 1 710 MPa exhibited typical gigacycle fatigue failure characteristics,whereas one with lower tensile strength of 1 010 MPa showed only traditional fatigue limit during the tests and no gigacycle failure could be found even when the specimen ran up to more than 10 8 cycles.Metallographic and fractographic analysis were carried out by an optical microscope （OM） and scanning electron microscope （SEM）.It showed two different crack initiation mechanisms that for the specimen with lower tensile strength the crack prefers surface initiation and for that with higher strength the crack initiates from subsurface inclusions revealed by a fish-eye like microstructure.

Nonlinear effects of line tension in adhesion of small droplets

Three-phase line tensions may become crucial in the adhesion of miero-nano or small droplets on solid planes. In this paper we study for the first time the nonlinear effects in adhesion spanning the full range of physically possible parameters of surface tension, line tension, and droplet size. It is shown that the nonlinear adhesion solution spaces can be characterized into four regions. Within each region the adhesion behaves essentially the same. Especially, inside the characteristic regions with violent nonlinearities, the co-existence of multiple adhesion states for given materials is disclosed. Besides, two common fixed points in the solution space are revealed. These new results are consistent with numerical analysis and experimental observations reported in the literatures.

Characteristic equation solution strategy for deriving fundamental analytical solutions of 3D isotropic elasticity

A simple characteristic equation solution strategy for deriving the fun- damental analytical solutions of 3D isotropic elasticity is proposed. By calculating the determinant of the differential operator matrix obtained from the governing equations of 3D elasticity, the characteristic equation which the characteristic general solution vectors must satisfy is established. Then, by substitution of the characteristic general solution vectors, which satisfy various reduced characteristic equations, into various reduced ad- joint matrices of the differential operator matrix, the corresponding fundamental analyt- ical solutions for isotropic 3D elasticity, including Boussinesq-Galerkin （B-G） solutions, modified Papkovich-Neuber solutions proposed by Min-zhong WANG （P-N-W）, and quasi HU Hai-chang solutions, can be obtained. Furthermore, the independence characters of various fundamental solutions in polynomial form are also discussed in detail. These works provide a basis for constructing complete and independent analytical trial func- tions used in numerical methods.

A direct proof of uniqueness of square-root of a positive semi-definite tensor

Understanding of the basic properties of the positive semi-definite tensor is a prerequisite for its extensive applications in theoretical and practical fields, especially for its square-root. Uniqueness of the square-root of a positive semi-definite tensor is proven in this paper without resorting to the notion of eigenvalues, eigenvectors and the spectral decomposition of the second-order symmetric tensor.

Shape gradient and classical gradient of curvatures:driving forces on micro/nano curved surfaces

Recent experiments and molecule dynamics simulations have shown that adhesion droplets on conical surfaces may move spontaneously and directionally. Besides, this spontaneous and directional motion is independent of the hydrophilicity and hydrophobicity of the conical surfaces. Aimed at this important phenomenon, a gen- eral theoretical explanation is provided from the viewpoint of the geometrization of micro/nano mechanics on curved surfaces. In the extrinsic mechanics on micro/nano soft curved surfaces, we disclose that the curvatures and their extrinsic gradients form the driving forces on the curved spaces. This paper focuses on the intrinsic mechanics on micro/nano hard curved surfaces and the experiment on the spontaneous and directional motion. Based on the pair potentials of particles, the interactions between an isolated particle and a micro/nano hard curved surface are studied, and the geometric foundation for the interactions between the particle and the hard curved surface is analyzed. The following results are derived： （a） Whatever the exponents in the pair potentials may be, the potential of the particle/hard curved surface is always of the unified curvature form, i.e., the potential is always a unified function of the mean curvature and the Gaussian curvature of the curved surface. （b） On the basis of the curvature-based potential, the geometrization of the micro/nano mechanics on hard curved surfaces may be realized. （c） Similar to the extrinsic mechanics on micro/nano soft curved surfaces, in the intrinsic mechanics on micro/nano hard curved surfaces, the curvatures and their intrinsic gradi- ents form the driving forces on the curved spaces. In other words, either on soft curved surfaces or hard curved surfaces and either in the extrinsic mechanics or the intrinsic mechanics, the curvatures and their gradients are all essential factors for the driving forces on the curved spaces. （d） The direction of the driving force induced by the hard curved surface is independent of the hydrophilieity and hydrophobicity of the curved surface, explaining the experimental phenomenon of the spontaneous and directional motion.

An elementary proof for the representation theorem of analytic isotropic tensor functions of a second-order tensor

Based on the Cayley-Hamilton theorem and fixed-point method,we provide an elementary proof for the representation theorem of analytic isotropic tensor functions of a second-order tensor in a three-dimensional(3D)inner-product space,which avoids introducing the generating function and Taylor series expansion.The proof is also extended to any finite-dimensional inner-product space.

Rotational dynamics of bottom-heavy rods in turbulence from experiments and numerical simulations

We successfully perform the three-dimensional tracking in a turbulent fluid flow of small axisymmetrical particles that are neutrally-buoyant and bottom-heavy,i.e.,they have a non-homogeneous mass distribu-tion along their symmetry axis.We experimentally show how a tiny mass inhomogeneity can affect the particle orientation along the preferred vertical direction and modify its tumbling rate.The experiment is complemented by a series of simulations based on realistic Navier-Stokes turbulence and on a point-like particle model that is capable to explore the full range of parameter space characterized by the gravi-tational torque stability number and by the particle aspect ratio.We propose a theoretical perturbative prediction valid in the high bottom-heaviness regime that agrees well with the observed preferential ori-entation and tumbling rate of the particles.We also show that the heavy-tail shape of the probability distribution function of the tumbling rate is weakly affected by the bottom-heaviness of the particles.

Specimen size effect on the fracture energy of architected stretchable materials

Architected stretchable materials with well-organized microarchitectures evolve very rapidly due to their potential in customizing mechanical properties and achieving exotic functions.In many applications,the architected stretchable materials are required to sustain large deformation,and their fracture is size-dependent.However,the size effect on the fracture of architected stretchable materials is still elusive.Here,we study this issue by experiment and finite element calculation.It is found that the fracture energy of architected stretchable materials increases with the specimen size ratio,H/h,within a range.When H/h reaches a transition ratio,Rt,the fracture energy approaches a plateau.This transition ratio differentiates the size-dependent and size-independent fracture behavior of architected stretchable materials.The mechanical properties of constituent material only have a minor effect on the transition ratio.The degree of constraint and stress concentration at the node,which are affected by the geometry of the unit-cell,dominate the specimen size effect.The result gives a practical guidance in choosing the specimen size to measure the steady state fracture energy of this class of materials.This work provides insights into the fracture of architected stretchable materials and design for fractureresistant architected stretchable devices.