Non-linear refined finite element elastic analysis of steel frames with generalized transverse member loads

In the finite element modelling of steel frames, external loads usually act along the members rather than at the nodes only. Conventionally, when a member is subjected to these transverse loads, they are converted to nodal forces which act at the ends of the elements into which the member is discretised by either lumping or consistent nodal load approaches. For a contemporary geometrically non-linear analysis in which the axial force in the member is large, accurate solutions are achieved by discretising the member into many elements, which can produce unfavourable consequences on the efficacy of the method for analysing large steel frames. Herein, a numerical technique to include the transverse loading in the non-linear stiffness formulation for a single element is proposed, and which is able to predict the structural responses of steel frames involving the effects of first-order member loads as well as the second-order coupling effect between the transverse load and the axial force in the member. This allows for a minimal discretisation of a frame for second-order analysis. For those conventional analyses which do include transverse member loading, prescribed stiffness matrices must be used for the plethora of specific loading patterns encountered. This paper shows, however, that the principle of superposition can be applied to the equilibrium condition, so that the form of the stiffness matrix remains unchanged with only the magnitude of the loading being needed to be changed in the stiffness formulation. This novelty allows for a very useful generalised stiffness formulation for a single higher-order element with arbitrary transverse loading patterns to be formulated. The results are verified using analytical stability function studies, as well as with numerical results reported by independent researchers on several simple structural frames.

Experimental and finite element analysis to identify the source of vibration of a coach

This paper presents an investigation on the cause of severe vibration problem of a coach with four-cylinder engine running at an idle state using vibration and impact hammer modal experiments to obtain the main vibration frequency components and the natural characteristics of the coach. The vibration results indicate that the main vibration component comes from the vibration transmitted from the engine to the chassis frame, which is closely related with the engine idle speed. Based on structural simulation analysis of the coach’s chassis frame and comparison with modal testing, the coach severe vibration problem was due to coupling resonance between the engine idle frequency and the fourth natural frequency of the chassis frame. The solution to eliminate the vibration problem is provided by changing the local structure stiffness of the chassis frame. The contribution of this paper lies in providing a solution to solve similar problems.

Finite element analysis of the in-plane shear behaviour of masonry panels confined with reinforced grouted cores

A combined experimental and numerical program was conducted to study the in-plane shear behaviour of hollow concrete masonry panels containing reinforced grout cores. This paper is focused on the numerical program. A two dimensional macromodelling strategy was used to simulate the behaviour of the confined masonry (CM) shear panels. Both the unreinforced masonry and the confining element were modelled using macromasonry properties and the steel reinforcement was modelled as an embedded truss element located within the grout using perfectly bonded constraint. The FE model reproduced key behaviours observed in the experiments, including the shear strength, the deformation and the crack patterns of the unconfined and confined masonry panels. The predictions of the validated model were used to evaluate the existing in-plane shear expressions available in the national masonry standards and research publications.

Determination of strain-rate-dependent mechanical behavior of living and fixed osteocytes and chondrocytes using atomic force microscopy and inverse finite element analysis

The aim of this paper is to determine the strain-rate-dependent mechanical behavior of living and fixed osteocytes and chondrocytes, in vitro. Firstly, Atomic Force Microscopy (AFM) was used to obtain the force-indentation curves of these single cells at four different strain-rates. These results were then employed in inverse finite element analysis (FEA) using Modified Standard neo-Hookean Solid (MSnHS) idealization of these cells to determine their mechanical properties. In addition, a FEA model with a newly developed spring element was employed to accurately simulate AFM evaluation in this study. We report that both cytoskeleton (CSK) and intracellular fluid govern the strain-rate-dependent mechanical property of living cells whereas intracellular fluid plays a predominant role on fixed cells’ behavior. In addition, through the comparisons, it can be concluded that osteocytes are stiffer than chondrocytes at all strain-rates tested indicating that the cells could be the biomarker of their tissue origin. Finally, we report that MSnHS is able to capture the strain-rate-dependent mechanical behavior of osteocyte and chondrocyte for both living and fixed cells. Therefore, we concluded that the MSnHS is a good model for exploration of mechanical deformation responses of single osteocytes and chondrocytes. This study could open a new avenue for analysis of mechanical behavior of osteocytes and chondrocytes as well as other similar types of cells.

Adaptive finite element analysis of elliptic problems based on bubble-type local mesh generation

A new mesh adaptivity algorithm that combines a posteriori error estimation with bubble-type local mesh generation (BLMG) strategy for elliptic differential equations is proposed. The size function used in the BLMG is defined on each vertex during the adaptive process based on the obtained error estimator. In order to avoid the excessive coarsening and refining in each iterative step, two factor thresholds are introduced in the size function. The advantages of the BLMG-based adaptive finite element method, compared with other known methods, are given as follows: the refining and coarsening are obtained fluently in the same framework; the local a posteriori error estimation is easy to implement through the adjacency list of the BLMG method; at all levels of refinement, the updated triangles remain very well shaped, even if the mesh size at any particular refinement level varies by several orders of magnitude. Several numerical examples with singularities for the elliptic problems, where the explicit error estimators are used, verify the efficiency of the algorithm. The analysis for the parameters introduced in the size function shows that the algorithm has good flexibility.

Finite element analysis of CFRP strengthened steel hollow sections under tension

This paper presents a nonlinear finite element (FE) model for the analysis of very high strength (VHS) steel hollow sections wrapped by high modulus carbon fibre rein forced polymer (CFRP) sheets. The bond strength of CFRP wrapped VHS circular steel hollow section under tension is investigated using the FE model. The three dimensional FE model by Nonlinear static analysis has been carried out by Strand 7 finite element software. The model is validated by the experimental data obtained from Fawzia et al [1]. A detail parametric study has been performed to examine the effect of number of CFRP layers, different diameters of VHS steel tube and different bond lengths of CFRP sheet. The analytical model developed by Fawzia et al. [1] has been used to determine the load carrying capacity of different diameters of CFRP strengthened VHS steel tube by using the capacity from each layer of CFRP sheet. The results from FE model have found in reasonable agreement with the analytical model developed by Fawzia et al [1]. This validation was necessary because the analytical model by Fawzia et al [1] was developed by using only one diameter of VHS steel tube and fixed (five) number of CFRP layers. It can be concluded that the developed analytical model is valid for CFRP strengthened VHS steel tubes with diameter range of 38mm to 100mm and CFRP layer range of 3 to 5 layers. Based on the results it can also be concluded that the effective bond length is consistent for different diameters of steel tubes and different layers of CFRP. Three layers of CFRP is considered most effective wrapping scheme due to the cost effectiveness. Finally the distribution of longitudinal and hoop stress has been determined by the finite element model for different diameters of CFRP strengthened VHS steel tube.

Novel non-linear elastic structural analysis with generalised transverse element loads using a refined finite element

In the finite element modelling of structural frames, external loads such as wind loads, dead loads and imposed loads usually act along the elements rather than at the nodes only. Conventionally, when an element is subjected to these general transverse element loads, they are usually converted to nodal forces acting at the ends of the elements by either lumping or consistent load approaches. In addition, it is especially important for an element subjected to the first- and second-order elastic behaviour, to which the steel structure is critically prone to; in particular the thin-walled steel structures, when the stocky element section may be generally critical to the inelastic behaviour. In this sense, the accurate first- and second-order elastic displacement solutions of element load effect along an element is vitally crucial, but cannot be simulated using neither numerical nodal nor consistent load methods alone, as long as no equilibrium condition is enforced in the finite element formulation, which can inevitably impair the structural safety of the steel structure particularly. It can be therefore regarded as a unique element load method to account for the element load nonlinearly. If accurate displacement solution is targeted for simulating the first- and second-order elastic behaviour on an element on the basis of sophisticated non-linear element stiffness formulation, the numerous prescribed stiffness matrices must indispensably be used for the plethora of specific transverse element loading patterns encountered. In order to circumvent this shortcoming, the present paper proposes a numerical technique to include the transverse element loading in the non-linear stiffness formulation without numerous prescribed stiffness matrices, and which is able to predict structural responses involving the effect of first-order element loads as well as the second-order coupling effect between the transverse load and axial force in the element. This paper shows that the principle of superposition can be applied to derive the generalized stiffness formulation for element load effect, so that the form of the stiffness matrix remains unchanged with respect to the specific loading patterns, but with only the magnitude of the loading (element load coefficients) being needed to be adjusted in the stiffness formulation, and subsequently the non-linear effect on element loadings can be commensurate by updating the magnitude of element load coefficients through the non-linear solution procedures. In principle, the element loading distribution is converted into a single loading magnitude at mid-span in order to provide the initial perturbation for triggering the member bowing effect due to its transverse element loads. This approach in turn sacrifices the effect of element loading distribution except at mid-span. Therefore, it can be foreseen that the load-deflection behaviour may not be as accurate as those at mid-span, but its discrepancy is still trivial as proved. This novelty allows for a very useful generalised stiffness formulation for a single higher-order element with arbitrary transverse loading patterns to be formulated. Moreover, another significance of this paper is placed on shifting the nodal response (system analysis) to both nodal and element response (sophisticated element formulation). For the conventional finite element method, such as the cubic element, all accurate solutions can be only found at node. It means no accurate and reliable structural safety can be ensured within an element, and as a result, it hinders the engineering applications. The results of the paper are verified using analytical stability function studies, as well as with numerical results reported by independent researchers on several simple frames.

Finite element analysis investigation of response of pumpkin tissue under uniaxial loading

Finite element analysis (FEA) models of uniaxial loading of pumpkin peel and flesh tissues were developed and validated using experimental results. The tensile model was developed for both linear elastic and plastic material models, the compression model was develop d only with the plastic material model. The outcomes of force versus time curves obtained from FEA models followed similar pattern to the experimental curves however the curve resulted with linear elastic material properties had a higher difference with the experimental curves. The values of predicted forces were determined and compared with the experimental curve. An error indicator was introduced and computed for each case and compared. Additionally Root Mean Square Error (RMSE) values were also calculated for each model and compared. The results of modelling were used to develop material model for peel and flesh tissues in FEA modelling of mechanical peeling of tough skinned vegetables.

Finite element analysis of deep transverse metatarsal ligaments mechanical response during landing

The deep transverse metatarsal ligaments (DTML) play an important role in stabilizing the metatarsal bones and manipulating foot transverse arch deformation. However, the biomechanical research about DTML in the foot maneuver is quite few. Due to the difficulties and lack of better measurement technology for these ligaments experimental monitor, the load transfer mechanism and internal stress state also hadn't been well addressed. The purpose of this study was to develop a detailing foot finite element model including DTML tissues, to investigate the mechanical response of DTML during the landing condition. The DTML was considered as hyperelastic material model was used to represent the nonlinear and nearly incompressible nature of the ligament tissue. From the simulation results, it is clearly to find that the peak maiximal principal stress of DTML was between the third and fourth metatarsals. Meanwhile, it seems the DTML in the middle position experienced higher tension than the sides DTML.

Utility of magnetic resonance imaging-based finite element analysis for the biomechanical stress analysis of hemorrhagic and non-hemorrhagic carotid plaques

Background: Biomechanical stress analysis has been used for plaque vulnerability assessment. The presence of plaque hemorrhage (PH) is a feature of plaque vulnerability and is associated with thromboembolic ischemic events. The purpose of the present study was to use finite element analysis (FEA) to compare the stress profiles of hemorrhagic and non-hemorrhagic profiles. Methods and Results: Forty-five consecutive patients who had suffered a cerebrovascular ischemic event with an underlying carotid artery disease underwent high-resolution magnetic resonance imaging (MRI) of their symptomatic carotid artery in a 1.5-T MRI system. Axial images were manually segmented for various plaque components and used for FEA. Maximum critical stress (M-CstressSL) for each slice was determined. Within a plaque, the maximum M-CstressSL for each slice of a plaque was selected to represent the maximum critical stress of that plaque (M-CstressPL) and used to compare hemorrhagic and non-hemorrhagic plaques. A total of 62% of plaques had hemorrhage. It was observed that plaques with hemorrhage had significantly higher stress (M-CstressPL) than plaques without PH (median [interquartile range]: 315 kPa [247-434] vs. 200 kPa [171-282], P=0.003). Conclusions: Hemorrhagic plaques have higher biomechanical stresses than non-hemorrhagic plaques. MRI-based FEA seems to have the potential to assess plaque vulnerability.

Finite element analysis of vulnerable atherosclerotic plaques: A comparison of mechanical stresses within carotid plaques of acute and recently symptomatic patients with carotid artery disease

Objectives: There is considerable evidence that patients with carotid artery stenosis treated immediately after the ischaemic cerebrovascular event have a better clinical outcome than those who have delayed treatment. Biomechanical assessment of carotid plaques using high-resolution MRI can help examine the relationship between the timing of carotid plaque symptomology and maximum simulated plaque stress concentration. Methods: Fifty patients underwent high-resolution multisequence in vivo MRI of their carotid arteries. Patients with acute symptoms (n=25) underwent MRI within 72 h of the onset of ischaemic cerebrovascular symptoms, whereas recently symptomatic patients (n=25) underwent MRI from 2 to 6 weeks after the onset of symptoms. Stress analysis was performed based on the geometry derived from in vivo MRI of the symptomatic carotid artery at the point of maximum stenosis. The peak stresses within the plaques of the two groups were compared. Results: Patient demographics were comparable for both groups. All the patients in the recently symptomatic group had severe carotid stenosis in contrast to patients with acute symptoms who had predominantly mild to moderate carotid stenosis. The simulated maximum stresses in patients with acute symptoms was significantly higher than in recently symptomatic patients (median (IQR): 313310 4 dynes/cm 2 (295 to 382) vs 2523104 dynes/cm 2 (236 to 311), p=0.02). Conclusions: Patients have extremely unstable, high-risk plaques, with high stresses, immediately after an acute cerebrovascular event, even at lower degrees of carotid stenoses. Biomechanical stress analysis may help us refine our risk-stratification criteria for the management of patients with carotid artery disease in future.

Correlation of carotid atheromatous plaque inflammation with biomechanical stress: Utility of USPIO enhanced MR imaging and finite element analysis

Objective: The aim of this study was to explore whether there is a relationship between the degree of MR-defined inflammation using ultra small super-paramagnetic iron oxide (USPIO) particles, and biomechanical stress using finite element analysis (FEA) techniques, in carotid atheromatous plaques. Methods and Results: 18 patients with angiographically proven carotid stenoses underwent multi-sequence MR imaging before and 36 h after USPIO infusion. T2 * weighted images were manually segmented into quadrants and the signal change in each quadrant normalised to adjacent muscle was calculated after USPIO administration. Plaque geometry was obtained from the rest of the multi-sequence dataset and used within a FEA model to predict maximal stress concentration within each slice. Subsequently, a new statistical model was developed to explicitly investigate the form of the relationship between biomechanical stress and signal change. The Spearman's rank correlation coefficient for USPIO enhanced signal change and maximal biomechanical stress was -0.60 (p = 0.009). Conclusions: There is an association between biomechanical stress and USPIO enhanced MR-defined inflammation within carotid atheroma, both known risk factors for plaque vulnerability. This underlines the complex interaction between physiological processes and biomechanical mechanisms in the development of carotid atheroma. However, this is preliminary data that will need validation in a larger cohort of patients.

Identifying vulnerable carotid plaques in vivo using high resolution magnetic resonance imaging-based finite element analysis

Object. Individuals with carotid atherosclerosis develop symptoms following rupture of vulnerable plaques. Biomechanical stresses within this plaque may increase vulnerability to rupture. In this report the authors describe the use of in vivo carotid plaque imaging and computational mechanics to document the magnitude and distribution of intrinsic plaque stresses. Methods. Ten (five symptomatic and five asymptomatic) individuals underwent plaque characterization magnetic resonance (MR) imaging. Plaque geometry and composition were determined by multisequence review. Intrinsic plaque stress profiles were generated from 3D meshes by using finite element computational analysis. Differences in principal (shear) stress between normal and diseased sections of the carotid artery and between symptomatic and asymptomatic plaques were noted. Results. There was a significant difference in peak principal stress between diseased and nondiseased segments of the artery (mean difference 537.65 kPa, p < 0.05). Symptomatic plaques had higher mean stresses than asymptomatic plaques (627.6 kPa compared with 370.2 kPa, p = 0.05), which were independent of luminal stenosis and plaque composition. Conclusions. Significant differences in plaque stress exist between plaques from symptomatic individuals and those from asymptomatic individuals. The MR imaging-based computational analysis may therefore be a useful aid to identification of vulnerable plaques in vivo.