正交各向异性旋转壳的极限分析

提出了一种轴对称壳极限分析的数值方法。该方法首先用拟协调元构造弹性应力场,然后用温度参数法构造残余应力场,将极限分析转化为一个数学规划问题,最后用序列二次规划(SQP)算法进行求解,得到了基于H ill屈服准则的正交各向异性旋转壳的下限值。同时由于von M ises屈服准则是H ill屈服准则的一个特例,因此该方法也可以用于各向同性轴对称壳的极限分析。分别计算了服从H ill屈服准则和von M ises屈服准则的受内压的组合壳的下限,计算结果可为各向异性压力容器的设计提供参考。     

基于Hoek-Brown准则的非常规态型近场动力学弹塑性模型

岩石介质在不同荷载环境下具有显著的弹塑性变形特征和断裂力学行为。为了同时描述岩石的弹塑性变形连续场以及断裂力学行为的非连续场,提出了一种基于Hoek-Brown强度准则的非常规态型近场动力学（non-ordinary state-based peridynamics, NOSBPD）弹塑性模型。首先,基于非局部变形框架下的材料对应性关系,建立了以Hoek-Brown强度准则为屈服准则的NOSBPD弹塑性本构方程,通过主应力空间的返回映射算法得到给定应变增量对应的应力增量,并给出了相应的增量模型积分算法;其次,提出了基于等效塑性应变的断裂准则,实现了岩石弹塑性断裂全过程力学行为的表征;最后,通过引入非均匀变形状态消除了NOSBPD模拟存在的零能模式问题。基于数值模拟结果与有限元结果以及试验数据进行对比分析,验证了所建立模型的正确性,可为岩石的弹塑性断裂力学行为研究提供有效分析手段。     

不同拉压模量厚壁球壳的弹塑性分析

基于理想弹塑性模型,对拉压模量不同的厚壁球壳进行弹塑性分析.在不同内压作用下,对厚壁球壳进行了弹性分析、弹塑性分析,最后进行塑性极限分析,并求解厚壁球壳各个状态下的应力及位移的大小.分别讨论了弹塑性边界位置不同时,球壳内壁开始屈服的弹性极限压力及球壳整个截面完全屈服的塑性极限压力.     

胞映射及其绘图方法

简要地说明了胞映射法分析了非线性系统全局性态的理论和方法；提出了一种绘图数据的压缩存储算法。这种算法可以节省大量存储空间和传送数据所需的时间，为胞映射法解决工程实际问题提供了条件，根据图形情况可节省存空间１０－１０００倍，文中附有两个数值实例。     

带接管高压球形容器塑性大应变的应力分析

本文在应力测试和爆破试验的基础上,采用轴对称问题塑性大应变有限元程序对带接管高压球形蓄力器进行弹塑性应力计算,分析带接管球壳承载后在连接部位的应力分布、最大应力点的位置及应力集中系数,按Mises屈服准则计算接管连接部位的局部屈服压力、连接部位以及球壳的全面屈服压力;进一步分析连接部位和球壳的塑性区域扩展过程,结构弹塑性和塑性阶段的应力分布和应力集中系数。本文又分析了容器的爆破试验和理论应力计算结果,探讨高压容器的破坏失效机理及爆破压力的计算方法。     

关于半隐式向后Euler本构积分算法的改进

在相关联流动法则与各向同性强化的条件下,对弹塑性有限元中半隐式向后Euler本构积分算法进行改进,改进后的算法,对较大的时间步长,调整后应力仍能返回到更新的屈服面上,从而加强了算法对时间步长的适应能力。对于特殊屈服准则,如D-P准则,本文算法本质上为半隐式径向返回。最后给出了算例说明。     

功能梯度圆柱壳的弹塑性屈曲

研究了功能梯度材料(FGM)圆柱壳在轴向均匀压缩载荷作用下的弹塑性屈曲行为.基于线性混合强化弹塑性模型,给出FGM圆柱壳的材料特性表达式和弹塑性本构方程.引入Hamilton原理,将FGM圆柱壳的弹塑性屈曲行为转化为求解辛空间的特征值问题.进一步利用分叉条件计算出正则方程广义特征值对应的屈曲临界载荷,联合屈服条件获得屈曲壳的弹塑性分界面位置.并讨论了材料梯度参数、结构几何参数对弹塑性屈曲临界载荷和弹塑性分界面的影响.     

叶片有限元分析中弹塑性过渡区应力奇异产生原因及解决方法

采用弹塑性有限元法、借助大型商业有限元软件对汽轮机叶片进行应力分析时,弹塑性过渡区应力的计算值有时会高于塑性区应力的计算值,即会产生应力奇异现象。为分析产生这一现象的原因,以8节点六面体单元为例,研究了有限元法计算应力的过程,并在理想弹塑性的条件下,采用有限元法和解析法计算了弹塑性过渡区单元节点应力。研究发现,有限元法通常采用高斯积分点应力值外推插值法得到单元节点应力,当单元一部分位于弹性区、另一部分位于塑性区时,这种外插算法会导致节点应力计算值高于结构的实际应力,甚至超出理想弹塑性材料的屈服极限,从而造成应力奇异。研究表明,在叶片弹塑性的有限元分析中,采用相邻高斯积分点应力加权平均的方法计算单元节点应力,可有效避免弹塑性过渡区应力产生奇异的现象。     

管节点弹塑性大挠度有限元分析

文中根据Von Mises屈服准则及Pradtl-Reuss增量关系,推导出了理想弹塑性材料、四节点板壳单元的增量本构关系,根据上述理论和利用有限元程序ANSYS实现了对管节点的弹塑性大挠度分析,并对三个管相贯节点作了弹塑性大挠度的全过程分析。     

稀疏矩阵文档图像处理

讨论基于稀疏矩阵的文档图像存储及处理方法 .采用三向量法或链表法表示稀疏图像 ,然后在稀疏域直接实现某些基于临域运算的图像处理算法 .分析表明 ,对于具有显著特征的文档图像能有效地节省存储空间并提高计算效率 .以卷积运算和一种文档图像处理运算为例 ,给出实验结果     

用弹塑性有限元素法解平压头下的压入问题

本文采用弹塑性有限元素法中的变刚度法和加权平均弹塑性矩阵,解光滑压头下平面应变的弹塑性硬化材料的压入问题。得到压力——位移曲线,屈服压力和塑性区的边界,并与理想刚塑性的滑移线场理论结果相比较。同时计算了在平面应变条件下,工业纯铝压缩时的冷变形抗力。计算结果与试验结果基本相符。     

二维弹塑性自然单元法算法实现

为了使自然单元法能够应用于土体等具有弹塑性本构关系的材料的分析计算,通过结合弹塑性理论及自然单元法自身特点,研究了在自然单元法中采用Von—Mises、Mohr—Coulomb和Drucker—Prager屈服准则解决二维弹塑性问题的算法,并利用面向对象的程序设计方法编制了相应的计算程序．通过算例验证了各屈服准则下算法的正确性,证明了自然单元法相对于常规有限元算法在精度上的优势．在自然单元法中实现了Mohr—Coulomb和Drucker—Prager屈服准则,拓展了自然单元法的适用范围．     

一种多波前并行处理的有限元算法

本文提出了一种求解大型有限元系统的新算法。该算法采用并行处理结构,首先将结构分成许多子区,然后利用多个波前在各个子区内并行地组集,消元,从而得到凝聚后的界面刚度阵和载荷阵。再串行组集和求解界面方程得界面位移,最后返回各个子区,并行求解内点应力和位移。从运算结果来看,该方法不但能有效地提高运算速度,减少计算时间,同时能有效地节省内存量,是一种求解大型结构有限元系统的有效途径。     

基于网格资源可用度的元任务映射算法

针对Min-Min启发式算法没有充分考虑网格的异构特性,导致资源负载不均衡的问题,提出一种基于资源可用度的启发式映射算法,在元任务与资源映射过程中,先计算资源对元任务的可用度,然后将元任务划分为优先映射与次优先映射两部分,再对这两组元任务分别应用经典Min-Min算法进行映射.采用仿真软件模拟网格异构环境,在元任务与资源一定的情况下,分别模拟这两种算法的元任务与资源的映射过程,得到这两种情况下的任务完成时间.实验结果表明,基于资源可用度的启发式算法在异构的网格环境中任务完成时间更短.     

边界元弹塑性分析的预测-校正迭代法

本文提出一种预测-校正迭代方法,使边界元初应变法的弹塑性分析可以适用于包括非强化材料在内的各种塑性流变特性的材料。其原理是利用结构因子即内部单元对其初应变的弹性响应来校正每步增量载荷作用下预测的塑性应变增量,使塑性区的等效应力逐步返回到屈服面上。在中型计算机上对均匀拉伸的双V缺口平板按理想弹塑性材料进行了弹塑性应力分析,计算结果与有限元法的结果十分一致,但计算时间大为缩短。在均匀内部应变的基础上,还给出一个区域型积分和边界型积分的转换关系,当采用简单几何形状的常数内部单元时,利用这个关系可以方便地得到内部单元区域积分的精确表达式。     

空间多样化约束下的移动k近邻查询

考虑为移动中的查询对象连续返回k个距离近并且满足空间多样化约束的对象，提出了空间多样化约束下的移动k近邻(SDC-MkNN)查询.在此，满足空间多样化约束代表对象之间的相互距离大于距离阈值.为了高效处理SDC-MkNN查询问题，提出了两种基于安全区域技术的算法.算法均通过减少重新计算查询结果的次数来提高查询效率.其中一种为精确算法EA，可连续返回精确的查询结果；另一种为近似算法ρAA，可连续返回具有近似率保障的近似查询结果.采用真实数据集验证了所提出算法的有效性.     

高频RCS预估中判别阴影区域的并行算法

利用物理光学方法预估高频雷达散射截面(RCS)时,需要判别目标表面的阴影区域。本文提出了一种基于树型结构的多重阴影判别算法,将判断面元是否被面元遮挡转变为判断面元是否被分组遮挡,极大地降低了计算复杂度。与此同时,引入并行技术,将基于树型结构的多重阴影遮挡算法并行化,进一步节省计算时间。数值算例表明,基于树型结构的多重阴影判别算法与传统的阴影判别方法相比,在保证计算精度一致的前提下,计算效率更加高效。     

一种提高铸件凝固过程数值计算效率的方法

基于显式差分法,提出了一种占用内存少交能提高砂型铸件凝固过程数值计算效率的方法,采用一个数据块存储铸件／砂型系统的温度场以节省内存,对砂和铸件区域采用不同的时间步攻分别进行计算以提高计算速度,并选择适当的方法处理个个计算域的耦合问题,通过一个实例证明,该算法可使模拟计算速度提高１倍,而计算所得的铸件温度场几乎与传统算法一致。     

一种适用于无线网络的椭圆曲线标量乘算法研究

单标量乘和多点标量乘是椭圆曲线密码体制中最主要的运算,在各种密码协议中起到核心作用。提出一种基于Shamir方法优化预计算的多标量乘算法,算法对标量从左到右进行编码,并将编码阶段与多标量乘的主计算阶段合并进行,节省了存储标量新编码的内存空间,更适合于内存受限的无线环境中应用。分析表明新算法在一般情况下效率可提高6％～50％。     

基于无向图转有向图的同构判别

无向同构图指的是在两个图中寻找顶点之间的映射关系,通过映射使原本形式各异的两图中的各条边保持对应的关系.为了有效提高寻找无向同构图的时间效率、简化操作,首先研究了无向图同构的矩阵存储方式,并针对性地提出了把无向图转换为有向图的同构算法.与矩阵存储算法相比,该判定算法的时间更为简短.最后给出了实现该算法的相关程序以及用该算法对无向图进行判定的过程和结果.