Probabilistic-Ellipsoid Hybrid Reliability Multi-Material Topology Optimization Method Based on Stress Constraint

This paper proposes a multi-material topology optimization method based on the hybrid reliability of the probability-ellipsoid model with stress constraint for the stochastic uncertainty and epistemic uncertainty of mechanical loads in optimization design.The probabilistic model is combined with the ellipsoidal model to describe the uncertainty of mechanical loads.The topology optimization formula is combined with the ordered solid isotropic material with penalization(ordered-SIMP)multi-material interpolation model.The stresses of all elements are integrated into a global stress measurement that approximates the maximum stress using the normalized p-norm function.Furthermore,the sequential optimization and reliability assessment(SORA)is applied to transform the original uncertainty optimization problem into an equivalent deterministic topology optimization(DTO)problem.Stochastic response surface and sparse grid technique are combined with SORA to get accurate information on the most probable failure point(MPP).In each cycle,the equivalent topology optimization formula is updated according to the MPP information obtained in the previous cycle.The adjoint variable method is used for deriving the sensitivity of the stress constraint and the moving asymptote method(MMA)is used to update design variables.Finally,the validity and feasibility of the method are verified by the numerical example of L-shape beam design,T-shape structure design,steering knuckle,and 3D T-shaped beam.

Stress tensor determination by modified hydraulic tests on pre-existing fractures:Method and stress constraints

The hydraulic testing of pre-existing fractures(HTPF)is one of the most promising in situ stress measurement methods,particularly for three-dimensional stress tensor determination.However,the stress tensor determination based on the HTPF method requires at least six tests or a minimum of 14-15 tests(under different conditions)for reliable results.In this study,we modified the HTPF method by considering the shear stress on each pre-existing fracture,which increased the number of equations for the stress tensor determination and decreased the number of tests required.Different shear stresses were attributed to different fractures by random sampling;therefore,the stress tensors were obtained by searching for the optimal solution using the least squares criterion based on the Monte Carlo method.Thereafter,we constrained the stress tensor based on the tensile strength criterion,compressive strength criterion,and vertical stress constraints.The inverted stress tensors were presented and analyzed based on the tensorial nature of the stress using the Euclidean mean stress tensor.Two stress-measurement campaigns in Weifang(Shandong Province,China)and Mercantour road tunnel(France)were implemented to highlight the validity and efficiency of the modified HTPF(M-HTPF)method.The results showed that the M-HTPF method can be applied for stress tensor inversion using only three to four tests on pre-existing fractures,neglecting the stress gradient.The inversion results were confined to relatively small distribution dispersions and were significantly reliable and stable due to the shear stresses on the fractures and the stress constraints employed.The M-HTPF method is highly feasible and efficient for complete stress tensor determination in a single borehole.

Dys-psychological Stress Effect on Expressions of P53 and NFκBp65 in Human Ovarian Carcinoma In Vivo

Objective： To investigate the dys-psychological stress effect on the growth of subcutaneous xenotransplanted tumor in nude mice bearing human epithelium ovarian carcinoma, and the influence on P53 and NFκBp65 expressions. Methods： The subcutaneous tumor xenografts were established by implanting human epithelium ovarian carcinoma tissues into nude mice and the dys-psychological stress model was established with restraint. The mice were randomized into the following four treatment groups with each group six mice respectively： tumor group （group A）, normal saline intraperitoneal injection; tumor with stress group （group B）, normal saline intraperitoneal injection; tumor therapy group （group C）, cisplatin intraperitoneal injection; and tumor therapy with stress group （group D）, cisplatin intraperitoneal injection. The expressions of P53 and NFκBp65 in tumor tissues were determined by Western blotting. Results： The expressions of P53 and NFκBp65 in each restraint group were enhanced compared with the control groups （P0.05）. Conclusion： The dys-psychological stress may induce the high expressions of P53 and NFκBp65 proteins and further promote tumor growth.

APPROACH FOR LAYOUT OPTIMIZATION OF TRUSS STRUCTURES WITH DISCRETE VARIABLES UNDER DYNAMIC STRESS, DISPLACEMENT AND STABILITY CONSTRAINTS

A mathematical model was developed for layout optimization of truss structures with discrete variables subjected to dynamic stress, dynamic displacement and dynamic stability constraints. By using the quasi-static method, the mathematical model of structure optimization under dynamic stress, dynamic displacement and dynamic stability constraints were transformed into one subjected to static stress, displacement and stability constraints. The optimization procedures include two levels, i.e., the topology optimization and the shape optimization. In each level, the comprehensive algorithm was used and the relative difference quotients of two kinds of variables were used to search the optimum solution. A comparison between the optimum results of model with stability constraints and the optimum results of model without stability constraint was given. And that shows the stability constraints have a great effect on the optimum solutions.

A METHOD FOR TOPOLOGICAL OPTIMIZATION OF STRUCTURES WITH DISCRETE VARIABLES UNDER DYNAMIC STRESS AND DISPLACEMENT CONSTRAINTS

A method for topological optimization of structures with discrete variables subjected to dynamic stress and displacement constraints is presented. By using the quasistatic method, the structure optimization problem under dynamic stress and displacement constraints is converted into one subjected to static stress and displacement constraints. The comprehensive algorithm for topological optimization of structures with discrete variables is used to find the optimum solution.

THE TOPOLOGICAL OPTIMIZATION FOR TRUSS STRUCTURES WITH STRESS CONSTRAINTS BASED ON THE EXIST-NULL COMBINED MODEL

A new exist-null combined model is proposed for the structural topology optimization. The model is applied to the topology optimization of the truss with stress constraints. Satisfactory computational result can be obtained with more rapid and more stable convergence as compared with the cross-sectional optimization. This work also shows that the presence of independent and continuous topological variable motivates the research of structural topology optimization.

Dynamic response characteristics of dry and watersaturated schist under impact loading

Many geological engineering hazards are closely related to the dynamic mechanical properties of rock materials.However,most existing studies on the dynamic mechanical properties of rock materials were conducted on the hard rocks such as sandstone,granite,limestone,and marble,whereas soft rocks,such as schist,are less studied.Therefore,in this study,a series of triaxial impact tests were conducted on dry and saturated schist by employing a modified triaxial split Hopkinson pressure bar system to reveal the coupling effects of water,strain rate,and triaxial confining pressure on the mechanical properties of schist.The results show that schist is a type of watersensitive rock and the stress-strain curve of saturated schist has apparent ductility.The effects of strain rate on dynamic strain,deformation modulus and peak stress were analyzed.The results also show that the dynamic peak stress is affected by the combined softening effect and viscous effect of water under impact loading.Finally,it was found that the failure mode of schist belongs to typical axial tensile failure under uniaxial impact tests,and shear failure is the main failure mode under triaxial impact tests.With the increase in confining pressure,the failure modes of schist change from tensile failure to shear failure.This research can provide useful parameters for geological engineering hazard prevention in mountain areas.

Topology Optimization with Aperiodic Load Fatigue Constraints Based on Bidirectional Evolutionary Structural Optimization

Because of descriptive nonlinearity and computational inefficiency,topology optimization with fatigue life under aperiodic loads has developed slowly.A fatigue constraint topology optimization method based on bidirectional evolutionary structural optimization(BESO)under an aperiodic load is proposed in this paper.In viewof the severe nonlinearity of fatigue damagewith respect to design variables,effective stress cycles are extracted through transient dynamic analysis.Based on the Miner cumulative damage theory and life requirements,a fatigue constraint is first quantified and then transformed into a stress problem.Then,a normalized termination criterion is proposed by approximatemaximum stress measured by global stress using a P-normaggregation function.Finally,optimization examples show that the proposed algorithm can not only meet the requirements of fatigue life but also obtain a reasonable configuration.

FORMING MICHELL TRUSS IN THREE-DIMENSIONS BY FINITE ELEMENT METHOD

The finite element method to form Michell truss in three_dimensions is presented.The orthotropic composite with fiber_reinforcement is employed as the material model to simulate Michell truss.The orientation and densities of fibers at nodes are taken as basic design variables.The stresses and strains at nodes are calculated by finite element method.An iteration scheme is suggested to adjust the orientations of fibers to be along the orientations of principal stresses, and the densities of fibers according to the strains in the orientations of fibers.The strain field satisfying Michell criteria and truss_like continuum are achieved after several iterations. Lastly, the Michell truss is showed by continuous lines, which are formed according to the orientations of fibers at nodes. Several examples are used to demonstrate the efficiency of the presented approach.

ALTERNATING OPTIMIZATION METHOD FOR ISOGEOMETRIC TOPOLOGY OPTIMIZATION WITH STRESS CONSTRAINTS

Topology optimization (TO) has developed rapidly recently. However, topology optimization with stress constraints still faces many challenges due to its highly non-linear properties which will cause inefficient computation, iterative oscillation, and convergence guarantee problems. At the same time, isogeometric analysis (IGA) is accepted by more and more researchers, and it has become one important tool in the field of topology optimization because of its high fidelity. In this paper, we focus on topology optimization with stress constraints based on isogeometric analysis to improve computation efficiency and stability. A new hybrid solver combining the alternating direction method of multipliers and the method of moving asymptotes (ADMM-MMA) is proposed to solve this problem. We first generate an initial feasible point by alternating direction method of multipliers (ADMM) in virtue of the rapid initial descent property. After that, we adopt the method of moving asymptotes (MMA) to get the final results. Several benchmark examples are used to verify the proposed method, and the results show its feasibility and effectiveness.

Effect of three dimensional stress state on unstable fracture condition and crack opening level in a new crack growth model

A new model for the analysis of fatigue crack growth in the metal structures was proposed. This model shows a promising capability of explaining various fatigue phenomena. The new crack growth model is further completed by a continuous empirical formula for estimating the value of variable fracture toughness during crack propagation and a modified continuous equation for the crack tip stress/strain constraint factor used to calculate the stress intensity factor at the opening level. The prediction results are proved to agree well with the observed phenomena in test.

Topology optimization of turbine disk considering maximum stress prediction and constraints

For the stress-constrained topology optimization of a turbine disk under centrifugal loads,the jagged boundaries of the mesh and the gray densities on the solid/void interfaces could make the calculated stress field inconsistent with the actual value.It may result in overestimating the maximum stress and thus affect the effectiveness of stress constraints.This paper proposes a new method for predicting the maximum stress to overcome the difficulty.In the process,a predicted density is newly defined to obtain stable boundaries with thin layers of gray elements,a transition factor is innovatively proposed to evaluate the effects of intermediate-density elements,two different stiffness penalty schemes are flexibly used to calculate the elastic modulus of elements,and a linear stress penalty is further adopted to relax the stress field of the structure.The proposed approach for predicting the maximum stress value is verified by the analysis of a structure with smooth boundaries and the topology optimization of a turbine disk.An updating scheme of the stress constraint in the topology optimization is also developed using the predicted maximum stress.Some key ingredients affecting the optimization results are discussed in detail.The results prove the effectiveness and efficacy of the proposed maximum stress prediction and developed stress constraint methods.

Topology Optimization of Truss-Like Structure with Stress Constraints Under Multiple-Load Cases

A new method for topology optimization of truss-like structures with stress constraints under multiple-load cases(MLCs)is presented.A spatial truss-like material model with three families of orthotropic members is adopted,in which the three families of members along three orthotropic directions are embedded continuously in a weak matrix.The densities and directions of the three families of members at the nodes are taken as the design variables.An optimality criterion is suggested based on the concept of directional stiffness.First,under each single-load case(SLC),the truss-like structure is optimized as per the fully stressed criterion.Accordingly,the directional stiffness of the optimal structure under an SLC at every node is obtained.Next,the directional stiffness of the truss-like structure under MLCs is determined by ensuring that the directional stiffness is as similar as possible to the maximum directional stiffness of the optimal structure under every SLC along all directions.Finally,the directions and densities of the members in the optimal truss-like structures under MLCs are obtained by solving the eigenvalue problems of the coefficient matrix of the directional stiffness at every node.Two examples are presented to demonstrate the effectiveness and efficiency of the method.

ELASTIC CALCULATION OF STRESSES IN RINGS USING AIRY STRESS FUNCTION

Stress calculation formulae for a ring have been obtained by using Airy stress function of the plane strain field with the decomposition of the solutions for normal stresses of Airy biharmonic equation into two parts when it is loaded under two opposite inside forces along a diameter. One part should fulfill a constraint condition about normal stress distribution along the circumference at an energy valley to do the minimum work. Other part is a stress residue constant. In order to verify these formulae and the computed results, the computed contour lines of equi-maximal shear stresses were plotted and quite compared with that of photo-elasticity test results. This constraint condition about normal stress distribution along circumference is confirmed by using Greens’ theorem. An additional compression exists along the circumference of the loaded ring, explaining the divorcement and displacement of singularity points at inner and outer boundaries.

Geometrically Nonlinear Topology Optimization of Continuum Structures Based on an Independent Continuous Mapping Method

A geometrically nonlinear topology optimization method for continuum structures is proposed based on the independent continuous mapping method.The stress constraint problem is studied due to the importance of structural strength in engineering applications.First,a topology optimization model is established for a lightweight structure with element stress as constraints.Second,the stress globalization method is adopted to convert local stress constraints into strain energy constraints,which overcomes the difficulties caused by local stress constraints,such as model establishment,sensitivity analysis,and massive solution calculations.Third,the sensitivity of the objective function and constraint function is analyzed,and the method of moving asymptotes is employed to solve the optimization model.In addition,the additive hyperelasticity technique is utilized to solve the numerical instability induced by structures undergoing large deformation.Numerical examples are given to validate the feasibility of the proposed method.The method provides a significant reference for geometrically nonlinear optimization design.

Topology Optimization of Geometrically Nonlinear Structures Under Thermal-Mechanical Coupling

A geometrically nonlinear topology optimization(GNTO)method with thermal–mechanical coupling is investigated.Firstly,the new expression of element coupling stress due to superimposed mechanical and thermal loading is obtained based on the geometrically nonlinear finite element analysis.The lightweight topology optimization(TO)model under stress constraints is established to satisfy the strength requirement.Secondly,the distortion energy theory is introduced to transform themodel into structural strain energy constraints in order to solve the implicit relationship between stress constraints and design variables.Thirdly,the sensitivity analysis of the optimization model is derived,and the model is solved by the method of moving asymptotes(MMA).Numerical examples show that temperature has a significant effect on the optimal configuration,and the TO method considering temperature load is closer to engineering design requirements.The proposed method can be extended to the GNTO design with multiple physical field coupling.