固体推进剂三维粘弹性Monte-Carlo随机有限元法

发展了一种粘弹性Monte-Carlo随机有限元法,并对固体火箭发动机药柱进行了三维随机分析。首先由Herrmann泛函出发导出了三维不可压和近似不可压粘弹性增量有限元列式。然后,同时考虑固体推进剂力学性能参数和载荷的随机性,结合Monte-Carlo技术对弹性约束的圆柱形中孔药柱进行了随机分析。为了克服直接Monte-Carlo法(DMCS)效率低的缺点,引入Latin超立方抽样技术(LHS)以提高计算效率。算例表明该方法收敛速度块、通用性强、易于工程应用。     

二维几何非线性问题的线有限元解法

提出了在线有限元方法中考虑几何非线性二维问题的理论公式和计算步骤，并以完全的拉格朗日表述为基础，首次导出了用线有限元方法进行结构分析中大变形问题载荷增量法的基本公式，对计算过程中的重要方面：包括对非线性微分方程线性化的方法，分步迭代的步骤等进行了详细讨论，以这些公式为基础，编制了相应的计算程序，并给出了算例，结果是令人满意的。     

微分求积单元法分析变曲率井内受径向约束管柱的非线性屈曲

基于Love的空间弯扭杆平衡方程,通过引入井壁约束条件导出了平面变曲率井内受径向约束管柱的平衡方程。采用微分求积单元和增量迭代法求解了曲率线性变化的井内管柱的非线性屈曲问题,通过与有限元计算结果的对比验证了所构建方法和编写程序的正确性,而且还表明微分求积单元法有方法简单、易于实施,计算量少、精度较高等优点。计算结果表明,等曲率井内管柱屈曲的临界载荷明显大于曲率线性变化井内管柱屈曲的临界载荷,另外,变曲率井眼的曲率变化对管柱弯矩、井壁约束力有显著的影响。     

旋翼／机身耦合系统的固有特性研究

直升机的旋翼与机身在桨毂处是机械耦合的，用旋翼阻抗和机身阻抗在桨毂处匹配的方法分析旋翼／机身耦合系统的固有特性。旋翼阻抗由旋翼在固定坐标系中的运动方程得到，这个运动方程由旋翼旋转坐标系中的桨叶模态方程导出，机身阻抗可由此可见 身的有限元分析计算或振动试验获得，还设计与加工了一个旋翼／梁耦合模型，对此模型进行了试验研究，试验结果和计算结果比较表明，分析旋翼／机身耦合系统固有特性的方法是可行的，且具有较高的精度。     

基于VUMAT的三维裂纹动态扩展有限元模拟

文章基于ABAQUS/EXPLICIT平台的VUMAT接口,直接对材料积分点添加自定义本构方程和失效准则,开发了三维裂纹动态扩展子程序.在网格适当加密的情况下就能完成三维裂纹动态扩展模拟.在加载策略上选用适应性较强的增量法,直接在软件中通过增量步控制求解过程和收敛性,计算结果与试验结果一致.本文方法基于商业有限元软件ABAQUS,仅对材料本构方程和积分点失效进行编码,无需调整和重新划分网格,操作简单,计算效率高,适合推广到复杂的本构关系和大型工程模型.     

粘弹性阻尼结构的动特性分析

本文研究了粘弹性阻尼结构动力特性的分析方法。根据分数导数模型表示的粘弹性材料本构方程,对一种阻尼材料用拟合法构造了它的力学性能解析表达式。建立了一种八自由度夹层梁有限元模型,从而给出了粘弹性阻尼结构的典型动力学方程。为分析粘弹性阻尼结构的动力特性,提出了一种非线性特征值问题的迭代解法。应用上述方法对悬臂约束阻尼层梁的复模态参数进行了分析,并进行了实验验证。结果表明,分析与实验的结果一致,说明本文所提出的方法是可靠的且有较好的精度和实用性。     

弹性容器内液体晃动的有限元法

本文给出了弹性容器内液体晃动的变分公式,通过单元离散和未知函数的插值导出了它的有限元方程,为解决方程的非对称性给出了一种方程对称化的方法,最后编制了计算程序,并与刚性壁容器内液体的晃动做了比较。     

高精度微分求积曲梁单元的建立与应用

首先由能量原理导出曲梁弯曲问题的控制微分方程，在此基础上应用微分求积法原理分别给出了曲梁内点和端点的微分求积方程，由此形成曲梁单元的刚度方程，从而建立了微分求积曲梁单元，并给出了曲梁结构刚度方程的边界条件。通过算例分析，得到了微分求积单元法结构离散时应使单元数量少的原则和求解精度与单元长度比基本无关的性质。与有限元方法的结果比较表明，本文导出的曲梁单元是一种具有很高求解精度的单元。》     

框、壳相互作用力q的弹塑性公式

本文根根载荷及结构的对称性,导出了柱壳平衡方程及挠度方程,并作了求解。文[1]弹性框的假定是不合理的,计算与试验结果相差很大。本文导出了框壳均进入塑性后q的计算公式,计算结果与试验数据相当符合,在工程上有直接使用价值。     

串级萃取最优萃取比精确方程

萃取比是优化串级萃取工艺设计中的重要参数。本文导出了计算最优萃取比的精确方程。相同条件下，最优萃取比精确议程的计算结果误差小于现有的最优萃取比方程。当分离系数小于3时，计算最优萃到比既可用最优萃取比方程，也可用最优萃取比精确方程。当分离系数大于3时，计算采用最优萃取比精确方程更为合适。     

固体推进剂药柱的近似不可压缩粘弹性随机有限元

采用Monte-Carlo粘弹性随机有限元法,分析了药柱的结构可靠度。针对固体火箭发动机药柱材料近似不可压缩的问题,基于Herrmann粘弹性有限元以提高计算的精度,考虑泊松比和松弛模量的随机性,采用拉丁超立方(LHS)技术计算了固体火箭发动机药柱结构的可靠度。算例表明该方法的有效性。     

最小二乘和奇异值分解法在边界元法中的应用

采用最小二乘与奇异值分解结合的方法，给出求解系数矩阵不满秩的线性代数方程组的数值方法，进而将此方法应用于边界无法中，处理给定外力的第一边值问题，特别地用于处理给定外力的三维裂纹问题。此外，本文还给出求解三维有限体裂纹问题的超奇异积分方程组，并使用有限部积分与边界元法为其建立了数值法。最后计算了若于典型例子的应力强度因子，数值结果与现有文献的相比，符合很好。     

稳定悬浮状态下无轴承开关磁阻电机电感特性

针对二维有限元和磁链法计算无轴承开关磁阻电机(Bearingless switched reluctance motor,BSRM)电感特性的局限性,建立了BSRM稳定悬浮时的三维有限元模型,采用双标量磁位法计算得到其磁场分布,应用一种仅需计算非线性系统能量增量即可计算电感的增强型能量增量法,计算其绕组电感.以-7.5 kW,12/8结构的BSRM为例,简单比较了双标量磁位法与二维有限元法、矢量磁位法和棱边有限元法分别求解BSRM磁场时的计算规模和平均计算时间,说明了双标量磁位法在大规模计算场合的优越性;研究了电感与转子径向位移的关系;计算获得了样机处于稳定悬浮状态下的主、悬浮绕组电感;研究了端部效应对电机磁场分布和绕组电感的影响.通过实验验证了电感计算的准确性,为优化设计和精确控制该类电机奠定了基础.     

有限元频响分析中网格划分技术研究

阐述了利用有限元进行频响分析过程中存在的一个突出问题；简单介绍了波动理论的一些基本公式，并提出有限元频响分析时划分网格的参考标准；首先通过对简单结构(包括梁、矩形板)的有限元计算以及与解析解的对比初步验证了该方法的有效性，然后给出了两个工程算例的应用情况。     

一种计算复合材料层合板面外应力的新型线性分析方法

复合材料层合板壳结构的面外应力对于层间分层的起始和扩展具有重要意义。首先简要地给出了经典的非协调位移有限元列式,然后基于修正的H-R变分原理导出了单元面外应力变量与位移变量的关系。面外应力的有限元线性系统为面外应力边界条件的引入提供了方便。同时,通过该线性系统可直接获得节点上唯一且连续的面外应力结果。与非协调位移有限元的结果比较,本文的方法可明显地提高数值结果的精度,且计算效率及资源要求均优于非协调辛元法和非协调广义混合元法。     

反复荷载下钢筋混凝土框架节点非线性分析

论述了反复荷载下分析、计算普通钢筋混凝土及轻骨料钢筋混凝土框架节点的非线性有限元分析方法 ,采用等效单轴应变的增量正交异性混凝土本构关系模型和反复荷载下应力应变关系的“焦点模型” ,编制二维非线性有限元分析程序 ,基于四个混凝土边节点的试验结果 ,对钢筋混凝土和轻骨料钢筋混凝土框架节点进行非线性有限元对比分析 ,计算的结果与试验吻合较好。     

动态有限元模型建立的逆方法

本文介绍一种用测量的特征对(归一化振型和固有频率)建立结构系统动态有限元模型的方法。有限元模型用质量和刚度参数表示,这些参数由正交条件和特征方程的加权欧氏范数取极小值来确定。对不完全测量特征对条件下解的唯一性和不完全测点集的处理方法进行了讨论。文中给出的实例表明,这里提出的方法是可行的。     

Fast SSED-MoM/FEM Analysis for Electromagnetic Scattering of Large-Scale Periodic Dielectric Structures

A hybrid method combining simplified sub entire domain basis function method of moment with finite el ement method （SSED-MoM/FEM） is accelerated for electromagnetic （EM） scattering analysis of large-scale peri- odic structures. The unknowns are reduced sharply with non-uniform mesh in FEM. The computational complexi ty of the hybrid method is dramatically declined by applying conjugate gradient-fast Fourier transform （CG-FFT） to the integral equations of both electric field and magnetic field. The efficiency is further improved by using OpenMP technique. Numerical results demonstrate that the SSEI）-MoM/FEM method can be accelerated for more than three thousand times with large scale periodic structures.