Crack path instabilities in DCDC experiments in the low speed regime

We studied the low speed fracture regime (10⁻⁴-10⁻⁹ m s⁻¹) in different glassy materials (soda-lime glass, glass ceramics) with variable but controlled length scale of heterogeneity. The chosen mechanical system enabled us to work in pure mode I (tensile) and at a fixed load on double cleavage drilled compression specimen. The internal residual stresses of studied samples were carefully relaxed by appropriate thermal treatment. By means of optical and atomic force microscopy techniques fracture surfaces have been examined. We have shown for the first time that the crack front line underwent an out-of-plane oscillating behavior as a result of a reproducible sequence of instabilities. The wavelength of such a phenomenon is in the micrometer range and its amplitude in the nanometer range. These features were observed for different glassy materials providing that a typical length scale characterizing internal heterogeneities was lower than a threshold limit estimated to few nanometers. This effect is the first clear experimental evidence of crack path instabilities in the low speed regime in a uniaxial loading experiment. This phenomenon has been interpreted by referring to the stability criterion for a straight crack propagation as presented by Adda-Bedia et al. [Phys. Rev. Lett. 76 (1996) 1497].     

20-relative neighborhood graphs are hamiltonian

For any two points p and q in the Euclidean plane, define LUN_pq = { v | v ∈ R², d_pv < d_pq and d_qv < d_pq}, where d_uv is the Euclidean distance between two points u and v . Given a set of points V in the plane, let LUN_pq(V) = V ∩ LUN_pq. Toussaint defined the relative neighborhood graph of V, denoted by RNG(V) or simply RNG, to be the undirected graph with vertices V such that for each pair p,q ∈ V, (p,q) is an edge of RNG(V) if and only if LUN_pq (V) = ?. The relative neighborhood graph has several applications in pattern recognition that have been studied by Toussaint. We shall generalize the idea of RNG to define the k-relative neighborhood graph of V, denoted by kRNG(V) or simply kRNG, to be the undirected graph with vertices V such that for each pair p,q ∈ V, (p,q) is an edge of kRNG(V) if and only if | LUN_pq(V) | < k, for some fixed positive number k. It can be shown that the number of edges of a kRNG is less than O(kn). Also, a kRNG can be constructed in O(kn₂) time. Let E_c = {e_pq| p ∈ V and q ∈ V}. Then G_c = (V,E_c) is a complete graph. For any subset F of E_c, define the maximum distance of F as max_epq∈Fd_pq. A Euclidean bottleneck Hamiltonian cycle is a Hamiltonian cycle in graph G_c whose maximum distance is the minimum among all Hamiltonian cycles in graph G_c. We shall prove that there exists a Euclidean bottleneck Hamiltonian cycle which is a subgraph of 20RNG(V). Hence, 20RNGs are Hamiltonian.     

Modeling hyperelastic behavior of rubber: A novel invariant‐based and a review of constitutive models

A constitutive phenomenological model completing the Gent‐Thomas concept is carried out to formulate laws governing the hyperelastic behavior of incompressible rubber materials. It is shown that the phenomenological Gent‐Thomas model (1958) and the constrained chain model (1992) give similar precise results at small to moderate deformation. On the other hand, comparisons of the outcome of the proposed model with that of the molecular model from the combined concepts of Flory‐Erman and Boyce‐Arruda (2000), and with those of the phenomenological models of Ogden (1982), Yeoh‐Fleming (1997), Pucci‐Saccomandi (2002) and Beda (2005) are made. Residual inconveniences raised by attractive continuum models in rubber elasticity literature have been successfully overcome. Results from both the statistical and phenomenological mechanics concepts are compared with the data of some useful classical materials (rubbers of Treloar, Rivlin‐Saunders, Pak‐Flory and Yeoh‐Fleming). The results permit one to see salient equivalence of the two theories for a more reliable prediction of stress‐stretch response for all states of any mode of deformation. A complete and exhaustive analysis of the Mooney plot that combines small and very large extension‐compression has been quite essential in assessing the validity of models. A method of identification of material parameters is presented and data of the simple tension suffice for the determination of the parameter values. It is shown that the ordinary identification procedures, such as the usual least squares, a very much used numerical method in materials investigation, can be unsuitable in some cases of hyperelastic modeling. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 1713–1732, 2007     

Finite element computations of two-dimensional arterial flow in the presence of a transverse magnetic field

A finite element solution of the Navier-Stokes equations for steady flow under the magnetic effect through a double-branched two-dimensional section of a three-dimensional model of the canine aorta is discussed. The numerical scheme involves transforming the physical co-ordinates to a curvilinear boundary-fitted co-ordinate system. The shear stress at the wall is calculated for a Reynolds number of 1000 with the branch-to-main aortic flow rate ratio as a parameter. The results are compared with earlier works involving experimental data and found to be in reasonable qualitative agreement. The steady flow, shear stress and branch flow under the effect of a magnetic field have been discussed in detail.     

Finite element method in the solution of the Euler and Navier-Stokes equations for internal flow

Finite element solutions of the Euler and Navier-Stokes equations are presented, using a simple dissipation model. The discretization is based on the weak-Galerkin weighted residual method and equal interpolation functions for all the unknowns are permitted. The nonlinearity is iterated upon using a Newton method and at each iteration the linear algebraic system is solved by a direct solver with all unknowns fully coupled. Results are presented for two-dimensional transonic inviscid flows and two- and three-dimensional incompressible viscous flows. Convergence of the algorithm is shown to be quadratic, reaching machine accuracy in very few iterations. The inviscid results demonstrate the existence of nonunique numerical solutions to the steady Euler equations.     

Interpretation of the polyelectrolyte and antipolyelectrolyte effects of poly(N‐vinylimidazole‐co‐sodium styrenesulfonate) hydrogels

Hydrogels of N‐vinylimidazole (VI) and sodium styrenesulfonate (SSS) were synthesized in aqueous solution by radical crosslinking copolymerization with N,N′‐methylene‐bis(acrylamide) as crosslinker. Swelling in several saline solutions was measured for hydrogel samples synthesized with different comonomer concentrations (C_T = 10, 25, or 40%) and with SSS mole fractions covering a broad range (f_SSS = 0–0.7), while the crosslinker ratio was 2 wt % in all cases. The degree of swelling in aqueous solution with a specific ionic strength (μ), plotted versus the SSS composition of the feed, shows a minimum for any set of samples synthesized with a fixed C_T. The dependence of swelling on μ shows both polyelectrolyte (f_SSS beyond the minimum) and antipolyelectrolyte behaviors (in the low f_SSS limit). It was found that the nonGaussian factor of the crosslinking density and the polymer‐solvent interaction parameter increase with f_SSS for any C_T. Moreover, in the low f_SSS limit, the osmotic swelling pressure is governed not only by the ionic contribution, but also by the polymer‐solvent mixing and, the concentration of mobile counterions inside the gel is not proportional to the net fixed charge but to the addition of cationic and anionic side groups, what discards the formation of ionic pairs. The antipolyelectrolyte effect is interpreted as due to the increasing protonation of VI as μ goes up. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 1683–1693, 2007