含压电片梁的静变形控制

首先借助阶梯函数，建立了含有任意分布的用作执行器的压电片的梁挠曲轴线方程。然后利用该方程，进行了梁的静变形控制的研究。最后结合实例，用“遗传＋配点”法对压电片的位置和尺寸进行了优化。结果表明，本文建立的挠曲方程形式简洁，求解简单：“遗传＋配点”法优化效果明显。     

多工况下结构鲁棒性拓扑优化设计

针对工程中存在多个随机不确定工况载荷作用的情况, 将鲁棒性设计思想引入到连续体结构拓扑优化设计, 发展和完善不确定性优化理论和计算方法. 基于概率模型和SIMP方法,提出以结构柔顺度标准差最小化为目标、具有体积约束的连续体鲁棒性拓扑优化数学模型.通过对目标函数及其灵敏度计算公式的推导, 采用数学规划法实现优化问题的求解. 数值算例验证了所提优化模型的正确性及算法的有效性, 并通过与确定性优化结果的比较,证明鲁棒性拓扑优化能够给出结构柔顺性变异更小的材料分布.      

多相材料结构拓扑优化: 体积约束还是质量约束?

针对拓扑优化设计中多年沿用的体积约束作为控制材料用量的思路, 研究了采用质量约束开展多相材料结构拓扑优化的新方法, 提出了相应的质量约束结构优化问题模型和材料插值模型. 通过研究广义SIMP模型和对等混合模型两类多相材料属性插值方法, 发现广义SIMP插值模型导致质量约束具有复杂的非线性特征, 给优化问题的全局寻优和数值求解带来困难;采用的线性对等混合材料插值策略保证质量约束函数具有线性和变量可分离特征, 为拓扑设计问题的数值求解, 特别是凸规划寻优提供了极大便利. 算例表明多相材料插值策略和结构优化模型的正确性和有效性; 优化设计结果的比较说明了使用质量约束相对于体积约束的合理性和必要性, 阐述了质量约束在解决结构轻量化布局设计问题的有效性以及所带来的显著收益.      

基于拓扑优化的声学结构材料分布设计

本文针对结构的声学设计问题进行研究,通过优化两种不同的材料在结构设计域内的拓扑分 布来最小化谐振结构所产生的声场中指定参考面/参考域内的声压。在研究中假定结构为线弹性小变 形结构,材料阻尼为Rayleigh阻尼,声学介质为无粘、可压缩、小扰动流体。对结构响应采用有限 元格式进行计算,对声场采用基于Helmholtz积分的边界元格式进行计算,由于声场在无穷远自由边 界的无反射条件在边界积分中能自动得到满足,该格式特别适合于具有开放边界的声场计算。建立 了结构有限元-声场边界元格式的耦合系统拓扑优化模型,导出了耦合系统敏感度分析的一般格式及 伴随格式。数值算例验证了所提出的结构-声学耦合系统优化方法的有效性和可靠性,并揭示了基于 声学准则的拓扑优化结果的有关特性 关键词边界积分,结构声学耦合系统,拓扑优化,敏感度分析,伴随方法     

道路约束条件下的移动机器人跟踪

为提高移动机器人的位置估计精度和跟踪效果,提出一种基于道路约束条件下的移动机器人鲁棒约束H∞滤波(CHF)跟踪算法.首先,将机器人移动的道路网络作为跟踪的约束条件,并利用当前统计模型对机器人的运动进行建模.其次,将道路约束条件作为机器人跟踪的非线性状态约束,利用最小协方差估计推导了鲁棒CHF递推方程.通过拉格朗日乘子法对非线性约束优化估计问题进行求解,并利用约束信息对CHF算法的状态更新过程进行了改进.最后,通过对CHF算法和无约束的H∞滤波算法的跟踪性能进行了对比分析和验证.仿真结果表明,该算法可以实现机器人的跟踪,且跟踪精度优于HF算法.     

混凝土复合介质爆炸的不确定优化设计及其高效求解方法

研究不确定因素对炸药在混凝土复合介质中爆炸的影响,并在不确定因素的影响下进行优化设计。用区间数模型来描述系统中的不确定参数,建立相应的区间数优化模型;通过高维代理模型的引入,以及对其中二层嵌套优化问题进行直接解耦,得到了一种高效的求解方法。研究表明：算法经过8个迭代步收敛到最优解,最优解得到的最佳炸药装填深度约为64cm,最佳姿态角约为30°;经过三个迭代步就能迅速定位到最优解附近。算法保证了在最大化毁伤面积的同时,尽量地追求不确定因素影响最小,即综合考虑了毁伤最大和偏差最小。该不确定优化设计方法能有效地推广到相关领域的优化设计,为处理爆炸系统的装药设计和爆炸初始状态参数设计等方面提供了新思路。     

杆件分组和网格密度对网壳截面优化结果影响

针对跨度为50\,m, 矢跨比为1/5的单层凯威特K6型球面网壳结构, 对6种杆件分组方式、 7种网格划分密度的42个网壳模型进行截面优化分析, 得出了与以往截面优化不同的、优化 效果更好的杆件分组方式和网格划分密度取值建议, 优化过程充分考虑了网壳结构强 度、刚度和整体稳定性, 使优化结果更具有工程实用价值. 可为工程设计人员提 供有益的参考.     

Enhanced multi‐objective optimization of a dimpled channel through evolutionary algorithms and multiple surrogate methods

Using a multi‐objective evolutionary algorithm (MOEA) and enhanced surrogate approximations, the present study demonstrates the numerical analysis and optimization of staggered‐dimple channels. Two surrogates, the response surface approximation (RSA) model and the Kriging (KRG) model, are applied in light of the surrogate fidelity of the approximate analysis. An enhanced Pareto‐optimal front is obtained by performing local resampling of the Pareto‐optimal front, which provides relatively more accurate Pareto‐optimal solutions in the design space for each surrogate model. Three dimensionless design variables are selected, which are related to geometric parameters, namely, the channel height, dimple print diameter, dimple spacing, and dimple depth. Two objective functions are selected that are related to the heat transfer and pressure loss, respectively. The objective‐function values are numerically evaluated through Reynolds‐averaged Navier–Stokes analysis at the design points that are selected through the Latin hypercube sampling method. Using these numerical simulations two surrogates, viz, the RSA and Kriging models, are constructed for each objective function and a hybrid MOEA is applied to obtain the Pareto‐optimal front. For the particular implementation of surrogate models, it is observed that Pareto‐optimal predictions of the RSA model are better than those of the KRG model, whereas the KRG model predicts equally well at the off‐Pareto‐region (region away from the Pareto‐optimal solutions), which is not the case with the RSA model. The local resampling of the Pareto‐optimal front increases the fidelity of the approximate solutions near the Pareto‐optimal region. The ratios of the channel height to the dimple print diameter and of the dimple print diameter to the dimple pitch are found to be more sensitive along the Pareto‐optimal front than the ratio of the dimple depth to the print diameter. The decrease of the ratio of the channel height to the dimple diameter and the increase of the ratio of the dimple print diameter to the pitch lead to greater heat transfer at the expense of the pressure loss, whereas the ratio of the dimple depth to the print diameter is rather insensitive to Pareto‐optimal solutions. Pareto‐optimal solutions at higher values of the Nusselt number are associated with higher values of the pressure loss due to the increased recirculation, mixing of fluid and vorticity generation. Copyright © 2010 John Wiley & Sons, Ltd.     

简单柔性机械臂的构造性逆动力学方法

本文提出了适于简单柔性机械臂的一种基于谱分析的构造性逆动力学方法。数值模拟和实验结果表明，这种方法提供的开环输入力矩，不仅可以实现机械臂特定的刚体转动，而且可以有效地减弱柔性对机械臂精确定位的影响。     

The automatic control of boundary-layer transition — Experiments and computation

In recent years there has been an increasing interest in the control of boundary-layer transition through the use of wall suction. In the current work suction is provided through one or more suction panels situated close to the leading edge of a plate. Experiments show that boundary-layer pressure fluctuation measurements can be used to identify the position of transition. Transition can be maintained at a desired location with minimum power consumption by employing an automatic adaptive feedback control loop which regulates the suction flow rates of two independent suction panels. This can be expressed as a constrained optimization problem. To allow the suction flow rates to be updated, a modified least mean squares algorithm is used within the control loop. Experimental measurements show that the control algorithm allows fast and stable convergence towards the optimum suction distribution for a double suction panel configuration. Numerical simulations have also been performed. The two-dimensional boundary layer was calculated allowing the viscous boundary layer to interact with the inviscid outer flow. Following linear stability theory the spatial growth rates are calculated by solving an Orr-Sommerfeld type eigenvalue problem, with the streamwise location of transition predicted via thee ^N -method. Applying the same optimization strategy as in the experiments, good qualitative agreement between computations and experiments was found. The optimization algorithm has been applied to computer models where the relation between suction flow rates and transition location is described by an empirical analytical function. This shows that the controller can in principle be applied to systems with more than two suction panels.Nomenclature b transition location with zero suction - d desired transition location - e(k) error signal - k iteration index - p rms pressure - p _ref reference rms pressure - r sum of the reference pressure - u streamwise velocity - u _e external velocity - inviscid external velocity - A wave amplitude - F( ) cost function - I identity matrix - N maximum amplification factor - P projection matrix - R Reynolds number - Re Reynolds number based on the boundary-layer thickness - R matrix of weights - Tu turbulence level - vector of suction flow rates - v normal velocity - v _wall suction velocity at the surface - x streamwise coordinates - x _m microphone location - x _T(k) measured transition location - y normal coordinate - y(k) sum of the measured pressures - w(k) noise - plate length - r +i _i - free stream velocity - * displacement thickness - gradient vector - Lagrange multiplier - controller gain - disturbance stream function - disturbance amplitude - wave frequency = complex wave number