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.

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.

Semi-analytical finite element analysis of a strip footing on an elastic reinforced soil

A plane strain elastic interaction analysis of a strip footing resting on a reinforced soil bed has been made by using a combined analytical and finite element method (FEM). In this approach the stiffness matrix for the footing has been obtained using the FEM, For the reinforced soil bed (halfplane) the stiffness matrix has been obtained using an analytical solution. For the latter, the reinforced zone has been idealised as (i) an equivalent orthotropic infinite strip (composite approach) and (ii) a multilayered system (discrete approach). In the analysis, the interface between the strip footing and reinforced halfplane has been assumed as (i) frictionless and (ii) fully bonded. The contact pressure distribution and the settlement reduction have been given for different depths of footing and scheme of reinforcement in soil. The load-deformation behaviour of the reinforced soil obtained using the above modelling has been compared with some available analytical and model test results. The equivalent orthotropic approach proposed in this paper is easy to program and is shown to predict the reinforcing effects reasonably well.

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.

Finite element analysis of laminated anisotropic thin-walled open-section beams

A finite element analysis of thin-walled open-section laminated anisotropic beams is presented herein. A two-noded, 8 degrees of freedom per node thin-walled open-section laminated anisotropic beam finite element has been developed and used. The displacements of the element reference axes are expressed in terms of one-dimensional first order Hermite interpolation polynomials and line member assumptions are invoked in the formulation of the stiffness matrix. The problems of: 1. (a) an isotropic material Z section straight cantilever beam, and 2. (b) a single-layer (0°) composite Z section straight cantilever beam, for which continuum solutions (exact/approximate) are possible, have been solved in order to evaluate the performance of the finite element. Its applicability has been shown by solving the following problems: 3. (c) a two-layer (45°/−45°) composite Z section straight cantilever beam, 4. (d) a three-layer (0°/45°/0°) composite Z section straight cantilever beam.

Finite element analysis of bimodulus composite stiffened thin shells of revolution

This paper is a sequel to the work published by the first and third authors[l] on stiffened laminated shells of revolution made of unimodular materials (materials having identical properties in tension and compression). A finite element analysis of laminated bimodulus composite thin shells of revolution, reinforced by laminated bimodulus composite stiffeners is reported herein. A 48 dot doubly curved quadrilateral laminated anisotropic shell of revolution finite element and it's two compatible 16 dof stiffener finite elements namely: (i) a laminated anisotropic parallel circle stiffener element (PCSE) and (ii) a laminated anisotropic meridional stiffener element (MSE) have been used iteratively. The constitutive relationship of each layer is assumed to depend on whether the fiberdirection strain is tensile or compressive. The true state of strain or stress is realized when the locations of the neutral surfaces in the shell and the stiffeners remain unaltered (to a specified accuracy) between two successive iterations. The solutions for static loading of a stiffened plate, a stiffened cylindrical shell. and a stiffened spherical shell, all made of bimodulus composite materials, have been presented.

Finite element analysis of moving contact in mechanically fastened joints

The Finite Element Method (FEM) has made a number of otherwise intractable problems solvable. An important aspect for achieving an economical and accurate solution through FEM is matching the formulation and the computational organisation to the problem. This was realised forcefully in the present case of the solution of a class of moving contact boundary value problems of fastener joints. This paper deals with the problem of changing contact at the pin-hole interface of a fastener joint. Due to moving contact, the stresses and displacements are nonlinear with load. This would, in general, need an interactive-incremental approach for solution. However, by posing the problem in an inverse way, a solution is sought for obtaining loads to suit given contact configuration. Numerical results are given for typical isotropic and composite plates with rigid pins. Two cases of loading are considered: (i) load applied only at the edges of the plate and (ii) load applied at the pin and reacted at a part of the edge of the plate. Load-contact relationships, compliance and stress-patterns are investigated. This paper clearly demonstrates the simplification achieved by a suitable formulation of the problem. The results are of significance to the design and analysis of fastener joints.

Finite element analysis of moving contact in mechanically fastened joints

The Finite Element Method (FEM) has made a number of otherwise intractable problems solvable. An important aspect for achieving an economical and accurate solution through FEM is matching the formulation and the computational organisation to the problem. This was realised forcefully in the present case of the solution of a class of moving contact boundary value problems of fastener joints. This paper deals with the problem of changing contact at the pin-hole interface of a fastener joint. Due to moving contact, the stresses and displacements are nonlinear with load. This would, in general, need an interactive-incremental approach for solution. However, by posing the problem in an inverse way, a solution is sought for obtaining loads to suit given contact configuration. Numerical results are given for typical isotropic and composite plates with rigid pins. Two cases of loading are considered: (i) load applied only at the edges of the plate and (ii) load applied at the pin and reacted at a part of the edge of the plate. Load-contact relationships, compliance and stress-patterns are investigated. This paper clearly demonstrates the simplification achieved by a suitable formulation of the problem. The results are of significance to the design and analysis of fastener joints.

Finite element analysis of curved anisotropic laminated hollow bars

Curved hollow bars of laminated anisotropic construction are used as structural members in many industries. They are used in order to save weight without loss of stiffness in comparison with solid sections. In this paper are presented the details of the development of the stiffness matrices of laminated anisotropic curved hollow bars under line member assumptions for two typical sections, circular and square. They are 16dof elements which make use of one-dimensional first-order Hermite interpolation polynomials for the description of assumed displacement state. Problems for which analytical or other solutions are available are first solved using these elements. Good agreement was found between the results. In order to show the capability of the element, application is made to carbon fibre reinforced plastic layered anisotropic curved hollow bars.

Finite element analysis of laminated anisotropic beams of bimodulus materials

This paper presents finite element analysis of laminated anisotropic beams of bimodulus materials. The finite element has 16 d.o.f. and uses the displacement field in terms of first order Hermite interpolation polynomials. As the neutral axis position may change from point to point along the length of the beam, an iterative procedure is employed to determine the location of zero strain points along the length. Using this element some problems of laminated beams of bimodulus materials are solved for concentrated loads/moments perpendicular and parallel to the layering planes as well as combined loads.