Evaluation of the substrate effect on indentation behavior of film/substrate system

A method to evaluate the substrate effect quantitatively in film indentation is proposed. For the thin film deposited on the substrate, the power function relationship is used to describe the loading curve of the film indentation behavior. The loading curve exponent of the power function which is the fitting parameter can reflect the substrate effect quantitatively. The finite element method is used to simulate the nanoindentation process of the film/substrate system. The loading curve exponent can be obtained from the simulation results. A substrate effect factor based on the loading curve exponent is defined to characterize the effect of the substrate on film indentation. Meanwhile, the dimensionless function of the loading curve exponent related with the material properties and indentation depth is obtained. The results can be helpful to the measurement of the mechanical properties of thin films by means of nanoindentation.     

光学系统透镜热形变特性仿真研究

光学系统透镜组件在环境温度变化时发生的形变,特别是透镜的厚度、面形以及透镜间空间位置的变化,对系统成像质量有着不可忽视的影响。分析计算了某型导引头光学系统透镜组件在温度载荷下的形变特性。计算结果表明,当环境温度从20℃降低到-40℃时,透镜表面位移呈现出显著的梯度变化,7个透镜中表面的最大位移达0.05mm,透镜厚度的变化率超过0.05%,透镜间距离的最大变化率为3.46%。利用数据重构技术,计算结果可以作为透镜镜面形变分析的参考数据。     

Influence of intraocular pressure on geometrical properties of a linear model of the eyeball: Effect of optical self-adjustment

In general, visual acuity does not change with variations in intraocular pressure. Experiments in vitro as well as our clinical findings lead us to hypothesise that the eyeball could possess certain mechanical properties enabling it to automatically produce a sharp image on the retina despite variations in intraocular pressure. Previously reported simple biomechanical models of the eye did not confirm this hypothesis. Here, we propose a generalised mechanical model of the eyeball in which we include an appropriate limbus ring that mimics the ciliary body and the iris. The Finite Element Method is used to model the eyeball and to test its behaviour. A set of geometrical and material parameters has been determined for the model so that the postulated function of the eye is preserved. Numerical simulations have confirmed the hypothesis. The anatomically justified inclusion of the limbus ring in the proposed model of the eyeball makes it more realistic than those previously reported.     

Correlation between the viscoelastic properties of a soft crystal and its microstructure

We investigate the rheological properties of a cubic fcc phase of micelles obtained by aggregation of a triblock copolymer (PEO)₁₂₇(PPO)₄₈(PEO)₁₂₇ in water as selective solvent. The resulting soft solid is submitted to a range of stresses varying from 20 to 800Pa in Couette geometry. Creep and flow behaviour can be distinguished and interpreted in terms of structural changes previously observed by SAXS under flow. Contrasting with other systems, no discontinuity in the flow behaviour is associated with the structural changes. The strong shear thinning is interpreted from the scattering data, as resulting from the nucleation of a new structure of hexagonal compact planes parallel to the Couette walls. This creates a lubricating domain in the gap, whose size grows with the applied shear rate. We argue moreover that the very existence of flow (as a steady state opposed to creep) is associated with this so-called layer-sliding structure in a fraction, however small, of the sample. Received on 4 June 1999 and Received in final form 6 September 1999     

Mechanical properties of silicon nanobeams with an undercut evaluated by combining the dynamic resonance test and finite element analysis

Mechanical properties of silicon nanobeams are of prime importance in nanoelectromechanical system applications.A numerical experimental method of determining resonant frequencies and Young’s modulus of nanobeams by combining finite element analysis and frequency response tests based on an electrostatic excitation and visual detection by using a laser Doppler vibrometer is presented in this paper.Silicon nanobeam test structures are fabricated from silicon-oninsulator wafers by using a standard lithography and anisotropic wet etching release process,which inevitably generates the undercut of the nanobeam clamping.In conjunction with three-dimensional finite element numerical simulations incorporating the geometric undercut,dynamic resonance tests reveal that the undercut significantly reduces resonant frequencies of nanobeams due to the fact that it effectively increases the nanobeam length by a correct value △L,which is a key parameter that is correlated with deviations in the resonant frequencies predicted from the ideal Euler-Bernoulli beam theory and experimentally measured data.By using a least-square fit expression including △L,we finally extract Young’s modulus from the measured resonance frequency versus effective length dependency and find that Young’s modulus of a silicon nanobeam with 200-nm thickness is close to that of bulk silicon.This result supports that the finite size effect due to the surface effect does not play a role in the mechanical elastic behaviour of silicon nanobeams with thickness larger than 200 nm.     

Mechanical properties an undercut evaluated resonance test and of silicon nanobeams with by combining the dynamic finite element analysis

Mechanical properties of silicon nanobeams are of prime importance in nanoelectromechanical system applications. A numerical experimental method of determining resonant frequencies and Young＇s modulus of nanobeams by combining finite element analysis and frequency response tests based on an electrostatic excitation and visual detection by using a laser Doppler vibrometer is presented in this paper. Silicon nanobeam test structures are fabricated from silicon-oninsulator wafers by using a standard lithography and anisotropic wet etching release process, which inevitably generates the undercut of the nanobeam clamping. In conjunction with three-dimensional finite element numerical simulations incorporating the geometric undercut, dynamic resonance tests reveal that the undercut significantly reduces resonant frequencies of nanobeams due to the fact that it effectively increases the nanobeam length by a correct value △L, which is a key parameter that is correlated with deviations in the resonant frequencies predicted from the ideal Euler-Bernoulli beam theory and experimentally measured data. By using a least-square fit expression including △L, we finally extract Young＇s modulus from the measured resonance frequency versus effective length dependency and find that Young＇s modulus of a silicon nanobeam with 200-nm thickness is close to that of bulk silicon. This result supports that the finite size effect due to the surface effect does not play a role in the mechanical elastic behaviour of silicon nanobeams with thickness larger than 200 nm.     

A novel method of rapidly modeling optical properties of actual photonic crystal fibres

<正>The flexible structure of photonic crystal fibre not only offers novel optical properties but also brings some difficulties in keeping the fibre structure in the fabrication process which inevitably cause the optical properties of the resulting fibre to deviate from the designed properties.Therefore,a method of evaluating the optical properties of the actual fibre is necessary for the purpose of application.Up to now,the methods employed to measure the properties of the actual photonic crystal fibre often require long fibre samples or complex expensive equipments.To our knowledge, there are few studies of modeling an actual photonic crystal fibre and evaluating its properties rapidly.In this paper,a novel method,based on the combination model of digital image processing and the finite element method,is proposed to rapidly model the optical properties of the actual photonic crystal fibre.Two kinds of photonic crystal fibres made by Crystal Fiber A/S are modeled.It is confirmed from numerical results that the proposed method is simple,rapid and accurate for evaluating the optical properties of the actual photonic crystal fibre without requiring complex equipment.     

金属半圆环/长板阵列的法诺共振特性

本文设计了金属半圆环/长板阵列,并应用有限元方法研究了该阵列的透射特性。研究表明:由于半圆环与长板之间的电场耦合,在该阵列中产生了法诺共振现象。法诺共振峰强烈地依赖于半圆环/长板的结构参数和相对位置,并且法诺共振峰对周围介质折射率有着较高的灵敏度,最高可以达到862.5 nm/RIU。这些结果有助于设计基于法诺共振的微纳光子学器件。     

The numerical prediction of planar viscoelastic contraction flows using the pom–pom model and higher-order finite volume schemes

This study investigates the numerical solution of viscoelastic flows using two contrasting high-order finite volume schemes. We extend our earlier work for Poiseuille flow in a planar channel and the single equation form of the extended pom–pom (SXPP) model [M. Aboubacar, J.P. Aguayo, P.M. Phillips, T.N. Phillips, H.R. Tamaddon-Jahromi, B.A. Snigerev, M.F. Webster, Modelling pom–pom type models with high-order finite volume schemes, J. Non-Newtonian Fluid Mech. 126 (2005) 207–220], to determine steady-state solutions for planar 4:1 sharp contraction flows. The numerical techniques employed are time-stepping algorithms: one of hybrid finite element/volume type, the other of pure finite volume form. The pure finite volume scheme is a staggered-grid cell-centred scheme based on area-weighting and a semi-Lagrangian formulation. This may be implemented on structured or unstructured rectangular grids, utilising backtracking along the solution characteristics in time. For the hybrid scheme, we solve the momentum-continuity equations by a fractional-staged Taylor–Galerkin pressure-correction procedure and invoke a cell-vertex finite volume scheme for the constitutive law. A comparison of the two finite volume approaches is presented, concentrating upon the new features posed by the pom–pom class of models in this context of non-smooth flows. Here, the dominant feature of larger shear and extension in the entry zone influences both stress and stretch, so that larger stretch develops around the re-entrant corner zone as Weissenberg number increases, whilst correspondingly stress levels decline.     

Thermal and mechanical finite element modeling of laser-generated ultrasound in coating–substrate system

Thermal and mechanical finite element modeling of laser-generated ultrasound in coating–substrate system is presented, which entails a numerical formulation for the transient responses in terms of the characteristics of the source that originated the waves in the thermoelastic regime. The formulation takes into account the temperature dependence of the material parameters and thermal diffusion from the source. Numerical results show that the temperature gradient field is equivalent to the body force source, namely, vertical forces propagate downward and horizontal forces outward in the specimen. Waveforms at the epicenter present fast oscillations immediately after the precursor, which are produced by multiple reflections of the longitudinal in the coating, the thickness of the coating has remarkable influence on the waveforms at the epicenter.     

Scattering properties of an individual metallic nano-spheroid by the incident polarized light wave

The scattering properties of a metallic nano-spheroid under the illumination of different polarized light waves are investigated using 3D boundary element method. The influences of different geometrical sizes of the nano-spheroid and incident directions of the illuminating light wave on the scattering spectrum are studied for different incident polarized light waves. The results show that the metallic nano-spheroid has two intrinsic resonant modes, corresponding to different polarization states and resonant wavelengths. The scattering enhancement, the resonant wavelength, and the location of the enhanced optical field are strongly dependent on the polarization properties of the illuminating light waves, and they can be modulated by appropriately choosing the polarization directions of the incident light wave.     

Investigation on periodic cracking of elastic film/substrate system by the extended finite element method

Multiple cracking behavior in a thin elastic film bonded to a thick elastic substrate is investigated by the extended finite element method. Stress and stress intensity factor are obtained using a periodic finite element model for the cracked film/substrate system. The influences of various parameters including crack length, film thickness, periodic crack spacing, and relative stiffness of the substrate on the stress and stress intensity factor are discussed in detail. It is demonstrated that the effects of geometric parameters are more sensitive than that of material property. In particular, the crack spacing has a saturation value due to interactions of neighboring cracks and relief of tensile stress in the film. The film/substrate couple with multiple periodic cracks can exhibit a positive potential in improving the durability of the film/substrate system.     

Simulation of the interface characterization of thin film on substrate system by bending creep tests

In the present paper, the finite element simulation of the bending creep tests of the thin film on substrate system is carried out. The purpose of the investigation is to understand the creep stress characterization of the thin film on substrate system with the three points bending creep test method, which plays an important role in the bending creep testing characterization, so as to provide some foundation on determination of interface properties of the thin film on substrate system by a bending creep testing. Finite element results shows that the influences of the thickness of thin film and the modulus ratio of thin film to substrate on stress distribution are important.     

空间太阳望远镜主镜支撑结构的优化设计

空间太阳望远镜主镜是有效口径为1m的抛物面镜,工作状态需要达到衍射极限,因此光学系统要求主镜面形误差小于λ/40(RMS),精度主要靠主镜支撑结构来保证。主镜支撑结构应满足地面调试、在轨及发射状态的需要。支撑结构试验样机已经加工完成,地面调试结果表明主镜的镜面变形满足整个光学系统的要求。试验样机强度和刚度还有较大余量,结构本身比较复杂。用有限元分析方法进行优化设计,优化后的主镜支撑结构满足地面调试、在轨及发射状态的需要,也能保证主镜的面形满足整个光学系统的要求,有效减轻仪器重量、简化支撑结构的同时,提高了整个仪器的可靠性。     

Reconstructing 3-D maps of the local viscoelastic properties using a finite-amplitude modulated radiation force

A modulated acoustic radiation force, produced by two confocal tone-burst ultrasound beams of slightly different frequencies (i.e. 2.0 MHz ± Δf/2, where Δf is the difference frequency), can be used to remotely generate modulated low-frequency (Δf ? 500 Hz) shear waves in attenuating media. By appropriately selecting the duration of the two beams, the energy of the generated shear waves can be concentrated around the difference frequency (i.e., Δf ± Δf/2). In this manner, neither their amplitude nor their phase information is distorted by frequency-dependent effects, thereby, enabling a more accurate reconstruction of the viscoelastic properties. Assuming a Voigt viscoelastic model, this paper describes the use of a finite-element-method model to simulate three-dimensional (3-D) shear-wave propagation in viscoelastic media containing a spherical inclusion. Nonlinear propagation is assumed for the two ultrasound beams, so that higher harmonics are developed in the force and shear spectrum. Finally, an inverse reconstruction algorithm is used to extract 3-D maps of the local shear modulus and viscosity from the simulated shear-displacement fields based on the fundamental and second-harmonic component. The quality of the reconstructed maps is evaluated using the contrast between the inclusion and the background and the contrast-to-noise ratio (CNR). It is shown that the shear modulus can be accurately reconstructed based on the fundamental component, such that the observed contrast deviates from the true contrast by a root-mean-square-error (RMSE) of only 0.38 and the CNR is greater than 30 dB. If the second-harmonic component is used, the RMSE becomes 1.54 and the corresponding CNR decreases by approximately 10–15 dB. The reconstructed shear viscosity maps based on the second harmonic are shown to be of higher quality than those based on the fundamental. The effects of noise are also investigated and a fusion operation between the two spectral components is applied to enhance the reconstruction quality. Finally, a modified shear-wave spectroscopy technique, shown to be more robust to noise, is described for the estimation of the viscoelastic properties inside and outside the spherical inclusion under conditions of increased noise.     

A numerical study of the target system of an ADSS with different flow guides

The mechanical design of the target module of an accelerator driven subcritical nuclear reactor system (ADSS) calls for an analysis of the related thermal-hydraulic issues because of large amount of heat deposition in the spallation region during the course of nuclear interactions with the molten lead bismuth eutectic (LBE) target. The LBE also should carry the entire heat generated as a consequence of the spallation reaction. The problem of heat removal by the LBE is a challenging thermal-hydraulic issue. For this, one has to examine the flows of low Prandtl number fluids (LBE) in a complex ADSS geometry. In this study, the equations governing the laminar flow and thermal energy are solved numerically using the streamline upwind Petrov-Galerkin (SUPG) finite element (FE) method. The target systems with a straight and a nozzle guide have been considered. The principal purpose of the analysis is to trace the flow and temperature distribution and thereby to check the suitability of the flow guide in avoiding the recirculation or stagnation zones in the flow space that may lead to hot spots.      

空间遥感器中窗口的热光学特性研究

光学窗口是空间光学遥感器应用系统中的重要组成部分 ,其热光学特性是其热控设计的指南。拟定某种温度边界条件后利用有限元法计算其稳态温度场分布 ,再将温度载荷和其它有关广义力载荷作为载荷 ,计算光学玻璃的变形。将变形值拟合为 Zemike多项式 ,输入到光学设计程序 Zmax TM中 ,计算出光学玻璃的变形所引入的波差。最终根据光学遥感器许用的系统波像差来给定光学玻璃的许用温度梯度     

Hyperspherical harmonics approach to the trinucleon system with Reid soft core potential: Calculation of geometrical structure coefficients

We present the full set of equations for the solution of the trinucleon problem by the hyperspherical harmonics expansion (HHE) method where nucleons interact via the Reid soft core (RSC) potential. The coupling potential matrix elements are expressed in terms of geometrical structure coefficients (GSC) and potential multipoles (PM). Introduction of GSC greatly simplifies the calculation of the potential matrix and makes the numerical algorithm efficient. A method for calculating all the twelve independent sets of GSC needed, by using the completeness property of the Jacobi polynomials has been presented. A convenient sum rule for each set of GSC has also been derived and precision of the calculated GSC has been checked by the sum rule. Such calculations of GSC are efficient and fast, in view of the complexity of the HHE equations.     

Numerical simulation of laser-generated surface acoustic waves in the transparent coating on a substrate by the finite element method

The optimum finite element model in the system consisting of a transparent coating and an opaque substrate is established based on the analysis of two important parameters: meshing size and time step, and the stability of solution. Taking into account the temperature dependence of material properties, the transient temperature and temperature gradient field are obtained. According to the thermoelastic theory, this temperature gradient field can be taken as a buried bulk source to generate ultrasonic wave. The surface acoustic waves (SAWs) are obtained. The influence of the coating thickness on the SAWs is analyzed. The model provides a useful tool for the determination of modes which are generated by a laser source in transparent coating on opaque substrate. The surface skimming longitudinal wave exists for the multiple oscillations and it charges from unipolar waveforms to dipolar.     

Effects of non-Gaussian noise on the dynamical properties of a logistic system

The dynamical properties of a tumor cell growth system described by the logistic system with coupling between non- Gaussian and Gaussian noise terms are investigated. The effects of the nonextensive index q on the stationary properties and the transient properties are discussed, respectively. The results show that the nonextensive index q can induce the tumor cell numbers to decrease greatly in the case of q 〉 1. Moreover, the switch from the steady stable state to the extinct state is speeded up as the increases of q, and the tumor cell numbers can be more obviously restrained for a large value of q. The numerical results are found to be in basic agreement with the theoretical predictions.