基于非线性变形场的柔性多体系统动力学

采用非线性变形场描述，将变形场表达为变形广义坐标的二阶小量形式，应用Ｋａｎｅ方法 建立了计及动力刚化项的柔性多体系统动力学模型，通过算例的数值仿真证明了本文模型的有效性和广泛的适用性。     

轴对称零件拉深过程中的最小无皱压边力

对轴对称零件的突缘变形区，动用能量法和二维折线式变位函数，计及材料常数ｎ、ｒ值和板厚变化，揭示了拉深中最小无皱压力的变化规律，指出其最大值明显滞后行突缘变形最大径向拉力的发生时刻，突破了两者基本上同时发生的传统结论计算结果与实验数据符合良好。     

空间可展开天线铰链热变形分析与试验验证

在轨热变形是影响空间可展开天线精度的关键因素之一,因此,空间可展开天线结构的在轨热变形分析具有十分重要的意义。随着空间探测任务要求的不断提高,铰链热变形对空间可展开天线指向精度的影响已经不可忽略。以某空间可展开天线的根部铰链和臂间铰链为研究对象,建立其有限元模型,计算铰链主弯方向和侧弯方向在温度每变化10℃时相对角度的变化。基于高低温箱和电子经纬仪构建铰链精度测试系统,开展铰链热变形试验。试验结果表明：根部铰链和臂间铰链在主弯方向和侧弯方向呈现纯结构热变形特性,仿真结果与试验分析结果基本一致,验证有限元模型和试验分析方法的正确性。根据铰链主弯和侧弯方向角度相对变化,评估热变形对铰链变形的影响。该建模和试验分析方法可为同类型空间可展开天线设计、分析和优化及类似可展开机构精度影响因素的分析提供参考。     

非对称偏馈环形网状天线索网预张力设计方法

提出一种考虑桁架柔性变形对索网平衡态影响的非对称偏馈环形网状天线索网预张力设计方法。该方法将天线整体结构分为内部索网与外部桁架两部分,通过内部索网预张力与外部桁架变形间的不断迭代,寻找满足索网形状及张力要求的天线平衡态,以实现网状天线索网预张力优化配置的目的。相比于常规基于刚性桁架假设的预张力设计,该方法可充分计及环形桁架与柔性索网协调变形,获得高精度反射面,且该方法适用于非对称偏馈天线,可有效避免后期张力调试,节约人力物力。算例分析验证了在进行预张力设计时计及桁架柔性的必要性及方法的正确性和有效性。     

大变形柔性铰链的静刚度分析及应用

大变形柔性铰链可由弹性较大的橡胶材料制成,而橡胶材料的应变-应力关系不符合胡克定律,其关系可由应变-能量函数关系式导出,为此,采用Mooney-Rivlin模型对形变能量函数进行简化,通过适当选取材料常数的比值,对其非线性的应力-应变关系进行分段线性化,导出圆角型柔性铰链在大变形时转角刚度的近似计算公式.通过有限元仿真对近似计算公式进行了验证.在此基础上,应用大变形柔性铰链构造扑翼式微型飞行器中的柔性胸腔,由转动电机作为驱动器,利用柔性铰链的变形实现翅膀的扇翅运动,实际应用效果表明可有效增大扇翅角并提高能量利用率.     

离心力场中的三维裂纹的J积分有限元计算

本文研究离心力场中的三维裂纹问题。在W.S.Blackburn给出的不计及离心力时的三维J积分公式的基础上,作者导出了离心力场中的三维裂纹的J积分公式为: J=integral from Γ(Wdy-T_iu_(i,1)ds)-integral from A((σ_(i3)u_(i,1),_3dA)-integral from A(F_iu_(i,1)dA) 本文将上式化为有限元表达式并研究了转换时的简化技巧;对中心含半圆表面裂纹的标准块在单向拉伸载荷作用下的J积分进行了计算和激光光弹实验;对某发动机涡轮盘的角隅裂纹进行了计及离心力的J积分计算。实践证明本文列出的方法是简便和切实可行的。     

变形反射面天线的效率

基于场相关原理,导出了适用于反射面缓变变形的反射面天线效率表达式,并讨论了它与Rutz 公式的关系,为反射面天线的设计提供了更为准确的电磁与结构关系式。     

基于考虑时间效应UH模型的一维固结分析

结合考虑时间效应的统一硬化(UH)模型得到一维固结下土的一维应力-应变关系,以此建立一维固结控制微分方程,可以统一考虑土的主、次固结变形.采用差分方法对单面排水固结问题进行了分析,加载初期不透水层处存在有效应力减小的松弛现象,推导了不透水层处有效应力随时间减小的公式并进行分析.随后,采用数值方法分析了次固结系数、超固结度、压缩回弹系数等因素对固结的影响规律,分析结果也验证了所得公式的结论.此外,通过室内试验对控制方程预测规律进行了验证分析.当主固结基本完成后,变形主要由次固结产生,基于这个条件推导了平均固结度与时间的关系,该公式可以用于计算固结基本完成需要的时间.      

工业机器人误差分析

本文讨论了工业机器人的静态误差和动态误差问题。导出了计算由杆件公差引起的手部静态误差的公式。用动态静力学方法分析了由杆件变形引起的手部动态误差,并以例子说明了计算过程。     

计及壁面传热的微型向心涡轮气动设计探讨

毫米级微型向心涡轮传热效应显著区别于常规向心涡轮,发展计及壁面传热效应的微型向心涡轮的气动设计方法对微型燃气轮机的设计具有重要意义.根据微型向心涡轮壁面非绝热的特点,通过理论分析,初步建立了计及固壁传热效应的微型向心涡轮速度三角形的分析模型.在此基础上发展了计及传热效应的微型向心涡轮的气动设计方法.分析得到了考虑传热影响的微型向心涡轮设计规律.利用三维数值模拟手段对模型和设计方法进行了验证.结果初步表明所发展计及壁面传热效应的模型和微型向心涡轮气动设计方法具有较好的可行性.     

蜂窝夹层结构剥离应力分析

为了研究蜂窝夹层结构受剥离力时面板的变形和蜂窝芯中的应力,对蜂窝夹层梁进行了数学分析和实验验证.将建筑中的弹性地基上梁的理论用于蜂窝夹层结构分析,视面板为弹性地基上的梁,视蜂窝芯为弹性地基,得出了蜂窝夹层结构受剥离力时的面板变形和夹芯中的应力分布的计算公式.通过分析和实验指出,将建筑中的弹性地基上梁的理论用于蜂窝夹层结构分析,具有计算量小、精确度较高和公式直观等优点.      

血管的新应变能函数(Ⅰ)

首先提出并且验证了可被应用于大变形计算的描述血管非线性应力-应变关系的新的本构理论.其推导的数学基础是假设应力-应变关系可以用张量多项式来表示.通过对比可以发现,拟合结果与试验数据完全吻合,并且优于已有的指数形式的表达式.为降低Taylor展开式的次数,对级数进行了修正,并且针对不同的拟合过程和方法进行了讨论.      

整体壁板压弯成形等效计算模型

分析了筋条内弯型整体壁板压弯成形过程中蒙皮和筋条的弹塑性变形分布,鉴于内弯筋条发生塑性变形而蒙皮主要为弹性变形的特点,提出了一种将蒙皮部分简化为虚拟材料筋条的几何等效方法.推导了基于惯性矩的虚拟材料等效系数计算公式,给出了复杂网格筋条壁板等效模型建立过程和虚拟材料系数计算方法.进行了Ⅰ型单筋条和网格筋条壁板的有限元分析和压弯实验,对比分析了完整模型和等效模型的计算效率和精度,并用实验数据进行了验证,表明建立的压弯成形有限元等效模型精度可靠、效率数倍提高,为大型整体壁板压弯成形分析和工艺参数优化提供了基础.     

遥感器球面反射镜蜂窝孔轻量化优化设计

针对遥感器球面反射镜背部蜂窝打孔减重结构,编制单元划分、节点生成程序,使之自动建立有限元模型文件;应用结构分析软件对该结构进行自重载况下的静力计算,变动各有关结构参数,得出一系列镜面最大变形值;处理这些数据,拟合曲线,进而导出该结构轻量化最优解的经验公式。     

血管的新应变能函数(Ⅱ)

用多项式表示高度非线性的本构关系时,多项式中往往含有很多高阶项,而要获得与这些高阶项对应的材料常数, 存在许多困难.为降低多项式的次数,应在多项式中加入修正项.对Taylor展开式进行修正的方法,以及在曲线拟合过程中需要使用的技巧有许多值得探讨之处.      

蜂窝结构异面变形动态塑性坍塌力学

在蜂窝结构静态塑性坍塌应力理论计算模型的基础上,建立了蜂窝结构在冲击载荷作用下的动态塑性坍塌应力的理论计算模型.考虑应变率对蜂窝结构动态力学性能的影响,使用Cowper-Symonds模型建立了静态塑性坍塌应力与动态塑性坍塌应力之间的关系;基于薄板弯曲理论,由平均冲击速度得到了应变率的计算公式,进而得出动态塑性坍塌应力关于冲击速度与静态塑性坍塌应力的函数关系.使用LS-DYNA动力学软件模拟了铝合金蜂窝结构在冲击载荷作用下的异面变形,数值计算结果与理论模型计算结果一致.     

ε-N曲线三参数幂函数公式

给出了一种描述ε-N曲线的三参数幂函数公式,与著名的Manson-Coffin公式相比,它能更好地拟合各种材料的应变-寿命试验数据,并且可以更方便地用于计算结构的损伤和寿命。     

Modeling and analysis of Earth’s gravity field measurement performance by inner-formation flying system

The Earth’s gravity field can be measured with high precision by constructing the purely gravitational orbit of the inner-satellite in Inner-formation Flying System (IFS), which is independently proposed by Chinese scholars and offers a new way to carry out gravity field measurement by satellite without accelerometers. In IFS, for the purpose of quickly evaluating the highest degree of recovered gravity field model and geoid error as well as analyzing the influence of system parameters on gravity field measurement, an analytical formula was established by spectral analysis method. The formula can reflect the analytical relationship between gravity field measurement performance and system parameters such as orbit altitude, the inner-satellite orbit determination error, the inner-satellite residual disturbances, data sampling interval and total measurement time. This analytical formula was then corrected by four factors introduced from numerical simulation of IFS gravity field measurement. By comparing computation results from corrected analytical formula and the actual gravity field measurement performance by CHAMP, the correctness and rationality of this analytical formula were verified. Based on this analytical formula, the influences of system parameters on IFS gravity field measurement were analyzed. It is known that gravity field measurement performance is a monotone decreasing function of orbit altitude, the inner-satellite orbit determination error, the inner-satellite residual disturbances, data sampling interval and the reciprocal of total measurement time. There is a match relationship between the inner-satellite orbit determination error and residual disturbances, in other words, the change rate of gravity field measurement performance with one of them is seriously restricted by their relative size. The analytical formula can be used to quantitatively evaluate gravity field measurement performance fast and design IFS parameters optimally. It is noted that the analytical formula and corresponding conclusions are applied to any gravity satellite which measures gravity field by satellite perturbation orbit.