不同模量弹性问题理论及有限元法研究进展

随着科学技术的日益发展，对材料力学性质的研究提出了更高的要求，研制新型的材料以及挖掘材料自身特性的潜力，已成为新的研究动向．简述了不同模量弹性问题理论及其有限元法的研究与发展．利用等效的概念，对Ambartsumyan有限元计算模型、Jones有限元计算模型、张允真等有限元计算模型、叶志明等计算模型进行改进和探讨．通过不同模量弹性理论及其有限元方法在实际工程构件问题分析中的应用表明，若沿用相同模量弹性理论或有限元对有关问题进行计算，其结果将与采用不同模量模型材料所得的刚度和强度有较大的偏差。     

界面对多相介质平均弹性性能和应力场的影响

在建立多相介质夹杂问题的自治模型基础上，应用自治有限元法研究了多相介质的平均弹性性能问题，通过数例讨论了在不同体分比下，界面对多相介质总体平均弹性性能的影响，结果表明，界面厚度对多相介质横向平均剪切模量和横向平均体积模量有明显的影响，而对纵向平均拉伸模量影响不大。为了进一步研究界面破坏特征，本文还分析了多相介质界面区域的应力场。     

基于光纤传感器的固体杨氏模量与泊松比测量方法研究

提出了一种基于光纤传感器的测量固体材料杨氏模量与泊松比的方法,给出了Sagnac光纤干涉仪的环形光纤作为点传感器检测声波及用于声速测量的原理,推导出固体中声波速度与材料的杨氏模量及泊松比的关系式.以光纤Sagnac干涉仪作为传感器贴附于三角铸铁导轨表面,分别测出声表面波与纵波的速度,根据理论推导,得到了铸铁的杨氏模量与泊松比.实验结果与理论数值基本吻合,表明了此测量方法的可行性,为固体材料性能评估提供了一种新的方法.     

DETERMINATION OF EFFECTIVE MODULI FOR FOAM PLASTICS BASED ON THREE PHASE SPHEROIDAL MODEL

The relations of bulk modulus,shear modulus,Young's modulus and the Poisson's ra-tio with porosity of foam plastics are determined by a three phase spheroidal model commonly used inComposite Mechanics.The results are compared with those using differential scheme.It is shownthat the material properties derived from the present model normally are larger than those obtained bydifferential scheme for foam plastics with identical porosity.The differences in shear moduli andYoung's moduli obtained by the two methods are small,but they are larger for bulk moduli of incom-pressible matrix and Poisson's ratios.The Young's moduli of high density foam plastics derived by thepresent model agree better with experimental ones.     

基于三相球模型确定泡沫塑料有效模量

通过复合材料三相球模型确定泡沫塑料的体积模量、剪切模量、杨氏模量和泊松比等材料等参数与材料孔隙比的关系，并将该文结果同微分法得到的结果进行了对比；结果表明：三相球模型确定的材料参数在相同孔隙比情况下，一般都高于微分法确定的结果，两种方法得到的有效剪切模量和杨氏模量差异较小，而体积模量和泊松比则差异较大，然而，由三相球模型确定的泡沫塑料杨氏模量仍可与实验较好地符合。此外，该文结果也同现有的理论及实验     

Optimal shape and non-homogeneity of a non-uniformly compressed column

The best possible distribution of Young's modulus and/or the cross-sectional area is found for a column which, for a given volume and length, carries the maximum possible axial loads which are non-uniformly distributed along its length and concentrated at the end-points. The column is elastically clamped at one end and free at the other, where the concentrated axial load is applied. The design variables are subject to upper and lower bounds. Sufficient optimality conditions are derived for a given function to be a solution of the optimization problem. The procedure to determine the optimal solutions is described. Numerical results are obtained by employing an iterative computational technique.     

Method for identification of the filled polymer material relaxation kernel in millisecond time range

The filled polymer materials exhibit viscoelastic properties in a wide time range including the millisecond range (∼10⁻²–10 ms) characteristic of different shock loadings of structures made of these materials. We propose a method for the identification of the filled polymer material relaxation kernel in the millisecond time range; this method is based on a shock loading test of a cylindrical sample made of this material. In this test, the disk indenter acceleration is measured by using a piezotransducer. The test scheme does not impose any rigid constraints on the sample dimensions. In particular, it is possible to use samples of typical dimensions of the order of 10 cm, for which the conditions that the sample material is representative of the structure material are necessarily satisfied. The relaxation kernel parameters are identified by numerical minimization of the theoretically predicted indenter velocity deviation from the velocity-time dependence obtained by integrating the acceleration transducer readings. The minimization problem is solved by using a genetic algorithm. The problem of theoretical prediction of the indenter velocity is solved numerically by using a reduced computational scheme whose parameters are chosen from the minimum condition for the deviation from the prediction obtained in the framework of the detailed computational scheme. The use of the reduced computational scheme permits decreasing the computational costs by 3–4 orders of magnitude compared with the detailed computational scheme, which is a necessary condition for the practical applicability of the genetic algorithm in identification problems. We present examples of relaxation kernel identification in the range of 0.1–10ms from the results of the test where the disk indenter raised to the height of 1m falls on the sample end surface.     

梯度弹性基础上正交异性薄板的屈曲分析

基于双参数弹性基础模型,研究了梯度弹性基础上正交异性薄板的屈曲问题. 首先,基于能量法与变分原理,给出了梯度弹性基础上正交异性薄板的屈曲控制方程,并得到了梯度弹性基础刚度系数K1 与K2的计算式;进而,通过将位移函数采用三角函数展开的方法,给出了单向压缩载荷作用下、四边简支正交异性弹性基础板屈曲载荷的计算式;在算例中,通过将该文的解退化到单纯的正交异性板,并与经典弹性解比较,证明了理论的正确性;最后,求解了弹性模量在厚度方向上呈幂律分布的梯度基础上的薄板屈曲问题,分析了基础上下表层材料弹性模量比与体积分数指数对屈曲载荷的影响.     

适用于几何非线性气动弹性分析的结构动力学降阶模型方法研究

几何非线性是壁板颤振和大展弦比机翼气动弹性等问题的一个主要特征,在进行数值仿真分析时往往需要采用商业非线性有限元求解器,存在计算量大和耦合迭代策略不易控制等问题。本文发展了一种适用于几何非线性的结构动力学降阶模型(CSD-ROM),利用广义坐标的非线性多项式表征非线性内力,采用参数识别方法获取多项式系数,并通过增加额外的线性模态来改善模型预测精度。基于此方法,分别针对壁板颤振、切尖三角翼的CFD/CSD-ROM非线性颤振问题开展了时域响应分析。计算结果表明,通过CSD-ROM计算出的壁板颤振速度为590 m/s,颤振频率为174 Hz,与有限元结果误差分别为0.8%和1.7%。马赫数0.879时切尖三角翼的颤振动压预测结果为2.25 psi,与非线性有限元相比的误差为3.8%。本文采用的非线性和线性模态基底组合方法,在保证计算精度的基础上可有效降低训练样本数量,一定程度上可替代非线性有限元开展气动弹性分析。     

非晶材料弛豫过程加杨氏模量的唯象研究

从晶态材料弛豫过程中模量的唯象描述的基本思想出发，同时考虑到近年来对非晶材料弛豫规律的种种研究，将传统的构造力学参数模型作为弛豫单元的作法与现在流行的双能级模型统一起来，对非晶材料弛豫过程中的杨氏模量进行了唯象处理，得到了与以往不同的弛豫规律。     

一种可用于机翼优化设计中的气动噪声预测模型

发展了一种可用于翼型/机翼外形设计中的气动噪声快速预测方法。相较于传统的半经验噪声预测方法,该方法以两方程非线性k-ε湍流模型模化雷诺应力的雷诺平均方程为背景,考虑了升力系数、三维流动效应以及机翼几何参数等因素对后缘噪声的影响。而相对于直接数值模拟或声类比拟方法,该方法虽不能准确预测噪声强度,但其计算量小,能给出不同翼型/机翼的相对总声压级,以及总声压级随升力系数的变化情况,易于应用于翼型/机翼气动外形优化设计中。通过计算分析二维NACA0012翼型几何参数或来流状态的改变所带来的气动噪声差异,与ANOPP软件及Brooks等计算结果进行对比,验证了该模型的可靠性。最后,计算分析二维、三维翼型/机翼气动噪声,凸显该方法在翼型/机翼气动外形优化设计中的应用价值。     

Computation of cascade flutter by uncoupled and coupled methods

A computational method for flutter prediction of turbomachinery cascades is presented. The flow through multiple blade passages is calculated using a time-domain approach with coupled aerodynamic and structural models. The unsteady Euler/Navier-Stokes equations are solved in quasi-three-dimensions using a second-order implicit scheme with dual time-stepping and a multigrid method. A structural model for the blades with bending and torsion degrees of freedom is integrated in time together with the flow field. Information between structural and aerodynamic models is exchanged until convergence in each real-time step. Computational results for a cascade are presented and compared with those obtained by the conventional energy method and with experimental and numerical data by other authors. Significant differences are found between the coupled and uncoupled methods at low mass ratios. A transonic test case with strong nonlinear phenomena is investigated with the fluid-structure coupled method. Results for inviscid flow are compared with results of Navier-Stokes computations.     

Effective elastic properties of porous materials: Homogenization schemes vs experimental data

In this paper, we focus on the prediction of elastic moduli of isotropic porous materials made of a solid matrix having a Poisson's ratio vm of 0.2. We derive simple analytical formulae for these effective moduli based on well-known Mean-Field Eshelby-based Homogenization schemes. For each scheme, we find that the normalized bulk, shear and Young's moduli are given by the same form depending only on the porosity p. The various predictions are then confronted with experimental results for the Young's modulus of expanded polystyrene (EPS) concrete. The latter can be seen as an idealized porous material since it is made of a bulk cement matrix, with Poisson's ratio 0.2, containing spherical mono dispersed EPS beads. The Differential method predictions are found to give a very good agreement with experimental results. Thus, we conclude that when vm=0.2, the normalized effective bulk, shear and Young's modulus of isotropic porous materials can be well predicted by the simple form (1 − p)² for a large range of porosity p ranging between 0 and 0.56.     

有限厚度模型校准二维MoS2的Bimodal-AFM杨氏模量测量

二维材料因其独特的晶体结构、新奇的物理特性和优异的力学性能, 在微纳机电系统、柔性电子器件等诸多领域有着广阔的应用前景. 弹性模量是二维材料的基本力学特性参量之一, 对其器件应用及应变调控有重要影响. 受限于二维结构和原子级厚度特征, 难以实现二维材料弹性模量的精确测量. 双模原子力显微镜的振幅调制-频率调制模式是一种高效测量二维材料杨氏模量的方法, 但刚性衬底对测量结果的影响不可忽视. 本工作通过双模原子力显微镜直接测得衬底与二维硫化钼的杨氏模量分布图, 并基于有限厚度模型对衬底效应进行修正, 得到了样品的本征杨氏模量值. 利用第一性原理计算得到了二维二硫化钼的弹性系数和杨氏模量, 对比发现实验和计算结果相当. 这说明双模原子力显微镜测量是一种可靠的二维材料杨氏模量直接测试方法, 且该方法无需制备悬空二维材料等繁琐步骤, 避免了常规测试中的不足. 本工作为大面积二维材料薄膜力学性能的程序化测试分析以及高通量力学实验数据的统计分析提供了可靠的实验基础.      

带制退器的膛口射流噪声数值模拟与实验研究

为了研究膛口装置对膛口噪声气动特性的影响,对带膛口制退器的某小口径武器的膛口射流噪声进行了数值模拟和实验研究。采用计算流体力学CFD (computational fluid dynamics)-计算气动声学CAA (computational aeroacoustics)耦合算法对膛口噪声进行数值模拟,即对膛口流场进行瞬态CFD模拟,获取流场数据,然后利用所得到的结果采用声学方程模拟声源信息求解声场。基于数值模拟结果,分析了膛口流场变化及噪声的指向性分布,并与实验结果进行了对比。研究表明：膛口制退器的安装改变了膛口流场结构,影响了膛口射流噪声的指向性分布。计算结果与实验结果的误差小于9%,验证了该计算方法的可行性。研究结果可为膛口射流噪声的预测及膛口制退器的设计提供一定的参考。

     

基于近场动力学的单晶冰弹性各向异性数值模拟方法

天然冰材料在变形与破坏行为上的各向异性特征是冰与结构相互作用中产生复杂载荷过程的关键诱因, 而天然冰各向异性的根源则在于单晶冰的各向异性. 目前, 学术界针对单晶冰各向异性的数值模拟方法研究仍较为缺乏. 为了准确再现天然冰材料的特殊力学性质, 本文基于近场动力学理论, 提出了一种单晶冰弹性各向异性的数值模拟方法. 该方法的核心思想是将单晶冰杨氏模量沿不同加载方向的变化规律引入到近场动力学力密度向量的影响函数中. 以前人实验测试得到的杨氏模量值为参考, 通过开展与C轴呈0°, 45°和90°三个加载方向的单晶冰单轴压缩数值模拟实验, 提出了针对该影响函数的修正和辅助参数标定方法, 最终在15°, 30°, 60°和75°等其他四个加载方向进行了验证. 结果表明: 本文提出的针对影响函数的修正与参数标定方法, 能够较为便捷地找到数值模型杨氏模量与参考杨氏模量相一致的影响函数最优解, 即本文提出的基于影响函数的近场动力学数值模拟方法, 能够合理、准确地模拟单晶冰的弹性各向异性行为. 本文研究成果可为后续多晶冰各向异性数值模拟方法的建立提供基础性参考.      

计算地球物理流体力学中的非线性不稳定问题(续)

四、算子差分方程和谱方法的非线性稳定问题§10 计算稳定性和能量守恒性与算子非负性的关系本节引进算子方程,并推广§6和§8的一些结果。对于非定常流体力学方程,往往都可以化为如下形式的“发展方程”:其中F=F(x,t)是待求函数,是一个非线性算子,X=X(x_1,…,x_k)是空间坐标,k是空间维数,t是时间坐标。例如,对于正压原始方程     

Modeling and prediction of bulk properties of open-cell carbon foam

The emerging ultralightweight material, carbon foam, was modeled with three-dimensional microstructures to develop a basic understanding in correlating microstructural configuration with bulk performance of open-cell foam materials. Because of the randomness and complexity of the microstructure of the carbon foam, representative cell ligaments were first characterized in detail at the microstructural level. The salient microstructural characteristics (or properties) were then correlated with the bulk properties through the present model. In order to implement the varying anisotropic nature of material properties in the foam ligaments, we made an attempt to use a finite element method to implement such variation along the ligaments as well as at a nodal point where the ligaments meet. The model was expected to provide a basis for establishing a process-property relationship and optimizing foam properties.The present model yielded a fairly reasonable prediction of the effective bulk properties of the foams. We observed that the effective elastic properties of the foams were dominated by the bending mode associated with shear deformation. The effective Young's modulus of the foam was strongly influenced by the ligament moduli, but was not influenced by the ligament Poisson's ratio. The effective Poisson's ratio of the foam was practically independent of the ligament Young's modulus, but dependent on the ligament Poisson's ratio. The effective Young's modulus of the carbon foam was dependent more on the transverse Young's modulus and the shear moduli of the foam ligaments, but less significantly on the ligament longitudinal Young's modulus. A parametric study indicated that the effective Young's modulus was significantly improved by increasing the solid modulus in the middle of the foam ligaments, but nearly invariant with that at the nodal point where the ligaments meet. Therefore, appropriate processing schemes toward improving the transverse and shear properties of the foam ligaments in the middle section of the ligaments rather than at the nodal points are highly desirable for enhancing the bulk moduli of the carbon foam.     

Fuzzy variational principle and its applications

Linear and non-linear peaky fuzzy numbers and their arithmetic operations are constructed for the analysis of engineering structures with fuzzy characteristic quantities. Fuzziness of the corresponding quantities is consistently incorporated into the functional of the total potential energy. A set of deterministic recursive equations is obtained as the alternative expressions of the fuzzy variational principle by means of the second-order perturbation technique. The fuzzy Ritz method and the fuzzy finite element method are presented as the applications of the fuzzy variational principle. Accordingly, the roundabout procedures frequently used in the formulations of the fuzzy finite element method are avoided. A benchmark problem of a bending beam with fuzzy Young's modulus under fuzzy external loading is solved by the developed fuzzy numerical methods. Numerical examples show that results determined by these two fuzzy methods are both little conservative, and are in good agreement with those obtained by the analytical method. Moreover, the fuzzy Ritz method or the fuzzy finite element method can provide more valuable information than the conventional deterministic methods.     

聚氨酯泡沫塑料压缩杨氏模量的理论预测

通过微分法导出了泡沫塑料剪切模量和体积模量所满足的微分方程组，再利用联系泡沫塑料泊松比和孔隙比的Ｋｅｒｎｅｒ－Ｒｕｓｃｈ经验关系及泡沫塑料弹性常数间满足足的关系，在基体材料不可压缩的假设下，确定了泡沫塑料的杨氏模量。本文针对几种密度的泡沫塑料，分别对它们的杨氏模量进行了理论预测和实验测定，结果表明：理论预测的模量在较高密度下与实验符合的很好，在低密度下也给出相当好的近似值。此外，本文的结果同其他理