载荷位置不确定条件下结构稳健性拓扑优化设计

研究了在载荷作用位置不确定条件下,连续体结构的稳健性拓扑优化设计,为外载荷的有效传递提供可靠的传力路径。首先基于非概率方法,将载荷位置的不确定性用区间变量表示。其次通过载荷移动的一阶和二阶导数信息,推导出载荷位置改变后结构柔顺度对拓扑设计变量的一阶导数公式。最后运用变密度法中的SIMP模型以及移动渐近线法(MMA),对连续体结构进行稳健性拓扑优化设计。所得结构与传统载荷位置确定条件下的优化结构相比,稳健性明显更高,充分展示了考虑载荷位置变化对结构拓扑构型设计的影响及其重要性。     

预应力索-桁架结构形状优化设计

就预应力索-桁架结构形状优化设计问题,考虑了施加预应力阶段、预应力与荷载共同作用阶段的性态约束条件,建立了设计变量包括截面尺寸、索力值、杆件及索节点坐标的形状优化数学模型;在求解方法上将设计变量分为两个子空间:第一子空间为索力值和截面尺寸优化设计空间,第二子空间为形状优化设计空间; 第一子空间给出新的求解方法以减轻结构重量,第二子空间用节点渐进法优化结构形状和布索位置以增加结构刚度.算例表明,该方法能使结构重量逐步减轻,结构刚度不减,形状逐步达到最优.     

声场-结构耦合系统声压约束下板重量优化设计研究

为了减小振动板重量并使其振动噪声满足设计要求,基于声场-结构耦合有限元模型,提出了使结构重量最小化、满足声压约束的优化设计模型。用有限元直接法计算耦合系统的声响应及域点声压对结构设计变量的灵敏度,用可行方向法对其进行了优化设计研究。以一矩形平板和立方体声场耦合系统为实例,以壳厚度为设计变量,给出了域点声压对结构设计变量的灵敏度、优化前和优化后声压级对比图及最优的设计变量值。经过优化设计,在1Hz~200Hz频段内域点声压最大值降低6分贝,结构重量减少1.88千克。结果表明,声压约束下结构重量最小化设计,通过对结构重量的重新分布,能同时满足结构轻量化及噪声指标的要求。     

基于区间因子法的随机-区间模型的拓扑优化设计

利用3σ准则将随机变量近似的转换为区间变量,构建了区间栽荷作用下的区问参数平面连续体结构的拓扑优化设计数学模型.以结构的形状拓扑信息为设计变量,结构总质量均值极小化为目标函数,满足单元应力非概率可靠性和结构总体积为约束条件,研究随机-区间模型的拓扑优化设计问题.基于区间因子法,导出了单元应力响应的均值和离差的数学计算表达式,并采用双方向渐进结构优化方法进行求解.通过假设确定性模型最优结构拓扑设计方案中的各个参数与随机-区间模型中参数具有相同分散性,对不同模型进行优化.计算结果表明:在相同载荷作用下,确定性模型下的最优解很可能是随机-区间模型下的不可行解.     

离散变量结构拓扑优化设计研究

研究了离散变量结构拓扑优化设计的若干问题，讨论了离散型优化模型的合理性，提出截面设计变量的离散程度和全局约束影响最优拓扑，是优化中不可忽视的因素，文中还提出了一种解离散变量桁架，刚架结构拓扑优化的启发式算法。     

质量对结构优化理性准则法计算结果的影响研究

在结构优化中,质量作为设计资源能改善结构性能,但质量会引起自重与惯性载荷等问题,因而质量还会导致结构性能劣化。虚功准则方程组忽略了载荷对设计变量的导数,不能考虑质量导致结构性能劣化的一面。对于航空航天器、高速运行的车辆、机械和高精度天线等以自重和惯性载荷为主的工程结构,这种忽略导致虚功法得到的解与最优解相差甚远。导重法是严密推导的理性准则法,克服了虚功法准则不准的缺陷,优化效果较虚功法大幅提高。但对于考虑自重载荷的三杆、五杆桁架等简单结构优化问题,虚功法也可以得到与导重法相同的最优解计算结果,从拉格朗日条件出发解释该问题,给出考虑载荷对设计变量的导数时虚功法也可求得最优解的条件,这对于结构理性准则法及其工程应用的深入研究均具有重要理论与实践意义。     

结构构型与结构间连接方式协同优化设计

通常的结构拓扑优化是在给定外部作用下的构型设计. 然而对于复杂结构,各部件由连接构件相连并通过它传递外部作用,所以连接方式影响部件承受载荷的分布,即连接构件的布局对部件的最优构型有重要影响. 研究结构部件构型与部件之间连接方式的协同优化问题,在连接区域引入一种特殊的材料模型,以材料的分布描述连接方式. 将承力结构与连接构件域的材料密度同时作为设计变量,提出了一种基于拓扑优化思想的连接方式与结构拓扑优化协同设计的优化模型和相应的求解方法. 给出的算例证明了这种设计方法与优化策略的有效性.      

热结构瞬态响应的耦合灵敏度分析方法与优化设计

研究结构瞬态热变形和热应力的灵敏度分析方法及其优化设计,灵敏度计算给出了直接法和伴随法两种算法.考虑了温度场的耦合作用,在直接法中需要计算温度场对设计变量的导数,在伴随法中需要计算热载荷对温度场的导数.数值算例验证了该方法的精度.伴随法在应用程序中的实现,为大型结构优化提供了高效率的灵敏度计算方法.     

带频率禁区的结构单元尺寸二阶灵敏度动力优化设计

结构在服役期间经常受到风、地震或车载等激振作用，优化调整结构设计参数来设计结构动力特性以使其避开外部激励对应的卓越频率或高能频带，是提升结构安全性的有效手段。本文针对尺寸约束、频率禁区约束和结构重量最小化的结构动力学优化设计问题，基于Kuhn-Tucker条件和频率对变量的灵敏度，采用泰勒二级展开式推导并建立了结构频率对变量的二阶灵敏度单因子迭代计算方法，利用MATLAB平台编制了计算程序。算例表明比一阶灵敏度算法计算效率高、收敛稳定性强，且修正因子在优化全过程中无需调整，操作更为便利，并给出了单因子的合理取值区间。发现并初步论证了在设计变量未全部到达约束上下界时，“重量不再降低”只能作为最优的必要条件，“高阶频率收敛于频率禁区上限” 则应是充要条件，更适合作为收敛判据，其可有力甄别“伪最优”状况。基于数据表征，初步揭示了优化变量的修改主要受频率梯度主导，并发现频率梯度值大小与变量修正幅度大体成反比的变化规律。本研究对结构的抗风、抗震和在役结构的动力加固改造设计均具有重要的理论指导价值和现实意义。     

热结构稳态响应的耦合灵敏度分析方法

研究结构稳态热变形和热应力的灵敏度分析方法，给出了直接法和伴随法两种算法。考虑了温度场的耦合作用，在直接法中需要计算温度场对设计变量的导数，在伴随法中需要计算热载荷对温度场的导数。对尺寸和形状两类设计变量的灵敏度分析算例，验证了本文方法的精度。伴随法在应用程序中的实现，为大型结构优化提供了高效率的灵敏度计算方法。     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

汶川地震诱发文家沟巨型滑坡 碎屑流基本特征及成因机制初步分析


2008年5月12日,汶川M80地震在四川省绵竹市清平乡文家沟内诱发一巨型滑坡。通过现场调查得知,滑坡前后缘高差455m,厚度20～30m,滑面为基岩层面,初始方量2750×107m3。滑体在运动中转化为碎屑流。滑坡-碎屑流总的水平运动距离为4022m,垂直运动距离为1443m,遗留的堆积物体积达5×107m3。滑坡距映秀—北川断裂仅36km,位于其下盘,地震烈度达XI度。滑坡导致文家沟中48人遇害,并形成一条完整的地震次生地质灾害链。初步分析表明滑坡启动速度快,滑坡向碎屑流转化过程明显、地点明确。碎屑流运动过程复杂,伴有强烈的“气垫效应”和“前缘气浪冲击效应”。作者认为,文家沟滑坡的高启动速度是长持时强烈地震动作用的结果,与山体的猛烈碰撞是导致滑体解体并转化为碎屑流的原因。     

Preface: Orbits around asteroid

One of the core issues in modern celestial mechanics is the orbital dynamics in the near-regime gravitational field of as- teroids, which provides deep insights into the mathematical nature of a class of nonlinear systems, and plays as a critical basis for in situ explorations of different science goals. Lots of efforts have been made to reveal the characteristics of orbital motion in the vicinity of asteroids, and to improve the skills of asteroid research in methodology.     

Effet de l'anisotropie élastique cristalline sur la distribution des facteurs de Schmid à la surface des polycristauxInfluence of crystalline elasticity anisotropy on Schmid factor distribution at the free surface of polycrystals

Experimental studies of the plasticity mechanisms of polycrystals are usually based on the Schmid factor distribution supposing crystalline elasticity isotropy. A numerical evaluation of the effect of crystalline elasticity anisotropy on the apparent Schmid factor distribution at the free surface of polycrystals is presented. Cubic elasticity is considered. Order II stresses (averaged on all grains with the same crystallographic orientation) as well as variations between averages computed on grains with the same crystallographic orientation but with different neighbour grains are computed. The Finite Element Method is used. Commonly studied metals presenting an increasing anisotropy degree are considered (aluminium, nickel, austenite, copper). Concerning order II stresses in strongly anisotropic metals, the apparent Schmid factor distribution is drifted towards small Schmid factor values (the maximum Schmid factor is equal to 0.43 instead of 0.5) and the slip activation order between characteristic orientations of the crystallographic standard triangle is modified. The computed square deviations of the stresses averaged on grains with the same crystallographic orientation but with different neighbour grains are a bit higher than the second order ones (inter-orientation scatter). Our numerical evaluations agree quantitatively with several observations and measures of the literature concerning stress and strain distribution in copper and austenite polycrystals submitted to low amplitude loadings. Hopefully, the given apparent Schmid factor distributions could help to better understand the observations of the plasticity mechanisms taking place at the free surface of polycrystals. To cite this article: M. Sauzay, C. R. Mecanique 334 (2006).     

Nonlinear behavior of matrix-inclusion composites under high confining pressure: application to concrete and mortar

This paper is devoted to a micromechanics-based simulation of the response of concrete to hydrostatic and oedometric compressions. Concrete is described as a composite made up of a cement matrix in which rigid inclusions are embedded. The focus is put on the role of the interface between matrix and inclusion which represent the interfacial transition zone (ITZ). A plastic behavior is considered for both the matrix and the interfaces. The effective response of the composite is derived from the modified secant method adapted to the situation of imperfect interfaces. To cite this article: T.H. Le et al., C. R. Mecanique 336 (2008).     

Dynamics and Breakup of Viscoelastic Liquids (A Review)

The phenomena of hydrodynamic breakup of liquid jets, drops, films, bridges, and filaments are reviewed for liquids with viscoelastic properties. The reasons for breakup are capillary instabilities, collisions with rigid obstacles, and other forms of dynamic action. The relationship between the properties of the liquids and the features of the breakup process is discussed.     

Flow stability analysis and excitation using pulsating jets

Classical flow stability applied to transition from laminar to turbulent flow may also describe the behavior of vorticity fluctuations created by a pulsating jet placed along a solid boundary. A numerical laminar flow experiment involving a pulsating jet placed along the surface of a duct with flow separation downstream, resulted in eliminating most part of the separated flow region. Applying the same approach to a turbulent flow, it was possible to develop a turbulent stability flow formulation and apply successfully turbulent pulsating jet flow separation control. To cite this article: D. Skamnakis, K. Papailiou, C. R. Mecanique 333 (2005).