Spherically isotropic,elastic spheres subject to diametral point load strength test

This paper presents an analytic solution for the stress concentrations within a spherically isotropic, elastic sphere of radius R subject to diametral point load strength test. The method of solution uses the displacement potential approach together with the Fourier–Legendre expansion for the boundary loads. For the case of isotropic sphere, our solution reduces to the solution by Hiramatsu and Oka, 1966 and agrees well with the published experimental observations by Frocht and Guernsey (1953) . A zone of higher tensile stress concentration is developed near the point loads, and the difference between this maximum tensile stress and the uniform tensile stress in the central part of the sphere increases with E/E′ (where E and E′ are the Youngs moduli governing axial deformations along directions parallel and normal to the planes of isotropy, respectively) , G′/G (where G and G′ are the moduli governing shear deformations in the planes of isotropy and the planes parallel to the radial direction) , and ν̄/ν′ (where ν̄ and ν′ are the Poissons ratios characterizing transverse reduction in the planes of isotropy under tension in the same plane and under radial tension, respectively) . This stress difference, in general, decreases with the size of loading area and the Poissons ratio.     

A model of dilute polymer solutions with hydrodynamic interaction and finite extensibility. II. Shear flows

The results from numerical calculations for steady shear, the start-up and cessation of steady shear, and the stress relaxation after a step shear stain are discussed in detail for a bead-spring chain model with consistently averaged hydrodynamic interaction between the beads and consistently averaged finite extensible springs. Calculations are made for a large range of spring stiffnesses 10 ⩽ b ⩽ ∞ and a hydrodynamic interaction strength h^* = 0.15 a value which has been estimated from experimental results. This model is found to satisfy the Hassager-Bird and the Lodge-Meissner relations.     

Turbulent Drag Reduction by a Near Wall Surface Tension Active Interface

In this work we study the turbulence modulation in a viscosity-stratified two-phase flow using Direct Numerical Simulation (DNS) of turbulence and the Phase Field Method (PFM) to simulate the interfacial phenomena. Specifically we consider the case of two immiscible fluid layers driven in a closed rectangular channel by an imposed mean pressure gradient. The present problem, which may mimic the behaviour of an oil flowing under a thin layer of different oil, thickness ratio h₂/h₁ =?9, is described by three main flow parameters: the shear Reynolds number Re _τ (which quantifies the importance of inertia compared to viscous effects), the Weber number We (which quantifies surface tension effects) and the viscosity ratio λ = ν₁/ν₂ between the two fluids. For this first study, the density ratio of the two fluid layers is the same (ρ₂ = ρ₁), we keep Re _τ and We constant, but we consider three different values for the viscosity ratio: λ =?1, λ =?0.875 and λ =?0.75. Compared to a single phase flow at the same shear Reynolds number (Re _τ =?100), in the two phase flow case we observe a decrease of the wall-shear stress and a strong turbulence modulation in particular in the proximity of the interface. Interestingly, we observe that the modulation of turbulence by the liquid-liquid interface extends up to the top wall (i.e. the closest to the interface) and produces local shear stress inversions and flow recirculation regions. The observed results depend primarily on the interface deformability and on the viscosity ratio between the two fluids (λ).     

Experimental study on the capillary thinning of entangled polymer solutions

The transient elongation behavior of entangled polymer and wormlike micelles (WLM) solutions has been investigated using capillary breakup extensional rheometry (CaBER). The transient force ratio X = 0.713 reveals the existence of an intermediate Newtonian thinning region for polystyrene and WLM solutions prior to the viscoelastic thinning. The exponential decay of X(t) in the first period of thinning defines an elongational relaxation time λ _x which is equal to elongational relaxation time λ _e obtained from exponential diameter decay D(t) indicating that the initial stress decay is controlled by the same molecular relaxation process as the strain hardening observed in the terminal regime of filament thinning. Deviations in true and apparent elongational viscosity are discussed in terms of X(t). A minimum Trouton ratio is observed which decreases exponentially with increasing polymer concentration leveling off at Tr_min = 3 for the solutions exhibiting intermediate Newtonian thinning and Tr_min ≈ 10 otherwise. The relaxation time ratio λ _e/ λ _s, where λ _s is the terminal shear relaxation time, decreases exponentially with increasing polymer concentration and the data for all investigated solutions collapse onto a master curve irrespective of polymer molecular weight or solvent viscosity when plotted versus the reduced concentration c[ η], with [ η] being the intrinsic viscosity. This confirms the strong effect of the nonlinear deformation in CaBER experiments on entangled polymer solutions as suggested earlier. On the other hand, λ _e ≈λ _s is found for all WLM solutions clearly indicating that these nonlinear deformations do not affect the capillary thinning process of these living polymer systems.     

The rheology of strongly-flocculated suspensions

The rheology of strongly-flocculated dispersions of colloidal particles has been investigated at particle concentrations where a continuous network is formed rather than a collection of discrete flocs. Such networks are shown to possess a true yield stress in both shear and in uniaxial compression (as realised in a centrifuge). Properties measured as a function of particle concentration and particle size include the yield stresses in shear (σ_y) and compression (P_y); the limiting and strain-dependent, instantaneous shear moduli G_O and G(γ); the elastic recovery at finite strains, and the rate of centrifugally-driven compaction. The yield stresses and moduli appear to show a power-law dependence on particle concentration with G_O and P_y, having the same power-law index and σ_y a somewhat lower one. The data are in part consistent with predictions based on the idea that the networks have a heterogeneous structure comprising a collection of interconnected fractal aggregates. The behaviour as a function of particle size and concentration is however not completely scaleable as might be expected on this basis. Thus, whereas the shear yield stress could be scaled to remove its dependence on particle radius a and volume fraction φ (over the measured range 0.25 μm ⩽ a ⩽ 3.4 μm; 0.05 ⩽ φ ⩽ 0.25) as could the strain dependent modulus (0.25 ⩽ a ⩽ 1.3 μm; 0.08 ⩽ 0.25), the particle-size and concentration dependence of P_y and G_O could only be scaled for particles with radii between 0.16 and 0.5 μm, smaller and larger particles having different and much higher power-law index in respect of their concentration dependencies. In the case of the smaller particles the failure of the scaling is thought to be due to an anomaly since these particles distort significantly under the influence of the strong van der Waals forces and this causes the aggregates to be more compact then they otherwise would be. The reasons for the failure at larger sizes is not clear.     

Large Strain Shear Compression Test of Sheet Metal Specimens

A hybrid experimental-computational procedure to establish accurate true stress-plastic strain curve of sheet metal specimen covering the large plastic strain region using shear compression test data is described. A new shear compression jig assembly with a machined gage slot inclined at 35° to the horizontal plane of the assembly is designed and fabricated. The novel design of the shear compression jig assembly fulfills the requirement to maintain a uniform volume of yielded material with characteristic maximum plastic strain level across the gage region of the Shear Compression Metal Sheet (SCMS) specimen. The approach relies on a one-to-one correlation between measured global load–displacement response of the shear compression jig assembly with SCMS specimen to the local stress-plastic strain behavior of the material. Such correlations have been demonstrated using finite element (FE) simulation of the shear compression test. Coefficients of the proposed correlations and their dependency on relative plastic modulus were determined. The procedure has been established for materials with relative plastic modulus in the range 5?×?10^?4?<?(E _p /E)?<?0.01. It can be readily extended to materials with relative plastic modulus values beyond the range considered in this study. Nonlinear characteristic hardening of the material could be established through piecewise linear consideration of the measured load–displacement curve. Validity of the procedure is established by close comparison of measured and FE-predicted load–displacement curve when the provisional hardening curve is employed as input material data in the simulation. The procedure has successfully been demonstrated in establishing the true stress-plastic strain curve of a demonstrator 0.0627C steel SCMS specimen to a plastic strain level of 49.2 pct.     

On the hyperbolic variation of elastic modulus in a non-homogeneous stratum

An exact formulation is presented for the governing dual integral equations representing the mixed boundary-value problem of the static stress distribution under a long rigid rectangular body lying on the free surface of a non-homogeneous stratum. The shear modulus of the stratum increases in the depth direction y from a value G₀ at the surface according to the hyperbolic variation.

$G\left(y\right)=\frac{G_{o}h}{h?y}$

It, therefore, simulates a practical soil, of arbitrary Poisson's ratio, which smoothly merges into a rigid bed at a depth h below the surface. The work shows that the problem is governed by both kinds of exponential integral function and that the effect of surface elastic properties is dominant in the solution of the governing equations. The limiting case of a homogeneous half-space is easily recovered from the general formulation as h tends to infinity.     

Analytical and experimental evaluation of double-notch shear specimens of orthotropic materials

A finite-difference analysis of the state of stress in a double-notch interlaminar shear strength specimen is developed. The effects of geometry and material parameters on the stress distributions are investigated. It has been found that, in agreement with previous determinations,^1–7 a uniform distribution of shear stress on the fracture plane does not exist. The shear stress distribution becomes more uniform for increased material anisotropy and for small (L/T) ratios, whereL is the distance between the notches andT is the specimen thickness. Also, it has been determined that the notch size (W) and the distance from the notches to the loaded ends of the specimen (h) do not influence the stress distributions significantly. The effects of variations in the (L/T) ratio, the notch size (W), and the length (h) were investigated experimentally. For a graphite/epoxy laminate of 0/90-deg square wave it has been found that the apparent shear strength determined by double-notch shear tests decreases significantly with an increase in (L/T) ratio. The decrease in the apparent shear strength with an increase inh, however, is very small. Also, the apparent shear strength is not affected significantly by increasing the notch sizeW.     

Existence and construction of quadratic invariants for the equation d2dξ2ν = λ(ξ)ddξν + f(ξ,ν)

If the equation ν_ζζ = λ(ξ)ν_ζ + F(ξ,ν) admits a first integral quadratic in the momentum c., the function F(ξ,ν) is defined; a new procedure is introduced to find the first integral. Among the examples is discussed the existence of quadratic invariants for the equation u″ λu′ = Q(u), where λ is a constant, and for the equation ν_ξ,ξ =βξ^σνⁿ. This last equation admits a first integral, quadratic in ν_ξ if σ = 0, or n = −2σ − 3, or n = −σ − 3. Particular eases are exhibited where no invariant quadratic in ν_ξ is found for the equation ν_ξξ = βξ^σνⁿ but it admits a first integral which is a polynomial of degree 4 or 6 in ν_ξ.     

On Near-Wall Treatment in (U)RANS-Based Closure Models

The present work is concerned with computational evaluation of a recently formulated near-wall relationship providing the value of the dissipation rate ε of the kinetic energy of turbulence k through its exact dependence on the Taylor microscale λ: ε = 10νk/λ ², (Jakirli? and Jovanovi?, J. Fluid Mech. 656:530–539, 2010). Dissipation rate determination benefits from the asymptotic behavior of the Taylor microscale resulting in its linear variation in terms of the wall distance (λ?∝?y) being valid throughout entire viscous sublayer. Accordingly, it can be applied as a unified near-wall treatment in all computational frameworks relying on a RANS-based model of turbulence (including also hybrid LES/RANS schemes) independent of modeling level—both main modeling concepts eddy-viscosity and Reynolds stress models can be employed. Presently, the feasibility of the proposed formulation was demonstrated by applying a conventional near-wall second-moment closure model based on the homogeneous dissipation rate ε _h ( ${\varepsilon_h =\varepsilon -0.5\partial \left( {{\nu \partial k}/ {\partial x_j }} \right)} / {\partial x_j }$ ; Jakirli? and Hanjali?, J. Fluid Mech. 539:139–166, 2002) and its instability-sensitive version, modeled in terms of the inverse turbulent time scale ω _h (ω _h ?=?ε _h /k; Maduta and Jakirli?, 2011), to a fully-developed channel flow with both flat walls and periodic hill-shaped constrictions mounted on the bottom wall in a Reynolds number range. The latter configuration is subjected to boundary layer separation from a continuous curved wall. The influence of the near-wall resolution lowering with respect to the location of the wall-closest computational node, coarsened even up to the viscous sublayer edge situated at $y_P^+ \approx 5$ in equilibrium flows, is analyzed. The results obtained follow closely those pertinent to the conventional near-wall integration with the wall-next node positioned at $y_P^+ \le 0.5$ .     

On Pressurized Curvilinearly Orthotropic Circular Disk, Cylinder and Sphere Made of Radially Nonuniform Material

A unified analysis is presented for the elastic response of a pressurized cylindrically anisotropic hollow disk under assumed conditions of plane stress, or a hollow cylinder under plane strain conditions, and a spherically anisotropic hollow sphere, made of material which is nonuniform in the radial direction according to the power law relationship. The solution for a cylinder under generalized plane strain is also presented. Two parameters play a prominent role in the analysis: the material nonuniformity parameter m, and the parameter ?? which accounts for the combined effects of material anisotropy, represented by the specified parameters (??, ??, ??), and material nonuniformity, represented by the parameter m. The radial and circumferential stresses are the linear combinations of two power functions of the radial coordinate, whose exponents (n ₁ and n ₂) depend on the parameters m and ??. New light is added to the stress amplification and shielding under combined effects of curvilinear anisotropy and radial nonuniformity. Different loading combinations are considered, including the equal pressure at both boundaries, and the uniform pressure at the inner or the outer boundary. While the stress state for the equal pressure loading is uniform in the case of isotropic uniform material (m=0, ??=1), and for one particular radially nonuniform and anisotropic material, it is strongly nonuniform for a general anisotropic or nonuniform material. If the aspect ratio of the inner and outer radii decreases (small hole in a large disk/cylinder or sphere), the magnitude of the circumferential stress at the inner radius increases for n ₁>0 (stress amplification), and decreases for n ₁<0 (stress shielding). Both can be achieved by various combinations of the material parameters m, ??, ??, and ??. While the stress amplification in the case of a pressurized external boundary occurs readily, it occurs only exceptionally in the case of a pressurized internal boundary. The effects of material parameters on the displacement response are also analyzed. The approximate character of the plane stress solution of a pressurized thin disk is discussed and the results are compared with those obtained by numerical solution of the exact three-dimensional disk model.     

State of stress at the vertex of a quarter-infinite crack in a half-space

Spherical coordinates are r, θ, φ. The half-space extends in θ < π/2. The crack occurs along φ = 0. The region to be investigated is the solid space-triangle (or cone) between the three planes θ = π/2, φ = +0 and φ = 2π ? 0, which planes are to be taken stress-free.In this space-angle a state of stress is considered in terms of the cartesian stress components σ_xx = r^λ?_xx(λ, θ, φ); σ_xy = r^λ(λ, θ, φ); etc. Possible values λ are determined from a characteristic (or eigenvalue) equation, expressing the condition that a determinant of infinite order is equal to zero. The root of λ which gives the most serious state of stress in the vertex region (the region r → 0) is the root closest to the limiting value Re λ > ?3/2. Knowledge of this state of stress, or at least of this value of λ is essential in the determination of the three-dimensional state of stress around a crack in a plate for distances of order of the plate thickenss.Along the front of the quarter-infinite crack (z-axis) the so called stress-intensity factor behaves like z^λ+½ (z → 0) and thus tends to zero, respectively to infinity, accordingly to Re λ being >?½ or <?½. But in the region z → 0 the notion stress-intensity factor loses its meaning. The required state of stress passes into the well-known state of plane strain around a crack tip if Poisson's ratio (v) tends to zero. The computed state of stress for the incompressible medium (v = ½) is confirmed by experiment.     

A shear-compression specimen for large strain testing

A new specimen geometry, the shear-compression specimen (SCS), has been developed for large strain testing of materials. The specimen consists of a cylinder in which two diametrically opposed slots are machined at 45° with respect to the longitudinal axis, thus forming the test gage section. The specimen was analyzed numerically for two representative material models, and various gage geometries. This study shows that the stress (strain) state in the gage, is three-dimensional rather than simple shear as would be commonly assumed. Yet, the dominant deformation mode in the gage section is shear, and the stresses and strains are rather uniform. Simple relations were developed and assessed to relate the equivalent true stress and equivalent true plastic strain to the applied loads and displacements. The specimen was further validated through experiments carried out on OFHC copper, by comparing results obtained with the SCS to those obtained with compression cylinders. The SCS allows to investigate a large range of strain rates, from the quasi-static regime, through intermediate strain rates (1–100 s⁻¹), up to very high strain rates (2×10⁴s⁻¹ in the present case).     

Extending the bounds of ‘steady’ RANS closures: Toward an instability-sensitive Reynolds stress model

The incapability of the conventional Unsteady RANS (Reynolds–Averaged Navier Stokes) models to adequately capture turbulence unsteadiness presents the prime motivation of the present work, which focuses on formulating an instability-sensitive, eddy-resolving turbulence model on the Second-Moment Closure level. The model scheme adopted, functioning as a ‘sub-scale’ model in the Unsteady RANS framework, represents a differential near-wall Reynolds stress model formulated in conjunction with the scale-supplying equation governing the homogeneous part of the inverse turbulent time scale ω_h (ω_h = ɛ_h/k). The latter equation was straightforwardly obtained from the model equation describing the dynamics of the homogeneous part of the total viscous dissipation rate ɛ, defined as ɛ_h = ɛ − 0.5ν∂²k/(∂x_j∂x_j) (Jakirlic and Hanjalic, 2002), by applying the derivation rules to the expression for ω_h. The model capability to account for vortex length and time scales variability was enabled through an additional term in the corresponding length-scale determining equation, providing a selective enhancement of its production, pertinent particularly to the highly unsteady separated shear layer region, modeled in terms of the von Karman length scale (comprising the second derivative of the velocity field) in line with the SAS (Scale-Adaptive Simulation) proposal (Menter and Egorov, 2010). The present model formulation, termed as SRANS model (Sensitized RANS), does not comprise any parameter depending explicitly on grid spacing. The predictive capabilities of the newly proposed length-scale determining model equation, solved in conjunction with Jakirlic and Hanjalic’s (2002) Reynolds stress model equation, are presently demonstrated by computing the flow configurations of increasing complexity featured by boundary layer separation from sharp-edged and continuous curved surfaces: backward-facing step flow, flow over a wall-mounted fence, flow over smoothly contoured periodically arranged hills and flow in a 3-D diffuser. The model performances are also assessed in capturing the natural decay of the homogeneous isotropic turbulence and the near-wall Reynolds stress anisotropy in a plane channel. In most cases considered the fluctuating velocity field was obtained starting from steady RANS results.     

Experimental research on the thermal performance of converging slot holes with different divergence angles

Thermal performances of two kinds of converging slot-hole (console) with different divergence angles have been measured using transient liquid crystal measurement technique which can process the nonuniform initial wall temperature. Four momentum ratios are tested. Consoles with different divergence angles produce different cooling effectiveness distributions in the upstream region. However, the cooling effectiveness distributions of the two consoles are similar in the downstream. The laterally averaged cooling effectiveness results show that the differences between the two consoles are very small and the best momentum ratio for both consoles’ cooling effectiveness distribution is around two. With the momentum ratio increasing, the normalized heat transfer coefficient h/h₀ of both consoles increases, but the h/h₀ value of small divergence case is higher and becomes progressively higher than that of large divergence case. Moreover, the effect of the couple vortices on the heat transfer coefficient distributions is more significant for the large divergence case. Both consoles provide the surface a certain degree of thermal protection, especially in the upstream region. The distributions of heat flux ratio q/q₀ are similar with those of η because the influence of η on q/q₀ is much larger than that of h/h₀ on q/q₀.     

Effect of curvature on dual solutions of flow through a curved circular tube

Dual solutions, i.e., two-vortex and four-vortex solutions, of the flow through a curved circular tube are numerically obtained by the spectral method for 0 ⩽ δ ⩽ 0.8 and 500 ⩽ Dn ⩽ 10000, where δ is the non-dimensional curvature of the tube and Dn the Dean number. It is found that the critical Dean number above which the four-vortex solution exists takes the lowest value 956 at δ = 0 and increases with δ. In terms of the Reynolds number of the flow, however, the critical Reynolds number decreases from infinity as δ increases from zero, takes the minimum value of about 250 at δ ≈ 0.52, and then increases again.     

Crack separation energy-rates for the DBCS model under biaxial modes of loading

A Modified version of the Dugdale-Bilby-Cottrell-Swinden (DBCS) model simulating the effect of plasticity at the tip of a crack in an infinite region was used by kfouri and rice (1978) to calculate the crack separation energy-rate G^Δ corresponding to a finite crack growth step Δa during plane strain mode I crack extension. The loading consisted of a remote uniaxial tension σp applied normally to the plane of the crack. Using Rice's path-independent integral J to characterize the applied load in the crack tip region, and assuming the length R of the crack tip plastic zone to be small compared with the length of the crack, an analytical expression was derived relating the ratios (J/G^Δ) and (2a/R) for small values of (2a/R). The analytical solution was incomplete in itself in that the value assumed in the plastic region of the DBCS model for the normal stress Y acting on the extending crack surfaces was the product of the yield stress in uniaxial tension σ_Y and an unknown parameter C, the value of which depends on the effect of the local hydrostatic stresses in the case of plane strain conditions. The analytical solution was compared with a numerical solution obtained from a plane strain elastic-plastic finite element analysis on a centre-cracked plate (CCP) of material obeying the von Mises yield criterion. The value used for the yield stress was 310 MN/m² and moderate isotropic linear hardening was applied with a tangent modulus of 4830 MN/m². A uniaxial tension σp was applied on the two appropriate sides of the plate. The comparisons showed that the analytical and finite element solutions were mutually consistent and they enabled the value of C to be established at 1.91. In the present paper similar comparisons are made between the analytical solution and the finite element solutions for the CCP of the same material under different biaxial modes of loading. By assuming further that the form of the analytical solution does not depend on the details of the geometry and of the loading at remote boundaries, a comparison has also been made with the results of a finite element analysis on a compact tension specimen (CTS) made of the same material as the CCP. The different values of C obtained in each case are consistent with investigations by other authors on the effect of load biaxiality on crack tip plasticity.     

Meltblown fibers: Influence of viscosity and elasticity on diameter distribution

Both melt viscosity (η_o) and elasticity (correlated here with the longest melt relaxation time λ₁) were found to control the diameter distribution of meltblown fibers. Fibers were formed by melt blowing binary polystyrene (PS) blends containing widely differing component molecular weights using a custom-built laboratory apparatus. Varying the concentration and molecular weight of a high molecular weight PS provided independent control over η_o and λ₁. These rheological parameters influence the average diameter (d_av) and the distribution of diameters (coefficient of variation, CV) of meltblown fibers in different ways. Increasing η_o leads to an increase in d_av but has little impact on CV. On the other hand, increasing λ₁ beyond a threshold value reduces CV while simultaneously increasing d_av. A one-dimensional slender-jet theoretical model with both upper convected Maxwell and Phan–Thien and Tanner constitutive equations was developed to investigate the influence of viscoelasticity and processing parameters on the properties of meltblown fibers. This model predicts a strong dependence of fiber diameter on the air shear stress and variations in fiber diameter with viscoelasticity that are in qualitative agreement with the experimental results. We believe these results suggest that carefully controlling the viscoelastic profile of polymers used in melt blowing is a viable approach for producing nanofibers with narrow fiber diameter distributions using current commercial equipment.     

Analytical solution of two intersecting cylindrical shells subjected to transverse moment on nozzle

In this paper a theoretical solution for two normally intersecting cylindrical shells subjected to transverse moment on the branch pipe is presented, which based on thin shell theory. The accurate shell equations, boundary conditions and calculating methods are adopted so that the solution presented can be applicable up to d/D⩽0.8 and λ=d/(DT)^1/2⩽8. The presented results are in very good agreement with experimental and numerical results for ORNL-1 Model. They are also in agreement with the results obtained by WRC Bulletin 297 when d/D is small.