线性压电材料的有限元分析

本文利用线性压电学理论，编制了线性压电材料四结点等参有限元程序，进行了校核，并对PZT-5A材料压电智能元件和压电材料标准断裂试进行了计算，计算包括：①采用实际工程应用的压电智能元件尺寸，计算了元件的压电响应；并针对元件内部电极尖端区域容易引起破坏的现象，计算了该区域的奇异应力、应变场及电场。②计算了加力和加电两种情况下压电材料标准断裂试件应力强度因子影响系数F_I和电位移强度因子影响系数F_D。裂纹面边界条件采用D-P条件，试件包括紧凑拉伸标准试件和三点弯曲标准试件。③应用Lagrange乘子法将Parton裂纹边界条件加于有限元程序中，计算了上述两种标准断裂试件的F_I和F_D，计算结果与采用D-P裂纹边界条件的计算结果有很大差异。

Numerical Study On Transient Flow In The Deep Naturally Fractured Reservoir With High Pressure

According to the experimental results and the characteristics of the pressure-sensitive fractured formation, a transient flow model is developed for the deep naturally-fractured reservoirs with different outer boundary conditions. The finite element equations for the model are derived. After generating the unstructured grids in the solution regions, the finite element method is used to calculate the pressure type curves for the pressure-sensitive fractured reservoir with different outer boundaries, such as the infinite boundary, circle boundary and combined linear boundaries, and the characteristics of the type curves are comparatively analyzed. The effects on the pressure curves caused by pressure sensitivity module and the effective radius combined parameter are determined, and the method for calculating the pressure-sensitive reservoir parameters is introduced. By analyzing the real field case in the high temperature and pressure reservoir, the perfect results show that the transient flow model for the pressure-sensitive fractured reservoir in this paper is correct.

Effects of rainfall infiltration on deep slope failure

With the finite element method and the limit equilibrium method, a numerical model has been estab-lished for examining the effects of rainfall infiltration on the stability of slopes. This model is able to reflect the variations in pore water pressure field in slopes, dead weight of the soil, and soil softening caused by rainfall infiltration. As a case study, an actual landslide located at the Nongji Jixiao in Chongqing was studied to analyze the effects of rainfall infiltration on the seepage field and slope sta-bility. The simulated results showed that a deep slope failure is prone to occur when rainfall infiltration leads to a remarkable variation in the seepage field, especially when the pore water pressure in slopes increases in a large range.

扩展有限元方法计算多夹杂问题时圆形夹杂与四边形单元的几何关系

用扩展有限元法XFEM(Extended Finite Element Method)解决夹杂问题时,夹杂与基质的界面把单元分成若干部分.求单元刚度矩阵时,需要分别在这各个部分求积分.找到便于程序编制的描述各积分区域几何形状的方法是亟待解决的问题.本文把各积分区域的形成过程看成是圆对四边形的多次切割.考虑切剩区域与圆的关系时,把不完整的边仍看作完整的边,把切剩区域看成是四边形或是切去一两条边的四边形.采用排列组合的方法,把它们与圆的所有位置关系列了出来.

On The Thermally Induced Cracking Of A Segmented Coating Deposited On The Outer Surface Of A Hollow Cylinder

In this work, the thermally induced cracking behavior of a segmented coating has been investigated. The geometry under consideration is a hollow cylinder with a segmented coating deposited onto its outer surface. The segmentation cracks are modeled as a periodic array of axial edge cracks. The finite element method is utilized to obtain the solution of the multiple crack problem and the Thermal Stress Intensity Factors (TSIFs) are calculated. Based on dimensional analysis, the main parameters affecting TSIFs are identified. It has been found that the TSIF is a monotonically increasing function of segmentation crack spacing. This result confirms that a segmented coating exhibits much higher thermal shock resistance than an intact counterpart, if only the segmentation crack spacing is narrow enough. The dependence of TSIF on some other parameters, such as normalized time, segmentation crack depth, convection severity as well as material constants, has also been discussed. (C) 2008 Elsevier B.V. All rights reserved.

Random response of nonlinear continuous systems

This thesis presents a technique for obtaining the stochastic response of a nonlinear continuous system. First, the general method of nonstationary continuous equivalent linearization is developed. This technique allows replacement of the original nonlinear system with a time-varying linear continuous system. Next, a numerical implementation is described which allows solution of complex problems on a digital computer. In this procedure, the linear replacement system is discretized by the finite element method. Application of this method to systems satisfying the one-dimensional wave equation with two different types of constitutive nonlinearities is described. Results are discussed for nonlinear stress-strain laws of both hardening and softening types.

EVALUATION OF FILM-SUBSTRATE ADHESION VIA IMPACT USING COATED BULLETS

A method was devised to evaluate the adhesion between a film and a substrate. A front-end coated bullet is accelerated by a gas gun and hits the substrate of the specimen under test. The impact generates a compressive stress pulse that propagates toward the film. After transmission through the interface, part of the pulse is reflected on the free surface of the film, and tensile stress arises at the film-substrate interface, possibly inducing debonding of the film. This dynamic process was demonstrated analytically and simulated numerically by the finite element method. The results validate the initial concept and lay the foundation for further optimization of this method.

飞秒激光在石英玻璃中诱导微爆炸的理论研究

利用有限元方法建立了二维模型，研究了飞秒激光作用下石英玻璃中导带电子的产生、激光能量的沉积、导带电子和能量扩散等微观过程．计算了导带电子扩散引起的局部净电荷及其形成的静电场分布，初步揭示了微爆炸的演化过程．

Thin metastructures with engineered thermal expansion

The geometry and constituent materials of metastructures can be used to engineer the thermal expansion coefficient. In this thesis, we design, fabricate, and test thin thermally stable metastructures consisting of bi-metallic unit cells and show how the coefficient of thermal expansion (CTE) of these metastructures can be finely and coarsely tuned by varying the CTE of the constituent materials and the unit cell geometry. Planar and three-dimensional finite element method modeling is used to drive the design and inform experiments, and predict the response of these metastructures. We demonstrate computationally the significance of out-of-plane effects in the metastructure response. We develop an experimental setup using digital image correlation and an infrared camera to experimentally measure full displacement and temperature fields during testing and accurately measure the metastructures’ CTE. We experimentally demonstrate high aspect ratio metastructures of Ti/Al and Kovar/Al which exhibit near-zero and negative CTE, respectively. We demonstrate robust fabrication procedures for thermally stable samples with high aspect ratios in thin foil and thin film scales. We investigate the lattice structure and mechanical properties of thin films comprising a near-zero CTE metastructure. The mechanics developed in this work can be used to engineer metastructures of arbitrary CTE and can be extended to three dimensions.

Dynamic response of structures with uncertain parameters

This thesis presents a technique for obtaining the response of linear structural systems with parameter uncertainties subjected to either deterministic or random excitation. The parameter uncertainties are modeled as random variables or random fields, and are assumed to be time-independent. The new method is an extension of the deterministic finite element method to the space of random functions.

First, the general formulation of the method is developed, in the case where the excitation is deterministic in time. Next, the application of this formulation to systems satisfying the one-dimensional wave equation with uncertainty in their physical properties is described. A particular physical conceptualization of this equation is chosen for study, and some engineering applications are discussed in both an earthquake ground motion and a structural context.

Finally, the formulation of the new method is extended to include cases where the excitation is random in time. Application of this formulation to the random response of a primary-secondary system is described. It is found that parameter uncertainties can have a strong effect on the system response characteristics.

Análise de fadiga em pontes rodoviárias de concreto armado

As pontes rodoviárias de concreto armado estão sujeitas às ações dinâmicas variáveis devido ao tráfego de veículos sobre o tabuleiro. Estas ações dinâmicas podem gerar o surgimento das fraturas ou mesmo a sua propagação na estrutura. A correta consideração destes aspectos objetivou o desenvolvimento de um estudo, de forma a avaliar os esforços do tráfego de veículos pesados sobre o tabuleiro. As técnicas para a contagem de ciclos de esforços e a aplicação das regras de dano acumulado foram analisadas através das curvas S-N de diversas normas estudadas. A ponte rodoviária investigada é constituída por quatro vigas longitudinais, três transversinas e por um tabuleiro de concreto armado. O modelo computacional, desenvolvido para a análise dinâmica da ponte, foi concebido com base no emprego de técnicas usuais de discretização através do método dos elementos finitos. O modelo estrutural da obra de arte rodoviária estudada foi simulado com base no emprego de elementos finitos sólidos tridimensionais. Os veículos são representados a partir de sistemas massa-mola-amortecedor. O tráfego dessas viaturas é considerado mediante a simulação de comboios semi-infinitos, deslocando-se com velocidade constante sobre o tabuleiro da ponte. As conclusões deste trabalho versam acerca da vida útil de serviço dos elementos estruturais de pontes rodoviárias de concreto armado submetidas às ações dinâmicas provenientes do tráfego de veículos pesados sobre o tabuleiro.

Análise da dissipação das tensões em dentes humanos restaurados com facetas laminadas de cerâmica, com três tipos de preparos, através do método dos elementos finitos

O objetivo deste trabalho é analisar in vitro a dissipação de tensões em incisivos centrais superiores humanos restaurados com facetas de cerâmica feldspática, através da análise do método dos elementos finitos, considerando cargas funcionais de mastigação e corte dos alimentos, em função de três tipos de preparos utilizados: sem proteção incisal; com proteção incisal em ângulo e com proteção incisal em degrau palatino. Foram utilizadas modelagens bidimensionais de um incisivo central superior e suas estruturas de suporte, simulando três situações: (Primeira modelagem) incisivo central superior com desgaste vestibular (em forma de janela); (Segunda modelagem) incisivo central superior com desgaste vestibular e proteção incisal em plano inclinado; (Terceira modelagem) incisivo central superior com desgaste vestibular, e proteção incisal com degrau palatino. Foi considerada uma carga (P=100N) com uma inclinação de 45 concentrada, simulando a região de contato do incisivo central inferior com o superior durante a mastigação e uma na região de contato topo a topo dos incisivos superior e inferior, simulando o corte dos alimentos. Após a análise dos dados obtidos pela distribuição de tensões, pode-se concluir que quanto à dissipação das tensões em todo o sistema proposto, com a aplicação de carga em 45, não foram observadas mudanças no estado tensional nos três diferentes preparos. Quando foi aplicada carga vertical, simulando o contato de topo, houve variação no estado tensional no sistema do dente com preparo em janela. Nas facetas, com a aplicação de carga em 45, nos preparos em janela e com proteção incisal em plano inclinado o resultado foi semelhante nos valores tensionais enquanto, nas facetas em dentes preparados com proteção incisal com degrau palatino, a distribuição foi mais homogênea tendo valores superiores, mostrando que o abraçamento do dente diminuiu a flexão.

Experimental, Numerical and Analytical Studies of the MHD-driven plasma jet, instabilities and waves

This thesis describes a series of experimental, numerical, and analytical studies involving the Caltech magnetohydrodynamically (MHD)-driven plasma jet experiment. The plasma jet is created via a capacitor discharge that powers a magnetized coaxial planar electrodes system. The jet is collimated and accelerated by the MHD forces.

We present three-dimensional ideal MHD finite-volume simulations of the plasma jet experiment using an astrophysical magnetic tower as the baseline model. A compact magnetic energy/helicity injection is exploited in the simulation analogous to both the experiment and to astrophysical situations. Detailed analysis provides a comprehensive description of the interplay of magnetic force, pressure, and flow effects. We delineate both the jet structure and the transition process that converts the injected magnetic energy to other forms.

When the experimental jet is sufficiently long, it undergoes a global kink instability and then a secondary local Rayleigh-Taylor instability caused by lateral acceleration of the kink instability. We present an MHD theory of the Rayleigh-Taylor instability on the cylindrical surface of a plasma flux rope in the presence of a lateral external gravity. The Rayleigh-Taylor instability is found to couple to the classic current-driven instability, resulting in a new type of hybrid instability. The coupled instability, produced by combination of helical magnetic field, curvature of the cylindrical geometry, and lateral gravity, is fundamentally different from the classic magnetic Rayleigh-Taylor instability occurring at a two-dimensional planar interface.

In the experiment, this instability cascade from macro-scale to micro-scale eventually leads to the failure of MHD. When the Rayleigh-Taylor instability becomes nonlinear, it compresses and pinches the plasma jet to a scale smaller than the ion skin depth and triggers a fast magnetic reconnection. We built a specially designed high-speed 3D magnetic probe and successfully detected the high frequency magnetic fluctuations of broadband whistler waves associated with the fast reconnection. The magnetic fluctuations exhibit power-law spectra. The magnetic components of single-frequency whistler waves are found to be circularly polarized regardless of the angle between the wave propagation direction and the background magnetic field.

Diseño estructural de un basculante de moto de competición

[ES]En este trabajo se plantea una propuesta de diseño estructural de un basculante de una moto de competición de 250cc y 4T adecuado a las normas de la competición MotoStudent entre universidades. En primer lugar se estudian los distintos tipos de basculante existentes y se ponderan sus ventajas e inconvenientes. Posteriormente se elige uno de los tipos en base a las características de la moto y se desarrolla. La herramienta utilizada para el desarrollo es el programa PTC Creo 2.0. Primero se realiza el diseño en CAD de la pieza y después se estudia el diseño con el método de elementos finitos y se modifica hasta alcanzar los requisitos necesarios. Finalmente se añaden elementos auxiliares necesarios y se realizan los planos de la pieza.

I. Numerical solution of exact pair equations. II. Numerical solution of first-order pair equations

Part I

Numerical solutions to the S-limit equations for the helium ground state and excited triplet state and the hydride ion ground state are obtained with the second and fourth difference approximations. The results for the ground states are superior to previously reported values. The coupled equations resulting from the partial wave expansion of the exact helium atom wavefunction were solved giving accurate S-, P-, D-, F-, and G-limits. The G-limit is -2.90351 a.u. compared to the exact value of the energy of -2.90372 a.u.

Part II

The pair functions which determine the exact first-order wavefunction for the ground state of the three-electron atom are found with the matrix finite difference method. The second- and third-order energies for the (1s1s)¹S, (1s2s)³S, and (1s2s)¹S states of the two-electron atom are presented along with contour and perspective plots of the pair functions. The total energy for the three-electron atom with a nuclear charge Z is found to be E(Z) = -1.125•Z² +1.022805•Z-0.408138-0.025515•(1/Z)+O(1/Z²)a.u.