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.

Strain concentration factor of heterogeneous materials and analytical influence functions based on Eshelby tensor

In this manuscript,Eshelby tensor is employed to assess the strain concentration that arises in the matrix phase at the interface,offering precise values and locations of maximum strain under specific loading conditions for both spherical and cylindrical inclusions.When compared to numerical simulation results,the analytical predictions grounded in the Eshelby tensor exhibit satisfactory accuracy.Then an analytical calculation method based on Eshelby tensor for the elastic strain influence functions of reduced-order homogenization(ROH)method is developed and adopted on particle-reinforced and fibrous composites,presenting its feasibility and advantage on off-line stage calculation of ROH method.The error analyses between analytical and numerical results are conducted.The numerical results also exhibit the necessity of finer interface partitioning to obtain the response on micro-scale with higher resolution.

Active tuning of anisotropic phonon polaritons in natural van der Waals crystals with negative permittivity substrates and its application in energy transport

Phonon polaritons(PhPs)exhibit directional in-plane propagation and ultralow losses in van der Waals(vdW)crystals,offering new possibilities for controlling the flow of light at the nanoscale.However,these PhPs,including their directional propagation,are inherently determined by the anisotropic crystal structure of the host materials.Although in-plane anisotropic PhPs can be manipulated by twisting engineering,such as twisting individual vdW slabs,dynamically adjusting their propagation presents a significant challenge.The limited application of the twisted bilayer structure in bare films further restricts its usage.In this study,we present a technique in which anisotropic PhPs supported by bare biaxial vdW slabs can be actively tuned by modifying their local dielectric environment.Excitingly,we predict that the iso-frequency contour of PhPs can be reoriented to enable propagation along forbidden directions when the crystal is placed on a substrate with a moderate negative permittivity.Besides,we systematically investigate the impact of polaritonic coupling on near-field radiative heat transfer(NFRHT)between heterostructures integrated with different substrates that have negative permittivity.Our main findings reveal that through the analysis of dispersion contour and photon transmission coefficient,the excitation and reorientation of the fundamental mode facilitate increased photon tunneling,thereby enhancing heat transfer between heterostructures.Conversely,the annihilation of the fundamental mode hinders heat transfer.Furthermore,we find the enhancement or suppression of radiative energy transport depends on the relative magnitude of the slab thickness and the vacuum gap width.Finally,the effect of negative permittivity substrates on NFRHT along the[001]crystalline direction ofα-MoO3 is considered.The spectral band where the excited fundamental mode resulting from the negative permittivity substrates is shifted to the first Reststrahlen Band(RB 1)ofα-MoO_(3) and is widened,resulting in more significant enhancement of heat flux from RB 1.We anticipate our results will motivate new direction for dynamical tunability of the PhPs in photonic devices.

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.

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.

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.

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 Hookes 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.

CHARACTERISTICS OF STRESS-INDUCED TRANSFORMATION AND MICROSTRUCTURE EVOLUTION IN Cu-BASED SMA

The mechanical behavior of shape memory alloys （SMAs） is closely related to the formation and evolution of its microstructures. Through theoretical analysis and experimental observations, it was found that the stress-induced martensitic transformation process of single crystal Cu-based SMA under uniaxial tension condition consisted of three periods： nucleation, mixed nucleation and growth, and merging due to growth. During the nucleation, the stress dropped rapidly and the number of interfaces increased very fast while the phase fraction increased slowly. In the second period, both the stress and the interface number changed slightly but the phase fraction increased dramatically. Finally, the stress and the phase fraction changed slowly while the number of interfaces decreased quickly. Moreover, it was found that the transformation could be of multi-stage： sharp stress drops at several strains and correspondingly, the nucleation and growth process occurred quasi-independently in several parts of the sample.

Effects of electroacupuncture on c-Fos expression in the spinal cord and brain of rats with chronic visceral hypersensitivity

BACKGROUND： Visceral hypersensitivity is the main cause of irritable bowel syndrome, c-Fos is a marker of visceral hypersensitivity in the central nervous system. Electroacupuncture can relieve chronic visceral hypersensitivity in rats, but the mechanism is still unknown. OBJECTIVE： To identify c-Fos expression in the spinal cord and cerebral cortex of rats with chronic visceral hypersensitivity, and to test the effects of electroacupuncture on pain sensitivity in rats with chronic visceral hypersensitivity. DESIGN, TIME AND SETTING： A randomized controlled animal experiment was performed at the Animal E：~perimental Center, Shanghai University of Traditional Chinese Medicine, from January to April, 2007. MATERIALS： A total of 24 neonatal, male, Sprague Dawley rats, aged five days old, were equally and randomly assigned into a normal group, a model group, and an electroacupuncture group. Rabbit anti-rat c-Fos antibody and Evision secondary antibody kits （Sigma, USA）, diaminobenzidine kit （Dako, Denmark）, and an LD202H electroacupuncture apparatus （Huawei, Beijing, China） were used in this study. METHODS： Neonatal rats from the model and electroacupuncture groups were used to establish rat models of chronic visceral hypersensitivity by the saccule stimulation method. After model establishment, 0.25 mm diameter electric needles were inserted into Tianshu （ST 25） and Shangjuxu （ST37） at a depth of approximately 0.5 cm, with an square wave （alternating current frequency at 100/20 Hz, amplitude ranged 0.2-0.6 ms, intensity at 1 mA） once for 20 minutes, once a day, for seven days. Rats in the normal and model groups were not treated. MAIN OUTCOME MEASURES： Following 7 days of treatment, c-Fos expression in the spinal cord and cerebral cortex was detected by immunohistochemistry. After the first electroacupuncture treatment, abdominal withdrawal reflex scores were investigated to evaluate the pain threshold for chronic visceral hypersensitivity in rats. RESULTS： Visceral hypersensitivity increased c-Fos staining （P 〈 0.05）, and electroacupuncture significantly decreased the number of these cells to near normal levels （P 〉 0.05）. Abdominal withdrawal reflex scores were significantly lower in the electroacupuncture and normal groups than in the model group （P 〈 0.05） and were similar between the electroacupuncture and normal groups （P 〉 0.05）. CONCLUSION： Electroacupuncture decreases c-Fos expression in the spinal cord and cerebral cortex and increases pain threshold in a chronic visceral hypersensitivity model in rats.

Quantitative calculation of local shear deformation in adiabatic shear band for Ti-6Al-4V

JOHNSON-COOK(J-C) model was used to calculate flow shear stress—shear strain curve for Ti-6Al-4V in dynamic torsion test. The predicted curve was compared with experimental result. Gradient-dependent plasticity(GDP) was introduced into J-C model and GDP was involved in the measured flow shear stress—shear strain curve, respectively, to calculate the distribution of local total shear deformation(LTSD) in adiabatic shear band(ASB). The predicted LTSDs at different flow shear stresses were compared with experimental measurements. J-C model can well predict the flow shear stress—shear strain curve in strain-hardening stage and in strain-softening stage where flow shear stress slowly decreases. Beyond the occurrence of ASB, with a decrease of flow shear stress, the increase of local plastic shear deformation in ASB is faster than the decrease of elastic shear deformation, leading to more and more apparent shear localization. According to the measured flow shear stress—shear strain curve and GDP, the calculated LTSDs in ASB are lower than experimental results. At earlier stage of ASB, though J-C model overestimates the flow shear stress at the same shear strain, the model can reasonably assess the LTSDs in ASB. According to the measured flow shear stress—shear strain curve and GDP, the calculated local plastic shear strains in ASB agree with experimental results except for the vicinity of shear fracture surface. In the strain-softening stage where flow shear stress sharply decreases, J-C model cannot be used. When flow shear stress decreases to a certain value, shear fracture takes place so that GDP cannot be used.

Entire deformational characteristics and strain localization of jointed rock specimen in plane strain compression

Shear band (SB), axial, lateral and volumetric strains as well as Poisson’s ratio of anisotropic jointed rock specimen (JRS) were modeled by Fast Lagrangian Analysis of Continua (FLAC). Failure criterion of rock was a composited Mohr-Coulomb criterion with tension cut-off. An inclined joint was treated as square elements of ideal plastic material beyond the peak strength. Several FISH functions were written to automatically find the addresses of elements in the joint and to calculate the entire deformational characteristics of plane strain JRS. The results show that for moderate joint inclination (JI), strain is only concentrated into the joint governing the behavior of JRS, leading to ideal plastic responses in axial and lateral directions. For higher JI, the post-peak stress-axial and lateral strain curves become steeper as JI increases owing to the increase of new SB’s length. Lateral expansion and precursor to the unstable failure are the most apparent, resulting in the highest Poisson’s ratio and even negative volumetric strain. For lower JI, the entire post-peak deformational characteristics are independent of JI. The lowest lateral expansion occurs, leading to the lowest Poisson’s ratio and positive volumetric strain all along. The present prediction on anisotropic strength in plane strain compression qualitatively agrees with the results in triaxial tests of rocks. The JI calculated by Jaeger’s formula overestimates that related to the minimum strength. Advantages of the present numerical model over the Jaeger’s model are pointed out.

DAMAGE AND FRACTURE EVALUATION OF GRANULAR COMPOSITE MATERIALS BY DIGITAL IMAGE CORRELATION METHOD

This paper presents the applications of digital image correlation technique to the mesoscopic damage and fracture study of some granular based composite materials including steel- fiber reinforced concrete,sandstone and crystal-polymer composite.The deformation fields of the composite materials resulted from stress localization were obtained by the correlation computation of the surface images with loading steps and thus the related damage prediction and fracture parameters were evaluated.The correlation searching could be performed either directly based on the gray levels of the digital images or from the wavelet transform(WT)coefficients of the transform spectrum.The latter was developed by the authors and showed higher resolution and sensitivity to the singularity detection. Because the displacement components came from the rough surfaces of the composite materials without any coats of gratings or fringes of optical interferometry,both surface profiles and the deformation fields of the composites were visualized which was helpful to compare each other to analyze the damage of those heterogeneous materials.

Numerical simulation of influence of shear dilatancy on deformation characteristics of shear band-elastic body system

The paper was numerically focused on investigation of deformation, failure and instability of shear band-surrounding elastic rock system in plane strain direct shear test considering shear dilatancy according to fast lagrangian analysis of continua (FLAC). The adopted failure criterion was a composite Mohr-Coulomb criterion with tension cut-off and post-peak constitutive relation of rock, i.e. linear strain-softening. Numerical results show that dilation angle affects the responses of elements, the number and the position of yielded elements. Increasing dilation angle results in higher load-carrying capability of elements, higher deformation or strain corresponding to peak stress, less brittle post-peak stress-deformation curve. Strain-hardening behavior can occur if dilation angle is high. Therefore, shear band-elastic rock body system tends to be stable and rock burst does not occur as dilation angle is increased. Moreover, the number of yielded elements is in- creased with dilation angle increase and two parallel plastic zones initially generated in the middle of two loading ends of specimen no longer remain collinear, reflecting increase in deformation resistant of the system. Evolution of volumetric strain rate was investigated based on three-dimensional curved surface diagrams. Approximately, volumetric strain rate concentration regions coincide with plastic zones. Volumetric strain rate in yielded elements is very evident. However, in elastic zones volumetric strain has not been nearly changed throughout the numerical tests.

The effects of compressibility and strength on penetration of long rod and jet

The approximate compressible model is adopted to study the effects of strength and compressibility on the penetration by WHA long rod and copper jet into semi-infinite target in detail. For WHA rod penetrating PMMA at 2 km/s <V <5 km/s, the compressibility has a significant effect on the penetration efficiency. We clarify how compressibility affects the penetration efficiency by changing the stagnation pressures of the rod and target. For WHA rod penetrating 4340 Steel and 6061-T6 Al at 2 km/s < V < 10 km/s, the effect of strength is strong and the effect of compressibility is negligible at lower impact velocity, whilst the effect of strength is weak and the effect of compressibility becomes stronger at higher impact velocity. For the copper jet penetrating 4030 Steel, 6061-T6 Al and PMMA. the virtual origin model is adopted, and the compressibility and strength are implicitly considered by the linear relation between the penetration velocity and impact velocity. The effects of compressibility and target resistance on penetration efficiency are studied. The results show that the target resistance has a significant effect on the penetration efficiency. Howver PMMA is much more compressible than copper and the huge difference of compressibility has a significant effect on the penetration by hypervelocity copper jet into PMMA.

Delamination strength of HTS tape under transverse tensile stress and its enhancement by using different Ag layer depositing temperatures

For the application of second generation high temperature superconducting coated conductors (CCs) with layered structures, thermal mismatch between different components and electromagnetic force exerted in superconducting layer in a working magnet can cause transverse tensile stress, which would result in delamination behavior. Therefore many research groups have designed experiments to measure the delamination strength and dedicate to improving that. However, the reason of the discrete distribution of measured data has still not get quantitatively studied, besides, there are lack of investigations on the method of changing depositing conditions to improve the delamination strength except by adding an additional metal layer. In this work, we adopt an anvil test device and obtain delamination strengths as 29.6 MPa of YBa2Cu3O7-x (YBCO)/buffer and 114.6 MPa of buffer/substrate by combing energy disperse spectroscopy (EDS) detection. The reason of discretized measurement data on the delamination strength is explained. Moreover, we find that different temperatures during Ag deposition determine the bonding force of Ag and YBCO layer. The Delamination strength between Ag and YBCO layer increases from 4.4 MPa to larger than 114.6 MPa with temperature elevated from 30℃ to 100℃. Hence we present a novel method for improving the delamination strength of YBCO CCs by setting an optimal temperature of Ag deposition.

Grain Size Effect on Electrical Conductivity and Giant Magnetoresistance of Bulk Magnetic Polycrystals

By solving the Boltzmann transport equation and considering the spin-dependent grain boundary scattering, the distribution of electrons in grains and the electrical transport properties in the applied magnetic field are studied. With regard to the dominant influence of grain boundary scattering which is taken as a boundary condition for the electrical transport, the grain size-dependent electrical conductivity is investigated. In addition, the reorientation of the relative magnetization between grains brings the change of the electron spin when the magnetonanocrystailine material is subjected to the magnetic field, resulting in the remarkable giant magnetoresistance effect.

Shear stress distribution and characteristics of deformation for shear band-elastic body system at pre-peak and post-peak

The distributed shear stress and the displacement across shear band, the evolution of plastic zones, and the load-carrying capacity of rock specimen were investigated in plane strain direct shear test according to Fast Lagrangian Analysis of Continua （FLAC）. And then the shear displacement distribution in normal direction of system composed of localized shear band and elastic rock was analyzed based on gradient-dependent plasticity. The adopted failure criterion was a composite of Mohr-Coulomb criterion, that is, the relation between tension cut-off and postpeak constitutive of rock was linear strain-softening. Numerical results show that shear stress field approximately undergoes three different stages. At first, shear stress is only concentrated in the middle of top and base of specimen. Next, shear stress in the middle of specimen tends to increase, owing to superposition of shear stresses. Interestingly, two peaks of shear stress appear far from the loading ends of specimen, and the peaks approach with the increase in timestep until elements at the center of specimen yield. Finally, relatively lower shear stress level is reached in large part of specimen except in the regions near the two ends. As flow stress decreases, the analytical shear displacement distribution in shear band based on gradient-dependent plasticity becomes steeps outside the band, it is linear and its slope tends to decrease. These theoretical results qualitatively agree with that of the present numerical predicted results. Main advantage of the analytical solution over the numerical results according to FLAC is that it is continuous, smooth and non-linear （except at elastic stage）.

Numerical simulation of solid tumor angiogenesis with Endostatin treatment:a combined analysis of inhibiting effect of anti-angiogenic factor and micro mechanical environment of extracellular matrix

To investigate the influence of anti-angiogenesis drug Endostatin on solid tumor angiogenesis, a mathematical model of tumor angiogenesis was developed with combined influences of local extra-cellular matrix mechanical environment, and the inhibiting effects of Angiostatin and Endostatin. Simulation results show that Angiostatin and Endostatin can effectively inhibit the process of tumor angiogenesis, and decrease the number of blood vessels in the tumor. The present model could be used as a valid theoretical method in the investigation of anti-angiogenic therapy of tumors.