Battery electrochemical-mechanical multi-physical field coupling method and device

Abstract

Embodiments of the present application provide a battery electrochemical-mechanical multi-physical field coupling method and device. The method comprises: constructing a gradient weighted smooth field based on the battery data; constructing an equation weak form of an electrochemical field and an equation weak form of a mechanical field based on the gradient weighted smooth field; and constructing a coupling relationship between the electrochemical field and the mechanical field based on the equation weak form of the electrochemical field and the equation weak form of the mechanical field. In this way, the numerical defects of low calculation precision of the finite element algorithm can be solved, and higher calculation precision can be obtained by using low-order tetrahedral elements, thereby greatly reducing the calculation time while ensuring the calculation precision.

Description

TECHNICAL FIELD

[0001] Embodiments of the present application relate to the field of data processing, and in particular to a battery electrochemical-mechanical multi-physical field coupling method and device. BACKGROUND

[0002] Computational accuracy and efficiency are important challenges faced by the electrochemical-mechanical multi-physical field coupling simulation analysis of nonlinear transient lithium batteries.

[0003] Currently, the nonlinear transient electrochemical-mechanical multi-field coupling simulation analysis for the field of lithium batteries usually adopts the traditional finite element method, and the coupling method usually adopts strong coupling based on the same unit stiffness matrix to assemble the electrochemical field and mechanics, or separately calculates the electrochemical field and the mechanical field, and then transmits the coupling variables in the form of external force for coupling. This coupling method not only suffers from the problems of low computational accuracy and computational efficiency inherent in the finite element method, but also suffers from the problems of nonlinear non-convergence and excessive iteration caused by strong coupling or weak coupling. SUMMARY

[0004] According to embodiments of the present application, a lithium battery electrochemical-mechanical multi-field coupling scheme based on a gradient weighting algorithm is provided, which can solve the numerical defects of low computational accuracy of the finite element algorithm, and can obtain higher computational accuracy by using low-order tetrahedral elements, thereby greatly reducing the computational time while ensuring computational accuracy; and can be applied to the electrochemical-mechanical coupling simulation scene in the electrode particle-cell multi-scale scene.

[0005] In a first aspect of the present application, a battery electrochemical-mechanical multi-physical field coupling method is provided. The method comprises:

[0006] obtaining battery data to be processed;

[0007] based on the battery data, constructing a gradient weighted smoothing domain;

[0008] based on the gradient weighted smoothing domain, constructing an equation weak form of an electrochemical field and an equation weak form of a mechanical field, respectively;

[0009] based on the equation weak form of the electrochemical field and the equation weak form of the mechanical field, constructing a coupling relationship between the electrochemical field and the mechanical field.

[0010] Further, the constructing a gradient weighted smoothing domain based on the battery data comprises:

[0011] based on the battery data, discretizing the battery solution domain into background grid elements of linear triangles and / or linear tetrahedrons;

[0012] Set the triangle and / or tetrahedron element in the background mesh as the primary element; set the adjacent element sharing the common edge or face with the primary element as the secondary element, and construct the gradient weighted smoothing domain.

[0013] Further, the constructing the weak form of the electrochemical field equation based on the gradient weighted smoothing domain comprises:

[0014] ;

[0015] ;

[0016] ;

[0017] ;

[0018] wherein,

[0019] N i represents the i-node shape function of the gradient weighted domain;

[0020] N j represents the j-node shape function of the gradient weighted domain;

[0021] c s and c l respectively represent the concentration field variable of the electrode region and the electrolyte region;

[0022] D s and D e respectively represent the concentration diffusion coefficient of the electrode region and the electrolyte region;

[0023] k eff represents the ionic conductivity of the electrolyte region;

[0024] and respectively represent the electrode potential and the electrolyte potential;

[0025] represents the gradient weighted smoothing integral domain;

[0026] represents the electrode region conductivity;

[0027] represents the material parameter at the electrode;

[0028] represents the external force coupling term of the electric field;

[0029] external force coupling term representing the concentration field;

[0030] a composite function partial derivative matrix representing the gradient weighted domain of node i;

[0031] a composite function partial derivative matrix representing the gradient weighted domain of node j;

[0032] a concentration field variable of the electrolyte region at the previous time step;

[0033] a material parameter at the electrolyte.

[0034] Further, the composite function partial derivative matrix can be calculated by:

[0035] ;

[0036] where w gr a weighting factor of the composite gradient weighted smoothing domain;

[0037] m represents the number of elements in the smoothing domain;

[0038] a gradient operator symbol;

[0039] N j a j-node shape function of the gradient weighted domain.

[0040] Further, the weighting factor of the composite gradient weighted smoothing domain can be calculated by:

[0041] ;

[0042] where, a main element volume;

[0043] an auxiliary element volume;

[0044] an auxiliary element number.

[0045] Further, the equation weak form of the mechanical field based on the gradient weighted smoothing domain comprises:

[0046] ;

[0047] where u s a discrete element displacement vector;

[0048] to be a mechanical material constitutive parameter.

[0049] Further, the constructing the coupling relationship between the electrochemical field and the mechanical field based on the weak form of the equation of the electrochemical field and the weak form of the equation of the mechanical field comprises:

[0050] constructing a single-field stiffness matrix of the electrochemical field and the mechanical field in the composite gradient smooth domain based on the weak form of the equation of the electrochemical field and the weak form of the equation of the mechanical field;

[0051] constructing the coupling relationship between the electrochemical field and the mechanical field based on the single-field stiffness matrix.

[0052] In a second aspect of the present application, a battery electrochemical-mechanical multi-physical field coupling device is provided. The device comprises:

[0053] an acquisition module configured to acquire battery data to be processed;

[0054] a first construction module configured to construct a gradient weighted smooth domain based on the battery data;

[0055] a second construction module configured to construct a weak form of an equation of an electrochemical field and a weak form of an equation of a mechanical field based on the gradient weighted smooth domain, respectively;

[0056] a coupling module configured to construct a coupling relationship between the electrochemical field and the mechanical field based on the weak form of the equation of the electrochemical field and the weak form of the equation of the mechanical field.

[0057] In a third aspect of the present application, an electronic device is provided. The electronic device comprises a memory and a processor, wherein the memory stores a computer program, and the processor implements the method as described above when executing the program.

[0058] In a fourth aspect of the present application, a computer readable storage medium is provided, which stores a computer program, and the program is executed by a processor to implement the method according to the first aspect of the present application.

[0059] The battery electrochemical-mechanical multi-physical field coupling method provided by the embodiments of the present application can solve the numerical defect of low calculation precision of the finite element algorithm by acquiring battery data to be processed, constructing a gradient weighted smooth domain based on the battery data, constructing a weak form of an equation of an electrochemical field and a weak form of an equation of a mechanical field based on the gradient weighted smooth domain, and constructing a coupling relationship between the electrochemical field and the mechanical field based on the weak form of the equation of the electrochemical field and the weak form of the equation of the mechanical field. Moreover, a higher calculation precision can be obtained by using a low-order tetrahedral element, which greatly reduces the calculation time while ensuring the calculation precision.

[0060] It is to be understood that the description in the summary is not intended to identify key or essential features of embodiments of the application, nor is it intended to limit the scope of the application. Other aspects of the application will be evident from the following description. BRIEF DESCRIPTION OF DRAWINGS

[0061] The above and other features, aspects, and advantages of embodiments of the present application will become more apparent from the following description in conjunction with the accompanying drawings. In the drawings:

[0062] Figure 1 a flow chart of a battery electro-chemo-mechanical multi-physics coupling method according to an embodiment of the present application;

[0063] Figure 2 a schematic diagram of a composite gradient-weighted smooth domain according to an embodiment of the present application; wherein, Figure 2 (a) is a schematic diagram of a background mesh; Figure 2 (b) is a schematic diagram of a gradient smooth domain;

[0064] Figure 3 a block diagram of a battery electro-chemo-mechanical multi-physics coupling device according to an embodiment of the present application;

[0065] Figure 4 a schematic diagram of a structure of a terminal device or a server suitable for implementing embodiments of the present application. DETAILED DESCRIPTION

[0066] In order to make the objects, technical solutions and advantages of the embodiments of the present disclosure clearer, the technical solutions in the embodiments of the present disclosure will be described clearly and completely below with reference to the drawings in the embodiments of the present disclosure. Obviously, the described embodiments are only some, but not all of the embodiments of the present disclosure. Based on the embodiments in the present disclosure, all other embodiments obtained by a person of ordinary skill in the art without creative work fall within the protection scope of the present disclosure.

[0067] In addition, the term "and / or" in this document is only a description of the association relationship of the associated objects, which means that there can be three relationships, for example, A and / or B can represent the three cases of A alone, A and B together, and B alone. In addition, the character " / " in this document generally represents an "or" relationship between the front and rear associated objects.

[0068] Figure 1 A flow chart of a battery electro-chemo-mechanical multi-physics coupling method according to an embodiment of the present disclosure is shown. The method comprises:

[0069] S110, obtaining battery data to be processed.

[0070] In some embodiments, the battery data to be processed can be acquired in a limited or wireless manner; the battery data includes battery solution domain and other related data.

[0071] In S120, a gradient-weighted smooth domain is constructed based on the battery data.

[0072] In some embodiments, as Figure 2 (a) shown, the battery solution domain is discretized into background mesh units in the shape of linear triangles and / or linear tetrahedrons, as Figure 2 (b) shown, the triangle or tetrahedron units in the background mesh are taken as primary units, and the adjacent units sharing edges or faces with the primary units are taken as auxiliary units, so as to construct a gradient-smooth domain with the primary units as the main part. The nodes in the smooth domain are formed by the nodes of the primary units and the auxiliary units, and thus the gradient-smooth domain node discrete information different from the background mesh can be formed.

[0073] In S130, the weak form of the electrochemical field equation and the weak form of the mechanical field equation are respectively constructed based on the gradient-weighted smooth domain.

[0074] In view of the defects of the traditional electrochemical field and mechanical field, i.e. the defects of constructing unit information and unit stiffness matrix by using the finite element method (the overall stiffness matrix of the solution domain is "too stiff", leading to low calculation precision of the linear tetrahedron unit), the present disclosure proposes an improved method for constructing the electrochemical field and the mechanical field based on the gradient-weighted unit algorithm.

[0075] Specifically, the weak form of the electrochemical field equation and the mechanical field equation is constructed based on the discrete information of the gradient-smooth domain; and the unit stiffness matrix of the electrochemical field and the mechanical field is constructed based on the shared gradient-smooth domain node shape function:

[0076] The weak form of the electrochemical field equation includes:

[0077] ;

[0078] ;

[0079] ;

[0080] ;

[0081] wherein,

[0082] N i represents the i node shape function of the gradient-weighted domain;

[0083] N j represents the j node shape function of the gradient-weighted domain;

[0084] c s and c l denote the concentration field variable in the electrode region and electrolyte region, respectively;

[0085] D s and D e denote the concentration diffusion coefficient in the electrode region and electrolyte region, respectively;

[0086] k eff denotes the ionic conductivity in the electrolyte region;

[0087] and denote the electrode potential and electrolyte potential, respectively;

[0088] denotes the gradient-weighted smooth integral domain;

[0089] denotes the electrode region conductivity;

[0090] denotes the material parameter at the electrode;

[0091] denotes the external force coupling term of the electric field;

[0092] denotes the external force coupling term of the concentration field;

[0093] denotes the composite function partial derivative matrix of the gradient-weighted domain of node i;

[0094] denotes the composite function partial derivative matrix of the gradient-weighted domain of node j;

[0095] denotes the concentration field variable in the electrolyte region at the previous time step;

[0096] denotes the material parameter at the electrolyte.

[0097] Further, the composite function partial derivative matrix in the above gradient smooth domain can be composed of the main element in the composite domain and the auxiliary element adjacent to the main element by surface or edge, that is:

[0098] ;

[0099] wherein m denotes the number of elements in the smooth domain;

[0100] denotes the gradient operator symbol;

[0101] N j denotes the j-node shape function of the gradient weighted domain;

[0102] w gr denotes the weighted factor of the composite gradient weighted smoothing domain, which can be expressed as:

[0103] ;

[0104] wherein, is the main unit volume;

[0105] is the auxiliary unit volume;

[0106] is the auxiliary unit number;

[0107] The weak form of the equation of the mechanical field includes:

[0108] ;

[0109] wherein, u s denotes the discrete unit displacement vector;

[0110] is the mechanical material constitutive parameter.

[0111] S140, based on the weak form of the equation of the electrochemical field and the weak form of the equation of the mechanical field, the coupling relationship between the electrochemical field and the mechanical field is constructed.

[0112] In some embodiments, after the electrochemical-mechanical correction coupling model is constructed by step S130 and the staggered iteration stabilization processing is constructed, the single field stiffness matrix of the electrochemical field and the mechanical field in the composite gradient smoothing domain is obtained; the mutual coupling relationship between the electrochemical field and the mechanical field can be further constructed.

[0113] Specifically, in the application scenario of actual electrochemical-mechanical coupling simulation analysis, the structural mechanical response speed is much faster than the electrochemical reaction. Therefore, it is necessary to first construct the key items of the coupling and influence of the mechanical field on the electrochemical field in the composite gradient smoothing domain. It mainly includes the influence and correction of the mechanical response on the concentration mass transfer and the influence and correction of the mechanical response on the electrochemical reaction:

[0114] ;

[0115] ;

[0116] wherein, and denote the electrode potential and the electrolyte potential, respectively;

[0117] η gr represents the modified potential in the gradient-weighted domain;

[0118] R , T and F are the molar gas constant, temperature and Faraday constant, respectively;

[0119] E represents the unit equilibrium potential;

[0120] α a , α c are the anode and cathode transfer numbers, respectively;

[0121] represents the exchange current density;

[0122] N represents the unit shape function;

[0123] D s represents the electrode area concentration diffusion coefficient;

[0124] n represents the number of electrons;

[0125] molar volume;

[0126] and represent the boundary domain and the gradient-smoothed domain, respectively;

[0127] represents the stress in the gradient-weighted domain;

[0128] Further, the coupling influence matrix of electrochemical field and mechanical structure deformation is constructed. Considering the volume change of the battery as a whole caused by lithium ion intercalation and deintercalation during the charging and discharging process of the lithium battery, the coupling relationship between the electrochemical field and the mechanical structure deformation can be constructed by the following way:

[0129] 1) Obtain the SOC state of the electrode area in the composite gradient-smoothed domain;

[0130] 2) Perform regional search to obtain the average particle size of the electrode particles or the average particle size of the electrode particles obtained based on the actual processing design process; slice the electrode area of the battery along the vertical direction of the thickness, height and width of the electrode area;

[0131] 3) Based on the slice region in step 2), all the cells in the region are cycled to reconstruct the gradient-weighted sub-domain in the slice region, and the SOC of the composite gradient-weighted domain is calculated through the composite gradient-weighted domain in the slice region to calculate the SOC of the composite gradient-weighted domain;

[0132] 4) Based on the SOC, the average strain along the vertical direction of the thickness, height and width of the cell in the slice region can be calculated by the following formula:

[0133] ;

[0134] wherein, ε ave represents the average strain along the vertical direction of the thickness, height and width of the cell in the slice region;

[0135] f i The test parameters can be modified or interpolated based on the test parameters;

[0136] 5) The average strain is applied to the mechanical field for coupling calculation.

[0137] In summary, the coupling relationship between the electrochemical field and the mechanical structure deformation can be obtained.

[0138] According to the embodiments of the present disclosure, the following technical effects are achieved:

[0139] The low-order linear unit is used to efficiently and accurately solve the electrochemical-mechanical multi-field coupling problem in lithium battery simulation analysis.

[0140] The mutual coupling and interaction relationship between the electrochemistry and the mechanics is constructed in the composite gradient-weighted smoothing domain, which has higher coupling accuracy and convergence than the traditional finite element method based on the construction of the unit.

[0141] By constructing the coupling between the stress and the concentration field mass transfer and electrode reaction kinetics (overpotential) based on the gradient-weighted smoothing domain, the calculation accuracy of the lithium battery electrochemical-mechanical simulation analysis is greatly improved.

[0142] By reconstructing the gradient-weighted sub-domain in the slice region through the cycle of all the cells in the region, and calculating the SOC of the composite gradient-weighted domain through the composite gradient-weighted domain in the slice region, the coupling effect of the electrochemical field on the mechanical field is obtained, and the electrochemical-mechanical coupling accuracy and the high-precision mapping from the electrode particle scale to the cell scale are improved.

[0143] It should be noted that, for the foregoing method embodiments, for the sake of simple description, they are all described as a combination of a series of actions, but those skilled in the art should know that the present application is not limited by the action sequence described, because according to the present application, certain steps can be performed in other sequences or simultaneously. Secondly, those skilled in the art should know that the embodiments described in the specification are all optional embodiments, and the actions and modules involved are not necessarily required by the present application.

[0144] The above is the introduction of the method embodiment, and the scheme described in the present application will be further described through the device embodiment.

[0145] Figure 3 A block diagram 300 of a battery electrochemical-mechanical multi-physical field coupling device according to an embodiment of the present application is shown, as shown in Figure 3 includes:

[0146] The acquisition module 310 is configured to acquire battery data to be processed.

[0147] The first construction module 320 is configured to construct a gradient weighted smooth field based on the battery data.

[0148] The second construction module 330 is configured to construct an equation weak form of an electrochemical field and an equation weak form of a mechanical field, respectively, based on the gradient weighted smooth field.

[0149] The coupling module 340 is configured to construct a coupling relationship between the electrochemical field and the mechanical field based on the equation weak form of the electrochemical field and the equation weak form of the mechanical field.

[0150] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working process of the described modules can refer to the corresponding process in the foregoing method embodiments, which will not be described here.

[0151] Figure 4 A structural schematic diagram of a terminal device or a server suitable for implementing an embodiment of the present application is shown.

[0152] As Figure 4 shown, the terminal device or the server includes a central processing unit (CPU) 401, which can perform various appropriate actions and processes according to programs stored in a read-only memory (ROM) 402 or programs loaded from a storage portion 408 into a random access memory (RAM) 403. In the RAM 403, various programs and data required for the operation of the terminal device or the server are also stored. The CPU 401, the ROM 402, and the RAM 403 are connected to each other through a bus 404. An input / output (I / O) interface 405 is also connected to the bus 404.

[0153] The following components are connected to the I / O interface 405: an input part 406 including a keyboard, a mouse, etc.; an output part 407 including a display such as a cathode ray tube (CRT), a liquid crystal display (LCD), etc., and a speaker, etc.; a storage part 408 including a hard disk, etc.; and a communication part 409 including a network interface card such as a LAN card, a modem, etc. The communication part 409 performs communication processing via a network such as the Internet. A drive 410 is also connected to the I / O interface 405 as necessary. A removable media 411 such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, etc. is attached to the drive 410 as necessary, so that a computer program read out therefrom is installed in the storage part 408 as necessary.

[0154] In particular, the above method flow steps can be implemented as a computer software program in accordance with embodiments of the present application. For example, embodiments of the present application include a computer program product comprising a computer program carried on a machine-readable medium, the computer program containing program code for executing the methods illustrated in the flowcharts. In such embodiments, the computer program can be downloaded and installed from a network via the communication part 409, and / or installed from the removable media 411. When the computer program is executed by the central processing unit (CPU) 401, the above-described functions defined in the system of the present application are executed.

[0155] It should be noted that the computer-readable medium shown in the application can be a computer-readable signal medium or a computer-readable storage medium or any combination of the two. The computer-readable storage medium may, for example, but is not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device or component, or any combination of the above. More specific examples of computer-readable storage media can include, but are not limited to, an electrical connection having one or more wires, a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above. In this application, the computer-readable storage medium can be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, device or component. In this application, the computer-readable signal medium can include a data signal carried in a baseband or as a carrier wave in a carrier wave, which carries computer-readable program code. Such a propagated data signal can take many forms, including but not limited to an electromagnetic signal, an optical signal, or any suitable combination of the above. The computer-readable signal medium can also be any computer-readable medium other than the computer-readable storage medium, which can send, propagate or transmit a program for use by or in conjunction with an instruction execution system, device or component. The program code contained on the computer-readable medium can be transmitted using any suitable medium, including but not limited to wireless, wire, optical cable, RF, etc., or any suitable combination of the above.

[0156] The flowcharts and block diagrams in the drawings illustrate the possible implementation architectures, functions and operations of the systems, methods and computer program products according to various embodiments of the application. In this regard, each block in the flowcharts or block diagrams can represent a module, a program segment, or a portion of code, which contains one or more executable instructions for implementing the specified logical functions. It should also be noted that in some alternative implementations, the functions noted in the blocks can occur in different orders than that shown in the drawings. For example, two blocks that are shown in succession can actually be executed substantially concurrently, or they can sometimes be executed in reverse order, depending on the involved functions. It should also be noted that each block in the block diagrams and / or flowcharts, and combinations of blocks in the block diagrams and / or flowcharts, can be implemented by a dedicated hardware-based system that performs the specified functions or operations, or can be implemented by a combination of special-purpose hardware and computer instructions.

[0157] The units or modules described in the embodiments of the present application can be implemented in the form of software or in the form of hardware. The units or modules described can also be arranged in a processor. In some cases, the names of the units or modules do not constitute a limitation on the units or modules themselves.

[0158] As another aspect, the present application also provides a computer readable storage medium, which can be included in the electronic device described in the above embodiments, or can exist separately without being assembled into the electronic device. The computer readable storage medium stores one or more programs, and the programs are used by one or more processors to execute the methods described in the present application.

[0159] The above description is merely preferred embodiments of the present application and a description of the principles of the technology used. Those skilled in the art should understand that the scope of the application described in the present application is not limited to the technical solutions formed by the specific combinations of the above technical features, and should also cover other technical solutions formed by the combinations of the above technical features or their equivalent features without departing from the above application concept. For example, the technical solutions formed by the mutual replacement of the above features and the technical features applied in the present application (but not limited to) having similar functions.

Claims

1. A battery electrochemical-mechanical multi-physical field coupling method, characterized in that, The method comprises: acquiring battery data to be processed; constructing a gradient-weighted smooth domain based on the battery data; constructing a weak form of an equation of an electrochemical field and a weak form of an equation of a mechanical field based on the gradient-weighted smooth domain respectively; wherein the constructing the weak form of the equation of the electrochemical field based on the gradient-weighted smooth domain comprises: ； ； ； ； wherein N i i-nodal shape function representing the gradient-weighted domain; N j jth nodal shape function representing the gradient-weighted domain; c s and denote the concentration field variable of the electrode region and the electrolyte region, respectively; D s and D e D and D represent the concentration diffusion coefficients of the electrode region and the electrolyte region, respectively; k eff represents the ionic conductivity of the electrolyte region; and are denoted as electrode potential and electrolyte potential, respectively; denotes the gradient-weighted smoothing integral domain; represents the electrode area conductivity; representing material parameters at the electrodes; external force coupling term representing an electric field; external force coupling term representing the concentration field; a matrix of partial derivatives of the composite function representing the gradient weighted domain of node i; a matrix of partial derivatives of the composite function representing the gradient weighted domain of node j; denotes the concentration field variable of the electrolyte region at the previous time step; denotes the material parameter at the electrolyte; the constructing the weak form of the equation of the mechanical field based on the gradient-weighted smooth domain comprises: ； wherein u s represents a discrete unit displacement vector; Mechanical material constitutive parameters; constructing a coupling relationship between the electrochemical field and the mechanical field based on the weak form of the equation of the electrochemical field and the weak form of the equation of the mechanical field; the constructing the coupling relationship between the electrochemical field and the mechanical field based on the weak form of the equation of the electrochemical field and the weak form of the equation of the mechanical field comprises: constructing a single-field stiffness matrix of the electrochemical field and the mechanical field in the composite gradient smooth domain based on the weak form of the equation of the electrochemical field and the weak form of the equation of the mechanical field; constructing the coupling relationship between the electrochemical field and the mechanical field based on the single-field stiffness matrix.

2. The method of claim 1, wherein, The constructing the gradient-weighted smooth domain based on the battery data comprises: discretizing a battery solving domain into background grid units in linear triangular and / or linear tetrahedral shapes based on the battery data; setting the triangular and / or tetrahedral units in the background grid as main units; setting adjacent units sharing edges or faces with the main units as auxiliary units to construct the gradient-weighted smooth domain.

3. The method of claim 2, wherein, The composite gradient-weighted smooth domain is calculated by the following method: ； where w gr denotes the weighting factor of the composite gradient weighted smoothing domain; m represents the number of units in the smooth domain; denotes the gradient operator symbol; N i a gradient-weighted domain is represented by i a nodal shape function.

4. The method of claim 3, wherein, The weighting factor of the composite gradient-weighted smooth domain is calculated by the following method: ； wherein, Vp is the volume of the primary unit; To assist unit volume; is the number of auxiliary units.

5. A battery electrochemical-mechanical multi-physical field coupling device, characterized in that, The method comprises: an acquisition module configured to acquire battery data to be processed; a first construction module configured to construct a gradient-weighted smooth domain based on the battery data; a second construction module configured to construct a weak form of an equation of an electrochemical field and a weak form of an equation of a mechanical field based on the gradient-weighted smooth domain respectively; wherein the constructing the weak form of the equation of the electrochemical field based on the gradient-weighted smooth domain comprises: ； ； ； ； wherein N i an i-node shaped function representing a gradient weighted domain; N j j nodal shape functions representing the gradient weighted domain; c s and denote the concentration field variable in the electrode region and in the electrolyte region, respectively; D s and D e Dc and De represent the concentration diffusion coefficients of the electrode region and the electrolyte region, respectively; k eff represents the ionic conductivity of the electrolyte region; and are denoted as electrode potential and electrolyte potential, respectively; denotes the gradient-weighted smoothing integral domain; represents the electrode area conductivity; representing material parameters at the electrodes; external force coupling term representing an electric field; external force coupling term representing the concentration field; a matrix of partial derivatives of the composite function representing the gradient weighted domain of node i; a matrix of partial derivatives of the composite function representing the gradient weighted domain of node j; concentration field variable representing the electrolyte region at the previous time step; represents the material parameter at the electrolyte; the constructing the weak form of the equation of the mechanical field based on the gradient-weighted smooth domain comprises: ； wherein u s represents a discrete unit displacement vector; Mechanical material constitutive parameters; a coupling module configured to construct a coupling relationship between the electrochemical field and the mechanical field based on the weak form of the equation of the electrochemical field and the weak form of the equation of the mechanical field; the constructing the coupling relationship between the electrochemical field and the mechanical field based on the weak form of the equation of the electrochemical field and the weak form of the equation of the mechanical field comprises: constructing a single-field stiffness matrix of the electrochemical field and the mechanical field in the composite gradient smooth domain based on the weak form of the equation of the electrochemical field and the weak form of the equation of the mechanical field; constructing the coupling relationship between the electrochemical field and the mechanical field based on the single-field stiffness matrix. 6.An electronic device comprising a memory and a processor, the memory having stored thereon a computer program, characterized in that, The processor executes the computer program to implement the method in any one of claims 1-4.

7. A computer-readable storage medium having stored thereon a computer program, characterized in that, The computer program is executed by the processor to implement the method in any one of claims 1-4.

Images