The present invention relates to a method of making one or more fibrils, the method comprising the steps of providing a material of manufacture of the fibrils; providing a random initial shape of the fibrils; calculating an adhesive force of the fibrils based on the material of manufacture of the fibrils and on the initial shape of the fibrils; adapting the random initial shape of the fibrils to maximize the adhesive force of the fibrils to form resultant shapes of the fibrils; selecting the resultant shape of the fibrils having the highest adhesive force of the fibrils; and producing having the selected resultant shape having the highest adhesive force of the fibrils. The invention further relates to a computer implemented method of simulating an adhesive force of one or more fibrils and to a fibril.
Geckos can climb almost any surface in nature swiftly and effortlessly. Such incredible climbing ability comes from the hair-like microstructures, setae, on their feet. The setae branch into many smaller hairs of seta with a spatulated tip ending, where these tips directly contact and adhere to surfaces leveraging mainly the most universal intermolecular force, the van der Waals force. The hierarchical and branching fibrillar structure of setae enable exploiting contact splitting and crack trapping and adhesion to even rough surfaces. Such fascinating adhesion evoked an inspiration for a synthetic elastomeric gecko pad. The research on synthetic gecko pad resulted in many follow-up studies on dry fibril adhesives, such as contact splitting and crack trapping, adhesion controllability, self-cleaning property, hierarchical structuring, multi-material composite fibrils, adhesion with any liquid super-repellency, adhesion to skin, soft robotic grippers, and robotic wall climbers.
The adhesion strength of synthetic gecko fibril arrays is studied typically on smooth surfaces, and it is a function of the size and shape of the fibril tip and stem in addition to the fibril density, fibril elastomeric material properties, and the fibril array backing layer. For a single fibril, the adhesion strength is especially sensitive to the geometry of the fibril structure. It is found that the tip geometry is crucial to the adhesion improvement, where bioinspired mushroom-shaped and T-shaped fibril tips improved the adhesion strength notably greater than the one from simple flat-punch fibrils. Also, multiple optimal mushroom- and T-shaped tip geometries have been proposed.
Limitations arise when optimizing such previously proposed tip shapes due to the nature of the current template-based design approach. Such approach is based on a pre-determined T-shaped or mushroom-shaped tip with a cylindrical stem, and suppresses the chances to explore and find out-of-the-box optimal tip and stem shapes.
For this reason, it is an object of the invention to provide a method of making one or more fibrils that results in improved fibril designs having a higher adhesion force than prior art fibrils of the same material. It is yet a further object of the invention to reduce the time period of manufacture and designing such fibrils.
This object is satisfied by the subject matter of the independent claims.
A method of making one or more fibrils, the method comprising the steps of:
In this way synthetic elastomeric dry microfibrillar adhesives can be created for a wide range of applications the fibers or arrays of fibers having the following unique properties, they are: residue-free, repeatable, tunable, controllable and silent adhesion; self-cleaning property; and breathability.
In contrast to the prior art where the design of such dry fibrillar adhesives is limited by a template-based design approach using a pre-determined bioinspired T-shaped or wedge-shaped mushroom tips, the present disclosure not only allow to vary the material of the fibrils, but also varies the shape and size of the fibrils in order to obtain fibrils of a given material having improved adhesion characteristics.
In this connection it should be noted that the step of adapting the initial random shape of the fibrils may take place a plurality of times, where each consecutive step may be based on the previous step. That is, the various adapting steps can be based on one another. Hence, the adapting of the initial—and therefore also the resultant shape—can also take place a plurality of times until the determined corresponding adhesive force of the resultant shapes converges.
The initial shape of the one or more fibrils may comprise an axisymmetric shape having an at least locally flat tip surface at an end thereof and the profile of the one or more fibrils are defined by a set of parameterizable curves, such as Bezier curves, spline curves, or polynomial curves.
This means that the initial shape can be an arbitrarily shaped fibril having no pre-defined specific shape. It may have a predefined shape such as one of a cylindrical shaped fibril, a T-shaped fibril, a mushroom shaped fibril, a wedge shaped fibril or the like. Preferably a shape is initially selected that has previously exhibited acceptable adhesion force characteristics. The profile of this shape is then varied in order to obtain the ideal shape for a given material and size of fibril.
The step of calculating the adhesive force of the one or more fibrils may take place by means of a computer implemented simulation. In this way the designing of the fibrils can be sped up significantly.
The normal stress occur at the interface determines the critical point of the separation between the fibril(s) and the object to be attached and the surface integral of the normal stress at the bottom of the stem of the fibril becomes the simulated adhesion. The calculation of the adhesion follows the standard finite element method and may be based on Cauchy's equations.
Thus, the method proposed herein teaches a machine learning (ML)-based computational design approach, where the 3D tip and stem shape of an adhesive fibril is searched in a broad design space and adapted automatically at a given fibril diameter scale and for a given material.
ML-based design optimizations have been widely used in various areas, such as optimal antenna design, composite design, mechanical property prediction, polymer dielectrics, aircraft airfoil design, mechanical structure design, and energy materials. Most of these studies rely on approximate models of finite-element-method (FEM)-based simulations by training artificial-neural-networks (ANNs) in order to accelerate the simulation evaluations. Improved evaluation speed of the approximated models makes the global optimization methods, such as genetic algorithms, feasible. However, the main drawbacks of such approach are high computational time and cost to create the training dataset for the ANNs and lower accuracy due to the approximation process.
In the present teaching a ML-based adhesive fibril design method with a data-efficient optimization framework and accurate FEM-based adhesion mechanics simulations is proposed. By employing a state-of-the-art global optimization method, Bayesian optimization, with comprehensive finite-element-method (FEM)-based adhesion mechanics simulations, the method produces optimal shapes of fibrils, which outperform previous adhesive fibril designs in both simulations and experiments.
Compared to a recently proposed machine learning approach, this disclosure provides the first experimental validation of fibrils designed by machine learning methods together with in-depth stress analyses. In addition, the proposed method is much more data efficient to generate a sub-optimal design in a reasonable duration (˜3 hours for each design in this study) and flexible for the changes in the simulation, where the governing physics of the simulation could be readily changed without any prior training of a model, allowing the creation of different adhesive fibrils with other adhesion mechanisms, such as suction. Finally, the method also optimizes the fibril 3D stem shape and deformation in addition to the tip shape, which is ignored in almost all previous template-based and in the previous machine learning-based fibril optimization method.
The step of adapting the initial shape of the one or more fibrils may take place a plurality of times, i.e. iteratively to obtain a wide parameter range from which to select the desired fibrils. In this connection it should be noted that the plurality of times may be selected between 2 and 1000000 iterations, in particular between 5 and 100000 iterations.
In some embodiments the step of adapting the initial shape of the one or more fibrils may take place until the determined adhesive forces of the resultant shapes converge, i.e. until the determined adhesive forces do not vary significantly anymore.
The step of adapting the initial shape of the one or more fibrils takes place by varying the smallest edge radius of the fibril tip with a size selected in the range of 1 nm to 10 μm. In this way the proposed method also takes into account real fabrication methods of the fibrils, i.e. the, preferably outer, shape and size may be determined by the fabrication process of the one or more fibrils. This is in stark contrast to the prior art where the edges are assumed to be of perfect 90° angles and hence leads to a distortion of the internal stresses of the fibrils.
The step of adapting the initial shape of the one or more fibrils may also take place by means of a computer implemented simulation. In this way the designing of the fibrils can be sped up significantly.
The shape adaptation may be based on probabilistic optimization methods, such as Bayesian optimization, which maximizes the objective function value, the adhesion of the one or more fibrils, while keeping the number of required finite element method (FEM) simulation evaluation low.
The step of producing the one or more fibrils may take place by means of at least one of additive manufacturing, two-photon polymerization, 3D printing, optical lithography, electron-beam lithography, focused ion beam machining, laser micro/nanomachining, mechanical or ultrasound micromachining, micro/nanoprinting, roll-to-roll replication, injection molding, compression molding, and polymer casting.
In this connection it should be noted that the fibrils can be formed directly at a backing substrate made of the same or a different material from the material of the one or more fibrils.
The material of manufacture of the one or more fibrils may be selected from the group of members consisting of organic, inorganic, polymers, rubbers, silicones, polyurethanes, biomaterials, biopolymers, composites, elastomers, liquid crystalline elastomers, thermoplastic elastomers, foams, fabric materials, particle materials, fibrous materials, or combinations of the foregoing. In this connection it is noted that the, preferably soft, material of the fibrils is preferably an elastomer, such as PDMS or a polyurethane elastomer.
It is preferred that the fibrils and/or the backing substrate further comprise a filler material selected from the group of members consisting of organics, inorganics, metals, alloys, ceramics, glass, polymers, rubbers, biomaterials, hydrogels, liquid materials, phase changing materials, composites, foams, fabric materials, particle materials, fibrous materials, micro- and nanoparticles and combinations of the foregoing. An example for the use of filler materials could be magnetic micro-particles, which enable the material to be movable, or metal particles, which can make the material switchable or electrically conductive so that it could e.g. be used as part of a sensor. It could also be possible to introduce a type of functionalized material in order to give the backing substrate of the one or more fibrils, for example, a further function apart from holding the fibrils together.
Thus, for example electrodes could form the backing substrate that then adheres to a patient's skin via the fibrils to make wearable medical devices.
The material of manufacture of the one or more fibrils may have a Young's modulus selected in the range of 0.01 to 10000 MPa, where the Young's modulus could be measured via tensile stress tests. Such comparably soft materials are well suited to be used for a wide range of adhesive applications in robotics, transfer printing, closures, consumer products, and medical devices etc.
In order to measure the Young's modulus of a material, tensile stress testing can be done with a commercial universal tensile testing machine such as an Instron tensile tester. Tensile testing, also known as tension testing, is a test, in which a material sample is subjected to a controlled tension until failure. The properties that are directly measured via a tensile test are ultimate tensile strength, breaking strength, maximum elongation and reduction in area. From these measurements—among other things—the Young's modulus can be determined.
The test process involves placing the test specimen in the tensile testing machine and slowly extending it until it fractures. During this process, the elongation of the gauge section is recorded against the applied force. The data can be manipulated such that it is no longer specific to the geometry of the test sample. The elongation measurement is used to calculate the engineering strain, ε, using the following equation:
ε=ΔL/L0=(L−L0)/L0
wherein ΔL is the change in gauge length, L0 is the initial gauge length, and L is the final length. The force measurement is used to calculate the engineering stress, σ, using the following equation:
σ=Fn/A
wherein F is the tensile force and A is the nominal cross-section of the specimen.
The Young's modulus can then be calculated via the following equation:
E≤σ/ε.
The method may further comprise the step of designing the fibril tip edge by defining a radius of curvature between one or more surfaces of the one or more fibrils.
The method may further comprise the steps of:
In this way the adhesion force at the boundary interface can be calculated in advance resulting in even more improved adhesion characteristics of the one or more fibrils to predefined surfaces.
The step of calculating an interfacial stress between the one or more fibrils and the boundary interface may take place by means of a computer implemented simulation.
The interfacial stress may be calculated directly by the Cauchy's equilibrium equations when the fibril in the simulation is gradually stretched. When the maximum stress at the interface exceeds the interfacial strength between the fibril material and the object material, which is measured experimentally.
According to a further aspect the present invention also relates to a computer implemented method of simulating an adhesive force of one or more fibrils to produce fibrils with an increased adhesive force, the method comprising the steps of:
In this connection the initial shape of the one or more fibrils can be a random initial shape that is arbitrarily selected. It may also be an initial shape that resembles the shape of a previously produced fibril, such as a T-shaped fibril, a wedge shaped fibril, a mushroom shaped fibril and the like.
Thus, a machine learning-based computational approach to optimize designs of adhesive fibrils, exploring a much broader design space is obtained. A combination of Bayesian optimization and finite element methods creates novel optimal designs of adhesive fibrils, which may be fabricated by two-photon-polymerization-based three-dimensional (3D) microprinting and double-molding-based replication out of polydimethylsiloxane (PDMS).
Such optimal elastomeric fibril designs outperform previously proposed designs by maximum 77% in the experiments in respect of dry adhesion performance on smooth surfaces. Furthermore, finite element analyses reveal that the adhesion of the fibrils is sensitive to the 3D fibril stem shape, tensile deformation, and fibril microfabrication limits, which contrast with the previous assumptions that mostly neglect the deformation of the fibril tip and stem, and focus only on the fibril tip geometry.
In this connection it should be noted that the step of adapting the initial random shape of the fibrils may take place a plurality of times, where each consecutive step may be based on the previous step. That is, the various adapting steps can be based on one another. Hence, the adapting of the initial—and therefore also the resultant shape—can also take place a plurality of times until the determined corresponding adhesive force of the resultant shapes converges.
In this connection it should further be noted that it is common knowledge that any information that is provided to and/or defined by and/or handled by a computer system is usually available as some sort of data set, such as, for example, a sequence of zeros and ones. Thus, it can also be possible that each input and/or output of a computer is given in the form of a data set.
The proposed computational fibril design method is data-efficient and flexible to be combined with other fibril adhesion mechanisms, such as suction and multimaterial-based composite dry fibril adhesion, which could help design future optimal fibrils for specific applications with less help from human intuition.
According to an embodiment of the invention the step of iteratively adapting the initial shape or the previously determined shape or the previously determined shape having the highest adhesive force so far of the one or more fibrils is carried out by varying a smallest edge radius of the surfaces of the one or more fibrils, in particular with said smallest edge radius comprising a size in the range of 1 nm to 10 μm.
According to another embodiment the step of iteratively adapting the initial shape of the one or more fibrils is based on probabilistic optimization methods, such as Bayesian optimization, with the random initial shape in particular being a random 3rd-order Bezier-curve-based side profile.
It may further be possible that the probabilistic optimization method comprises the step of maximizing an objective function value, i.e. the adhesion of the one or more fibrils, while keeping the number of a required finite element method, FEM, simulation evaluation low, i.e. below 1000 simulation evaluations.
In this connection it should be noted that the optimization method may be defined by
where f(·) is the adhesion, in particular calculated from FEM, x is a parameter which constitutes profile curves of a set of parameterizable curves, such as Bezier curves, spline curves, or polynomial curves, and A is a simple set which constitutes a hyper-rectangle of the parameter boundaries of x.
Additionally, it can be possible that f(x) is a Gaussian Processes, GP, which is a collection of random variables such that a joint distribution of every finite subset of random variables is a multivariate Gaussian defined by:
f˜GP(μ,k), in particular f=GP(μ,k)
wherein μ(x) and k(x) are a mean and a covariance function, and wherein the GP surrogates the relationship between the parameter x and f(x) using black-box representation of the adhesion simulated in the FEM.
According to another embodiment of the invention the initial shape of the one or more fibrils comprises an axisymmetric shape having an at least locally flat tip surface at an end thereof and the profile of the one or more fibrils are defined by a set of parameterizable curves, such as Bezier curves, spline curves, or polynomial curves.
According to yet another embodiment the initial shape of the one or more fibrils resembles the shape of a previously produced shape, such as a T-shaped fibril, a wedge-shaped fibril, a mushroom-shaped fibril and/or combinations of the foregoing.
It may further be possible that the surfaces of the one or more fibrils are arranged such that the initial shape of the one or more fibril comprises an axisymmetric shape having an at least locally flat tip surface at an end thereof and the profile of the one or more fibrils are defined by a set of parameterizable curves, such as Bezier curves, spline curves, or polynomial curves.
The surfaces of the one or more fibrils may also be arranged such that the initial shape of the one or more fibrils resembles the shape of a previously produced shape, such as a T-shaped fibril, a wedge-shaped fibril, a mushroom-shaped fibril and/or combinations of the foregoing.
According to another embodiment the step of calculating the adhesive force of the one or more fibrils is carried out by standard finite element method.
In this connection it should be noted that the standard finite element method may be based on Cauchy's equation.
It may further be possible that the step of adapting the initial shape of the one or more fibrils takes place a plurality of times with the plurality of times being selected between 2 and 1000000 iterations, i.e. iteratively.
According to an embodiment the adapting the initial shape of the one or more fibrils stops when the calculated adhesive forces of the resultant shapes of the one or more fibrils converge.
According to a further embodiment the material of manufacture of the one or more fibrils is selected from the group of members consisting of organic, inorganic, polymers, rubbers, silicones, polyurethanes, biomaterials, biopolymers, composites, elastomers, liquid crystalline elastomers, thermoplastic elastomers, foams, fabric materials, particle materials, fibrous materials, or combinations of the foregoing. In this connection it is noted that the soft material is preferably an elastomer, such as PDMS, a silicone elastomer, vinylsiloxane, or a polyurethane elastomer.
In this connection it should be noted that the material of manufacture of the one or more fibrils may have a Young's modulus selected in the range of 0.01 to 10000 MPa.
According to a further embodiment of the invention the computer implemented method further comprises the step of designing the fibril tip edge by defining a radius of curvature between the surfaces of the one or more fibrils.
According to another embodiment the computer implemented method further comprises the steps of: defining a boundary interface to which the one or more fibrils are attachable, defining a material of the boundary interface; and calculating an interfacial stress between the one or more fibrils and the boundary interface.
In this connection it should be noted that the interfacial stress may be calculated by the Cauchy's equilibrium equations.
According to a further aspect the present invention also relates to a fibril obtainable by a method of making the one or more fibrils and/or simulated by the computer implemented method, the fibril optionally having at least one of an adhesive force selected in the range of 0.000001 mN to 1000 mN, a shape selected from the group of members comprising a profile composed of freeform curves; and interfacial stress to a boundary surface selected in the range of 0.001 kPa to 10000 MPa and a material selected from the group of members consisting of organic, inorganic, polymers, rubbers, silicones, polyurethanes, biomaterials, biopolymers, composites, elastomers, liquid crystalline elastomers, thermoplastic elastomers, foams, fabric materials, particle materials, fibrous materials, or combinations of the foregoing.
Such fibril designs can be used for a wide range of adhesive applications in robotics, transfer printing, closures, consumer products, and medical devices etc.
Further embodiments of the invention are described in the following description of the Figures and/or the dependent claims attached to this specification. The invention will be explained in the following in detail by means of embodiments and with reference to the drawing in which is shown:
In the following the same reference numerals will be used for parts having the same or equivalent function. Any statements made having regard to the direction of a component are made relative to the position shown in the drawing and can naturally vary in the actual position of application.
The present invention relates to a method of making one or more fibrils 10 (see
It is preferred that the fibrils 10 and/or a backing substrate 24 (see
Following the selection of the material, a random initial shape of the one or more fibrils 10 is selected. This initial shape of the one or more fibrils 10 comprises an axisymmetric shape having an at least locally flat tip surface at an end thereof and the profile of the one or more fibrils 10 are defined by a set of parameterizable curves, such as Bezier curves, spline curves, or polynomial curves. By way of example, the fibrils 10 shown in
Thereafter, an adhesive force of the one or more fibrils 10 is calculated. This adhesive force is based on the material of manufacture of the one or more fibrils 10 and on the initial shape of the one or more fibrils 10. The calculation of the adhesive force of the one or more fibrils 10 takes place by means of a computer implemented simulation.
This computer implemented simulation is based on the normal stress that occurs at the interface and that determines the critical point of the separation between the fibril(s) 10 and an object to be attached and the surface integral of the normal stress at the bottom of the stem of the fibril 10 becomes the simulated adhesion. The calculation of the adhesion is based on the standard finite element method which follows Cauchy's equations of equilibrium as
where σxx, σyy, σzz are normal stresses and τxy, τyz, τxz are shear stresses on an infinitesimal cube in x-, y-, and z-directional normal surfaces; x1, x2, x3 are coordinate system; and Fx, Fy, Fz are the body forces acting on the infinitesimal cube. The calculation of the equation could be performed using a weighted residual method in weak formulation.
The calculation of the adhesion is further processed from the stress calculations from above. The stresses calculated from the Cauchy's equations of equilibrium provide the interfacial stress between the fibril of interest and the adhesion interface. The integration of the normal stress along the interface or along the root (base and/or stem) of the fibril becomes the adhesion force of the fibril during the stretch of the fibril. When the normal stress at the interface exceeds the interfacial strength of the interface, the adhesion force at that moment becomes the adhesion of the fibril. The calculation of these could be formutated in commercial finite element analysis tools.
In order to carry out this step of adaption reliably the initial shape of the one or more fibrils 10 is adapted between 5 and 1000000 iterations.
Specifically, the step of adapting the initial shape of the one or more fibrils 10 takes place by varying the smallest edge radius of the fibril tip with a size selected in the range of 1 nm to 10 μm. This shape and size may be determined by the fabrication process of the one or more fibrils 10. For example, if the fibrils 10 are cast in a mold, then the radius may be of the order of 500 to 10000 nm, whereas fibrils 10 produced e.g. by means of polymer casting may have a radius of 500 nm to 10 μm and fibrils 10 produced e.g. by means of electron-beam lithography may produce fibrils with radii in the range of 1 nm to 5 μm.
The method may further comprise the step of designing the fibril tip edge by defining a radius of curvature between one or more surfaces of the one or more fibrils.
In this connection it should be noted that the representation of the profile may be performed in a parameterized vector representation. Because of vectorized representation of the profile, the designed profile is not subject to the resolution of the design system. By changing the parameters of the profile, the profile could be further adapted, where the optimal trials of the parameters are given by the Bayesian optimizer.
Thereafter, the shape and size of the random initial shape of the one or more fibrils 10 is adapted, i.e. varied in order to maximize the adhesive force of the one or more fibrils 10 to form resultant shapes of the one or more fibrils 10.
This step of adapting the initial shape of the one or more fibrils 10 also takes place by means of a computer implemented simulation. In order to carry out this step of adaption reliably the initial shape of the one or more fibrils 10 is adapted between 5 and 1000000 iterations in order to obtain different adhesives forces for the different fibril shapes.
The shape adaptation is based on probabilistic optimization methods, such as Bayesian optimization, which maximizes the objective function value, the adhesion of the one or more fibrils 10, while keeping the number of required finite element method (FEM) simulation evaluation low. The optimization process solves the general optimization represented as
where f(·) is the adhesion calculated from FEM, x is the parameter which constitutes the profile curves of the set of parameterizable curves, such as Bezier curves, spline curves, or polynomial curves, A is a simple set which constitutes a hyper-rectangle of the parameter boundaries of x.
The f(x) is represented by Gaussian Processes (GP) which is a collection of random variables such that the joint distribution of every finite subset of random variables is multivariate Gaussian:
f˜GP (μ,k), in particular f=GP(μ,k)
where μ(x) and k(x) are the mean and covariance function. The GP surrogates the relationship between the parameter x and f(x) using black-box representation of the adhesion simulated in FEM. By choosing different kinds of acquisition functions of the Bayesian optimizer, the shape of the fibrils is optimized based on the parameters generated by Bayesian optimizer.
Once the shape has been varied for the one or more fibrils 10 and their respective adhesion force has been simulated, one proceeds to select the resultant shape of the one or more fibrils 10 having the highest adhesive force of the one or more fibrils 10. Also, this step can be implemented in a computer in order to speed up the manufacturing process. Thereby, for a specific material of manufacture that is desired for a specific application of the fibrils 10 the best shape for the desired adhesion force of the fibrils 10 can be selected. For example, the highest adhesion force may be desired for the fibrils 10. Applications may also exist, where a pre-defined adhesion force may be desired at which the fibrils 10 automatically detach from a surface as a kind of pre-determined breaking point.
It should also be noted in this connection that for different sizes of fibrils 10 of the same material different shapes may yield the best adhesion force for the fibrils 10, i.e. smaller fibrils 10 of the same material may have a shape that differs from the shape of larger fibrils 10 of the same material in order to yield the optimum adhesion force for this size and material of fibril 10.
Following the selection of the resultant shape of the one or more fibrils 10, these are produced with the selected resultant shape having the highest adhesive force of the one or more fibrils 10 for the desired material of manufacture.
The step of producing the one or more fibrils takes place by means of at least one of additive manufacturing, two-photon polymerization, 3D printing, optical lithography, electron-beam lithography, focused ion beam machining, laser micro/nanomachining, mechanical or ultrasound micromachining, micro/nanoprinting, roll-to-roll replication, injection molding, compression molding, and polymer casting.
By way of example, the one or more fibrils can be produced using the following method: In a first step a rigid master 12 (see e.g.
In a second step the master 12 including its patterns 14 are first covered with a non-stick coating (not shown) and then with a first soft material 16. Thus, the patterns 14 are replicated with the first soft material 16 in order to produce a negative mold 18 comprising cavities 20 with an inner shape corresponding to the outer shape of the patterns 14, e.g. the double re-entrant shape, T-shaped, punch shaped, wedge shaped, mushroom shaped, arrow shaped, and mixed shaped fibrils of the fibrils to be produced. The first soft material 16 usually comprises a Young's modulus in the range of 10 kPa to 5 MPa. Also, higher values are possible as long as they still lie within a range where the material can be regarded as “soft”.
The “non-stick coating” is supposed to reduce the surface energy of the cavities 20 in order to simplify the further step of demolding the master 12 from the first soft material 16 (the mold 18). In commonly known methods for manufacturing double re-entrant shapes out of rigid materials a non-stick coating is not necessary since they can be produced with additive manufacturing processes. Here it is the case that a double re-entrant shape is supposed to be produced out of soft materials, which is why a negative mold 18 is needed.
Therefore, the surface of the master is first activated by providing oxygen molecules to the master. Said oxygen molecules are then illuminated with UV radiation which causes the oxygen to react with, for example, the carbon dioxide in the air to produce ozone. These ozone molecules are able to reach every edge and corner of the patterns to react and activate the surface. In a further step, passivation is initiated by providing a silane, preferably fluorosilane, or another chemical treatment to the activated surface, which can react with the activated surface in order to give it a sensation of a “non-sticky coating”, i. e. to help in a further step of demolding the master from the first material.
Another possibility would be to provide an oxygen plasma to the master surface to activate the surface before passivating it with fluorosilane or another chemical treatment.
In another embodiment of the invention, the surface and the patterns of the master are simply covered with a polymer coating. When using this method, a provision of UV radiation is not necessary.
In a further step, regardless of which “coating”-method has been chosen, before the master 12 is demolded, the first soft material 16 is cured since the material 16 is usually in a liquid phase when it is applied to the master 12. Depending on the exact material used for the first soft material 16 and depending on the ambient temperature said curing can be done within 10 minutes, when the temperature is raised to 150° C., within an hour, when the temperature is raised up to 90° C., or can even take up to 48 hours when the temperature is kept at room temperature (example of PDMS with a standard 10:1 monomer to crosslinker ratio). For softer formulations like, for example, a 20:1 ratio room temperature curing takes even longer.
In order to demold the master 12, without destroying it, from the first material 16 it can be mechanically peeled off by using conventional tweezers. Since the surface of the master 12 has been prepared with the activation and passivation steps or with the polymer coating, as described before, it can be peeled off rather easily from the first material 16 without destroying the delicate structures of the cavities 20 or the patterns 14. Other demolding techniques like chemical dissolving may generally be possible depending on the used materials but are nevertheless infeasible since they would destroy the master and make it not reusable.
After the mold 18 is completely separated from the master (which is only possible because of the previous “coating”-steps), first the surface of the mold 18, and especially the cavities 20, are coated by using one of the coating techniques described above (activating/passivating or polymer coating).
Then, vacuum is applied to the cavities 20 in order to remove air (not shown) from the cavities 20. This will help in the next step where the cavities 18 are filled with a second soft material 22. The second soft material 22 can be the same material as the first soft material 16. In general, it is also possible to use different materials, but it has been shown that soft materials for both the first and the second material 16, 22 are preferred since they guarantee the production of the delicate fibril structures. In the depicted case of
In order to fill the cavities 20 completely, it might be necessary to tilt the negative mold 18 to a certain degree in order to let the second material 22 flow into the cavities 20, especially also in the undercut parts of the cavities 20 (see e. g.
One possibility would also be to provide the second material 22 at the mold 18 without the second material initially coming into contact with the cavities 20, before applying the vacuum to the mold 18. Then, after the vacuum has been applied, the mold 18 is moved, in particular tilted, relative to the second material such that the second material 22 flows into and completely fills the degassed cavities 20 to form the fibrils within the cavities 20.
In this connection it is noted that it is either possible to only fill the cavities 20 with the second soft material 22 or to cover the whole mold 18 with the material 22 in order to also produce a backing substrate 24, which is supposed to hold the array of fibrils 10 together. Most common techniques use the latter and additionally provide a backing substrate 24. Said backing substrate 24 can have a thickness of about 0.5 to 1 mm or even up to 10 cm. For certain applications it is also possible to attach the backing substrate 24 to a payload (not shown) which is generally not limited in size and dimensions.
The next step comprises curing the second soft material 22 to finally produce the fibrils 10. Optionally some of the cured second material 22 can be kept in order to provide said backing substrate 24. In some embodiments it is also possible to remove the parts of the cured backing substrate 24 with a razor blade or the like (not shown) in order to only keep the fibrils 10 at the end. Furthermore, it could also be possible to first add the backing substrate 24 and then add a third material, which could act as a “new” backing substrate.
Lastly the first soft material 16 is demolded from the second soft material 22 to arrive at the produced fibrils 10 (and optionally the backing substrate 24).
Thus, in summary, for the elastomeric fibrillar surface fabrication, first a two-photon polymerization method was used to 3D-print the master 12 fibril array (patterns 14) on a glass substrate. The master 12 was then fluorosilanized to allow subsequent demolding, and PDMS (first soft material 16) was then cast and cured on the fibril array 14. The negative replica (negative mold 18) was peeled off and fluorosilanized and the subsequent molding 18 yielded a PDMS replica of the original master 12. A softer PDMS composition (20:1 monomer to crosslinker ratio) for the negative replica was employed to facilitate demolding of the final replica, for which standard PDMS (second soft material 22) was used (10:1 monomer to crosslinker ratio). The final replicas were characterized as prepared without any further processing.
Examples of the fibrils 10, which can be produced with the above disclosed method are shown in
The top surface of the tip 28 can be flat and smooth typically. Furthermore, the tip 28 can comprise a side wall with a given height and with a given sharpness and/or curvature at a tip edge that can be determined with the fabrication process resolution.
Said side wall can then comprise an almost 90 degree turn and an undercut cavity under the tip surface. The undercut area combines with the fibril tip 28 at the neck area of the fibril body 26 with a given curvature.
The method may further comprise the steps of:
The step of calculating an interfacial stress between the one or more fibrils and the boundary interface takes place by means of a computer implemented simulation. The interfacial stress is calculated directly by the Cauchy's equilibrium equations denoted above when the fibril 10 in the simulation is gradually stretched. When the maximum stress at the interface exceeds the interfacial strength between the fibril material and the object material, which is measured experimentally.
The present disclosure also relates to a computer implemented method of simulating an adhesive force of one or more fibrils to produce fibrils 10 with an increased adhesive force, the method comprising the steps of:
A fibril 10 as obtained by using a method as defined above and/or simulated on the basis of the methods described herein typically has at least one of an adhesive force selected in the range of 0.000001 mN to 1000 mN, a shape selected from the group of members comprising a profile composed of freeform curves, i.e. a fibril 10 having a base 32, a body 26 and a tip 28 that are generally axissysmmetric; and interfacial stress to a boundary surface selected in the range of 0.001 kPa to 10000 MPa and a material selected from the group of members consisting of organic, inorganic, polymers, rubbers, silicones, polyurethanes, biomaterials, biopolymers, composites, elastomers, liquid crystalline elastomers, thermoplastic elastomers, foams, fabric materials, particle materials, fibrous materials, or combinations of the foregoing of polymers.
The proposed computational single fibril design framework has two components: FEM-based adhesion mechanics simulation and Bayesian optimization (
In the simulation, the shape of the fibril 10 is constructed by design parameters with the following assumptions. First, the fibril 10 is assumed to be axisymmetric and the attached surface is assumed to be smooth and locally flat, as in the case of typical T- and mushroom-shaped fibril 10 designs. Additionally, the aspect ratio of the fibril is set to be 1:1 in order to compare the fibril 10 with other fibril designs at the same aspect ratio, which removes the effect of having different aspect ratios on adhesion tests.
Next, the side profile of the fibril 10 is constructed by a Bezier curve, which connects the top edge of the fibril 10 to the base of the fibril 10 smoothly with a given number of control points (
The effects of Bezier curve order are used to represent the fibril profile on adhesion (
Additionally, by predicting the maximum adhesion of MLFs 10 with higher order Bezier curves,
Finally, the simulation takes the fabrication constraints into account, where an infinite sharp edge is not possible to fabricate in practice. Therefore, the edge 30 of the tip 10 is rounded by a certain fillet radius, estimated from the experimental fabrication test results by the two-photon-polymerization-based direct 3D microprinting technique used in this study.
Once the shape of a fibril 10 is determined at a given iteration step, the adhesion of such a fibril 10 is simulated in a FEM simulation environment. Here, the simulation considers nonlinear deformation of a fibril during adhesion tests by incorporating the Mooney-Rivlin hyperelastic model, and a realistic detachment of the fibril 10 from a flat-punch probe using a Dugdale-Barenblatt cohesive zone model (
When a crack forms at the tip-substrate contact interface due to a normal stress exceeding the interfacial theoretical strength, the adhesion quickly reaches its maximum and decays subsequently. Here, the material properties and fabrication constraints determining the achievable edge sharpness are characterized experimentally and calibrated before the design of the fibrils 10.
The aforementioned simulation framework is directly linked to a Bayesian optimization framework. Bayesian optimization is a state-of-the-art global optimization method for problems where it is vital to optimize a performance criterion while keeping the number of evaluations of the system small. This method becomes an ideal candidate for minimizing the number of FEM-based adhesion evaluations without compromising the simulation accuracy. Here, the optimization goal is to maximize the adhesive force. The optimizer takes the Bezier-curve control points as optimizable variables and reaches to the optimal design in a given iteration number (
A previous study used a trained deep-neural-network (DNN), which simulates the adhesion profile at the interface. The difference lies in the simulation accuracy and data efficiency. First, the DNN-based approach does not compute the adhesion directly from the simulation, rather it estimates the adhesion indirectly with an infinitesimal fibril stretch. This approach would not be accurate when the fibril 10 is deformed due to the elasticity of the material, where the stress profile drastically changes.
Soft elastomers, such as polydimethylsiloxane (PDMS) used in our experimental fibrils 10 and in all previous T- and wedge-shaped mushroom fibril studies, do stretch significantly during the pull-off process resulting in hyperelastic regime. In contrast, the present method estimates adhesion directly right before the detachment considering the hyperelastic body deformation of the fibrils 10. Secondly, the present method is data efficient by orders of magnitude more than the DNN-based approach. More specifically, the present method allows to find an optimal shape within 150 function evaluations, whereas the DNN-based method requires more than 200,000 function evaluations to create a training dataset. Furthermore, the DNN-based approach is not flexible because a single change in the simulation might require another whole set of data causing a significant computational time.
Once the optimal fibril design is obtained, such fibril 10 is fabricated using two-photon polymerization (2PP) process and a subsequent double-molding-based replication technique (
First, computer-aided-design (CAD) files of the fibrils 10 are 3D-printed using a high-resolution rigid photoresist using 2PP to create a master template, and then fluorosilanized to allow easy replication. Afterwards, PDMS is casted on the master template to form a negative replica, cured, peeled off, and then fluorosilanized. Finally, the fluorosilanized PDMS negative mold 18 is replicated again to create a PDMS fibril 10 replica of the original master 12 (
The machine learning-based optimally designed fibrils (hereafter MLFs) with various tip diameters (30, 50, 70, 90 μm) are designed and fabricated (
A triplicate of each design was generated forming a tripod for accurate alignments with a smooth spherical glass probe (the diameter of the probe (8 mm) is much larger than the diameter of the fibril (30-90 μm); therefore, a flat-flat contact interface can be assumed locally), while keeping the number of fibrils minimum to reduce the effect of multiple fibrils on a soft backing layer (
The MLFs 10 result in unique stem 26 and tip 28 shapes compared to the previous fibrils. The shape of the MLFs 10 features a sharp tip edge 30 with a large wedge, a narrow neck, and a thickening body towards the base 32 of the fibril 10 (
Compared to the previously reported narrow neck designs to minimize any misalignment errors, these new designs utilize the stem deformation to alter stress profile at the interface along the stretch of the fibril 10 (see
In general, the predicted adhesion shows agreement with the experiments. In most of the cases the trend is captured, where the flat-punch fibrils 10 show the lowest and the MLFn 10 fibrils show the highest adhesion followed by T-shaped fibrils 10. Even the absolute adhesion values are predicted close to the experiments. The increasing discrepancy of the adhesion in the large fibrils 10 is considered to be attributed from the relatively increasing curvature of the measurement probe compared to the size of the tip of the fibrils 10, where flat-flat interfacial contact is compromised.
MLFn fibrils 10 show overall superior adhesion compared to the fibrils proposed in the prior art (
The present further finite element analysis (FEA) on the interfacial stress and fibril stretch (
In case of MLFn 10 (70 μm tip diameter), when a fibril 10 is insignificantly stretched (25% strain in
Fabrication constraints impact the optimal design and adhesion. The fabrication constraints worsen the adhesion performance due to the edge singularity induced by a rounded edge 30. The presented design methodology reflects such fabrication constraint, where it changes the optimal shape to cope with the edge singularity. To understand how a fabrication constraint changes the optimal design proposed by the algorithm, a set of design optimizations is performed. The edge fillet radius, re, (not the radius of the fibril) is chosen as a control variable, and the algorithm performs design optimizations for each re. Here the fibril diameter is set as 70 μm as a representative case.
The present simulation analysis suggests that the fabrication constraint greatly impacts the optimal shapes of the fibrils 10, showing that such constraints must be considered in the design processes. When such constraint is negligible compared to the size of the fibril 10, the optimal shape shows sharp edge with relatively thick neck and stem (
The design changes when the fabrication constraint becomes non-negligible. As the constraint increases, the fibril tends to show much narrower neck compared to the designs with less constraint. Also, the body shape changes by having a narrower body 26. These could be interpreted as to minimize deformation at the top body, by transferring most of the deformation to the lower body, which might reduce the risk of increasing stress peak at the tip-surface interface. Nonetheless, as the constraint increases, the adhesion deteriorates multifold due to the increased peak stress at the edge. Additionally, the analysis on the effect of the size of the fibril tip 28 with a given fabrication limit is shown, where such fabrication limit affects the optimal shape similar to above analysis (
The effect of the fibril body shape is further analyzed to understand its role in the fibril adhesion. It is found that the adhesion significantly depends on the fibril stem shape 26, whereas most of the previous studies focused only on the tip geometry. For example, small differences in the stem shapes of MLFv and MLFn (70 μm tip diameter,
The change in the radial control points shows high sensitivity in the adhesion performance. As the radial control points move towards the axis, the neck shape gets narrower resulting in the change of the curvature of the lower body, while the general wedge shape remains steady. This change of the stem shape is directly reflected in the adhesion, by having lower adhesion with a thinner stem (
The interfacial stress analysis of these variants (
The additional simulation study suggests that concave shape of the top surface resulted in a dramatic decrease of the adhesion, which is the worst-case scenario in dry adhesion optimization due to the increased edge singularity (
The study proposes a design framework of adhesive fibrils 10 based on Bayesian optimization and FEM-based adhesion mechanics simulations, its experimental validation, and analyses of the newly designed optimal shapes. While most of the previous studies focused on designs based on bioinspired T- and wedge-shaped tips and human intuition, this disclosure provides another perspective purely based on a machine learning algorithm. This new perspective discovered novel optimal designs and improved the adhesion performance significantly in the experiments. Compared to the previous optimal design frameworks, the proposed framework allows to incorporate full adhesion FEM simulations in the optimization loop for accurate results. The optimal designs computed by the proposed framework are validated by the simulations and experiments, and provide the evidences that 1) the large axial stem deformation causes the fibril tip 28 stress profile change, 2) fabrication constraints are important to be incorporated in the design process, and 3) fibril stem 26 shape is also important to the adhesion performance.
The design method is compared with the other approaches by allowing hyperelastic deformation of a fibril 10 in simulations. There are numerous studies to predict or to optimize adhesion using simulations. However, most of the studies assume that the adhesive fibril does not change shape during the pull-off or deforms very small staying in a linear deformation regime. However, soft elastomers, such as PDMS, stretches significantly during the pull-off process (more than 100% strain before detachment in our experiments, see
This disclosure also provides an insight to infer the role of shear stress on the normal adhesion. The trial by generating shear minimized MLFv reveals that indeed the shear stress may not affect the normal adhesion adversely, rather there is a positive correlation between them. Comparing the simulation and experimental results (
The FEM simulation is designed to estimate adhesion based on the design parameters (
Material characterization experiments are performed before the simulation. This is because soft elastomer (PDMS) fibril 10 might behave differently in small scales than in bulk in terms of mechanical properties. For this calibration, first, different shapes of fibrils 10 (flat punch, wedge-shaped, T-shaped in 30, 50, 70, 90 μm diameter) are fabricated using with 2PP and double-molding-based replication method (
A shape of a fibril 10 is determined by two sets of parameters: fixed parameter set and optimizable variable set. The fixed parameter set includes aspect ratio, tip diameter, edge fillet radius, the Mooney-Rivlin parameters and the theoretical interfacial strength. The optimizable parameter sets include two Bezier control points and the fibril base 32 diameter. The geometry of the fibril is axisymmetric. Once the shape is determined, a fibril 10 is meshed with more than 10,000 triangular elements, with extremely fine meshes towards the edge of the fibril to handle the interfacial stress singularity.
Bayesian optimization. Once the FEM simulation is constructed, the optimizable variable set, x, is directly handled by the Bayesian optimizer to produce the optimal adhesive fibril 10. Here, x is a 5-dimensional vector including two Bezier control points in r- and z-coordinate, and the diameter of the fibril base 32. The optimization problem is formulated as
where f(·) is the estimated adhesive force from the simulation, x* is the optimal variable set, p is the fixed parameter set, xmin is the lower bound of x, and xmax is the upper bound of x. This optimization problem is solved with a prebuilt Bayesian optimization function (bayesopt.m) in a proprietary programming language (MatLAB R2020a, The MathWorks, Inc.) with the expected improvement acquisition function. In this disclosure, the adhesion converges within 150 iterations. The communication between the optimizer and the simulator is established via LiveLink™ (COMSOL Multiphysics 5.4, COMSOL Inc.). The total of 150 iterations combined with optimization and simulation took less than 3 hours to produce the final design of the fibrils 10 using a performance desktop computer (20 cores of Intel Xeon CPU E5-2680 v2, 2.80 GHz (40 logical cores), 192 GB RAM, NVIDIA Quadro K5000 graphics card). The process of the Bayesian optimization and convergence is shown in
The fabrication process is similar to the fabrication method presented in PCT/EP2020/084461 whose disclosure regarding the precise steps of manufacture are hereby incorporated by reference, with modification in the fabrication parameters for sharper edge definitions. First, the fibrils 10 are designed using a CAD software (SolidWorks 2020, Dassault Systemes Co.), and then converted to a 2PP raster file (DeScribe, Nanoscribe GmbH), and printed using a 2PP 3D microprinter (Photonic Professional GT, Nanoscribe GmbH). The masters 12 are printed on an indium tin oxide (ITO)-coated glasses using a rigid IP-S commercial photoresist (Nanoscribe GmbH). To enhance the adhesion between the ITO-coated glass and the fibrils 10 to be printed, the ITO-coated glass is exposed to 1,1,1,3,3,3-hexamethyldisilazan (HMDS) in a desiccator more than 6 hours before the 3D print. Afterwards, the master 12 is developed in propylene glycol monomethyl ether acetate (PGMEA, Sigma-Aldrich Inc.) for 40 min and rinsed in isopropyl alcohol (IPA) for 3 mins. After the master fibril array (4×5) is fabricated, the array is exposed to ozone plasma (ozone plasma chamber, Zepto, Diener electronic GmbH) in 80° C. for 30 mins for surface activation, and subsequently placed in a vacuum desiccator with a glass vial containing 0.1 mL of trichloro (1H,1H,2H,2H-perfluoroctyl) silane (Sigma-Aldrich Inc.) for 35 mins. Finally, the master 12 is placed in 90° C. oven for 35 mins.
Degassed PDMS mixture (20:1 monomer to crosslinker ratio) is casted on top of the fluorosilanized master array and degassed for 30 mins in a vacuum desiccator. Then the master 12 and the cast are baked in a 90° C. oven for 35 mins. Afterwards, the negative mold 18 is peeled off. To silanize the negative mold 18, a similar procedure is applied. The negative mold is exposed to ozone plasma in 80° C. for 30 mins, and placed in a vacuum desiccator with trichlorosilane for 40 mins. Afterwards, the negative mold 18 is baked in 65° C. for 3 hours. After the negative mold 18 is peeled off, PDMS mixture (10:1 monomer to crosslinker ratio) is casted on top of the negative mold 18, degassed for 30 mins, and cured for 3 hours in a 65° C. preset oven. The positive replica is ready after peeled off from the negative mold 18. The surface roughness of the fibril surface is Rq of 0.42 μm with cutoff frequencies of λs=0.25 μm and λc=0.25 mm. Young's modulus of PDMS (65° C. cured for 3 hours) is 1.6 MPa. As the material behave different at the microscale, adhesion test of the T-shape fibril 10 (
Adhesion measurements. The adhesion of the fibrils 10 is tested in a customized adhesion setup. A load cell (GSO-25, Transducer Techniques LLC) is attached to a vertically moving (axial direction to the fibrils) high-precision motorized piezo stage (LPS-65 2″, Physik Instrumente GmbH) with 5 nm positioning resolution. Horizontal positions (xy-position) are adjusted by manual xy-stage with an additional motorized piezo stage. The motion control of the motorized stages and the force measurement from the loadcell are performed simultaneously in a custom-built program in LabVIEW environment (National Instruments Co.) with 1 KHz data acquisition rate. The raw data is filtered to reduce noise with 10 points median filter. An inverted optical microscope (Axio Observer A1, Carl Zeiss AG) is placed under the setup to visualize the contact interface for alignments and adhesion analyses.
A spherical glass probe with a diameter of 4 mm is used as a probe directly connected to the load cell. The surface roughness of the glass probe is Rq of 1.04 μm with cutoff frequencies of λs=0.25 μm and λc=0.25 mm. The xy-position of the probe is carefully adjusted to place the probe at the center of the triplicate of the fibrils 10. The approach and retraction speeds are 10% of the fibril size per second (e.g., 7 μm s−1 for 70 μm fibril) to minimize the viscoelastic effects. The probe indents until it reaches a certain preload (150 kPa for all sizes; e.g., 0.577 mN for 70 μm tip diameter fibril 10), waits for 30 seconds for the relaxation, and is retracted (
Number | Date | Country | Kind |
---|---|---|---|
21162253.5 | Mar 2021 | EP | regional |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2022/056341 | 3/11/2022 | WO |