INVERSE DDA OF A SYSTEM OF BLOCKS IN SPHERE COORDINATES

1　IntroductionDiscontinuousDeformationAnalysis(DDA)developedbyShiin 1 980’s[1 ,2 ] maybeusedtocalculatethedeformationanddisplacementsinamulti blocksystemandtoanalyzetherelationbe tweenforcesanddisplacementsinthesystem .Ifthesurveydisplacementofeverypointordirectionisknown,alldisplacementsanddeformationsofblocks,eventheglide ,theclosureandthestretchoftheboundary planeofblocks ,canbecalculated .TheresultisoptimallyfittedaccordingtotheLeastSquarePrinciple.DDAismainlyusedinrockblockssystems…     

井中激电地-井方式井旁球体正反演

以激发极化法中的地-井方式中的体极化球形体作为研究对象,采用解析法对地下半空间的电位场进行求解,利用提取出的井中二次场异常电位差进行反演,得出极化球体的球心离井的距离和球心的埋深。通过计算机进行正反演模拟,验证了进行正反演计算的解析式的正确性,并从物理意义上对正反演过程中出现的现象进行分析。最后给出该方法的适用条件以及注意事项。研究对激发极化法地下勘探目标的定量分析具有明确意义。     

不等电极距情况下高密度电阻率法的反演结果分析

通过高密度电阻率法的模型实验,得出不等电极距情况对直立高阻薄板和低阻球体反演结果的影响:缺失电极对它正下方的ρ_s影响最大,而异常形态主要在水平方向有明显的变化,异常体随电极缺失的位置有趋势变化,电极缺失的数目越多,异常形态变化越大,但是从整体来看反演剖面图基本相似。     

1966年以来华北地区一系列七级大震破裂过程的数值模拟

用非连续变形分析方法(DDA+FEM)数值模拟，在华北地区各地块相互制约的块体系统环境中，地块边界断层上发生1966年邢台地震、 1969年渤海地震、 1975年海城地震、 1976年唐山地震等不同类型破裂模式大震的破裂过程. 数值模拟结果给出每个大震释放的主应力场，最大剪应力变化等值线图，地震前后位移变化矢量图，发震断层滑移随时间变化以及走滑错距和应力降等震源参数. 这些结果与地震的震源机制，用地震波资料研究得到的震源参数,宏观等震线,地表观测的水平位移矢量图基本一致. 其中1969年渤海地震正交破裂模式的结果与宏观等震线及小震分布图像更接近. 1976年唐山地震复杂震源模式与该震早期余震分布图像更相符. 表明用DDA+FEM方法进行数值模拟研究，能较好地模拟地震破裂过程.     

Applying disk-based discontinuous deformation analysis (DDA) to simulate Donghekou landslide triggered by the Wenchuan earthquake

In this paper, for the first time, disk-based discontinuous deformation analysis (DDA) is applied to simulate a real landslide triggered by an earthquake. For this purpose, the kinematic behaviour of the Donghekou landslide triggered by the Wenchuan earthquake is simulated and the results obtained using disk-based DDA are compared with actual data. The comparisons show that there is a good agreement between the results obtained using disk-based DDA and observed data. The simulation results provided an understanding of the failure behaviour and temporal evolution of the landslide. This study shows that disk-based DDA is a practical numerical tool that can be used to simulate the post-failure behaviour of landslides triggered by an earthquake.     

卫星地面站雷电直击效应防护设计探讨

卫星地面站雷电防护工程通常采用接闪杆作为天线主体的雷电直击效应防护措施,接闪杆的高度不仅影响其保护范围,而且影响其截闪概率。通过计算年预计雷击次数对截闪概率进行量化分析,结果表明,截闪概率近似与接闪杆高度平方成正比,采用过高的接闪杆将增大地面站遭受雷击电磁脉冲损坏的风险。为优化接闪杆设计,定义了保护体积的概念,并利用折线法与滚球法分别计算了三维立体空间内的保护范围。结果表明:接闪杆的保护范围与高度呈现非线性相关,当接闪杆超过一定高度(折线法超过30m,滚球法超过0.8倍滚球半径)后,对保护范围的影响十分有限;当接闪杆高度低于0.4倍滚球半径时,滚球法保护范围较大,反之则折线法保护范围较大。对接闪杆接闪瞬间周边的磁场强度分析结果表明,无屏蔽环境下地面站电子系统与接闪杆的常规距离远小于两者的理论安全距离,实际工程中难以通过增大接闪杆与卫星地面站的间距消除雷击电磁脉冲危害。为降低这一风险,卫星地面站直击雷防护应优先采用天线自带接闪杆的方式,条件不具备时也应尽量避免采用单支高大接闪杆,可选取适当的计算方法,采用多支较低接闪杆共同防护的方案。     

Investigation of highly efficient algorithms for solving linear equations in the discontinuous deformation analysis method

Large‐scale engineering computing using the discontinuous deformation analysis (DDA) method is time‐consuming, which hinders the application of the DDA method. The simulation result of a typical numerical example indicates that the linear equation solver is a key factor that affects the efficiency of the DDA method. In this paper, highly efficient algorithms for solving linear equations are investigated, and two modifications of the DDA programme are presented. The first modification is a linear equation solver with high efficiency. The block Jacobi (BJ) iterative method and the block conjugate gradient with Jacobi pre‐processing (Jacobi‐PCG) iterative method are introduced, and the key operations are detailed, including the matrix‐vector product and the diagonal matrix inversion. Another modification consists of a parallel linear equation solver, which is separately constructed based on the multi‐thread and CPU‐GPU heterogeneous platforms with OpenMP and CUDA, respectively. The simulation results from several numerical examples using the modified DDA programme demonstrate that the Jacobi‐PCG is a better iterative method for large‐scale engineering computing and that adoptive parallel strategies can greatly enhance computational efficiency. Copyright © 2015 John Wiley & Sons, Ltd.     

变比例尺地图投影系统

本文对编制城市游览图提出了变比例尺地图投影系统。通过由普通城市平面图向辅助球面作逆投影A,再由辅助球面向平面作非A投影,构成了变比例尺地图的数学基础。由不同性质的投影的组合,能起到适应不同城市街区结构的特点。使本系统具有相当的灵活性。文中还讨论了辅助球适宜的大小和不同方位投影之间的变换公式。文末试作了北京市的变比例尺地图。     

一种利用积分球原理的日照计新方案

该文提出了一种利用积分球原理设计的新型日照传感器。该传感器由限角挡环、积分球体和光电探测器组成。积分球为一有5个针孔的空心正球体。针孔可阻挡散射光进入球内，而射入球内的直射光，经内壁多次漫反射后，其亮度呈均匀分布。与日照阈值下相应的亮度相比较可确定日照之有无。该日照传感器具有无转动部件、结构简单、无需调整及便于操作等优点。