Computer Simulation Of Random Sphere Packing In An Arbitrarily Shaped Container

Most simulations of random sphere packing concern a cubic or cylindric container with periodic boundary, containers of other shapes are rarely studied. In this paper, a new relaxation algorithm with pre-expanding procedure for random sphere packing in an arbitrarily shaped container is presented. Boundaries of the container are simulated by overlapping spheres which covers the boundary surface of the container. We find 0.4 similar to 0.6 of the overlap rate is a proper value for boundary spheres. The algorithm begins with a random distribution of small internal spheres. Then the expansion and relaxation procedures are performed alternately to increase the packing density. The pre-expanding procedure stops when the packing density of internal spheres reaches a preset value. Following the pre-expanding procedure, the relaxation and shrinking iterations are carried out alternately to reduce the overlaps of internal spheres. The pre-expanding procedure avoids the overflow problem and gives a uniform distribution of initial spheres. Efficiency of the algorithm is increased with the cubic cell background system and double link data structure. Examples show the packing results agree well with both computational and experimental results. Packing density about 0.63 is obtained by the algorithm for random sphere packing in containers of various shapes.

Degenerate four-wave-mixing signals from dc- and ac-driven semiconductor superlattice

Within the one-dimensional tight-binding model;rnd chi-3 approximation, we have calculated four-wave-mixing (FWM) signals for a semiconductor superlattice in the presence of both static and high-frequency electric fields. When the exciton effect is negligible, the time-periodic field dynamically delocalizes the otherwise localized Wannier-Stark states, and accordingly quasienergy band structures are formed, and manifest in the FWM spectra as a series of equally separated continua. The width of each continuum is proportional to the joint width of the valence and conduction minibands and is independent of the Wannier-Stark index. The realistic homogeneous broadening blurs the continua into broad peaks, whose line shapes, far from the Lorentzian, vary with the delay time in the FWM spectra. The swinging range of the peaks is just the quasienergy bandwidth. The dynamical delocalization (DDL) also induces significant FWM signals well beyond the excitation energy window. When the Coulomb interaction is taken into account, the unequal spacing between the excitonic Wannier-Stark levels weakens the DDL effect, and the FWM spectrum is transformed into groups of discrete lines. Strikingly, the groups are evenly spaced by the ac field frequency, reflecting the characteristic of the quasienergy states. The homogeneous broadening again smears out the line structures, leading to the excitonic FWM spectra quite similar to those without the exciton effect. However, all these features predicted by the dynamical theory do not appear in a recent experiment [Phys. Rev. Lett. 79, 301 (1997)], in which, by using the static approximation the observed Wannier-Stark ladder with delay-time-dependent spacing in the FWM spectra is attributed to a temporally periodic dipole field, produced by the Bloch oscillation of electrons in real space. The contradiction between the dynamical theory and the experiments is discussed. In addition, our calculation indicates that the dynamical localization coherently enhances the time-integrated FWM signals. The feasibility of using such a technique to study the dynamical localization phenomena is shown. [S0163-1829(99)10607-6].

深基坑土钉支护机理的现场测试及三维数值模拟研究

Soil nailing is a sort of new support system developed in recent year, and it has already been extensive application in the supporting system of deep-foundation pit . Based on the in situ measured data of hairun mansion deep-foundation pit excavation in Beijing, which include steel bar stress and lateral deformation, the following categories are systematically studied: (1) Magnitude of lateral deformation and displacement distribution vary with time at different excavation depth of deep-foundation pit excavation. Magnitude and distribution of soil nailing force vary with time at different excavation depth of deep-foundation pit excavation. By contract with the in situ measured data, we process the three-dimension numerical simulation of deep-foundation pit excavation with soil nailing support by Fast Lagrangian Analysis of Continua 3D(FLAC~(3D)). By above in situ measurement and studying, we found that the mechanism of soil nailing support is: by stress sharing, stress transmission and diffusion, and displacement restriction, composite soil nailing system can improve the strength and rigidity of soil and decrease the deformation and failure state by enhancing the whole stability of slope.

青藏高原壳幔形变与层圈耦合数值模拟

Based on multi-principle (such as structures, tectonics and kinematics) exploratory data and related results of continental dynamics in the Tibetan plateau, the author reconstructed the geological-geophysical model of lithospherical structure and tectonic deformation, and the kinetics boundary conditions for the model. Then, the author used the numerical scheme of Fast Lagrangian Analysis of Continua (FLAC), to stimulate the possible process of the stress field and deformational field in the Tibetan plateau and its adjacent area, since the convergence-collision between the Indian continent and Eurasia continent about 50Ma ago. With the above-mentioned results, the author discussed the relationship between crustal movement in shallow layer and the deformational process in interior layers, and its possible dynamic constraints in deep. At the end of the paper, an integrative model has been put forward to explain the outline images of crust-mantle deformation and coupling in the Tibetan Plateau. (1) The characteristics of crust-mantle structure of the Tibetan plateau have been shown to be very complex, and vertical and horizontal difference is significant. The general characteristics of crust-mantle of the Tibetan plateau may be that it's layering in depth direction, and shows blocking from south to north and belting from east to west, mainly according to the results of about 20 seismic sections, such as wide-angle seismic profiles, CMP, seismic tomography and so on. (2) The crust had shortened about 2200km, while the shortening is different for different block from south to north in the Tibetan plateau. It is about 11.5mm/a in Himalayan block, about 9.0mm/a in Lhas-Gangdese block, about 7.0mm/a in Qiangtang block and Songpan-Ganzi-Kekexili block, about 8.0mm/a in Kunlun-Qaidam, and about ll.Omm/a in Qilian block, since the convergence-collision between the Indian continent and Eurasia continent about 50Ma ago. Which - in demonstrates the shortening rate decreases from south to north, but this rate increases near the north edge of the Tibetan plateau. The crust thickening rate is about 0.4mm/a in the whole Tibetan plateau; and this rate is about 0.5mm/a in Himalayan block, about 0.4mm/a in Lhas-Gangdese block, about 0.3mm/a in Qiangtang block, about 0.2mm/a in Songpan-Ganzi-Kekexili block and about O.lmm/a in Kunlun-Qaidam-Qilian block, since the convergence-collision between the Indian continent and Eurasia continent about 50Ma ago. This implies that the thickening rate decreases in the blocks of the Tibetan plateau. From south to north, the displacement of eastern boundary in the Tibetan plateau is about 37mm/a in Himalayan block, about 45mm/a in Lhas-Gangdese block, about 47mm/a in Qiangtang block, about 43mm/a in Songpan-Ganzi-Kekexili block, and about 35mm/a in Kunlun-Qaidam-Qilian block, since the collision-matching between the Indian continent and Eurasia continent had happened about 50Ma ago. This implies that the rate of eastward displacement is biggest in the middle of plateau, and decreases to both sides. The transition of S-N compression stress field in Tibetan Plateau, since about 28Ma+ ago, may be caused by two reasons: On one hand, the movement direction of Eurasia continent changed from northward to southward about 28Ma± ago in the northern plateau. On the other hand, the front belt that is located between India continent's and Eurasia continent's convergence-collision, had moved southward to high Himalayan from Indus-Brahmaputra suture almost at the same time in southern plateau. Affected by the stress field, the earlier tectonics rotated clockwise, NE and NW conjugate strike-slip faults developed, and the SN rift formed. This indicated that the EW movement started. The ratio between upper crust and lower crust of different blocks from south to north in the Tibetan plateau during the process of deformation are as following: about 3.5~5:1 in Himalayan block, about 1~5: 3-4 (which is about 1:3o--4 in south and about 4~5:3 in north) in Lhas-Gangdese block, about 1:3~447mm/a in these blocks: Which is located to the north of Banggong-nujiang suture.