多个裁剪参数曲面的有限元混合网格剖分法

根据冲压模具形状特点,对由CAD/CAM系统建立的冲压模具的多裁剪参数曲面数学模型,提出了一种符合冲压成形有限元分析需要的三角形和四边形混合网格剖分算法.按几何离散控制参数即最长边、最短边、距离误差和角度误差,将多裁剪参数曲面离散为多边形集.利用曲面曲线查找裁剪参数曲面相邻关系,建立B-Rep表示的曲面几何模型,并更新多边形集.然后由多边形集得到三角形和四边形混合网格,结果网格单元是相容的,在相邻裁剪参数曲面边界处无裂缝和覆盖.      

基于自动有限元建模的民机机翼结构布局优化

针对民用飞机中广泛采用的双梁式机翼结构,以机翼结构纵向构件(长桁)和横向构件(翼肋)的数量为参数生成机翼CAD模型,然后基于Patran的PCL(Patran Command Language)语言,根据整体模型分割、局部模型编号、局部网格控制、整体有限元生成的思路,实现了以纵向和横向构件的数量变化为基础自动进行机翼结构有限元模型的构建及分析,并将上述流程集成到机翼结构的布局优化中.最后参照某大型民用飞机的机翼外形尺寸建立了CAD模型,以机翼的静强度、刚度和蒙皮稳定性为约束,以机翼结构质量最低为目标,对上、下蒙皮长桁和翼肋数量进行了优化.优化后的布局使机翼结构质量降低了10.1%,表明了布局优化方法的有效性和可行性.     

板料成形数值模拟接触搜索模型的研究

基于参数曲面描述法和有限元网格描述法,提出模具几何形状模型的曲面单元描述法,采用了逐级更新的接触搜索算法,推导出求交迭代初值优化选择方案及满足无穿透原理的接触判断准则,并建立了板料成形数值模拟接触搜索模型,此模型不仅改善了数值模拟的收敛性,而且提高了计算效率和稳定性.     

任意曲面上射线的寻迹方法

为计算任意曲面上天线间的隔离度,提出一种射线寻迹的新算法——最小夹角法.首先对任意曲面表面进行三角形网格剖分,并按照文中要求的格式生成模型表面网格数据;用文中提出的“最小夹角法”对上述网格数据进行处理,从而达到快速准确地找出任意曲面上、任意两点间的短程线.求解了位于圆柱体、椭球体及任意组合体上任意两点间的短程线长度,验证了“最小夹角法”的准确性.     

任意曲面上射线的寻迹方法

为计算任意曲面上天线间的隔离度,提出一种射线寻迹的新算法——最小夹角法.首先对任意曲面表面进行三角形网格剖分,并按照文中要求的格式生成模型表面网格数据;用文中提出的"最小夹角法"对上述网格数据进行处理,从而达到快速准确地找出任意曲面上、任意两点间的短程线.求解了位于圆柱体、椭球体及任意组合体上任意两点间的短程线长度,验证了"最小夹角法"的准确性.     

飞机复合材料构件设计、分析和制造数据集成

实施复合材料构件数字化定义是实现复合材料构件设计、分析和制造数据集成的基础.复合材料构件数字化定义包括几何建模和材料结构建模.总结了材料结构建模过程和材料铺放数据的组成,指出当前复合材料构件分析无法基于最终的真实纤维路径进行,采用离散三角形单元的真实纤维路径来近似铺层纤维的最终状态,并利用可扩展标记语言(XML,Extensive Markup Language)文件对材料铺放数据进行了描述,探讨了通过纤维路径数据到有限元网格的自动映射实现材料铺放数据的自动传递.最后,研究了设计数据到排样下料系统和激光投影系统的集成和接口数据组织方法.      

飞行器RCS计算前置处理中裁剪曲面剖分算法

在雷达散射截面(RCS)计算中,首先需要将物体的表面进行三角剖分,称为前置处理,三角片单元的边长通常为0.1个波长(一般1个波长取3 cm).飞行器外表形状复杂,在建立数学化模型中,使用了大量的裁剪曲面,曲面之间的关系非常复杂,而3 mm边长的三角片剖分,使得飞行器的三角片单元的数目巨大,给前置处理带来难度.对此,提出了一种裁剪曲面三角网格剖分的方法.利用飞行器外形设计的特点,在以拱高为逼近误差的前提下,把曲面离散一系列等参数线,等参数线上等弦长布点,两两参数值相邻的参数线生成网格单元.根据网格单元与裁剪区域的位置关系,将落在裁剪区域外的网格单元进行裁剪处理.然后对曲面间边界处的网格进行相容性处理,最终实现多张裁剪曲面的三角网格剖分.      

半规则三角网格模型细分曲面重构

针对逆向工程中的三角网格重构问题,提出了一种保持尖锐特征的半规则三角网格模型细分曲面重构算法,以充分利用细分曲面的多分辨特性.首先经提取尖锐特征和删除最大独立点集得到基网格,然后利用插值Loop细分和最近点法向投影法不断调整半规则网格得到重采样网格,最后运用渐进插值(PI,Progressive Interpolation)算法生成半规则细分曲面.实例表明,重构后获得的细分曲面连接性好,可以直接进行细分小波分析.     

基于边界保持的非闭合细分曲面小波变换

针对非闭合Loop细分曲面,提出一种小波变换的方法.该方法利用网格边界所在的平面作为映射镜面,将非闭合的网格曲面转换为具有对称形状的闭合网格曲面.在此情况下,原网格边界上的顶点变为新生成闭合网格的内部顶点,在对其进行小波变换后将镜像复制的部分去除即可得到新网格及其边界.该算法同时能够保证网格曲面边界在进行小波变换过程中始终保持在同一平面内,为基于多分辨技术的细分曲面数控加工定位面的确定带来方便.     

基于造型系统的曲面及二维域有限元网格剖分

介绍了一个基于几何造型系统的有限元分析前置处理系统．其主要功能是对几何造型系统的ＮＵＲＢＳ曲面以及二维任意形体有限元网格划分．该系统可作为几何造型系统与有限元分析程序之间的接口．     

高超音速热流数值计算中的误差匹配原则

通过对粘性边界层的数值模拟和误差分析,提出了差分格式的误差匹配原则.在边界层内的计算网格一般都具有很大的长宽比,而计算误差与网格长宽比密切相关.为了使不同坐标方向上的误差保持在同一量级上,需根据误差匹配原则确定计算网格的尺寸.对绕钝锥的高超音速流动进行了数值模拟,进而计算了钝锥表面的热流分布.数值计算表明,参照误差匹配原则进行网格划分是有益的,可在保持流场数值计算精度和提高热流计算准确性的同时,使网格的数量和分布更为合理.     

空域空间网格标识原理及应用

在对空域空间有限元的相关研究中，构建空域空间的网格标识系统，是将信息技术应用于空域管理、推进空域数字化的重要基础。在美军全球区域参考系统和通用地理位置参考系统基础上，提出了一种同时具备标识平面和立体空间位置能力的参考系统。通过在地球表面以上特定高度之下的空间内，建立一种递归剖分策略方法，将空域空间划分为若干空间网格，不同尺度下的网格既是一种地理空间位置参考系统，也是一种组成空域的基本空域体，以此实现空间位置参考与空间位置标定的混合表达，为开展空域数字化管理提供底层模型。      

Loop细分曲面精加工刀具轨迹生成

细分曲面既能表示连续的几何设计模型也能表示离散的加工模型,避免了模型转换的复杂中间过程.细分曲面除了对于构造具有任意拓扑结构的复杂零件具有巨大的优势外,对于数字化制造也极具发展潜力.因此,对基于Loop细分曲面的精加工刀具轨迹生成算法进行了研究.首先利用基于弦长误差的自适应插值Loop细分得到精加工模型;然后利用等斜率跟踪法将精加工模型分割为平坦区域和非平坦区域.对各个区域依次进行处理,不同区域实施不同的刀具轨迹规划.据此,既避免了因细分过程数据量过大而导致的曲面精度不足或表面质量降低的问题,又可顺利计算整张曲面相对均匀残留高度的刀具轨迹.最后,为验证该算法的可行性进行铣削加工实验.      

Geomagnetic cutoffs: a review for space dosimetry applications.

The earth's magnetic field acts as a shield against charged particle radiation from interplanetary space, technically described as the geomagnetic cutoff. The cutoff rigidity problem (except for the dipole special case) has "no solution in closed form". The dipole case yields the Stormer equation which has been repeatedly applied to the earth in hopes of providing useful approximations of cutoff rigidities. Unfortunately the earth's magnetic field has significant deviations from dipole geometry, and the Stormer cutoffs are not adequate for most applications. By application of massive digital computer power it is possible to determine realistic geomagnetic cutoffs derived from high order simulation of the geomagnetic field. Using this technique, "world-grids" of directional cutoffs for the earth's surface and for a limited number of satellite altitudes have been derived. However, this approach is so expensive and time consuming it is impractical for most spacecraft orbits, and approximations must be used. The world grids of cutoff rigidities are extensively used as lookup tables, normalization points and interpolation aids to estimate the effective geomagnetic cutoff rigidity of a specific location in space. We review the various options for estimating the cutoff rigidity for earth-orbiting satellites.     

B样条应用于超限插值的网格生成方法

采用B样条曲线曲面生成技术与传统网格生成方法中超限插值技术,并将两者结合起,形成一种计算网格的生成方法.它发挥了两者的优点和长处:能充分反映复杂外形的完全信息,提高对外形的适应程度;可控制切矢的边界条件,提高生成网格的质量.算例表明了上述特点.      

地质构造的三维可视化

地质构造三维可视化可以将地质勘探数据用直观的图形方式显示,是正确认识地质构造的重要手段,从而为油藏描述提供科学的依据.通过根据已知地质数据的三角剖分和插值,建立了地层层面、断面三角网格模型.研究了地质层面、地质整体模型、地质立体剖切和动态模型等多种三维地质模型可视化流程.提出了建立地质立体模型的方法,即由缝合相邻地层层面、相邻断线分别形成的地质体边界面和断层三角形环状曲面,以及顶地层三角网格层面,共同围成三维地质立体模型.在这一过程中,实现了最近邻优先重构、地层层面三角网格模型的外边界闭合圈多边形追踪及平面与三角网格求交线等关键算法.实验结果提供了从整体到局部的多种地质构造显示手段,从而帮助人们准确快速掌握地质构造.      

基于目标特性的航天器太阳光压面积计算方法

太阳光压是影响深空探测航天器轨道确定与预报精度最主要的摄动力.针对实际任务需求，采用了一种基于目标特性的光压面积建模与计算方法，根据航天器形状、尺寸、表面材料以及材料光学特性等信息，实现了分析型光压模型的建立与求解，提高了计算效率和精度，可快速计算目标在光照方向上的光压面积、投影面积以及光压比例因子等参数.通过长方体光压面积理论值与仿真值的对比，验证了该方法的准确性和有效性.针对复杂结构探测器开展了光压面积计算，可为深空探测航天器精密定轨中的光压模型解算、定轨及预报提供参考.      

探空火箭定时定点入轨的一种制导控制规律

针对特定探测天体,给出了特殊用途的探空火箭与其实现空间交会的时刻与地点的计算方法.根据特定天体的运行轨道,发射前算出标称交会飞行轨道,装订在箭载计算机内.火箭发射后,利用箭载惯性导航系统确定自身当前的位置与速度,比对标称飞行轨道参数得出飞行偏差,通过控制火箭推力偏斜调整飞行轨道,使探空火箭在交会时刻到达交会点,并在交会时刻相对与惯性空间的速度为0.定义了研究所用的各种坐标系,建立了火箭飞行动力学方程.研究了标称飞行轨道最优交会点选取,交会时间与发射时间计算等问题.给出了发射后动力飞行段的制导控制规律,核心思想是将控制信号分解为时间控制、当地水平面上的海拔高度控制、南北控制与东西控制,通过设置偏置量减小关机后轨道摄动因素引起的漂移.利用计算机数值仿真验证了这种制导控制规律的可行性.     

基于陆态网的区域电离层TEC空间变化特性分析

为研究中国陆态网区域电离层TEC在空间小尺度、高分辨率情况下的变化特性及适用精度范围，利用陆态网260个GNSS连续运行观测站数据，解算并生成2016-2017年731天陆态网区域电离层RIM格网，并进行精度验证.在同一RIM格网中，分别在经度和纬度方向上对间隔不同经纬度的TEC格网点作差分析.结果表明:陆态网区域内经度方向上TEC最大变化率和平均变化率分别为0.30TECU·（°）-1和0.11TECU·（°）-1；经度间隔1°时，TEC差值小于2TECU，且随着经度间隔的增大，其TEC差值也随之增大，并表现出一定的半年和周年变化规律；纬度方向上TEC最大变化率和平均变化率分别为1.7TECU·（°）-1和0.46TECU·（°）-1；陆态网区域内电离层TEC随纬度减小而增大，纬度间隔1°时，99.4%的TEC差值小于4TECU，且随着纬度间隔的增大，其TEC差值也随之增大，并表现出一定的半年和周年变化规律；间隔相同情况下，纬度方向上TEC的变化比经度方向上大.      

超音速轴对称有旋流特征线法的计算程序

编写了一个ＦＯＲＴＲＡＮ７７程序，计算了Ｍａ０为２．０￣５．０的超速轴对称有旋流Ｃａｕｃｈｙ问题的一些特征线网，但推进方向是从初值线在物面，这样，可以对物形进行反设计，在设计科波机时，可以选用这种流场中的流面作为其升力面的组成部分。