NONLINEAR TRANSIENT RESPONSE OF FUNCTIONALLY GRADED SHALLOW SPHERICAL SHELLS SUBJECTED TO MECHANICAL LOAD AND UNSTEADY TEMPERATURE FIELD

This paper is concerned with the transient deformation of functionally graded （FG） shallow spherical shells subjected to time-dependent thermomechanical load. Based on Timoshenko- Mindlin hypothesis and yon Karman nonlinear theory, a set of nonlinear governing equations of motion for FG shallow spherical shells in regard to transverse shear deformation and all the inertia terms are established using Hamilton＇s principle. The collocation point method and Newmark- beta scheme in conjunction with the finite difference method are adopted to solve the governing equations of motion and the unsteady heat conduction equation numerically. In the numerical examples, the transient deflection and stresses of FG shallow spherical shells with various material properties under different loading conditions are presented.     

Thermal and mechanical analysis of material response to non-steady ramp and steady shock wave loading

Ramp wave experiments on the Sandia Z accelerator provide a new approach to study the rapid compression response of materials at pressures, temperatures and stress or strain rates not attainable in conventional shock experiments. Due to its shockless nature, the ramp wave experiment is often termed as an isentropic (or quasi-isentropic) compression experiment (ICE). However, in reality there is always some entropy produced when materials are subjected to large amplitude compression even under shockless loading. The entropy production mechanisms that cause deformation to deviate from the isentropic process can be attributed to mechanical and thermal dissipations. The former is due to inelasticity associated with various deformation mechanisms and the rate effect that is inherent in all the deformation processes and the latter is due to irreversible heat conduction. The main purpose of the current study is to gain insights into the effects of ramp and shock loading on the entropy production and thermomechanical responses of materials. Another purpose is to investigate the role of heat conduction in the material response to both the non-steady ramp wave and steady shock.Numerical simulations are used to address the aforementioned research objectives. The thermomechanical response associated with a steady shock wave is investigated first by solving a set of nonlinear ordinary differential equations. Using the steady wave solutions as the reference, the material responses under non-steady ramp waves are then studied with numerical wave propagation simulation. It is demonstrated that the material response to ramp and shock loading is essentially a manifestation of the interaction between the time scale associated with the loading and the intrinsic time scales associated with mechanical deformation and heat transfer. At lower loading rates as encountered in ramp loading, the loading path is closer to an isentrope and results in lower entropy production. The reasonable ramp rate to obtain a quasi-isentropic state depends on the intrinsic time scales of the dissipation mechanisms which are strongly material dependent. Thus shockless loading does not necessarily produce an isentropic response. Between two equilibrium states, heat conduction was shown to have significant effect on the temperature history but it contributes little to the overall temperature change if the specific heat remains constant. It also affects the history of entropy, but only the irreversible part of heat conduction contributes to the net entropy change. The various types of thermomechanical responses of materials would manifest themselves more significantly in terms of the thermal history than the mechanical history. Thus temperature measurement appears to be an important experimental tool in distinguishing the various mechanisms for the thermomechancial responses of the materials.     

Kármán-type equations for a higher-order shear deformation plate theory and its use in the thermal postbuckling analysis

Kármán-type nonlinear large deflection equations are derived occnrding to the Reddy’s higher-order shear deformation plate theory and used in the thermal postbuckling analysis The effects of initial geometric imperfections of the plate areincluded in the present study which also includes th thermal effects.Simply supported,symmetric cross-ply laminated plates subjected to uniform or nomuniform parabolictemperature distribution are considered. The analysis uses a mixed GalerkinGolerkinperlurbation technique to determine thermal buckling louds and postbucklingequilibrium paths.The effects played by transverse shear deformation plate aspeclraio, total number of plies thermal load ratio and initial geometric imperfections arealso studied.     

纵横载荷作用下功能梯度梁的弯曲和过屈曲

基于一阶非线性梁理论和物理中面概念,导出了纵横向载荷作用下功能梯度材料(FGM)梁非线性弯曲和过屈曲问题的控制方程,并获得了该问题的精确解;据此解研究了梯度材料性质、外载荷、横向剪切变形以及边界条件等因素对功能梯度材料梁非线性力学行为的影响,分析中假设功能梯度材料性质只沿梁厚度方向,并按成分含量的幂指数函数形式变化。结果表明:纵横载荷共同作用下,功能梯度梁的弯曲构形将有无限多个;随着梯度指数的增大,梁的变形减小,临界载荷升高;随着长高比的增大,横向剪切变形的影响减小。     

Effects of deformed volume, volume fraction and particle size on the deformation behaviour of W/Cu composites

Tungsten/copper (W/Cu) particle reinforced composites were used to investigate the scaling effects on the deformation and fracture behaviour. The effects of the volume fraction and the particle size of the reinforcement (tungsten particles) were studied. W/Cu-80/20, 70/30 and 60/40 wt.% each with tungsten particle size of 10 μm and 30 μm were tested under compression and shear loading. Cylindrical compression specimens with different volumes (D_S = H) were investigated with strain rates between 0.001 s⁻¹ and about 5750 s⁻¹ at temperatures from 20 °C to 800 °C. Axis-symmetric hat-shaped shear specimens with different shear zone widths were examined at different strain rates as well. A clear dependence of the flow stress on the deformed volume and the particle size was found under compression and shear loading. Metallographic investigation was carried out to show a relation between the deformation of the tungsten particles and the global deformation of the specimens. The size of the deformed zone under either compression or shear loading has shown a clear size effect on the fracture of the hat-shaped specimens.The quasi-static flow curves were described with the material law from Swift. The parameters of the material law were presented as a function of the temperature and the specimen size. The mechanical behaviour of the composite materials were numerically computed for an idealized axis-symmetric hat-shaped specimen to verify the determined material law.     

Nonlinear deformation and buckling of elastic inhomogeneous shells under thermomechanical loads

The paper outlines the fundamentals of the method of solving static problems of geometrically nonlinear deformation, buckling, and postbuckling behavior of thin thermoelastic inhomogeneous shells with complex-shaped midsurface, geometrical features throughout the thickness, or multilayer structure under complex thermomechanical loading. The method is based on the geometrically nonlinear equations of three-dimensional thermoelasticity and the moment finite-element scheme. The method is justified numerically. Results of practical importance are obtained in analyzing poorely studied classes of inhomogeneous shells. These results provide an insight into the nonlinear deformation and buckling of shells under various combinations of thermomechanical loads     

A new non-linear higher-order shear deformation theory for large-amplitude vibrations of laminated doubly curved shells

A consistent higher-order shear deformation non-linear theory is developed for shells of generic shape, taking geometric imperfections into account. The geometrically non-linear strain-displacement relationships are derived retaining full non-linear terms in the in-plane displacements; they are presented in curvilinear coordinates in a formulation ready to be implemented. Then, large-amplitude forced vibrations of a simply supported, laminated circular cylindrical shell are studied (i) by using the developed theory, and (ii) keeping only non-linear terms of the von Kármán type. Results show that inaccurate results are obtained by keeping only non-linear terms of the von Kármán type for vibration amplitudes of about two times the shell thickness for the studied case.     

Mechanical response and texture evolution of AZ31 alloy at large strains for different strain rates and temperatures

In order to study the behavior of material under finite deformation at various strain rates, the responses of AZ31 Mg sheet are measured under uniaxial (tension and compression) and multiaxial (simple shear) loadings along rolling direction (RD), 45° to rolling direction (DD), 90° to rolling direction (TD), and normal to the sheet (ND) to large strains. The material exhibits positive strain rate sensitivity (SRS) at room and elevated temperatures; the SRS is more pronounced at high temperatures and lower strain rates. The r-value of the material under tensile loading at room temperatures is higher in TD at lower strain rate. Texture measurements on several failed specimens are reported under tension and simple shear after finite plastic deformation of about 20% equivalent strain. The as-received material exhibits a strong fiber with equal fractions of grains having the c-axis slightly tilted away from the sheet normal towards both +RD and −RD. Pole figures obtained after tensile loading along the rolling direction (RD) show that the texture of the material strengthens even at low strains, with c-axis perpendicular to the sheet plane and prism planes lining up in a majority of grains. However, the tensile loading axis along TD does not lead to similar texture strengthening; the c-axis distribution appears to be virtually unchanged from the virgin state. The pole figures obtained after in-plane compression along RD brings the c-axes of the grains parallel to the loading direction. The pole figures after simple shear loading show that the c-axis rotates to lie on the sheet plane consistent with a compression axis 45° away on the sheet plane.     

An <Emphasis Type="Italic">In Situ</Emphasis> Ultrasonic Technique for Simultaneous Measurement of Longitudinal and Shear Wave Speeds in Solids

Ultrasonic wave speed measurements are widely used to infer the elastic properties of solids. In the standard method, longitudinal and shear transducers are used separately to measure the corresponding wave speeds in a material. A new experimental method is introduced for simultaneously measuring the longitudinal and the shear wave speeds using a single set of longitudinal or shear transducers. The method can also be used to measure the wave speeds in situ during deformation by placing the transducers along the loading axis. The transducers are housed in a specially designed fixture such that they are not subjected to loading. The technique is demonstrated by performing uniaxial compression experiments on fully dense isotropic solids (where the wave speeds are not expected to change during deformation) and in polymeric foams (where the wave speeds are affected by damage).     

Shear failure of tungsten alloy at mesoscale modeled by Movable cellular Automata

Shear failure and deformation localization are studied for tungsten alloy undergoing impact failure and “self-sharpening” in ordnance. Movable cellular Automata (MCA) method is employed to simulate the interior-structure of tungsten alloy such that shear deformation failure at mesoscale level under shear loading can be analyzed. The results show that the velocity of every automaton would vary continuously with deformation. A large velocity vortex sweeping throughout the specimen occurs at the initial loading stage. This corresponds to uniform deformation. Subsequently, the deformation adopts a pattern with more localized vortices. The vortex-morphology would change continually to sustain the global deformation field. The contact between the local vortices induces stress concentration. This means that the deformation of certain automaton pairs would not accommodate the equilibrium of vortex rotation. The different mechanisms of micro-damage mechanisms along the interface between the W grain and matrix are developed quickly throughout the whole specimen. Such shear instability leads to failure. It is also shown that the W grain size effect is sufficient to cause premature shear deformation instability of the tungsten alloy, even at low shear stress.     

Loading direction-dependent shear behavior at different temperatures of single-layer chiral graphene sheets

The loading direction-dependent shear behavior of single-layer chiral graphene sheets at different temperatures is studied by molecular dynamics(MD) simulations.Our results show that the shear properties(such as shear stress–strain curves, buckling strains, and failure strains) of chiral graphene sheets strongly depend on the loading direction due to the structural asymmetry. The maximum values of both the critical buckling shear strain and the failure strain under positive shear deformation can be around 1.4 times higher than those under negative shear deformation. For a given chiral graphene sheet, both its failure strain and failure stress decrease with increasing temperature. In particular, the amplitude to wavelength ratio of wrinkles for different chiral graphene sheets under shear deformation using present MD simulations agrees well with that from the existing theory. These findings provide physical insights into the origins of the loading direction-dependent shear behavior of chiral graphene sheets and their potential applications in nanodevices.     

Inelastic behavior of an AS4/PEEK composite under combined transverse compression and shear. Part I: experiments

The nonlinear behavior in shear and transverse compression of unidirectional AS4/PEEK and their interaction are investigated experimentally. The composite is rate dependent even at room temperature and its rate exponent is similar to that of neat PEEK. The material is tested under pure shear, pure compression and under biaxial loading histories. The biaxial tests are performed in a custom facility on thin strips of the material. The facility allows freedom to choose the loading path in the biaxial stress and strain spaces of interest. Tests are performed for three biaxial loading paths. In the first, the specimen is sheared then compressed while the shear stress is held constant; in the second, the specimen is compressed then sheared while the compressive stress is held constant; and in the third, the specimen is loaded simultaneously by proportional amounts of compression and shear. It was found that the induced deformation is influenced significantly by the loading history followed. Also, initial loading in shear or compression has only a modest effect on subsequent loading of the other type. An unorthodox yielding behavior for the composite results from this lack of interaction. Finally, the stresses at failure are found to trace an elliptical path in the shear–transverse compression plane, but the failure stress state is not significantly affected by the loading path followed.     

不同加载状态下TA2钛合金绝热剪切破坏响应特性

一般认为绝热剪切现象在宏观上表现为材料动态本构失稳,即热软化大于应变硬化.本文采用帽型受迫剪切试样研究TA2钛合金的动态力学特性和本构失稳过程.首先对剪切区加载应力状态进行理论和数值分析,通过合理设计帽型试样,剪切区变形可近似按剪切状态处理;结合二维数字图像相关法(two-dimensional digital image correlation,DIC-2D)直接测试试样剪切区应变演化,给出帽型受迫剪切实验的等效应力-应变响应曲线.进一步,利用Hopkinson压杆对TA2钛合金开展动态压缩及帽型剪切对比试验研究,比较压缩、剪切试验得到的等效应力-应变曲线,采用"冻结"试样方法分析试样中绝热剪切局域化演化过程,探讨不同加载状态下TA2钛合金的绝热剪切破坏现象及其动态力学响应特性.实验结果表明,在塑性变形初始阶段,动态压缩及剪切加载下的等效应力-应变曲线符合较好,但随塑性损伤发展及绝热剪切带形成,两者出现分离,表明损伤及绝热剪切演化过程与应力状态相关.剪切试样实验得到的本构"软化"特性能够反映绝热剪切带起始、破坏演化过程的力学响应特性,而在动态压缩实验中,即使试样中已出现双锥形的绝热剪切带及局部裂纹分布,其表观等效应力-应变曲线并不出现软化特征,动态压缩实验无法得到关于绝热剪切起始、发展以及破坏的本构软化响应特性.     

Lüders bands propagation of 1045 steel under multiaxial stress state

Inhomogeneous plastic deformation of 1045 steel under monotonic loading was experimentally studied. Thin-walled tubular specimens were used in the experiments and custom-made small strain gages were bonded on the specimen surface to characterize the local deformation. Experiments were conducted under tension, torsion, and combined tension–torsion. During the propagation of Lüders bands, the local deformation experienced two-stage deformation: an abrupt plastic deformation stage followed by a slower deformation process. In some area of the gage section of the specimen, a small amount of initial plastic deformation occurred before the Lüders front reached. During the propagation of Lüders bands, multiple Lüders fronts can be formed. Under tension, torsion, and combined tension–torsion with a constant axial load, the Lüders front was approximately parallel to the material plane of maximum shear stress. When the combined axial-torsion followed a proportional fashion, the stress–extensometer strain responses were dependent on the axial/torsional loading ratio, and the Lüders fronts were oriented differently and propagated along the specimen axis at a different velocity. The local strain was inhomogeneous even at the work-hardening stage. The relationships between the equivalent stress and the equivalent plastic strain were found to be practically identical for all the loading cases studied.     

功能梯度中厚圆/环板轴对称弯曲问题的解析解

基于一阶剪切变形板理论，导出了热/机载荷作用下，位移形式的功能梯度 中厚圆/环板轴对称弯曲问题的控制方程，获得了问题的位移和内力的一般解析解. 作为特 例，分别研究了边缘径向固定和可动的夹紧和简支的4种实心功能梯度圆板，给出了它们的 解，并分析了热/机载荷作用下解的形态，讨论了横向剪切变形、材料梯度常数和边界条件， 对板的轴对称弯曲行为的影响.     

Multiaxial and non-proportional loading responses,anisotropy and modeling of Ti–6Al–4V titanium alloy over wide ranges of strain rates and temperatures

Results from a series of multiaxial loading experiments on the Ti–6Al–4V titanium alloy are presented. Different loading conditions are applied in order to get the comprehensive response of the alloy. The strain rates are varied from the quasi-static to dynamic regimes and the corresponding material responses are obtained. The specimen is deformed to large strains in order to study the material behavior under finite deformation at various strain rates. Torsional Kolsky bar is used to achieve shear strain rates up to 1000 s⁻¹. Experiments are performed under non-proportional loading conditions as well as dynamic torsion followed by dynamic compression at various temperatures. The non-proportional loading experiments comprise of an initial uniaxial loading to a certain level of strain followed by biaxial loading, using a channel-type die at various rates of loadings. All the non-proportional experiments are carried out at room temperature. Experiments are also performed to investigate the anisotropic behavior of the alloy. An orthotropic yield criterion [proposed by Cazacu, O., Plunkett, B., Barlat, F., 2005. Orthotropic yield criterion for hexagonal closed packed metals. International Journal of Plasticity 22, 1171–1194.] for anisotropic hexagonal closed packed materials with strength differential is used to generate the yield surface. Based on the definition of the effective stress of this yield criterion, the observed material response for the different loading conditions under large deformation is modeled using the Khan–Huang–Liang (KHL) equation assuming isotropic hardening. The model constants used in the present study, were pre-determined from the extensive uniaxial experiments presented in the earlier paper [Khan, A.S., Suh, Y.S., Kazmi R., 2004. Quasi-static and dynamic loading responses and constitutive modeling of titanium alloys. International Journal of Plasticity 20, 2233–2248]. The model predictions are found to be extremely close to the observed material response.     

Compression response of a thermoplastic elastomer gel tissue surrogate over a range of strain-rates

The mechanical response of a thermoplastic elastomer synthetic ballistic gel is studied over a range of strain-rates. Experiments were conducted at room temperature under uniaxial stress compression at rates ranging from 0.001 to 500/s. Low-rate experiments (<1/s) were conducted with a servo-hydraulic load frame. High-rate experiments (>100/s) were conducted with a polymeric Kolsky bar, along with several modifications to improve data quality. These modifications include the use of a commercial force transducer, a normal displacement interferometer, and a line laser extensometer. Because of the low shear strength and comparatively high compressibility of these materials, inertial effects are very pronounced. Specimen size is varied in an effort to study inertial effects at various loading rates. High speed photography is also used to demonstrate the presence of non-uniform deformation, due to both inertia and friction between the specimen and the loading surfaces. Finally, numerical simulation is used to verify trends observed in the experiments and further validate the data. It is concluded that this material is rate sensitive, with an almost three-fold increase in stiffness over the range of strain-rates studied.     

Analysis of anisotropy and plastic spin effects on localization phenomena

Summary The main objective of the paper is the investigation of the influence of the anisotrophy and plastic spin effects on criteria for adiabatic shear band localization of plastic deformation. A theory of thermoplasticity is formulated within a framework of the rate-type covariance material structure with a finite set of internal state variables. The theory takes into consideration such effects as plastic non-normality, plastic-induced anisotropy (kinematic hardening), micro-damage mechanism, thermomechanical coupling and plastic spin. The next objective of the paper is to focus attention on cooperative phenomena in presence of the plastic spin, and the discussion on the influence of synergetic effects on localization criteria. A particular constitutive law for the plastic spin is assumed. The necessary condition for a localized plastic deformation region to be formed is obtained. This condition is accomplished by the assumption that some eigenvalues of the instantaneous adiabatic acoustic tensor vanish. A procedure has been developed which allows us to discuss two separate groups of effects on the localization phenomenon along a shear band. Plastic spin, spatial covariance and kinematic hardening effects are investigated at an isothermal process in an undamaged solid. In the second case, an adiabatic process in a damaged solid is discussed when the spatial covariance terms and the plastic spin are neglected. Here the thermomechanical coupling, micro-damage mechanism and kinematic hardening effects are examined. For both cases, the criteria for adiabatic shear band localization are obtained in an exact analytical form. Particular attention is focused on the analysis of the following effects: (i) plastic non-normality; (ii) plastic spin; (iii) covariant terms; (iv) plastic strain-induced anisotropy; (v) micro-damage mechanism; (vi) thermomechanical couplings. Cooperative phenomena are considered, and synergetic effects are investigated. A discussion of the influence of the plastic spin, kinematic hardening and covariant terms on the shear band localization conditions is presented. A numerical estimation of the effects discussed is given. Received 24 April 1997; accepted for publication 23 December 1997     

The growth of shear bands in composite microstructures

Shear band formation in materials with inhomogeneous and composite microstructures is influenced by factors that usually do not come into play in monolithic materials. Experiments and calculations have shown that inhomogeneities in material properties enhance the localization of deformation. This investigation concerns the propagation of shear bands in a two-phase tungsten composite under the conditions of nominally pure shear deformation. Finite element calculations are carried out to delineate the effects of different grain–matrix morphologies. In the numerical models, the initiation of shear bands is triggered by a notch, simulating the effect of defects such as microcracks and microvoids in materials. Calculations demonstrate that phase morphology, particle size and the relative location of initiation site have significant influences on the development of localized deformation. The work and energy evolutions are tracked for each constituent phase in the microstructures. In addition, the exchange of thermal energy through heat flow between the phases is analyzed. The results show that a strong correlation exists between the course of shear band propagation and the thermomechanical coupling between microscopic phases.