3D Finite Element Research Articles

Finite element analysis and optimization studies on tibia implant of SS 316L steel and Ti6Al4V alloy

Tibial fractures account for approximately 15% of all fractures, typically resulting from high-energy trauma. A critical surgical approach to treat these fractures involves the fixation of the tibia using a plate with minimally invasive osteosynthesis. The selection and fixation of the implant plate are vital for stabilizing the fracture. This selection is highly dependent on the plate's stability, which is influenced by factors like the stresses generated in the plate due to the load on the bone, as well as the plate's length, thickness, and number of screw holes. Minimizing these stresses is essential to reduce the risk of implant failure, ensuring optimal stress distribution and promoting faster, more effective bone healing. In the present work, the finite element and statistical approach was used to optimize the geometrical parameters of the implant plate made of SS 316L steel and Ti6Al4V alloy. A 3D finite element model was developed for analyzing the stresses and deformation, and implant plates were manufactured to validate the results with the help of an experiment conducted on the universal testing machine. A strong correlation was observed between the experimental and predicted results, with an average error of 8.6% and 8.55% for SS316L and Ti6Al4V alloy, respectively. Further, using the signal-to-noise ratio for the minimum stress condition was applied to identify the optimum parameters of the plate. Finally, regression models were developed to predict the stresses generated in SS316L and Ti6Al4V alloy plates with different input conditions. The statistical model helps us to develop the relation between different geometrical parameters of the Tibia implant plate. As determined by the present work, the parameter most influencing is implant plate length. This outcome will be used to select the implant for a specific patient, resulting in a reduction in implant failure post-surgery.

Finite element simulations of quartz crystal microbalances (QCM): from Sauerbrey to fractional viscoelasticity under water

2D finite element simulations are performed on QCM working in the thickness-shear mode and loaded with different homogeneous films. They include a purely elastic film, a viscoelastic Maxwellian liquid, viscoelastic-Voigt solid, and the fractional viscoelastic (power-law) version of each case. Unlike single-relaxation kind models, fractional viscoelasticity considers the relaxation-time spectrum often found in polymeric materials. The films are tested in air or covered with liquids of different viscosities. Two substrate thicknesses are tested: 100 nm and 500 nm, the latter being close to the condition that promotes the resonance of the adsorbed film. In all cases the simulations are compared with small-load approximation theory (SLA). The 100 nm films follow the theory closely, although significant deviations of the SLA are observed as the overtone number n increases, even in purely elastic films. We also show that it is possible to identify the viscoelastic ‘fingerprint’ of the 100 nm films in air using raw data and Sauerbrey’s equivalent thickness obtained with the QCM in the 3 < n < 13 range. These numerical data are validated by experimental measurements of crosslinked polydimethylsiloxane films with thicknesses ∼150 nm. In contrast, the 500 nm films deviate notoriously from the SLA, for all viscoelastic models and overtones, with the largest deviation observed in the elastic film. When a liquid layer covers the QCM without an adsorbed film, the only overtone that numerically reproduces the theoretical value is the fundamental, n = 1. For n > 1, strong coupling between the solid and liquid is detected, and the original vibration modes of the crystal are altered by the presence of the liquid. Finally, the numerical simulations suggest that it is possible to detect whether a viscoelastic film is formed under a liquid layer using only the information from n = 1. In these film/liquid systems we also observe the so-called missing-mass effect, although the theory and simulations exhibit different levels of impact of such effect when the liquid viscosity is high.

Research on Performance of Interior Permanent Magnet Synchronous Motor with Fractional Slot Concentrated Winding for Electric Vehicles Applications

The fractional-slot, concentrated-winding, interior permanent magnet synchronous motor (FSCW IPMSM) has advantages, such as reducing motor copper consumption, improving flux-weakening capability, and motor fault tolerance, and has certain development potential in application fields such as electric vehicles. However, fractional-slot concentrated-winding motors often contain rich harmonic components due to their winding characteristics, leading to increased motor losses and back electromotive force harmonics, thereby affecting the efficiency and constant power speed regulation range of the motor. Based on this, this article first uses the winding function method to explore the inductance and saliency ratio of the interior permanent magnet synchronous motor with different slot pole combinations in the fractional-slot concentrated- winding of electric vehicles. Secondly, this article will establish a 2D finite element parameterized model to analyze and compare the performance of fractional-slot concentrated-winding motors with different slot pole combinations, including air gap magnetic density, back electromotive force distortion rate, overload multiple, and torque. The structural parameters of the motor were optimized with the objective of minimizing the torque ripple under the constraint of minimizing the average torque reduction. The motor slot width, permanent magnet angle, and permanent magnet pole arc angle were analyzed and optimized. The simulation results showed that 12 slots and 8 poles were the optimal design schemes, providing a theoretical basis for the selection of slot pole coordination in the fractional-slot concentrated-winding interior permanent magnet synchronous motor for electric vehicles.

The role of underground archaeological remains on the seismic response of a historic tower

The Italian territory features plenty of historic masonry buildings, many of which have undergone reconstruction, often atop ancient ruins. The presence of archaeological remains in the shallow layers of the ground is always disregarded in the seismic assessment of masonry buildings. To address this gap, the paper delves into the influence of underground masonry ruins on the seismic response of a bell tower built above a historic soil deposit. A comprehensive 3D finite-element model encompassing the tower, foundation, soil, and underground remains was developed. A parametric analysis was conducted, involving modifications to the geometry, mechanical properties, and location of the buried ruins. From the parametric study it emerged that the buried walls could alter the seismic signal at the base of the structure. The most critical scenario occurs when the buried archaeological ruins have a grid spacing comparable to the characteristic size of the foundation, especially in soft soils. Analysis of displacements, accelerations, and spectral accelerations at different elevations along the tower height revealed that the presence of archaeological ruins could modify the seismic response of the structure built above them.

Strain-gradient and damage failure behavior in particle reinforced heterogeneous matrix composites

A strain gradient plasticity model with strengthening effect and damage effect (SGSED) is developed under the continuum medium framework for particle-reinforced heterogeneous matrix composites (PRHMCs). Boron carbide (B4C) particle-reinforced ultrafine-grained (UFG)/fine-grained (FG) heterogeneous matrix composites (B4Cp/(UFG/FG)) are fabricated, in which the UFG region consisted of carbon nanotubes (CNTs)/6061 aluminum (Al) flakes with grains in the ultrafine range and the FG region is processed via 6061Al. The SGSED model is written into the user subroutines using commercial finite element (FE) calculation software, and the three-dimensional (3D) FE representative volume element (RVE) for B4Cp/(UFG/FG) composites is established, from which the distribution of the interface-affected-zone (IAZ) formed of the strain gradient caused by the uncoordinated deformation of the UFG-FG heterogeneous matrix and reinforced phase-matrix is calculated. The evolution of the strain gradient in the deformation process of composites and the influence of the strain gradient on the progressive damage and crack evolution of composites are analyzed, and the strain gradient strengthening-toughening mechanism of composites is revealed. It is found that the IAZ has a considerable strengthening-toughening effect on the composites, which can reduce stress concentration at the interface between the reinforced phase and the matrix, and slow down the crack propagation of the matrix.

Experimental and numerical investigation of the damage propagation in regularly arrayed short fiber reinforced composite laminates under low‐velocity impact

AbstractThe outstanding mechanical performance and flowability of unidirectionally arrayed chopped strands (UACSs) give them an advantage in the manufacture of engineering structures with complex geometry. For the application of practical structures, the impact responses and damage evolution of material under low‐velocity impact must be investigated in advance. In this study, UACS laminates and continuous carbon fiber laminates with a stacking sequence of [0/90]4s and a thickness of 2 mm were fabricated for low‐velocity impact tests at 7, 11, 15, and 20 J. The impact responses and postimpact intralaminar damage area were analyzed according to the experimental results, including impact load responses and ultrasound C‐scan inspections. Moreover, to predict the damage evolution of the internal structure, the 3D finite element models were constructed in ABAQUS using a progressive damage model (PDM) through a user‐material subroutine VUMAT. Compared with continuous carbon fiber laminates, the dissipated energy of UACS laminates increases by approximately 10.64% and 57.37% for the 15 and 20 J, respectively. However, the intralaminar damage area of UACS laminates decreased by 29.96% and 28.16% at 15 and 20 J, respectively, since the discontinuous slits in UACS laminates can guide damage paths and suppress damage propagation.Highlights The regularly arrayed short fiber reinforced composite was studied. Revealed the low‐velocity impact responses and predicted the damage evolution. UACS and CFRP laminates have comparable impact performance. Illustrated the slits design can suppress the damage to a relatively small area.

A study on solid-shell finite element formulations applied to nonlinear thermoelastic analysis of thin-walled structures

The solid-shell elements have shown great advantages in three-dimensional (3D) finite element (FE) simulation of thin-walled structures. To overcome the numerical lockings for 3D element with a large span-thickness ratio, the assumed natural strain (ANS) method combined by either the hybrid stress (HS) or enhanced assumed strain (EAS) formulations are commonly used. In this work, two types of finite element formulations of the solid-shell elements are developed for nonlinear thermoelastic analysis of thin-walled structures, which are termed as the “ANS+HS” and “ANS+EAS” formulations. Accordingly, two eight-node solid-shell elements (CSSH8 and CSSE8) are developed based on the “ANS+HS” and “ANS+EAS” formulations, respectively. The Green–Lagrange displacement–strain relation is applied to involve the geometrical nonlinearities, and both the thermal expansion and temperature-dependent material properties are considered. Nonlinear thermoelastic equilibrium equations are constructed in the framework of the two types of finite element formulations using the Hellinger–Reissner and conventional variational principles, respectively. The proposed method can easily realize three different coupling analyses for thermal–mechanical loads by modifying the parameters of nonlinear thermoelastic equilibrium equations. Numerical examples demonstrate that the proposed method using CSSH8 and CSSE8 elements traces the nonlinear thermoelastic response of the isogrid stiffened panel as accurately as the shell, solid-shell, and solid elements of ABAQUS. Furthermore, the path-following capability of the proposed method using the CSSH8 with “ANS+HS” formulation is slightly superior to that using the CSSE8 with “ANS+EAS” formulation, but both better than ABAQUS.

Simulation and experimental thermal analysis of ultrasonic vibration-assisted high-efficiency deep grinding of γ-TiAl blade tenon

Gamma titanium-aluminum intermetallic compounds (γ-TiAl) have attracted considerable attention in the aerospace industry because of their outstanding thermal stability and broad range of properties. To analyze the grinding temperature during high-efficiency deep grinding (HEDG) and ultrasonic vibration-assisted high-efficiency deep grinding (UVHEDG) of γ-TiAl materials, both an analytical thermal model and a finite element simulation model were developed. Subsequent comparison trials between HEDG and UVHEDG were conducted to validate the accuracy of the simulation. The results show that the final prediction error can be effectively maintained within 15 %, and the temperature measurement curve fits well when an ultrasonic heat influence factor is incorporated into the 3D finite element model. The introduction of ultrasonic vibrations into the HEDG reduced the maximum grinding temperature by 39.1 %, effectively preventing grinding burns. The highest temperature consistently occurred in the top zone of the blade tenon tooth due to different heat conduction paths caused by the complex profile. This led to a 32.4 % and 31.9 % reduction in the grinding temperature in the bottom zone under HEDG and UVHEDG, respectively. The surface roughness of the machined blade tenon tooth reached 0.664 μm at the top and 0.721 μm at the bottom under UVHEDG, whereas it was 0.735 μm and 0.843 μm under HEDG.

Crushing behavior of GFRP composite-reinforced PVC tubes: Experimental testing and numerical simulation

This paper introduces glass fiber reinforced polymer (GFRP)-reinforced Polyvinyl Chloride (PVC) tubes, both corrugated and non-corrugated, designed as energy absorber devices. The PVC tubes were externally and internally reinforced with GFRP composite oriented at ±45∘ and subjected to quasi-static axial compression tests. Results indicated that all reinforced tubes exhibited significantly higher load-bearing capacity, energy absorption (EA) capability, and crushing force efficiency (CFE) compared to standard PVC tubes. Among the tested specimens, externally reinforced corrugated tubes demonstrated the highest specific energy absorption (SEA), surpassing other configurations by 17.5 kJ/kg when considering both pre- and post-crushing stages combined. However, these corrugated specimens showed instability during crushing, reflected in poor instantaneous crush force efficiency (iCFE) and the lowest iCFE among the composite tubes, with an average decrease of 43.59%.The corrugation notably increased the initial peak load, enhancing energy absorption in the pre-crushing stage without compromising the stability of crush force efficiency. Additionally, the combination of external and internal reinforcement significantly improved CFE and iCFE. Consequently, the PVC tubes combining corrugation with both external and internal reinforcement emerged as the best-performing configuration among all tested tubes.Furthermore, a 3D Finite Element (FE) model was developed using ABAQUS FE code with user-defined subroutines to simulate the crushing process. The constitutive models and numerical procedures employed are detailed. The FE model’s predictions showed a satisfactory correlation with experimental results, providing valuable insights into the crushing mechanics and offering a predictive tool for future design optimizations.

A 3D finite element analysis of biomechanical effects on teeth and bone during true intrusion of anteriors using miniscrews

The primary objective of this study was to investigate the biomechanical effects and stresses on bone, periodontal ligament (PDL), cementum and displacement along X-, Y- and Z-axis during true intrusion of incisors using mini-implants with Finite Element Analysis; the secondary objective of the study was to find out the best method for anterior intrusion in clinical practice to treat anterior deep bite malocclusions. A 3D finite element method was used to simulate true anterior intrusion with sliding mechanics using mini-implants. Two groups were modelled with mini-implants placed distal to lateral incisors for Model1, and below the anteriornasalspine (ANS) for Model2, to achieve intrusion. von Mises stress, principal stress on PDL and alveolar bone, displacements in all 3 planes were determined. Amongst the modalities compared in the present study, the stresses on bone and PDL were showing best behavior for mini-implants placed distal to lateral incisors (Model1). The teeth showed controlled tooth movement in Model1. Maximum stress was found in the cortical bone and in the PDL. Nature of the stress changed from compressive to tensile from cervical area to root apex, concentrating mainly at the apical area. Amongst the modalities compared, the best controlled tooth movements for anterior intrusion to treat anterior deep bite malocclusions, was for mini-implants placed distal to lateral incisors (Model1).

Numerical simulation of thermal and stress fields for multilayer and multi-pass weaving WAAM of magnesium alloy

Purpose A 3D finite element (FE) model based on the double ellipsoidal heat source was developed to investigate the evolution of temperature and stress fields during the multilayer and multi-pass wire and arc additive manufacturing (WAAM) process. This paper aims to investigate the evolution of temperature and stress fields during the multilayer and multi-pass wire and arc additive manufacturing (WAAM) process by developing a 3D finite element (FE) model based on the double ellipsoidal heat source. Design/methodology/approach Experimental thermal cycle curves and residual stresses were obtained by thermocouples and X-ray diffraction, respectively. The validity of the model was verified by the corresponding experimental results. Findings The deposition process of the upper pass led to the partial remelting of the lower deposited pass. The thermal process of the current-deposited pass alleviated the stress concentration in the previous-formed passes. A more uniform temperature distribution could be obtained by using the reciprocating deposition path. Compared to the reciprocating deposition path, the peak values of the transverse and longitudinal tensile residual stresses of the deposited sample under the unidirectional deposition path were reduced by 15 MPa and increased by 13 MPa, respectively. The heat conduction in the deposited passes could be improved by extending the inter-pass cooling time appropriately. With an increase in the inter-pass cooling time, the longitudinal residual stress in the middle region of sample along longitudinal and transverse directions showed increase and decrease–increase trends, respectively, while the transverse residual stress exhibited decrease trend. Originality/value This study enhances the understanding of temperature and stress fields evolution during the multilayer and multi-pass cold metal transfer-WAAM processes of magnesium alloy and provides the reference for parameter optimization.

Hybrid Continuum and Multi‐Particle Simulation of RDX Powder Subjected to Drop Impact

AbstractA variety of energetic material (EM) devices could be designed more efficiently with a more rapid protocol of modeling and testing if computational tools are more predictive. Drop weight impact tester (DWIT) is widely used to experimentally study energetic material response to impact. Drop weight impact on a thin layer of RDX powder is studied here and computationally modeled to compute the likelihood of EM powder ignition from low‐speed impact. The overall objective is to develop a thermomechanical‐based FEM methodology that can qualitatively predict ignition‐likelihood to study iginition mechanics in EM powders and understand how anvil properties and striker impact conditions influence it. The mechanical confinement conditions offered by an anvil and striker type, influence energy transmission and heat localization within the energetic powder. Efficiently modeling both the bulk behavior and the mesoscopic behavior, such as, localized shear and particle‐to‐particle frictional heating, is necessary to design insensitive EM‐devices. A 3D finite element modeling (FEM) approach is described that more accurately models the impact conditions than previously developed 2D models. The results are compared with the experimental results and the 2D simulations. A continuum‐based explicit FEM can only simulate the bulk‐behavior of the powder. To accurately predict hot‐spot ignition mechanisms, meso‐scale multi‐particle models of the EM is needed. For computational efficiency, hybrid models that combines continuum and multi‐particle FEM (MPFEM) were utilized. The computed temperature profiles are compared for two anvil types and demonstrates that an assessment of ignition likelihood is viable with the hybrid methodologies proposed.

Development of a new soil–structure contact stress sensor for underground construction applications

This paper describes the design, development, calibration and validation of a novel soil–structure contact stress sensor. The new sensor design combines a novel operating principle, fibre Bragg grating (FBG) strain sensing and data-driven mapping techniques to create a multi-axis contact stress sensor that is both economical and suitably robust for deployment in underground construction applications. The instrumentation process is informed using a ‘virtual twin’ of the sensor in which synthetic data is generated by extracting and interpolating virtual FBG strains obtained from a large number of 3D finite-element calculations. A physical prototype is subsequently developed to demonstrate proof of concept. Results from laboratory validation tests give confidence in the sensor's ability to provide accurate contact stress measurements in typical soil–structure interface shear applications. In particular, the novel sensor structure and operating principle was shown to achieve excellent measurement of effective normal stress. The new sensor design harnesses many of the inherent benefits of FBGs including immunity to electromagnetic noise and water ingress, and the use a single lightweight cable and connector, which significantly simplifies installation on site compared to electrical multi-axis sensors.

Damage Analysis of Rolling Stock Automatic Coupler under Cyclic Loads

To ensure the operational safety of the automatic couplers of a rolling stock subjected to cyclic loading, reliability and residual service life must be determined. This paper involved the prediction of the durability of the automatic coupler SA3 by analytical summation of cyclic damage. During the investigation, the cyclic characteristics and damage modes of the automatic coupler material that influence the accumulation of damage are determined. The stress and strain state assessment model was developed using a 3D finite element method for the automatic coupler housing as a geometrically complex component. A methodology was used to assess damage per load cycle that is applicable to any sequence of automatic coupler load cycles. For this purpose, the authors used low-cycle stationary load dependencies that account for quasi-static and low-cycle fatigue damage. The investigation showed that a coupler may develop a crack due to accumulated quasi-static and fatigue damage. For damage summation, the dependences of low-cycle stationary stress-controlled load accounting for low-cycle, quasi-static fatigue damage are proposed in view of the variation of the load on the automatic coupler during operation depending on the weight of the rolling stock, velocity, and railway relief. The proposed methodology is applicable to the calculation of other housings under similar loading conditions.

Spur gear tooth root stress analysis by a 3D flexible multibody approach and a full-FE contact-based formulation

This paper proposes an original method to determine the gear tooth root stresses from a 3D finite element (FE) flexible multibody approach and a full-FE contact-based formulation. The contact problem is dealt with an augmented Lagrangian formulation whereas the analysis is performed by a preconditioned gradient solver (PCG). Tooth flank modifications are directly introduced within the 3D model. This one is thus able to take into account straightforwardly tooth bending and Hertzian-like deformations as well as the micro-geometry effect. Simulations are performed for several mesh periods, without making any assumptions about load distribution, tooth and gear blank flexibilities, and possible premature or delayed contacts between tooth pairs in quasi-static conditions. A precise distribution of tooth root stresses associated with instantaneous contacts conditions is then computed. For this study, a single stage spur gear with micro-geometry modifications corresponding to an arc-shaped profile crowning is modeled. Several output torques are considered. The obtained results are compared to those obtained using a 2D FE ISO-based model, where external forces are applied along the theoretical line of action.

Lifing Assessment of Gas Turbine Blade Root Affected by Out-of-Tolerances.

Current and future heavy-duty gas turbines (GTs) are being developed as an alternative or support to renewable energy sources (RESs). Therefore, GTs are subjected to several instances of being switched on and off; thus, the material fatigue limit can be reached in a short time. In such a scenario, possible out-of-tolerances (OoTs) in critical components must be considered. In this paper, OoTs related to critical parameters in the attachment geometry of a rotor blade are considered to estimate their impact on component life through a 2D finite element (FE) analysis. First, a mesh refinement is performed to obtain mesh-independent results; second, OoT geometries are simulated to determine stresses and strains at the blade attachment and disc groove. The mesh refinement process is critical to ensure model accuracy for both nominal and OoT geometries. The results show that OoTs can lead to important reductions in the life of the intended component, both on the blade and disc sides. These results could be useful in updating the maintenance plan for components and could be used for future insights, with further work extending the study to 3D geometry, for example, and evaluating the effect of other geometries.

Development of acoustic impedance tube system using 3D printing and finite element analysis

This study aims to demonstrate and provide a practical guide on correctly designing, fabricating, and verifying the acoustic impedance tube system. Many parts required to build this tube system were uniquely designed and 3D printed, such as the driver case, specimen holder, microphone holder, etc. In addition, this acoustic tube is designed to meet three functions: standard sound absorption measurement, in-situ measurement, and microphone mismatch measurement. Various types of specimens are tested to verify the accuracy of the measurement and how it compares to the poroelastic theoretical models and FEM models. To estimate the poroelastic parameters, measurement and FEM models were used to run the optimization process. The validation study shows good agreement in all three cases: measurement, the poroelastic model, and the FEM model.

Numerical investigation and seismic performance evaluation of an innovative prefabricated structural system

This paper focuses on an innovative structural system with separated gravity and lateral resisting systems (SGLR system), which is suitable for multi-story prefabricated buildings and demonstrates different seismic performance from the traditional rigid frame system (RF system). Elaborate three-dimensional (3D) finite element models for RF and SGLR systems are established and verified by test results. Based on the numerical models, the load-displacement relationship, force mechanism, energy dissipation capacity, local buckling, and stress and strain distribution of SGLR system are investigated in detail in comparison with RF system. In order to conduct the seismic performance evaluation, four groups of RF and SGLR substructures are designed with different lateral stiffness. Pushover analysis and response spectrum analysis are carried out, and safety factors are calculated on the basis of the incremental N2 method. It is indicated that SGLR system has lower loading capacity, smaller safety margin but larger deformation capability than RF system with the same lateral stiffness. Considering that the inter-story drift ratio is a common design index, more stringent design requirements should be proposed for SGLR system to obtain seismic performance equivalent to that of RF system.

Assessment of damage prediction models for composite laminates under single and repeated low-velocity impacts

The aim is to explore the influence of damage initiation criteria and evolution methods on the damage prediction for composite laminates under single and repeated low-velocity impacts. A 3D finite element model is established to analyze the impact behavior of laminate, and a damage analysis process including the initiation determination, progressive evolution, and constitutive relationship is designed. Hashin criterion based on tensile strain for matrix is modified by an empirical assumption method. Besides, the quantitative evolution of damage area during impact is studied, which provides a new perspective for elucidating the damage mechanism. The numerical model is first validated against the available experimental data. The influence of initiation criteria and evolution methods on damage prediction is then performed in two-phase research. The first phase discusses the prediction capabilities of different combinations under single impact. In the second phase, a method for removing residual vibrations under repeated impacts is proposed. The differences caused by the removal and non-removal of residual vibrations are analyzed, and the sensitivity of these differences to various combinations is discussed. Results show that the numerical prediction agrees well with the experiment in dynamic mechanical response curve. Combining Hashin-Strain criterion with linear equivalent strain method and Puck criterion with exponential equivalent displacement method yields optimal results. Moreover, the sensitivity to differences caused by the residual vibrations effects is related to the initiation criterion, and the evolution method.

Coupling model between building post structure and floors with the spectral element method

With the densification of urban area and public transport, the structure-borne noise disturbance due to railway traffic is a major issue. In order to mitigate it, ground-borne vibration propagation in building has to be investigated. The case of post-beam constructions necessitates 3D modeling with choices between complex, precise and computationally time-consuming approaches (3D finite element method for example), and fast but less accurate approaches like the statistical energy analysis (SEA) which need a strong modal density. For such post-beam constructions, the spectral element method (SEM) appears to be a good compromise. This method allows reducing the number of elements to represent the truss structure compared to the finite element method. However, the building walls and floor are difficult to consider with the SEM. In this paper, the challenge of coupling post-beam structures and floors is examined, a method is proposed and results are discussed. The aim is to develop a method applicable to different buildings while keeping advantage of SEM low discretization.