在CorelDRAW平台上进行Fry法有限应变测量的新技术

有限应变测量是韧性变形定量估算的主要手段。Fry法有限应变测量以变形颗粒,如变形砾石、鲕粒、近等轴矿物等为应变标志体,地质实践中较容易获得,是岩石有限应变测量的最普遍手段。其传统测量过程极为繁琐,难以满足大范围应变测量的要求,数据精度也较低,限制了Fry法的应用。前人曾根据Fry法原理设计了计算机程序,但其程序相关编程语言如今已基本不能使用。本研究基于目前地质上较普遍使用的制图软件CorelDRAW平台的VBA编程,设计出一个基于CorelDRAW平台的Fry法应变测量宏程序FRY-1。在新图层上将所有变形颗粒的中心点标出后,运行FRY-1宏程序即可快速得到精确的Fry图解,并计算出相对精确的R值。通过实例与“最邻近心对心法”的测量结果进行的对比,证明该方法简便易操作,并可大幅提高有限应变测量的效率和有效减小误差。     

利用惯量椭圆进行岩石有限应变分析

提出了利用惯量椭圆对任意形状矿物颗粒进行描述的方法，并利用惯量椭圆理论，计算了岩石薄片中任意形状矿物颗粒的惯量椭圆，通过颗粒面积对椭圆参数进行标准化，得到每一矿物颗粒的等效应变椭圆。等效应变椭圆能够有效地反映对应矿物颗粒的优选方位以及变形特征，进而利用椭圆的矩阵参数形式对等效应变椭圆进行统计分析，获得岩石的有限应变椭圆；同时给出了相应的数值计算方法，编制了软件Straindesk，并得到了成功的应用。该方法克服了先前应变测量中的局限性，方便实现计算机的自动分析，具有较强的有效性和广泛的适用性。     

利用变形花岗岩体中的长石矿物进行有限应变测量初探

岩石有限应变测量的方法较多，相比之下，“长短轴法”最简单、快捷，野外易于操作，但是对所选的对象变形前要求近圆形，或者椭圆形，且随机分布。花岗岩区很难找到满足此条件的颗粒。在《1：25万玉门镇幅》调研的过程中，试以长石颗粒为研究对象，进行长短轴法变形测量，将其结果与相邻变形砾石、杏仁体和包体测量结果进行对比，发现它们具有相似的应变椭圆轴率，证实在变形花岗岩区可以用变形长石颗粒研究区域变形的特点。     

岩石有限应变测量反向轮法的计算机CSD软件设计

用计算机实现反向轮法测岩石有限应变的过程, 可提高有限应变测量的准确性和效率.开发出了用Visual Basic 5.0编写的CSD软件, 利用岩石薄片图像或者矿物颗粒轮廓图测出岩石有限应变的大小及应变椭圆长轴方向, 同时测出矿物颗粒分布的优选方位.其操作相当简便, 首先统计出图像中0°~180°各方向上矿物颗粒的边界数目; 然后在方位-边界坐标系中投点, 并利用最小二乘法进行数据点的多项式曲线拟合; 最后求出曲线的极值点坐标, 并根据坐标绘制相应的应变椭圆.软件运行中, 统计出的边界数据和拟合曲线以及应变椭圆图形都是可视的, 并能进行相应的保存.      

3D-DIC系统在岩石力学试验中的应用

介绍了全场3D-DIC系统,在试样周围均匀布置三组3D采集元,利用散斑图像的采集,尽可能在较大范围内捕捉试样表面的变形并进行分析。为了验证3D-DIC系统测量结果的可靠性,以江持安山岩为试验材料,进行了单轴压缩条件下的恒定应变速率和交替应变速率试验。通过3D-DIC系统的测量结果与应变片测量结果的对比发现,两者吻合度较好,表明3D-DIC系统能够满足岩石力学试验非接触式、可视化变形测量的需求。另外以田下凝灰岩为试验材料,通过图像分析得到了加载过程中岩石表面的全场变形信息,发现岩石变形的不均匀性和应变局部化现象,表明利用3D-DIC系统对岩石破裂过程变形测量比传统变形测量方法更具有优越性,为岩石力学试验研究提供了一种新的、便捷的测试方法。     

华北克拉通基底花岗质片麻岩变形和流变学研究——以辽西寺儿堡地区为例

辽西寺儿堡镇新太古代花岗质片麻岩内发育的宏观、微观构造变形特征表明该地区曾遭受了强烈的韧性变形改造。花岗质岩石变形程度在初糜棱岩–糜棱岩之间,岩石经历了SWW向左行剪切作用改造。岩石中石英有限应变测量判别结果表明,构造岩类型为L-S型,为平面应变。岩石的剪应变平均值为1.43,运动学涡度值为0.788~0.829,指示岩石形成于以简单剪切为主的一般剪切变形中。此外,石英颗粒以亚颗粒旋转重结晶和颗粒边界迁移重结晶作用为主,长石颗粒塑性拉长,部分发生膨凸式重结晶作用;石英组构特征(EBSD)揭示石英以中–高温柱面滑移为主;石英颗粒边界具有明显的分形特征,分形维数值为1.151~1.201,指示了中高温变形条件。综合石英、长石的变形行为、石英组构特征以及分形法Kruhl温度计的判别结果,推断辽西寺儿堡镇新太古代花岗质片麻岩经历过480~600℃的中高温变形,其同构造变质相为高绿片岩相-低角闪岩相。花岗质岩石的古差异应力为10.62~12.21 MPa,估算的应变速率为10~(–11.67)~10~(–13.34) s~(–1),即缓慢的变形,可能记录早期中高温、低应变速率的韧性变形过程,反映华北克拉通基底中下部地壳变形特征。     

骨节状石香肠构造应变计初探

考虑两期纯剪切变形成因的骨节状石香肠。假设骨节缝处香肠层的厚度近似等于香肠层的原始厚度,再根据等面积法原理恢复骨节状石香肠层的原始长度,求得石香肠层的拉伸有限应变,并采用物理模拟方法对骨节状石香肠等面积法进行了检验,误差范围分别为3.23%,2.84%,3.94%和6.82%。将该种方法应用于湖北铁山含骨节状石香肠岩石的有限应变测量,得到石香肠层的伸长有限应变值为33.80%~48.22%。再利用矩形断裂石香肠对骨节状石香肠应变测量结果进行检验,取得比较一致的结果。并将两期构造变形的岩石有限应变进行了分离。     

大理岩常规三轴压缩下强度和变形特性的试验研究

通过在伺服试验机上对中、粗颗粒大理岩进行常规三轴试验,基于试验结果,研究了两种颗粒大理岩的强度和变形特性。结果表明,低围压时岩样内部材料并未均匀化,岩石表现为应变软化特性;而高围压时岩样内部材料强度由低到高逐渐屈服,变形趋于均匀,岩石出现塑性流动特性。岩石的峰值应变与围压显著成正线性关系;中颗粒和粗颗粒岩样之间确实存在差异,但内摩擦系数可以表征材料力学性质,与粒径没有关系;岩石残余强度对围压的敏感性显著高于峰值强度。     

赣西北中元古界双桥山群岩石有限应变分析

对赣西北中元古界双山群岩石进行系统的有限应变测量,结果表明,双桥山群岩石的主要应变型为平面应主为,其对数Flinn参数K值在0．74－1．44之间,若石应变强度T值在2．59－4．60之间,平均为3．15,表明岩石应变强度较大,岩石有限应变强度具有南强北北,东强西弱的规律,有限应变的椭球X轴优势方位为近东西向,反映双桥山群岩石在变形作用过程中,近南北向压应力作用占主导地位。     

应变测量及其地质意义

利用岩石中的某些标志(如鲕、砾石等)分析自然变形岩石中的有限应变,对了解变形岩石构造的几何特性、形成方式及其应变历史,都有很重要的意义.早在19世纪中期,西方一些学者已开始利用变形化石作为应变标志研究地质构造(特别是劈理)的成因.此后,不少学者利用岩石中的应变标志对变形岩石中的构造以及区域的应变历史进行了研究.其中尤以克鲁斯(E.Cloos,1947)利用鲕石的变形研究美国马里兰州南山的岩石变形,是一次系统的尝试.在《岩石的褶皱和断裂》("Folding and fracturing of rocks",1967)一书中,兰姆     

Strain analysis using quartz deformation bands

Quartz deformation bands are kink bands in quartz crystals. A deformation band develops as a region of localized crystal-plastic deformation with boundaries perpendicular to the slip plane and slip direction, which usually is along an -axis in the basal plane. Under cross-polarized light, the difference in crystallographic orientation between a deformation band and its host is indicated by a difference in extinction positions. The displacement between the c axis in a deformation band and the c axis in the host represents the angular shear of the deformation band in the direction of the c axis in the host grain. Assuming the deformation is homogeneous at the grain scale, the angular shear of the grain (the gauge) is calculated by multiplying the angular shear of the deformation band by the ratio of the sheared part to the whole grain. Using the strain-gauge method for three-dimensional infinitesimal strain analysis, a minimum number of five grains measured on universal stage is needed to solve for the deviatoric strain components of the aggregate if the strain is homogeneous in the aggregate. Data from more than five grains are used to find the best-fit strain components by a least-squares method. The principal strains and their orientations are found from these strain components by calculating the eigenvalues and eigenvectors. A 3-D strain ellipsoid also is obtained from strain ellipses in three perpendicular planes determined from the two-dimensional flat-stage measurements by the Wellman method. Both the strain-gauge method and the Wellman method are tested by using synthetic data sets and applied to a naturally deformed sample. Both methods give similar results; the established Wellman method thus confirms the strain-gauge calculation.     

Pure shear to simple shear-dominated transpression during the Eburnean major D2 deformation,Daléma sedimentary basin,eastern Senegal

Field studies in the Palaeoproterozoïc Daléma basin, Kédougou-Kéniéba Inlier, reveal that the main tectonic feature comprises alternating large shear zones relatively well-separated by weakly deformed surrounding rock domains. Analysis of the various structures in relation to this major D₂ phase of Eburnean deformation indicates partitioning of sinistral transpressive deformation between domains of dominant transcurrent and dominant compressive deformation. Foliation is mostly oblique to subvertical and trending 0–30° N, but locally is subhorizontal in some thrust-motion shear zones. Foliation planes of shear zones contain a superimposed subhorizontal stretching lineation which in places cross-cuts a steeply plunging stretching lineation which is clearly expressed in the metasedimentary rocks of weakly deformed surrounding domains. In the weakly deformed domains, the subhorizontal lineation is absent, whereas the oblique to subvertical lineation is more fully developed. Finite strain analyses of samples from surrounding both weakly deformed and shearing domains, using finite strain ratio and the Fry method, indicate flattened ellipsoid fabrics. However, the orientation of the long axis (X) of the finite strain ellipsoid is horizontal in the shear zones and oblique within the weakly deformed domains. Exceptionally, samples from some thrust zones indicate a finite strain ellipsoid in triaxial constriction fabrics with a subhorizontal long axis (X). In addition, the analysis of the strain orientation starting from semi-ductile and brittle structures indicates that a WNE–ESE (130° N to 110° N) orientation of strain shortening axis occurred during the Eburnean D₂ deformation.     

古气候记录的时间标尺及相位关系计算——以塔吉克斯坦Chashmanigar黄土剖面为例

轨道调谐法是确定第四纪地层时间标尺的主要手段,但人工进行轨道调谐确定时间标尺具有可操作性差、客观性差、工作量大及很难得到最佳结果等缺陷,本文描述了利用计算机自动进行轨道调谐以厘定古气候记录时间标尺的方法步骤.对于不同的气候替代性指标时间序列之间在主要周期频率上的相位关系,作者提出了利用演化凝聚函数法进行计算的新方法.作为实例,利用上述方法对塔吉克斯坦Chashmanigar黄土记录进行了时间标尺和相位关系演化的计算.     

The relation between microfabric and strain in a progressively deformed quartzite sequence from Central Australia

In the Ormiston Nappe Complex, west of Alice Springs, central Australia, a deformed zone up to 0.7 km thick is developed in the sedimentary Heavitree Quartzite. The deformed zone is adjacent to a major thrust fault and is defined by mylonitic foliation, which is parallel to the thrust plane and by isoclinal folds. Recognition of original detrital quartz grains allows strain ellipsoids to be measured across the zone. The strain generally plots in the flattening field and many specimens show pure flattening strain. The mylonitic foliation is an axial-plane structure to the folds and is parallel to the XY-plane of the strain ellipsoid. A quartz elongation lineation may be present within the foliation and is parallel to the principal extension direction (the X-axis) of the strain ellipsoid.Strain is accommodated principally by intracrystalline plastic deformation of the quartz grains. In detail the strain is not homogeneous and may vary even between adjacent grains of the same specimen. Quartz optic axis fabrics reflect this strain inhomogeneity. If the strain ellipsoid is an oblate spheroid, c-axes lie in small-circle girdles about the principal shortening axis (the Z-axis). With general triaxial strain, the c-axes lie in a great-circle girdle or girdles which intersect the foliation parallel to the intermediate strain axis (the Y-axis) and lie symmetrically about the Z-axis. A random population of grains from a specimen often shows a composite c-axis pattern between these two types.With approach to the thrust there is an increase in the amount of strain within the specimens. The increasing strain correlates with an increase in the degree of c-axis preferred orientation of the deformed detrital grains and in the amount of new strain-free grains present in the deformed quartzite. Adjacent to the thrust the quartzite is completely composed of polygonal new grains. The new grains probably formed under syntectonic conditions caused by movement along the thrust. The bulk of the new grains developed by increasing misorientation between the subareas of an initially polygonized old grain. There is no evidence of any marked host control on new-grain orientation, but new grain c-axis plots are generally similar to the corresponding old-grain plots from the same specimen.     

Deformed graptolites, finite strain and volume loss during cleavage formation in rocks of the taconic slate belt, New York and Vermont, U.S.A.

An analysis of graptolites in the Taconic slate belt of eastern New York and western Vermont (U.S.A.) shows that they are nearly ideal strain markers. For the three species used in this study, Orthograptus whitfieldii, Orthograptus calcaratus and Climacograptus bicornis, the spacing of thecae is constant except for the first five or so thecae in the proximal part of the fossil rhabdosome. Further, the thecal apertures are perpendicular to the long axis of the stipe. Observations of the thecal spacing in deformed rocks leads to a determination of extension (e=(l_f−l₀)/l₀) and a measure of the angle between thecal aperture and stipe axis yields a direct determination of angular shear strain. In practice, we find it is most straightforward to use length changes to determine the magnitude of principal strains. In Taconic slates, e₁ ranges from 1.0 to 0.24, e₂ ranges from 0.23 to −0.43 and e₃ ranges from −0.56 to −0.74. Thus, we find that the absolute finite strain in these slates is constrictional at three sites, plane strain at another and a true flattening at only one site. An examination of volume changes based on strain results in determinations of between 81% volume loss and 7% volume gain, with volume losses between 28% and 81% in 9 out of 10 calculations. These conclusions are in accord with previous determinations of volume loss based on reduction spot analyses and are consistent with the observation that pressure dissolution was a common grain scale deformation process in cleavage formation but that these slates lack abundant veins, fibrous overgrowths or other identifiable sites of reprecipitation.     

Strain analysis from objects with a random distribution: A generalized center-to-center method

Existing methods of strain analysis such as the center-to-center method and the Fry method estimate strain from the spatial relationship between point objects in the deformed state. They assume a truncated Poisson distribution of point objects in the pre-deformed state. Significant deviations occur in nature and diffuse the central vacancy in a Fry plot, limiting the its effectiveness as a strain gauge. Therefore, a generalized center-to-center method is proposed to deal with point objects with the more general Poisson distribution, where the method outcomes do not depend on an analysis of a graphical central vacancy. This new method relies upon the probability mass function for the Poisson distribution, and adopts the maximum likelihood function method to solve for strain. The feasibility of the method is demonstrated by applying it to artificial data sets generated for known strains. Further analysis of these sets by use of the bootstrap method shows that the accuracy of the strain estimate has a strong tendency to increase either with point number or with the inclusion of more pre-deformation nearest neighbors. A poorly sorted, well packed, deformed conglomerate is analyzed, yielding strain estimate similar to the vector mean of the major axis directions of pebbles and the harmonic mean of their axial ratios from a shape-based strain determination method. These outcomes support the applicability of the new method to the analysis of deformed rocks with appropriate strain markers.     

A phenomenological numerical approach for investigating grain size evolution in ductiley deforming rocks

The sizes of recrystallised grains in exhumed ductile shear zones are often used to infer conditions of deformation (i.e. stress, strain rate and temperature). Here we present a simple numerical method of calculating the dynamic evolution of grain size during ductile deformation. Our phenomenological method is based on the fact that the dynamic competition between grain growth and recrystallisation will drive grains towards a steady-state size. At each time increment, grain growth and reduction contributions are calculated, with magnitudes which depend on the difference between the current grain size and a desired steady-state grain size. In our models we use a recrystallised grain size piezometer to calculate the steady-state grain size for a given stress. Our numerical routine is incorporated into the SULEC finite element package, allowing us to explore spatial and temporal changes in grain size.As a test, we compare model results to measured grain sizes in quartz layers thinned and recrystallised around rigid garnet porphyroclasts under simple shear dominated deformation in the Alpine Fault Zone of New Zealand. Numerical models are able to replicate observed grain size variations, with boundary conditions consistent with those constrained for the central Alpine Fault Zone.     

Implicit and explicit procedures for the yield vertex non-coaxial theory

The yield vertex non-coaxial model is different from classical elastoplastic models, in that there is an additional plastic strain rate tangential to yield surfaces, as well as the plastic strain rate normal to yield surfaces, when orientations of principal stress change. This feature raises concerns on its finite element implementations. In nonlinear finite element numerical iterations, a large tangential plastic strain rate is likely to make the trial total strain rate direct inside a yield surface, which entails convergence difficulty. Some modifications are introduced on the non-coaxial model itself to make numerical convergence easier in the work published in Yang and Yu (2010) [20]. This paper is an extension of the previous work. Instead of modifying the non-coaxial model itself, this paper concerns the use of finite element explicit procedure, which is suitable for highly discontinuous problems. The simulations of shallow foundation load-settlement responses indicate that the finite element explicit procedure, assisted with a robust and explicit automatic substepping integration scheme of the non-coaxial model, does not encounter numerical difficulty. In addition, the overall trends of implicit and explicit simulations are similar.