浅析三维动画在水利工程设计中的应用及前景

三维动画是近年来随着计算机软硬件技术的发展而产生的新兴技术,由于其精确性、真实性和无限的可操作性,被广泛应用于医学、教育、军事、娱乐等诸多领域。通过对Bentley软件各模块的运用,对水电站工程水工建筑物及施工区地形、地貌建立可实时显示的仿真模型,使项目人员能够直观、清晰地从任意位置多视角地查看工程场景,对水利工程的设计和施工起到了十分重要的作用。     

应用沉积学数据库进行图解式交流的连接

在现代信息技术中，应用数据库进行交流存在着一个极为重要的问题，即用户的观点差异问题。为了能够应用信息系统进行交流，传统的系统要求用户具有专门的训练。这种专门训练往往是长期的，甚至必须具有良好的计算机基础知识和勤奋好学的性格。因此，无论这些信息系统多么先进和有益，上述原因足以使一些人望而却步，放弃使用计算机信息资源。基于这一问题，贝尔格莱德大学矿业和地质学院开发了一个连接软件，它可弥补上述传统用户交流手段的不足。本文介绍了这一软件的结构，即沉积学数据的录入方法和应用数据库进行图解式交流的连接软件的原理。     

Development of methods for constructing models for intra-year irregularity in the Earth’s rotation

The development of methods for the construction of stochastic, dynamical models for intra-year irregularity of the Earth’s rotation is considered. A correlational model based on harmonically additive and parametrically random, colored and broadband, gravitational-tidal perturbations from the Sun and Moon is developed. One-dimensional and multi-dimensional characteristic functions are found for the case of Gaussian and non-Gaussian colored and broadband fluctuations in the irregularity of the Earth’s rotation. Examples of computer modeling of the irregularity in the Earth’s rotation based on a priori and a posteriori IERS data are presented.     

The computation of deformations

The mathematics of irrotational deformation are simplified by presentation in matrix form. Matrix equations are easily programmed and are easily interpreted in geometrical terms. Graphical operations commonly carried out on orientation nets such as rotation of data can be translated into simple matrix equations for use with a computer. If the shape and orientation of the deformation ellipsoid for a pure shear are known, a matrix can be constructed for use as a deformation matrix. This can be used to deform other deformation ellipsoids to obtain a resultant ellipsoid. It can also be used to deform geological structures such as lineations and planes. The matrix equations for these operations are very simple, but their numerical solution often requires a computer.     

The computation of deformations

The mathematics of irrotational deformation are simplified by presentation in matrix form. Matrix equations are easily programmed and are easily interpreted in geometrical terms. Graphical operations commonly carried out on orientation nets such as rotation of data can be translated into simple matrix equations for use with a computer. If the shape and orientation of the deformation ellipsoid for a pure shear are known, a matrix can be constructed for use as a deformation matrix. This can be used to deform other deformation ellipsoids to obtain a resultant ellipsoid. It can also be used to deform geological structures such as lineations and planes. The matrix equations for these operations are very simple, but their numerical solution often requires a computer.     

Numerical modelling of groundwater vulnerability: the example Namibia

The UNSAT-H, HELP3 and MACRO4.3 computer codes, which have been developed for simulating the water balance of the unsaturated zone of soils and unconsolidated sediments, are also capable simulating water flow in low-porosity media such as fractured rock. The codes can be used to model the ability of rocks and overlying soils in the vadose zone to protect the groundwater in the uppermost aquifer. The net infiltration rates simulated by the different codes are compared against recharge determined by the chloride-balance method in Namibia. The dual-permeability code MACRO4.3 was found to produce more realistic estimates of net infiltration than the UNSAT-H and HELP3 codes, which are based on a single-permeability or effective-continuum method. The net infiltration rate together with the water storage in the unsaturated zone and the groundwater depth are used to calculate the residence time of pore water in the unsaturated zone. This parameter determines the intrinsic vulnerability of the aquifer.     

CAN总线与LPT通信接口设计的与实现

现场总线是用于自动化系统最底层的现场设备或仪表互联的通信网络。CAN是现场总线中的一种,具有较高的性价比和较强的实时处理能力。微机LPT并行口具有增强工作方式(EPP),可以进行双向快速数据传输。用CAN控制器SJA1000设计实现计算机并口LPT与CAN总线通信,介绍了并行口EPP方式的信号定义,寄存器工作时序及利用EPP模式进行数据快速双向传输的方法。给出了通信接口电路的设计方法和CAN总线与并行口EPP模式的通信程序,并对通信程序进行了的调试与验证。     

Principles for selecting a list of reference radio sources for a celestial coordinate system

Time series of the coordinates of radio sources defining the celestial coordinate frame are analyzed. Methods for selecting such sources so as to enhance the stability of the frame are considered. Some of these sources, including so-called “defining” sources, demonstrate significant proper motions. Since the stability of the celestial coordinate frame is determined by an absence of global rotation relative to the defining sources (no net rotation), variation in their coordinates will lead to a rotation of the axes of the celestial coordinate frame. The parameters of this rotation are calculated for two physical models for the motions of extragalactic radio sources. The motions displayed by the first group of sources are linear and uniform. Since the apparent speeds of radio sources are often close to, and sometimes exceed, the speed of light, it is supposed that such radio sources have relativistic jets or plasma clouds that move with speeds roughly equal to the speed of light. The observed uniform, linear motion can then be explained by precession of the jet. The second group of sources display non-linear motions, interpreted as a manifestation of the acceleration of matter by the jet. It is assumed that a cloud of particles that moves into the path of the jet is accelerated to relativistic speeds by the jet. A list of sources that should form a very stable coordinate system for several decades into the future is composed based on these two models.     

基于BP神经网络的土体细观力学参数反演分析

利用离散元方法对颗粒材料的细观力学特性研究, 目前确定数值计算模型的细观力学参数大多数通过反复调试获取, 效率低、可重复性差。本文采用开源的颗粒离散元程序LMGC开展了土体双轴压缩数值试验, 通过25组土体细观力学参数计算得到相应的宏观力学参数, 建立了BP人工神经网络反演系统。利用土体物理试验得到的土体宏观力学参数, 输入BP神经网络, 反演得到土体的细观力学参数。将所得细观力学特性参数输入所建立的土体数值计算模型, 得到土体破坏过程中的应力-应变关系曲线, 以及土体颗粒的力链图和旋转变形云图。所建立的土体数值试验模型能够较好地模拟土体变形破坏过程, 利用BP神经网络反演细观力学参数以及数值模型计算得到的土体宏观力学参数与物理试验吻合较好, 误差在10%左右, 土颗粒间力链云图以及旋转变形云图较好地揭示了土体变形破坏的机理。     

在AutoCAD200O中应用VBA技术自动绘制等高线

应用AutoCAD进行地形图的计算机辅助制图能极大地提高工作效率，但因其没有根据数字地形分析（DTM模型）自动产生等高线的功能，用VBA编程技术，将外业采集的离散点通过建立三角网，在三角形边上内插等高点的方式，能在AutoCAD中实现等高线的自动绘制。     

金则石粒径测试方法的研究

利用数字摄像和计算机技术，对金钢石的粒径等参数进行了测试，对数据的分析揭示了金刚石粒径，等效短径和平均网径之间的内在关系，提出了以平均网径为基础，建立金刚石粒径标准的方法。     

Making ends meet in Gondwana: retracing the transforms of the Indian Ocean and reconnecting continental shear zones

Interpretation of the detailed patterns of ocean-floor transforms revealed by satellite altimetry enables the creation of the Indian Ocean to be described quantitatively as four consecutive plate-tectonic regimes separated at 200, 136, 89 and 43 Ma. Each regime is reversed in turn by keeping transform termini coincident and colinear until conjugate points on the margins of pre-existing plates regain their pre-regime integrity. Progressive elimination of the Indian Ocean, demonstrable as a smooth computer animation ( http://www.kartoweb.itc.nl/gondwana ), leads to a refined re-assembly of the continental fragments of central Gondwana that is substantiated by new geological data. A sequence of Euler interval poles that describes the dispersal of the Gondwana fragments, time-calibrated against available magnetic anomaly data, is given. The model requires a mid-Cretaceous position for India's southern tip about 1000 km south of Madagascar, prior to India's rapid northward migration.     

2008年汶川地震破裂的滑移向量

估计同震滑移向量对于认识和理解破裂方式和破裂过程具有重要意义。2008年汶川大地震在青藏高原东缘龙门山推覆构造带的中央断裂和前山断裂上各形成了一条长250 km和72 km的地表破裂带。地震发生后至今,已经发表了大量有关同震位错沿破裂带分布的论文和报告,但绝大部分都仅仅是破裂的走向位错和垂直位错,极少有同震滑移向量的报道。这不仅是因为野外难以直接测量到水平缩短量(或拉张量),而且还因为这些走滑位错实际上是视走滑位错,部分或全部来自水平缩短或拉张。因此,仅仅根据视走滑同震位错和垂直同震位错估计的同震总滑移量肯定包含了相当大的误差。尝试利用据不同走向参考线测量到的一组(两个以上)视走滑位错来计算水平滑移向量的这一新方法,获得了中央破裂带上的7个水平同震滑移向量,并结合垂直位错量进一步计算了走滑、倾滑和水平缩短三个同震滑移分量以及断层倾角和破裂面上的同震滑移向量,综合出露破裂面的擦痕所指示的滑移向量,并对比根据矩张量解获得的震源深度的滑移向量,得出以下认识:(1)破裂南段的地表滑移向量的方位角明显小于震源深度滑移向量的方位角,表明在破裂从震源向地表传播过程中破裂面上的滑移向量发生了逆时针旋转;(2)滑移方位角向北东方向逐渐增大,表明地平面上水平滑移向量表现出顺时针旋转的趋势,而且在破裂向北东方向传播过程中近地表的走滑分量逐渐减小而倾滑分量逐渐增大;(3)几乎在每一个观测点倾滑分量都大于走滑分量,表明汶川地震的破裂方式在任何地点都是以逆冲运动为主;(4)破裂面倾角在10.4°～64.7°,平均值为41°,与天然破裂露头和探槽揭示的结果基本一致;(5)滑移向量沿破裂带的分布显示,走滑分量中段大而两端小,倾滑分量则相反,中段小两端大。     

Microplate motions in the hotspot reference frame

We have computed motions of the major plates (seven large plates and seven medium‐sized plates) and 38 microplates relative to the hotspot reference frame, and present velocities of these 52 plates. Moreover, using updated plate boundaries for the present, we have computed new geometrical factors for plates and microplates, useful for kinematic calculations and to obtain the net‐rotation of the lithosphere and plate velocities in the mean‐lithosphere reference frame. Instead of a continuum or gradational distribution of the plates by size, the plates clearly partition into three groups each having their own characteristics. For the seven large plates, rotation poles generally lie in high latitudes; the seven medium‐sized plates have rotation poles in a restricted equatorial area; the 38 small plates show the greatest scatter. Moreover subsets of the 52 plates reveal differing fractal behaviour: the large, middle and small groupings each have a characteristic fractal dimension, suggestive of microplate clustering. The highest angular velocities occur for some of the smallest plates, with the location of their rotation poles closeby. Terra Nova, 18, 276–281, 2006     

宝石折射仪测试结果的计算机模拟与分析

本文在分析了宝石折射仪原理的基础上,提出可以用光率体在垂直于折射仪长轴方向的椭圆切面来判断所测宝石折射率的大小。据此原理,建立了分析折射仪测试结果的数学模型,用计算机模拟出了可能存在的30种折射率曲线型式。并对这些曲线型式的特点进行了分析。     

Calculation of finite incremental deformations in ductile geological materials and structural models

Although the total continuous deformation of rocks can be the net result of several tectonic pulses (finite increments), it will be assumed that there are only two. It is further assumed that differential displacements are unknown, but that two generations of rotation gauges are available in a rock body, together with two respective generations of strain gauges. A mathematical procedure is outlined herein for calculating the first or second deformation increment in natural structures and putty models with elliptical markers. To illustrate their practical value, one set of equations is used to determine the pattern of large incremental strain associated with the lateral spreading of a mature model diapir. This model pattern may help to interpret igneous fabric patterns of undeformed granitoid plutons.     

Folding by cataclastic flow at shallow crustal levels in the Canyon Range,Sevier orogenic belt,west-central Utah

Folds form by ductile deformation typically involving continuous flow. In the elastico-frictional regime, such deformation may be accomplished by cataclastic flow involving collective movement on a population of fractures and zones. The Canyon Range (CR) syncline, part of the CR thrust sheet in west-central Utah, developed in this regime. The CR syncline is composed of thick-bedded quartzite units with a small material contrast between layers, limiting limb rotation by flexural slip alone. Thus, fracture populations developed to accommodate fold tightening by limb rotation and thinning, and the formation of transverse zones across the fold. Several generations of fracture and deformation zone (DZ) networks are recognized from mesoscopic and microscopic evidence, and can be related to stages of folding. The net result of the large number of distributed fractures and deformation zones is a continuous deformation that is homogeneous at the scale of the outcrop. All these lines of evidence suggest that large-scale cataclastic flow accommodated folding by allowing rigid mesoscopic blocks to slide along bounding DZs.Along its length, the CR syncline consists of several segments bounded by transverse zones with different mechanisms accommodating fold tightening in adjacent segments. In one segment, fold tightening progressed by limb rotation, and then out-of-the-core thrusting. In contrast, fold tightening in the adjoining segments occurred by rotation and thinning of one limb and possible hinge migration, with the steeply dipping to overturned limb showing progressive thinning of units on a megascopic scale and progressive increase in the thickness and density of deformation zones at all scales.     

The geometry of planar domino-style normal faults above a dipping basal detachment

Geometrical analysis of planar domino-style normal faults rooted into a dipping basal detachment fault allows derivation of equations which relate: (1) horizontal extension within the upper plate; (2) dip of the detachment; (3) final fault dips; (4) rotation that faults and beds undergo; and (5) net slip on domino-style faults. Past geometrical models have focused on extremely idealized and non-unique geometries, in which domino-style faults are parallel and similar cut-off points are all at the same elevation after faulting. This corresponds to evenly spaced domino-style faults above a horizontal detachment. Considering non-parallel faults and dipping detachments introduces unique geometries, which allows calculation, for example, of permissible depth ranges to the detachment. Horizontal extension varies significantly: (1) with dip of the detachment; and (2) for synthetic and antithetic cases, in which the domino-style faults dip in the same or opposite direction, respectively, as the detachment. When other factors (e.g. rotation and fault dips) are held constant, horizontal extension greatly increases for moderate detachment dips and moderately decreases for antithetic detachment dips, as compared to the horizontal detachment case. This is important because the synthetic case has been widely reported.     

Horizontal displacements in Japan

Recent horizontal displacements in Japan are discussed by using the results of the first-order triangulation surveys. Large horizontal displacements are found in East Hokkaido, Tohoku, South Kanto, Tokai and Nankai districts (Fig.1). Crustal activity in these districts is also briefly discussed. The original triangulation survey in Japan neglected Laplace (azimuth) observations and consequently rotation and divergence around the assumed fixed point often appear as horizontal displacements. This is especially true in the central Japan where horizontal displacements result from cancelling the apparent rotation and divergence around the assumed fixed point (Fig.2).On the other hand, strain measurements avoid such shortcoming. Considering with the results of the first triangulation, the horizontal earth-strain can be calculated for every subsequent triangulation net. It is interesting that the velocity of the maximum shear strain is almost 2–3 · 10^/t7/year throughout Japan (Fig.3), even though the seismic activities show large regional discrepancies.     

Neotectonic deformation in the western sector of tectonic escape in Anatolia: palaeomagnetic study of the Afyon region, central Turkey

Following final closure of the Neotethyan Ocean during the late Miocene, deformation in central Turkey has led to crustal thickening and uplift to produce the Anatolian Plateau followed by westward extrusion of terranes by strike–slip. Widespread volcanism has accompanied this latter (neotectonic) phase, and palaeomagnetic study of the volcanism shows a coherent record of differential block rotations, indicating that the Anatolian region is not a plate (or ‘platelet’) sensu stricto but is undergoing distributed internal deformation. To evaluate the scale of neotectonic rotations in the transition zone near the western limit of tectonic escape and the border of the extensional domain in central-west Turkey, we have studied the palaeomagnetism at 82 sites in volcanic suites distributed along a 140-km lineament with north–south trend and ranging in age from 18 to 8 Ma. Comparable deflection of magnetic remanence from the present field direction is identified along the full length of the lineament. A mean clockwise rotation of 12.3±4.2° is determined for this western sector of the Anatolian strike–slip province. Since similar rotations are observed in the youngest and oldest units, this cumulative rotation occurred after the late Miocene. When interpreted together with results elsewhere in Anatolia, it is inferred that the rotation is later than crustal thickening and uplift of the Anatolian Plateau and entirely a facet of the tectonic escape. Inclinations are mostly 10° shallower than the predicted Miocene field and are considered to reflect the presence of a persistent inclination anomaly in the Mediterranean region. Larger rotations departing from the regional trend are also observed within the study region, but are confined to the vicinity of major faults, notably those bounding the Afyon-Ak ehir Graben.The pattern of neotectonic declinations across Anatolia identifies strong anticlockwise rotation in the east near the Arabian pincer with progressive reduction in the amount of rotation towards the west; it becomes zero or slightly clockwise at the western extremity of the accreted terrane collage. Rotations also appear to become generally younger towards the south. Crustal deformation has therefore been distributed, and the net effect of terrane extrusion to the west and south has been to expand the curvature of the Tauride Arc. The westward radial expansion of the extruded terranes is inferred to combine with backroll on the Hellenic Arc to produce the contemporary extensional province in western Turkey.