Three-dimensional topology optimization for geotechnical foundations in granular soil

This article presents three-dimensional structural optimization in geotechnical engineering for foundations in granular soil. The general design (topology) of a shallow foundation is optimized with respect to its deformational behaviour within the service limit state. The SIMP (solid isotropic material with penalization) method is applied to optimize the distribution of foundation material. The soil is modelled as a hypoplastic material with a constitutive model suitable for optimization using finite element analysis. Two load cases are examined. The optimized topology is validated against two-dimensional optimization and 1g-model test results. The present study proves the applicability and shows the potential of topology optimization in geotechnical engineering.     

3D seismic analysis of gravity-driven and basement influenced normal fault growth in the deepwater Otway Basin,Australia

We use three-dimensional (3D) seismic reflection data to analyse the structural style and growth of a normal fault array located at the present-day shelf-edge break and into the deepwater province of the Otway Basin, southern Australia. The Otway Basin is a Late Jurassic to Cenozoic, rift-to-passive margin basin. The seismic reflection data images a NW-SE (128–308) striking, normal fault array, located within Upper Cretaceous clastic sediments and which consists of ten fault segments. The fault array contains two hard-linked fault assemblages, separated by only 2 km in the dip direction. The gravity-driven, down-dip fault assemblage is entirely contained within the 3D seismic survey, is located over a basement plateau and displays growth commencing and terminating during the Campanian-Maastrichtian, with up to 1.45 km of accumulated throw (vertical displacement). The up-dip normal fault assemblage penetrates deeper than the base of the seismic survey, but is interpreted to be partially linked along strike at depth to major basement-involved normal faults that can be observed on regional 2D seismic lines. This fault assemblage displays growth initiating in the Turonian-Santonian and has accumulated up to 1.74 km of throw.Our detailed analysis of the 3D seismic data constraints post-Cenomanian fault growth of both fault assemblages into four evolutionary stages: [1] Turonian-Santonian basement reactivation during crustal extension between Australia and Antarctica. This either caused the upward propagation of basement-involved normal faults or the nucleation of a vertically isolated normal fault array in shallow cover sediments directly above the reactivated basement-involved faults; [2] continued Campanian-Maastrichtian crustal extension and sediment loading eventually created gravitational instability on the basement plateau, nucleating a second, vertically isolated normal fault array in the cover sediments; [3] eventual hard-linkage of fault segments in both fault arrays to form two along-strike, NW-SE striking fault assemblages, and; [4] termination of fault growth in the latest Maastrichtian. We document high variability of throw along-strike and down-dip for both fault assemblages, thereby providing evidence for lateral and vertical segment linkage. Our results highlight the complexities involved in the growth of both gravity-driven normal fault arrays (such as those present in the Niger Delta and Gulf of Mexico) and basement-linked normal fault arrays (such as those present in the North Sea and Suez Rift) with the interaction of an underlying and reactivating basement framework. This study provides an excellent example of spatial variability in growth of two normal fault assemblages over relatively short spatial scales (∼2 km separation down-dip).     

Influence of syn-sedimentary faults on orogenic structures in a collisional belt: Insights from the inner zone of the Northern Apennines (Italy)

This paper discusses the possible influence of syn-sedimentary structures on the development of orogenic structures during positive tectonic inversion in the inner Northern Apennines (Italy). Examples from key areas located in southern Tuscany provided original cartographic, structural and kinematics data for Late Oligocene-Early Miocene thrusts, organized in duplex systems, verging in the opposite direction of the foreland propagation (back-thrusts), which affected the Late Triassic-Oligocene sedimentary succession of the Tuscan Domain, previously affected by pre-orogenic structures. These latter consist of mesoscopic-to cartographic-scale Jurassic syn-sedimentary normal faults and extensional structures, which gave rise to effective stratigraphic lateral variation and mechanical heterogeneities. Structural analysis of both syn-sedimentary faults and back-thrusts were therefore compared in order to discuss the possible role of the pre-existing anisotropies in influencing the evolution of the back-thrusts. As a result, it can be reasonably proposed that back-thrusts trajectories and stacking pattern were controlled by relevant syn-sedimentary normal faults; these latter were reactivated, in some cases, if properly oriented. Such an issue adds new inputs for discussing the potential role of structural inheritance during tectonic inversions, and helps to better understand the processes suitable for the development of back-thrusts in the inner zones of orogenic belts, as it is the case of the inner Northern Apennines.     

Parasitic folds with wrong vergence: How pre-existing geometrical asymmetries can be inherited during multilayer buckle folding

Parasitic folds are typical structures in geological multilayer folds; they are characterized by a small wavelength and are situated within folds with larger wavelength. Parasitic folds exhibit a characteristic asymmetry (or vergence) reflecting their structural relationship to the larger-scale fold. Here we investigate if a pre-existing geometrical asymmetry (e.g., from sedimentary structures or folds from a previous tectonic event) can be inherited during buckle folding to form parasitic folds with wrong vergence. We conduct 2D finite-element simulations of multilayer folding using Newtonian materials. The applied model setup comprises a thin layer exhibiting the pre-existing geometrical asymmetry sandwiched between two thicker layers, all intercalated with a lower-viscosity matrix and subjected to layer-parallel shortening. When the two outer thick layers buckle and amplify, two processes work against the asymmetry: layer-perpendicular flattening between the two thick layers and the rotational component of flexural flow folding. Both processes promote de-amplification and unfolding of the pre-existing asymmetry. We discuss how the efficiency of de-amplification is controlled by the larger-scale fold amplification and conclude that pre-existing asymmetries that are open and/or exhibit low amplitude are prone to de-amplification and may disappear during buckling of the multilayer system. Large-amplitude and/or tight to isoclinal folds may be inherited and develop type 3 fold interference patterns.     

滇西高黎贡变质岩带南段淡色花岗岩脉年代学特征及构造意义

高黎贡构造带作为藏东三江地区大地构造上的一条重要构造带,在印度-欧亚板块碰撞过程中起着重要的调节作用,并被认为是该地区发生"大陆逃逸"与地块旋转的西边界。本文以高黎贡构造带中的深变质岩为研究对象,通过对高黎贡变质岩带南段进行详细地构造解析,发现其内部广泛发育一系列或顺岩石面理、或切割面理分布的淡色花岗岩脉;并选取了5个淡色花岗岩脉样品进行锆石U-Pb测年,这5个样品的锆石~(206)Pb/~(238)U年龄值主要分布于二个区间:68~40Ma和24~21Ma,并在45Ma左右形成一个巨大的峰值。基于锆石的晶形、晶内结构和微量元素组成特征等方面进行锆石成因学分析,认为本文所测试的5个淡色花岗岩脉样品的锆石为深熔作用过程中形成的锆石;从而推测高黎贡变质岩带南段发生深熔作用并形成淡色花岗岩脉的时代主要集中于68~40Ma和24~21Ma。     

基于遥感技术的汪清县森林资源动态监测研究与实践

以吉林省汪清县为研究区,根据地理国情普查成果应用及监测工作通知要求,探讨了遥感技术在汪清县森林资源监测应用上的主要内容和方法。利用吉林省地理国情普查、基础测绘成果数据,结合Landsat、SPOT、ALOS等卫星影像数据,进行汪清县土地覆盖信息提取及动态监测、森林类型提取及动态监测、植被覆盖度提取及动态监测、森林参数模型反演及动态监测、驱动因子选取等内容进行了详细地说明,同时,针对汪清县森林资源各项的动态变化情况进行了驱动力分析。     

GPS/BDS卫星摄动力数值分析

由精密星历利用拉格朗日插值公式求二次导数的方法计算了卫星在J2000.0惯性坐标系下的总加速度;利用现有的力模型计算了地球中心引力,地球非球形摄动力,太阳、月球和其他行星的摄动力,地球固体潮摄动力,相对论效应摄动力对GPS/BDS卫星所产生的加速度数值大小;利用G-file里的BERNE太阳光压模型参数计算了GPS卫星太阳光压摄动加速度大小;对GPS/BDS卫星所受的不同摄动力进行了数值分析,对同一摄动力对不同类型卫星的影响进行了数值分析比较。结果表明,现有力模型与GPS/BDS卫星所受的实际作用力仍有一定的差距,不同类型卫星所受摄动力有明显差异,在精密定轨的实际应用中应根据不同类型卫星建立合适的力学模型。     

基于蜂群算法的城市土地利用变化建模

通过引入人工蜂群算法用于构建土地利用变化的驱动力模型,分析土地利用变化的驱动力机制。算法原理通过模仿蜜蜂采蜜行为,自动搜索和提取土地利用变化样本中不同土地变化类型所对应的驱动力分类规则。分类规则的构建采用“IF…THEN”形式,并选取3种不同的适应度函数分别进行模拟验证。研究案例基于UCI实验数据集和美国纽卡斯尔市真实土地利用变化数据集。由实验结果可知,采用蜂群算法模型的总体精度和Kappa系数评价优于其它算法,表明蜂群算法应用于土地利用变化建模具有可行性。     

三维物体上二阶波浪绕射的实时模拟

采用二阶时域理论对非线性波浪在任意三维物体周围的绕射问题进行了研究，对自由表面边界条件进行Taylor级数展开，应用摄动展开可以建立相应的边值问题，而且此边值问题的计算域不随时间变化，运用基于B－样条的边界元方法求解每一时刻的波浪场，二阶自由表面边界条件在时间上进行数值积分，在自由表面加了一个人工阻尼层以避免波浪的反射，速度势分解为已知的入射势和未知的散射势，初始条件采用二阶Stokes波浪场，通过加入物体表面边界条件，得到散射势在时间和空间上的发展，本文对圆柱所受规则波的二阶波浪力和波浪爬高进行了计算，数值结果表明此理论计算准确，效率高，数值稳定。     

Practical spreading of driving forces for wave-driven currents

Profiles of wave-induced current on plane beaches are generally smooth. Previously, the smooth nature was often explained by involving lateral mixing concept. This paper describes a new spreading approach instead of the previous lateral mixing approach. The spreading of the driving forces leads to similar results on the velocity profiles of the previous approach. The results of the new method show good agreement with some laboratory measurements for a plane beach. The present method was then extended for two-dimensional plan problems. The present approach was tested against a laboratory experiment, and the computed flow field agrees reasonably well with the measured flow field. A merit of the present method is that the computational time increment is not restricted by the lateral mixing.