膜结构找形方法——改进力密度法

在膜结构找形中,力密度法引入了"力密度"的概念,使问题简化.但力密度法需要多次试算后才能确定合适的力密度取值,很大程度上依赖个人经验,给找形过程带来一定困难.本文针对力密度法的找形过程进行改进,避开其中"力密度"的概念,直接引入膜面应力和索拉力作为初始条件,以节点不平衡力作为控制误差,避免了传统力密度法需要反复试算力密度取值的弊端,使找形计算过程简洁高效.据此编制找形程序,通过复杂算例验证了该算法的正确性和普适性.     

张拉膜结构力密度法混合找形分析

论述了张拉膜结构力密度法混合找形的基本理论;所谓力密度法混合找形,即部分单元力密度控制,部分单元弹性控制;力密度控制采用线性求解,弹性控制采用非线性求解,通过迭代计算混合找形求出各结点的坐标值。据此编制了相应的计算软件;对工程实例进行了验算,结果表明,本文给出的计算结果与德国著名软件easy的计算结果相吻合。     

考虑面力密度的力法找形和预应力精确确定

随着膜结构应用领域的拓宽,实践对膜结构形状的构造提出了各种不同的要求。考虑膜单元的面力密度的力法找形方法,通过调整面力密度能够获得不同形状的曲面。对索膜结构找力分析的优化复位平衡法,能够对给定的拓扑曲面求解获得精确的自平衡预应力分布。把面力密度的力法找形和优化的复位平衡法的找力分析结合起来,能够完成形状构造的全过程,丰富了索膜结构的形状构造理论。     

力密度找形分析方法及计算机实现

本文通过对力密度法的分析研究，概括了力密度法的基本原理和找形分析的过程，用矩阵形式给出了力密度的基本列式。通过对各个矩阵基本元素分析，找到了力密度矩阵的组装规律和一维变带宽存储的规则及实现算法和过程，并给出了详细的力密度矩阵的装配过程。通过基于C 语言的系统框架的程序，实现了该基于一维变带宽存储的力密度找形分析方法。最后，通过几个数值算例表明了，力密度法找形理论的先进性，及本文提出的算法过程实现的可行性和正确性。由于本文的算法是基于效率比较高的一维变带宽压缩存储格式，因此，可以实现大型索膜结构的找形分析。     

张拉整体三棱柱几何作图法找形与找力

力学和几何学是密不可分的,本文给出了一种几何作图法进行端面平行时张拉整体三棱柱的找形方法和作图确定自平衡内力大小的找力方法,并推导出找形后张拉整体三棱柱自平衡内力的力密度计算公式,通过算例验证了几何作图法找形与找力的正确性。几何作图法找形和找力方法操作简单,直观可控,可以在CAD软件中实现,甚至可以纸上手绘实现;自平衡力密度公式简单,表达直观,是力密度法的一种几何实现。     

张拉膜结构的有限元分析

应用有限元程序ANSYS进行了薄膜结构的找形分析和荷载分析．张拉结构受力后将产生较大的位移,因此需要采用非线性理论．通过一系列计算,可以发现影响找形结果的因素有：单元的划分,初始弹性模量的选取,膜的预张力,边索和脊索的预张拉力值及其比值．在实际工程中,由于膜结构建筑要求和边界条件的限制,找到具有极小曲面性质的初始平衡形状较难,因此,一般把应力的差值控制在一定的范围内．工程实例计算结果表明：在风荷载作用下膜结构的竖向位移反应比较大,风吸力通常是索膜结构设计的控制因素．     

结构力学中拟无矩状态的计算问题

 通过两个算例介绍了结构力学中拟无矩状态的特点及其计算中应注意的问题,以及无矩理论 解在计算中的应用.     

钢筋混凝土圆球壳基于双剪屈服准则无矩理论的极限分析

利用已建立的基于双剪屈服准则的钢筋混凝土壳体的屈服条件,研究了圆球壳屈服时无矩理论下内力状态,导出了壳体各区域的本构关系和塑性屈服条件,并用机动法和静力平衡法计算了圆球壳在竖向均布载下的极限载荷.计算结果与用已有的塑性屈服线理论计算结果相符.     

索穹顶结构的找力分析方法--不平衡力迭代法

索穹顶的刚度和承载能力都是通过施加预应力来获取的,预应力设计是对索穹顶进行任何其它分析的基础,所以它的找力分析(Force finding)是一项很重要的工作。本文提出一种找力分析方法———不平衡力迭代法,它可以克服整体可行预应力法的一些不足:即使在索穹顶构件分组发生错误时仍可找到满足已知形状的预应力分布,并且该方法具有良好的稳定性和运算效率,算例表明不平衡力迭代法非常适合大型复杂索穹顶结构的找力分析。     

Conforming radial point interpolation method for spatial shell structures on the stress-resultant shell theory

The implementation of the conforming radial point interpolation method (CRPIM) for spatial thick shell structures is presented in this paper. The formulation of the discrete system equations is derived from a stress-resultant geometrically exact theory of shear flexible shells based on the Cosserat surface. A discrete singularity-free mapping between the five degrees of freedom of the Cosserat surface and the normal formulation with six degrees of freedom is constructed by exploiting the geometry connection between the orthogonal group and the unit sphere. A radial basis function is used in both the construction of shape functions based on arbitrarily distributed nodes as well as in the surface approximation of general spatial shell geometries. The major advantage of the CRPIM is that the shape functions possess a delta function property and the interpolation function obtained passes through all the scattered points in the influence domain. Thus, essential boundary conditions can be easily imposed, as in finite element method. A range of shape parameters is studied to examine the performance of CRPIM for shells, and optimal values are proposed. The phenomena of shear locking and membrane locking are illustrated by presenting the membrane and shear energies as fractions of the total energy. Several benchmark problems for shells are analyzed to demonstrate the validity and efficiency of the present CRPIM. The convergence rate of the results using a Gaussian (EXP) radial basis is relatively high compared to those using a multi-quadric (MQ) radial basis for the shell problems.     

A method to solve the efficiency-accuracy trade-off of multi-harmonic balance calculation of structures with friction contacts

The steady-state nonlinear forced response of systems with frictional damping can be computed in the frequency domain through the Harmonic Balance Method (HBM). A critical point is the selection of the number of harmonic terms used to represent the solution. In linear systems, this number is easily determined by the harmonic content of the forcing function (e.g. mono-harmonic). However, if nonlinearities are present, higher order harmonics may need to be included to ensure a proper representation of friction forces and displacements, with a detrimental effect on the computational time.The paper presents a novel method to solve the efficiency-accuracy trade-off of harmonic selection for nonlinear systems. This method warns the user whenever the number of retained harmonic terms is inadequate. As a result, it enables the user to run the simulation with a low number of retained harmonics, e.g. designers are typically interested in the first harmonic of the solution. The calculation is repeated with a larger harmonic support only when strictly necessary to keep the error below a user-defined threshold.The method is first compared to existing adaptive HBM techniques, highlighting its novel contributions. It is then carefully validated against high-order multi-harmonic calculations and against direct time integration. Its performance in terms of accuracy vs. computational time is highlighted. The method is then implemented in a state-of-the-art numerical tool for the design of underplatform dampers for turbine blades. Finally, its outcome is compared with experimental results.     

基于IGA-SIMP法的连续体结构应力约束拓扑优化

建立了一种IGA-SIMP框架下的连续体结构应力约束拓扑优化方法。基于常用的SIMP模型,将非均匀有理B样条（NURBS）函数用于几何建模、结构分析和设计参数化,实现了结构分析和优化设计的集成统一。利用高阶连续的NURBS基函数,等几何分析（IGA）提高了结构应力及其灵敏度的计算精度,增加了拓扑优化结果的可信性。为处理大量局部应力约束,提出了基于稳定转换法修正的P-norm应力约束策略,以克服拓扑优化中的迭代振荡和收敛困难。通过几个典型平面应力问题的拓扑优化算例表明了本文方法的有效性和精确性。应力约束下的体积最小化设计以及体积和应力约束下的柔顺度最小化设计的算例表明,基于稳定转换法修正的约束策略可以抑制应力约束体积最小化设计中的迭代振荡现象,获得稳定收敛的优化解;比较而言,体积和应力约束下的柔顺度最小化设计的迭代过程更加稳健,适合采用精确修正的应力约束策略。     

Improving convergence of incremental harmonic balance method using homotopy analysis method

We have deduced incremental harmonic balance an iteration scheme in the （IHB） method using the harmonic balance plus the Newton-Raphson method. Since the convergence of the iteration is dependent upon the initial values in the iteration, the convergent region is greatly restricted for some cases. In this contribution, in order to enlarge the convergent region of the IHB method, we constructed the zeroth-order deformation equation using the homotopy analysis method, in which the IHB method is employed to solve the deformation equation with an embedding parameter as the active increment. Taking the Duffing and the van der Pol equations as examples, we obtained the highly accurate solutions. Importantly, the presented approach renders a convenient way to control and adjust the convergence.     

Piecewise constrained optimization harmonic balance method for predicting the limit cycle oscillations of an airfoil with various nonlinear structures

A method is proposed to calculate the periodic solutions of piecewise nonlinear systems. The method is based on analytical derivation of nonlinear multi-harmonic equations of motion. Since periodic variations of nonlinear forces are characterized by different states, the vibration cycle is broken into sequential transition intervals according to the instant sets of state transitions. Analytical formulations of the harmonic coefficients of the nonlinear forces and its derivatives with respect to the harmonic coefficients of displacements are developed. Sensitivities of the harmonic coefficients of periodic solutions are determined for constructing explicit expressions for vibration amplitude levels as a function of structural parameters. Numerical investigations of the limit cycle oscillations and its sensitivities of an airfoil with different piecewise nonlinearities have been performed. The results show that the developed method is capable of determining the periodic solutions and its sensitivities with respect to the structural parameters. In order to guarantee time continuity of the nonlinear force, for the hysteresis model it is not right to track the periodic solutions by using the preload or freeplay as the continuation parameters.     

Modified method of characteristics for solving population balance equations

This paper presents a new numerical method for solving the population balance equation using the modified method of characteristics. Aggregation and break‐up are neglected but the density function variations in the three‐dimensional space and its dependence on the external fields are accounted for. The method is an interpretation of the Lagrangian approach. Based on a pre‐specified grid, it follows the particles backward in time as opposed to forward in the case of traditional method of characteristics. Unlike the direct marching method, the inverse marching method uses a fixed grid thus, making it compatible with other numerical schemes (e.g. finite‐volume, finite elements) that may be used to solve other coupled equations such as the mass, momentum, and energy conservation equations. The numerical solutions are compared with the exact analytical solutions for simple one‐dimensional flow cases. Very good agreement between the numerical and the theoretical solutions has been obtained confirming the validity of the numerical procedure and the associated computer program. Copyright © 2003 John Wiley & Sons, Ltd.     

统一骨架与连续体的结构拓扑优化的ICM理论与方法

技术了基于ＩＣＭ方法的结构拓扑优化新模型并应用于骨架与连续体结构。ＩＣＭ方法意指独立、连续变量与映射及其反演。新模型将两种结构统一建立了具有重量目标函数和多工况下应力与位移约束下的优化问题,提出的过滤函数是ＩＣＭ方法的关键技术之一。说明了优化策略与算法。     

基于自然单元法的轴对称结构极限上限分析

自然单元法是一种以自然邻近插值为试函数的新兴无网格数值方法,其形函数的计算不涉及矩阵求逆,也不需要任何人为参数。为了充分发挥自然单元法的优势,本文基于极限分析上限定理建立了轴对称结构极限上限分析的整套求解算法。轴对称结构的位移场由自然邻近插值构造,并且采用罚函数法处理材料的不可压条件。为了消除目标函数非光滑所引起的数值困难,采用逐步识别刚性区和塑性区,并对两者用不同方法进行处理。数值算例结果表明,本文提出的轴对称结构极限上限分析方法是行之有效的。