旋转拟稳定状态下的耦合热应力分析

本文引用旋转动坐标系，以导出在旋转拟稳定状态下的耦合热弹性方程，并应用以贝塞尔函数为坐标函数的加权余量法求解该方程。对数值结果进行了分析，导出了一些具有实际意义的结论     

一种用于计算含裂纹旋转壳动态特性的模型

本文用“含裂纹单元”模型（在垂直于裂纹面方向的单元尺度为零），求解含裂纹旋转壳的振动问题。同时给出了相应结构的实验数据，数值计算的结果和实验分析的结果具有很好的一致性。     

旋转三角形绳系卫星编队系统动态稳定性分析

本文研究了平动点附近的旋转三角形绳系卫星编队系统的动态稳定性问题.为了便于分析,假设系绳始终保持张紧状态,且不计系绳的质量.为了建立更加符合实际的数学模型,考虑了姿态运动和轨道运动的相互影响,即姿轨耦合,并在此基础上推导出了考虑姿轨耦合的非线性动力学方程.最后用数值方法模拟了在不同的旋转速度情况下系统的动态稳定性,模拟结果表明旋转速率越大系统稳定性越好.     

刚粘塑性强化球形薄壳在撞击体作用下的大变形动力分析

本文研究了在平头圆柱体的冲击下，薄球形壳的动力大变形。通过等度量变换，给出了球壳的变形模态。在Ｐｅｒｚｙｎａ和Ｓｙｍｏｎｄｓ的粘塑性本构模型基础上，设材料从Ｍｉｓｅｓ屈服条件，并假定弯矩的膜力不耦合。根据刚粘塑性强化材料的变分原理。给出球壳的运动方程结果与实验进行了比较，二者吻合较好。     

复合圆柱壳冲击压缩数值模拟及稳定性研究

针对复合圆柱壳在炸药爆轰作用下的动力学响应及在此过程中伴随的失稳问题,研究了其制造工艺中可能出现的缺陷以及圆柱壳中铜线螺旋角和直径对复合圆柱壳稳定性产生的影响。采用SPH-FEM耦合算法,建立了复合圆柱壳二维细节模型,并提出了一种基于圆柱壳内壁粒子速度历史的失稳判据,计算了在不同参数条件下复合圆柱壳的失稳时间及对应的压缩率,对影响复合圆柱壳稳定性的因素进行了评估。分析结果表明,在复合圆柱壳制备过程中存在的折返层缺陷和铜线直径对复合圆柱壳的稳定性有较大影响,而螺旋角度对其稳定性影响不大。     

基于约束关系的下颈椎运动测量新方法

下颈椎各个关节准确稳定的运动参数是医生对颈椎病进行诊断的重要依据。本文提出了一种满足颈椎耦合约束关系,基于双目视觉原理的下颈椎运动参数测量方法。首先,根据下颈椎围绕瞬时旋转中心运动的特点和下颈椎耦合运动的规律,建立了满足耦合约束关系的下颈椎运动学模型;其次,在下颈椎解剖学特征位置(棘突和横突)上布置标志点,借助多台摄像机,基于双目视觉原理和三维数字图像相关(3D-DIC)技术,直接测量不考虑耦合约束关系的下颈椎前屈、后伸、左旋、右旋及侧屈等运动时的角度;最后,利用考虑耦合约束关系的运动学模型对下颈椎运动角度变化曲线进行了校正。波动范围的分析结果表明测量的稳定性得到了提高。利用本方法实现了对下颈椎关节运动角度的实时测量,提供了一种对下颈椎关节运动进行稳定测量的有效手段。     

旋转双曲面自然通风冷却塔塔筒筒壁在集中载荷作用下的应力分析

对于双曲线冷却塔这样一种大型薄壳结构,通常的设计和计算都是在假定壳体承受分布荷载的前提下进行的,而在实际当中有时会遇到塔体受集中荷载作用的情况.本文针对这一问题,着重叙述了在计算中如何对集中荷载进行数学上的简化处理,利用现有的旋转壳应力分析程序给出一实际的数值算例,以及对计算结果的分析.     

高温外压下薄壁筒蠕变屈曲实验研究

本文介绍了作者研制的圆柱壳外压蠕变屈曲实验装置,用1Cr18N_i9Ti薄壁式筒在650℃下进行了外压失稳及蠕变屈曲实验,所得结果与理论分析有较好的一致性。     

在均布外压作用下复合材料三角形网格加筋圆锥壳体总体稳定性分析

本文在文献[1]的基础上加以推广,分析了在均布外压作用下复合材料三角形网格加筋圆锥壳体总体稳定性,考虑了蒙皮的正交拉弯耦合效应和加筋筋条的偏心效应.基于外压实验观察结果,通过选择恰当的位移分量表达式,本文利用Galerkin法获得了临界外压的近似表达式,并对某一C/E复合材料三角形内网格加筋圆锥壳体的临界外压值进行了理论计算,所得结果与实验值相当吻合.文中所得结果可供工程应用参考.     

无伞末敏弹气动特性风洞实验研究

为实现末敏弹的无伞稳态扫描运动,设计了一种短圆柱轴向非对称尾翼子弹气动外形。设计了模型和实验装置,进行了大迎角低速风洞实验。获得了模型在有/无减旋翼、固定和自由旋转条件下的气动力数据,测量了模型在气动力作用下的转速。实验结果表明:设计的尾翼能够为模型提供较大的阻力,旋转时模型的阻力系数有所减小;在有减旋翼的情况下,模型能够保持较低的稳定转速;模型的静稳定性较差;模型旋转时稳定性有较大的改善。需要进一步改进尾翼,提高静稳定性,为扫描角的稳定提供保证。实验结果可为无伞末敏弹的设计和改进提供参考。     

TSP问题的单元划分法

提出一种利用单元划分法求解ＴＳＰ问题的新求法，该法计算量小，且计算结果理想     

Analogy solutions of axisymmetrical twist problems by axisymmetrical finite element program

It used to be considered that an axisymmetrical problem and a twist problem of an axisymmetrical body cannot be simulated by each other, because the number of unknown variables in an axisymmetrical problem is greater than that in a twist problem, and the governing equations are not the same. This paper proposes a degenerated analogy method, by which the twist problems of axisymmetrical bodies can be simulated by axisymmetrical problems with finite element programs.An ordinary structural analysis method can be used to analyze an axisymmetrical problem, but a twist problem of axisymmetrical bodies is treated as a 3-dimensional problem usually. According to the method proposed in this paper, the analysis of a twist problem can be simulated by the analysis of an axisymmetrical body with a structural analysis problem. The example of analysis computation is also given. Thecomputed result is in agreement with the theoretical result.In this paper, the constitutive relation of the degenerated analogy problem is given.The authors suggest that a twist problem of a body made of any materials is simulated by an axisymmetrical problem of a body made of orthotropic material. If you have to use some program for the axisymmetrical problem to be limited to isotropic materials the penalty coefficient method can be used to solve the problem.     

流体饱和多孔隙介质波动方程多尺度反演

基于多尺度的思想,将小波多分辨分析和多尺度方法相结合,应用于流体饱和多孔隙介质孔隙率的反演。利用小波变换,将原始反问题分解为不同尺度上的一系列子反问题,并按照尺度从粗到细的顺序依次求解。在每一个尺度上,都采用稳定、收敛快的正则化高斯牛顿法求解,次一级尺度上求出的“全局最优解”作为上一级的初始解,依此类推,直到求出原始问题的真正的全局最优解。通过与传统的正则化高斯牛顿法相比较,显示了小波多尺度法是一个大范围收敛、能够有效节省计算量的方法,数值模拟的结果也表明了方法的有效性。     

一种改进格林元方法及在渗流问题中的应用

采用格林公式和基本解推导出直接边界积分方程来求解渗流问题.边界积分方程数值离散基于格林元方法(Green element methond),改进了原方法中压力和压力导数的求解方法,命名为混合边界元方法(Mixed boundary element method).相较于格林元类方法,该方法显式考虑了求解节点的外法向流量值和压力值,并使求得的数值解在求解区域上能够连续,符合实际的物理过程,在不增加额外未知数的情况下提高了计算精度.分析了不同网格类型对模拟计算结果的影响,并对稳定渗流问题、非稳定(瞬态)渗流问题和非稳态问题进行了实例计算,结果显示改进方法提高了计算精度,并对各类渗流问题有较好的适应性.     

Meshless methods for the inverse problem related to the determination of elastoplastic properties from the torsional experiment

The problem of determining the elastoplastic properties of a prismatic bar from the given experimental relation between the torsional moment M and the angle of twist per unit length of the rod’s length θ is investigated as an inverse problem. The proposed method to solve the inverse problem is based on the solution of some sequences of the direct problem by applying the Levenberg-Marquardt iteration method. In the direct problem, these properties are known, and the torsional moment is calculated as a function of the angle of twist from the solution of a non-linear boundary value problem. This non-linear problem results from the Saint-Venant displacement assumption, the Ramberg–Osgood constitutive equation, and the deformation theory of plasticity for the stress–strain relation. To solve the direct problem in each iteration step, the Kansa method is used for the circular cross section of the rod, or the method of fundamental solutions (MFS) and the method of particular solutions (MPS) are used for the prismatic cross section of the rod. The non-linear torsion problem in the plastic region is solved using the Picard iteration.     

Numerical method for optimum motion of undulatory swimming plate in fluid flow

A numerical method for the optimum motion of an undulatory swimming plate is presented. The optimum problem is stated as minimizing the power input under the condition of fixed thrust. The problem is singular for the invisible modes, and therefore the commonly used Lagrange multiplier method cannot predict an optimum solution but just a saddle point. To eliminate the singularity, an additional amplitude inequality constraint is added to the problem. A numerical optimization code with a sequential quadratic programming method is used to solve the problem. The method is applied to several cases of the motion of two-dimensional and three-dimensional undulatory plates, and the optimum results are obtained.     

变系数渗流场的正交各向异性边界元分析

本文给出了计算复杂渗流场的简化方法.首先由加权余量法得到变系数渗流问题的积分方程.再将区域及边界剖分,对各单元使用形心处的渗透系数而化为正交各向异性问题.最后用正交各向异性问题的基本解,得到对应于边界单元法的线性代数方程组.算例表明方法可靠.     

轴对称体扭转问题的轴对称有限元模拟解

轴对称体的轴对称问题与扭转问题一向被认为是两个互相不能模拟的问题.前者的未知量与方程多于后者,形式也不相同.本文提出一种退化模拟方法.能够把扭转问题模拟为轴对称问题的一类特殊情况来解.一般的结构分析程序都能够分析轴对称问题,但轴对称体的扭转问题通常作为三维问题处理.按本文提出的方法,可用结构分析程序的轴对称分析功能模拟扭转分析.本文还给出模拟计算的算例.计算结果表明与理论解完全一致.本文对退化模拟的材料本构关系进行了研究,建议在数值计算时以各向异性材料的轴对称问题模拟任何材料的扭转问题.当限定用各向同性材料的轴对称问题来模拟时,采用了罚系数法.     

Difference inversion model of wave equation

A numerical iterative model was derived from the difference method and a perturbation assumption to calculate the coefficient function of a wave equation. The method was used to solve the disaccord problem of numerical precision between the direct problem model and inverse problem model, and its serial problems using the old method. Numerical simulation calculation shows that the method is feasible and effective.     

Stress analysis near the tips of a transverse crack in an elastic semi-strip

The plane elastic problem for a semi-strip with a transverse crack is investigated. The initial problem is reduced to a one-dimensional continuous problem by use of an integral transformation method with a generalized scheme. The one-dimensional problem is first formulated as a vector boundary problem, and then reduced to a system of three singular integral equations(SIEs). The system is solved by use of an orthogonal polynomial method and a special generalized method. The contribution of this work is the consideration of kernel fixed singularities in solving the system. The crack length and its location relative to the semi-strip's lateral sides are investigated to simplify the problem's statement. This simplification reduces the initial problem to a system of two SIEs.