Strain gradient plasticity solution for an internally pressurized thick-walled cylinder of an elastic linear-hardening material

An elastic-plastic solution is presented for an internally pressurized thick-walled plane strain cylinder of an elastic linear-hardening plastic material. The solution is derived in a closed form using a strain gradient plasticity theory. The inner radius of the cylinder enters the solution not only in non-dimensional forms but also with its own dimensional identity, which differs from that in classical plasticity based solutions and makes it possible to capture the size effect at the micron scale. The classical plasticity solution of the same problem is recovered as a special case of the current solution. To further illustrate the newly derived solution, formulas and numerical results for the plastic limit pressure are provided. These results reveal that the load-carrying capacity of the cylinder increases with decreasing inner radius at the micron scale. It is also seen that the macroscopic behavior of the pressurized cylinder can be well described by using classical plasticity based solutions.     

基于应变梯度理论的粘塑性厚壁圆筒和球壳极限内压分析

基于应变梯度塑性理论,分析了内压作用下厚壁圆筒和球壳的塑性极限荷载．结果表明:圆筒内径在微米量级时,存在尺度效应现象,内径减小,其尺度效应增强;变形越大,影响越大;应变速率敏感指数越大,尺度效应越明显．经典塑性理论结果是当前解的特例．     

Strain gradient solutions of half-space and half-plane contact problems

General solutions for the problems of an elastic half-space and an elastic half-plane, respectively, subjected to a symmetrically distributed normal force of arbitrary profile are analytically derived using a simplified strain gradient elasticity theory (SSGET) that contains one material length scale parameter. Mindlin’s potential function method and Fourier transforms are employed in the formulation, and the half-space and half-plane contact problems are solved in a unified manner. The specific solutions for the problems of a half-space/plane subjected to a concentrated normal force or a uniformly distributed normal force are obtained by directly applying the general solutions, which recover the existing classical elasticity-based solutions of the Flamant and Boussinesq problems as special cases. In addition, the indentation problems of an elastic half-space indented by a flat-ended cylindrical punch, a spherical punch, and a conical punch, respectively, are solved using the general solutions, leading to hardness formulas that are indentation size- and material microstructure-dependent. Numerical results reveal that the displacement and stress fields in a half-space/plane given by the current SSGET-based solutions are smoother than those predicted by the classical elasticity-based solutions and do not exhibit the discontinuity and/or singularity displayed by the latter. Also, the indentation hardness values based on the newly obtained half-space solution are found to increase with decreasing indentation radius and increasing material length scale parameter, thereby explaining the microstructure-dependent indentation size effect.     

理想塑性固体中逐步扩展裂纹的弹塑性场

本文从三维的塑性流动理论出发,导出了关于理想塑性固体平面应变问题的基本方程。利用这些方程,分析了不可压缩理想塑性固体的逐步扩展裂纹顶端的弹塑性场。得到了关于应力和速度的一阶渐近场。分析了弹性卸载区的演变过程和中心扇形区的发展过程。预示了出现二次塑性区的可能性。最后给出了关于应力场二阶渐近分析。     

FEM solutions for plane stress mode-I and mode-II cracks in strain gradient plasticity

The strain gradient plasticity theory is used to investigate the crack-tip field in a power law hardening material. Numerical solutions are presented for plane-stress mode I and mode II cracks under small scale yielding conditions. A comparison is made with the existing asymptotic fields. It is found that the size of the dominance zone for the near-tip asymptotic field, recently obtained by Chen et al., is on the order 5% of the intrinsic material lengthI. Remote from the dominance zone, the computed stress field tends to be the classical HRR field. Within the plastic zone only force-stress dominated solution is found for either mode I or mode II crack.     

Plastic flow in a conical channel: Qualitative features of the solutions under different yield conditions

The problem of the convergence of the solutions of problems of plasticity theory, with a yield condition which depends on the hydrostatic stress, to solutions based on classical plasticity theory with von Mises or Tresea conditions is considered, with a particular choice of the parameters of the material model. For the case of axisymmetric flow of material in a channel with converging and diverging walls, solutions according to two plasticity theories with a yield condition which depends on the hydrostatic stress are compared with the classical solution. It is shown that only the solution using Spencer's model shows all the main features of the classical solution. As the internal criterion of the choice of the preferred plasticity theory when examining a special class of problems, it is suggested that the criterion of the convergence of the solutions to the solutions of classical plasticity theory should be used.     

Quasistatic evolution for Cam-Clay plasticity: a weak formulation via viscoplastic regularization and time rescaling

Cam-Clay nonassociative plasticity exhibits both hardening and softening behaviour, depending on the loading. For many initial data the classical formulation of the quasistatic evolution problem has no smooth solution. We propose here a notion of generalized solution, based on a viscoplastic approximation. To study the limit of the viscoplastic evolutions we rescale time, in such a way that the plastic strain is uniformly Lipschitz with respect to the rescaled time. The limit of these rescaled solutions, as the viscosity parameter tends to zero, is characterized through an energy-dissipation balance, that can be written in a natural way using the rescaled time. As shown in Dal Maso and DeSimone (Math Models Methods Appl Sci 19:1?C69, 2009) and Dal Maso and Solombrino (Netw Heterog Media 5:97?C132, 2010), the proposed solution may be discontinuous with respect to the original time. Our formulation allows us to compute the amount of viscous dissipation occurring instantaneously at each discontinuity time.     

权函数法研究高速旋转厚壁筒的应力强度因子

采用厚壁筒在内压作用下的应力强度因子作为参考载荷的应力强度因子,通过权函数方法推导出了内壁带二维径向边裂纹的旋转厚壁筒的应力强度因子公式．这些公式可用于计算旋转厚壁筒在不同裂纹深度、转速、材料和尺寸情况下的应力强度因子．算例表明该文的公式具有良好的精度．同时还研究了旋转圆筒应力强度因子随裂纹深度和内外径比之间的变化规律,方便了工程应用．     

Axisymmetric flow over a backward-facing step in an annular pipe

The incompressible flow of a Newtonian fluid over a backward-facing step is investigated numerically. The geometry is an annular pipe in which the radius of the inner cylinder decreases suddenly. Keeping the radial expansion ratio fixed axisymmetric flows are computed for outlet radius ratios from 0.1 to 1 (ratio of the inner to the outer outlet radius). The Reynolds number at which the flow separates from the outer cylinder decreases as the outlet radius ratio decreases for constant inlet geometry. The growth with Reynolds number of the recirculation zone on the inner outlet cylinder just behind the step is strongly reduced when the recirculation zone on the outer cylinder is established. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Incompatibility Induced Nonlocal Crystal Plasticity

Motivated by the fact that locally inhomogeneous elastic or plastic deformations may result in incompatibilities of the fictitious intermediate configuration a strain gradient crystal plasticity model is developed. Thereby incompatibilities can be accounted for and scale dependent material behavior, as also observed experimentally, is predictable. A nonlocal extension of existing local formulations is proposed which does not require additional boundary conditions and thus maintains the classical BVP structure. On the numerical side key developments are an extended FE-formulation for rate-(in)dependent strain gradient plasticity and a local FE-formulation which bases the gradient computation on an operator split combined with a smoothing algorithm. Comparative numerical studies for classical examples proove the superior efficiency of the second approach. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A semi-analytical solution for shallow tunnels with radius-iterative-approach in semi-infinite space

This study presents a semi-analytical elastic-plastic solution for a shallow tunnel subjected to ground loss in the strain-softening surrounding rock. The most important contribution is the radius-iterative-approach in which the initial plastic radius is first determined by the strain continuity boundary condition on the elastic-plastic interface and then corrected to the precise one. The corrected approach follows three steps: (1) Applying the radius increment technique to semi-infinite space (2) Carrying out the plastic radius correction by using iteration method from the elastic-plastic interface to the tunnel wall. (3) If the calculated convergence value is equal to the convergence value on the tunnel wall, the accurate determination of the plastic region, stresses, and displacements, of the whole half plane, can be derived consequently. All the results compare favorably with numerical simulation results. The study completes the theoretical framework for addressing the fundamental problem of shallow tunnels excavated in the semi-infinite space and also provides a useful theoretical tool for potential application on the tunnel and underground engineering problems.     

Viscoelastic flow in a curved channel: A similarity solution for the Oldroyd-B fluid

A similarity solution is used to analyse the flow of the Oldroyd fluid B, which includes the Newtonian and Maxwell fluids, in a curved channel modelled by the narrow annular region between two circular concentric cylinders of large radius. The solution is exact, including inertial forces. It is found that the non-Netonian kinematics are very similar to the Newtonian ones, although some stress components can become very large. At high Reynolds number a boundary layer is developed at the inner cylinder. The structure of this boundary layer is asymptotically analysed for the Newtonian fluid. Non-Newtonian stress boundary layers are also developed at the inner cylinder at large Reynolds numbers.     

Effects of slip on the non-linear flows of a third grade fluid

The present analysis considers the non-linear problems of steady flow of a third grade fluid between the concentric cylinders. A complete analysis of mathematical modeling is made when no-slip condition is no longer valid. Exact analytic solutions of the following two non-linear problems are derived: (i) when inner cylinder moves and outer cylinder remains stationary and (ii) for inner cylinder at rest and outer cylinder in motion. Graphical results are presented to illustrate the analytic solutions. The corresponding results of no-slip condition are deduced as the limiting cases when the slip parameter is equal to zero.     

Creeping flow of a sphere nearby a cylinder

We present a three-dimensional solution of a sphere nearby an infinite cylinder at low Reynolds number. We utilize the Lamb’s general solution based on spherical harmonics and develop a framework based on cylindrical harmonics to solve the flow field around the sphere and outside the cylinder, respectively. The solution is solved semi-analytically by considering geometrical parameters, including sphere radius, sphere velocity, separation distance and cylinder radius. The drag force coefficients of the sphere which are dependent on the distance between the cylinder surface and the sphere, as well as the velocity contours in the vicinity of the sphere, are analyzed. We also provide an analytical formula to calculate the drag force. The analytical formula has good quantitative agreement with the semi-analytical solution when the radius of the cylinder is smaller than the sphere. Such analysis can give insights into the details of the complex interaction between the sphere and cylinder.     

Thermal analysis of moving induction heating of a hollow cylinder with subsequent spray cooling: Effect of velocity,initial position of coil,and geometry

In this paper, temperature analysis of the complete process of moving induction heat treatment is performed using numerical methods. A non-linear and transient magneto-thermal coupled problem with a moving coil which is considered as moving heat source, is investigated by an efficient finite-element procedure. A vertical hollow circular cylinder is heated by the moving coil at a given velocity along it, and the heated parts then quenched by a moving water–air spray. The effects of natural convection with air on the both inner and outer surfaces of cylinder, and also radiation of outer surface of cylinder with ambient are taken into account. For quenching of work-piece, a specific kind of atomized spray cooling which utilizes a mixture of water and air with different mass fractions is used. This procedure includes moving boundary conditions, temperature-dependent properties, and change in magnetic permeability of specified alloy at the Curie temperature. Obtained numerical results have been verified by comparison with analytical solutions using Green’s function methods. Also, the effect of velocity, initial position of inductor and inner to outer radius ratio on temperature distribution are investigated.     

The limit equilibrium of an anisotropic medium under the general plasticity condition

A theory of the limit equilibrium of an anisotropic medium under the general plasticity condition in the plane strain state is developed. The proposed yield criterion (the limit equilibrium condition) is obtained by combining the von Mises–Hill yield criterion of an ideally plastic anisotropic material and Prandtl's limit equilibrium condition for a medium under the general plasticity law. It is shown that the problem is statically determinate, i.e., if the boundary conditions are specified in stresses, the stress state in plastic region can only be obtained using equilibrium equations. It is established that the equations describing the stress state are hyperbolic and have two families of characteristic curves that intersect at variable angles. In deriving the equations describing the velocity field, the material is assumed to be rigid plastic, and the associated law of flow is applied. It is shown that the equations for the velocities are also hyperbolic, and their characteristic curves are identical with those of the equations for stresses. However, the directions of the principal values of the stress and strain rate tensors are different due to the anisotropy of the material. The characteristic directions differ from the isotropic case in that the normal and tangential components of the stress tensor do not satisfy the limit conditions. It is established that the equations obtained allow of partial solutions, and in this case, at least one family of characteristic curves consists of straight lines. The conditions along the lines of discontinuity of the velocity are investigated, and it is shown that, as in the isotropic case, these are characteristic curves of the system of governing equations. In the anisotropic formulation, the well-known Rankine problem of the limit state of a ponderable layer is solved. From an analysis of the velocity field it is shown that plastic flow of the entire layer is possible only for a slope angle equal to the angle of internal friction. For slope angles less than the angle of internal friction, the solutions obtained are solutions of problems of the pressure of the medium on the retaining walls. The change in this pressure as a function of the parameters of anisotropy is investigated, and turns out to be significant.     

Analytical solution of the quasi-static thermoelasticity problem in a pressurized thick-walled cylinder subjected to transient thermal loading

Thermoelasticity problem in a thick-walled cylinder is solved analytically using the finite Hankel transform. Time-dependent thermal boundary conditions are assumed to act on the inner surface of the cylinder. For the mechanical boundary conditions two different cases are assumed: Traction–displacement problem (traction is prescribed on the inner surface and the fixed displacement boundary condition on the outer one) and Traction–Traction problem (tractions are prescribed on both the inner and outer surfaces of the hollow cylinder). The quasi-static solution of the thermoelasticity problem is derived analytically, i.e., the transient thermal response of the cylinder is derived and then, quasi-static structural problem is solved and closed form relations are extracted for the thermal stresses in the two problems. The results show to be in accordance with that cited in the literature in the special cases.     

Dynamics of Vortex Rossby Waves in Tropical Cyclones,Part 1: Linear Time‐Dependent Evolution on an f‐Plane

Analyses of observational data on hurricanes in the tropical atmosphere indicate the existence of spiral rainbands which propagate outward from the eye and affect the structure and intensity of the hurricane. These disturbances may be described as vortex Rossby waves. This paper describes the evolution of barotropic vortex Rossby waves in a cyclonic vortex in a two‐dimensional configuration where the variation of the Coriolis force with latitude is ignored. The waves are forced by a constant‐amplitude boundary condition at a fixed radius from the center of the vortex and propagate outward. The mean flow angular velocity profile is taken to be a quadratic function of the radial distance from the center of the vortex and there is a critical radius at which it is equal to the phase speed of the waves. For the case of waves with steady amplitude, an exact solution is derived for the steady linearized equations in terms of hypergeometric functions; this solution is valid in the outer region away from the critical radius. For the case of waves with time‐dependent amplitude, asymptotic solutions of the linearized equations, valid for late time, are obtained in the outer and inner regions. It is found that there are strong qualitative similarities between the conclusions on the evolution of the vortex waves in this configuration and those obtained in the case of Rossby waves in a rectangular configuration where the latitudinal gradient of the Coriolis parameter is taken into account. In particular, the amplitude of the steady‐state outer solution is greatly attenuated and there is a phase change of across the critical radius, and in the linear time‐dependent configuration, the outer solution approaches a steady state in the limit of infinite time, while the amplitude of the inner solution grows on a logarithmic time scale and the width of the critical layer approaches zero.     

Towards variational constitutive updates for finite strain plasticity models based on non-associative evolution equations

In this contribution, first steps towards variational constitutive updates for finite strain plasticity theory based on non-associative evolution equations are presented. These schemes allow to compute the unknown state variables such as the plastic part of the deformation gradient, together with the deformation mapping, by means of a fully variational minimization principle. Therefore, standard optimization algorithms can be applied to the numerical implementation leading to a very robust and efficient numerical implementation. Particularly, for highly non-linear, singular or nearly ill-posed physical models like that corresponding to crystal plasticity showing a large number of possible active slip planes, this is a significant advantage compared to standard constitutive updates such as the by now classical return-mapping algorithm. While variational constitutive updates have been successfully derived for associative plasticity models, their extension to more complex constitutive laws, particularly to those featuring non-associative evolution equations, is highly challenging. In the present contribution, a certain class of non-associative finite strain plasticity models is discussed and recast into a variationally consistent format. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)