A fixed-grid b-spline finite element technique for fluid-structure interaction

We present a fixed-grid finite element technique for fluid-structure interaction problems involving incompressible viscous flows and thin structures. The flow equations are discretised with isoparametric b-spline basis functions defined on a logically Cartesian grid. In addition, the previously proposed subdivision-stabilisation technique is used to ensure inf-sup stability. The beam equations are discretised with b-splines and the shell equations with subdivision basis functions, both leading to a rotation-free formulation. The interface conditions between the fluid and the structure are enforced with the Nitsche technique. The resulting coupled system of equations is solved with a Dirichlet-Robin partitioning scheme, and the fluid equations are solved with a pressure-correction method. Auxiliary techniques employed for improving numerical robustness include the level-set based implicit representation of the structure interface on the fluid grid, a cut-cell integration algorithm based on marching tetrahedra and the conservative data transfer between the fluid and structure discretisations. A number of verification and validation examples, primarily motivated by animal locomotion in air or water, demonstrate the robustness and efficiency of our approach. © 2013 John Wiley & Sons, Ltd.

An important pitfall of pseudo-static finite element analysis

Finite Element (FE) pseudo-static analysis can provide a good compromise between simplified methods of dynamic analysis and time domain analysis. The pseudo-static FE approach can accurately model the in situ, stresses prior to seismic loading (when it follows a static analysis simulating the construction sequence) is relatively simple and not as computationally expensive as the time domain approach. However this method should be used with caution as the results can be sensitive to the choice of the mesh dimensions. In this paper two simple examples of pseudo-static finite element analysis are examined parametrically, a homogeneous slope and a cantilever retaining wall, exploring the sensitivity of the pseudo-static analysis results on the adopted mesh size. The mesh dependence was found to be more pronounced for problems with high critical seismic coefficients values (e.g. gentle slopes or small walls), as in these cases a generalised layer failure mechanism is developed simultaneously with the slope or wall mechanism. In general the mesh width was found not to affect notably the predicted value of critical seismic coefficient but to have a major impact on the predicted movements. © 2012 Elsevier Ltd.

A discrete dislocation analysis of fatigue crack growth

Cyclic loading of a plane strain mode I crack under small scale yielding is analyzed using discrete dislocation dynamics. The dislocations are all of edge character, and are modeled as line singularities in an elastic solid. At each stage of loading, superposition is used to represent the solution in terms of solutions for edge dislocations in a half-space and a non-singular complementary solution that enforces the boundary conditions, which is obtained from a linear elastic, finite element solution. The lattice resistance to dislocation motion, dislocation nucleation, dislocation interaction with obstacles and dislocation annihilation are incorporated into the formulation through a set of constitutive rules. An irreversible relation between the opening traction and the displacement jump across a cohesive surface ahead of the initial crack tip is also specified, which permits crack growth to emerge naturally. It is found that crack growth can occur under cyclic loading conditions even when the peak stress intensity factor is smaller than the stress intensity required for crack growth under monotonic loading conditions; however below a certain threshold value of ΔKI no crack growth was seen.

Plastic limit solution for the response of plate anchors under combined translational and torsional undrained loading

Plate anchors are increasingly being used to moor large floating offshore structures in deep and ultradeep water. These facilities impart substantial vertical uplift loading to plate anchors. However, extreme operating conditions such as hurricane loading often result in partial system failures, with significant change in the orientation of the remaining intact mooring lines. The purpose of this study is to investigate the undrained pure translational (parallel to plate) and torsional bearing capacity of anchor plates idealized as square and rectangular shaped plates. Moreover, the interaction response of plate anchors under combined translational and torsional loading is studied using a modified plastic limit analysis (PLA) approach. The previous PLA formulation which did not account for shear-normal force interaction on the vertical end faces of the plate provides an exact solution to the idealized problem of an infinitely thin plate but only an approximate solution to the problem of a plate of finite thickness. This is also confirmed by the three-dimensional finite element (FE) results, since the PLA values exceed FE results as the thickness of the plate increases. By incorporating the shear-normal interaction relationship in the modified solution, the torsional bearing capacity factors, as well as the plate interaction responses are enhanced as they show satisfactory agreement with the FE results. The interaction relationship is then obtained for square and rectangular plates of different aspect ratios and thicknesses. The new interaction relationships could also be used as an associated plastic failure locus for combined shear and torsional loading to predict plastic displacements and rotations in translational and torsional loading modes as well. Copyright © 2011 by the International Society of Offshore and Polar Engineers (ISOPE).