[go: up one dir, main page]

CN106777768B - Optimal design method for eliminating tensile wrinkles of film structure - Google Patents

Optimal design method for eliminating tensile wrinkles of film structure Download PDF

Info

Publication number
CN106777768B
CN106777768B CN201710011463.1A CN201710011463A CN106777768B CN 106777768 B CN106777768 B CN 106777768B CN 201710011463 A CN201710011463 A CN 201710011463A CN 106777768 B CN106777768 B CN 106777768B
Authority
CN
China
Prior art keywords
film structure
design
constraint
optimization
thin film
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201710011463.1A
Other languages
Chinese (zh)
Other versions
CN106777768A (en
Inventor
罗阳军
亢战
李明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN201710011463.1A priority Critical patent/CN106777768B/en
Publication of CN106777768A publication Critical patent/CN106777768A/en
Application granted granted Critical
Publication of CN106777768B publication Critical patent/CN106777768B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Moulding By Coating Moulds (AREA)
  • Separation Using Semi-Permeable Membranes (AREA)

Abstract

The invention discloses an optimization design method for eliminating tensile wrinkles of a film structure, and solves the problem that a macroscopic film structure and a graphene structure are easy to wrinkle under the action of tensile load. On the basis of plane finite element analysis, the minimum principal stress of each unit is restrained to be a positive value, the principal stress distribution of the film structure is regulated, the film material distribution is designed by adopting a topological optimization technology, and then the structural form with a curve boundary or a hole is obtained, so that the aim of completely eliminating wrinkles is fulfilled. The mode ensures that wrinkles are completely eliminated, the accurate positioning of the appearance and the holes of the complex film structure can be realized, the automation degree of the optimized design is high, and the research and development efficiency of the wrinkle-free innovative design of the film structure can be ensured.

Description

一种用于消除薄膜结构拉伸褶皱的优化设计方法An Optimal Design Method for Eliminating Tensile Wrinkles in Film Structures

技术领域technical field

本发明属于航空航天薄膜结构设计和石墨烯纳米材料结构设计技术领域,提出了一种以消除褶皱现象为目的的薄膜结构优化设计方法。本发明采用一种新的优化模型及方法对薄膜边界进行裁剪和对内部进行孔洞设计,实现拉伸载荷下薄膜结构主应力分布的调控,从而有效消除薄膜在拉伸状态下可能产生的褶皱,完全实现薄膜类结构的无皱化。The invention belongs to the technical field of aerospace thin film structure design and graphene nanomaterial structure design, and proposes a thin film structure optimization design method aiming at eliminating wrinkle phenomenon. The invention adopts a new optimization model and method to cut the boundary of the film and design the holes in the interior, so as to realize the regulation of the principal stress distribution of the film structure under the tensile load, so as to effectively eliminate the possible wrinkles of the film in the stretched state. Completely achieve wrinkle-free film-like structures.

背景技术Background technique

随着科技的发展和人类文明的进步,作为一种典型的结构型式,空间薄膜结构越来越多地被用于航空航天结构,如太阳帆、充气天线、薄膜反射镜等。这些结构充分利用了薄膜易折叠\易展开、重量轻、体积小等优势,可以解决火箭运载对体积、质量的限制与不断增加的大尺寸、大口径使用要求之间的矛盾,因此极具诱人的应用前景。但是,由于弯曲刚度非常小,薄膜即使在张拉状态下也很容易发生面外屈曲现象,即褶皱。在航空航天应用中,经常要求这类薄膜必须保持表面光滑。例如薄膜反射镜,固定边界条件容易导致薄膜产生褶皱,进而影响表面光的反射从而降低成像的精确性。太阳帆的四个角受到集中力也易产生褶皱,从而影响光子反射角度及太阳光子压力方向。并且,大的褶皱还可能导致光子能量的局部集中引起局部高温,发生蠕变现象影响薄膜寿命。因此,如何有效消除薄膜褶皱现象在航空航天领域显得尤其重要。With the development of science and technology and the progress of human civilization, as a typical structural type, space thin-film structures are increasingly used in aerospace structures, such as solar sails, inflatable antennas, and thin-film mirrors. These structures make full use of the advantages of easy folding and unfolding of films, light weight and small size, which can solve the contradiction between the limitation of rocket carrying on volume and quality and the increasing requirements for large size and large caliber use, so they are very attractive. human application prospects. However, due to the very small bending stiffness, the film is prone to out-of-plane buckling, or wrinkling, even in tension. In aerospace applications, such films are often required to maintain a smooth surface. For example, thin-film mirrors, fixed boundary conditions can easily lead to wrinkles in the thin film, which in turn affects the reflection of surface light and reduces the accuracy of imaging. The four corners of the solar sail are also prone to wrinkles due to the concentrated force, which affects the photon reflection angle and the direction of the solar photon pressure. In addition, large wrinkles may also lead to local concentration of photon energy and cause local high temperature, and the creep phenomenon will affect the life of the film. Therefore, how to effectively eliminate the phenomenon of film wrinkles is particularly important in the aerospace field.

另外,在纳米材料领域,作为目前发现的最薄、强度最大、导电导热性能最强的一种新型纳米材料,石墨烯被科学家称为可能“彻底改变21世纪”的“新材料之王”,在移动设备、航空航天、新能源电池领域极具发展潜力。石墨烯(Graphene)是一种只有一个原子层厚度的准二维材料结构,厚度大约为0.335nm,具有非常类似于平面薄膜的力学特性,在拉伸荷载作用下也会产生面外褶皱现象,将影响其优良的电学、力学等性能。因此,在很多应用中有必要消除石墨烯结构中的褶皱现象。In addition, in the field of nanomaterials, as a new type of nanomaterial with the thinnest, the strongest, and the strongest electrical and thermal conductivity found so far, graphene is called by scientists as the "king of new materials" that may "completely change the 21st century". It has great potential for development in the fields of mobile devices, aerospace, and new energy batteries. Graphene is a quasi-two-dimensional material structure with a thickness of only one atomic layer, with a thickness of about 0.335 nm. It will affect its excellent electrical and mechanical properties. Therefore, it is necessary to eliminate the wrinkling phenomenon in the graphene structure in many applications.

薄膜褶皱是一个高度几何非线性的后屈曲现象,在建筑领域通常可以通过改变外荷载或约束边界条件来引入双轴应力加以消除。然而,在航天结构或者纳米材料结构中,由于空间展开、轻质量、纳米制造技术等方面的限制,上述的方法在航天薄膜结构和石墨烯结构中几乎不可能实现。因此,如何在不改变外荷载和约束边界条件的情况下,只通过改变结构自身的拓扑和形状达到消除褶皱的目的,无疑是一个非常重要但具有挑战性的问题。在一些已有研究和工程应用中,多采取经验或者试验方法对部分简单薄膜结构进行外形设计,无法进行推广和应用。对于具有复杂荷载或约束边界的薄膜结构,亟需开展一种以完全消除褶皱为目的的普适性优化设计方法,自动、准确、高效地寻找创新性拓扑形式,实现薄膜类结构的无皱化。Membrane wrinkling is a highly geometrically nonlinear post-buckling phenomenon that can usually be eliminated in the architectural field by introducing biaxial stresses by changing external loads or restraining boundary conditions. However, in aerospace structures or nanomaterial structures, due to the limitations of space expansion, light weight, and nanofabrication technology, the above methods are almost impossible to achieve in aerospace thin film structures and graphene structures. Therefore, how to eliminate wrinkles only by changing the topology and shape of the structure itself without changing the external loads and constraint boundary conditions is undoubtedly a very important but challenging problem. In some existing research and engineering applications, experience or experimental methods are mostly used to design the shape of some simple thin-film structures, which cannot be popularized and applied. For membrane structures with complex loads or constrained boundaries, it is urgent to develop a universal optimization design method for the purpose of completely eliminating wrinkles, to automatically, accurately and efficiently find innovative topological forms, and to achieve wrinkle-free membrane structures. .

发明内容SUMMARY OF THE INVENTION

本发明主要解决薄膜结构在拉伸荷载作用下容易产生褶皱现象的问题,提出一种以完全消除褶皱为目的的薄膜结构优化设计方法。在平面有限元分析基础上,通过控制每个单元的最小主应力为正值,调节薄膜的应力状态,采用拓扑优化技术对薄膜材料分布进行设计,进而获得带有曲线边界或孔洞的结构形式,以达到完全消除褶皱的目的。这种方式保证了对褶皱的完全消除,可实现外形和孔洞的精确定位,优化设计的自动化程度高,将保证薄膜结构无皱设计的研发效率。The invention mainly solves the problem that the film structure is prone to wrinkles under the action of tensile load, and proposes an optimal design method for the film structure for the purpose of completely eliminating wrinkles. On the basis of plane finite element analysis, by controlling the minimum principal stress of each element to be a positive value, the stress state of the film is adjusted, and topology optimization technology is used to design the film material distribution, and then the structural form with curved boundaries or holes is obtained, In order to achieve the purpose of completely eliminating wrinkles. This method ensures the complete elimination of wrinkles, can achieve precise positioning of shapes and holes, and has a high degree of automation in the optimized design, which will ensure the research and development efficiency of the wrinkle-free design of the film structure.

为了达到上述目的,本发明的技术方案为:In order to achieve the above object, the technical scheme of the present invention is:

一种用于消除薄膜结构拉伸褶皱的优化设计方法,具体包括以下步骤:An optimized design method for eliminating tensile folds of a film structure, which specifically includes the following steps:

第一步,对薄膜结构进行无皱化拓扑优化First step, wrinkle-free topology optimization of the thin film structure

(1)根据结构的尺寸要求和实际加载情况确定设计域,建立薄膜结构拓扑优化初始设计;在设计域中施加荷载和约束边界,划分有限元单元网格;(1) Determine the design domain according to the size requirements of the structure and the actual loading situation, and establish the initial design of the topology optimization of the membrane structure; apply loads and constraint boundaries in the design domain, and divide the finite element mesh;

(2)建立薄膜结构无皱化拓扑优化模型,使薄膜结构的整体刚度最大化或者整体柔顺性最小化(2) Establish a wrinkle-free topology optimization model of the film structure to maximize the overall stiffness or minimize the overall flexibility of the film structure

a)约束一:每个有限元单元最小主应力大于零,即

Figure BDA0001204869340000031
其中,e为有限元单元编号,σ1为最大主应力,σ2为最小主应力;a) Constraint 1: The minimum principal stress of each finite element element is greater than zero, namely
Figure BDA0001204869340000031
Among them, e is the finite element element number, σ 1 is the maximum principal stress, σ 2 is the minimum principal stress;

b)约束二:确定薄膜面积用量,作为面积约束上限;所述的薄膜面积用量为设计域面积的60%-90%。b) Constraint 2: Determine the area dosage of the film as the upper limit of the area constraint; the film area dosage is 60%-90% of the design domain area.

c)设计变量:设计域内有限元单元的相对密度ρe,ρe取值为α和1之间,代表单元处薄膜材料的分布;其中,α是一个远小于1的正数;c) Design variables: the relative density ρ e of the finite element element in the design domain, ρ e is between α and 1, representing the distribution of the film material at the element; where α is a positive number far less than 1;

(3)根据步骤(2)建立的拓扑优化模型,对约束一等效变换为

Figure BDA0001204869340000032
其中I1和J2分别为有限元单元的应力第一和第二不变量;(3) According to the topology optimization model established in step (2), the constraint-equivalent transformation is
Figure BDA0001204869340000032
where I 1 and J 2 are the stress first and second invariants of the finite element element, respectively;

(4)对步骤(3)变换后的约束采用约束松弛处理,避免应力奇异解现象;采用SIMP惩罚策略和优化算法进行迭代求解,得到薄膜结构最优材料分布;(4) Constraint relaxation processing is adopted for the transformed constraints in step (3) to avoid the phenomenon of stress singular solution; SIMP penalty strategy and optimization algorithm are used for iterative solution to obtain the optimal material distribution of the thin film structure;

第二步,对薄膜结构进行形状优化设计The second step is to optimize the shape of the thin film structure

在第一步(4)得到的薄膜结构拓扑形式的基础上,对薄膜结构边界及孔洞的具体几何参数进行优化,考虑最小主应力约束,获得更加详细和准确的结构形状参数。得到的结构形式满足最小主应力为正的要求和最大化整体刚度,构型比较简单,易于加工制造,经有非线性后屈曲限元分析和试验验证,实现真正的无皱化。On the basis of the topological form of the thin film structure obtained in the first step (4), the specific geometric parameters of the thin film structure boundary and holes are optimized, and the minimum principal stress constraints are considered to obtain more detailed and accurate structural shape parameters. The obtained structural form satisfies the requirement that the minimum principal stress is positive and maximizes the overall stiffness, the configuration is relatively simple, and it is easy to manufacture.

本发明的有益效果为:The beneficial effects of the present invention are:

在薄膜结构优化前,该类结构在拉伸载荷作用下存在显著的大面积褶皱现象。采用本发明的“无皱化”结构形式的薄膜后,通过数值模拟和试验考核,整个薄膜的最小主应力均保证为大于零的正值,完全观测不到褶皱产生。结构容易加工制造,只需要进行简单的裁剪和开孔。同时,本发明所建立的无皱化优化方法避免了在优化过程中进行复杂的后屈曲计算,所耗费的工作量极小,将显著提高设计效率,有望成为航空航天领域、微纳米领域中薄膜结构创新性设计的有效方法。Before the optimization of the film structure, this kind of structure has a significant large-area wrinkling phenomenon under the action of tensile load. After adopting the film of the present invention in the form of "wrinkle-free" structure, through numerical simulation and test evaluation, the minimum principal stress of the entire film is guaranteed to be a positive value greater than zero, and no wrinkles are observed at all. The structure is easy to process and manufacture, only simple cutting and opening are required. At the same time, the wrinkle-free optimization method established by the present invention avoids the complicated post-buckling calculation in the optimization process, and consumes a very small amount of work, which will significantly improve the design efficiency, and is expected to become a thin film in the aerospace field and micro-nano field. Effective methods for the innovative design of structures.

附图说明Description of drawings

图1为本发明实施例提供的一种四角拉伸结构设计域。FIG. 1 is a design domain of a four-corner tensile structure provided by an embodiment of the present invention.

图2(a)为薄膜面积比为70%时四角拉伸结构最优设计图。Figure 2(a) shows the optimal design of the four-corner tensile structure when the film area ratio is 70%.

图2(b)为薄膜面积比为80%时四角拉伸结构最优设计图。Figure 2(b) shows the optimal design of the four-corner tensile structure when the film area ratio is 80%.

图3为本发明实施例提供的一种中间具有硬块的两端拉伸结构设计域。FIG. 3 is a design domain of a tensile structure at both ends with a hard block in the middle according to an embodiment of the present invention.

图4为薄膜面积比为80%时中间具有硬块的两端拉伸结构最优设计图。Figure 4 is the optimal design diagram of the tensile structure at both ends with a hard block in the middle when the area ratio of the film is 80%.

具体实施方式Detailed ways

以下结合技术方案和附图详细叙述本发明的具体实施例。The specific embodiments of the present invention will be described in detail below with reference to the technical solutions and the accompanying drawings.

第一步,对薄膜结构进行无皱化拓扑优化First step, wrinkle-free topology optimization of the thin film structure

(1)根据结构的尺寸要求和实际加载情况确定设计域,建立薄膜结构拓扑优化初始设计;在设计域中施加荷载和约束边界并划分有限元网格。图1为四角拉伸结构的设计域,划分有限元网格数N=6400个,图3为中间具有刚性硬块的两端拉伸结构设计域,划分有限元网格数N=5000个。两种初始结构在拉伸荷载作用下均存在明显的褶皱行为。(1) Determine the design domain according to the size requirements of the structure and the actual loading situation, and establish the initial design of the topology optimization of the membrane structure; apply loads and constraint boundaries in the design domain and divide the finite element mesh. Figure 1 shows the design domain of the four-corner tensile structure, with N=6400 finite element meshes, and Figure 3 shows the design domain of the tensile structure at both ends with rigid blocks in the middle, with N=5000 finite element meshes. Both initial structures have obvious folding behavior under tensile load.

(2)建立薄膜结构无皱化拓扑优化模型:(2) Establish a topological optimization model of thin film structure without wrinkling:

(a)目标:薄膜结构的整体刚度最大化或者整体柔顺性最小化。(a) Objective: Maximize the overall stiffness or minimize the overall compliance of the membrane structure.

(b)约束一:要求每个有限元单元的最小主应力大于零,即

Figure BDA0001204869340000041
其中e为有限元单元编号,σ1为最大主应力,σ2为最小主应力。(b) Constraint 1: The minimum principal stress of each finite element element is required to be greater than zero, namely
Figure BDA0001204869340000041
where e is the finite element element number, σ 1 is the maximum principal stress, and σ 2 is the minimum principal stress.

(c)约束二:给定薄膜面积用量,作为面积约束上限。所述的薄膜面积用量为整个设计域面积的70%和80%。(c) Constraint 2: The amount of film area is given as the upper limit of the area constraint. The area usage of the film is 70% and 80% of the entire design domain area.

(d)设计变量:设计域内单元的相对密度ρe,取值为α=0.001和1之间,代表单元处薄膜材料的分布。(d) Design variable: the relative density ρ e of the element in the design domain, which is between α=0.001 and 1, representing the distribution of the thin film material at the element.

(3)根据步骤(2)建立的拓扑优化模型,对约束一等效变换为

Figure BDA0001204869340000051
其中I1和J2分别为应力第一和第二不变量。(3) According to the topology optimization model established in step (2), the constraint-equivalent transformation is
Figure BDA0001204869340000051
where I 1 and J 2 are the stress first and second invariants, respectively.

(4)对步骤(3)变换后的约束采用约束松弛处理,避免应力奇异解现象。采用cosine-type松弛方法,其中cosine-type松弛函数的表达式为θ=(1-cos(ρe·π))/2。(4) Constraint relaxation is used for the transformed constraints in step (3) to avoid the phenomenon of stress singular solution. The cosine-type relaxation method is adopted, wherein the expression of the cosine-type relaxation function is θ=(1-cos(ρ e ·π))/2.

(5)优化模型采用SIMP惩罚策略和优化算法进行迭代求解,得到薄膜结构最优材料分布,分别见图2和图4。(5) The optimization model adopts the SIMP penalty strategy and the optimization algorithm to iteratively solve, and obtains the optimal material distribution of the thin film structure, as shown in Figure 2 and Figure 4, respectively.

第二步,对薄膜结构进行细节形状优化设计The second step is to carry out the detailed shape optimization design of the thin film structure

在第一步(5)得到的薄膜结构拓扑的基础上,对薄膜结构边界和孔洞的具体几何参数进行优化,考虑最小主应力约束,获得更加详细和准确的结构形状参数。得到的构型比较简单,易于加工制造,经有非线性后屈曲限元分析和全尺寸试验验证,整个薄膜结构没有褶皱发生,验证了本发明所提方法的有效性。On the basis of the thin film structure topology obtained in the first step (5), the specific geometric parameters of the thin film structure boundaries and holes are optimized, and the minimum principal stress constraints are considered to obtain more detailed and accurate structural shape parameters. The obtained configuration is relatively simple, easy to manufacture, and has been verified by nonlinear post-buckling finite element analysis and full-scale test that no wrinkles occur in the entire film structure, which verifies the effectiveness of the method proposed in the present invention.

本发明的实质是采用拓扑优化方法得到具有曲边边界或孔洞的构型,以对整个薄膜最小主应力进行控制和优化,达到无皱化目的。其对前述各实施例所记载的优化模型、方法、方案进行修改,或者对其中部分或者全部方法特征进行等同替换(例如采用水平集或者显式曲线对结构边界与孔洞进行描述、采用其他拓扑优化方法、改变目标函数或约束具体形式等),并不使相应方法与方案的本质脱离本发明各实施例方法与方案的范围。The essence of the present invention is to use the topology optimization method to obtain a configuration with curved edge boundaries or holes, so as to control and optimize the minimum principal stress of the entire film to achieve the purpose of no wrinkling. It modifies the optimization models, methods, and schemes described in the foregoing embodiments, or performs equivalent replacements on some or all of the method features (for example, using level sets or explicit curves to describe structural boundaries and holes, using other topology optimization methods) method, changing the objective function or constraining the specific form, etc.), does not make the essence of the corresponding method and scheme deviate from the scope of the method and scheme of each embodiment of the present invention.

Claims (4)

1.一种用于消除薄膜结构拉伸褶皱的优化设计方法,其特征在于以下步骤:1. an optimized design method for eliminating stretch folds of film structure is characterized in that the following steps: 第一步,对薄膜结构进行无皱化拓扑优化First step, wrinkle-free topology optimization of the thin film structure (1)根据结构的尺寸要求和实际加载情况确定设计域,建立薄膜结构拓扑优化初始设计;在设计域中施加荷载和约束边界,划分有限元单元网格;(1) Determine the design domain according to the size requirements of the structure and the actual loading situation, and establish the initial design of the topology optimization of the membrane structure; apply loads and constraint boundaries in the design domain, and divide the finite element mesh; (2)建立薄膜结构无皱化拓扑优化模型,使薄膜结构的整体刚度最大化或者整体柔顺性最小化(2) Establish a wrinkle-free topology optimization model of the film structure to maximize the overall stiffness or minimize the overall flexibility of the film structure a)约束一:每个有限元单元最小主应力大于零,即
Figure FDA0003815383010000011
其中,e为有限元单元编号,σ2为最小主应力;
a) Constraint 1: The minimum principal stress of each finite element element is greater than zero, namely
Figure FDA0003815383010000011
Among them, e is the number of the finite element element, and σ 2 is the minimum principal stress;
b)约束二:确定薄膜面积用量,作为面积约束上限;所述的薄膜面积用量为设计域面积的60%-90%;b) Constraint 2: Determine the amount of film area as the upper limit of the area constraint; the amount of film area is 60%-90% of the design domain area; c)设计变量:设计域内有限元单元的相对密度ρe,ρe取值为α和1之间,代表单元处薄膜材料的分布;其中,α是一个远小于1的正数;c) Design variables: the relative density ρ e of the finite element element in the design domain, ρ e is between α and 1, representing the distribution of the film material at the element; where α is a positive number far less than 1; (3)根据步骤(2)建立的拓扑优化模型,对约束一等效变换为
Figure FDA0003815383010000012
其中I1和J2分别为有限元单元的应力第一和第二不变量;
(3) According to the topology optimization model established in step (2), the constraint-equivalent transformation is
Figure FDA0003815383010000012
where I 1 and J 2 are the stress first and second invariants of the finite element element, respectively;
(4)对步骤(3)变换后的约束采用约束松弛处理,避免应力奇异解现象;采用SIMP惩罚策略和优化算法进行迭代求解,得到薄膜结构最优材料分布;(4) Constraint relaxation processing is adopted for the transformed constraints in step (3) to avoid the phenomenon of stress singular solution; SIMP penalty strategy and optimization algorithm are used for iterative solution to obtain the optimal material distribution of the thin film structure; 第二步,对薄膜结构进行形状优化设计The second step is to optimize the shape of the thin film structure 在得到的薄膜结构拓扑的基础上,对薄膜结构边界及孔洞的具体几何参数进行优化,考虑最小主应力约束,获得更加详细和准确的结构形状参数。On the basis of the obtained thin film structure topology, the specific geometric parameters of the thin film structure boundary and holes are optimized, and the minimum principal stress constraints are considered to obtain more detailed and accurate structural shape parameters.
2.根据权利要求1所述的优化设计方法,其特征在于,第一步中所述的约束松弛处理包括ε松弛方法或qp松弛方法。2 . The optimization design method according to claim 1 , wherein the constraint relaxation process in the first step includes an ε relaxation method or a qp relaxation method. 3 . 3.根据权利要求1所述的优化设计方法,其特征在于,所述的约束松弛处理包括cosine-type松弛方法,其中cosine-type松弛函数的表达式为θ=(1-cos(ρe·π))/2。3. The optimization design method according to claim 1, wherein the constraint relaxation processing comprises a cosine-type relaxation method, wherein the expression of the cosine-type relaxation function is θ=(1-cos(ρ e · π))/2. 4.根据权利要求1或2或3所述的优化设计方法,其特征在于,所述的优化算法为准则法、MMA算法、ESO方法或Level set方法。4. The optimization design method according to claim 1, 2 or 3, wherein the optimization algorithm is a criterion method, an MMA algorithm, an ESO method or a Level set method.
CN201710011463.1A 2017-01-09 2017-01-09 Optimal design method for eliminating tensile wrinkles of film structure Expired - Fee Related CN106777768B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710011463.1A CN106777768B (en) 2017-01-09 2017-01-09 Optimal design method for eliminating tensile wrinkles of film structure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710011463.1A CN106777768B (en) 2017-01-09 2017-01-09 Optimal design method for eliminating tensile wrinkles of film structure

Publications (2)

Publication Number Publication Date
CN106777768A CN106777768A (en) 2017-05-31
CN106777768B true CN106777768B (en) 2022-09-30

Family

ID=58951193

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710011463.1A Expired - Fee Related CN106777768B (en) 2017-01-09 2017-01-09 Optimal design method for eliminating tensile wrinkles of film structure

Country Status (1)

Country Link
CN (1) CN106777768B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108133097B (en) * 2017-12-21 2021-01-19 大连理工大学 Space film structure clamp shape optimization design method for inhibiting wrinkles
CN111285336B (en) * 2018-12-10 2023-02-03 北京清正泰科技术有限公司 Method for eliminating two-dimensional material wrinkles
CN112800642B (en) * 2020-12-29 2022-07-19 中国人民解放军国防科技大学 Film wrinkle strength evaluation method and system based on nonlinear finite element
CN112784457B (en) * 2021-01-19 2022-05-20 中国人民解放军国防科技大学 Layout optimization method and system of thin-film structure reinforcement strips based on genetic algorithm
CN113221201B (en) * 2021-05-12 2023-08-04 西安电子科技大学 Folding design method of small-curvature paraboloid-of-revolution film

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5746446A (en) * 1994-05-23 1998-05-05 Automotive Technologies International, Inc. Plastic film airbag
WO2010001401A1 (en) * 2008-07-03 2010-01-07 Cadent Ltd. Method, apparatus and system for use in dental procedures
WO2011050161A1 (en) * 2009-10-21 2011-04-28 The Trustees Of Columbia University In The City Of New York Systems and methods for self-assembling ordered three-dimensional patterns by buckling of thin films bonded to curved compliant substrates
CN102043883A (en) * 2010-12-29 2011-05-04 长沙理工大学 Material breakage constraint-based continuum structure topology design modeling and optimization design method
CN103020377A (en) * 2012-12-26 2013-04-03 国电联合动力技术有限公司 Bolt strength analysis method for directly-driven wind-driven generator
CN103020361A (en) * 2012-12-17 2013-04-03 华南理工大学 Method for extracting no-checkerboard topological diagram from compliant mechanism
US8655481B2 (en) * 2010-04-09 2014-02-18 Victor Shi-Yueh Sheu IMR (in-mold roller or in-mold release)/IMF (in-mold forming) making method using a digital printer printing and pre-forming technique
WO2014124044A1 (en) * 2013-02-06 2014-08-14 The Board Of Trustees Of The University Of Illinois Self-similar and fractal design for stretchable electronics
EP3096195A1 (en) * 2015-05-19 2016-11-23 Benecke-Kaliko AG Method for controlling a machine tool
CN106299583A (en) * 2016-08-15 2017-01-04 西安电子科技大学 The spaceborne Electrostatic deformation film reflector face anti-discharge method of attaching of deployable antenna electrode

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014011222A1 (en) * 2012-07-13 2014-01-16 Massachusetts Institute Of Technology Thin films with micro-topologies prepared by sequential wrinkling

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5746446A (en) * 1994-05-23 1998-05-05 Automotive Technologies International, Inc. Plastic film airbag
WO2010001401A1 (en) * 2008-07-03 2010-01-07 Cadent Ltd. Method, apparatus and system for use in dental procedures
WO2011050161A1 (en) * 2009-10-21 2011-04-28 The Trustees Of Columbia University In The City Of New York Systems and methods for self-assembling ordered three-dimensional patterns by buckling of thin films bonded to curved compliant substrates
US8655481B2 (en) * 2010-04-09 2014-02-18 Victor Shi-Yueh Sheu IMR (in-mold roller or in-mold release)/IMF (in-mold forming) making method using a digital printer printing and pre-forming technique
CN102043883A (en) * 2010-12-29 2011-05-04 长沙理工大学 Material breakage constraint-based continuum structure topology design modeling and optimization design method
CN103020361A (en) * 2012-12-17 2013-04-03 华南理工大学 Method for extracting no-checkerboard topological diagram from compliant mechanism
CN103020377A (en) * 2012-12-26 2013-04-03 国电联合动力技术有限公司 Bolt strength analysis method for directly-driven wind-driven generator
WO2014124044A1 (en) * 2013-02-06 2014-08-14 The Board Of Trustees Of The University Of Illinois Self-similar and fractal design for stretchable electronics
EP3096195A1 (en) * 2015-05-19 2016-11-23 Benecke-Kaliko AG Method for controlling a machine tool
CN106299583A (en) * 2016-08-15 2017-01-04 西安电子科技大学 The spaceborne Electrostatic deformation film reflector face anti-discharge method of attaching of deployable antenna electrode

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
An enhanced aggregation method for topology optimization with local stress constraints;Yangjun luo;《Computer Methods in Applied Mechanics and Engineering》;20130228;31-41 *
向量式有限元膜单元及其在膜结构褶皱分析中的应用;赵阳;《建筑结构学报》;20150120;127-135 *
薄膜结构褶皱研究述评;李作为等;《中国安全科学学报》;20040730(第07期);全文 *
铝合金温室拱形ETFE膜屋面抗风性能分析;乔克等;《农业工程学报》;20120715(第14期);全文 *

Also Published As

Publication number Publication date
CN106777768A (en) 2017-05-31

Similar Documents

Publication Publication Date Title
WO2018126465A1 (en) Optimization design method for removing tensile wrinkles from thin-film structure
CN106777768B (en) Optimal design method for eliminating tensile wrinkles of film structure
Du et al. Connecting microstructures for multiscale topology optimization with connectivity index constraints
Ye et al. Topology optimization of conformal structures on manifolds using extended level set methods (X-LSM) and conformal geometry theory
Luo et al. Wrinkle-free design of thin membrane structures using stress-based topology optimization
Hao et al. Design of manufacturable fiber path for variable-stiffness panels based on lamination parameters
CN109657284A (en) A kind of equal geometry Topology Optimization Method towards Meta Materials
Zhang et al. Geometrically nonlinear large deformation analysis of functionally graded carbon nanotube reinforced composite straight-sided quadrilateral plates
Kwok et al. Four-dimensional printing for freeform surfaces: design optimization of origami and kirigami structures
CN107220461B (en) A high-efficiency optimization method for composite shell structures with variable stiffness
CN110222432B (en) A Genetic Algorithm-Based Parameter Optimal Design Method for Locally Constrained Damping Plates
CN107145657B (en) Non-uniform geometric variable block method for reactor neutron diffusion equation
Albanesi et al. An optimization method based on the evolutionary and topology approaches to reduce the mass of composite wind turbine blades
Albanesi et al. Application of the inverse finite element method to design wind turbine blades
CN111523264B (en) Multiphase material microstructure topology optimization method with ultimate elastic performance
CN105426592B (en) A kind of Electrostatic deformation film reflector surface antenna analysis method
Timme et al. Towards three-dimensional global stability analysis of transonic shock buffet
Lee et al. Finite element analysis of wrinkling membrane structures with large deformations
Naumenko et al. A layer-wise theory of shallow shells with thin soft core for laminated glass and photovoltaic applications
CN107590325A (en) A kind of fiber-reinforced composite materials structures optimization method based on Shepard interpolation
Vigliotti et al. Analysis and design of lattice materials for large cord and curvature variations in skin panels of morphing wings
Nguyen-Van et al. Geometrically nonlinear analysis of composite plates and shells via a quadrilateral element with good coarse-mesh accuracy
Gu et al. Form-finding design of electrostatically controlled deployable membrane antenna based on an extended force density method
Iwasa Approximate estimation of wrinkle wavelength and maximum amplitude using a tension-field solution
Nie et al. Efficient prediction of the effective nonlinear properties of porous material by FEM-Cluster based Analysis (FCA)

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220930

CF01 Termination of patent right due to non-payment of annual fee