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 PDFInfo
- 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
Links
- 239000004744 fabric Substances 0.000 title claims abstract description 113
- 238000004088 simulation Methods 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000003860 storage Methods 0.000 title claims abstract description 14
- 230000005484 gravity Effects 0.000 claims abstract description 53
- 238000004458 analytical method Methods 0.000 claims abstract description 29
- 238000004364 calculation method Methods 0.000 claims description 36
- 238000005452 bending Methods 0.000 claims description 23
- 238000004590 computer program Methods 0.000 claims description 15
- 238000013016 damping Methods 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 10
- 239000006185 dispersion Substances 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 6
- 238000010008 shearing Methods 0.000 claims description 5
- 108010001267 Protein Subunits Proteins 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 9
- 238000011160 research Methods 0.000 abstract description 8
- 239000000463 material Substances 0.000 abstract description 5
- 230000007547 defect Effects 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 230000005489 elastic deformation Effects 0.000 description 3
- 238000005034 decoration Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000001413 cellular effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000000053 physical method Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling 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
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
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,wherein,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,Ewpand EwtRespectively the Young modulus of the fabric in the warp and weft directions, H is the rigidity modulus,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,as the bending energy at the present time is,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,is a matrix wbA Hilbert-Schmidt norm of the weighted curvature tensor variations;
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
Calculating inertial internal forceViscous damping forceWeight of fabricAnd air resistanceWherein,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
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,wherein,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,Ewpand EwtRespectively the Young modulus of the fabric in the warp and weft directions, H is the rigidity modulus,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,wherein,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,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
The stress calculation subunit calculates the inertial internal forceViscous damping forceWeight of fabricAnd air resistanceWherein,p is the density of the fabric,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,
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 asS432 calculates the stretch-shear energy at the current time from the initial riemann gauge number and the riemann gauge number at the current time,wherein,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,Ewpand EwtRespectively the Young modulus of the fabric in the warp and weft directions, H is the rigidity modulus,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,wherein,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,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
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 asDifferential 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: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 calculatedS442 calculating the inertial internal forceViscous damping forceWeight of fabricAnd air resistanceWherein,p is the density of the fabric,mu is a damping coefficient; s443, according to the first dynamic principle, obtaining a dynamic equation reflecting the movement of the fabric,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 energyCalculating inertial internal forceViscous damping forceWeight of fabricAnd air resistanceThe 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
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,wherein,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,Ewpand EwtRespectively the Young modulus of the fabric in the warp and weft directions, H is the rigidity modulus,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,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,is a matrix wbWeighted curvature tensor variationalThe hilbert-schmitt norm of (d);
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
Calculating inertial internal forceViscous damping forceWeight of fabricAnd air resistanceWherein,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.
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
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,wherein,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,Ewpand EwtRespectively the Young modulus of the fabric in the warp and weft directions, H is the rigidity modulus,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,wherein,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,is a matrix wbA Hilbert-Schmidt norm of the weighted curvature tensor variations;
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
The stress calculation subunit calculates the inertial internal forceViscous damping forceWeight of fabricAnd air resistanceWherein,p is the density of the fabric,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,
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.
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)
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)
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 |
-
2019
- 2019-10-12 CN CN201910969172.2A patent/CN110717297B/en active Active
Patent Citations (5)
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)
Title |
---|
LI,SS 等: "Construction and initial analysis of five Fosmid libraries of mitochondrial genomes of cotton (Gossypium)", 《CHINESE SCIENCE BULLETIN》 * |
吴明珠 等: "基于零曲率的轻柔织物建模与动态仿真研究", 《丝绸》 * |
朱俊平等: "基于谱多流形学习的聚类方法研究", 《数学的实践与认识》 * |
陈瑛 等: "基于黎曼零曲率的织物静态造型研究", 《计算机仿真》 * |
Cited By (2)
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 |
---|---|---|
Li et al. | Space-time editing of elastic motion through material optimization and reduction | |
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 | |
US7791606B2 (en) | Goal-directed cloth simulation | |
EP2069986A2 (en) | Physical simulations on a graphics processor | |
US20190108300A1 (en) | Methods for realistic and efficient simulation of moving objects | |
CN111488670B (en) | Nonlinear mass point spring soft tissue deformation simulation method | |
CN109002630B (en) | Rapid simulation method for super-elastic material | |
US8392154B2 (en) | System and method for real-time cloth simulation | |
CN104360896A (en) | Parallel fluid simulation acceleration method based on GPU (Graphics Processing Unit) cluster | |
US20230061175A1 (en) | Real-Time Simulation of Elastic Body | |
JP4936522B2 (en) | Image processing method and image processing apparatus | |
CN112862942B (en) | Physical special effect simulation method, device, electronic equipment and storage medium | |
CN109002571B (en) | Cloth dynamic simulation method based on isogeometric spring particle model | |
CN105677950A (en) | Fabric deformation simulation method | |
Bender et al. | Adaptive cloth simulation using corotational finite elements | |
CN110717297B (en) | Method and system for establishing fabric simulation model, storage medium and electronic equipment | |
Bender et al. | Efficient Cloth Simulation Using an Adaptive Finite Element Method. | |
CN115455753A (en) | Collision simulation method and device for soft tissue and rigid ground | |
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 | |
Gu et al. | Constraint solving order in position based dynamics | |
Yao et al. | Realistic and stable animation of cloth |
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 |