Research on Bearing Capacity Calculation Method of Large Diameter Cylinder in Soft Clay

A large diameter cylinder inserted in soils is a new type of engineering structures used in offshore and port engineering. The mechanism of its bearing capacity and the analysis of its stability are important to its design and applications. In this paper, the finite element method is used to analyze the reacting forces of the soft soil foundation on the structure under the wave action. A simplified method is proposed, based on the plastic limit method, for the safety and stability analysis. Our analysis shows that the assumptions made in this paper and the mechanism used are reasonable, and the results obtained are appropriate. The calculation method is very efficient and can be used to evaluate main parameters of the structure in its preliminary designs.

Computation of flow over a rotating body on unstructured Chimera mesh

Flow around moving boundary is ubiquitous in engineering applications. To increse the efficienly of the algorithm to handle moving boundaries is still a major challenge in Computational Fluid Dynamics (CFD). The Chimera grid method is one type of method to handle moving boundaries. A concept of domain de-composition has been proposed in this paper. In this method, sub-domains are meshed independently and governing equations are also solved separately on them. The Chimera grid method was originally used only on structured (curvilinear) meshes. However, in a problem which involves both moving boundary and complex geometry, the number of sub-domains required in a traditional (structured) Chimera method becomes fairly large. Thus the time required in the interior boundary locating, link-building and data exchanging also increases. The use of unstructured Chimera grid can reduce the time consumption significantly by the reduction of domain(block) number. Generally speaking, unstructured Chimera grid method has not been developed. In this paper, a well-known pressure correction scheme - SIMPLEC is modified and implemented on unstructured Chimera mesh. A new interpolation scheme regarding the pressure correction is proposed to prevent the possible decoupling of pressure. A moving-mesh finite volume approach is implemented in an inertial reference frame. This approach is then used to compute incompressible flow around a rotating circular and elliptic cylinder. These numerical examples demonstrate the capability of the proposed scheme in handling moving boundaries. The numerical results are in good agreement with other experimental and computational data in literature. The method proposed in this paper can be efficiently applied to more challenge cases such as free-falling objects or heavy particles in fluid.

A simple method to simulate shrinkage-induced cracking in cement-based composites by lattice-type modeling

A new numerical procedure is proposed to investigate cracking behaviors induced by mismatch between the matrix phase and aggregates due to matrix shrinkage in cement-based composites. This kind of failure processes is simplified in this investigation as a purely spontaneous mechanical problem, therefore, one main difficulty during simulating the phenomenon lies that no explicit external load serves as the drive to propel development of this physical process. As a result, it is different from classical mechanical problems and seems hard to be solved by using directly the classical finite element method (FEM), a typical kind of "load -> medium -> response" procedures. As a solution, the actual mismatch deformation field is decomposed into two virtual fields, both of which can be obtained by the classical FEM. Then the actual response is obtained by adding together the two virtual displacement fields based on the principle of superposition. Then, critical elements are detected successively by the event-by-event technique. The micro-structure of composites is implemented by employing the generalized beam (GB) lattice model. Numerical examples are given to show the effectiveness of the method, and detailed discussions are conducted on influences of material properties.

A surface energy model and application to mechanical behavior analysis of single crystals at sub-micron scale

Size effects of mechanical behaviors of materials are referred to the variation of the mechanical behavior due to the sample sizes changing from macroscale to micro-/nanoscales. At the micro-/nanoscale, since sample has a relatively high specific surface area (SSA) (ratio of surface area to volume), the surface although it is often neglected at the macroscale, becomes prominent in governing the energy effect, although it is often neglected at the macroscale, becomes prominent in governing the mechanical behavior. In the present research, a continuum model considering the surface energy effect is developed through introducing the surface energy to total potential energy. Simultaneously, a corresponding finite element method is developed. The model is used to analyze the axial equilibrium strain problem for a Cu nanowire at the external loading-free state. As another application of the model, from dimensional analysis, the size effects of uniform compression tests on the microscale cylinder specimens for Ni and Au single crystals are analyzed and compared with experiments in literatures. (C) 2009 Elsevier B.V. All rights reserved.

An Instrumented Indentation Method for Young's Modulus Measurement with Accuracy Estimation

The relationships between indentation responses and Young's modulus of an indented material were investigated by employing dimensional analysis and finite element method. Three representative tip bluntness geometries were introduced to describe the shape of a real Berkovich indenter. It was demonstrated that for each of these bluntness geometries, a set of approximate indentation relationships correlating the ratio of nominal hardness/reduced Young's modulus H (n) /E (r) and the ratio of elastic work/total work W (e)/W can be derived. Consequently, a method for Young's modulus measurement combined with its accuracy estimation was established on basis of these relationships. The effectiveness of this approach was verified by performing nanoindentation tests on S45C carbon steel and 6061 aluminum alloy and microindentation tests on aluminum single crystal, GCr15 bearing steel and fused silica.

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.

Crank-Nicholson method for rate equations in powder random lasers

3rd International Conference on Mathematical Modeling in Physical Sciences (IC-MSQUARE 2014)

Otimização de torre de aço para aerogerador eólico.

Diversas formas de geração de energia vêm sendo desenvolvidas com o objetivo de oferecer alternativas ecologicamente corretas. Neste contexto, a energia eólica vem se destacando na região Nordeste do Brasil, devido ao grande potencial dos ventos da região. As torres, que representam parcela significativa do custo total do sistema, tendem a crescer buscando ventos mais fortes e permitindo assim a utilização de aerogeradores com maior capacidade de geração de energia. Este trabalho tem como objetivo formular um modelo de otimização de torres tubulares de aço, para aerogeradores eólicos. Busca-se minimizar o volume total (custo, indiretamente), tendo como variáveis de projeto as espessuras da parede da torre. São impostas restrições relativas à frequência natural e ao comportamento estrutural (tensão e deslocamento máximo de acordo com recomendações da norma Europeia). A estrutura da torre é modelada com base no Método dos Elementos Finitos e o carregamento atuante na estrutura inclui os pesos da torre, do conjunto de equipamentos instalados no topo (aerogerador), e o efeito estático da ação do vento sobre a torre. Para verificação das tensões, deslocamentos e frequências naturais, foram utilizados elementos finitos de casca disponíveis na biblioteca do programa de análise ANSYS. Os modelos de otimização foram também implementados no modulo de otimização do programa ANSYS (design optimization), que utiliza técnicas matemáticas em um processo iterativo computadorizado até que um projeto considerado ótimo seja alcançado. Nas aplicações foram usados os métodos de aproximação por subproblemas e o método de primeira ordem. Os resultados obtidos revelam que torres para aerogeradores merecem atenção especial, em relação à concepção do projeto estrutural, sendo que seu desempenho deve ser verificado através de metodologias completas que englobem além das análises clássicas (estáticas e dinâmicas), incluam também as análises de otimização.

Application of domain decomposition and transmission line modelling techniques to 2D, time-domain, finite element problems

In this paper a recently published finite element method, which combines domain decomposition with a novel technique for solving nonlinear magnetostatic finite element problems is described. It is then shown how the method can be extended to, and optimised for, the solution of time-domain problems. © 1999 IEEE.

A new Lagrangian-Eulerian shell-fluid coupling algorithm based on level sets

We propose a computational method for the coupled simulation of a compressible flow interacting with a thin-shell structure undergoing large deformations. An Eulerian finite volume formulation is adopted for the fluid and a Lagrangian formulation based on subdivision finite elements is adopted for the shell response. The coupling between the fluid and the solid response is achieved via a novel approach based on level sets. The basic approach furnishes a general algorithm for coupling Lagrangian shell solvers with Cartesian grid based Eulerian fluid solvers. The efficiency and robustness of the proposed approach is demonstrated with a airbag deployment simulation. It bears emphasis that in the proposed approach the solid and the fluid components as well as their coupled interaction are considered in full detail and modeled with an equivalent level of fidelity without any oversimplifying assumptions or bias towards a particular physical aspect of the problem.

Towards jet flow LES of conceptual nozzles for acoustic predictions

Progress in simulating chevron nozzle jet flows using ILES/RANS-ILES approaches and using the Ffowcs Williams and Hawkings (FW-H) surface integral method to predict the radiated far field sound is presented in this paper. With the focus on the realistic chevron geometries, SMC001 and SMC006, coarse and fine meshes are generated in the range of 3∼13 million mesh cells. Throughout this work, to minimize numerical dissipation introduced by mesh quality issues, the hexahedral cell type is used. Numerical simulations are then carried out with cell-vertex and cell-centered codes. Despite the modest grids, mean velocities and turbulent statistics are found to be in reasonable accord with measurements. Also, far field sound levels predicted by the FW-H post processor are encouraging. Copyright © 2008 by the American Institute of Aeronautics and Astronautics, Inc.

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.

Subdivision-stabilised immersed b-spline finite elements for moving boundary flows

An immersed finite element method is presented to compute flows with complex moving boundaries on a fixed Cartesian grid. The viscous, incompressible fluid flow equations are discretized with b-spline basis functions. The two-scale relation for b-splines is used to implement an elegant and efficient technique to satisfy the LBB condition. On non-grid-aligned fluid domains and at moving boundaries, the boundary conditions are enforced with a consistent penalty method as originally proposed by Nitsche. In addition, a special extrapolation technique is employed to prevent the loss of numerical stability in presence of arbitrarily small cut-cells. The versatility and accuracy of the proposed approach is demonstrated by means of convergence studies and comparisons with previous experimental and computational investigations.

Active control of high-frequency vibration: Optimisation using the hybrid modelling method

This work presents active control of high-frequency vibration using skyhook dampers. The choice of the damper gain and its optimal location is crucial for the effective implementation of active vibration control. In vibration control, certain sensor/actuator locations are preferable for reducing structural vibration while using minimum control effort. In order to perform optimisation on a general built-up structure to control vibration, it is necessary to have a good modelling technique to predict the performance of the controller. The present work exploits the hybrid modelling approach, which combines the finite element method (FEM) and statistical energy analysis (SEA) to provide efficient response predictions at medium to high frequencies. The hybrid method is implemented here for a general network of plates, coupled via springs, to allow study of a variety of generic control design problems. By combining the hybrid method with numerical optimisation using a genetic algorithm, optimal skyhook damper gains and locations are obtained. The optimal controller gain and location found from the hybrid method are compared with results from a deterministic modelling method. Good agreement between the results is observed, whereas results from the hybrid method are found in a significantly reduced amount of time. © 2012 Elsevier Ltd. All rights reserved.

Finite element analysis of stress and strain distributions in InAs/GaAs quantum dots

In this paper, we perform systematic calculations of the stress and strain distributions in InAs/GaAs truncated pyramidal quantum dots (QDs) with different wetting layer (WL) thickness, using the finite element method (FEM). The stresses and strains are concentrated at the boundaries of the WL and QDs, are reduced gradually from the boundaries to the interior, and tend to a uniform state for the positions away from the boundaries. The maximal strain energy density occurs at the vicinity of the interface between the WL and the substrate. The stresses, strains and released strain energy are reduced gradually with increasing WL thickness. The above results show that a critical WL thickness may exist, and the stress and strain distributions can make the growth of QDs a growth of strained three-dimensional island when the WL thickness is above the critical value, and FEM can be applied to investigate such nanosystems, QDs, and the relevant results are supported by the experiments.