复合材料弹塑性多尺度分析模型与算法

对材料非线性多尺度分析的计算模型与算法进行研究.在构建周期分布单胞分析算法的基础上,发展针对复合材料结构材料非线性多尺度分析的一般有限元方法.方法的特点是将所建立的单胞分析过程作为有限元分析的子程序嵌入到总程序系统当中,完成对应的高斯点应力计算,因而使所发展的方法具有实现方便的特点.给出数值计算结果,验证了方法与所发展的多尺度有限元分析程序的正确与有效性.     

饱和土壤固结分析的算法研究

对饱和土壤固结分析的广义Biot方程求解,提出一个新的算法.文中给出了算法的精度及稳定性条件.     

从规划法求解看有摩擦接触解的不唯一性

以单点接触问题为例研究了有摩擦接触问题的几个特点，给出了这类问题出现不唯一解的条件。进而从二次规划求解角度分析了解的不稳定性与不唯一性的数学原理，文中给出了规划法求解有摩擦接触问题应注意的问题     

基于氮氧自由基的新型稀土配合物的合成、结构及光谱性质

合成了一种新型的氮氧自由基-稀土配合物[Dy（hfac）3（H2O）2]·（NITMePy）(1, CCDC:846691, hfac=六氟乙酰丙酮,NITMePy=2-（2-吡啶基-3-甲基）-4,4,5,5-四甲基氧化咪唑啉-3氧化-1-氧基自由基),其结构和光谱性质经元素分析、红外光谱和紫外光谱表征。结果表明：配...     

解Drucker-Prager塑性问题的二阶锥互补法

基于经典弹塑性理论中多数屈服准则具有凸锥数学结构的事实,将在大规模计算中更具潜力的锥规划法引入弹塑性分析。考虑到弹塑性流动理论有关联与非关联之分,本文提出利用锥型互补法求解弹塑性问题。具体以Drucker-Prager弹塑性模型为例,首先利用最大塑性功耗散原理和圆锥对偶理论等工具,建立了弹塑性本构方程的等价二阶锥互补模型;然后,基于参变量变分原理和有限元技术,建立了弹塑性增量分析的二阶锥线性互补模型;最后,利用一类半光滑Newton算法求解。数值算例表明了本文方法的有效性。     

颗粒材料内部滑动变形局部化的离散元模拟

颗粒材料在加载过程中表现出非常复杂的变形性质.在微观尺度上,单个颗粒的运动、单个孔隙胞元内的变形都是非均匀的,但也不是完全随机、没有规律可循的,而是呈现一定结构性的分布,如微带、剪切带等.本文用基于孔隙胞元的离散元方法对颗粒体进行双轴加载数值试验,模拟了以滑动变形表征的变形局部化现象.数值结果发现,对应加载过程中的不同...     

冠脉支架纵向柔顺性数值模拟

柔顺性是冠脉支架重要的力学性能之一。本文针对一种具有非对称结构的冠脉支架,利用有限元方法深入研究与弯曲方向有关的支架柔顺性,在柔顺性数值模拟中使该支架获得较大的弯曲曲率并且考虑了接触因素。在一个周期结构内的多个弯曲方向上分别进行模拟,详细研究了支架抗弯刚度的各向异性。计算结果显示：在支架自接触现象发生前,支架的柔顺性表...     

基于SiPESC平台的模态综合法研究及实现

模态综合法作为主流的动力子结构方法,可在有效的精度下缩减整体结构的自由度,特别是对大型复杂结构的动力计算具有较高的效率,已被集成到许多商业CAE软件中,因此把它集成到SiPESC平台上是一项很有意义的工作.本文实现了两种模态综合法的集成,分别是Craig-Bampton方法、多层Craig-Bampton方法.研究工作同时提出了一种高效的子结构界面点搜寻算法,所实现的Craig-Bampton方法考虑了子结构的旋转及平移,多层Craig-Bampton方法可以考虑各种多层子结构的划分方式.SiPESC平台是由大连理工大学开发的新一代的有限元数值计算软件平台,它利用了新型的软件架构和程序设计思想,因此在SiPESC平台上进行程序设计开发是十分方便的.通过构建数值算例同Nastran的计算结果进行对比,数值算例验证了本文实现方法的正确性和有效性.     

基于哈密顿动力系统新变分原理的保辛算法之二：算法保辛性质证明

文献[1]给出了哈密顿系统的一个新的变分原理,并基于此变分原理,通过选择一个时间步长两端不同广义位移或广义动量为独立变量,给出了四种不同类型的求解哈密顿动力系统的数值方法。本文将分别证明这四类数值方法都是保辛的数值方法。     

Improved precise integration method for differential Riccati equation

An improved precise integration method(IPIM) for solving the differential Riccati equation(DRE) is presented.The solution to the DRE is connected with the exponential of a Hamiltonian matrix,and the precise integration method(PIM) for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix.The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE.Based on the error analysis,the criterion for choosing two parameters of the PIM is given.Three kinds of IPIMs for solving the DRE are proposed.The numerical examples show that the IPIM is stable and gives the machine accuracy solutions.