Modeling Method of C/C-ZrC Composites and Prediction of Equivalent Thermal Conductivity Tensor Based on Asymptotic Homogenization

This article proposes a modeling method for C/C-ZrC composite materials.According to the superposition of Gaussian random field,the original gray model is obtained,and the threshold segmentation method is used to generate the C-ZrC inclusion model.Finally,the fiber structure is added to construct the microstructure of the three-phase plain weave composite.The reconstructed inclusions can meet the randomness of the shape and have a uniform distribution.Using an algorithm based on asymptotic homogenization and finite element method,the equivalent thermal conductivity prediction of the microstructure finite element model was carried out,and the influence of component volume fraction on material thermal properties was explored.The sensitivity of model parameters was studied,including the size,mesh sensitivity,Gaussian complexity,and correlation length of the RVE model,and the optimal calculation model was selected.The results indicate that the volume fraction of the inclusion phase has a significant impact on the equivalent thermal conductivity of the material.As the volume fraction of carbon fiber and ZrC increases,the equivalent thermal conductivity tensor gradually decreases.This model can be used to explore the impact of materialmicrostructure on the results,and numerical simulations have studied the relationship between structure and performance,providing the possibility of designing microstructure based on performance.

Theoretical and experimental investigation of the resonance responses and chaotic dynamics of a bistable laminated composite shell in the dynamic snap-through mode

The dynamic model of a bistable laminated composite shell simply supported by four corners is further developed to investigate the resonance responses and chaotic behaviors.The existence of the 1:1 resonance relationship between two order vibration modes of the system is verified.The resonance response of this class of bistable structures in the dynamic snap-through mode is investigated,and the four-dimensional(4D)nonlinear modulation equations are derived based on the 1:1 internal resonance relationship by means of the multiple scales method.The Hopf bifurcation and instability interval of the amplitude frequency and force amplitude curves are analyzed.The discussion focuses on investigating the effects of key parameters,e.g.,excitation amplitude,damping coefficient,and detuning parameters,on the resonance responses.The numerical simulations show that the foundation excitation and the degree of coupling between the vibration modes exert a substantial effect on the chaotic dynamics of the system.Furthermore,the significant motions under particular excitation conditions are visualized by bifurcation diagrams,time histories,phase portraits,three-dimensional(3D)phase portraits,and Poincare maps.Finally,the vibration experiment is carried out to study the amplitude frequency responses and bifurcation characteristics for the bistable laminated composite shell,yielding results that are qualitatively consistent with the theoretical results.

Editorial: Computational mechanics of granular materials

Most of granular materials are highly heteroge- neous, composed of voids and particles with different sizes and shapes. Geological matter, soil and clay in nature, geo-structure, concrete, etc. are practical ex- amples among them. From the microscopic view, a lo- cal region in the medium is occupied by particles with small but finite sizes and granular material is naturally modeled as an assembly of discrete particles in contacts On the other hand, the local region is identified with a material point in the overall structure and this discon- tinuous medium can then be represented by an effective continuum on the macroscopic level

Editorial:Mechanics of intelligent materials and structures

Recently, intelligent or smart materials and structures have been received more and more attention due to their distinguished multi-field coupling properties and wide applications in aerospace, automobiles, civil structures, medical devices, information storage, energy harvesting and so on. It is of academic challenge to fully understand the complex multi-field coupling behaviors of various smart materials and structures, and of engineering sig- nificance to enhance the performance and reliability of these materials and structures in industrial applications. The papers in the special topic of Mechanics of Intelligent Materials and Structures focus on the understanding of the electromechanical, magneto-elastic, and magneto-rheological coupling behav- iors and properties of smart materials and structures for applications in vibration control, resonators, and various functional devices.

Slip boundary effect on the critical Reynolds number of subcritical transition in channel flow

In this letter,the effect of slip boundary on the origin of subcritical transition in two-dimensional chan-nel flows is studied numerically and theoretically.It is shown that both the positive and the negative slip lengths will increase the critical Reynolds number of localized wave packet and hence postpone the transition.By applying a variable transformation and expanding the variables about a small slip length,it is illustrated that the slip boundary effect only exists in the second and higher order modulations of the no-slip solution,and hence explains the power law found in simulations,i.e.the relative increment of the critical Reynolds number due to the slip boundary is proportional to the square of the slip length.

Theoretical analysis of surface waves in piezoelectric medium with periodic shunting circuits

The investigations of surface waves in the piezoelectric medium bring out great possibility in designing smart surface acoustic wave(SAW)devices.It is important to study the dispersion properties and manipulation mechanism of surface waves in the semi-infinite piezoelectric medium connected with periodic arrangement of shunting circuits.In this study,the extended Stroh formalism is developed to theoretically analyze the dispersion relations of surface waves under different external circuits.The band structures of both the Rayleigh wave and the Bleustein-Gulyaev(BG)wave can be determined and manipulated with proper electrical boundary conditions.Furthermore,the electromechanical coupling effects on the band structures of surface waves are discussed to figure out the manipulation mechanism of adjusting electric circuit.The results indicate that the proposed method can explain the propagation behaviors of surface waves under the periodic electrical boundary conditions,and can provide an important theoretical guidance for designing novel SAW devices and exploring extensive applications in practice.

Numerical simulation of the mechanical behavior of superconducting tape in conductor on round core cable using the cohesive zone model

Cables composed of rare-earth barium copper oxide(REBCO)tapes have been extensively used in various superconducting devices.In recent years,conductor on round core(CORC)cable has drawn the attention of researchers with its outstanding current-carrying capacity and mechanical properties.The REBCO tapes are wound spirally on the surface of CORC cable.Under extreme loadings,the REBCO tapes with layered composite structures are vulnerable,which can lead to degradation of critical current and even quenching of superconducting devices.In this paper,we simulate the deformation of CORC cable under external loads,and analyze the damage inside the tape with the cohesive zone model(CZM).Firstly,the fabrication and cabling of CORC are simulated,and the stresses and strains generated in the tape are extracted as the initial condition of the next step.Then,the tension and bending loads are applied to CORC cable,and the damage distribution inside the tape is presented.In addition,the effects of some parameters on the damage are discussed during the bending simulations.

Controllable rectification on the thermal conductivity of porous YBa_(2)Cu_(3)O_(7−x) superconductors from 3D-printing

Superconducting YBa_(2)Cu_(3)O_(7−x)(YBCO)bulks have promising applications in quasi-permanent magnets,levitation,etc.Recently,a new way of fabricating porous YBCO bulks,named direct-ink-writing(DIW)3D-printing method,has been reported.In this method,the customized precursor paste and programmable shape are two main advantages.Here,we have put forward a new way to customize the YBCO 3D-printing precursor paste which is doped with Al_(2)O_(3)nanoparticles to obtain YBCO with higher thermal conductivity.The great rheological properties of precursor paste after being doped with Al_(2)O_(3)nanoparticles can help the macroscopic YBCO samples with high thermal conductivity fabricated stably with high crystalline and lightweight properties.Test results show that the peak thermal conductivity of Al_(2)O_(3)-doped YBCO can reach twice as much as pure YBCO,which makes a great effort to reduce the quench propagation speed.Based on the microstructure analysis,one can find that the thermal conductivity of Al_(2)O_(3)-doped YBCO has been determined by its components and microstructures.In addition,a macroscopic theoretical model has been proposed to assess the thermal conductivity of different microstructures,whose calculated results take good agreement with the experimental results.Meanwhile,a microstructure with high thermal conductivity has been found.Finally,a macroscopic YBCO bulk with the presented high thermal conductivity microstructure has been fabricated by the Al_(2)O_(3)-doped method.Compared with YBCO fabricated by the traditional 3D-printed,the Al_(2)O_(3)-doped structural YBCO bulks present excellent heat transfer performances.Our customized design of 3D-printing precursor pastes and novel concept of structural design for enhancing the thermal conductivity of YBCO superconducting material can be widely used in other DIW 3D-printing materials.

Mathematical modeling of mixed convective MHD Falkner-Skan squeezed Sutterby multiphase flow with non-Fourier heat flux theory and porosity

In a wide variety of mechanical and industrial applications,e.g.,space cooling,nuclear reactor cooling,medicinal utilizations(magnetic drug targeting),energy generation,and heat conduction in tissues,the heat transfer phenomenon is involved.Fourier’s law of heat conduction has been used as the foundation for predicting the heat transfer behavior in a variety of real-world contexts.This model’s production of a parabolic energy expression,which means that an initial disturbance would immediately affect the system under investigation,is one of its main drawbacks.Therefore,numerous researchers worked on such problem to resolve this issue.At last,this problem was resolved by Cattaneo by adding relaxation time for heat flux in Fourier’s law,which was defined as the time required to establish steady heat conduction once a temperature gradient is imposed.Christov offered a material invariant version of Cattaneo’s model by taking into account the upper-connected derivative of the Oldroyd model.Nowadays,both models are combinedly known as the Cattaneo-Christov(CC)model.In this attempt,the mixed convective MHD Falkner-Skan Sutterby nanofluid flow is addressed towards a wedge surface in the presence of the variable external magnetic field.The CC model is incorporated instead of Fourier’s law for the examination of heat transfer features in the energy expression.A two-phase nanofluid model is utilized for the implementation of nano-concept.The nonlinear system of equations is tackled through the bvp4c technique in the MATLAB software 2016.The influence of pertinent flow parameters is discussed and displayed through different sketches.Major and important results are summarized in the conclusion section.Furthermore,in both cases of wall-through flow(i.e.,suction and injection effects),the porosity parameters increase the flow speed,and decrease the heat transport and the influence of drag forces.

The Boundary Element Method for Ordinary State-Based Peridynamics

The peridynamics(PD),as a promising nonlocal continuum mechanics theory,shines in solving discontinuous problems.Up to now,various numerical methods,such as the peridynamic mesh-free particlemethod(PD-MPM),peridynamic finite element method(PD-FEM),and peridynamic boundary element method(PD-BEM),have been proposed.PD-BEM,in particular,outperforms other methods by eliminating spurious boundary softening,efficiently handling infinite problems,and ensuring high computational accuracy.However,the existing PD-BEM is constructed exclusively for bond-based peridynamics(BBPD)with fixed Poisson’s ratio,limiting its applicability to crack propagation problems and scenarios involving infinite or semi-infinite problems.In this paper,we address these limitations by introducing the boundary element method(BEM)for ordinary state-based peridynamics(OSPD-BEM).Additionally,we present a crack propagationmodel embeddedwithin the framework ofOSPD-BEM to simulate crack propagations.To validate the effectiveness of OSPD-BEM,we conduct four numerical examples:deformation under uniaxial loading,crack initiation in a double-notched specimen,wedge-splitting test,and threepoint bending test.The results demonstrate the accuracy and efficiency of OSPD-BEM,highlighting its capability to successfully eliminate spurious boundary softening phenomena under varying Poisson’s ratios.Moreover,OSPDBEMsignificantly reduces computational time and exhibits greater consistencywith experimental results compared to PD-MPM.

INFLUENCE OF NOSE PERTURBATIONS ON BEHAVIORS OF ASYMMETRIC VORTICES OVER SLENDER BODY

The influence of nose perturbations on the behaviors of asymmetric vortices over a slender body with a three-caliber ogive nose is studied in this paper. The tests of a nose-disturbed slender body with surface pressure measurement were conducted at a low speed wind tunnel with subcritical Reynolds number of 1×105 at angle of attack α=50°. The experiment results show that the behaviors and structure of asymmetric vortices over the slender body are mainly controlled by manual perturbation on the nose of body as compared with geometrical minute irregularities on the test model from the machining tolerances. The effect of the perturbation axial location on asymmetric vortices is the strongest if its location is near the model apex. There are four sensitive circumferential locations of manual perturbation at which bistable vortices over the slender body are switched by the perturbation. The flowfield near the reattachment line of lee side is more sensitive to the perturbation, because the saddle point to saddle point topological structure in this reattachment flowfield is unstable. Various types of perturbation do not change the perturbation effect on the behaviors of bistable asymmetric vortices.

Space-time evolution rules of acoustic emission location of unloaded coal sample at different loading rates

By using MTS815 rock mechanics test system,a series of acoustic emission(AE) location experiments were performed under unloading confining pressure,increasing the axial stress.The AE space-time evolution regularities and energy releasing characteristics during deformation and failure process of coal of different loading rates are compared,the influence mechanism of loading rates on the microscopic crack evolution were studied,combining the AE characteristics and the macroscopic failure modes of the specimens,and the precursory characteristics of coal failure were also analyzed quantitatively.The results indicate that as the loading rate is higher,the AE activity and the main fracture will begin earlier.The destruction of coal body is mainly the function of shear strain at lower loading rate and tension strain at higher rate,and will transform from brittleness to ductility at critical velocities.When the deformation of the coal is mainly plasticity,the amplitude of the AE ringing counting rate increases largely and the AE energy curves appear an obvious ''step'',which can be defined as the first failure precursor point.Statics of AE information shows that the strongest AE activity begins when the axial stress level was 92-98%,which can be defined as the other failure precursor point.As the loading rate is smaller,the coal more easily reaches the latter precursor point after the first one,so attention should be aroused to prevent dynamic disaster in coal mining when the AE activity reaches the first precursor point.

UNSTEADY FLOWS OF A GENERALIZED SECOND GRADE FLUID WITH THE FRACTIONAL DERIVATIVE MODEL BETWEEN TWO PARALLEL PLATES

The fractional calculus approach in the constitutive relationship model of a generalized second grade fluid is introduced.Exact analytical solutions are obtained for a class of unsteady flows for the generalized second grade fluid with the fractional derivative model between two parallel plates by using the Laplace transform and Fourier transform for fractional calculus.The unsteady flows are generated by the impulsive motion or periodic oscillation of one of the plates.In addition,the solutions of the shear stresses at the plates are also determined.

Effects of constitutive parameters on adiabatic shear localization for ductile metal based on JOHNSON-COOK and gradient plasticity models

By using the widely used JOHNSON-COOK model and the gradient-dependent plasticity to consider microstructural effect beyond the occurrence of shear strain localization,the distributions of local plastic shear strain and deformation in adiabatic shear band(ASB)were analyzed.The peak local plastic shear strain is proportional to the average plastic shear strain,while it is inversely proportional to the critical plastic shear strain corresponding to the peak flow shear stress.The relative plastic shear deformation between the top and base of ASB depends on the thickness of ASB and the average plastic shear strain.A parametric study was carried out to study the influence of constitutive parameters on shear strain localization.Higher values of static shear strength and work to heat conversion factor lead to lower critical plastic shear strain so that the shear localization is more apparent at the same average plastic shear strain.Higher values of strain-hardening exponent,strain rate sensitive coefficient,melting point,thermal capacity and mass density result in higher critical plastic shear strain,leading to less apparent shear localization at the same average plastic shear strain.The strain rate sensitive coefficient has a minor influence on the critical plastic shear strain,the distributions of local plastic shear strain and deformation in ASB.The effect of strain-hardening modulus on the critical plastic shear strain is not monotonous.When the maximum critical plastic shear strain is reached,the least apparent shear localization occurs.

Analysis of localized shear deformation of ductile metal based on gradient-dependent plasticity

Shear localization in linear strain softening heterogeneous material under simple shear was investigated analytically. The closed-form solutions obtained based on gradient plasticity theory considering interactions and interplaying among microstructures due to heterogeneity of metal material show that in the normal direction of shear band, elastic shear displacement is linear; while plastic and total shear displacement are non-linear. Elastic shear strain in the band is uniform and the non-uniformity of total shear displacement stems from localized plastic shear displacement. In the center of the band, plastic and total shear displacement all reach their maximum values. In strain-softening process, elastic displacement decreases as flow shear stress decreases. Contrarily, plastic and total shear displacement increase and manifest shear localization occurs progressively. Under the same shear stress level, plastic and total shear displacement increase as strain softening modulus and elastic shear modulus decrease. The present analytical solutions were compared with many experimental results and the agreement is good.

Level set method for numerical simulation of a cavitation bubble,its growth, collapse and rebound near a rigid wall

A level set method of non-uniform grids is used to simulate the whole evolution of a cavitation bubble, including its growth, collapse and rebound near a rigid wall. Single-phase Navier-Stokes equation in the liquid region is solved by MAC projection algorithm combined with second-order ENO scheme for the advection terms. The moving inter-face is captured by the level set function, and the interface velocity is resolved by ＂one-side＂ velocity extension from the liquid region to the bubble region, complementing the second-order weighted least squares method across the interface and projection inside bubble. The use of non-uniform grid overcomes the difficulty caused by the large computational domain and very small bubble size. The computation is very stable without suffering from large flow-field gradients, and the results are in good agreements with other studies. The bubble interface kinematics, dynamics and its effect on the wall are highlighted, which shows that the code can effectively capture the ＂shock wave＂-like pressure and velocity at jet impact, toroidal bubble, and complicated pressure structure with peak, plateau and valley in the later stage of bubble oscillating.

Volume Change of Heterogeneous Quasi-brittle Materials in Uniaxial Compression

The volumetric strain was categorized into elastic and plastic parts. The farmer camposed of axial and lateral strains is uniform and determined by Hooke＇s law ; however, the latter consisting of axial and lateral strains is a fuaction af thickness af shear band determined by grndieat-dependeat plasticity by cansidering the heterngeneity of quasi- brittle materials. The non- uniform lateral strain due to the fact that shear band was farmed in the middle of specimen was averaged within specimen to precisely assess the volumetric strain. Then, the analytical expression for volumetric strain was verified by comparison with two earlier experimental results for concrete and rack. Finally, a detailed parametric study was carried out to investigate effects of constitutive parameters （ shear band thickness, elastic and softening rnoduli ） and geometrical size of specimen（ height and width of specimen ） on the volume dilatancy.

Calculation of temperature distribution in adiabatic shear band based on gradient-dependent plasticity

A method for calculation of temperature distribution in adiabatic shear band is proposed in terms of gradient-dependent plasticity where the characteristic length describes the interactions and interplaying among microstructures. First, the increment of the plastic shear strain distribution in adiabatic shear band is obtained based on gradient-dependent plasticity. Then, the plastic work distribution is derived according to the current flow shear stress and the obtained increment of plastic shear strain distribution. In the light of the well-known assumption that 90% of plastic work is converted into the heat resulting in increase in temperature in adiabatic shear band, the increment of the temperature distribution is presented. Next, the average temperature increment in the shear band is calculated to compute the change in flow shear stress due to the thermal softening effect. After the actual flow shear stress considering the thermal softening effect is obtained according to the Johnson-Cook constitutive relation, the increment of the plastic shear strain distribution, the plastic work and the temperature in the next time step are recalculated until the total time is consumed. Summing the temperature distribution leads to rise in the total temperature distribution. The present calculated maximum temperature in adiabatic shear band in titanium agrees with the experimental observations. Moreover, the temperature profiles for different flow shear stresses are qualitatively consistent with experimental and numerical results. Effects of some related parameters on the temperature distribution are also predicted.

Analysis of damage localization for ductile metal in process of shear band propagation

Distribution of localized damage in shear band can’ t be predicted theoretically based on classical elastoplastic theory. The average damage variable in shear band was considered to be a non-local variable. Based on non-local theory, an analytical expression for the localized damage in strain-softening region of shear band in the process of shear band propagation was presented using boundary condition and symmetry of local damage variable, etc. The results show that dynamic shear softening modulus, dynamic shear strength and shear elastic modulus influence the distribution of the localized damage in shear band. Internal length of ductile metal only governs the thickness of shear band. In the strain-softening region of shear band, the local damage variable along shear band’s tangential and normal directions is non-linear and highly non-uniform. The non-uniformities in the normal and tangential directions of shear band stem from the interactions and interplaying among microstructures and the non-uniform distribution of shear stress, respectively. At the tail of the strain-softening region, the maximum value of local damage variable reaches 1. This means that material at this position fractures completely. At the tip of shear band and upper as well as lower boundaries, no damage occurs. Local damage variable increases as dynamic shear softening modulus decreases or shear elastic modulus increases, leading to difficulty in identification or detection of damage for less ductile metal material at higher strain rates.

Analysis of progressive failure of pillar and instabilitycriterion based on gradient-dependent plasticity

A mechanical model for strain softening pillar is proposed considering the characteristics of progressive shear failure and strain localization. The pillar undergoes elastic, strain softening and slabbing stages. In the elastic stage, vertical compressive stress and deformation at upper end of pillar are uniform, while in the strain softening stage there appears nonuniform due to occurrence of shear bands, leading to the decrease of load-carrying capacity. In addition, the size of failure zone increases in the strain softening stage and reaches its maximum value when slabbing begins. In the latter two stages, the size of elastic core always decreases. In the slabbing stage, the size of failure zone remains a constant and the pillar becomes thinner. Total deformation of the pillar is derived by linearly elastic Hookes law and gradient-dependent plasticity where thickness of localization band is determined according to the characteristic length. Post-peak stiffness is proposed according to analytical solution of averaged compressive stress-average deformation curve. Instability criterion of the pillar and roof strata system is proposed analytically (using) instability condition given by Salamon. It is found that the constitutive parameters of material of pillar, the geometrical size of pillar and the number of shear bands influence the stability of the system; stress gradient controls the starting time of slabbing, however it has no influence on the post-peak stiffness of the pillar.