Computer simulation of the damage evolution of fiber-reinforced polypropylene-matrix composites with matrix defects

The damage evolution of fiber-reinforced polypropylene-matrix composites with matrix defects was studied via a Monte Carlo technique combined with a finite element method. A finite element model was constructed to predict the effects of various matrix defect shapes on the stress distributions. The results indicated that a small matrix defect had almost no effect on fiber stress distributions other than interfacial shear stress distributions. Then, a finite element model with a statistical distribution of the fiber strength was constructed to investigate the influences of the spatial distribution and the volume fraction of matrix defects on composite failure. The results showed that it was accurate to use the shear-lag models and Green's function methods to predict the tensile strength of composites even though the axial stresses in the matrix were neglected.

Kinetic isotope effects of proton transfer reaction for amines by ultrasonic relaxation methods: Unexpected evidence from the results in ethylamine solutions

Ultrasonic absorption coefficients for ethylamine in heavy water (D2O) and in light water (H2O) have been measured in the frequency range from 0.8 to 220 MHz at 25 degrees C. A single relaxational process has been observed in these two kinds of solutions. From the concentration dependence of the ultrasonic relaxation parameters, and following the reaction mechanism proposed by Eigen et al. for ethylamine in H2O, the causes of the relaxations have been attributed to a perturbation of an equilibrium associated with a deuteron or proton transfer reaction. The rate and equilibrium constants have been estimated from deuterioxide or hydroxide ion concentration dependence of the relaxation frequency, and the kinetic isotope effects have been determined. In addition, the standard volume changes of the reactions have been calculated from the concentration dependence of the maximum absorption per wavelength, and the adiabatic compressibility has also been determined from the density and sound velocity for ethylamine in D2O and in H2O, respectively. These results are compared with those for propylamine and butylamine and are discussed in relation to the different kinetic properties between D2O and H2O, the reaction radii derived by Debye theory, and the structural properties of the reaction intermediate.

Isotope effect in proton-transfer and association-dissociation reactions for butylamine by ultrasonic relaxation methods

Ultrasonic absorption coefficients were measured for butylamine in heavy water (D2O) in the frequency range from 0.8 to 220 MHz and at concentrations from 0.0278 to 2.5170 mol dm(-3) at 25 degrees C; two kinds of relaxation processes were observed. One was found in relatively dilute solutions (up to 0.5 mol dm(-3)), which was attributed to the hydrolysis of butylamine. In order to compare the results, absorption measurements were also carried out in light water (H2O). The rate and thermodynamic parameters were determined from the concentration dependence of the relaxation frequency and the maximum absorption per wavelength. The isotope effects on the diffusion-controlled reaction were estimated and the stability of the intermediate of the hydrolysis was considered while comparing it with the results for propylamine in H2O and D2O. Another relaxation process was observed at concentrations greater than 1 mol dm(-3) in D2O. In order to examine the solution characteristics, proton NMR measurements for butylamine were also carried out in D2O. The chemical shifts for the gamma- and delta-proton in butylamine molecule indicate the existence of an aggregate. From profiles of the concentration dependence of the relaxation frequency and the maximum absorption per wavelength of sound absorption, the source of the relaxation was attributed to an association-dissociation reaction, perhaps, associated with a hydrophobic interaction. The aggregation number, the forward and reverse rate constants and the standard volume change of the reaction were determined. It was concluded from a comparison with the results in H2O that the hydrophobic interaction of butylamine in D2O is stronger than that in H2O. Also, the isotope effect on this reaction was interpreted in terms of the solvent structure.

基于广义Lippmann-Schwinger方程的地震波场模拟算法与应用

In the last several decades, due to the fast development of computer, numerical simulation has been an indispensable tool in scientific research. Numerical simulation methods which based on partial difference operators such as Finite Difference Method (FDM) and Finite Element Method (FEM) have been widely used. However, in the realm of seismology and seismic prospecting, one usually meets with geological models which have piece-wise heterogeneous structures as well as volume heterogeneities between layers, the continuity of displacement and stress across the irregular layers and seismic wave scattering induced by the perturbation of the volume usually bring in error when using conventional methods based on difference operators. The method discussed in this paper is based on elastic theory and integral theory. Seismic wave equation in the frequency domain is transformed into a generalized Lippmann-Schwinger equation, in which the seismic wavefield contributed by the background is expressed by the boundary integral equation and the scattering by the volume heterogeneities is considered. Boundary element-volume integral method based on this equation has advantages of Boundary Element Method (BEM), such as reducing one dimension of the model, explicit use the displacement and stress continuity across irregular interfaces, high precision, satisfying the boundary at infinite, etc. Also, this method could accurately simulate the seismic scattering by the volume heterogeneities. In this paper, the concrete Lippmann-Schwinger equation is specifically given according to the real geological models. Also, the complete coefficients of the non-smooth point for the integral equation are introduced. Because Boundary Element-Volume integral equation method uses fundamental solutions which are singular when the source point and the field are very close，both in the two dimensional and the three dimensional case, the treatment of the singular kernel affects the precision of this method. The method based on integral transform and integration by parts could treat the points on the boundary and inside the domain. It could transform the singular integral into an analytical one both in two dimensional and in three dimensional cases and thus it could eliminate the singularity. In order to analyze the elastic seismic wave scattering due to regional irregular topographies, the analytical solution for problems of this type is discussed and the analytical solution of P waves by multiple canyons is given. For the boundary reflection, the method used here is infinite boundary element absorbing boundary developed by a pervious researcher. The comparison between the analytical solutions and concrete numerical examples validate the efficiency of this method. We thoroughly discussed the sampling frequency in elastic wave simulation and find that, for a general case, three elements per wavelength is sufficient, however, when the problem is too complex, more elements per wavelength are necessary. Also, the seismic response in the frequency domain of the canyons with different types of random heterogeneities is illustrated. We analyzed the model of the random media, the horizontal and vertical correlation length, the standard deviation, and the dimensionless frequency how to affect the seismic wave amplification on the ground, and thus provide a basis for the choice of the parameter of random media during numerical simulation.