scholarly journals Topology optimization of quasistatic contact problems

Author(s):  
Andrzej Myśliński

Topology optimization of quasistatic contact problemsThis paper deals with the formulation of a necessary optimality condition for a topology optimization problem for an elastic contact problem with Tresca friction. In the paper a quasistatic contact model is considered, rather than a stationary one used in the literature. The functional approximating the normal contact stress is chosen as the shape functional. The aim of the topology optimization problem considered is to find the optimal material distribution inside a design domain occupied by the body in unilateral contact with the rigid foundation to obtain the optimally shaped domain for which the normal contact stress along the contact boundary is minimized. The volume of the body is assumed to be bounded. Using the material derivative and asymptotic expansion methods as well as the results concerning the differentiability of solutions to quasistatic variational inequalities, the topological derivative of the shape functional is calculated and a necessary optimality condition is formulated.

Author(s):  
Guang Dong ◽  
Zheng-Dong Ma ◽  
Gregory Hulbert ◽  
Noboru Kikuchi ◽  
Sudhakar Arepally ◽  
...  

Efficient and reliable sensitivity analyses are critical for topology optimization, especially for multibody dynamics systems, because of the large number of design variables and the complexities and expense in solving the state equations. This research addresses a general and efficient sensitivity analysis method for topology optimization with design objectives associated with time dependent dynamics responses of multibody dynamics systems that include nonlinear geometric effects associated with large translational and rotational motions. An iterative sensitivity analysis relation is proposed, based on typical finite difference methods for the differential algebraic equations (DAEs). These iterative equations can be simplified for specific cases to obtain more efficient sensitivity analysis methods. Since finite difference methods are general and widely used, the iterative sensitivity analysis is also applicable to various numerical solution approaches. The proposed sensitivity analysis method is demonstrated using a truss structure topology optimization problem with consideration of the dynamic response including large translational and rotational motions. The topology optimization problem of the general truss structure is formulated using the SIMP (Simply Isotropic Material with Penalization) assumption for the design variables associated with each truss member. It is shown that the proposed iterative steps sensitivity analysis method is both reliable and efficient.


Author(s):  
Guang Dong ◽  
Yulan Song

The topology optimization method is extended to solve a single phase flow in porous media optimization problem based on the Two Point Flux Approximation model. In particular, this paper discusses both strong form and matrix form equations for the flow in porous media. The design variables and design objective are well defined for this topology optimization problem, which is based on the Solid Isotropic Material with Penalization approach. The optimization problem is solved by the Generalized Sequential Approximate Optimization algorithm iteratively. To show the effectiveness of the topology optimization in solving the single phase flow in porous media, the examples of two-dimensional grid cell TPFA model with impermeable regions as constrains are presented in the numerical example section.


2015 ◽  
Vol 2015 ◽  
pp. 1-20 ◽  
Author(s):  
Rudolf Pernis ◽  
Tibor Kvackaj

The calculation of average material contact pressure to rolls base on mathematical theory of rolling process given by Karman equation was solved by many authors. The solutions reported by authors are used simplifications for solution of Karman equation. The simplifications are based on two cases for approximation of the circular arch: (a) by polygonal curve and (b) by parabola. The contribution of the present paper for solution of two-dimensional differential equation of rolling is based on description of the circular arch by equation of a circle. The new term relative stress as nondimensional variable was defined. The result from derived mathematical models can be calculated following variables: normal contact stress distribution, front and back tensions, angle of neutral point, coefficient of the arm of rolling force, rolling force, and rolling torque during rolling process. Laboratory cold rolled experiment of CuZn30 brass material was performed. Work hardening during brass processing was calculated. Comparison of theoretical values of normal contact stress with values of normal contact stress obtained from cold rolling experiment was performed. The calculations were not concluded with roll flattening.


Author(s):  
Hae Chang Gea ◽  
Xing Liu ◽  
Euihark Lee ◽  
Limei Xu

In this paper, topology optimization under multiple independent loadings with uncertainty is presented. In engineering practice, load uncertainty can be found in many applications. From the literature, researchers have focused mainly on problems containing only a single uncertain external load. However, such idealistic problems may not be very useful in engineering practice. Problems involving multi-loadings with uncertainty are more commonly found in engineering applications. This paper presents a method to solve a system which contains multiple independent loadings with load uncertainty. First, a two-level optimization problem is formulated. The upper level problem is a typical topology optimization problem to minimize the mean compliance in the design using the worst case conditions. The lower level optimization problem is to solve for the worst loadings corresponding to the critical structure response. At the lower level formulation, an unknown-but-bounded model is used to define uncertain loadings. There are two challenges in finding the worst loading case: non-convexity and multi-loadings. The non-convexity problem is addressed by reformulating the problem as an inhomogeneous eigenvalue problem by applying the KKT optimality conditions and the multi-uncertain loadings problem is solved by an iterative method. After the worst loadings are generated, the upper level problem can be solved by a general topology optimization method. The effectiveness of the proposed method is demonstrated by numerical examples.


Author(s):  
Behzad Majdi ◽  
Arash Reza

The present study aims at providing a topology optimization of multi-material compliant mechanisms using solid isotropic material with penalization (SIMP) approach. In this respect, three multi-material gripper, invertor, and cruncher compliant mechanisms are considered that consist of three solid phases, including polyamide, polyethylene terephthalate, and polypropylene. The alternating active-phase algorithm is employed to find the distribution of the materials in the mechanism. In this case, the multiphase topology optimization problem is divided into a series of binary phase topology optimization sub-problems to be solved partially in a sequential manner. Finally, the maximum displacement of the multi-material compliant mechanisms was validated against the results obtained from the finite element simulations by the ANSYS Workbench software, and a close agreement between the results was observed. The results reveal the capability of the SIMP method to accurately conduct the topology optimization of multi-material compliant mechanisms.


Author(s):  
Qi Xia ◽  
Tielin Shi

Load-normalized strain energy increments between consecutive load steps are aggregated through the Kreisselmeier-Steinhauser (KS) function, and the KS function is proposed as a stiffness criterion of geometrically nonlinear structures. A topology optimization problem is defined to minimize the KS function together with the perimeter of structure and a volume constraint. The finite element analysis is done by remeshing, and artificial weak material is not used. The topology optimization problem is solved by using the level set method. Several numerical examples in two dimensions are provided. Other criteria of stiffness, i.e., the end compliance and the complementary work, are compared.


Sign in / Sign up

Export Citation Format

Share Document