Physical and numerical modelling of the inherent variability of shear strength in soil mechanics
-
- Geomechanics and Engineering
- /
- v.17 no.1
- /
- pp.31-45
- /
- 2019
In this study the spatial variability of soils is substantiated physically and numerically by using random field theory. Heterogeneous samples are fabricated by combining nine homogeneous soil clusters that are assumed to be elements of an adopted random field. Homogeneous soils are prepared by mixing different percentages of kaolin and bentonite at water contents equivalent to their respective liquid limits. Comprehensive characteristic laboratory tests were carried out before embarking on direct shear experiments to deduce the basic correlations and properties of nine homogeneous soil clusters that serve to reconstitute the heterogeneous samples. The tests consist of Atterberg limits, and Oedometric and unconfined compression tests. The undrained shear strength of nine soil clusters were measured by the unconfined compression test data, and then correlations were made between the water content and the strength and stiffness of soil samples with different consistency limits. The direct shear strength of heterogeneous samples of different stochastic properties was then evaluated by physical and numerical modelling using FISH code programming in finite difference software of
When a building suffers damages under moderate to severe loading condition, its physical properties such as damping and stiffness parameters will change. There are different practical methods besides various numerical procedures that have successfully detected a range of these changes. Almost all the previous proposed methods used to work with translational components of mode shapes, probably because extracting these components is more common in vibrational tests. This study set out to investigate the influence of using both rotational and translational components of mode shapes, in detecting damages in 3-D steel structures elements. Three different sets of measured components of mode shapes are examined: translational, rotational, and also rotational/translational components in all joints. In order to validate our assumptions two different steel frames with three damage scenarios are considered. An iterative model updating program is developed in the MATLAB software that uses the OpenSees as its finite element analysis engine. Extensive analysis shows that employing rotational components results in more precise prediction of damage location and its intensity. Since measuring rotational components of mode shapes still is not very convenient, modal dynamic expansion technique is applied to generate rotational components from measured translational ones. The findings indicated that the developed model updating program is really efficient in damage detection even with generated data and considering noise effects. Moreover, methods which use rotational components of mode shapes can predict damage's location and its intensity more precisely than the ones which only work with translational data.
In this paper, semi-rigid elasto-plastic post-buckling analysis of a space frame was performed using various explicit arc-length methods. Various explicit arc-length methodsand a large-deformation and small-strain elasto-plastic 3D space frame element with semi-rigid connections and plastic hinges were developed. This element can be appliedto both explicit and implicit numerical algorithms. In this study, the Dynamic Relaxation method was adopted in the predictor and corrector processesto formulate an explicit arc-length algorithm. The developed "explicit-predictor" or "explicit-corrector" were used in the elasto-plastic post-buckling analysis. The Eulerian equations for a beam-column with finite rotation, which considers the bowing effects, were adopted for the elastic system and extended to theinelastic system with a plastic hinge concept. The derived tangent stiffness matrix was asymmetrical due to the finite rotation. The joint connection elements were introduced for semi-rigidity using a static condensation technique. Semi-rigid elasto-plastic post-buckling analyses were carried out to demonstrate the potential of the developed explicit arc-length method and advanced space frame element in terms of accuracy and efficiency.
Oral implants must fulfill certain criteria arising from special demands of function, which include biocompatibility, adequate mechanical strength, optimum soft and hard tissue integration, and transmission of functional forces to bone within physiological limits. And one of the critical elements influencing the long-term uncompromise functioning of oral implants is load distribution at the implant- bone interface, Factors that affect the load transfer at the bone-implant interface include the type of loading, material properties of the implant and prosthesis, implant geometry, surface structure, quality and quantity of the surrounding bone, and nature of the bone-implant interface. To understand the biomechanical behavior of dental implants, validation of stress and strain measurements is required. The finite element analysis (FEA) has been applied to the dental implant field to predict stress distribution patterns in the implant-bone interface by comparison of various implant designs. This method offers the advantage of solving complex structural problems by dividing them into smaller and simpler interrelated sections by using mathematical techniques. The purpose of this study was to evaluate the stresses induced around the implants in bone using FEA, A 3D FEA computer software (SOLIDWORKS 2004, DASSO SYSTEM, France) was used for the analysis of clinical simulations. Two types (external and internal) of implants of 4.1 mm diameter, 12.0 mm length were buried in 4 types of bone modeled. Vertical and oblique forces of lOON were applied on the center of the abutment, and the values of von Mises equivalent stress at the implant-bone interface were computed. The results showed that von Mises stresses at the marginal. bone were higher under oblique load than under vertical load, and the stresses were higher at the lingual marginal bone than at the buccal marginal bone under oblique load. Under vertical and oblique load, the stress in type I, II, III bone was found to be the highest at the marginal bone and the lowest at the bone around apical portions of implant. Higher stresses occurred at the top of the crestal region and lower stresses occurred near the tip of the implant with greater thickness of the cortical shell while high stresses surrounded the fixture apex for type N. The stresses in the crestal region were higher in Model 2 than in Model 1, the stresses near the tip of the implant were higher in Model 1 than Model 2, and Model 2 showed more effective stress distribution than Model.
This study investigates the biomechanical efficacies of various cement augmentation techniques with or without pressurization for varying degrees of osteoporotic femur. For this study, a biomechanical analysis using a finite element method (FEM) was undertaken to evaluate surgical procedures, Simulated models include the non-cemented(i.e., hip screw only, Type I), the cement-augmented(Type II), and the cemented augmented with pressurization(Type III) models. To simulate the fracture plane and other interfacial regions, 3-D contact elements were used with appropriate friction coefficients. Material properties of the cancellous bone were varied to accommodate varying degrees of osteoporosis(Singh indices, II∼V). For each model. the following items were analyzed to investigate the effect surgical procedures in relation to osteoporosis of varying degrees : (a) von Mises stress distribution within the femoral head in terms of volumetric percentages. (b) Peak von Mises stress(PVMS) within the femoral head and the surgical constructs. (c) Maximum von Mises strain(MVMS) within the femoral head, (d) micromotions at the fracture plane and at the interfacial region between surgical construct and surrounding bone. Type III showed the lowest PVMS and MVMS at the cancellous bone near the bone-construct interface regardless of bone densities. an indication of its least likelihood of construct loosening due to failure of the host bone. Particularly, its efficacy was more prominent when the bone density level was low. Micromotions at the interfacial surgical construct was lowest in Type III. followed by Type I and Type II. They were about 15-20% of other types. which suggested that pressurization was most effective in limiting the interfacial motion. Our results demonstrated the cement augmentation with hip screw could be more effective when used with pressurization technique for the treatment of intertrochanteric fractures. For patients with low bone density. its effectiveness can be more pronounced in limiting construct loosening and promoting bone union.
Fuzzy logic based Control Theory has gained much interest in the industrial world, thanks to its ability to formalize and solve in a very natural way many problems that are very difficult to quantify at an analytical level. This paper shows a solution for treating membership function inside hardware circuits. The proposed hardware structure optimizes the memoried size by using particular form of the vectorial representation. The process of memorizing fuzzy sets, i.e. their membership function, has always been one of the more problematic issues for the hardware implementation, due to the quite large memory space that is needed. To simplify such an implementation, it is commonly [1,2,8,9,10,11] used to limit the membership functions either to those having triangular or trapezoidal shape, or pre-definite shape. These kinds of functions are able to cover a large spectrum of applications with a limited usage of memory, since they can be memorized by specifying very few parameters ( ight, base, critical points, etc.). This however results in a loss of computational power due to computation on the medium points. A solution to this problem is obtained by discretizing the universe of discourse U, i.e. by fixing a finite number of points and memorizing the value of the membership functions on such points [3,10,14,15]. Such a solution provides a satisfying computational speed, a very high precision of definitions and gives the users the opportunity to choose membership functions of any shape. However, a significant memory waste can as well be registered. It is indeed possible that for each of the given fuzzy sets many elements of the universe of discourse have a membership value equal to zero. It has also been noticed that almost in all cases common points among fuzzy sets, i.e. points with non null membership values are very few. More specifically, in many applications, for each element u of U, there exists at most three fuzzy sets for which the membership value is ot null [3,5,6,7,12,13]. Our proposal is based on such hypotheses. Moreover, we use a technique that even though it does not restrict the shapes of membership functions, it reduces strongly the computational time for the membership values and optimizes the function memorization. In figure 1 it is represented a term set whose characteristics are common for fuzzy controllers and to which we will refer in the following. The above term set has a universe of discourse with 128 elements (so to have a good resolution), 8 fuzzy sets that describe the term set, 32 levels of discretization for the membership values. Clearly, the number of bits necessary for the given specifications are 5 for 32 truth levels, 3 for 8 membership functions and 7 for 128 levels of resolution. The memory depth is given by the dimension of the universe of the discourse (128 in our case) and it will be represented by the memory rows. The length of a world of memory is defined by: Length = nem (dm(m)+dm(fm) Where: fm is the maximum number of non null values in every element of the universe of the discourse, dm(m) is the dimension of the values of the membership function m, dm(fm) is the dimension of the word to represent the index of the highest membership function. In our case then Length=24. The memory dimension is therefore 128*24 bits. If we had chosen to memorize all values of the membership functions we would have needed to memorize on each memory row the membership value of each element. Fuzzy sets word dimension is 8*5 bits. Therefore, the dimension of the memory would have been 128*40 bits. Coherently with our hypothesis, in fig. 1 each element of universe of the discourse has a non null membership value on at most three fuzzy sets. Focusing on the elements 32,64,96 of the universe of discourse, they will be memorized as follows: The computation of the rule weights is done by comparing those bits that represent the index of the membership function, with the word of the program memor . The output bus of the Program Memory (μCOD), is given as input a comparator (Combinatory Net). If the index is equal to the bus value then one of the non null weight derives from the rule and it is produced as output, otherwise the output is zero (fig. 2). It is clear, that the memory dimension of the antecedent is in this way reduced since only non null values are memorized. Moreover, the time performance of the system is equivalent to the performance of a system using vectorial memorization of all weights. The dimensioning of the word is influenced by some parameters of the input variable. The most important parameter is the maximum number membership functions (nfm) having a non null value in each element of the universe of discourse. From our study in the field of fuzzy system, we see that typically nfm 3 and there are at most 16 membership function. At any rate, such a value can be increased up to the physical dimensional limit of the antecedent memory. A less important role n the optimization process of the word dimension is played by the number of membership functions defined for each linguistic term. The table below shows the request word dimension as a function of such parameters and compares our proposed method with the method of vectorial memorization[10]. Summing up, the characteristics of our method are: Users are not restricted to membership functions with specific shapes. The number of the fuzzy sets and the resolution of the vertical axis have a very small influence in increasing memory space. Weight computations are done by combinatorial network and therefore the time performance of the system is equivalent to the one of the vectorial method. The number of non null membership values on any element of the universe of discourse is limited. Such a constraint is usually non very restrictive since many controllers obtain a good precision with only three non null weights. The method here briefly described has been adopted by our group in the design of an optimized version of the coprocessor described in [10].
Statement of problem & Purpose: The purpose of this study was to investigate the effect of a mouth guard material properties on the skull and brain when they were under impact loads on mandible. Material and methods: Two customized mouth protectors having different material propeerst ieach other were made for a female Korean who had no history of brain trauma, no cerebral diseases, nomal occlusion and natural dentition. The 3D finite element model of human skull and brain scanned by means of computed tomography was constructed. The FEM model of head was composed of 407,825 elements and 82,138 nodes, including skull, brain, maxilla, mandible, articular disc, teeth and mouth guard. The stress concentrations on maxillary teeth, maxilla and skull with two mouth guards were evaluated under oblique impact load of 800N onto mandibular 3 loading points for 0.1sec. And the brain relative displacement was compared in two different mouth guard materials under same condition. Result and Conclusion: The results were as follows; 1. In comparison of von Mises stress on maxillary teeth, a soft mouth guard material had significantly lower stress values on measuring point than a hard mouth protector materials (P < .05). 2. In comparison of von Mises stress on maxilla and skull, A soft mouth protector material had significantly lower stress values on measuring point than a hard mouth protector materials (P < .05). 3. For impact loads on mandible, there were more stress concentrated area on maxilla and skull with hard mouth guard than soft with mouth protector. 4. For impact loads on mandible, brain relative displacement had little relation with mouth guard material properties. In results of this study, soft mouth guard materials were superior to hard mouth guard materials for mandible impact loads for prevention of sports injuries. Although the results of this study were not enough to figure out the roles of needed mouth guard material properties for a human head, we got some knowledge of the pattern about stress concentration and distribution on maxilla and skull for impact loads with soft or hard mouth protector. More studies are needed to substantiate the relationship between the mouth guard materials and sports injuries.
The wall shear stress in the vicinity of end-to end anastomoses under steady flow conditions was measured using a flush-mounted hot-film anemometer(FMHFA) probe. The experimental measurements were in good agreement with numerical results except in flow with low Reynolds numbers. The wall shear stress increased proximal to the anastomosis in flow from the Penrose tubing (simulating an artery) to the PTFE: graft. In flow from the PTFE graft to the Penrose tubing, low wall shear stress was observed distal to the anastomosis. Abnormal distributions of wall shear stress in the vicinity of the anastomosis, resulting from the compliance mismatch between the graft and the host artery, might be an important factor of ANFH formation and the graft failure. The present study suggests a correlation between regions of the low wall shear stress and the development of anastomotic neointimal fibrous hyperplasia(ANPH) in end-to-end anastomoses. 30523 T00401030523 ^x Air pressure decay(APD) rate and ultrafiltration rate(UFR) tests were performed on new and saline rinsed dialyzers as well as those roused in patients several times. C-DAK 4000 (Cordis Dow) and CF IS-11 (Baxter Travenol) reused dialyzers obtained from the dialysis clinic were used in the present study. The new dialyzers exhibited a relatively flat APD, whereas saline rinsed and reused dialyzers showed considerable amount of decay. C-DAH dialyzers had a larger APD(11.70
The wall shear stress in the vicinity of end-to end anastomoses under steady flow conditions was measured using a flush-mounted hot-film anemometer(FMHFA) probe. The experimental measurements were in good agreement with numerical results except in flow with low Reynolds numbers. The wall shear stress increased proximal to the anastomosis in flow from the Penrose tubing (simulating an artery) to the PTFE: graft. In flow from the PTFE graft to the Penrose tubing, low wall shear stress was observed distal to the anastomosis. Abnormal distributions of wall shear stress in the vicinity of the anastomosis, resulting from the compliance mismatch between the graft and the host artery, might be an important factor of ANFH formation and the graft failure. The present study suggests a correlation between regions of the low wall shear stress and the development of anastomotic neointimal fibrous hyperplasia(ANPH) in end-to-end anastomoses. 30523 T00401030523 ^x Air pressure decay(APD) rate and ultrafiltration rate(UFR) tests were performed on new and saline rinsed dialyzers as well as those roused in patients several times. C-DAK 4000 (Cordis Dow) and CF IS-11 (Baxter Travenol) reused dialyzers obtained from the dialysis clinic were used in the present study. The new dialyzers exhibited a relatively flat APD, whereas saline rinsed and reused dialyzers showed considerable amount of decay. C-DAH dialyzers had a larger APD(11.70