生物力学与基因-献给周培源教授诞辰100周年

生物界包罗万象,其中有力的作用,所以有生物力学.自Galileo,Harvey, Boreli, Hooke, Euler,Young等创始以来,生物力学阐明了鸟飞鱼游,人体运动,血液循环,人工脏器等,对人世社会,有所贡献.生物力学的基础是质点力学,传统地用连续体力学的概念来简化.但近年做生物组织在应力的作用下改造的问题,引起了必须更改传统连续体力学的几个公理的问题.我们将仔细讨论这些公理,然后指出新公理存在的理由,是由于基因在细胞里的日常工作.基因不单主宰遗传,变异;并且忙着控制日常生活.不过,现在仅见其端倪.详细的情形,要等将来来阐发了.      

A program-control unit for fatigue-testing machines

Though equipment is available to produce programmed tests on fatigue specimens, such equipment is usually expensive to purchase and complicated to maintain. The tendency, therefore, has been to use constant-amplitude tests and, when designing a member where the loading spectrum is known, to combine the test results according to various theories, thus estimating the fatigue life. This article describes a control unit for an electromagnetic vibrator. The unit, being basically simple in concept, has enabled loads of varying magnitude to be applied to a fatigue specimen in accord with instructions on a punched input tape.     

Mean stress dependent nonlinear hyperelasticity coupled with damage stiffness degradation. A thermodynamical approach

Porous materials, such as geomaterials, exhibit a behaviour dependent on the confining pressure. The aim of this paper is to study the degradation of the elastic stiffness of mean stress dependent materials, due to the deterioration of the microstructure during loading.Continuum damage mechanics offers a framework to model this rigidity deterioration. In addition to the concept of effective stress, a choice has to be made between two widely used hypotheses, the principle of strain equivalence and the principle of equivalent elastic energy, in order to build a complete modelling framework.A mean stress dependent hyperelastic formulation is used to ensure a conservative behaviour, and associated to the two previous damage modelling assumptions, whose effects are compared. This allows for mean stress dependent elasticity to be reproduced, with elastic moduli increasing with mean stress while decreasing with damage.     

Development of a tyre traction testing facility

An experimental setup comprising up of an indoor soil bin, a single wheel tester (SWT), a soil processing trolley, a drawbar pull loading device and an instrumentation unit was developed to perform traction tests in the soil bin to study the effect of soil, tyre and system parameters on the performance of tyres. The design of the single wheel tester was such that the dynamic weight reaction force is equal to that measured statically. It is a simple wheeled device, capable of testing tyres of up to 1.5 m in diameter, vertical force up to 19 kN, net pull up to 7.2 kN, torque up to 5.5 kN m, and speed up to 3.5 km/h.     

Indirect boundary element method for unsteady linearized flow over prolate and oblate spheroids and hemispheroidal protuberances

The indirect boundary element method was used to study the hydrodynamics of oscillatory viscous flow over prolate and oblate spheroids, and over hemispheroidal bodies hinged to a plate. Analytic techniques, such as spheroidal coordinates, method of images, and series representations, were used to make the numerical methods more efficient. A novel method for computing the hydrodynamic torque was used, since for oscillatory flow the torque cannot be computed directly from the weightings. Instead, a Green's function for torque was derived to compute the torque indirectly from the weightings. For full spheroids, the method was checked by comparing the results to exact solutions at low and high frequencies, and to results computed using the singularity method. For hemispheroids hinged to a plate, the method for low frequencies was checked by comparing the results to previous results, and to exact solutions at high frequencies. Copyright © 2004 John Wiley & Sons, Ltd.     

The effect of variations in the creep exponent on the buckling of circular cylindrical shells

The senior author solved the problem of axially symmetrical creep buckling of thin circular cylindrical shells subjected to uniform axial compression. In that analysis the constitutive equation was a power law, and the exponent was taken to be equal to three. The purpose of this work was to extend the solution to a range of values of the creep exponent, n. To cope with the increasing algebraic complexity, a digital computer was employed in two ways: to generate the set of equations symbolically, and then to solve these equations. The computer programs were used to generate numerical solutions for the cases in which n was equal to 3, 5, 7 and 9. Two simple extrapolation techniques were then employed to obtain approximate solutions to the critical time problem for values of n up to 29.     

开放式柔性材料力学实验平台的建设与实践

提出建设以机器编号为标识,以计算机网络通信技术为支撑的柔性力学实验平台,在这个平台上,学生自主进行实验,而教师在后台监控学生实验,既提高了学生实验的积极性,也减轻了教师的工作量,使实验设备得到很好的保护和利用.     

Estimating deflection of a simple beam model using fiber optic bragg-grating sensors

In this paper, we develop a method to estimate the bridge deflection using fiber optic Bragg-grating (FBG) strain sensors. For most structural evaluation of bridge integrity, it is very important to measure the geometric profile, which is a major factor representing the global behavior of civil structures, especially bridges. In the past, because of the lack of appropriate methods to measure the deflection curve of bridges on site, the measurement of deflection has been restricted to just a few discrete points along the bridge, and the measuring points have been limited to the locations installed with displacement transducers. However, by applying classical beam theory, a formula is rearranged to estimate the continuously deflected profile by using strains measured directly from several points. In addition, FBG strain sensors, which are electromagnetic, noise-free and multipoint measurable, are employed to obtain flexural strains more accurately and stably. The regression analysis is performed to obtain a strain function from the measured strain data. Finally, the deflection curve can be estimated by applying the strain function to the formula. An experimental test has also been carried out to verify the developed method.     

扰动法在结构分枝失稳分析中的应用

对结构进行平衡路径的跟踪分析，是全面了解该结构的受力性能所必须进行的一项工作。目前的结构非线性稳定分析技术一般仅对极值点失稳型问题较为有效，而对分枝点失稳型问题则困难较多。对于具有缺陷敏感性的结构，如拱结构、壳体等，在普通荷载作用下，其失稳路径常包含分枝点。文献[1]提出并认为位移扰动法和力扰动法在分析结构分枝失稳时具有很好的效果。本文采用多个不同类型的算例，对扰动法在结构分枝失稳问题中的应用进行了分析比较，表明该方法具有较强的跟踪能力。最后，就扰动法在结构分枝失稳问题中的应用提出几点建议。     

Load-unload response ratio——An interplay between earthquake prediction and mechanics

地震预测是世界性科学难题,地震现象虽然复杂其物理实质倒是明确：地震就是地壳块体的快速剪切脆断,相应地,地震的孕育过程就是震源区介质的损伤、演化,并最终导致破坏的过程,这一过程主要是力学过程,抓住这一点就抓住了问题的物理本质,但是,在研究地震预测时遇到的力学问题和通常的（工程）力学问题有所不同,根据地震问题的特点,紧扣地震孕育过程的物理本质,提出了加卸载响应比这一地震预测新思路.
文中介绍了加卸载响应比理论的基本科学问题,包括如何对地壳加载/卸载,如何选择适当的地球物理参数作为响应量,以及怎样定义加卸载响应比,用实验研究、数值模拟和理论分析3 种基础研究手段,揭示了地震孕育过程中加卸载响应比共同的演化规律：孕震初期加卸载响应比在1 附近涨落,之后上升至峰值点,地震不在峰值点发生,而是在下降过程中发生,从峰值点到地震发生这段滞后时间称为T₂,T₂ 和震级有关,为了预测地震,必须在全国范围内作加卸载响应比的时空扫描,在时空扫描基础上,结合量纲分析,更多地考虑当地的地球物理条件,拟定了全面预测未来地震的时、空、强的思路,回顾了用这种思路进行地震预测实践的历程,     

An Integral Method to Calculate Water Saturation in a Centrifuge Experiment to Determine Capillary Pressure

The centrifuge method is commonly used to determine the capillary pressure of a porous medium, and the original approximating method for data analysis developed by Hassler and Brunner is still being used. Its limitations are, however, not well understood. Application to analyze experiments where one of the assumptions was obviously violated had been given in the literature. While the result appeared to be quite reasonable, it was not clear how close was it to reality. One of the objectives of this paper is to review the assumptions that is required to develop this method, so that the experimental condition in which it is applicable can be established. The other objective is to derive a completely different solution technique to this problem. There is no need to assume that the ratio of the inlet radius to the exit radius of the core to the center of the centrifuge be close to 1. With the freedom from this limitation it is, therefore, possible to construct machines at lower cost and to improve on the data quality by allowing longer cores.     

An icing physics study by using lifetime-based molecular tagging thermometry technique

An experimental investigation was conducted to quantify the unsteady heat transfer and phase changing process within small icing water droplets in order to elucidate underlying physics to improve our understanding of the important micro-physical process of icing phenomena. A novel, lifetime-based molecular tagging thermometry (MTT) technique was developed and implemented to achieve temporally-and-spatially resolved temperature distribution measurements to reveal the time evolution of the unsteady heat transfer and dynamic phase changing process within micro-sized water droplets in the course of icing process. It was found that, after a water droplet impinged onto a frozen cold surface, the liquid water at the bottom of the droplet would be frozen and turned to solid ice rapidly, while the upper portion of the droplet was still in liquid state. As the time goes by, the interface between the liquid phase water and solid phase ice was found to move upward continuously with more and more liquid water within the droplet turned to solid ice. Interestingly, the averaged temperature of the remaining liquid water within the small icing droplet was found to increase, rather than decrease, continuously in the course of icing process. The temperature increase of the remaining liquid water is believed to be due to the heat release of the latent heat during solidification process. The volume expansion of the water droplet during the icing process was found to be mainly upward to cause droplet height growth rather than radial to enlarge the contact area of the droplet on the test plate. As a result, the spherical-cap-shaped water droplet was found to turn to a prolate-spheroid-shaped ice crystal with cusp-like top at the end of the icing process. The required freezing time for the water droplets to turn to ice crystals completely was found to depend on the surface temperature of the test plate strongly, which would decrease exponentially as the surface temperature of the frozen cold test plate decreases.     

Sloshing and energy dissipation in an egg: SPH simulations and experiments

Tall structures, such as towers and bridges, can oscillate at excessive magnitudes when subjected to wind and earthquake loads. Liquid sloshing absorbers can be used to suppress these excessive oscillations by tuning the frequency of the sloshing to the critical frequency of the structure. Sloshing absorbers are simple structures consisting of a partially full container of liquid with a free surface. Tuning ensures that significant amounts of harmful energy can be extracted from the structure to the sloshing liquid. However, there needs to be a rapid means of dissipating this energy to avoid its returning back to the structure (then back to the liquid periodically).A hen׳s egg seems to have evolved to efficiently dissipate energy to protect its embryo using sloshing of its liquid content. Hence, the potential to implement the egg׳s unique properties as a sloshing absorber for structural control, is the main focus of this study. Numerical simulations, using Smoothed Particle Hydrodynamics (SPH), and experimental comparisons are presented in this paper. One objective is to demonstrate the ability of SPH to simulate complex free surface behaviour in three dimensions. Such a tool is then useful to identify different dissipation modes. Effects of fill volume and viscosity on the rate of dissipation, are also investigated.     

Notes on the X-ray Theory

My aim is to better compare thermoelasticity theory with a continuum theory of X-ray observations of crystals, to be described, to make it easier to adapt ideas and techniques from one to the other. To this end, I will present some different ways of formulating the X-ray theory. Also, I will adapt symmetry considerations from the former, partly to analyze problems of a kind not considered in thermoelasticity theory. This revised version was published online in August 2006 with corrections to the Cover Date.     

Application of simple,periodic homogenization techniques to non-linear heat conduction problems in non-periodic,porous media

Often, detailed simulations of heat conduction in complicated, porous media have large runtimes. Then homogenization is a powerful tool to speed up the calculations by preserving accurate solutions at the same time. Unfortunately real structures are generally non-periodic, which requires unpractical, complicated homogenization techniques. We demonstrate in this paper, that the application of simple, periodic techniques to realistic media, that are just close to periodic, gives accurate, approximative solutions. In order to obtain effective parameters for the homogenized heat equation, we have to solve a so called “cell problem”. In contrast to periodic structures it is not trivial to determine a suitable unit cell, which represents a non-periodic media. To overcome this problem, we give a rule of thumb on how to choose a good cell. Finally we demonstrate the efficiency of our method for virtually generated foams as well as real foams and compare these results to periodic structures.     

Sedimentation of an elliptical object in a two-dimensional foam

The sedimentation of an elliptical object in a dry two-dimensional, monodisperse foam is simulated. The calculations are quasi-static, allowing the identification and separation of the elastic, plastic and viscous response of the foam to the motion. In addition to its weight, the forces on the ellipse are due to the network of soap films and the pressures in the bubbles. These give rise to non-zero torque, lift and drag forces, causing the motion of the ellipse to deviate from a vertical path. Highly stretched films are formed in the wake of the ellipse and asymmetry in the flow field, with bubbles moving from the front to the back of the ellipse along only one side, causes it to rotate from a metastable state with its major axis perpendicular to gravity into a stable orientation with its major axis parallel to the direction of gravity. When the orientation is intermediate between these two limits, there is a significant lift force which causes the ellipse to move laterally. A larger, more eccentric, ellipse rotates more quickly.     

Fast prediction of the fatigue behavior of short-fiber-reinforced thermoplastics based on heat build-up measurements: application to heterogeneous cases

Short-fiber-reinforced thermoplastics components for structural applications are usually very complex parts as stiffeners, ribs and thickness variations are used to compensate the quite low material intrinsic stiffness. These complex geometries induce complex local mechanical fields but also complex microstructures due to the injection process. Accounting for these two aspects is crucial for the design in regard to fatigue of these parts, especially for automotive industry. The aim of this paper is to challenge an energetic approach, defined to evaluate quickly the fatigue lifetime, on three different heterogeneous cases: a classic dog-bone sample with a skin-core microstructure and two structural samples representative of the thickness variations observed for industrial components. First, a method to evaluate dissipated energy fields from thermal measurements is described and is applied to the three samples in order to relate the cyclic loading amplitude to the fields of cyclic dissipated energy. Then, a local analysis is detailed in order to link the energy dissipated at the failure location to the fatigue lifetime and to predict the fatigue curve from the thermomechanical response of one single sample. The predictions obtained for the three cases are compared successfully to the Wöhler curves obtained with classic fatigue tests. Finally, a discussion is proposed to compare results for the three samples in terms of dissipation fields and fatigue lifetime. This comparison illustrates that, if the approach is leading to a very relevant diagnosis on each case, the dissipated energy field is not giving a straightforward access to the lifetime cartography as the relation between fatigue failure and dissipated energy seems to be dependent on the local mechanical and microstructural state.     

功能梯度双材料弱/微间断界面的冲击断裂分析

提出强间断、弱间断、微间断和全连续界面的概念与分类,建立功能梯度弹性双材料弱间断 界面冲击断裂问题的力学模型,采用积分变换法推导问题的Cauchy奇异积分方程,并用配 点法求得数值解. 分析表明,弱/微间断性对于FGMs界面裂纹应力强度因子有着重要影响, 而且微间断性是优于弱间断性的一种界面力学性能连接关系. 以FGMs界面某一侧 的力学性能函数在界面处的Taylor展开式的低阶项作为界面另一侧的力学性能函数,便可 以使FGMs界面成为``微间断'界面. 界面的一阶微间断对应力强度因子的减小作用较为明 显,而高阶(二阶及以上)微间断对应力强度因子的影响较小. 减小界面的弱间断程度或使 FGMs界面具备``微间断性',都将利于提高功能梯度双材料界面抗冲击断裂能力,在一定 程度上达到界面增韧的目的.     

浅谈复合型工程力学专业人才的培养

作为一名一线的力学教学科研人员, 通过自身的教学体会及深入相关行业了解和调研, 结合 实际情况, 客观地提出当前工程力学专业人才培养中存在的问题, 并思考当前高校工程力学 专业人才培养中需要做出的改革.     

区域无线电定位系统基站的最优几何布局

地面基站的几何布局方式是影响区域无线电定位系统定位精度的重要因素,水平几何精度递变因子（HDOP）是衡量地面基站几何布局对系统定位精度影响的基本指标。为求解区域无线电定位系统基站的最优几何布局,给出了区域无线电定位系统HDOP的计算方法,建立了求解区域无线电定位系统地面基站最优几何布局的数学模型;提出了应用惩罚函数对数学模型进行求解的方法,将约束优化问题转化为无约束优化问题;并针对求解过程中目标函数难于求导的问题,提出了采用模拟退火算法进行全局寻优,使求得的结果逼近全局最优解。仿真结果表明,这种结合惩罚函数法和模拟退火算法进行全局寻优的方法,对于区域无线电定位系统求解地面基站的最优布局方案具有工程实用价值。