Lattice type of fracture model for concrete

Concrete is usually described as a three-phase material, where matrix, aggregate and interface zones are distinguished. The beam lattice model has been applied widely by many investigators to simulate fracture processes in concrete. Due to the extremely large computational effort, however, the beam lattice model faces practical difficulties. In our investigation, a new lattice called generalized beam (GB) lattice is developed to reduce computational effort. Numerical experiments conducted on a panel subjected to uniaxial tension show that the GB lattice model can reproduce the load-displacement curves and crack patterns in agreement to what are observed in tests. Moreover, the effects of the particle overlay on the fracture process are discussed in detail. (C) 2007 Elsevier Ltd. All rights reserved.

Toward an understanding of the turbidity measurement of fleterocoagulation rate constants of dispersions containing particles of different sizes

Our previous studies have shown that the determination of coagulation rate constants by turbidity measurement becomes impossible for a certain operating wavelength (that is, its blind point) because at this wavelength the change in the turbidity of a dispersion completely loses its response to the coagulation process. Therefore, performing the turbidity measurement in the wavelength range near the blind point should be avoided. In this article, we demonstrate that the turbidity measurement of the rate constant for coagulation of a binary dispersion containing particles of two different sizes (heterocoagulation) presents special difficulties because the blind point shifts with not only particle size but also with the component fraction. Some important aspects of the turbidity measurement for the heterocoagulation rate constant are discussed and experimentally tested. It is emphasized that the T-matrix method can be used to correctly evaluate extinction cross sections of doublets formed during the heterocoagulation process, which is the key data determining the rate constant from the turbidity measurement, and choosing the appropriate operating wavelength and component fraction are important to achieving a more accurate rate constant. Finally, a simple scheme in experimentally determining the sensitivity of the turbidity changes with coagulation over a wavelength range is proposed.

Lattice instability at a fast moving crack tip

A molecular dynamics method is used to analyze the dynamic propagation of an atomistic crack tip. The simulation shows that the crack propagates at a relatively constant global velocity which is well below the Rayleigh wave velocity. However the local propagation velocity oscillates violently, and it is limited by the longitudinal wave velocity. The crack velocity oscillation is caused by a repeated process of crack tip blunting and sharpening. When the crack tip opening displacement exceeds a certain critical value, a lattice instability takes place and results in dislocation emissions from the crack tip. Based on this concept, a criterion for dislocation emission from a moving crack tip is proposed. The simulation also identifies the emitted dislocation as a source for microcrack nucleation. A simple method is used to examine this nucleation process. (C) 1996 American Institute of Physics.

A new expression of hardening coefficients for fcc-crystal and calibration of the material constants

In order to describe the effect of latent hardening on the macro-plastic behavior of foc-crystal, a new expression for hardening coefficient is proposed in which there are 12 material constants, each having clear physical meaning. And a method of material constant calibration is suggested and used to determine the material constants of copper and aluminum crystal. The simulated load-elongation curves along various crystallographic orientations are comparable with the experimental ones.

Algorithm for simulating fracture processes in concrete by lattice modeling

To simulate fracture behaviors in concrete more realistically, a theoretical analysis on the potential question in the quasi-static method is presented, then a novel algorithm is proposed which takes into account the inertia effect due to unstable crack propagation and meanwhile requests much lower computational efforts than purely dynamic method. The inertia effect due to load increasing becomes less important and can be ignored with the loading rate decreasing, but the inertia effect due to unstable crack propagation remains considerable no matter how low the loading rate is. Therefore, results may become questionable if a fracture process including unstable cracking is simulated by the quasi-static procedure excluding completely inertia effects. However, it requires much higher computational effort to simulate experiments with not very high loading rates by the dynamic method. In this investigation which can be taken as a natural continuation, the potential question of quasi-static method is analyzed based on the dynamic equations of motion. One solution to this question is the new algorithm mentioned above. Numerical examples are provided by the generalized beam (GB) lattice model to show both fracture processes under different loading rates and capability of the new algorithm.

A simple method to simulate shrinkage-induced cracking in cement-based composites by lattice-type modeling

A new numerical procedure is proposed to investigate cracking behaviors induced by mismatch between the matrix phase and aggregates due to matrix shrinkage in cement-based composites. This kind of failure processes is simplified in this investigation as a purely spontaneous mechanical problem, therefore, one main difficulty during simulating the phenomenon lies that no explicit external load serves as the drive to propel development of this physical process. As a result, it is different from classical mechanical problems and seems hard to be solved by using directly the classical finite element method (FEM), a typical kind of "load -> medium -> response" procedures. As a solution, the actual mismatch deformation field is decomposed into two virtual fields, both of which can be obtained by the classical FEM. Then the actual response is obtained by adding together the two virtual displacement fields based on the principle of superposition. Then, critical elements are detected successively by the event-by-event technique. The micro-structure of composites is implemented by employing the generalized beam (GB) lattice model. Numerical examples are given to show the effectiveness of the method, and detailed discussions are conducted on influences of material properties.

Numerical investigation of crack growth in concrete subjected to compression by the generalized beam lattice model

The beam lattice-type models, such as the Euler-Bernoulli (or Timoshenko) beam lattice and the generalized beam (GB) lattice, have been proved very effective in simulating failure processes in concrete and rock due to its simplicity and easy implementation. However, these existing lattice models only take into account tensile failures, so it may be not applicable to simulation of failure behaviors under compressive states. The main aim in this paper is to incorporate Mohr-Coulomb failure criterion, which is widely used in many kinds of materials, into the GB lattice procedure. The improved GB lattice procedure has the capability of modeling both element failures and contact/separation of cracked elements. The numerical examples show its effectiveness in simulating compressive failures. Furthermore, the influences of lateral confinement, friction angle, stiffness of loading platen, inclusion of aggregates on failure processes are respectively analyzed in detail.

Scattering by two degenerate anisotropic modes in square-lattice dielectric photonic crystals

We systematically investigate the square-lattice dielectric photonic crystals that have been used to demonstrate flat slab imaging experimentally. A right-handed Bloch mode is found in the left-handed frequency region by using the plane wave expansion method to analyze the photonic band structure and equifrequency contours. Using the multiple scattering theory, numerical simulations demonstrate that the left-handed mode and the right-handed mode are excited simultaneously by a point source and result in two kinds of transmitted waves. Impacted by the evanescent waves, superposition of these transmitted waves brings on complicated near field distributions such as the so-called imaging and its disappearance.

Optical properties of a square-lattice photonic crystal within the partial bandgap

Optical properties of a two-dimensional square-lattice photonic crystal are systematically investigated within the partial bandgap through anisotropic characteristics analysis and numerical simulation of field pattern. Using the plane-wave expansion method and Hellmann-Feynman theorem, the relationships between the incident and refracted angles for both phase and group velocities are calculated to analyze light propagation from air to photonic crystals. Three kinds of flat slab focusing are summarized and demonstrated by numerical simulations using the multiple scattering method. (c) 2007 Optical Society of America

Negative refraction in a honeycomb-lattice photonic crystal

In this paper we theoretically investigate a photonic crystal with dielectric rods in a honeycomb lattice. Two left-handed frequency regions are found in the second and third photonic band by using the plane wave expansion method to analyze the photonic band structure and equifrequency contours. Subwavelength imaging by the photonic crystal flat lens are systematically studied by numerical simulations using the multiple scattering method. Different from the photonic crystals with noncircular dielectric rods in air, this structure is almost isotropic at the optimal frequency for superlensing. As a comparison, flat slab focusing is also demonstrated at other frequencies in the two left-handed regions. (c) 2006 Elsevier Ltd. All rights reserved.

Determination and analysis of the optical constants of thin films of nickel(II) and copper(II) hydrazone complexes by spectroscopic ellipsometry

Thin films of four nickel(II) and copper(II) hydrazone complexes, which will hopefully be used as recording layers for the next-generation of high-density recordable disks, were prepared by using the spin-coating method. Absorption spectra of the thin films on K9 optical glass substrates in the 300-700 nm wavelength region were measured. Optical constants (complex refractive indices N) and thickness d of the thin films prepared on single-crystal silicon substrates in the 275-675 nm wavelength region were investigated on a rotating analyzer-polarizer scanning ellipsometer by fitting the measured ellipsometric angles (Psi(lambda) and Delta(lambda)) with a 3-layer model (Si/dye film/air). The dielectric functions epsilon and absorption coefficients alpha as a function of the wavelength were then calculated. Additionally, a design to achieve high reflectivity and optimum dye film thickness with an appropriate reflective layer was performed with the Film Wizard software using a multilayered model (PC substrate/reflective layer/dye film/air) at 405 nm wavelength.

Band structure of photonic crystal fibers with silica rings in triangular lattice

The characteristics of the cladding band structure of air-core photonic crystal fibers with silica rings in triangular lattice are investigated by using a standard plane wave method. The numerical results show that light can be localized in the air core by the photonic band gaps of the fiber. By increasing the air-filling fraction, the band gap edges of the low frequency photonic band gaps shift to shorter wavelength.. whereas the band gap width decreases linearly. In order to make a specified light fall in the low frequency band gaps of the fiber, the interplay of the silica ring spacing and the air-filling fraction is also analyzed. It shows that the silica ring spacing increases monotonously when the air-filling fraction is increased, and the spacing range increases exponentially. This type fiber might have potential in infrared light transmission. (c) 2006 Elsevier B.V. All rights reserved.

Design of a compact polarizing beam splitter based on a photonic crystal ring resonator with a triangular lattice

We propose and simulate a new kind of compact polarizing beam splitter (PBS) based on a photonic crystal ring resonator (PCRR) with complete photonic bandgaps. The two polarized states are separated far enough by resonant and nonresonant coupling between the waveguide modes and the microring modes. Some defect holes are utilized to control the beam propagation. The simulated results obtained by the finite-difference time-domain method show that high transmission (over 95%) is obtained and the polarization separation is realized with a length as short as 3.1 mu m. The design of the proposed PBS can be flexible, thanks to the advantages of PCRRs.

The impact of imperfect symmetry on band edge modes of a two-dimensional photonic crystal with square lattice

A theoretical analysis has been performed by means of the plane-wave expansion method to examine the dispersion properties of photons at high symmetry points of an InP based two-dimensional photonic crystal with square lattice. The Q factors are compared qualitatively. The mechanism of surface-emitting is due to the photon manipulation by periodic dielectric materials in terms of Bragg diffraction. A surface-emitting photonic crystal resonator is designed based on the phenomenon of slow light. Photonic crystal slabs with different unit cells are utilized in the simulation. The results indicate that the change of the air holes can affect the polarization property of the modes. So we can find a way to improve the polarization by reducing the symmetry of the structure.