瓷修复体界面断裂行为的模拟实验研究

本文利用云纹干涉法和云纹干涉－－有限元混合法，对瓷修复体的模拟双材料模型界面断裂问题进行了实验研究。用云纹干涉和数字错位云纹干涉法测量带边裂纹的双材料四点简支梁在剪切作用下界面表面的剪应变分布及界面两侧局部表面的位移场，实验表明，由于界面两两侧材料力学性质不同，表现出界面剪切断裂问题的非称性和裂尖附近复合型断裂的特点；用云纹干涉法和有限元法相结合的混合法对粘接界面角点应力奇异性进行研究，并对角点附近应力应变场作了分析，得到了应力奇异指数与边界楔角，载荷的关系，证明了用界面应力强度因子Kf来描述界面端部区域应力分布的公式，并得到了双材料界面端部区域的应力应变分布情况。本文的实验结果为进一步研究口腔金瓷修复体界面的优化设计提供了基础，同时也说明云纹干涉法对于双材料界面断裂行为的研究是有效的。     

Investigation of the Nonlinear Behavior of Beams*

ABSTRACT A system of equations that describes the nonlinear behavior of teams is presented. The differential equations are solved by using the Galerkin method, and a system of nonlinear algebric equations is obtained. A method to deal with discontinuities in structural properties and load distribution along a beam is presented. The derivation includes general expressions for gravity-type loads. Two methods of calculating the resultant moment along the beam are described. One method incorporates differentiation of the expressions for displacements, while the second method is based on integration of loads along the beam. The numerical model is used in order to investigate displacements and moment resultants along a cantilevered beam. Different problems that are associated with nonlinear behavior are presented and discussed. Good agreement is obtained between theoretical and experimental results.     

考虑尺度效应的微梁静态弯曲数值分析

基于修正偶应力理论及表面弹性理论,本文提出了一种新的双曲线剪切变形梁模型,用于均匀微尺度梁的静态弯曲分析。该理论可以直接利用本构关系获得横向剪切应力,满足梁顶部和底部的无应力边界条件,避免了引入剪切修正因子。根据广义Young-Laplace方程建立了梁的内部与表面层的应力连续性条件,单一的变量场可以描述梁的位移模式。通过在位移场中考虑表面层厚度以及表面层的应力连续条件,可以使新模型能够更准确地预测微尺寸和表面能相关的尺度效应。通过Hamilton原理推导出了梁的控制方程和边界条件。应变能除了考虑经典弹性理论,还要考虑微结构效应和表面能。Navier-type的解析解适用于简支边界条件,而基于拉格朗日插值的微分求积法(DQEM)可以研究在不同边界条件下的力学响应。把该数值解与Navier方法得出的解析解作了对比,得出:微尺度梁在考虑表面能或微尺寸效应、不同载荷和梁高变化下的响应一致;当不考虑微结构相关性和表面能效应时,该模型退化为经典的欧拉梁模型。     

Nondestructive damage detection in structures using changes in compliance

A methodology to identify damage in a structure is presented in this paper. The method utilizes a new form of damage index based on the changes in the distribution of the compliance of the structure due to damage. The changes in the compliance distribution are obtained using the mode shapes of the pre-damaged and the post-damaged state of the structure. The validity of the method is demonstrated using numerically generated data from beam structures and experimental data from a free–free beam structure with inflicted damage. In the numerical and experimental examples, the damage identification performance of the proposed method is compared with that of the existing strain-energy-based method. The results of the numerical and experimental studies indicate that the proposed compliance-based damage index method can be used in damage identification of the structure.     

Modeling and bifurcation analysis of the centre rigid-body mounted on an external Timoshenko beam

Ａｓｗｅｋｎｏｗ,ｉｔｉｓｔｈｅｃｈａｒａｃｔｅｒｉｓｔｉｃｏｆｔｈｅｆｌｅｘｉｂｌｅｍｕｌｔｉｂｏｄｙｓｙｓｔｅｍｔｈａｔｔｈｅｒｉｇｉｄｍｏｔｉｏｎｉｓｃｏｕｐｌｅｄｗｉｔｈｔｈｅｆｌｅｘｉｂｌｅｄｅｆｏｒｍａｔｉｏｎ［１～３］．Ｔｈｅｒｅｆｏｒｅ,ｆｏｒｔｈｅｆｌｅｘｉｂｌｅｍｕｌｔｉｂｏｄｙｓｙｓｔｒｍ,ｉｔｉｓｎｅｃｅｓｓａｒｙｔｏｉｎｖｅｓｔｉｇａｔｅｔｈｅｒｅｇｕｌａｒｏｆｔｈｅｄｅｓｔａｂｉｌｉｚａｔｉｏｎａｎｄｂｉｆｕｒｃａｔｉｏｎｏｆｉｔｓｃｏｎｆｉｇｕｒａｔｉｏｎｂｅｆｏｒ…     

Pull-in analysis of electrically actuated viscoelastic microbeams based on a modified couple stress theory

A new beam model is developed for the viscoelastic microbeam based on a modified couple stress model which contains only one material length scale parameter. The governing equations of equilibrium together with initial conditions and boundary conditions are obtained by a combination of the basic equations of modified couple stress theory and Hamilton’s principle. This new beam model is then used for an electrically actuated microbeam-based MEMS structure. The dynamic and quasi-static governing equations of an electrically actuated viscoelastic microbeam are firstly given where the axial force created by the midplane stretching effect is also considered. Galerkin method is used to solve above equation and this method is also validated by the finite element method (FEM) when our model is reduced into an elastic case. The numerical results show that the instantaneous pull-in voltage, durable pull-in voltage and pull-in delay time predicted by this newly developed model is larger (longer) than that predicted by the classical beam model. A comparison between the quasi-static model results and the dynamic model results is also given.     

Boundary element method for model analysis of 2-D composite structure

The boundary element method is used for the modal analysis of free vibration of 2-D composite structures in this paper. Since the particular solution method is used to treat the terms of body forces (inertial forces) in the equation of motion, only static fundamental solutions are needed in solving the problem. For an isotropic cantilever beam, the numerical results obtained by using the BEM presented in this paper are in good agreement, with, those of using FEM or other BEM, but this BEM can also be used to analyze problems for anisotropic materials. For simply supported composite laminated beams, the comparisons of the numerical reslts obtained by this method with the analytical results obtained by 1-D laminated beam theory indicate that if the ratio of length/thickness is greater than 20, the results of the two methods are in good agreement, but if the ratio of length/thickness is less than 20, big errors will occur for 1-D laminated beam theory.     

大变形柔性多体系统非线性动力学数值仿真与实验研究

提出了一种作大范围运动柔性梁的非接触动态测试技术.在基于位移的柔性多体系统几何精确建模及非线性有限元分析技术的基础上,利用EAGLE-500运动分析系统及其相应的分析软件对作大范围运动钛合金柔性梁作了实验研究,并且利用之前提出的几何精确梁理论进行数值仿真.数值仿真结果与实验结果完全吻合,验证了作者所提的几何精确梁理论及...     

腰椎力学分析的数值模拟与实验研究

采用数值模拟和实验测试技术对两种不同内固定法的腰椎模型进行应力和变形分析,基于CT图像建立L4-S1的三维数值模型,经ANSYS计算分析得出五种工况下的终板应力值;在实验中采用了一种薄膜压力测试传感器结合图像处理的方法,提高测试椎间盘压力分布的精度;同时采用数字图像相关技术对腰椎骨上下关节突在承载情况下的空间位移进行了测量,获得了腰椎间盘（L3-L4）在承受轴压、前屈后伸和侧弯情况下的压力分布,以及对应的关节突的位移迹线。结果表明：本研究采用的数值分析技术和实验开发的测试技术可操作性强,精度满足要求,有望在类似的生物力学分析中得到应用。     

屈曲型数字微镜结构的机电动力学分析

设计一种基于压杆屈曲原理的数字微镜装置，针对该装置建立机电动力学模型，应用龙格-库塔法和有限元数值方法求解其动力学模型中的隐函数。通过对该模型的完整机电动力学微分方程数值求解所作的运动仿真，得到了微镜基板的连续角位移曲线，表明该装置具有良好的光束稳定特性。     

A combined experimental and numerical method for the solution of generalized elasticity problems

This work presents the method for the investigation of three-dimensionally stressed bodies with arbitrary shape which are under the action of an outside system of arbitrary forces. The combined method is based on syntheses of photoelastic experimental methods (other experimental methods may also be used) and digital methods of discrete analysis. Experimental procedures are used for defining superfluous boundary conditions. The boundary-value problem with such boundary conditions is solved by numerical methods. This approach qualitatively changes the very essence of experimental methods and essentially widens their range. It reduces the amount of measurements required and, at the same time, allows one to obtain complete stress fields throughout a body in a short time. In comparison with numerical methods, the combined method increases the accuracy of problem solutions and, at the same time, reduces the time required for complete investigations.     

基于辛弹性力学解析本征函数的有限元应力磨平方法

在实际工程结构的结构强度与优化等力学数值分析中,应力计算结果的精度是非常重要的。有限元法是得到最广泛应用的一类数值方法,并形成了众多通用的有限元程序系统。这些程序系统采用的几乎都是基于最小总势能的位移法,虽然其分析给出的有限元位移场具有较高的精度,但所得到的有限元应力场的精度较位移场大大降低。基于极坐标辛对偶体系所提供的平面弹性力学的解析辛本征展开解,并借用有限元程序系统所得到的节点位移,本文提出了一个应力分析的改进方法。数值结果表明,本方法给出的应力分析精度得到大幅提高,并具有良好的数值稳定性,可用于有限元程序系统的后处理,以提高应力尤其是关键区域应力的分析精度。     

一维六方准晶层合简支梁自由振动与屈曲的精确解

伪Stroh型公式能够将多场耦合材料的控制方程转化为线性特征系统来求解,从而获得多层结构简支边界条件的精确解.本文利用伪Stroh型公式,研究一维六方准晶层合简支梁的自由振动和屈曲问题,通过传递矩阵法,获得准晶层合梁自由振动固有频率与临界屈曲载荷的精确解.通过与已有梁的剪切变形理论结果比较,验证了本文伪Stroh型公式的正确性和有效性.通过数值算例,分析由两种不同准晶材料组成的三明治层合梁的叠层方式、高跨比、层厚比及层数对梁的固有频率、临界屈曲载荷及其模态的影响规律.结果表明,叠层顺序和梁的高跨比、层厚比对准晶层合梁的自由振动固有频率和临界屈曲载荷有很大影响,可通过调整梁的几何尺寸和叠层顺序得到准晶层合梁的最佳固有频率和临界屈曲载荷.本文给出的精确解可为工程上研究准晶梁的各种数值解法和实验方法提供理论参考.     

Photoelastic determination ofK I andK II : A numerical study on experimental data

The most efficient photoelastic methods to obtain stress intensity factors are those based on stress functions series expansions. The coefficients of these expansions are fitted to the experimental isochromatic pattern using an overdeterministic Newton-Raphson least squares method. In this paper, a study has been carried out to analyze the influence on the results of several numerical and experimental factors. It is shown that accurate values of the stress intensity factorsK _I andK _II can be obtained by following some recommendations given in the text and summarized in the conclusions at the end of the paper.     

基于单元的子区间摄动有限元方法研究

针对区间相关性导致区间扩张的问题,探讨了区间数之间的相关性并给出了降低区间扩张的子区间摄动方法。文中给出了基于单元的子区间摄动有限元计算公式和子区间划分数目的近似计算公式,同时文中讨论了区间有限元计算精度问题,给出了可提高计算效率的一些措施。对桁架结构和平面应力问题梁结构算例分析结果表明文中方法可以达到一定的计算精度,并且是合理可行的。     

均布荷载作用下悬臂磁电弹性梁的解析解

对磁电弹性平面问题进行了研究，给出了用拟调和位移函数表达的通解，进而以试凑法按平面应力问题推导出了均布荷载作用下悬臂磁电弹性粱的解析解，所得解有易于理解、便于校对、形式统一简洁的特点。本文还将计算结果与压电材料和弹性材料相应结果进行了分析、比较，为验证各种数值计算方法提供了参考依据。     

Dynamic Analysis of Multi-Beam MEMS Structures for the Extraction of the Stress-Strain Response of Thin Films

Freestanding MEMS structures made of two long connected beams from different materials are fabricated and released in order to extract the stress-strain properties of thin films. The first material, named actuator, contains a high internal tensile stress component and, when released, pulls on the other beam. The strain in the beams is calculated based on the measurement of the displacement with respect to the reference configuration using scanning electron microscopy. The stress is estimated using two different methods. The first method, already reported, is based on the displacement of the actuator and the knowledge of its internal stress. The method which constitutes the novelty of the present study is based on the dynamic analysis of the multi-beam structures, and the determination of the stress value that corresponds to the measured resonance frequencies. The dynamic analysis is performed via two different methods: (i) the modified Rayleigh–Ritz technique and (ii) the Euler–Bernoulli beam dynamics. Results are provided for palladium thin films which deform plastically and for monocrystalline silicon thin films, exhibiting a purely elastic behavior. The results show the higher accuracy of the dynamic measurements for the estimation of the stress compared to the static method. The dynamic measurements also show that the Rayleigh–Ritz technique tends to give a higher bound for the resonance frequencies compared to the Euler–Bernoulli technique. This dynamic method extends the potential of this on-chip material testing technique which can also be adapted to stress controlled sensors applications.     

基于p型有限元法和围线积分法计算复合型应力强度因子

传统有限元法由于采用低阶插值计算应力强度因子时,需要划分的网格数较多,收敛速度较慢,得到的应力强度因子精度不足。p型有限元法在网格确定时通过增加插值多项式的阶数来提高计算精度,具有网格划分少、收敛速度快、精度高、自适应能力强等特点。本文采用基于p型有限元法的有限元计算软件StressCheck计算得到应力场和位移场,并由围线积分法导出混合型应力强度因子(SIFs)。通过几个经典算例,分析了围线的选择对计算精度的影响,计算了不同裂纹长度、不同裂纹角度和裂纹在应力集中区域不同位置时的应力强度因子。并将数值结果、理论解与文献中其他数值计算方法所得的部分结果进行了对比分析,结果表明自由度数不大于7000时,导出的应力强度因子相对误差最大不超过1.2%,数值解表现出较高的精度及数值稳定性。     

Nonlinear resonances in an axially excited beam carrying a top mass: simulations and experiments

In this paper, nonlinear resonances in a coupled shaker-beam-top mass system are investigated both numerically and experimentally. The imperfect, vertical beam carries the top mass and is axially excited by the shaker at its base. The weight of the top mass is below the beam’s static buckling load. A semi-analytical model is derived for the coupled system. In this model, Taylor-series approximations are used for the inextensibility constraint and the curvature of the beam. The steady-state behavior of the model is studied using numerical tools. In the model with a single beam mode, parametric and direct resonances are found, which affect the dynamic stability of the structure. The model with two beam modes not only shows an additional second harmonic resonance, but also reveals some extra small resonances in the low-frequency range, some of which can be identified as combination resonances. The experimental steady-state response is obtained by performing a (stepped) frequency sweep-up and sweep-down with respect to the harmonic input voltage of the amplifier-shaker combination. A good correspondence between the numerical and experimental steady-state responses is obtained.     

基于一种广义梁理论的弹性地基FGM简支梁自由振动解析解

基于一种扩展的n阶广义剪切变形梁理论(n-GBT),应用Hamilton原理,建立了以轴向位移、横向位移及转角为未知函数的Winkler-Pasternak弹性地基功能梯度材料(FGM)梁的自由振动方程,采用Navier法获得了弹性地基FGM简支梁自由振动的精确解。与多种梁理论预测结果进行比较,讨论并给出了GBT阶次n的理想取值;分析了梯度指标、跨厚比及地基刚度对FGM梁频率的影响。结果表明:本文方法有效且适用范围广,若采用高阶剪切梁理论模型,宜取n≥3的奇数;FGM梁的自振频率随材料梯度指标的增大而减小;随跨厚比的增加而增大,但当跨厚比大于20,跨厚比增加对频率的影响很小;随地基刚度的增加而增大,地基刚度足够大时,频率趋于收敛。