模拟弹塑性碰撞的修正SPH方法

从基本的无网格光滑粒子法SPH(Smoothed Particle Hydrodynamics)近似出发,修正了模拟固体力学中大变形弹塑性碰撞的SPH方法.在边界处采用修正的边界条件,弹塑性分析过程中采用增量理论计算应力,迭代过程中用守恒光滑法进行滤波修正消除拉力不稳定.对SPH方法进行了程序实现,给出了杆弹塑性碰撞的算例.计算分析表明,SPH方法节点的影响域较大、精度较相同节点间距有限元法的结果有一定差距,但是通过增加粒子数量可以提高SPH的精度,保持了其简单性和计算大变形的特性.      

空间碎片超高速碰撞数值模拟的SPH方法

利用光滑质点动力学SPH(Smoothed Particle Hydrodynamics)方法对Whipple防护结构在空间碎片超高速碰撞下的物理过程进行了数值模拟.在数值模拟中,为了充分发挥SPH方法和有限元方法FEM(Finite Element Method)的优点,利用有限元单元和SPH节点混合建模,将有限元单元和SPH节点(SPH nodes)通过定义接触条件相结合,在大变形和飞溅区域采用SPH节点建模,而小变形区域则采用有限元单元建模,从而大大节省求解时间,提高计算效率.计算结果表明,弹丸在穿透前板后,形成二次碎片,碎片云经膨胀和拉长,对后板造成轻微的损伤,这和文献的相关试验数据是符合的.利用SPH方法对空间碎片的超高速碰撞过程进行数值模拟,不仅很好地预测了Whipple防护结构的破坏情况,而且对整个碰撞过程,包括碎片云的形成、膨胀和拉长过程都有形象的描述,符合超高速碰撞的试验测试结果.     

球坐标系六片网格下三维定态行星际太阳风模拟

采用二阶MacCormack差分格式, 利用稳态的磁流体(MHD)方程组在球坐标系六片网格下模拟研究了行星际太阳风. 六片网格系统能有效避免极区奇性和网格收敛性. 迭代按径向方向推进求解, 很大程度上减少了计算量, 节约了计算时间. 内边界条件根据太阳与行星际观测确定, 比较测试了5种内边界条件, 模拟给出了1922卡林顿周的背景太阳风结构. 几种内边界条件所得模拟结果与行星际观测基本吻合. 太阳风速度采用McGregor 等的经验公式给出, 磁场由水平电流片(HCCS)模型得到, 密度和温度分别根据动量守恒和气压守恒得到, 研究表明采用这样的边界条件模拟结果最佳.      

试验测定DFR的升降法方法

提出了试验测定结构件细节疲劳额定值(DFR)的升降法方法,可以充分估计并节约总试验时间.通过数值模拟比较了升降法数据处理的半子样极大似然法和马尔可夫方法.为满足可靠度要求,对半子样极大似然法估计的疲劳强度标准差进行了修正,并给出了保证95％可靠度的疲劳强度的置信度要求的最少试件数.最后给出了升降法确定DFR的表达式.      

蜂窝夹芯板超高速碰撞仿真

蜂窝夹芯板是一种多用于航天器的特殊结构,采用光滑粒子流体动力学算法(SP H)对蜂窝夹芯板进行超高速碰撞研究时主要存在2个问题:蜂窝芯的几何建模过程复杂;蜂 窝芯层非常薄,为保证仿真精度采用相等的粒子间距必然生成超过计算容忍的粒子数量.结 合参数化程序设计语言(APDL)、用户界面设计语言(UIDL)进行结构建模,编制了APDL生成蜂 窝芯的程序,采用UIDL将该程序集成到ANSYS交互式界面中,能够方便快捷的生成任意尺度 的正六边形蜂窝芯.针对此方法得到的模型,采用SPH单元与壳单元相结合的方法,其中蜂 窝芯采用壳单元,入射体、蜂窝夹芯板的前、后面板采用SPH单元,对蜂窝夹芯板进行了超 高速碰撞仿真研究.仿真结果与相关试验进行了对比,仿真得到的后板破口尺寸与试验比较 接近,证实了该方法适用于蜂窝夹芯板的超高速仿真.      

基于SPH方法的球形贮箱液体晃动动力学分析

针对发生液体破碎现象的航天器球形贮箱液体晃动问题,基于光滑粒子流体动力学方法（SPH方法）对N S方程进行离散近似，分别对核函数、应力张量、边界问题、邻域粒子搜索等方面进行选择确定或计算处理，最终得到基于SPH方法的流体动力学模型。通过建立包括不同燃料液体体积及不同挡板安装位置的仿真模型，分析了航天器球形贮箱的受力及局部压强。通过分析贮箱受力情况及设置测量点压强情况，得到运动参数对晃动的影响规律，为航天器球形贮箱的设计提供了理论依据。     

非结构三角形网格的一个ENO型有限体积方法

基于双曲型守恒律方程,对非结构三角形网格给出了一种ENO(Essentially Nonoscillatory Scheme)型有限体积格式,方法的主要思想是先对每一个三角形单元构造一个加权的二次插值多项式,而在计算交界面的流通量时采用了两点高斯积分公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法.最后给出了该格式收敛的数值阶,并对前台阶问题进行了计算.      

辛算法的纠飘研究

辛算法较RK(Runge-Kutta)方法,保持辛结构不变或保持哈密顿系统规律性不变是突出的优点,但点态数值精度并不理想.推导出了三阶、四阶辛算法的漂移量计算公式,通过补偿漂移量就能提高点态数值精度,既保辛结构又保证点态数值高精度,适合于工程应用.建立了分数步对称辛算法(即FSJS算法)的纠漂公式,制定了漂移的约束标准.相关算例的数值结果表明:三阶FSJS算法漂移量最小,点态数值精度更高.     

磁层顶通量传输事件的经验重构

在地球磁层顶附近观测到的通量传输事件（Flux Transfer Event,FTE）一般被认为是瞬态局域磁重联的产物,是太阳风质量、动量和能量进入地球内磁层的重要通道.重构FTE的磁场结构可促进对其形成、演化过程及其与周围等离子体环境相互作用的理解.Grad-Shafranov重构法和磁通量绳拟合法等传统磁场重构方法适用于满足特定物理条件的磁场结构.基于平面线性插值原理,设计了一种不受具体物理条件限定的二维FTE磁场结构重构法.模型测试以及对THEMIS和Cluster卫星簇分别观测到的两个FTE的实际应用表明,在合适的多卫星位形条件下,该方法能快速有效重构出FTE的磁场空间分布,有助于推测FTE的磁场线位形,理解卫星测量数据的时间变化,以及分析等离子体物理量相对于FTE的磁场空间分布特征.      

交互式多模型粒子滤波优化重采样算法

针对标准交互式多模型粒子滤波(IMMPF)算法中存在粒子退化及多样性匮乏问题,提出了交互式多模型粒子滤波优化重采样(IMMPFOR)算法,利用线性优化理论改善模型中具有小权值的粒子精度。该算法的新颖性体现在给定量测信息条件下,利用线性优化方法及模型交互概率将每个模型中拥有小权值的粒子替换成新的粒子。新的粒子既包含本模型中粒子信息,又包含了本模型与其他模型交互后的粒子信息。目标跟踪的仿真结果证明:每个模型新产生的粒子集合可以准确地近似真实状态后验概率密度函数,系统的估计精度与标准IMMPF算法相比有较大提升。     

Comparison of Ionospheric Total Electron Content (TEC) over Sonmiani (Pakistan) with NeQuick-2 and IRI-2012 during July 2014 – June 2015

In this study, Total Electron Content (TEC) observations acquired by a GNSS receiver installed at Sonmiani (Geog. Coord. 25.19°N, 66.74°E, Geomag. Coord. 17.62°N, 141.5°E) are being reported for the first time. The data utilized is hourly instantaneous TEC values during 10 International Quiet Days (IQDs) per month from Jul-14 to Jun-15, totaling 120 observation days for monitoring nominal TEC. The findings confirm the semi-annual trend of TEC over Sonmiani, which lies at the northern crest of Equatorial Ionization Anomaly (EIA) region. The TEC measurements are then compared with NeQuick-2 and International Reference Ionosphere (IRI-2012) models. It was found that the TEC values derived from NeQuick-2 are in better agreement with GNSS measurements than those from IRI-2012. The TEC measurements also show seasonal variation which is largest during Equinox months. The TEC value in Dec solstice is higher than the Jun solstice, which confirms that the seasonal anomaly is playing a major role in this region during the course of study.     

chirp-scaling算法中运动补偿问题研究

由于避免了插值计算,在大前视角和大范围成像时,chirp-scaling算法(CSA)的性能优于距离-多普勒算法(RDA).自聚焦算法在SAR成像处理中通常是必不可少的,但是现有的自聚焦算法几乎都是与RDA相结合的.本文提出了能够将自聚焦算法与CSA相结合的方法,使得CSA更具有实用性.本文同时也提出了将基于运动传感器的运动补偿和基于雷达数据的自聚焦算法与CSA结合的方法.该方法非常适合于高分辨率机载SAR成像处理.      

NVNA谐波相位检验件设计及定标

本文介绍了一种用于检验非线性网络分析仪谐波相位指标的检验件设计方法，并搭建了基于实时示波器的测量系统对谐波相位进行定标。通过高速实时示波器对检验件的输出信号进行采样，利用傅立叶变换提取谐波相位信息。为了得到高准确度的测量结果，采用去缠绕、最小二乘等算法对相位结果进行优化处理。实验结果表明，该方法可以作为NVNA谐波相位指标的验证手段对NVNA的谐波相位进行验证。     

GPS多接收机组合的共视时间比对技术研究

为提高GPS共视时间比对的性能，本文介绍了利用中科院国家授时中心（NTSC）和德国物理技术研究院（PTB）两个守时实验室的多台接收机，基于多接收机组合技术原理，构建了NTSC和PTB各自的多接收机组合系统，采用数据融合技术，计算得出两个多接收机系统的观测数据并进行了时间比对试验及性能分析。结果表明，多接收机组合不仅提高了GPS共视的可靠性和稳定性，且提高了时间比对链路的精度。多接收机链路共视时间比对结果的标准偏差STDEV为1.36ns，比单接收机链路时间比对结果的STDEV值平均提高了19.4%，日稳可达3.2×10^-14。     

同步三星目标运动状态快速检测方法

针对星载无源定位系统难以从定位结果及时准确地检测目标运动状态的问题,利用同步三星时差定位体制的定位误差特性,根据静止假设下同步双星时/频差与三星定位结果应相匹配的原理,提出了基于定位一致性的目标运动状态快速检测算法.先通过定位结果的误差分析掌握目标实际位置所在区域,然后在静止假设下统计该区域的频差范围,并将定位时刻的频差观测值与之比较,实现目标运动状态检测,并且推导了检测概率的理论表达式.仿真和试验结果验证了算法的有效性,时/频差估计误差标准差为10 μs/100 mHz时,对速度大于40 m/s的目标,检测概率高于99%.实测数据表明该算法显著提高了系统目标运动状态检测的准确性和时效性.     

基于主元分析和BP神经网络对刀具VB值预测

对声发射信号进行5层小波分解提取6个频段的能量值,把它与切削速度、切削深度、进给量和切削时间一起作为刀具状态的特征向量.通过主元分析进行降维、消除特征向量间的相关性后,把得到的主元作为BP(Back Propagation)神经网络的输入向量.BP神经网络应用改进的LM(Levenberg-Marquart)算法进行学习,利用输入向量对网络进行训练后,实现对刀具后刀面磨损量VB的预测.实验结果显示:基于主元分析和LM算法改进的BP神经网络建立的预测系统,网络输出与实测VB值的误差0.03以内;根据预测VB值的范围可判别出刀具的不同状态.      

Near real-time PPP-based monitoring of the ionosphere using dual-frequency GPS/BDS/Galileo data

Ionosphere delay is very important to GNSS observations, since it is one of the main error sources which have to be mitigated even eliminated in order to determine reliable and precise positions. The ionosphere is a dispersive medium to radio signal, so the value of the group delay or phase advance of GNSS radio signal depends on the signal frequency. Ground-based GNSS stations have been used for ionosphere monitoring and modeling for a long time. In this paper we will introduce a novel approach suitable for single-receiver operation based on the precise point positioning (PPP) technique. One of the main characteristic is that only carrier-phase observations are used to avoid particular effects of pseudorange observations. The technique consists of introducing ionosphere ambiguity parameters obtained from PPP filter into the geometry-free combination of observations to estimate ionospheric delays. Observational data from stations that are capable of tracking the GPS/BDS/GALILEO from the International GNSS Service (IGS) Multi-GNSS Experiments (MGEX) network are processed. For the purpose of performance validation, ionospheric delays series derived from the novel approach are compared with the global ionospheric map (GIM) from Ionospheric Associate Analysis Centers (IAACs). The results are encouraging and offer potential solutions to the near real-time ionosphere monitoring.     

Influence of land surface parameters on the spatio-seasonal land surface temperature regime in rural West Bengal,India

The magnitude and causes of changes in the land surface temperature of rural areas have not been extensively studied. The thermal band of Landsat imagery is taken to extract winter, summer, and monsoon season land surface temperature (LST) and relate it to surface parameters over a 30-year period. From the extracted parameters constructed a prospective surface temperature (PST) model using Multivariate Adaptive Regression Splines. The Chandrabhaga river basin in West Bengal of the lateritic Rarh Tract at the Chota Nagpur Plateau fringe was chosen as the study area because it is far from urban influences, to avoid the well-known heat island effect. Over the study period, summer and winter average LST increased linearly by 0.085?°C/y and 0.016?°C/y respectively. These results were validated with air temperature (RMSE?=?x and y, respectively). Over time more of the area is in the higher temperature zones, e.g., in April 2011, 4% area exceeded >32°, whereas in 2015 this proportion reached 52%. PST models of all the seasons were moderate to highly correlate (0.57–0.87) with actual LST, showing the value of this model. It also revealed the relative importance of the regional factors. Based on this information factor management is a scientific step to restrict or minimize the temperature rise effect.     

Solar and magnetic control on night-time enhancement in TEC near the crest of the Equatorial Ionization Anomaly

The ionospheric Total Electron Content (TECs), derived by dual frequency signals from the Global Positioning System (GPS) recorded near the Indian equatorial anomaly region, Bhopal (23.2°N, 77.4°E, Geomagnetic 14.2°N) were analyzed for the period of January, 2005 to February, 2008. The work deals with monthly, diurnal, solar and magnetic activity variations on night-time enhancement in TEC. From a total of 157 night-time enhancements, 75 occur during pre-midnight and 82 post-midnight hours. The occurrence of night-time enhancement in TEC is utmost during summer months, followed by equinox and winter months. The occurrence of night-time enhancement in TEC decreases with increase in solar and magnetic activities. We observed that peak size and half amplitude duration are positively correlated, while time of occurrence of night-time enhancement in TEC and time of peak enhancement are negatively correlated with solar activity. The peak size, half amplitude duration, time of peak enhancement and time of occurrence of night-time enhancement in TEC shows negative correlation with magnetic activity. The results have been compared with the earlier ones and discussed in terms of possible source mechanism responsible for the enhancement at anomaly crest region.     

改进的SSDA并行算法及其在TEM中的应用

介绍了一种把图像匹配中的序贯相似性检测算法(SSDA,Sequential Similarity Detection Algorithm)应用到地形匹配中的方法,把地形匹配系统中的数字高程模型、高程数据和实时剖面数据分别视为SSDA算法中的搜索图、灰度值和模板,同时提出了动态门限序列和分组选取随机点的方法,最后将改进的SSDA算法进行并行化并用消息传递接口(MPI, Message Passing Interface)实现.实验结果表明,改进的SSDA算法能有效地提高地形匹配的速度和计算精度,而相应的并行程序也能取得较好的加速比,从而一定程度上解决了传统的地形匹配算法时间复杂度高、实时性差的问题.