CN110717297A - Method and system for establishing fabric simulation model, storage medium and electronic equipment - Google Patents

Method and system for establishing fabric simulation model, storage medium and electronic equipment Download PDF

Info

Publication number
CN110717297A
CN110717297A CN201910969172.2A CN201910969172A CN110717297A CN 110717297 A CN110717297 A CN 110717297A CN 201910969172 A CN201910969172 A CN 201910969172A CN 110717297 A CN110717297 A CN 110717297A
Authority
CN
China
Prior art keywords
tensor
initial
fabric
calculating
moment
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.)
Granted
Application number
CN201910969172.2A
Other languages
Chinese (zh)
Other versions
CN110717297B (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.)
Wuhan Textile University
Original Assignee
Wuhan Textile University
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 Wuhan Textile University filed Critical Wuhan Textile University
Priority to CN201910969172.2A priority Critical patent/CN110717297B/en
Publication of CN110717297A publication Critical patent/CN110717297A/en
Application granted granted Critical
Publication of CN110717297B publication Critical patent/CN110717297B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides a method and a system for establishing a fabric simulation model, a storage medium and electronic equipment, wherein the method comprises the following steps: obtaining a fabric curved surface at an initial moment, and dispersing the fabric curved surface into quadrilateral mesh surface patches with different sizes according to different material attributes; analyzing initial coordinates of the gravity center of each quadrilateral mesh surface patch; calculating an initial Riemann degree tensor and an initial Riemann curvature tensor according to the initial coordinates; taking the initial moment as a start, calculating internal deformation energy according to the Riemann scale tensor and the Riemann curvature tensor at the current moment in sequence, and calculating the coordinate position of the gravity center at the next moment by combining stress analysis; and establishing a fabric simulation motion model according to the coordinate position. The invention researches the fabric in the curved Riemannian manifold space, and the invention can improve the existing fabric modeling theory and method by adopting the non-European geometric theory, and makes up the defects of complex fabric modeling, slow convergence speed, poor stability, poor simulation effect and the like in the European space.

Description

Method and system for establishing fabric simulation model, storage medium and electronic equipment
Technical Field
The invention relates to the technical field of fabric simulation, in particular to a method and a system for establishing a fabric simulation model, a storage medium and electronic equipment.
Background
Current methods for studying fabric models and simulations are generally divided into three categories: geometric methods, physical methods and hybrid methods. The fabric simulation technology based on geometric modeling only focuses on the external expression form of the fabric, and is more like an advanced drawing tool, and the simulation speed is high; research on methods based on physical modeling has been dominant in recent years, and can be divided into the following according to the way of simulating fabric structures: discrete models including spring-mass models, mesh models, particle models; continuous models, including elastic deformation models, aerodynamic models, finite element models, and the like; the method based on the hybrid modeling adopts a geometric modeling method to model, considers the physical properties of the fabric in the simulation process, and considers both the simulation speed and the reality sense of the simulation, but the simulation effect is still unsatisfactory. The existing research depends on a fabric model established in a Euclidean space, namely improvements on simulation models or methods of different fabrics in the Euclidean space are made from the aspects of response speed or effect of simulation, but the problems of complex modeling, large calculation amount, difficult simulation effect and the like still exist.
Disclosure of Invention
The invention aims to provide a method and a system for establishing a fabric simulation model, a storage medium and electronic equipment, so that a fabric can be put in a curved Riemannian manifold space for research, the research result can improve the conventional fabric modeling theory and simulation method, and the defects of complexity, low convergence speed, poor stability, poor simulation effect and the like of fabric modeling in an Euclidean space are overcome.
The technical scheme provided by the invention is as follows:
the invention provides a method for establishing a fabric simulation model, which comprises the following steps:
acquiring a fabric curved surface at an initial moment in a Riemann manifold space, dispersing the fabric curved surface into a plurality of quadrilateral mesh surface patches, and equivalently using all the quadrilateral mesh surface patches as a dispersion unit to perform analysis processing;
analyzing the gravity center of each quadrilateral mesh surface patch to obtain the initial coordinate of the gravity center;
calculating an initial Riemannian degree tensor and an initial Riemannian curvature tensor of each quadrilateral grid surface patch at an initial moment according to the initial coordinates;
calculating the internal deformation energy of each quadrilateral mesh patch at the current moment sequentially according to the Riemannian scale tensor at the current moment and the Riemannian curvature tensor at the current moment by taking the initial moment as a start, and calculating the coordinate position of the gravity center of each quadrilateral mesh patch at the next moment by combining the stress analysis of the quadrilateral mesh patch at the current moment and the internal deformation energy, wherein the difference between the next moment and the current moment is a preset time step length;
and establishing a fabric simulation motion model according to the coordinate position.
Further, calculating the internal deformation energy of each quadrilateral mesh patch at the current time according to the riemann scale tensor at the current time and the riemann curvature tensor at the current time, and calculating the coordinate position of the gravity center of each quadrilateral mesh patch at the next time by combining the stress analysis of the quadrilateral mesh patch at the current time and the internal deformation energy specifically includes:
acquiring an instantaneous vector of the gravity center of any quadrilateral mesh patch at the current moment;
calculating the Riemannian degree tensor of any quadrilateral mesh patch at the current moment and the Riemannian curvature tensor of the current moment according to the instantaneous vector;
calculating the internal deformation energy of each quadrilateral mesh patch at the current time according to the initial Riemannian degree tensor, the initial Riemannian curvature tensor, the Riemannian degree tensor at the current time and the Riemannian curvature tensor at the current time;
and calculating the coordinate position of the next moment by combining the stress analysis of the quadrilateral mesh surface patch and the internal deformation energy.
Further, calculating the internal deformation energy of each quadrilateral mesh patch at the current time according to the riemann measure tensor at the current time and the riemann curvature tensor at the current time specifically includes:
defining a region of the any one quadrilateral mesh patch as Ω, Ω { (l)1,l2)∈R2Define the instantaneous vector of the center of gravity as
Figure BDA0002231512640000031
Calculating the stretching and shearing energy at the current moment according to the initial Riemannian degree quantity and the Riemannian degree quantity at the current moment,
Figure BDA0002231512640000032
wherein the content of the first and second substances,for the stretch-shear energy at the present moment, GtIs the Riemann scale number, G, of the current time0Is the initial number of the Riemann scale,
Figure BDA0002231512640000034
Ewpand EwtRespectively the Young modulus of the fabric in the warp and weft directions, H is the rigidity modulus,
Figure BDA00022315126400000317
is a matrix wsA Hilbert-Schmitt norm of the weighted measure tensor variabilities;
calculating the bending energy of the current moment according to the initial Riemann curvature tensor and the Riemann curvature tensor of the current moment,
Figure BDA0002231512640000035
wherein the content of the first and second substances,
Figure BDA0002231512640000036
as the bending energy at the present time is,Rtis the curvature tensor of the current time, R0Is the initial riemann curvature tensor,
Figure BDA0002231512640000037
Rfwpand RfwtBending stiffness of the fabric in the warp and weft directions respectively,
Figure BDA00022315126400000318
is a matrix wbA Hilbert-Schmidt norm of the weighted curvature tensor variations;
calculating the internal deformation energy from the tensile shear energy and the bending energy
Figure BDA0002231512640000038
Further, calculating the coordinate position of the next moment by combining the stress analysis of the quadrilateral mesh patch and the internal deformation energy specifically comprises:
calculating the elastic internal force of the quadrilateral mesh patch according to the internal deformation energy
Figure BDA00022315126400000310
Calculating inertial internal force
Figure BDA00022315126400000311
Viscous damping force
Figure BDA00022315126400000312
Weight of fabric
Figure BDA00022315126400000313
And air resistanceWherein the content of the first and second substances,
Figure BDA00022315126400000315
p is the density of the fabric,mu is a damping coefficient;
according to the first dynamic principle, a dynamic equation reflecting the motion of the fabric is obtained,
and calculating the coordinate position of the next moment according to the dynamic equation.
The invention also provides a system for establishing the fabric simulation model, which comprises the following steps:
the dispersion module is used for obtaining a fabric curved surface at an initial moment in a Riemannian manifold space, dispersing the fabric curved surface into a plurality of quadrilateral mesh surface patches, and equivalently using all the quadrilateral mesh surface patches as dispersion units for analysis processing;
the gravity center calculation module analyzes the gravity center of each quadrilateral mesh surface patch obtained by the dispersion module and obtains the initial coordinate of the gravity center;
the parameter calculation module is used for calculating an initial Riemannian degree tensor and an initial Riemannian curvature tensor of each quadrilateral mesh surface patch at an initial moment according to the initial coordinates obtained by the gravity center calculation module;
the coordinate calculation module is used for calculating the internal deformation energy of each quadrilateral mesh patch at the current moment by taking the initial moment as a start according to the Riemannian scale tensor at the current moment and the Riemannian curvature tensor at the current moment which are obtained by the parameter calculation module in sequence, and calculating the coordinate position of the gravity center of each quadrilateral mesh patch at the next moment by combining the stress analysis of the quadrilateral mesh patch at the current moment and the internal deformation energy, wherein the difference between the next moment and the current moment is a preset time step;
and the model establishing module is used for establishing a fabric simulation motion model according to the coordinate position obtained by the coordinate calculating module.
Further, the coordinate calculation module specifically includes:
the initial definition unit is used for acquiring the instantaneous vector of the gravity center of any quadrilateral mesh patch at the current moment;
the tensor calculation unit is used for calculating the Riemannian degree tensor of any quadrilateral mesh patch at the current moment and the Riemannian curvature tensor at the current moment according to the instantaneous vector defined by the primary definition unit;
a parameter calculation unit which calculates internal deformation energy of each quadrilateral mesh patch at the current time according to the initial Riemannian degree tensor, the initial Riemannian curvature tensor, the Riemannian degree tensor at the current time and the Riemannian curvature tensor obtained by the tensor calculation unit;
and the coordinate calculation unit is used for calculating the coordinate position of the next moment by combining the stress analysis of the quadrilateral mesh surface patch and the internal deformation energy obtained by the parameter calculation unit.
Further, the parameter calculating unit specifically includes:
defining a sub-unit, defining the region of any quadrilateral mesh patch as omega, omega { (l)1,l2)∈R2Define the instantaneous vector of the center of gravity as
Figure BDA00022315126400000510
A parameter calculating subunit that calculates the stretch-shear energy at the current time from the initial riemannian tensor obtained by the tensor calculating unit, the riemannian tensor at the current time, and the parameters defined by the defining subunit,
Figure BDA0002231512640000051
wherein the content of the first and second substances,
Figure BDA0002231512640000052
for the stretch-shear energy at the present moment, GtIs the Riemann scale number, G, of the current time0Is the initial number of the Riemann scale,
Figure BDA0002231512640000053
Ewpand EwtRespectively the Young modulus of the fabric in the warp and weft directions, H is the rigidity modulus,
Figure BDA0002231512640000054
is a matrix wsA Hilbert-Schmitt norm of the weighted measure tensor variabilities;
the parameter calculating subunit calculates the bending energy at the current time according to the initial Riemann curvature tensor obtained by the tensor calculating unit, the Riemann curvature tensor at the current time and the parameters defined by the defining subunit,
Figure BDA0002231512640000055
wherein the content of the first and second substances,
Figure BDA0002231512640000056
is the bending energy at the present moment, RtIs the curvature tensor of the current time, R0Is the initial riemann curvature tensor,
Figure BDA0002231512640000057
Rfwpand RfwtBending stiffness of the fabric in the warp and weft directions respectively,
Figure BDA0002231512640000058
is a matrix wbA Hilbert-Schmidt norm of the weighted curvature tensor variations;
the parameter calculation subunit calculates the internal deformation energy from the tensile shear energy and the bending energy
Further, the coordinate calculation unit specifically includes:
a stress calculating subunit for calculating the elastic internal force of the quadrilateral mesh patch according to the internal deformation energy
Figure BDA0002231512640000061
The stress calculation subunit calculates the inertial internal force
Figure BDA0002231512640000062
Viscous damping force
Figure BDA0002231512640000063
Weight of fabric
Figure BDA0002231512640000064
And air resistance
Figure BDA0002231512640000065
Wherein the content of the first and second substances,
Figure BDA0002231512640000066
p is the density of the fabric,
Figure BDA0002231512640000067
mu is a damping coefficient;
a stress analysis subunit for obtaining a dynamic equation reflecting the motion of the fabric according to the first dynamic principle and the force calculated by the stress calculation subunit,
Figure BDA0002231512640000068
and the coordinate calculation subunit calculates the coordinate position of the next moment according to the dynamic equation obtained by the stress analysis subunit.
The invention also provides a storage medium having stored thereon a computer program which, when executed by a processor, implements any of the methods described above.
The invention also provides an electronic device comprising a memory and a processor, wherein the memory stores a computer program running on the processor, and the processor implements any one of the methods described above when executing the computer program.
By the method and the system for establishing the fabric simulation model, the storage medium and the electronic equipment, the fabric is put in a curved Riemann manifold space for research, and the modeling and simulation problems of the fabric are solved by a manifold surface theory in the space. The method can make up the defects of complex fabric modeling, low convergence speed, poor stability, poor simulation effect and the like in Euclidean space, and can improve the simulation time performance while ensuring the vivid fabric simulation effect. The method not only theoretically improves the fabric modeling and simulation theory, but also has wide industrial application prospect, and the result can be directly applied to the application fields of virtual garment animation simulation, movie and television advertisement, game and entertainment and the like, and especially can generate great promotion effect on the development of the current virtual garment design, production and sale technology.
Drawings
The above features, technical features, advantages and implementations of the method, system, storage medium and electronic device for building a fabric simulation model will be further described in the following detailed description of preferred embodiments in a clearly understandable manner with reference to the accompanying drawings.
FIG. 1 is a flow diagram of one embodiment of a method of building a fabric simulation model of the present invention;
FIG. 2 is a flow diagram of another embodiment of a method of building a fabric simulation model of the present invention;
fig. 3 is a schematic structural diagram of an embodiment of the system for building a fabric simulation model according to the present invention.
Detailed Description
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following description will explain specific embodiments of the present invention with reference to the drawings of the specification. It is obvious that the drawings in the following description are only some examples of the invention, from which other drawings and embodiments can be derived by a person skilled in the art without inventive effort.
For the sake of simplicity, only the parts relevant to the present invention are schematically shown in the drawings, and they do not represent the actual structure as a product. In addition, in order to make the drawings concise and understandable, components having the same structure or function in some of the drawings are only schematically illustrated or only labeled. In this document, "one" means not only "only one" but also a case of "more than one".
In the embodiment of the present invention, as shown in fig. 1, a method for building a fabric simulation model includes:
s100, obtaining a fabric curved surface at an initial moment in a Riemannian manifold space, dispersing the fabric curved surface into a plurality of quadrilateral mesh surface patches, and performing analysis processing by equivalently using all the quadrilateral mesh surface patches as discrete units;
s200, analyzing the gravity center of each quadrilateral mesh surface patch to obtain the initial coordinate of the gravity center;
s300, calculating an initial Riemannian degree tensor and an initial Riemannian curvature tensor of each quadrilateral grid surface patch at an initial moment according to the initial coordinates;
s400, starting from an initial moment, sequentially calculating internal deformation energy of each quadrilateral mesh patch at the current moment according to the Riemannian degree tensor at the current moment and the Riemannian curvature tensor at the current moment, and calculating the coordinate position of the gravity center of each quadrilateral mesh patch at the next moment by combining the stress analysis of the quadrilateral mesh patch at the current moment and the internal deformation energy, wherein the difference between the next moment and the current moment is a preset time step length;
s500, establishing a fabric simulation motion model according to the coordinate position.
Specifically, in the present embodiment, the riemann manifold is a differential manifold having the riemann tensor, has a simple structure, has the maximum symmetry and isotropic properties, and is suitable for studying the deformed form of the flexible fabric, the structure of which is also isotropic. The elastic deformation of the curved surface of the fabric is related to shearing, bending, and the like occurring inside the fabric, and the deformation of the curved surface of the fabric is expressed by the riemann tensor and the riemann curvature tensor. The fabric dynamic simulation based on the Riemann manifold aims to research the material attribute of the fabric and the deformation rule of the curved surface of the fabric on the basis of the Riemann space manifold theory.
Firstly, initializing a series of parameters to obtain a fabric curved surface at an initial moment, and dispersing the fabric curved surface in a Riemannian manifold space into a plurality of quadrilateral mesh surface patches according to different material properties, wherein the quadrilateral mesh surface patches can be regular quadrilaterals, such as rectangles, rhombuses, parallelograms and the like, or irregular quadrilaterals, and the four sides are not equal. The quadrilateral mesh patches with the same material property have the same size, and the quadrilateral mesh patches with different material properties can be the same or similar or completely different. In order to increase the calculation speed, the fabric curved surface is usually discretized into the same quadrilateral mesh patches in the practical application process, but since the fabric curved surface itself may be irregular, the fabric curved surface is discretized into the same or similar quadrilateral mesh patches as much as possible.
Since each parameter, for example, coordinates, of the fabric curved surface at the initial time are known, the center of gravity of each quadrilateral mesh patch is analyzed, initial coordinates of the center of gravity are obtained, and an initial riemann tensor and an initial riemann curvature tensor of each quadrilateral mesh patch at the initial time are calculated according to the initial coordinates. Relative positions of the centers of gravity of the quadrilateral mesh patches are unchanged, but due to deformation of the curved surface of the fabric, coordinates of the centers of gravity of the quadrilateral mesh patches in a world coordinate system or other coordinate systems are changed, so that a fabric simulation motion model can be established by determining the coordinate positions of the centers of gravity of the quadrilateral mesh patches at different moments.
And then, starting from the initial time, sequentially calculating the internal deformation energy of each quadrilateral mesh patch at the current time according to the Riemannian degree tensor at the current time and the Riemannian curvature tensor at the current time, and calculating the coordinate position of the gravity center of each quadrilateral mesh patch at the next time by combining the stress analysis of the quadrilateral mesh patch at the current time and the internal deformation energy, wherein the difference between the next time and the current time is a preset time step length.
The method comprises the steps that the initial moment is the current moment, the coordinate position of the gravity center of each quadrilateral mesh patch at the next moment is calculated by combining stress analysis of the quadrilateral mesh patches according to the initial Riemannian degree tensor and the initial Riemannian curvature tensor, and the fabric curved surface moves according to the coordinate position of the gravity center of each quadrilateral mesh patch at the next moment to obtain a simulation model. And then the next moment is a new current moment, and the coordinate position of the gravity center of each quadrilateral mesh patch at the new next moment is calculated. And by analogy, calculating the coordinate position of the gravity center of each quadrilateral mesh patch at the new next moment at one time, and then adjusting the curved surface of the fabric to finally obtain the fabric simulated motion model.
The invention researches the fabric in a curved Riemannian manifold space, solves the modeling and simulation problems of the fabric by using a manifold surface theory in the space, ensures the simulation effect of the fabric to be vivid and can also improve the simulation time performance.
Another embodiment of the present invention is an optimized embodiment of the foregoing embodiment, and as shown in fig. 2, the main improvement of this embodiment is that, compared with the foregoing embodiment, the S400 calculates internal deformation energy of each quadrilateral mesh patch at the present time according to the riemann scale tensor at the present time and the riemann curvature tensor at the present time, and the calculating, by combining the force analysis of the quadrilateral mesh patch at the present time and the internal deformation energy, a coordinate position of a center of gravity of each quadrilateral mesh patch at the next time specifically includes:
s410, acquiring an instantaneous vector of the gravity center of any quadrilateral mesh patch at the current moment;
s420, calculating a Riemannian degree tensor of the any one quadrilateral mesh patch at the current moment and a Riemannian curvature tensor of the any one quadrilateral mesh patch at the current moment according to the instantaneous vector;
s430, calculating internal deformation energy of each quadrilateral mesh patch at the current time according to the initial Riemannian degree tensor, the initial Riemannian curvature tensor, the Riemannian degree tensor at the current time and the Riemannian curvature tensor at the current time;
and S440, calculating the coordinate position of the next moment by combining the stress analysis of the quadrilateral mesh patch and the internal deformation energy.
Specifically, in this embodiment, the curved surface of the fabric in the riemann manifold space is discretized into a plurality of quadrilateral mesh patches, and an initial riemann measure tensor and an initial riemann curvature tensor of each quadrilateral mesh patch at an initial time are calculated. Any quadrilateral mesh patch is selected to calculate the coordinate position of the gravity center of the quadrilateral mesh patch at the current moment, the calculation methods of the coordinate positions of the gravity centers of the quadrilateral mesh patches at the current moment are the same, and one of the quadrilateral mesh patches is taken as an example for description in the embodiment.
And acquiring an instantaneous vector of the gravity center of any quadrilateral mesh patch at the current moment, wherein the instantaneous vector is obtained by calculating the coordinate of the gravity center at the current moment according to the related parameters of the previous moment. And calculating the Riemannian scale tensor of the corresponding quadrilateral mesh patch at the current moment and the Riemannian curvature tensor of the current moment according to the instantaneous vector, and calculating the internal deformation energy of the quadrilateral mesh patch corresponding to the current moment by combining the initial Riemannian scale tensor, the initial Riemannian curvature tensor, the Riemannian scale tensor at the current moment and the Riemannian curvature tensor at the current moment. And then establishing a dynamic equation reflecting the motion of the fabric based on the stress analysis of the quadrilateral mesh surface patches, and solving by using the quadrilateral mesh surface patches as finite element units to obtain the instantaneous positions of the gravity centers of all the quadrilateral mesh surface patches.
Further, the step S430 of calculating the internal deformation energy of each quadrilateral mesh patch at the current time according to the initial riemann scale tensor, the initial riemann curvature tensor, the riemann scale tensor at the current time, and the riemann curvature tensor at the current time specifically includes: s431 defines an area of the any one quadrilateral mesh patch as Ω, Ω { (l)1,l2)∈R2Define the instantaneous vector of the center of gravity as
Figure BDA0002231512640000103
S432 calculates the stretch-shear energy at the current time from the initial riemann gauge number and the riemann gauge number at the current time,
Figure BDA0002231512640000101
wherein the content of the first and second substances,
Figure BDA0002231512640000102
for the stretch-shear energy at the present moment, GtIs the Riemann scale number, G, of the current time0Is the initial number of the Riemann scale,
Figure BDA0002231512640000111
Ewpand EwtRespectively the Young modulus of the fabric in the warp and weft directions, H is the rigidity modulus,
Figure BDA0002231512640000112
is a matrix wsA Hilbert-Schmitt norm of the weighted measure tensor variabilities; s433 calculates a bending energy at the current time from the initial riemann curvature tensor and the riemann curvature tensor at the current time,
Figure BDA0002231512640000113
wherein the content of the first and second substances,
Figure BDA0002231512640000114
is the bending energy at the present moment, RtIs the curvature tensor of the current time, R0Is the initial riemann curvature tensor,
Figure BDA0002231512640000115
Rfwpand RfwtBending stiffness of the fabric in the warp and weft directions respectively,
Figure BDA0002231512640000116
is a matrix wbA Hilbert-Schmidt norm of the weighted curvature tensor variations; s434 calculating the internal deformation energy from the tensile shear energy and the bending energy
Figure BDA0002231512640000117
Specifically, the region of any one quadrilateral mesh patch is defined as Ω, Ω { (l)1,l2)∈R2In which l1,l2Parameterized respectively for the warp and weft edges of a quadrilateral mesh patch, R2Representing that the quadrilateral mesh patch belongs to a two-dimensional space, and simultaneously defining the instantaneous vector of the gravity center of the quadrilateral mesh patch at the current moment as
Figure BDA0002231512640000118
Differential manifold MnThe higher Riemann scale tensor g is an inner product given on a tangent space at a point p on M<,>pAnd the inner product is smooth. The fabric surface is considered as a differential manifold, then the riemann gauge tensor can be computed as:
Figure BDA0002231512640000119
riemann curvature tensor: r ═ RijklI, j, k, l ═ 1,2 …, n, the riemann curvature tensor RijklIs a 4 th order tensor, 4 indices can all vary from 1 to n. And calculating the stretching shearing energy at the current moment according to the initial Riemannian degree tensor and the Riemannian degree tensor at the current moment, calculating the bending energy at the current moment according to the initial Riemannian curvature tensor and the Riemannian curvature tensor at the current moment, and finally obtaining the internal deformation energy at the current moment.
Further, the step S440 of calculating the coordinate position of the next moment by combining the force analysis of the quadrilateral mesh patch and the internal deformation energy specifically includes: s441, according to the internal deformation energy, the elastic internal force of the quadrilateral mesh surface patch is calculated
Figure BDA0002231512640000121
S442 calculating the inertial internal force
Figure BDA0002231512640000122
Viscous damping force
Figure BDA0002231512640000123
Weight of fabric
Figure BDA0002231512640000124
And air resistanceWherein the content of the first and second substances,
Figure BDA0002231512640000126
p is the density of the fabric,
Figure BDA0002231512640000127
mu is a damping coefficient; s443, according to the first dynamic principle, obtaining a dynamic equation reflecting the movement of the fabric,
Figure BDA0002231512640000128
s444 calculates a coordinate position of the next time according to the dynamic equation.
Specifically, the elastic deformation of the curved surface of the fabric is related to shearing, bending and the like generated in the fabric, and the elastic internal force of the quadrilateral mesh surface patch is calculated through the internal deformation energy
Figure BDA0002231512640000129
Calculating inertial internal force
Figure BDA00022315126400001210
Viscous damping force
Figure BDA00022315126400001211
Weight of fabric
Figure BDA00022315126400001212
And air resistance
Figure BDA00022315126400001213
The fabric is stressed by internal force and external force, wherein the internal force comprises inertial internal force, viscous damping force and elastic internal force, and the external force comprises gravity and air resistance. Thus, according to the first kinetic principle, a dynamic equation reflecting the motion of the fabric is obtained. Solving partial differential equations by adopting a finite element method, taking quadrilateral patches of the fabric as finite element units, and solving the partial differential equations (dynamic equations) to obtain the instantaneous positions of the gravity centers of all quadrilateral mesh patches so as to determine the coordinate positions of the gravity centers.
All or part of the flow of the method according to the embodiments of the present invention may be implemented by a computer program, which may be stored in a computer-readable storage medium and used by a processor to implement the steps of the embodiments of the method. Wherein the computer program comprises computer program code, which may be in the form of source code, object code, an executable file or some intermediate form, etc. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, usb disk, removable hard disk, magnetic disk, optical disk, computer Memory, Read-only Memory (ROM), Random Access Memory (RAM), electrical carrier wave signals, telecommunications signals, software distribution medium, etc. It should be noted that the computer readable medium may contain content that is subject to appropriate increase or decrease as required by legislation and patent practice in jurisdictions, for example, in some jurisdictions, computer readable media does not include electrical carrier signals and telecommunications signals as is required by legislation and patent practice.
An embodiment of the present invention further provides an electronic device, which includes a memory and a processor, wherein the memory stores a computer program running on the processor, and the processor executes the computer program to implement all or part of the method steps in the first embodiment.
The Processor may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), an off-the-shelf Programmable Gate Array (FPGA) or other Programmable logic device, discrete Gate or transistor logic, discrete hardware components, etc. The general purpose processor may be a microprocessor or the processor may be any conventional processor or the like which is the control center for the computer device and which connects the various parts of the overall computer device using various interfaces and lines.
The memory may be used to store the computer programs and/or modules, and the processor may implement various functions of the computer device by running or executing the computer programs and/or modules stored in the memory and invoking data stored in the memory. The memory may mainly include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required by at least one function (such as a sound playing function, an image playing function, etc.), and the like; the storage data area may store data (such as audio data, video data, etc.) created according to the use of the cellular phone, etc. In addition, the memory may include high speed random access memory, and may also include non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), at least one magnetic disk storage device, a Flash memory device, or other volatile solid state storage device.
It should be noted that the above embodiments can be freely combined as necessary. The foregoing is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, various modifications and decorations can be made without departing from the principle of the present invention, and these modifications and decorations should also be regarded as the protection scope of the present invention.

Claims (10)

1. The method for establishing the fabric simulation model is characterized by comprising the following steps:
acquiring a fabric curved surface at an initial moment in a Riemann manifold space, dispersing the fabric curved surface into a plurality of quadrilateral mesh surface patches, and equivalently using all the quadrilateral mesh surface patches as a dispersion unit to perform analysis processing;
analyzing the gravity center of each quadrilateral mesh surface patch to obtain the initial coordinate of the gravity center;
calculating an initial Riemannian degree tensor and an initial Riemannian curvature tensor of each quadrilateral grid surface patch at an initial moment according to the initial coordinates;
calculating the internal deformation energy of each quadrilateral mesh patch at the current moment sequentially according to the Riemannian scale tensor at the current moment and the Riemannian curvature tensor at the current moment by taking the initial moment as a start, and calculating the coordinate position of the gravity center of each quadrilateral mesh patch at the next moment by combining the stress analysis of the quadrilateral mesh patch at the current moment and the internal deformation energy, wherein the difference between the next moment and the current moment is a preset time step length;
and establishing a fabric simulation motion model according to the coordinate position.
2. The method for building a fabric simulation model according to claim 1, wherein the step of calculating the internal deformation energy of each quadrilateral mesh patch at the current time according to the riemann measure tensor at the current time and the riemann curvature tensor at the current time, and the step of calculating the coordinate position of the center of gravity of each quadrilateral mesh patch at the next time by combining the stress analysis of the quadrilateral mesh patch at the current time and the internal deformation energy specifically comprises the steps of:
acquiring an instantaneous vector of the gravity center of any quadrilateral mesh patch at the current moment;
calculating the Riemannian degree tensor of any quadrilateral mesh patch at the current moment and the Riemannian curvature tensor of the current moment according to the instantaneous vector;
calculating the internal deformation energy of each quadrilateral mesh patch at the current time according to the initial Riemannian degree tensor, the initial Riemannian curvature tensor, the Riemannian degree tensor at the current time and the Riemannian curvature tensor at the current time;
and calculating the coordinate position of the next moment by combining the stress analysis of the quadrilateral mesh surface patch and the internal deformation energy.
3. The method for building a fabric simulation model according to claim 2, wherein the calculating the internal deformation energy of each quadrilateral mesh patch at the current time according to the riemann measure tensor at the current time and the riemann curvature tensor at the current time specifically includes:
defining a region of the any one quadrilateral mesh patch as Ω, Ω { (l)1,l2)∈R2Define the instantaneous vector of the center of gravity as
Figure FDA0002231512630000021
Calculating the stretching and shearing energy at the current moment according to the initial Riemannian degree quantity and the Riemannian degree quantity at the current moment,
Figure FDA0002231512630000022
wherein the content of the first and second substances,
Figure FDA0002231512630000023
for the stretch-shear energy at the present moment, GtIs the Riemann scale number, G, of the current time0Is the initial number of the Riemann scale,
Figure FDA0002231512630000024
Ewpand EwtRespectively the Young modulus of the fabric in the warp and weft directions, H is the rigidity modulus,
Figure FDA0002231512630000025
is a matrix wsA Hilbert-Schmitt norm of the weighted measure tensor variabilities;
calculating the bending energy of the current moment according to the initial Riemann curvature tensor and the Riemann curvature tensor of the current moment,wherein the content of the first and second substances,is the bending energy at the present moment, RtIs the curvature tensor of the current time, R0Is the initial riemann curvature tensor,
Figure FDA0002231512630000028
Rfwpand RfwtBending stiffness of the fabric in the warp and weft directions respectively,is a matrix wbWeighted curvature tensor variationalThe hilbert-schmitt norm of (d);
calculating the internal deformation energy from the tensile shear energy and the bending energy
Figure FDA00022315126300000210
Figure FDA00022315126300000211
4. The method for building the fabric simulation model according to claim 3, wherein the step of calculating the coordinate position of the next moment by combining the force analysis of the quadrilateral mesh patch and the internal deformation energy specifically comprises the following steps:
calculating the elastic internal force of the quadrilateral mesh patch according to the internal deformation energy
Figure FDA0002231512630000031
Figure FDA0002231512630000032
Calculating inertial internal force
Figure FDA0002231512630000033
Viscous damping force
Figure FDA0002231512630000034
Weight of fabric
Figure FDA0002231512630000035
And air resistance
Figure FDA0002231512630000036
Wherein the content of the first and second substances,
Figure FDA0002231512630000037
p is the density of the fabric,
Figure FDA0002231512630000038
mu is a damping coefficient;
according to the first dynamic principle, a dynamic equation reflecting the motion of the fabric is obtained,
and calculating the coordinate position of the next moment according to the dynamic equation.
5. The system for establishing the fabric simulation model is characterized by comprising the following steps:
the dispersion module is used for obtaining a fabric curved surface at an initial moment in a Riemannian manifold space, dispersing the fabric curved surface into a plurality of quadrilateral mesh surface patches, and equivalently using all the quadrilateral mesh surface patches as dispersion units for analysis processing;
the gravity center calculation module analyzes the gravity center of each quadrilateral mesh surface patch obtained by the dispersion module and obtains the initial coordinate of the gravity center;
the parameter calculation module is used for calculating an initial Riemannian degree tensor and an initial Riemannian curvature tensor of each quadrilateral mesh surface patch at an initial moment according to the initial coordinates obtained by the gravity center calculation module;
the coordinate calculation module is used for calculating the internal deformation energy of each quadrilateral mesh patch at the current moment by taking the initial moment as a start according to the Riemannian scale tensor at the current moment and the Riemannian curvature tensor at the current moment which are obtained by the parameter calculation module in sequence, and calculating the coordinate position of the gravity center of each quadrilateral mesh patch at the next moment by combining the stress analysis of the quadrilateral mesh patch at the current moment and the internal deformation energy, wherein the difference between the next moment and the current moment is a preset time step;
and the model establishing module is used for establishing a fabric simulation motion model according to the coordinate position obtained by the coordinate calculating module.
6. The system for building a fabric simulation model according to claim 5, wherein the coordinate calculation module specifically comprises:
the initial definition unit is used for acquiring the instantaneous vector of the gravity center of any quadrilateral mesh patch at the current moment;
the tensor calculation unit is used for calculating the Riemannian degree tensor of any quadrilateral mesh patch at the current moment and the Riemannian curvature tensor at the current moment according to the instantaneous vector defined by the primary definition unit;
a parameter calculation unit which calculates internal deformation energy of each quadrilateral mesh patch at the current time according to the initial Riemannian degree tensor, the initial Riemannian curvature tensor, the Riemannian degree tensor at the current time and the Riemannian curvature tensor obtained by the tensor calculation unit;
and the coordinate calculation unit is used for calculating the coordinate position of the next moment by combining the stress analysis of the quadrilateral mesh surface patch and the internal deformation energy obtained by the parameter calculation unit.
7. The system for building a fabric simulation model according to claim 6, wherein the parameter calculation unit specifically includes:
defining a sub-unit, defining the region of any quadrilateral mesh patch as omega, omega { (l)1,l2)∈R2Define the instantaneous vector of the center of gravity as
Figure FDA0002231512630000041
A parameter calculating subunit that calculates the stretch-shear energy at the current time from the initial riemannian tensor obtained by the tensor calculating unit, the riemannian tensor at the current time, and the parameters defined by the defining subunit,
Figure FDA0002231512630000042
wherein the content of the first and second substances,
Figure FDA0002231512630000043
for the stretch-shear energy at the present moment, GtAs the Riemann scale number of the current time,G0Is the initial number of the Riemann scale,
Figure FDA0002231512630000044
Ewpand EwtRespectively the Young modulus of the fabric in the warp and weft directions, H is the rigidity modulus,
Figure FDA0002231512630000045
is a matrix wsA Hilbert-Schmitt norm of the weighted measure tensor variabilities;
the parameter calculating subunit calculates the bending energy at the current time according to the initial Riemann curvature tensor obtained by the tensor calculating unit, the Riemann curvature tensor at the current time and the parameters defined by the defining subunit,
Figure FDA0002231512630000051
wherein the content of the first and second substances,
Figure FDA0002231512630000052
is the bending energy at the present moment, RtIs the curvature tensor of the current time, R0Is the initial riemann curvature tensor,Rfwpand RfwtBending stiffness of the fabric in the warp and weft directions respectively,
Figure FDA0002231512630000054
is a matrix wbA Hilbert-Schmidt norm of the weighted curvature tensor variations;
the parameter calculation subunit calculates the internal deformation energy from the tensile shear energy and the bending energy
Figure FDA0002231512630000055
8. The system for building a fabric simulation model according to claim 7, wherein the coordinate calculation unit specifically includes:
a stress calculating subunit for calculating the elastic internal force of the quadrilateral mesh patch according to the internal deformation energy
Figure FDA0002231512630000056
The stress calculation subunit calculates the inertial internal force
Figure FDA0002231512630000057
Viscous damping forceWeight of fabric
Figure FDA0002231512630000059
And air resistance
Figure FDA00022315126300000510
Wherein the content of the first and second substances,
Figure FDA00022315126300000511
p is the density of the fabric,
Figure FDA00022315126300000512
mu is a damping coefficient;
a stress analysis subunit for obtaining a dynamic equation reflecting the motion of the fabric according to the first dynamic principle and the force calculated by the stress calculation subunit,
Figure FDA00022315126300000513
and the coordinate calculation subunit calculates the coordinate position of the next moment according to the dynamic equation obtained by the stress analysis subunit.
9. A storage medium having a computer program stored thereon, characterized in that: the computer program, when executed by a processor, implements the method of any of claims 1 to 4.
10. An electronic device comprising a memory and a processor, the memory having stored thereon a computer program that runs on the processor, characterized in that: the processor, when executing the computer program, implements the method of any of claims 1 to 4.
CN201910969172.2A 2019-10-12 2019-10-12 Method and system for establishing fabric simulation model, storage medium and electronic equipment Active CN110717297B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910969172.2A CN110717297B (en) 2019-10-12 2019-10-12 Method and system for establishing fabric simulation model, storage medium and electronic equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910969172.2A CN110717297B (en) 2019-10-12 2019-10-12 Method and system for establishing fabric simulation model, storage medium and electronic equipment

Publications (2)

Publication Number Publication Date
CN110717297A true CN110717297A (en) 2020-01-21
CN110717297B CN110717297B (en) 2022-07-12

Family

ID=69211479

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910969172.2A Active CN110717297B (en) 2019-10-12 2019-10-12 Method and system for establishing fabric simulation model, storage medium and electronic equipment

Country Status (1)

Country Link
CN (1) CN110717297B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113435098A (en) * 2021-06-30 2021-09-24 西南交通大学 Method for accurately simulating appearance of deformed fabric thin-layer soft substance

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102262781A (en) * 2011-05-11 2011-11-30 浙江工业大学 Method for registration of ink-jet printing texture image based on unit decomposition optical flow field
CN106530338A (en) * 2016-10-31 2017-03-22 武汉纺织大学 Method and system for matching MR image feature points before and after nonlinear deformation of biological tissue
CN108647411A (en) * 2018-04-25 2018-10-12 桂林电子科技大学 A kind of truss stress model modeling method expressed based on the spaces Grassmann and recurrence surface
CN109002571A (en) * 2018-05-08 2018-12-14 杭州电子科技大学 Cloth dynamic emulation method based on equal geometry mass-spring modeling
WO2019037657A1 (en) * 2017-08-23 2019-02-28 湖南碳康生物科技有限公司 Carbon fiber composite artificial bone and preparation method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102262781A (en) * 2011-05-11 2011-11-30 浙江工业大学 Method for registration of ink-jet printing texture image based on unit decomposition optical flow field
CN106530338A (en) * 2016-10-31 2017-03-22 武汉纺织大学 Method and system for matching MR image feature points before and after nonlinear deformation of biological tissue
WO2019037657A1 (en) * 2017-08-23 2019-02-28 湖南碳康生物科技有限公司 Carbon fiber composite artificial bone and preparation method
CN108647411A (en) * 2018-04-25 2018-10-12 桂林电子科技大学 A kind of truss stress model modeling method expressed based on the spaces Grassmann and recurrence surface
CN109002571A (en) * 2018-05-08 2018-12-14 杭州电子科技大学 Cloth dynamic emulation method based on equal geometry mass-spring modeling

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
LI,SS 等: "Construction and initial analysis of five Fosmid libraries of mitochondrial genomes of cotton (Gossypium)", 《CHINESE SCIENCE BULLETIN》 *
吴明珠 等: "基于零曲率的轻柔织物建模与动态仿真研究", 《丝绸》 *
朱俊平等: "基于谱多流形学习的聚类方法研究", 《数学的实践与认识》 *
陈瑛 等: "基于黎曼零曲率的织物静态造型研究", 《计算机仿真》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113435098A (en) * 2021-06-30 2021-09-24 西南交通大学 Method for accurately simulating appearance of deformed fabric thin-layer soft substance
CN113435098B (en) * 2021-06-30 2022-11-15 西南交通大学 Method for accurately simulating appearance of deformed fabric thin-layer soft substance

Also Published As

Publication number Publication date
CN110717297B (en) 2022-07-12

Similar Documents

Publication Publication Date Title
Müller et al. Position based dynamics
CN101496028A (en) Method of simulating deformable object using geometrically motivated model
Duan et al. Volume preserved mass–spring model with novel constraints for soft tissue deformation
Li et al. Space-time editing of elastic motion through material optimization and reduction
US8666712B2 (en) Physical simulations on a graphics processor
US7791606B2 (en) Goal-directed cloth simulation
CN109002630B (en) Rapid simulation method for super-elastic material
US8392154B2 (en) System and method for real-time cloth simulation
US20230061175A1 (en) Real-Time Simulation of Elastic Body
CN104360896A (en) Parallel fluid simulation acceleration method based on GPU (Graphics Processing Unit) cluster
JP4936522B2 (en) Image processing method and image processing apparatus
CN112862942B (en) Physical special effect simulation method, device, electronic equipment and storage medium
CN111488670B (en) Nonlinear mass point spring soft tissue deformation simulation method
CN110717297B (en) Method and system for establishing fabric simulation model, storage medium and electronic equipment
CN109002571B (en) Cloth dynamic simulation method based on isogeometric spring particle model
Bender et al. Efficient Cloth Simulation Using an Adaptive Finite Element Method.
Cetinaslan ESPEFs: Exponential spring potential energy functions for simulating deformable objects
CN114912303A (en) Three-dimensional fluid display method and device, electronic equipment and storage medium
CN109118561B (en) Hierarchical dynamic simulation method based on position
Choi et al. Image morphing using mass-spring system
Cetinaslan et al. Energy Embedded Gauss-Seidel Iteration for Soft Body Simulations
Gu et al. Constraint solving order in position based dynamics
Mandal et al. Physics-based Mesh Deformation with Haptic Feedback and Material Anisotropy
Song et al. Cutting and fracturing models without remeshing
Choi et al. A fast overlap control in image morphing using mass-spring system

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