基于静力响应面的结构有限元模型修正方法

提出了基于静力响应面的结构有限元模型修正方法.运用响应面方法,将结构静力响应和结构参数之间复杂的隐式关系用显式函数近似表达出来;在此响应面模型(函数)基础上,通过优化计算对结构有限元模型参数进行修正.阐述了基于静力响应面的结构有限元模型修正方法的基本理论和一般实现过程.对两跨连续梁结构的静力模型修正数值算例分析结果表明:基于静力响应面的有限元模型修正方法可以减少结构有限元计算的次数、提高模型修正的优化效率,结构有限元模型修正结果具有可接受的精度.     

基于敏度分析的拉压不同模量桁架问题的数值分析

利用光滑函数技术,提出光滑化的拉压不同弹性模量问题的应力应变关系,与有限元方法相结合,建立了拉压不同模量一维连续体与桁架结构的数值求解模型,推导了敏度计算公式,采用Newton-Raphson算法进行求解.数值结果表明,本文算法具有较高的计算精度和收敛速度.     

水力劈裂问题的态型近场动力学建模

将态型近场动力学理论引入水力劈裂问题的模拟。构建了能反映岩土类材料准脆性断裂特征的态型近场动力学本构模型,并在物质点间相互作用力模型中加入等效水压力项,以实现在新生裂纹面上跟踪施加水压力。同时,考虑裂纹面间的接触,引入物质点间的短程排斥力作用,并设计了相应的接触算法。通过自编程序将模型和算法应用于含初始裂纹、不含初始裂纹以及含坝基软弱结构面的混凝土重力坝在高水头作用下的水力劈裂过程模拟,并与扩展有限元等模拟结果对比,验证了本文模型和算法的可行性和准确性。     

RTM充模过程数值模拟的隐式有限元算法

建立了基于欧拉方法描述树脂传递模塑(RTM)工艺充模过程的基本数学方程,并采用有限元隐式时间积分方法对基本方程进行了数值求解．编制了基于隐式有限元算法及传统有限元控制体算法的程序,通过具体算例比较了这两种算法的优缺点．与传统的有限元控制体法相比,该文提出的隐式有限元算法能节省计算时间,特别适合于单元、节点数目多的情况．隐式有限元算法是一种纯有限元方法,不需要使用控制体积技术,采用该算法计算出的流动前沿与时间步长无关。     

球形弹对金属靶板侵彻问题的数值模拟

基于球形空穴膨胀模型(SCE),采用ABAQUS有限元商业软件并结合自主开发的ABAQUS用户子程序对球形弹侵彻金属靶进行了有限元3D数值模拟。根据空穴膨胀理论,靶体对侵彻弹体的影响可以用一个作用在弹体表面的力函数代替,这样在进行数值模拟时就无需划分靶体网格,也避免了复杂的接触问题,从而使模拟大大简化。模拟过程中考虑到弹体的可变形性和入射时的微小偏航角,并且考虑了弹体在运动过程中和靶体的接触分离效应。模拟结果与文献中的实验结果进行了比较,模拟结果与实验结果吻合得很好。     

基于CUDA的有限元矩阵并行装配算法研究

构建航天飞行器的结构有限元模型是准确模拟飞行仿真、完成飞行器在轨飞行阶段结构故障监测和诊断的基础。采用细长体飞行器简化梁模型,提出新的基于CUDA(Compute Unified Device Architecture)的有限元单元刚度矩阵生成和总刚度矩阵组装算法。依据梁单元矩阵的对称性,结合GPU硬件架构提出并行生成算法并进行改进。为有效减少装配时间,在装配过程中采用着色算法,提出了基于GPU(Graphics Processing Unit)共享内存的非零项组装策略,通过在不同计算平台下算例对比,验证了新算法的快速性。数值算例表明,本文算法的求解效率较高,针对一定计算规模内的模型可满足快速计算与诊断的实时性要求。     

基于PAM-CRASH的鸟撞飞机风挡动响应分析

结合显式动力分析有限元软件PAM-CRASH及其提供的SPH算法,建立了鸟撞飞机风挡数值分析模型。对某飞机风挡进行了动响应分析,计算了风挡中轴线上四点位移时间曲线;鸟撞击过程的仿真结果表明,SPH鸟体模型能有效模拟撞击时鸟体溅射成碎片的情形;建立了鸟撞击作用下风挡破坏判据,对风挡在试验条件下是否破坏进行了模拟计算。计算结果和试验结果吻合较好,表明本文建立的鸟撞飞机风挡数值模型是有效的。     

针对有砟道床动力特性分析的嵌入式离散元-有限元耦合方法

在离散元-有限元耦合方法中,离散元和有限元交界面处的耦合方式对整体有砟道床的力学行为影响显著.采用基于球形单元的镶嵌单元或粘结单元模拟有砟道床时,由于球形单元和有限单元表面的自锁能力较差,使道砟层在列车载荷作用下容易产生侧向滑移,导致数值模型不稳定.此外,在实际铁路道床中,底部道砟均不同程度地嵌入路堤.为此,发展了一种嵌入式离散元-有限元耦合方法,通过设置一层嵌入地基有限元模型中的球形颗粒传递离散元域和有限元域间的力学参数,实现离散元和有限元方法的耦合.数值结果表明,嵌入式离散元-有限元耦合模型能够有效降低有砟道床的侧向位移,数值结果更加稳定,在处理与有砟道床类似的连续介质与散体介质的耦合问题时推荐采用嵌入式耦合算法.     

结构动力学有限元模型修正的目标函数及算法

结构动力学有限元模型修正是结构动力学领域的一个热点问题,回顾了结构动力学有限元模型修正研究的发展历史和现状,简要评述了结构动力学有限元模型修正所使用的设计变量,着重阐述了各种结构动力学有限元模型修正方法中所使用的目标函数及修正算法,讨论了工程结构动力学有限元模型修正的一些策略,最后对结构动力学有限元模型修正技术的发展进行了总结和展望.     

基于CUDA的有限元矩阵并行装配算法研究

构建航天飞行器的结构有限元模型是准确模拟飞行仿真、完成飞行器在轨飞行阶段结构故障监测和诊断的基础。采用细长体飞行器简化梁模型，提出新的基于CUDA（Compute Unified Device Architecture）的有限元单元刚度矩阵生成和总刚度矩阵组装算法。依据梁单元矩阵的对称性，结合GPU硬件架构提出并行生成算法并进行改进。为有效减少装配时间，在装配过程中采用着色算法，提出了基于GPU（Graphics Processing Unit）共享内存的非零项组装策略，通过在不同计算平台下算例对比，验证了新算法的快速性。数值算例表明，本文算法的求解效率较高，针对一定计算规模内的模型可满足快速计算与诊断的实时性要求。     

RESEARCH DEVELOPMENTS OF SMOOTHED PARTICLE HYDRODYNAMICS METHOD AND ITS COUPLING WITH FINITE ELEMENT METHOD

拉格朗日型的有限元法和光滑粒子法在模拟材料大变形问题时各存优缺点, 而有限元与光滑粒子耦合算法实现了在小变形区域采用有限元法计算, 在局部的大变形区域采用光滑粒子法计算, 有效地综合了有限元法计算效率高和光滑粒子法能够自然地模拟材料大变形问题的特点.重点论述了有限元法、光滑粒子法以及有限元与光滑粒子耦合算法的研究现状及应用进展, 并讨论了各方法中需要进一步解决的问题.最后通过算例对3种方法的计算精度和计算效率进行了分析, 供研究人员参考.     

基于结构动力修改的桩基检测方法

采用结构动力修改技术研究桩身损伤检测，将其转化为在一定约束条件下拟合桩顶频响函数的优化问题。用遗传算法加以求解。文中给出的数值模拟试验和工程应用实例，显示该方法能在较大范围内可靠地反演出桩身横截面积、桩长、波速等参数，达到检测桩身缺陷的目的，实现测试和分析一体化，使桩基检测易于进行并且结果实用可靠。     

三维厚度机翼上亚临界定常流的有限元计算

本文首次使用了曲六面体带节点导数等参元,通过变分法,计算了三维机翼亚临界定常位势流中的压力分布,计算结果与试验结果符合较好,文中还例用了有限元法中的局部线性化理论,插值出元素中心点密度去处理迭代过程中元素系数阵的计算,其结果与通常的有限元计算相吻合,但计算时间却明显减少了。文章最后以有限元法的局部线性化理论和局部线性化的有限元法做了阐述。     

Phenomenological and numerical modelling of short fibre reinforced cementitious composites

In this paper, a constitutive model for short fibre reinforced cementitious composites will be presented. This model is based on the St. Venant–Kirchhoff model, which is a special case of a hyperelastic material. This model is refined to include the fibre orientation distribution. Numerical FEM simulations with the developed constitutive model and fracture simulations using the discrete element method are presented. The outcomes of these numerical methods demonstrate how important it is to monitor and further to control the fibre orientation distribution during the manufacturing process. As the manufacturing process might involve casting, as, e.g., in the case of steel fibre reinforced concrete, an outlook on simulations of the manufacturing process in order to predict and to control the fibre orientation distribution is given.     

二维高亚音速Laval喷管流场的有限元计算

本文使用了八节点曲四边形等参元，通过变分有限元法，对二维高亚音速Ｌａｖａｌ喷管位势流场进行了计算，结果是满意的，并和一维结果进行了对比。文中还使用了有限元法特有的局部线化理论，取单元中心点密度去处理迭代过程中单元系数阵的计算，其结果与通常的有限元法计算相一致，但计算时间却大大减少了。     

Photoelastic analysis of the singular stress field in a bimaterial wedge

This paper presents an experimental investigation of the singular stress field near the vertex of a bimaterial wedge using a digital photoelastic technique. Special attention is given to the casting of bimaterial wedge specimens and analysis technique for extracting stress intensity factors from photoelastic samples. Different bimaterial wedge specimens are made of two different photoelastic materials bonded through a special casting procedure and loaded in simple tension. A new multiple-parameter method is developed to obtain the stress intensity factor reliably from the isochromatic fringe patterns and the series representation of the stress field at the vertex of the wedge. Experimental results are compared with finite element predictions, and good agreement is observed.     

有限元网格快速消稳的画家算法

消隐功能对于有限元分析模型的显示是相当重要的。本文针对有限元离散网格的特点,研究并实现了有限元网格隐藏面快速消隐的画家算法,这一算法先将所有单元面分成前向面和后向面集合,通过对各面顶点集合排序,可以极迅速地找到整个模型的边界前后面,对边界前向面采用一种新的动态排序画家算法实现了迅速高效的自动消隐。     

三维连续与非连续变形分析

将石根华博士所提出的二维非连续变形分析——Discontinuous Deformation Analysis（DDA）方法扩展到三维情况,并对三维不连续块体进行有限元网格剖分,即块体之间的接触采用DDA描述,块体内部的位移场和应力场则采用有限单元法描述,从而将三维DDA与有限元方法结合起来,增强了DDA方法与有限元方法解决实际工程问题的能力,实现了三维连续与非连续变形分析.给出了基本公式的推导过程和各子矩阵的形式.典型接触、碰撞算例证明了所提出方法的有效性和正确性.     

Neumann随机有限元的一种推广形式

针对Monte—Carlo随机有限元,研究在多样本环境下如何提高单样本计算效率的问题。证明了关于线性问题Neumann随机有限元法的算法与线性系统的等刚度迭代法是一致的;在此基础上,将迭代效率更高的算法一预处理共轭梯度法引入线性Monte—Carlo随机有限元系统,建立了相应的有限元列式;最后,利用两个算例比较了本方法与Neumann随机有限元法,结果显示随着随机变量离散度的增大,Neumann法的计算效率急剧下降,而本方法的计算效率表现稳定;对于离散度不太大的随机变量,Neumann法已经不收敛,而本方法在随机变量离散度很大的情况下,仍然保持很好的收敛性。