Normal shock boundary layer control with various vortex generator geometries

Various vortex generators which include ramp, split-ramp and a new hybrid concept "ramped-vane" are investigated under normal shock conditions with a diffuser at Mach number of 1.3. The dimensions of the computational domain were designed using Reynolds Average Navier-Stokes studies to be representative of the flow in an external-compression supersonic inlet. Using this flow geometry, various vortex generator concepts were studied with Implicit Large Eddy Simulation. In general, the ramped-vane provided increased vorticity compared to the other devices and reduced the separation length downstream of the device centerline. In addition, the size, edge gap and streamwise position respect to the shock were studied for the ramped-vane and it was found that a height of about half the boundary thickness and a large trailing edge gap yielded a fully attached flow downstream of the device. This ramped-vane also provided the largest reduction in the turbulent kinetic energy and pressure fluctuations. Additional benefits include negligible drag while the reductions in boundary layer displacement thickness and shape factor were seen compared to other devices. © 2011 Elsevier Ltd.

Multiphase Flow Measurement in the Offshore Oil Industry

In recent years, considerable research has been conducted into the development of a three-phase flowmeter suitable for use in an offshore environment, and oil/gas/water three-phase metering becomes an important aspect in multiphase flow measurement. This paper discusses the importance of three-phase flow measurement in offshore oil industry, describes the current development in this area, and points out the principal strategies which may be used to meter three-phase flow.

I. Rigid body penetration into brittle materials. II. Phase change effect on shock wave propagation

Part I.

We have developed a technique for measuring the depth time history of rigid body penetration into brittle materials (hard rocks and concretes) under a deceleration of ~ 10⁵ g. The technique includes bar-coded projectile, sabot-projectile separation, detection and recording systems. Because the technique can give very dense data on penetration depth time history, penetration velocity can be deduced. Error analysis shows that the technique has a small intrinsic error of ~ 3-4 % in time during penetration, and 0.3 to 0.7 mm in penetration depth. A series of 4140 steel projectile penetration into G-mixture mortar targets have been conducted using the Caltech 40 mm gas/ powder gun in the velocity range of 100 to 500 m/s.

We report, for the first time, the whole depth-time history of rigid body penetration into brittle materials (the G-mixture mortar) under 10⁵ g deceleration. Based on the experimental results, including penetration depth time history, damage of recovered target and projectile materials and theoretical analysis, we find:

1. Target materials are damaged via compacting in the region in front of a projectile and via brittle radial and lateral crack propagation in the region surrounding the penetration path. The results suggest that expected cracks in front of penetrators may be stopped by a comminuted region that is induced by wave propagation. Aggregate erosion on the projectile lateral surface is < 20% of the final penetration depth. This result suggests that the effect of lateral friction on the penetration process can be ignored.

2. Final penetration depth, P_max, is linearly scaled with initial projectile energy per unit cross-section area, e_s , when targets are intact after impact. Based on the experimental data on the mortar targets, the relation is P_max(mm) 1.15e_s (J/mm² ) + 16.39.

3. Estimation of the energy needed to create an unit penetration volume suggests that the average pressure acting on the target material during penetration is ~ 10 to 20 times higher than the unconfined strength of target materials under quasi-static loading, and 3 to 4 times higher than the possible highest pressure due to friction and material strength and its rate dependence. In addition, the experimental data show that the interaction between cracks and the target free surface significantly affects the penetration process.

4. Based on the fact that the penetration duration, t_max, increases slowly with e_s and does not depend on projectile radius approximately, the dependence of t_max on projectile length is suggested to be described by t_max(μs) = 2.08e_s (J/mm² + 349.0 x m/(πR²), in which m is the projectile mass in grams and R is the projectile radius in mm. The prediction from this relation is in reasonable agreement with the experimental data for different projectile lengths.

5. Deduced penetration velocity time histories suggest that whole penetration history is divided into three stages: (1) An initial stage in which the projectile velocity change is small due to very small contact area between the projectile and target materials; (2) A steady penetration stage in which projectile velocity continues to decrease smoothly; (3) A penetration stop stage in which projectile deceleration jumps up when velocities are close to a critical value of ~ 35 m/s.

6. Deduced averaged deceleration, a, in the steady penetration stage for projectiles with same dimensions is found to be a(g) = 192.4v + 1.89 x 10⁴, where v is initial projectile velocity in m/s. The average pressure acting on target materials during penetration is estimated to be very comparable to shock wave pressure.

7. A similarity of penetration process is found to be described by a relation between normalized penetration depth, P/P_max, and normalized penetration time, t/t_max, as P/P_max = f(t/t_max, where f is a function of t/t_max. After f(t/t_max is determined using experimental data for projectiles with 150 mm length, the penetration depth time history for projectiles with 100 mm length predicted by this relation is in good agreement with experimental data. This similarity also predicts that average deceleration increases with decreasing projectile length, that is verified by the experimental data.

8. Based on the penetration process analysis and the present data, a first principle model for rigid body penetration is suggested. The model incorporates the models for contact area between projectile and target materials, friction coefficient, penetration stop criterion, and normal stress on the projectile surface. The most important assumptions used in the model are: (1) The penetration process can be treated as a series of impact events, therefore, pressure normal to projectile surface is estimated using the Hugoniot relation of target material; (2) The necessary condition for penetration is that the pressure acting on target materials is not lower than the Hugoniot elastic limit; (3) The friction force on projectile lateral surface can be ignored due to cavitation during penetration. All the parameters involved in the model are determined based on independent experimental data. The penetration depth time histories predicted from the model are in good agreement with the experimental data.

9. Based on planar impact and previous quasi-static experimental data, the strain rate dependence of the mortar compressive strength is described by σ_f/σ⁰_f = exp(0.0905(log(έ/έ_0) ^1.14, in the strain rate range of 10^-7/s to 10³/s (σ⁰_f and έ are reference compressive strength and strain rate, respectively). The non-dispersive Hugoniot elastic wave in the G-mixture has an amplitude of ~ 0.14 GPa and a velocity of ~ 4.3 km/s.

Part II.

Stress wave profiles in vitreous GeO₂ were measured using piezoresistance gauges in the pressure range of 5 to 18 GPa under planar plate and spherical projectile impact. Experimental data show that the response of vitreous GeO₂ to planar shock loading can be divided into three stages: (1) A ramp elastic precursor has peak amplitude of 4 GPa and peak particle velocity of 333 m/s. Wave velocity decreases from initial longitudinal elastic wave velocity of 3.5 km/s to 2.9 km/s at 4 GPa; (2) A ramp wave with amplitude of 2.11 GPa follows the precursor when peak loading pressure is 8.4 GPa. Wave velocity drops to the value below bulk wave velocity in this stage; (3) A shock wave achieving final shock state forms when peak pressure is > 6 GPa. The Hugoniot relation is D = 0.917 + 1.711u (km/s) using present data and the data of Jackson and Ahrens [1979] when shock wave pressure is between 6 and 40 GPa for ρ₀ = 3.655 gj cm³ . Based on the present data, the phase change from 4-fold to 6-fold coordination of Ge⁺⁴ with O^-2 in vitreous GeO₂ occurs in the pressure range of 4 to 15 ± 1 GPa under planar shock loading. Comparison of the shock loading data for fused SiO₂ to that on vitreous GeO₂ demonstrates that transformation to the rutile structure in both media are similar. The Hugoniots of vitreous GeO₂ and fused SiO₂ are found to coincide approximately if pressure in fused SiO₂ is scaled by the ratio of fused SiO₂to vitreous GeO₂ density. This result, as well as the same structure, provides the basis for considering vitreous Ge0₂ as an analogous material to fused SiO₂ under shock loading. Experimental results from the spherical projectile impact demonstrate: (1) The supported elastic shock in fused SiO₂ decays less rapidly than a linear elastic wave when elastic wave stress amplitude is higher than 4 GPa. The supported elastic shock in vitreous GeO₂ decays faster than a linear elastic wave; (2) In vitreous GeO₂ , unsupported shock waves decays with peak pressure in the phase transition range (4-15 GPa) with propagation distance, x, as α 1/x^-3.35 , close to the prediction of Chen et al. [1998]. Based on a simple analysis on spherical wave propagation, we find that the different decay rates of a spherical elastic wave in fused SiO₂ and vitreous GeO₂ is predictable on the base of the compressibility variation with stress under one-dimensional strain condition in the two materials.

A new ground motion intensity measure, peak filtered acceleration (PFA), to estimate collapse vulnerability of buildings in earthquakes

In this thesis, we develop an efficient collapse prediction model, the PFA (Peak Filtered Acceleration) model, for buildings subjected to different types of ground motions.

For the structural system, the PFA model covers modern steel and reinforced concrete moment-resisting frame buildings (potentially reinforced concrete shear wall buildings). For ground motions, the PFA model covers ramp-pulse-like ground motions, long-period ground motions, and short-period ground motions.

To predict whether a building will collapse in response to a given ground motion, we first extract long-period components from the ground motion using a Butterworth low-pass filter with suggested order and cutoff frequency. The order depends on the type of ground motion, and the cutoff frequency depends on the building’s natural frequency and ductility. We then compare the filtered acceleration time history with the capacity of the building. The capacity of the building is a constant for 2-dimentional buildings and a limit domain for 3-dimentional buildings. If the filtered acceleration exceeds the building’s capacity, the building is predicted to collapse. Otherwise, it is expected to survive the ground motion.

The parameters used in PFA model, which include fundamental period, global ductility and lateral capacity, can be obtained either from numerical analysis or interpolation based on the reference building system proposed in this thesis.

The PFA collapse prediction model greatly reduces computational complexity while archiving good accuracy. It is verified by FEM simulations of 13 frame building models and 150 ground motion records.

Based on the developed collapse prediction model, we propose to use PFA (Peak Filtered Acceleration) as a new ground motion intensity measure for collapse prediction. We compare PFA with traditional intensity measures PGA, PGV, PGD, and Sa in collapse prediction and find that PFA has the best performance among all the intensity measures.

We also provide a close form in term of a vector intensity measure (PGV, PGD) of the PFA collapse prediction model for practical collapse risk assessment.

Structural and clumped-isotope constraints on the mechanisms of displacement along low-angle detachments

Despite years of research on low-angle detachments, much about them remains enigmatic. This thesis addresses some of the uncertainty regarding two particular detachments, the Mormon Peak detachment in Nevada and the Heart Mountain detachment in Wyoming and Montana.

Constraints on the geometry and kinematics of emplacement of the Mormon Peak detachment are provided by detailed geologic mapping of the Meadow Valley Mountains, along with an analysis of structural data within the allochthon in the Mormon Mountains. Identifiable structures well suited to constrain the kinematics of the detachment include a newly mapped, Sevier-age monoclinal flexure in the hanging wall of the detachment. This flexure, including the syncline at its base and the anticline at its top, can be readily matched to the base and top of the frontal Sevier thrust ramp, which is exposed in the footwall of the detachment to the east in the Mormon Mountains and Tule Springs Hills. The ~12 km of offset of these structural markers precludes the radial sliding hypothesis for emplacement of the allochthon.

The role of fluids in the slip along faults is a widely investigated topic, but the use of carbonate clumped-isotope thermometry to investigate these fluids is new. Faults rocks from within ~1 m of the Mormon Peak detachment, including veins, breccias, gouges, and host rocks, were analyzed for carbon, oxygen, and clumped-isotope measurements. The data indicate that much of the carbonate breccia and gouge material along the detachment is comminuted host rock, as expected. Measurements in vein material indicate that the fluid system is dominated by meteoric water, whose temperature indicates circulation to substantial depths (c. 4 km) in the upper crust near the fault zone.

Slip along the subhorizontal Heart Mountain detachment is particularly enigmatic, and many different mechanisms for failure have been proposed, predominantly involving catastrophic failure. Textural evidence of multiple slip events is abundant, and include multiple brecciation events and cross-cutting clastic dikes. Footwall deformation is observed in numerous exposures of the detachment. Stylolitic surfaces and alteration textures within and around “banded grains” previously interpreted to be an indicator of high-temperature fluidization along the fault suggest their formation instead via low-temperature dissolution and alteration processes. There is abundant textural evidence of the significant role of fluids along the detachment via pressure solution. The process of pressure solution creep may be responsible for enabling multiple slip events on the low-angle detachment, via a local rotation of the stress field.

Clumped-isotope thermometry of fault rocks associated with the Heart Mountain detachment indicates that despite its location on the flanks of a volcano that was active during slip, the majority of carbonate along the Heart Mountain detachment does not record significant heating above ambient temperatures (c. 40-70°C). Instead, cold meteoric fluids infiltrated the detachment breccia, and carbonate precipitated under ambient temperatures controlled by structural depth. Locally, fault gouge does preserve hot temperatures (>200°C), as is observed in both the Mormon Peak detachment and Heart Mountain detachment areas. Samples with very hot temperatures attributable to frictional shear heating are present but rare. They appear to be best preserved in hanging wall structures related to the detachment, rather than along the main detachment.

Evidence is presented for the prevalence of relatively cold, meteoric fluids along both shallow crustal detachments studied, and for protracted histories of slip along both detachments. Frictional heating is evident from both areas, but is a minor component of the preserved fault rock record. Pressure solution is evident, and might play a role in initiating slip on the Heart Mountain fault, and possibly other low-angle detachments.

Avaliação de Correlações e Equações de Estado para Determinação de Fatores de Compressibilidade de Gás Natural

O fator de compressibilidade (Z) de gás natural é utilizado em vários cálculos na engenharia de petróleo (avaliação de formações, perda de carga em tubulações, gradiente de pressão em poços de gás, cálculos de balanço de massa, medição de gás, compressão e processamento de gás). As fontes mais comuns de valores de Z são medições experimentais, caras e demoradas. Essa propriedade também é estimada por correlações empíricas, modelos baseados no princípio dos estados correspondentes ou equações de estado (EOS). Foram avaliadas as capacidades das EOS de Soave-Redlich-Kwong (SRK), Peng-Robinson (PR), Patel-Teja (PT), Patel-Teja-Valderrama (PTV), Schmidt-Wenzel (SW), Lawal-Lake-Silberberg (LLS) e AGA-8 para previsão desta propriedade em aproximadamente 2200 pontos de dados experimentais. Estes pontos foram divididos em quatro grupos: Grupo 1 (Presença de frações C7+, Grupo 2 (temperaturas inferiores a 258,15 K), Grupo 3 (pressões superiores a 10000 kPa) e Grupo 4 (pressões inferiores a 10000 kPa). Os cálculos utilizando as equações de estado sob diferentes esquemas de previsão de coeficientes binários de interação foram cuidadosamente investigados. Os resultados sugerem que a EOS AGA-8 apresenta os menores erros para pressões de até 70000 kPa. Entretanto, observou-se uma tendência de aumento nos desvios médios absolutos em função das concentrações de CO2 e H2S. As EOS PTV e a EOS SW são capazes de predizer o fator de compressibilidade (Z) com desvios médios absolutos entre os valores calculados e experimentais com precisão satisfatória para a maioria das aplicações, para uma variada faixa de temperatura e pressão. Este estudo também apresenta uma avaliação de 224 métodos de cálculo de Z onde foram utilizadas 8 correlações combinadas com 4 regras de mistura para estimativa de temperaturas e pressões pseudorreduzidas das amostras, junto com 7 métodos de caracterização das propriedades críticas da fração C7+, quando presente na composição do gás. Em função dos resultados são sugeridas, para diferentes tipos de sistemas, as melhores combinações de correlações com regras de mistura capazes de predizer fatores de compressibilidade (Z) com os menores erros absolutos médios relativos

Reflex impairment as a measure of vitality and survival potential of Atlantic cod (Gadus morhua)

A ssur ing the v itality and survival potential of live-caught Atlantic cod (Gadus morhua) is important for improving the sorting of fish before net penning operations designed to hold fish for growth and later market. When Atlantic cod are captured by Danish seine, the most commonly used fishing gear for live-caught fish, they undergo stressors such as forced swimming, net abrasion, and air exposure. Laboratory experiments (at an air temperature of 9°C and water temperature of 8°C) were conducted with the aim of constructing a RAMP (reflex action mortality predictor) curve for prediction of vitality and survival potential in Atlantic cod captured in Danish seines, by varying the levels of these stressors. Atlantic cod exposed to increased duration in air (5–20 min) showed increased reflex impairment and mortality, with 75% mortality at 10 minutes of air exposure. Forced swimming in combination with net abrasion and air exposure did not increase reflex impairment or mortality above that associated with air exposure alone. The Atlantic cod RAMP curves indicated that fish with reflex impairment less than 50% would not show mortality and would likely recover from capture stress.

Identificação da formação de titanato de bário a partir de mistura reacional calcinada em diferentes temperaturas via difração de raios X, espectroscopia fotoacústica e análise térmica

Os materiais ferroelétricos têm sido utilizados em muitas áreas da tecnologia e da ciência, pois possuem um grande número de aplicações, como: sensores; transdutores; capacitores; dispositivos ópticos; dentre outras. A busca por novos materiais cerâmicos ferroelétricos tem sido grande. Um dos materiais cerâmicos ferroelétricos mais estudados é o titanato de bário (BT). São vários os métodos de produção e caracterização do titanato de bário. Neste trabalho, pós cerâmicos de titanato de bário foram obtidos por reação do estado sólido a partir de misturas reacionais calcinadas em diferentes temperaturas entre 400C e 900C. Foram três as misturas reacionais: não dopadas; dopadas com 1%; e dopadas com 5% de dióxido de cério (CeO2). A identificação da formação do BT, nos pós cerâmicos produzidos, foi feita a partir de três técnicas de caracterização: difração de raios X (DRX); espectroscopia fotoacústica (PAS); e técnicas de análise térmica. Com a técnica DRX, difratogramas mostraram que a plena formação do titanato de bário ocorreu a partir da temperatura de calcinação de 700C. Para a amostra não dopada com cério e calcinada a 800C, houve deslocamento de todos os picos de difração. Nas amostras dopadas com dióxido de cério houve deslocamento de todos os picos de difração, em relação as amostras não dopadas. Observou-se também que nas amostras dopadas com 5% de CeO2, e calcinadas a 700C e 800C, resíduos de dióxido de cério foram observados nos difratogramas. Com a técnica PAS, espectros de absorção foram obtidos. Foi possível observar uma grande diferença de absorção da amostra calcinada a 600 e 630C, indicando a formação do titanato de bário a partir da temperatura de 630C, nas amostras sem a dopagem dióxido de cério. Houve um alargamento nas bandas de absorção a partir da temperatura de 600C, quando o dióxido de cério entrou na matriz. Foi também possível determinar as energias de band-gap das amostras utilizando o método de Tauc. Com as técnicas de análise térmica, em especial através da técnica termogravimétrica (TG/DTG), foi comprovado que até 400C não havia formação de titanato de bário. Visto que nesta temperatura de calcinação houve a maior perda de massa durante a rampa de aquecimento. O início da formação do titanato de bário foi observado a partir da temperatura de calcinação de 500C, assim como nas técnicas DRX e PAS. Portanto, com os resultados apresentados, foi demonstrada a identificação da formação do titanato de bário nas misturas reacionais calcinadas, com auxílio das potencialidades das três técnicas utilizadas.

Efeito da estimulação transcraniana por corrente contínua sobre o desempenho físico: análise das respostas cerebrais, metabólicas e perceptivas durante e após exercícios aeróbio e de força

Estimulação transcraniana por corrente contínua (ETCC) sobre áreas corticais pré-selecionadas, tem aumentado o desempenho físico de diferentes populações. Porém, lacunas persistem no tocante aos mecanismos subjacentes à estes efeitos. Assim, a presente tese objetivou: a) investigar os efeitos da ETCC anódica (aETCC) e placebo (Sham) no córtex motor (CM) de indivíduos saudáveis sobre o desempenho de força máxima; b) comparar os efeitos da ETCC sobre a produção de força máxima e estabilidadade da força durante exercícios máximo e submáximo em sujeitos hemiparéticos e saudáveis; c) investigar o efeito da ETCC sobre a conectividade funcional inter-hemisférica (coerência eletroencefalográfica cEEG) do córtex pré-frontal (CPF), desempenho aeróbio e dispêndio energético (EE) durante e após exercício máximo e submáximo. No 1 estudo, 14 adultos saudáveis executaram 2 sessões de exercício máximo de força (EMF) dos músculos flexores e extensores do joelho dominante (3 séries de 10 rep máximas), precedidos por aETCC ou Sham (2mA; 20 mim). aETCC não foi capaz de aumentar o trabalho total e pico de torque (PT), resistência à fadiga ou atividade eletromiográfica durante o EMF. No 2 estudo, 10 hemiparéticos e 9 sujeitos saudáveis receberam aETCC e Sham no CM. O PT e a estabilidade da força (coeficiente de variação - CV) foram avaliados durante protocolo máximo e submáximo de extensão e flexão unilateral do joelho (1 série de 3 reps a 100% do PT e 2 séries de 10 reps a 50% do PT). Nenhuma diferença no PT foi observada nos dois grupos. Diminuições no CV foram obervadas durante a extensão (~25-35%, P<0.001) e flexão de joelho (~22-33%, P<0.001) após a aETCC comparada com Sham nos hemiparéticos, entretanto, somente o CV na extensão de joelhos diminuiu (~13-27%, P<0.001) nos saudáveis, o que sugere que aETCC pode melhorar o CV, mas não o PT em sujeitos hemiparéticos. No 3 estudo, 9 adultos saudáveis realizaram 2 testes incrementais máximos precedidos por aETCC ou Sham sobre o CPF com as respostas cardiorrespiratórias, percepção de esforço (PSE) e cEEG do CPF sendo monitoradas. O VO2 de pico (42.64.2 vs. 38.23.3 mL.kg.min-1; P=0,02), potência total (252.776.5 vs. 23773.3 W; P=0,05) e tempo de exaustão (531.1140 vs. 486.7115.3 seg; P=0,04) foram maiores após aETCC do que a Sham. Nenhuma diferença foi encontrada para FC e PSE em função da carga de trabalho (P>0,05). A cEEG do CPF aumentou após aETCC vs. repouso (0.700.40 vs. 0.380.05; P=0,001), mas não após Sham vs. repouso (0.360.49 vs. 0.330.50; P=0,06), sugerindo que a aETCC pode retardar a fadiga aumentando a conectividade funcional entre os hemisférios do CPF e desempenho aeróbio durante exercício exaustivo. No 4 estudo, o VO2 e EE foram avaliados em 11 adultos saudáveis antes, durante a aETCC ou Sham no CPF e 30 min após exercício aeróbio submáximo isocalórico (~200kcal). Diferenças não foram observadas no VO2 vs. repouso durante aETCC e Sham (P=0.95 e P=0.85). Porém, a associação entre exercício e aETCC aumentou em ~19% o EE após ao menos, 30 min de recuperação após exercício quando comparada a Sham (P<0,05).

Impact of high-k dielectrics on breakdown performances of SiC and diamond Schottky diodes

This paper presents a comparison between SiC and diamond Schottky barrier diodes using the oxide ramp termination. The influences of the dielectric thickness and relative permittivity on the diode's electrical performance are investigated. Typical commercial drift layer parameters are used for this study. The extension of the space charge area throughout the drift region and the current distribution at breakdown are shown. The efficiency of the termination is also evaluated for both SiC and diamond diodes. © (2009) Trans Tech Publications, Switzerland.

On-line topographic analysis in the SEM

A novel method for on-line topographic analysis of rough surfaces in the SEM has been investigated. It utilises a digital minicomputer configured to act as a programmable scan generator and automatic focusing unit. The computer is coupled to the microscope through digital-to-analogue converters which enable it to generate ramp waveforms allowing the beam to be scanned over a small sub-region of the field under program control. A further digital-to-analogue converter regulates the current supply to the objective lens of the microscope. The video signal is sampled by means of an analogue-to-digital converter and the resultant binary code stored in the computer's memory as an array of numbers describing relative image intensity. Computations based on the intensity gradient of the image allow the objective lens current to be found for the in-focus condition, which may be related to the working distance through a previous calibration experiment. The sensitivity of the method for detecting small height changes is theoretically of the order of 1 μm. In practice the operator specifies features of interest by means of a mobile spot cursor injected into the SEM display screen, or he may scan the specimen at sub-regions corresponding to pre-determined points on a regular grid defined by him. The operation then proceeds under program control. | A novel method for on-line topographic analysis of rough surfaces in the SEM has been investigated. It utilizes a digital minicomputer configured to act as a programmable scan generator and automatic focusing unit. A further digital-to-analog converter regulates the current supply to the objective lens of the microscope. The video signal is sampled by means of an analog-to-digital converter and the resultant binary code stored in the computer's memory as an array of numbers describing relative image intensity. The sensitivity of the method for detecting small height changes is theroretically of the order of 1 mu m.

Topographical analysis in the SEM using an automatic focusing technique

A novel technique for automated topographical analysis in the SEM has been investigated. It utilizes a 16-bit minicomputer arranged to act as an automatic focusing unit. The computer is coupled to the objective lens of the microscope, by means of a digital to analogue converter, and may regulate the excitation of the lens under program control. Further digital-to-analogue converters allow the computer to act as a programmable scan generator by applying ramp waveforms to the scan amplifiers, permitting the beam to be swept over a small sub-region of the field of interest. The video signal is sampled and applied to an analogue-to-digital converter; the resultant binary numbers are stored in computer memory as an array of values representing relative image intensities within a subregion. A differencing algorithm applied to the collected data allows the level of objective lens excitation to be found at which the sharpness of the image is optimized, and the excitation may be related to the working distance for that subregion through a previous calibration experiment. The sensitivity of the method for detecting small height changes is theoretically of the order of 1 μm.

Analytical techniques for indentation of viscoelastic materials

Indentation of linearly viscoelastic materials is explored using elastic-viscoelastic correspondence analysis for both conical-pyramidal and spherical indentation. Boltzmann hereditary integrals are used to generate displacement-time solutions for loading at constant rate and creep following ramp loading. Experimental data for triangle- and trapezoidal-loading are examined for commercially-available polymers and compared with analytical solutions. Emphasis is given to the use of multiple experiments to test the fidelity and predictive capability of the obtained material creep function. Plastic deformation occurs in sharp indentation of glassy polymers and is found to complicate the viscoelastic analysis. A new method is proposed for estimating a material time-constant from peak displacement or hardness data obtained in pyramidal indentation tests performed at different loading rates.

Nanoindentation measurements of bone viscoelasticity as a function of hydration state

Bone is an anisotropic material, and its mechanical properties are determined by its microstructure as well as its composition. Mechanical properties of bone are a consequence of the proportions of, and the interactions between, mineral, collagen and water. Water plays an important role in maintaining the mechanical integrity of the composite, but the manner in which water interacts within the ultrastructure is unclear. Dentine being an isotropic two-dimensional structure presents a homogenous composite to examine the dehydration effects. Nanoindentation methods for determining the viscoelastic properties have recently been developed and are a subject of great interest. Here, one method based on elastic-viscoelastic correspondence for 'ramp and hold' creep testing (Oyen, J. Mater. Res., 2005) has been used to analyze viscoelastic behavior of polymeric and biological materials. The method of 'ramp and hold' allows the shear modulus at time zero to be determined from fitting of the displacement during the maximum load hold. Changes in the viscoelastic properties of bone and dentine were examined as the material was systematically dehydrated in a series of water:solvent mixes. Samples of equine dentine were sectioned and cryo-polished. Shear modulus was obtained by nanoindentation using spherical indenters with a maximum load hold of 120s. Samples were tested in different solvent concentrations sequentially, 70% ethanol to 50% ethanol, 70 % ethanol to 100% ethanol, 70% ethanol to 70% methanol to 100% methanol, and 70% ethanol to 100% acetone, after storage in each condition for 24h. By selectively removing and then replacing water from the composite, insights in to the ultrastructure of the tissue can be gained from the corresponding changes in the experimentally determined moduli, as well as an understanding of the complete reversibility of the dehydration process. © 2006 Materials Research Society.

Estimating material elasticity by spherical indentation load-relaxation tests on viscoelastic samples of finite thickness

A two-step viscoelastic spherical indentation method is proposed to compensate for 1) material relaxation and 2) sample thickness. In the first step, the indenter is moved at a constant speed and the reaction force is measured. In the second step, the indenter is held at a constant position and the relaxation response of the material is measured. Then the relaxation response is fit with a multi-exponential function which corresponds to a three-branch general Maxwell model. The relaxation modulus is derived by correcting the finite ramp time introduced in the first step. The proposed model takes into account the sample thickness, which is important for applications in which the sample thickness is less than ten times the indenter radius. The model is validated numerically by finite element simulations. Experiments are carried out on a 10% gelatin phantom and a chicken breast sample with the proposed method. The results for both the gelatin phantom and the chicken breast sample agree with the results obtained from a surface wave method. Both the finite element simulations and experimental results show improved elasticity estimations by incorporating the sample thickness into the model. The measured shear elasticities of the 10% gelatin sample are 6.79 and 6.93 kPa by the proposed finite indentation method at sample thickness of 40 and 20 mm, respectively. The elasticity of the same sample is estimated to be 6.53 kPa by the surface wave method. For the chicken breast sample, the shear elasticity is measured to be 4.51 and 5.17 kPa by the proposed indentation method at sample thickness of 40 and 20 mm, respectively. Its elasticity is measured by the surface wave method to be 4.14 kPa. © 2011 IEEE.