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.

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.

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.

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

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.

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.

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.

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.

Vibration characteristics of rotor system with tilting-pad journal bearing of elastic and damped pivots

In view of an entire dynamic model of tilting-pad journal bearing(TPJB) in which the pads swing and vibrate along geometric direction of preload, a TPJB of elastic and damped pivots was designed and manufactured. Vibration experiments were carried out under the conditions of different rotor bending stiffness and oil supply pressure to find out the relationship between the new bearing's vibration depression effect and other dynamic parameters of the rotor. The result shows that critical amplitudes can be efficaciously reduced while system's stability can be remarkably improved by this bearing. Besides, the bearing's effect of vibration depression weakens as the rotor bending stiffness increases, but heightens it as the oil supply pressure increases.

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.

State-vector equation with damping and vibration analysis of laminates

Based on the modified mixed Hellinger-Reissner（H-R） variational principle for elastic bodies with damping, the state-vector equation with parameters is directionally derived from the principle. A new solution for the harmonic vibration of simply supported rectangular laminates with damping is proposed by using the precise integration method and Muller method. The general solutions for the free vibration of underdamping, critical damp and overdamping of composite laminates are given simply in terms of the linear damping vibration theory. The effect of viscous damping force on the vibration of composite laminates is investigated through numerical examples. The state-vector equation theory and its application areas are extended.

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.