Multi-sound field synergistic copper alloy calendering material modulation control method and system

Abstract

The invention provides a copper alloy calendering material modulation control method and system based on multiple sound field collaboration, and relates to the field, and the method comprises the steps of obtaining tissue state information, dividing sound field action partitions, calculating phase regulation parameters based on propagation time difference, establishing a reverse mapping relation between tissue density and sound field energy, and collecting dislocation line orientation distribution. And calculating an included angle concentration degree index, and optimizing the phase difference and energy ratio between the excitation sources of the adjacent sound fields according to the included angle concentration degree index to realize accurate regulation and control of dislocation orientation of the calendered material, so that the calendering performance of the copper alloy is improved.

Description

TECHNICAL FIELD

[0001] The present application relates to metal calendering technology, and in particular to a copper alloy calendering material modulation control method and system with multiple sound fields. BACKGROUND

[0002] With the development of high-end equipment manufacturing and electronic information technology, higher requirements are put forward for the performance of copper alloy calendering materials. As a non-contact energy transmission method, sound field has shown unique advantages in the field of material processing, providing a new technical approach for the modulation of copper alloy calendering materials.

[0003] The current sound field modulation technology for copper alloy calendering materials has obvious deficiencies. The existing sound field action mode mainly uses single sound source excitation, which cannot implement differential modulation according to the organizational gradient difference in the thickness direction of the calendering material, resulting in uneven distribution of sound field energy and unsatisfactory modulation effect. The existing technology fails to establish an accurate mapping relationship between dislocation orientation distribution and sound field parameters, making it difficult to achieve directional modulation of dislocation orientation, limiting the performance improvement of calendering materials, and failing to meet the precise control requirements of high-end applications for the anisotropic performance of copper alloy calendering materials. SUMMARY

[0004] The embodiments of the present application provide a copper alloy calendering material modulation control method and system with multiple sound fields, which can solve the problems in the prior art.

[0005] In a first aspect, the embodiments of the present application provide a copper alloy calendering material modulation control method with multiple sound fields, comprising:

[0006] Obtaining the organizational state information of the copper alloy material to be calendered and the target calendering performance index;

[0007] Dividing the calendering material along the thickness direction into multiple sound field action zones according to the organizational state information, and arranging sound field excitation sources in each zone;

[0008] Detecting the propagation time of each sound field excitation source at different depth layers, calculating the propagation time difference, and converting the propagation time difference into a phase control parameter;

[0009] Extracting the organizational density distribution curve of the calendering material, establishing an inverse mapping relationship between the organizational density distribution curve and the output energy of the sound field excitation source, and calculating the energy distribution coefficient of each depth layer;

[0010] Driving each sound field excitation source to apply a cooperative sound field action to the calendering material according to the phase control parameter and the energy distribution coefficient, collecting the dislocation line orientation distribution of the calendering material under the cooperative sound field action, counting the angle distribution between the dislocation line and the calendering direction, and calculating the angle concentration index;

[0011] Determine the target strengthening coefficient of dislocation orientation based on the target calendering performance index, calculate the sound wave superposition gain according to the angle concentration index, compare the sound wave superposition gain with the target strengthening coefficient, optimize the phase difference and energy ratio between adjacent sound field excitation sources based on the comparison result, and realize the coordinated regulation of dislocation orientation of the calendering material.

[0012] According to the organization state information, the calendering material is divided into multiple sound field action partitions along the thickness direction, and the sound field excitation sources are arranged in each partition, which includes:

[0013] According to the organization state information, the organization parameters of the calendering material at different thickness positions are extracted, the organization difference value between adjacent thickness positions is calculated, the organization mutation position is determined along the thickness direction of the calendering material according to the organization difference value, and the calendering material is divided into multiple sound field action partitions with the organization mutation position as the boundary, and the organization characteristic value of each sound field action partition is calculated;

[0014] Based on the organization characteristic value and the size of the sound field action partition, the number of sound field excitation sources required in each sound field action partition is calculated, the sound field excitation sources are subjected to sound field superposition test, the action radius of the sound field excitation sources is obtained, and the arrangement interval of the sound field excitation sources is determined according to the action radius and the number of sound field excitation sources;

[0015] According to the arrangement interval and the size of the sound field action partition, the arrangement coordinates of the sound field excitation sources and the corresponding sound field intensity distribution are calculated, the arrangement coordinates satisfying the preset sound field coverage threshold are selected, and the sound field excitation sources are arranged at the selected arrangement coordinates.

[0016] Detect the propagation time of each sound field excitation source at different depth layers, calculate the propagation time difference, and convert the propagation time difference into phase regulation parameters, which includes:

[0017] Collect the sound wave propagation signals of each sound field excitation source at different depth layers, and extract the sound wave frequency drift and sound wave amplitude attenuation rate through time-frequency domain joint analysis of the sound wave propagation signals;

[0018] Invert the dislocation density distribution of each depth layer through the sound wave frequency drift, invert the grain boundary density distribution of each depth layer through the sound wave amplitude attenuation rate, calculate the acoustic impedance correction value of each depth layer based on the dislocation density distribution and the grain boundary density distribution, and correct the propagation time of each depth layer according to the acoustic impedance correction value to obtain the corrected propagation time;

[0019] Calculate the difference between the corrected propagation times to obtain the propagation time difference, perform phase conversion on the propagation time difference and the working frequency of each sound field excitation source to obtain the excitation phase difference value, and construct a space-time distribution map of the excitation phase difference value changing with the depth layer position;

[0020] extract a spatial gradient vector field of the excitation phase difference values from the space-time distribution map, identify a depth layer position with a maximum gradient amplitude in the spatial gradient vector field, and take an excitation phase difference value of a sound field excitation source corresponding to the depth layer position with the maximum gradient amplitude as a reference phase difference value;

[0021] calculate a phase deviation between an excitation phase difference value of another depth layer sound field excitation source and the reference phase difference value as a phase control parameter.

[0022] extract a tissue density distribution curve of the calendered material, establish an inverse mapping relationship between the tissue density distribution curve and an output energy of the sound field excitation source, and calculate an energy distribution coefficient of each depth layer, including:

[0023] obtain tissue density values of different depth layers of the calendered material along a thickness direction, arrange the tissue density values according to depth layer positions, and construct a tissue density distribution curve;

[0024] calculate a change rate of the tissue density values between adjacent depth layers, identify a depth layer with a change rate deviating from a mean value and mark it as a tissue density abnormal area;

[0025] read a tissue density value of the tissue density abnormal area, establish a sound wave propagation impedance distribution, deduce a sound wave propagation path according to the sound wave propagation impedance distribution, and calculate a sound wave energy cumulative attenuation amount;

[0026] take the sound wave energy cumulative attenuation amount as an output energy increment requirement of the sound field excitation source of each depth layer, and establish an inverse mapping relationship between the tissue density distribution curve and the output energy of the sound field excitation source;

[0027] read the output energy increment requirement of the sound field excitation source of each depth layer, and construct a depth transmission link of the output energy increment requirement of the sound field excitation source;

[0028] identify an anchor depth layer with a mutation in the depth transmission link, perform bidirectional diffusion correction on the output energy increment requirement of the sound field excitation source of adjacent depth layers from the anchor depth layer as a starting point, normalize the corrected output energy increment requirement of the sound field excitation source, and obtain an energy distribution coefficient of each depth layer.

[0029] drive each sound field excitation source to exert a cooperative sound field action on the calendered material according to the phase control parameter and the energy distribution coefficient, collect a dislocation line orientation distribution of the calendered material under the cooperative sound field action, count an angle distribution between the dislocation line and the calendering direction, and calculate an angle concentration index, including:

[0030] convert the phase control parameter into a frequency control amount of the sound field excitation source, convert the energy distribution coefficient into a power control amount of the sound field excitation source, and generate a sound wave according to the frequency control amount and the power control amount;

[0031] The propagation time and pressure amplitude of the collected sound wave in the calendered material are calculated to obtain the propagation time difference and pressure amplitude difference of adjacent sound field excitation sources, and the propagation characteristics of the sound wave in the depth direction are obtained.

[0032] The beam pointing angle of the sound field excitation source is adjusted according to the propagation characteristics, the beam pointing angle and the pressure amplitude difference are combined to determine the superposition position of the sound wave, a cooperative sound field is formed, the stress state of the calendered material in the action area of the cooperative sound field is collected, the driving force of the dislocation motion is calculated according to the stress state, and the starting point of the dislocation slip is determined.

[0033] The position coordinates of the dislocation line are recorded from the starting point, the spatial trajectory of the dislocation line is constructed according to the position coordinates, the node information of the dislocation line is extracted, the local direction of the dislocation line is calculated according to the node information, the local direction of the dislocation line is projected to the calendering direction, and the dislocation orientation distribution is obtained.

[0034] The included angle values in the dislocation orientation distribution are counted, the number distribution of the included angle values is calculated, and the included angle concentration index is obtained by normalizing the number distribution.

[0035] The target enhancement coefficient of the dislocation orientation is determined based on the target calendering performance index, the sound wave superposition gain is calculated according to the included angle concentration index, the sound wave superposition gain is compared with the target enhancement coefficient, and the phase difference and energy ratio between adjacent sound field excitation sources are optimized based on the comparison result to realize the cooperative regulation of the dislocation orientation of the calendered material.

[0036] The stress requirement and texture requirement in the target calendering performance index are obtained, the crystal orientation distribution and dislocation slip distribution of the calendered material are calculated, the critical shear force of the dislocation motion is determined according to the crystal orientation distribution and dislocation slip distribution, and the critical shear force is converted into the target enhancement coefficient of the dislocation orientation.

[0037] The sound wave propagation signals of adjacent sound field excitation sources are collected, the sound wave propagation time difference and pressure amplitude difference are calculated, the included angle value distribution in the included angle concentration index is extracted, the dynamic response function of the included angle value distribution and the sound wave propagation signal is established, the output power of the adjacent sound field excitation source is adjusted according to the dynamic response function, and the adjusted sound wave propagation signal is converted into the sound wave superposition gain.

[0038] The difference between the sound wave superposition gain and the target enhancement coefficient is taken as the regulation deviation of the dislocation orientation, the phase compensation value and the power distribution ratio of the adjacent sound field excitation source are calculated according to the regulation deviation, and the output parameters of the sound field excitation source are adjusted based on the phase compensation value and the power distribution ratio to realize the cooperative regulation of the dislocation orientation of the calendered material.

[0039] The dynamic response function of the included angle value distribution and the sound wave propagation signal is established, and the output power of the adjacent sound field excitation source is adjusted according to the dynamic response function.

[0040] Statistical analysis is performed on the included angle value distribution data to obtain the peak position and distribution width of the included angle value, the peak position is taken as the dominant direction parameter of dislocation orientation, and the distribution width is taken as the dispersion degree parameter of dislocation orientation;

[0041] The sound wave propagation signals inside the calendered material of adjacent sound field excitation sources are collected, the sound wave arrival time and the sound wave peak amplitude are extracted by time domain analysis on the sound wave propagation signals, and the sound wave propagation time difference and the sound wave peak amplitude difference between adjacent sound field excitation sources are calculated.

[0042] The dominant direction parameter and the dispersion degree parameter are set as output variables of the dynamic response function, the sound wave propagation time difference and the sound wave peak amplitude difference are set as input variables of the dynamic response function, and the dynamic response function of the included angle value distribution and the sound wave propagation signal is constructed by fitting the functional relationship between the input variables and the output variables.

[0043] The partial derivatives of the dynamic response function with respect to the sound wave propagation time difference and the sound wave peak amplitude difference are calculated, the sound wave propagation time difference target value and the sound wave peak amplitude difference target value are converted into the output power adjustment amount of the adjacent sound field excitation source according to the partial derivatives, and the output power adjustment amount is applied to the corresponding sound field excitation source to complete the adjustment of the output power of the adjacent sound field excitation source.

[0044] In a second aspect of the embodiment of the present application, a multi-sound field cooperative copper alloy calendering material modulation control system is provided, comprising:

[0045] A first unit is configured to obtain the organizational state information of the copper alloy material to be calendered and the target calendering performance index.

[0046] A second unit is configured to divide the calendering material into multiple sound field action sub-zones along the thickness direction according to the organizational state information, and arrange sound field excitation sources in each sub-zone.

[0047] A third unit is configured to detect the propagation time of each sound field excitation source at different depth layers, calculate the propagation time difference, and convert the propagation time difference into a phase control parameter.

[0048] A fourth unit is configured to extract the organizational density distribution curve of the calendering material, establish the inverse mapping relationship between the organizational density distribution curve and the output energy of the sound field excitation source, and calculate the energy distribution coefficient of each depth layer.

[0049] A fifth unit is configured to drive each sound field excitation source to apply a cooperative sound field action to the calendering material according to the phase control parameter and the energy distribution coefficient, collect the dislocation line orientation distribution of the calendering material under the cooperative sound field action, statistically analyze the included angle distribution between the dislocation line and the calendering direction, and calculate an included angle concentration index.

[0050] The sixth unit is configured to determine a target enhancement coefficient of dislocation orientation based on a target calendering performance index, calculate a sound wave superposition gain according to the included angle concentration index, compare the sound wave superposition gain with the target enhancement coefficient, and optimize the phase difference and energy ratio between adjacent sound field excitation sources based on the comparison result to realize the coordinated regulation of the dislocation orientation of the calendering material.

[0051] In a third aspect, the present application provides an electronic device, comprising:

[0052] a processor;

[0053] a memory for storing processor-executable instructions;

[0054] The processor is configured to invoke the instructions stored in the memory to execute the method described above.

[0055] In a fourth aspect, the present application provides a computer-readable storage medium having computer program instructions stored thereon, wherein the computer program instructions are executed by a processor to implement the method described above.

[0056] In the embodiment, the sound field is precisely controlled at different depth layers by partitioning the calendering material, and the utilization efficiency of the sound field energy is improved. The sound field energy of each depth layer can be accurately configured according to the reverse mapping relationship established based on the tissue density distribution, and the uniformity of the sound field effect is ensured. The dislocation orientation of the calendering material can be coordinated and regulated by using the concentration index of the included angle between the dislocation orientation distribution and the calendering direction, and dynamically optimizing the sound wave superposition gain, so as to improve the uniformity of the tissue and the performance stability of the calendering material, and finally obtain a copper alloy calendering product with excellent performance. BRIEF DESCRIPTION OF DRAWINGS

[0057] Figure 1 A flowchart of the copper alloy calendering material modulation control method of the present application is shown in the figure.

[0058] Figure 2 A flowchart of the energy distribution correction algorithm of the present application is shown in the figure. DETAILED DESCRIPTION

[0059] To make the objectives, technical solutions, and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be described below in conjunction with the accompanying drawings of the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative work fall within the scope of protection of the present application.

[0060] The technical solutions of the present application are described in detail below with specific examples. The following specific examples can be combined with each other, and some examples may not be described in detail for the same or similar concepts or processes.

[0061] Figure 1 The flowchart of the copper alloy calendering material modulation control method for the multi-sound field coordination of the embodiments of the present application is shown in Figure 1 The method comprises:

[0062] Obtaining the organizational state information of the copper alloy material to be calendered and the target calendering performance index;

[0063] According to the organizational state information, the calendering material is divided into multiple sound field action partitions along the thickness direction, and sound field excitation sources are arranged in each partition;

[0064] Detecting the propagation time of each sound field excitation source at different depth layers, calculating the propagation time difference, and converting the propagation time difference into a phase control parameter;

[0065] Extracting the organizational density distribution curve of the calendering material, establishing the inverse mapping relationship between the organizational density distribution curve and the output energy of the sound field excitation source, and calculating the energy distribution coefficient of each depth layer;

[0066] According to the phase control parameter and the energy distribution coefficient, each sound field excitation source is driven to exert a coordinated sound field action on the calendering material, the dislocation line orientation distribution of the calendering material under the coordinated sound field action is collected, the angle distribution between the dislocation line and the calendering direction is counted, and the angle concentration index is calculated;

[0067] Based on the target calendering performance index, a target enhancement coefficient of dislocation orientation is determined, the sound wave superposition gain is calculated according to the angle concentration index, the sound wave superposition gain is compared with the target enhancement coefficient, the phase difference and the energy ratio between adjacent sound field excitation sources are optimized based on the comparison result, and the coordinated regulation of the dislocation orientation of the calendering material is realized.

[0068] According to the organizational state information, the calendering material is divided into multiple sound field action partitions along the thickness direction, and sound field excitation sources are arranged in each partition;

[0069] According to the organizational state information, the organizational parameters of the calendering material at different thickness positions are extracted, the organizational difference value between adjacent thickness positions is calculated, the organizational mutation position is determined along the thickness direction of the calendering material according to the organizational difference value, the calendering material is divided into multiple sound field action partitions with the organizational mutation position as the boundary, and the organizational characteristic value of each sound field action partition is calculated;

[0070] The number of sound field excitation sources required in each sound field action subzone is calculated based on the size of the sound field action subzone and the tissue characteristic value, a sound field superposition test is performed on the sound field excitation sources, the action radius of the sound field excitation sources is obtained, and the arrangement spacing of the sound field excitation sources is determined according to the action radius and the number of sound field excitation sources;

[0071] The arrangement coordinates of the sound field excitation sources and the corresponding sound field intensity distribution are calculated according to the arrangement spacing and the size of the sound field action subzone, the arrangement coordinates that meet the preset sound field coverage threshold are selected according to the sound field intensity distribution, and the sound field excitation sources are arranged at the selected arrangement coordinates.

[0072] The organization parameters of the calendered material at different thickness positions are extracted according to the organization state information, the organization difference value between adjacent thickness positions is calculated, and then the organization mutation position is determined and the sound field action subzone is divided. The grain size, grain boundary distribution, second phase distribution and other organization parameters of each thickness position are obtained. For a copper alloy calendered material, a measurement point can be taken every 0.5 mm in the thickness direction, and the average grain size at this position is recorded. The difference between the grain sizes of the adjacent two measurement points is the organization difference value, and when the organization difference value exceeds the preset threshold, it is determined as the organization mutation position. Taking a brass calendered material with a thickness of 10 mm as an example, 20 measurement points are set in the thickness direction, and it is found that the organization difference values at positions of 2.5 mm and 7.5 mm are 5 microns and 6 microns respectively, both of which exceed the preset threshold of 4 microns, so the two positions are determined as the organization mutation positions, and the calendered material is divided into three sound field action subzones: a surface subzone of 0-2.5 mm, a middle subzone of 2.5-7.5 mm and a bottom subzone of 7.5-10 mm. The organization characteristic value of each subzone is calculated, including the weighted average value of the parameters such as the average grain size, the standard deviation of the grain size, and the grain boundary area ratio in the subzone. The calculation results show that the organization characteristic values of the three subzones are 0.68, 0.82 and 0.75 respectively.

[0073] The number of sound field excitation sources needed in each sound field action partition is calculated based on the tissue characteristic value and the size of the sound field action partition. The number of sound field excitation sources is related to the tissue characteristic value and the size of the partition. The larger the tissue characteristic value, the more heterogeneous the tissue, the more sound field excitation sources needed. The larger the size of the partition, the more sound field excitation sources needed. Sound field superposition test is performed on the sound field excitation sources. Multiple sound field excitation sources are placed on the calendering material surface, the distance between them is adjusted, the sound field intensity distribution is measured, and the effective range of the sound field is determined. The action radius of the sound field excitation source is defined as the distance corresponding to the sound field intensity decaying to 50% of the initial value. For the above brass calendering material, the action radius of the ultrasonic transducer used is determined to be 25 mm through the sound field superposition test. The layout spacing of the sound field excitation source is determined according to the action radius and the number of sound field excitation sources needed. The layout spacing should be less than twice the action radius to ensure the continuity of the sound field coverage. For the three partitions of the surface layer, the middle layer and the bottom layer, the number of sound field excitation sources needed is calculated to be 4, 6 and 5 respectively according to the respective tissue characteristic value and the size of the partition. Considering that the action radius is 25 mm, the layout spacing is determined to be 40 mm, ensuring that the sound field of adjacent excitation sources has enough overlapping area.

[0074] The layout coordinates of the sound field excitation source and the corresponding sound field intensity distribution are calculated according to the layout spacing and the size of the sound field action partition. The layout coordinates of the sound field excitation source are calculated by the partition size and the layout spacing. Taking a rectangular calendering material as an example, the length is 500 mm, the width is 300 mm, 4 sound field excitation sources are laid out in the surface layer partition, and the layout coordinates are (100, 75), (100, 225), (400, 75) and (400, 225) respectively, with the coordinate unit being millimeter and the coordinate origin being located at one corner of the calendering material. The sound field intensity distribution at each layout position is calculated, and the sound field intensity decreases with the increase of the distance from the excitation source. The sound field intensity distribution can be calculated by the sound field propagation model, considering the parameters such as the acoustic impedance and absorption coefficient of the material. The layout coordinates that meet the preset sound field coverage threshold are selected. The preset sound field coverage threshold is 30% of the initial sound field intensity, that is, the sound field intensity at any position in the partition is required to be not less than 30% of the sound field intensity at the excitation source.

[0075] The calculation can be made that the 4 layout coordinates of the surface layer partition can make 95% of the area meet the coverage threshold requirement, the 6 layout coordinates of the middle partition can make 96% of the area meet the coverage threshold requirement, and the 5 layout coordinates of the bottom layer partition can make 94% of the area meet the coverage threshold requirement. After the layout coordinates are determined, the sound field excitation sources are laid out at the selected positions. Different frequencies and powers of ultrasonic transducers can be selected as sound field excitation sources for different partitions. The transducer with a frequency of 20 kHz and a power of 500 W is selected for the surface layer partition, the transducer with a frequency of 25 kHz and a power of 600 W is selected for the middle partition, and the transducer with a frequency of 30 kHz and a power of 700 W is selected for the bottom layer partition.

[0076] The present application can divide the calendering material into multiple sound field action partitions according to the differences in the organizational state of the calendering material in the thickness direction, and arrange appropriate sound field excitation sources according to the organizational characteristics of each partition to realize accurate regulation of the sound field. Compared with the traditional single sound field action mode, the non-uniformity of the internal organization of the calendering material is fully considered, and through the synergistic action of multiple sound fields, the utilization efficiency of sound field energy is significantly improved, and the modulation effect of the sound field on the material organization is enhanced.

[0077] In an optional embodiment, the propagation time of each sound field excitation source at different depth layers is detected, the propagation time difference is calculated, and the propagation time difference is converted into a phase control parameter, which includes:

[0078] The sound wave propagation signals of each sound field excitation source at different depth layers are collected, and the sound wave frequency drift and the sound wave amplitude attenuation rate are extracted by joint time-frequency domain analysis of the sound wave propagation signals;

[0079] The dislocation density distribution of each depth layer is inverted by the sound wave frequency drift, the grain boundary density distribution of each depth layer is inverted by the sound wave amplitude attenuation rate, the acoustic impedance correction value of each depth layer is calculated based on the dislocation density distribution and the grain boundary density distribution, and the propagation time of each depth layer is corrected according to the acoustic impedance correction value to obtain the corrected propagation time;

[0080] The difference between the corrected propagation times is calculated to obtain the propagation time difference, the propagation time difference and the working frequency of each sound field excitation source are phase converted to obtain the excitation phase difference value, and a space-time distribution map of the excitation phase difference value changing with the depth layer position is constructed;

[0081] The spatial gradient vector field of the excitation phase difference value is extracted from the space-time distribution map, the depth layer position with the maximum gradient amplitude in the spatial gradient vector field is identified, and the excitation phase difference value of the sound field excitation source corresponding to the depth layer position with the maximum gradient amplitude is taken as the reference phase difference value;

[0082] The phase deviation between the excitation phase difference values of the sound field excitation sources at other depth layers and the reference phase difference value is calculated as the phase control parameter.

[0083] In the embodiment, the sound wave propagation signals of each sound field excitation source at different depth layers are collected, and the sound wave frequency drift and the sound wave amplitude attenuation rate are extracted by time-frequency domain joint analysis. During the collection process, multiple ultrasonic transducers are arranged on the surface of the copper alloy calendered material as sound field excitation sources, and sound wave receiving probes are arranged at different depth positions inside the calendered material. Taking a brass calendered material with a thickness of 10 mm as an example, a receiving point is set every 1 mm in the depth direction, and a total of 10 depth layers are set. The pulse-echo method is used to measure the sound wave propagation signals at each depth layer, and the ultrasonic flaw detector is used to collect the ultrasonic echo signals, with a sampling frequency of 100 MHz and a sampling time of 100 microseconds. After the collected sound wave propagation signals are preprocessed, the time-frequency domain joint analysis is performed on the signals by short-time Fourier transform, and the sound wave frequency drift and the sound wave amplitude attenuation rate are extracted. In the time-frequency spectrum analysis process, the Hanning window is selected as the window function, the window length is set to 512 points, and the overlap rate is 50%. By analyzing the change of the main frequency peak value in the time-frequency spectrum, the sound wave frequency drift at each depth layer is extracted. For the sound wave with an initial frequency of 20 kHz, the frequency drift measured at a depth of 5 mm is -320 Hz. By calculating the natural logarithm of the ratio of the amplitude of the echo signal at each depth layer to the amplitude of the initial transmitted signal, the sound wave amplitude attenuation rate is obtained. Similarly, the amplitude attenuation rate measured at a depth of 5 mm is 0.42 dB / mm.

[0084] The dislocation density distribution at each depth layer is inverted by the sound wave frequency drift, the grain boundary density distribution at each depth layer is inverted by the sound wave amplitude attenuation rate, the sound impedance correction value at each depth layer is calculated based on the dislocation density distribution and the grain boundary density distribution, and the corrected propagation time is obtained by correcting the propagation time at each depth layer according to the sound impedance correction value. The frequency drift is affected by the dislocations in the microstructure of the material, and the frequency drift is related to the dislocation density. According to the physical principle of sound wave propagation in metal materials, the frequency drift is approximately proportional to the dislocation density, and the frequency drift can be converted into dislocation density by a calibration coefficient. For the measured brass calendered material, the dislocation density corresponding to the frequency drift of -320 Hz is about 8×10 10 / cm 2 . The sound wave amplitude attenuation is mainly affected by the grain boundary scattering, and the amplitude attenuation rate is closely related to the grain boundary density. By establishing the corresponding relationship between the amplitude attenuation rate and the grain boundary density, the grain boundary density distribution at each depth layer can be inverted. For the measured amplitude attenuation rate of 0.42 dB / mm, the corresponding grain boundary density is about 4×10 4 / mm 2The dislocation density and the grain boundary density both affect the acoustic impedance of the material. By comprehensively considering the two parameters, the acoustic impedance correction value can be calculated. When calculating the acoustic impedance correction value, the dislocation density and the grain boundary density are respectively given a weight of 0.6 and 0.4. The acoustic impedance correction value at a depth of 5 mm is calculated to be 0.87. The corrected propagation time is obtained by multiplying the original propagation time by the acoustic impedance correction value. For a depth layer with an original propagation time of 6.5 microseconds, the corrected propagation time is 5.66 microseconds.

[0085] The difference between the corrected propagation times is calculated to obtain the propagation time difference. The propagation time difference is phase-converted with the operating frequency of each sound field excitation source to obtain the excitation phase difference value. A space-time distribution map of the excitation phase difference value changing with the depth layer position is constructed. The difference between the corrected propagation times of adjacent depth layers is the propagation time difference. For two depth layers with depths of 4 mm and 5 mm, the corrected propagation times are 4.53 microseconds and 5.66 microseconds, respectively, and the propagation time difference is 1.13 microseconds. The propagation time difference can be converted into a phase difference by the operating frequency. The calculation formula is: the phase difference is equal to the propagation time difference multiplied by the operating frequency and then multiplied by 360 degrees. For a sound field excitation source with an operating frequency of 20 kHz, a propagation time difference of 1.13 microseconds is converted into a phase difference of 8.14 degrees. The data of all depth layers of the calendered material are processed to obtain the excitation phase difference values of different depth layers, and a space-time distribution map of the excitation phase difference value changing with the depth is constructed. The space-time distribution map is represented in the form of a heat map, with the horizontal coordinate representing the depth position and the vertical coordinate representing the time.

[0086] The spatial gradient vector field of the excitation phase difference value is extracted from the space-time distribution map, and the depth layer position with the maximum gradient amplitude in the spatial gradient vector field is identified. The excitation phase difference value of the sound field excitation source corresponding to the depth layer position with the maximum gradient amplitude is taken as the reference phase difference value. The spatial gradient vector field represents the rate of change of the excitation phase difference value in space, which can be obtained by calculating the difference between the phase difference values of adjacent depth layers. For the space-time distribution map of the brass calendered material, the calculated spatial gradient vector field shows that the gradient amplitude is maximum at a depth of 7 mm, reaching 4.3 degrees / mm, and the phase difference value at this position is 12.5 degrees, which is taken as the reference phase difference value.

[0087] The phase deviation between the excitation phase difference value of the other depth layer sound field excitation source and the reference phase difference value is calculated as a phase control parameter. The difference between the phase difference value and the reference phase difference value of each depth layer is the phase deviation, which is the phase control parameter of the sound field excitation source. For a depth of 2 mm, the phase difference value is 3.6 degrees, and the phase deviation is -8.9 degrees; for a depth of 5 mm, the phase difference value is 8.14 degrees, and the phase deviation is -4.36 degrees. These phase control parameters are used to adjust the initial phase of each sound field excitation source to enable the sound field to form enhanced interference at the target depth layer, thereby achieving precise modulation of the specific depth layer. In specific implementation, the phase of the driving signal of each sound field excitation source is controlled by a digital signal processor to achieve precise phase control.

[0088] In the embodiment, the synergistic effect of multiple sound fields in the copper alloy rolling material can be effectively improved. By accurately detecting the propagation characteristics of sound waves in different depth layers and combining the propagation time correction with the material microstructure parameters, the sound field energy is accurately focused at the target depth layer. This multiple sound field synergistic method based on phase control breaks through the limitation of traditional sound field control methods in the depth direction, and can realize directional energy transfer to the specific depth layer according to the distribution characteristics of the internal structure of the rolling material, thereby accurately controlling the material structure evolution process at the microscale.

[0089] As shown in Figure 2 , the energy distribution correction algorithm flow of the embodiment is shown.

[0090] In an optional implementation, the organization density distribution curve of the rolling material is extracted, an inverse mapping relationship between the organization density distribution curve and the output energy of the sound field excitation source is established, and the energy distribution coefficient of each depth layer is calculated, including:

[0091] The organization density values of the rolling material at different depth layers along the thickness direction are obtained, the organization density values are arranged according to the depth layer positions, and an organization density distribution curve is constructed;

[0092] The change rate of the organization density values between adjacent depth layers is calculated, the depth layers with a change rate deviating from the mean value are identified and marked as organization density abnormal regions;

[0093] The organization density values of the organization density abnormal regions are read, a sound wave propagation impedance distribution is established, a sound wave propagation path is deduced according to the sound wave propagation impedance distribution, and a sound wave energy cumulative attenuation amount is calculated;

[0094] The sound wave energy cumulative attenuation amount is taken as the output energy increment requirement of the sound field excitation source of each depth layer, and an inverse mapping relationship between the organization density distribution curve and the output energy of the sound field excitation source is established;

[0095] reading the sound field excitation source output energy increment requirement of each depth layer, constructing a depth transfer link of the sound field excitation source output energy increment requirement;

[0096] identifying the anchor point depth layer where the mutation occurs in the depth transfer link, performing bidirectional diffusion correction on the sound field excitation source output energy increment requirement of the adjacent depth layer from the anchor point depth layer as the starting point, and normalizing the corrected sound field excitation source output energy increment requirement to obtain the energy distribution coefficient of each depth layer.

[0097] In the embodiment, the organization density values of the calendered material at different depth layers along the thickness direction are obtained, and the organization density values are arranged according to the depth layer positions to construct an organization density distribution curve. The organization density represents the distribution of discontinuities such as pores and defects in the material, and is an important indicator for evaluating the material organization quality. The organization density values are obtained by combining metallographic microscope observation and image analysis. For example, metallographic samples are taken from different depths of a 12 mm thick H62 brass calendered material, and after grinding, polishing and etching, an organization photo is taken under a metallographic microscope. The photo resolution is 2048x1536 pixels, and the magnification is 200 times. The photo is binarized by using an image analysis software, and the ratio of the black area (pores, defects, etc.) to the entire field area is calculated. The organization density is obtained by subtracting the ratio from 1. The calendered material is sampled every 1 mm from the surface to the bottom to obtain the organization density values of 12 depth layers. The organization density of the surface layer (0-1 mm) is 0.985, the organization density of the middle layer (5-6 mm) is 0.962, and the organization density of the bottom layer (11-12 mm) is 0.978. These data points are arranged according to the depth, and a continuous organization density distribution curve is generated by cubic spline interpolation. The horizontal coordinate of the curve is the depth, and the vertical coordinate is the organization density.

[0098] The change rate of the organization density values between adjacent depth layers is calculated, and the depth layer whose change rate deviates from the mean value is identified and marked as an organization density abnormal area. The change rate of the organization density between adjacent depth layers is calculated by dividing the difference between the organization densities of adjacent two points by the depth interval. For the H62 brass calendered material, 11 change rate values are calculated, and the average value of these change rates is 0.002 / mm. The change rate deviating from the mean value by more than 2 times the standard deviation is set as an abnormal value, and the standard deviation is calculated as 0.003 / mm. Therefore, the depth layer whose absolute value of the change rate is greater than 0.008 / mm is marked as an organization density abnormal area. Through analysis, the change rate at the depth of 3-4 mm is -0.009 / mm, and the change rate at the depth of 8-9 mm is 0.011 / mm. These two areas are marked as organization density abnormal areas.

[0099] The tissue density value of the abnormal region of the tissue density of the read tissue is read, the sound wave propagation impedance distribution is established, the sound wave propagation path is deduced according to the sound wave propagation impedance distribution, and the sound wave energy cumulative attenuation amount is calculated. The sound wave propagation impedance is directly related to the tissue density of the material, and the two are approximately linearly related. The sound wave propagation impedance calculation formula is: the sound wave propagation impedance is equal to the density of the material multiplied by the sound wave propagation speed in the material. For H62 brass, the basic sound wave propagation impedance is about 30×10 6 kg / m 2 ·s, and on this basis, the change is corrected according to the change of the tissue density. The sound wave propagation impedance increases by about 0.5×10 6 kg / m 2 ·s for each 0.01 decrease in tissue density. According to this relationship, the sound wave propagation impedance at a depth of 3-4 mm is calculated to be 30.9×10 6 kg / m 2 ·s, and the sound wave propagation impedance at a depth of 8-9 mm is calculated to be 29.7×10 6 kg / m 2 ·s. During the propagation of the sound wave, reflection and transmission occur due to the change in the sound impedance, resulting in energy attenuation. The sound wave energy attenuation amount is proportional to the square of the difference in the sound impedance. Through the sound propagation model calculation, the sound wave energy cumulative attenuation amount from the surface to a depth of 3-4 mm is about 2.1 dB, and the sound wave energy cumulative attenuation amount from the surface to a depth of 8-9 mm is about 4.3 dB.

[0100] The sound wave energy cumulative attenuation amount is taken as the energy increment requirement of the sound field excitation source output at each depth layer, and the reverse mapping relationship between the tissue density distribution curve and the sound field excitation source output energy is established. The energy increment requirement of the sound field excitation source output is proportional to the sound wave energy cumulative attenuation amount. In order to offset the attenuation of the sound wave during the propagation, the output energy of the sound field excitation source needs to be adjusted accordingly. The reverse mapping relationship is established by a fitting function, and for the H62 brass calendering material, an exponential function is used: the energy increment requirement of the sound field excitation source output is equal to the basic energy multiplied by the exponential function of the sound wave energy cumulative attenuation amount. The basic energy is set to 100 W, and the exponential coefficient is 0.15. According to this relationship, the energy increment requirement of the sound field excitation source output at a depth of 3-4 mm is calculated to be 36.7 W, and the energy increment requirement of the sound field excitation source output at a depth of 8-9 mm is calculated to be 90.2 W.

[0101] The energy increment requirement of the sound field excitation source output of each depth layer is read, and a depth transfer link of the energy increment requirement of the sound field excitation source output is constructed. The depth transfer link is a data structure representing the transfer relationship between energy requirements at different depth layers. The construction method is to arrange the energy increment requirements of each depth layer in depth order to form a one-dimensional array, and then calculate the transfer coefficient through the difference between adjacent elements. For the 12 depth layers of the H62 brass calendering material, the constructed depth transfer link contains 12 nodes, and each node records the energy increment requirement of the depth layer and the transfer coefficient with the adjacent depth layer.

[0102] An anchor point depth layer where a mutation occurs in the depth transfer link is identified, and the energy increment requirement of the sound field excitation source output of the adjacent depth layer is bidirectionally diffused and corrected from the anchor point depth layer as the starting point. The energy distribution coefficient of each depth layer is obtained by normalizing the corrected energy increment requirement of the sound field excitation source output. In the depth transfer link, the node where the transfer coefficient between adjacent nodes changes by more than 50% is identified as a mutation point, and the corresponding depth layer is the anchor point depth layer. For the H62 brass calendering material, two anchor point depth layers of 3-4 mm and 8-9 mm are identified. The energy increment requirement of the adjacent depth layer is bidirectionally diffused and corrected from the anchor point depth layer as the starting point. The diffusion correction adopts an exponential decay model, and the correction amplitude decreases with increasing distance. After correction, the energy increment requirement of the 12 depth layers is normalized to make the sum equal to 1, and the energy distribution coefficient of each depth layer is obtained. After normalization, the energy distribution coefficient of the 3-4 mm depth layer is 0.12, the energy distribution coefficient of the 8-9 mm depth layer is 0.19, and the energy distribution coefficients of the remaining depth layers are between 0.05 and 0.10.

[0103] In this embodiment, the copper alloy calendering material modulation control process under the synergistic action of multiple sound fields can be effectively guided. The non-uniformity of the internal organization of the calendering material is fully considered, and precise energy compensation is performed for the abnormal organization density area, overcoming the inherent defect of the sound wave energy attenuation with depth in the traditional sound field processing method. Through accurate allocation of sound field energy in the depth direction, the sound field energy is evenly distributed in each depth layer of the calendering material, significantly improving the uniformity and effectiveness of the sound field processing.

[0104] In an alternative embodiment, each sound field excitation source is driven to exert a synergistic sound field effect on the calendering material according to the phase control parameter and the energy distribution coefficient, the dislocation line orientation distribution of the calendering material under the synergistic sound field effect is collected, the angle distribution between the dislocation line and the calendering direction is counted, and the angle concentration index is calculated, including:

[0105] The phase control parameter is converted into a frequency control quantity of the sound field excitation source, the energy distribution coefficient is converted into a power control quantity of the sound field excitation source, and a sound wave is generated according to the frequency control quantity and the power control quantity;

[0106] The propagation time and pressure amplitude of the collected sound wave in the calendered material are obtained, the propagation time difference and pressure amplitude difference of adjacent sound field excitation sources are calculated, and the propagation characteristics of the sound wave in the depth direction are obtained;

[0107] The beam pointing angle of the sound field excitation source is adjusted according to the propagation characteristics, the beam pointing angle and the pressure amplitude difference are combined to determine the superposition position of the sound wave, a cooperative sound field is formed, the stress state of the calendered material in the action area of the cooperative sound field is collected, the driving force of the dislocation motion is calculated according to the stress state, and the starting point of the dislocation slip is determined;

[0108] The position coordinates of the dislocation line are recorded from the starting point, the spatial trajectory of the dislocation line is constructed according to the position coordinates, the node information of the dislocation line is extracted, the local direction of the dislocation line is calculated according to the node information, the local direction of the dislocation line is projected to the calendering direction, and the dislocation orientation distribution is obtained.

[0109] The included angle values in the dislocation orientation distribution are counted, the number distribution of the included angle values is calculated, and the included angle concentration index is obtained by normalizing the number distribution.

[0110] In the embodiment, the phase control parameter is converted into the frequency control quantity of the sound field excitation source, the energy distribution coefficient is converted into the power control quantity of the sound field excitation source, and the sound wave is generated according to the frequency control quantity and the power control quantity. The phase control parameter is converted into the frequency control quantity of the sound field excitation source through frequency modulation. For the depth layer with a phase deviation of-8.9 degrees, the frequency control quantity is calculated as the basic frequency plus the frequency offset. The frequency offset is equal to the basic frequency multiplied by the phase deviation divided by 360 degrees. Taking a basic frequency of 20 kHz as an example, the frequency offset corresponding to a phase deviation of-8.9 degrees is-0.49 kHz, so the frequency control quantity of the sound field excitation source of this depth layer is 19.51 kHz. The energy distribution coefficient is directly converted into the power control quantity of the sound field excitation source. For the depth layer with an energy distribution coefficient of 0.12, the power control quantity is equal to the total system power multiplied by the coefficient. Taking the total system power of 1000 W as an example, the power control quantity of the sound field excitation source of this depth layer is 120 W. According to the calculated frequency control quantity and power control quantity, a control signal is generated through a digital signal processor to drive an ultrasonic power amplifier to output a corresponding excitation signal, and finally the sound wave is converted by an ultrasonic transducer. For H62 brass calendered material, 8 sound field excitation sources are set, respectively corresponding to different depth layers, and independently controlled according to the respective frequency control quantity and power control quantity.

[0111] The propagation time and pressure amplitude of the sound wave in the rolled material are collected, the propagation time difference and pressure amplitude difference of adjacent sound field excitation sources are calculated, and the propagation characteristics of the sound wave in the depth direction are obtained. An ultrasonic receiving probe array is arranged on the surface of the rolled material to receive sound wave signals from each depth layer. The probe array is composed of 16 piezoelectric sensors uniformly distributed in an area of 100 mm x 100 mm. The collected sound wave signals include two key parameters: propagation time and pressure amplitude. The propagation time is obtained by measuring the time interval from the emission to the reception of the sound wave, and the pressure amplitude is obtained by measuring the voltage peak value of the received signal. For H62 brass rolled material, the sound wave propagation time at a depth of 3 mm is 2.16 microseconds, and the pressure amplitude is 0.85 MPa; the sound wave propagation time at a depth of 4 mm is 2.87 microseconds, and the pressure amplitude is 0.72 MPa. The propagation time difference and pressure amplitude difference between adjacent depth layers are calculated, and the propagation time difference between 3-4 mm depth layers is 0.71 microseconds, and the pressure amplitude difference is 0.13 MPa. Through data analysis of all depth layers, the propagation characteristic curve of the sound wave in the depth direction is obtained, which represents the variation law of sound wave velocity and attenuation with depth.

[0112] The beam pointing angle of the sound field excitation source is adjusted according to the propagation characteristics, the beam pointing angle and the pressure amplitude difference are combined to determine the sound wave superposition position, a cooperative sound field is formed, the stress state of the rolled material in the cooperative sound field action area is collected, the driving force of dislocation motion is calculated according to the stress state, and the starting point of dislocation slip is determined. The beam pointing angle of the sound field excitation source is adjusted through phased array technology. Based on the propagation characteristics of the sound wave in the depth direction, the required beam pointing angle of each depth layer is calculated. For the depth layer that needs sound field enhancement, the beam pointing angle is adjusted to the direction of maximum sound wave energy. For the 3-4 mm depth layer in the H62 brass rolled material, the beam pointing angle is adjusted to 5.3 degrees; for the 8-9 mm depth layer, the beam pointing angle is adjusted to 7.8 degrees. The beam pointing angle and the pressure amplitude difference are combined to determine the optimal sound wave superposition position. The sound wave superposition position is selected in the area where the pressure amplitude difference is maximum and the beam pointing angle is appropriate, so as to realize the maximum focusing of sound field energy. In the determined cooperative sound field action area, the in-situ transmission electron microscope is used to observe the microstructure change of the rolled material, and the stress state data is collected. The stress state is obtained by measuring the degree of lattice distortion. The greater the degree of lattice distortion, the higher the local stress. According to the stress state data, the driving force of dislocation motion is calculated. The driving force of dislocation motion is equal to the applied shear stress minus the critical shear stress. When the driving force is greater than zero, the dislocation begins to slip. For the H62 brass rolled material, the measured shear stress at a depth of 4 mm is 58 MPa, the critical shear stress is 45 MPa, and the calculated driving force of dislocation motion is 13 MPa, which determines the position as the starting point of dislocation slip.

[0113] The position coordinates of the dislocation line are recorded from the starting point, the spatial trajectory of the dislocation line is constructed according to the position coordinates, the node information of the dislocation line is extracted, the local direction of the dislocation line is calculated according to the node information, the local direction of the dislocation line is projected to the calendering direction, and the dislocation orientation distribution is obtained. Using a high-resolution transmission electron microscope, the position coordinates are recorded every 50 nm along the extension direction of the dislocation line from the starting point of dislocation slip to form a spatial trajectory dataset of the dislocation line. For a typical dislocation line in the H62 brass calendering material, 20 position coordinate points are recorded to form a complete spatial trajectory. A three-dimensional spatial curve is constructed according to the position coordinates, and the node information of the dislocation line is extracted. The node information includes node position, curvature and tangent direction. For the recorded dislocation line trajectory, the plane determined by the adjacent three coordinate points is calculated, and the normal vector of the plane is the local direction of the dislocation line at the node. The local direction of the dislocation line is projected to the calendering direction, the included angle between the two is calculated, and the dislocation orientation distribution is obtained. For the observed dislocation line, the included angle between the local direction and the calendering direction is distributed between 15 degrees and 75 degrees, and the average included angle is 42 degrees.

[0114] The included angle values in the dislocation orientation distribution are counted, the number distribution of the included angle values is calculated, and the included angle concentration index is obtained by normalizing the number distribution. The included angle between the dislocation line and the calendering direction is divided into 9 intervals: 0-10 degrees, 10-20 degrees,..., 80-90 degrees, the number of dislocation line nodes in each interval is counted, and a number distribution histogram of the included angle values is formed. For the H62 brass calendering material, the included angle distribution shows a bimodal characteristic, mainly concentrated in the 30-40 degree and 60-70 degree intervals, accounting for 28% and 23% of the total number of nodes respectively. The number distribution is normalized to calculate the included angle concentration index. The included angle concentration index is defined as the proportion of the number of nodes in the dominant included angle interval to the total number of nodes. Among them, the 30-40 degree interval is determined as the dominant included angle interval, and the included angle concentration index is 0.28. The higher the included angle concentration index, the more concentrated the dislocation line orientation, and the stronger the anisotropy of the material structure; the lower the index, the more uniform the dislocation line orientation distribution, and the better the isotropy of the material structure. The change of the included angle concentration index before and after the action of the synergistic sound field can evaluate the effectiveness of the sound field treatment on the organization regulation of the calendering material.

[0115] In this embodiment, precise regulation of the microstructure of the calendering material can be achieved, especially effective adjustment of the dislocation line orientation distribution. Through precise control of the phase and energy, a synergistic sound field in the depth direction is formed, effectively overcoming the defect of uneven energy attenuation in the depth direction of a single sound field. Under the action of the sound field, the slip and rearrangement of the dislocation line are directionally regulated, making the dislocation line orientation tend to be reasonably distributed, effectively reducing the stress concentration in the calendering material and improving the uniformity of the material structure and the stability of the performance.

[0116] In an optional embodiment, a target enhancement coefficient of dislocation orientation is determined based on a target calendering performance index, a sound wave superposition gain is calculated according to an included angle concentration index, the sound wave superposition gain is compared with the target enhancement coefficient, and a phase difference and an energy ratio between adjacent sound field excitation sources are optimized based on a comparison result to realize collaborative regulation of dislocation orientation of the calendering material, including:

[0117] A stress requirement and a texture requirement in the target calendering performance index are obtained, a crystal orientation distribution and a dislocation slip distribution of the calendering material are calculated, a critical shear force of dislocation movement is determined according to the crystal orientation distribution and the dislocation slip distribution, and the critical shear force is converted into a target enhancement coefficient of dislocation orientation;

[0118] Sound wave propagation signals of adjacent sound field excitation sources are collected, a sound wave propagation time difference and a pressure amplitude difference are calculated, an included angle value distribution in the included angle concentration index is extracted, a dynamic response function of the included angle value distribution and the sound wave propagation signals is established, an output power of the adjacent sound field excitation sources is adjusted according to the dynamic response function, and the adjusted sound wave propagation signals are converted into a sound wave superposition gain;

[0119] A difference between the sound wave superposition gain and the target enhancement coefficient is taken as a regulation deviation of dislocation orientation, a phase compensation value and a power distribution ratio of the adjacent sound field excitation sources are calculated according to the regulation deviation, output parameters of the sound field excitation sources are adjusted based on the phase compensation value and the power distribution ratio, and collaborative regulation of dislocation orientation of the calendering material is realized.

[0120] In the rolling process, the stress requirement and texture requirement in the target rolling performance index are obtained, the crystal orientation distribution and dislocation slip distribution of the rolled material are calculated, the critical shear stress of dislocation motion is determined according to the crystal orientation distribution and dislocation slip distribution, and the target enhancement coefficient of dislocation orientation is converted from the critical shear stress. The target rolling performance index is usually determined by the product application requirement, including mechanical performance indexes such as tensile strength, yield strength and elongation, and characteristic indexes such as conductivity and texture degree. For the C2680 brass rolling strip for electronic use, the target tensile strength is 430 MPa, the yield strength is 340 MPa, the elongation is not less than 22%, the conductivity is not less than 28% IACS, and the texture requirement is that the proportion of cubic texture is not less than 65%. The crystal orientation distribution function is obtained by analyzing the crystal orientation of the rolled material by electron backscatter diffraction technology. For the C2680 brass rolling strip, the crystal orientation distribution is mainly concentrated in two texture components of {110} <112> and {112} <111>, accounting for 45% and 25% respectively. According to the crystal orientation distribution, the possible slip system is determined, and the dislocation slip distribution is calculated. The dislocation slip distribution refers to the distribution of dislocation density on different slip systems. Using dislocation statistical model analysis, it is obtained that the dislocation density on the main slip system is 8*10^10 / m², and the dislocation density on the secondary slip system is 3*10^10 / m². According to the crystal orientation distribution and dislocation slip distribution, the critical shear stress of dislocation motion is calculated by combining the Schmidt factor rule, and the critical shear stress represents the minimum stress required for dislocation to start moving. For the C2680 brass rolling strip, the critical shear stress of the main slip system is 68 MPa. The target enhancement coefficient of dislocation orientation is converted from the critical shear stress, and the conversion relationship is established based on the dislocation kinematics model. The target enhancement coefficient of dislocation orientation is defined as the growth ratio of dislocation density on the main slip system after the action of acoustic field relative to the original state. For the case of critical shear stress of 68 MPa, the target enhancement coefficient of dislocation orientation is calculated to be 1.35, which means that the dislocation density on the main slip system needs to be increased by 35%.

[0121] The sound wave propagation signals of adjacent sound field excitation sources are collected, the sound wave propagation time difference and the pressure amplitude difference are calculated, the included angle value distribution in the included angle concentration index is extracted, the dynamic response function of the included angle value distribution and the sound wave propagation signal is established, the output power of the adjacent sound field excitation source is adjusted according to the dynamic response function, and the adjusted sound wave propagation signal is converted into sound wave superposition gain. The piezoelectric sensor array is used to collect the sound wave propagation signals of adjacent sound field excitation sources, the sampling frequency is set to 100 MHz, and the sampling time is 200 μs. For the C2680 brass rolled strip, the adjacent sound field excitation sources arranged at the depths of 4 mm and 5 mm, the collected sound wave propagation times are 2.63 μs and 3.28 μs respectively, and the calculated propagation time difference is 0.65 μs; the collected pressure amplitudes are 0.78 MPa and 0.65 MPa respectively, and the calculated pressure amplitude difference is 0.13 MPa. The transmission electron microscope is used to observe the dislocation line orientation, and the included angle value distribution between the dislocation line and the rolling direction is counted. For the C2680 brass rolled strip, the included angle values are mainly concentrated in the two intervals of 25°-35° and 55°-65°, and the included angle concentration index is 0.32, indicating that the dislocations in the dominant included angle interval account for 32% of the total dislocations. The dynamic response function of the included angle value distribution and the sound wave propagation signal is established, which describes the influence law of the sound field effect on the dislocation orientation distribution. The dynamic response function adopts a nonlinear regression model, the input variable is the sound wave propagation time difference and the pressure amplitude difference, and the output variable is the change rate of the included angle concentration index. According to the established dynamic response function, the output power of the adjacent sound field excitation source is adjusted. When the included angle concentration index is lower than the target value, the output power of the sound field excitation source corresponding to the depth layer in the dominant included angle interval is increased; when the index is higher than the target value, the output power of the sound field excitation source corresponding to the depth layer in the non-dominant interval is increased. For the C2680 brass rolled strip, the output power of the sound field excitation source at the depth of 4 mm is adjusted from the original 180 W to 210 W, and the output power of the sound field excitation source at the depth of 5 mm is adjusted from 150 W to 135 W. The adjusted sound wave propagation signal is converted into sound wave superposition gain through the energy superposition model. The sound wave superposition gain is defined as the ratio of the adjusted sound field energy to the original sound field energy, and the calculation result is 1.28, indicating that the sound field energy is increased by 28%.

[0122] The difference between the sound wave superposition gain and the target enhancement coefficient is taken as a regulation deviation of dislocation orientation, the phase compensation value and the power distribution ratio of adjacent sound field excitation sources are calculated according to the regulation deviation, the output parameters of the sound field excitation sources are adjusted based on the phase compensation value and the power distribution ratio, and the synergistic regulation of dislocation orientation is realized. The difference between the sound wave superposition gain 1.28 and the target enhancement coefficient 1.35 is 0.07, which is the regulation deviation of dislocation orientation. A positive regulation deviation indicates that the current sound field intensity is insufficient and needs to be further enhanced, and a negative regulation deviation indicates that it is too strong and needs to be appropriately weakened. The phase compensation value of adjacent sound field excitation sources is calculated according to the regulation deviation. The phase compensation value is proportional to the regulation deviation, and the proportional coefficient is determined through experiments. For a regulation deviation of 0.07, the calculated phase compensation value is 12.6°. The phase compensation value is used to adjust the phase difference between adjacent sound field excitation sources to optimize the sound wave interference effect. The power distribution ratio is calculated according to the regulation deviation, and the power distribution ratio is defined as the ratio of the output power of adjacent sound field excitation sources. For a regulation deviation of 0.07, the calculated power distribution ratio is 1.65, indicating that the output power of the sound field excitation source at a depth of 4mm should be 1.65 times that of the sound field excitation source at a depth of 5mm. Based on the phase compensation value and the power distribution ratio, the output parameters of the sound field excitation sources are adjusted through an ultrasonic wave control system. The phase compensation value of 12.6° is converted into a phase modulation control signal, and the phase adjustment is realized through a digital phase shifter; the power distribution ratio of 1.65 is converted into a power modulation control signal, and the power adjustment is realized through a programmable power amplifier. The sound field excitation sources are driven to work according to the adjusted parameters to form an optimized synergistic sound field. The change of dislocation orientation distribution is monitored in real time through an electron microscope, and when the included angle concentration index reaches the set threshold, the sound field processing is completed. For C2680 brass rolled strip, after 20 minutes of synergistic sound field processing, the included angle concentration index increases from the original 0.32 to 0.44, the dislocations are mainly concentrated in the interval with an included angle of 30°±5° with the rolling direction, the dislocation arrangement is more ordered, the crystal orientation is more consistent, and the synergistic regulation of dislocation orientation is realized.

[0123] In the present embodiment, by establishing the correspondence between the rolling performance index and the dislocation orientation distribution, the macroscopic performance requirement is converted into the microscopic organization regulation target, and the precise control of the sound field parameters is used to realize the precise matching of the target enhancement coefficient. Compared with the traditional heat treatment or mechanical processing method, the internal microstructure of the material can be adjusted directionally without changing the chemical composition and surface state of the material. Through the synergistic optimization of the phase difference and the energy ratio, the limitations of single sound field regulation are effectively overcome, and the utilization efficiency and regulation effect of the sound field energy are significantly improved.

[0124] In an alternative embodiment, a dynamic response function of the included angle value distribution and the sound wave propagation signal is established, and the output power of adjacent sound field excitation sources is adjusted according to the dynamic response function.

[0125] The peak position and the distribution width of the included angle value are obtained by statistical analysis of the included angle value distribution data, the peak position is taken as the dominant direction parameter of the dislocation orientation, and the distribution width is taken as the dispersion degree parameter of the dislocation orientation;

[0126] The sound wave propagation signals inside the calendering material of adjacent sound field excitation sources are collected, the sound wave arrival time and the sound wave peak amplitude are extracted by time domain analysis of the sound wave propagation signals, and the sound wave propagation time difference and the sound wave peak amplitude difference between adjacent sound field excitation sources are calculated;

[0127] The dominant direction parameter and the dispersion degree parameter are set as output variables of the dynamic response function, the sound wave propagation time difference and the sound wave peak amplitude difference are set as input variables of the dynamic response function, and the dynamic response function of the included angle value distribution and the sound wave propagation signal is constructed by fitting the functional relationship between the input variables and the output variables;

[0128] The partial derivatives of the dynamic response function with respect to the sound wave propagation time difference and the sound wave peak amplitude difference are calculated, the sound wave propagation time difference target value and the sound wave peak amplitude difference target value are converted into the output power adjustment amount of the adjacent sound field excitation source according to the partial derivatives, the output power adjustment amount is applied to the corresponding sound field excitation source, and the adjustment of the output power of the adjacent sound field excitation source is completed.

[0129] The peak position and the distribution width of the included angle value are obtained by statistical analysis of the included angle value distribution data, the peak position is taken as the dominant direction parameter of the dislocation orientation, and the distribution width is taken as the dispersion degree parameter of the dislocation orientation. 1500 included angle data are obtained by collecting the orientation data of 300 dislocation lines from different regions of the calendering material, and measuring 5 points of each dislocation line. The included angle values in the range of 0° to 90° are divided into 18 intervals, the data frequency in each interval is counted, and a distribution histogram is drawn. The peak position and the distribution width of the included angle distribution are extracted by using the kernel density estimation method to smooth the distribution. The included angle distribution presents a double-peak characteristic with a main peak position at 32° and a secondary peak position at 58°, the main peak distribution width is 7.5°, the secondary peak distribution width is 10.2°, and the main peak data accounts for 42% of the total data. The main peak position 32° is determined as the dominant direction parameter of the dislocation orientation, and its distribution width 7.5° is determined as the dispersion degree parameter of the dislocation orientation. These two parameters comprehensively describe the orientation characteristics of the dislocation lines in the calendering material.

[0130] The sound wave propagation signals inside the calendered material of adjacent sound field excitation sources were collected. The sound wave arrival time and the sound wave peak amplitude were extracted by time domain analysis of the sound wave propagation signals. The sound wave propagation time difference and the sound wave peak amplitude difference between adjacent sound field excitation sources were calculated. Eight piezoelectric sensors were arranged on the surface of the calendered material in a 5mm x 5mm grid to receive the sound wave signals emitted by adjacent sound field excitation sources with depths of 3mm and 6mm respectively. The collected sound wave signals were filtered to remove environmental noise and high-frequency interference, and the signal envelope was extracted using Hilbert transform. The time when the amplitude of the sound wave signal is 3 times the background noise is defined as the sound wave arrival time, and the maximum value of the signal envelope is the sound wave peak amplitude. For the 3mm deep sound field excitation source, the sound wave arrival time is 1.86μs and the peak amplitude is 0.87MPa; for the 6mm deep sound field excitation source, the sound wave arrival time is 3.72μs and the peak amplitude is 0.63MPa. The sound wave propagation time difference between adjacent sound field excitation sources is 1.86μs, and the sound wave peak amplitude difference is 0.24MPa. These parameters reflect the space-time characteristics of sound wave propagation inside the calendered material.

[0131] The advantage direction parameter and the dispersion degree parameter were set as the output variables of the dynamic response function, and the sound wave propagation time difference and the sound wave peak amplitude difference were set as the input variables of the dynamic response function. The dynamic response function of the angle value distribution and the sound wave propagation signal was constructed by fitting the functional relationship between the input variables and the output variables. To establish the relationship between the input and output variables, 20 sets of sound field processing experiments with different parameter combinations were designed. The frequency difference and power ratio of adjacent sound field excitation sources were changed to produce different sound wave propagation time differences (1.2μs to 2.4μs) and sound wave peak amplitude differences (0.15MPa to 0.35MPa) while keeping other conditions consistent. The dislocation orientation distribution of the samples after each experimental treatment was measured, and the changes in the advantage direction parameter and the dispersion degree parameter were recorded. The dynamic response function was established using multivariate regression analysis method. Considering the possible nonlinear relationship between variables, a quadratic polynomial model was used for fitting. The fitted dynamic response function shows that when the sound wave propagation time difference increases, the advantage direction parameter approaches the calendering direction (the value decreases), and the dispersion degree parameter decreases, indicating that the dislocation orientation is more concentrated; when the sound wave peak amplitude difference increases, the advantage direction parameter moves away from the calendering direction (the value increases), and the dispersion degree parameter increases, indicating that the dislocation orientation distribution is more dispersed. The goodness of model fitting is 0.89, indicating that the established dynamic response function can well describe the relationship between sound field parameters and dislocation orientation distribution.

[0132] The partial derivatives of the dynamic response function of the sound wave propagation time difference and the sound wave peak amplitude difference are calculated, the sound wave propagation time difference target value and the sound wave peak amplitude difference target value are converted into the output power adjustment amount of the adjacent sound field excitation source according to the partial derivatives, the output power adjustment amount is applied to the corresponding sound field excitation source, and the output power adjustment of the adjacent sound field excitation source is completed. Based on the established dynamic response function, the partial derivatives of the function at the current working point (1.86 μs, 0.24 MPa) with respect to the two input variables are calculated. The calculation results show that the partial derivative of the advantage direction parameter with respect to the sound wave propagation time difference is -8.7° / μs, and the partial derivative of the advantage direction parameter with respect to the sound wave peak amplitude difference is 12.3° / MPa; the partial derivative of the dispersion degree parameter with respect to the sound wave propagation time difference is -3.2° / μs, and the partial derivative of the dispersion degree parameter with respect to the sound wave peak amplitude difference is 5.6° / MPa. According to the product performance requirements, the target parameters of the dislocation orientation are determined: the advantage direction parameter is 25°, and the dispersion degree parameter is 5°. There is a difference between the current measured value and the target value: the difference of the advantage direction parameter is 7°, and the difference of the dispersion degree parameter is 2.5°. By solving the equation set, the parameter differences are converted into the required sound wave propagation time difference adjustment amount and the sound wave peak amplitude difference adjustment amount. It is calculated that the sound wave propagation time difference needs to be increased by 0.42 μs, and the sound wave peak amplitude difference needs to be reduced by 0.13 MPa. There is a corresponding relationship between the sound wave propagation parameters and the output power of the sound field excitation source: the propagation time difference is mainly controlled by the frequency difference, and the peak amplitude difference is mainly controlled by the power ratio. Through the sound wave propagation model calculation, it is necessary to increase the output power of the 3mm depth sound field excitation source from the original 220W to 265W, and to decrease the output power of the 6mm depth sound field excitation source from the original 180W to 160W. The power adjustment amount calculated by the ultrasonic control system is applied to the corresponding sound field excitation source, the change of the sound wave propagation parameter is monitored in real time until the target value is reached, and the accurate adjustment of the output power of the adjacent sound field excitation source is completed.

[0133] In the embodiment, the output power of the sound field excitation source is adjusted by the method of establishing the dynamic response function, which provides a scientific basis and implementation path for the modulation control of the copper alloy calendering material in multiple sound field cooperation. Through the fine analysis of the sound wave propagation characteristics and the scientific construction of the dynamic response function, the influence of the sound field parameter change on the dislocation orientation can be accurately predicted, and the output power configuration of the sound field excitation source can be optimized accordingly. The optimized and controlled sound field action can effectively adjust the slip direction and arrangement of the dislocation lines, so that the dislocation orientation is more reasonable and the internal stress distribution is more uniform.

[0134] In a second aspect of the embodiment of the present application, a copper alloy calendering material modulation control system in multiple sound field cooperation is provided, and the system comprises:

[0135] The first unit is used for acquiring the organization state information of the copper alloy material to be calendered and the target calendering performance index.

[0136] The second unit is configured to divide the calendering material into a plurality of sound field action sub-zones in the thickness direction according to the organization state information, and arrange sound field excitation sources in each sub-zone;

[0137] The third unit is configured to detect the propagation time of each sound field excitation source at different depth layers, calculate the propagation time difference, and convert the propagation time difference into a phase control parameter;

[0138] The fourth unit is configured to extract an organization density distribution curve of the calendering material, establish an inverse mapping relationship between the organization density distribution curve and the output energy of the sound field excitation source, and calculate an energy distribution coefficient of each depth layer;

[0139] The fifth unit is configured to drive each sound field excitation source to exert a cooperative sound field action on the calendering material according to the phase control parameter and the energy distribution coefficient, collect dislocation line orientation distribution of the calendering material under the cooperative sound field action, count the angle distribution between the dislocation line and the calendering direction, and calculate an angle concentration index;

[0140] The sixth unit is configured to determine a target enhancement coefficient of the dislocation orientation based on a target calendering performance index, calculate a sound wave superposition gain according to the angle concentration index, compare the sound wave superposition gain with the target enhancement coefficient, and optimize the phase difference and the energy ratio between adjacent sound field excitation sources based on the comparison result, so as to realize cooperative regulation of the dislocation orientation of the calendering material.

[0141] In a third aspect, an electronic device is provided, including:

[0142] a processor;

[0143] a memory for storing processor-executable instructions;

[0144] The processor is configured to invoke the instructions stored in the memory to execute the method described above.

[0145] In a fourth aspect, a computer-readable storage medium is provided, which stores computer program instructions, and the computer program instructions are executed by a processor to implement the method described above.

[0146] The present application can be a method, device, system and / or computer program product. The computer program product can include a computer readable storage medium having computer readable program instructions loaded thereon for performing various aspects of the present application.

[0147] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present application, and are not intended to limit the present application; although the present application has been described in detail with reference to the above embodiments, those skilled in the art should understand that the technical solutions recorded in the above embodiments can be modified, or some or all of the technical features can be replaced by equivalents; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the scope of the technical solutions of the embodiments of the present application.

Claims

1. A method for controlling the temper of a copper alloy rolled material in conjunction with a plurality of sound fields, characterized by, The method comprises the following steps: obtaining the organizational state information of the copper alloy material to be calendered and target calendering performance indicators; dividing the calendering material along the thickness direction into multiple sound field action sub-zones according to the organizational state information, and arranging sound field excitation sources in each sub-zone; detecting the propagation time of each sound field excitation source at different depth layers, calculating the propagation time difference, and converting the propagation time difference into a phase control parameter; extracting the organizational density distribution curve of the calendering material, establishing a reverse mapping relationship between the organizational density distribution curve and the output energy of the sound field excitation source, and calculating the energy distribution coefficient of each depth layer; driving each sound field excitation source to apply a cooperative sound field action to the calendering material according to the phase control parameter and the energy distribution coefficient, collecting the dislocation line orientation distribution of the calendering material under the cooperative sound field action, counting the angle distribution between the dislocation line and the calendering direction, and calculating an angle concentration index; determining a target enhancement coefficient of the dislocation orientation based on the target calendering performance indicators, calculating a sound wave superposition gain according to the angle concentration index, comparing the sound wave superposition gain with the target enhancement coefficient, optimizing the phase difference and energy ratio between adjacent sound field excitation sources based on the comparison result, and realizing the cooperative control of the dislocation orientation of the calendering material.

2. The method of claim 1, wherein, The method of dividing the calendering material along the thickness direction into multiple sound field action sub-zones according to the organizational state information and arranging sound field excitation sources in each sub-zone comprises the following steps: extracting the organizational parameters of the calendering material at different thickness positions according to the organizational state information, calculating the organizational difference value between adjacent thickness positions, determining the organizational mutation position along the thickness direction of the calendering material according to the organizational difference value, dividing the calendering material into multiple sound field action sub-zones based on the organizational mutation position, and calculating the organizational characteristic value of each sound field action sub-zone; calculating the number of sound field excitation sources required in each sound field action sub-zone based on the organizational characteristic value and the size of the sound field action sub-zone, conducting a sound field superposition test on the sound field excitation sources, obtaining the action radius of the sound field excitation sources, and determining the arrangement spacing of the sound field excitation sources according to the action radius and the number of sound field excitation sources; calculating the arrangement coordinates of the sound field excitation sources and the corresponding sound field intensity distribution according to the arrangement spacing and the size of the sound field action sub-zone, selecting the arrangement coordinates that satisfy a preset sound field coverage threshold, and arranging the sound field excitation sources at the selected arrangement coordinates.

3. The method of claim 1, wherein, The method of detecting the propagation time of each sound field excitation source at different depth layers, calculating the propagation time difference, and converting the propagation time difference into a phase control parameter comprises the following steps: collecting the sound wave propagation signals of each sound field excitation source at different depth layers, and extracting the sound wave frequency drift and the sound wave amplitude attenuation rate through time-frequency domain joint analysis of the sound wave propagation signals; inverting the dislocation density distribution of each depth layer through the sound wave frequency drift, inverting the grain boundary density distribution of each depth layer through the sound wave amplitude attenuation rate, calculating the acoustic impedance correction value of each depth layer based on the dislocation density distribution and the grain boundary density distribution, and correcting the propagation time of each depth layer according to the acoustic impedance correction value to obtain the corrected propagation time; calculating the difference between the corrected propagation times to obtain the propagation time difference, performing phase conversion on the propagation time difference and the working frequency of each sound field excitation source to obtain an excitation phase difference value, and constructing a space-time distribution map of the excitation phase difference value varying with the depth layer position; extract a spatial gradient vector field of the excitation phase difference values from the space-time distribution map, identify a depth layer position with a maximum gradient amplitude in the spatial gradient vector field, and take an excitation phase difference value of a sound field excitation source corresponding to the depth layer position with the maximum gradient amplitude as a reference phase difference value; calculate a phase deviation between an excitation phase difference value of another depth layer sound field excitation source and the reference phase difference value as a phase control parameter.

4. The method of claim 1, wherein, extract a structure density distribution curve of the calendered material, establish an inverse mapping relationship between the structure density distribution curve and an output energy of the sound field excitation source, and calculate an energy distribution coefficient of each depth layer, including: obtain structure density values of different depth layers of the calendered material along a thickness direction, arrange the structure density values according to depth layer positions, and construct a structure density distribution curve; calculate a change rate of the structure density values between adjacent depth layers, identify a depth layer with a change rate deviating from a mean value and mark the depth layer as a structure density abnormal region; read a structure density value of the structure density abnormal region, establish a sound wave propagation impedance distribution, deduce a sound wave propagation path according to the sound wave propagation impedance distribution, and calculate a sound wave energy cumulative attenuation amount; take the sound wave energy cumulative attenuation amount as an output energy increment requirement of the sound field excitation source of each depth layer, establish an inverse mapping relationship between the structure density distribution curve and the output energy of the sound field excitation source; read the output energy increment requirement of the sound field excitation source of each depth layer, and construct a depth transmission link of the output energy increment requirement of the sound field excitation source; identify an anchor depth layer with a mutation in the depth transmission link, and bidirectionally diffuse and correct the output energy increment requirement of the sound field excitation source of adjacent depth layers from the anchor depth layer as a starting point, normalize the corrected output energy increment requirement of the sound field excitation source, and obtain the energy distribution coefficient of each depth layer.

5. The method of claim 1, wherein, drive each sound field excitation source to exert a cooperative sound field action on the calendered material according to the phase control parameter and the energy distribution coefficient, collect dislocation line orientation distribution of the calendered material under the cooperative sound field action, count an angle distribution between the dislocation line and the calendering direction, and calculate an angle concentration index, including: convert the phase control parameter into a frequency control amount of the sound field excitation source, convert the energy distribution coefficient into a power control amount of the sound field excitation source, and generate a sound wave according to the frequency control amount and the power control amount; collect a propagation time and a pressure amplitude of the sound wave in the calendered material, calculate a propagation time difference and a pressure amplitude difference between adjacent sound field excitation sources, and obtain a propagation characteristic of the sound wave in the depth direction; adjust a beam pointing angle of the sound field excitation source according to the propagation characteristic, determine a sound wave superposition position by combining the beam pointing angle and the pressure amplitude difference, form a cooperative sound field, collect a stress state of the calendered material in an action area of the cooperative sound field, calculate a driving force of dislocation motion according to the stress state, and determine a starting point of dislocation slip; record a position coordinate of the dislocation line from the starting point, construct a spatial trajectory of the dislocation line according to the position coordinate, extract node information of the dislocation line, calculate a local direction of the dislocation line according to the node information, project the local direction of the dislocation line to the calendering direction, and obtain a dislocation orientation distribution; The angle value in the dislocation orientation distribution is counted, the number distribution of the angle value is calculated, and the number distribution is normalized to obtain the angle concentration index.

6. The method of claim 1, wherein, Based on the target calendering performance index, a target enhancement coefficient of dislocation orientation is determined, a sound wave superposition gain is calculated according to the angle concentration index, the sound wave superposition gain is compared with the target enhancement coefficient, and based on the comparison result, the phase difference and energy ratio between adjacent sound field excitation sources are optimized to realize the synergistic regulation of the dislocation orientation of the calendering material, including: The stress requirement and texture requirement in the target calendering performance index are obtained, the crystal orientation distribution and dislocation slip distribution of the calendering material are calculated, the critical shear stress of dislocation movement is determined according to the crystal orientation distribution and dislocation slip distribution, and the critical shear stress is converted into a target enhancement coefficient of dislocation orientation; The sound wave propagation signals of adjacent sound field excitation sources are collected, the sound wave propagation time difference and pressure amplitude difference are calculated, the angle value distribution in the angle concentration index is extracted, the dynamic response function of the angle value distribution and the sound wave propagation signal is established, the output power of the adjacent sound field excitation sources is adjusted according to the dynamic response function, and the adjusted sound wave propagation signal is converted into a sound wave superposition gain; The difference between the sound wave superposition gain and the target enhancement coefficient is taken as a regulation deviation of the dislocation orientation, the phase compensation value and the power distribution ratio of the adjacent sound field excitation sources are calculated according to the regulation deviation, and the output parameters of the sound field excitation sources are adjusted based on the phase compensation value and the power distribution ratio to realize the synergistic regulation of the dislocation orientation of the calendering material.

7. The method of claim 6, wherein, The dynamic response function of the angle value distribution and the sound wave propagation signal is established, and the output power of the adjacent sound field excitation sources is adjusted according to the dynamic response function, including: Statistical analysis is performed on the angle value distribution data to obtain the peak position and distribution width of the angle value, the peak position is taken as the dominant direction parameter of the dislocation orientation, and the distribution width is taken as the dispersion degree parameter of the dislocation orientation; The sound wave propagation signals of adjacent sound field excitation sources inside the calendering material are collected, the sound wave arrival time and the sound wave peak amplitude are extracted by time domain analysis of the sound wave propagation signal, and the sound wave propagation time difference and the sound wave peak amplitude difference between adjacent sound field excitation sources are calculated; The dominant direction parameter and the dispersion degree parameter are set as the output variables of the dynamic response function, the sound wave propagation time difference and the sound wave peak amplitude difference are set as the input variables of the dynamic response function, and the dynamic response function of the angle value distribution and the sound wave propagation signal is constructed by fitting the functional relationship between the input variables and the output variables; The partial derivative of the dynamic response function with respect to the sound wave propagation time difference and the sound wave peak amplitude difference is calculated, the sound wave propagation time difference target value and the sound wave peak amplitude difference target value are converted into the output power adjustment amount of the adjacent sound field excitation sources according to the partial derivative, and the output power adjustment amount is applied to the corresponding sound field excitation source to complete the adjustment of the output power of the adjacent sound field excitation sources.

8. A multi-acoustic field coordinated copper alloy calendering material conditioning control system for implementing the method of any of the preceding claims 1-7, characterized by, including: A first unit for obtaining the organizational state information of the copper alloy material to be calendered and the target calendering performance index; A second unit for dividing the calendering material into multiple sound field action sub-zones along the thickness direction according to the organizational state information, and arranging sound field excitation sources in each sub-zone; The third unit is configured to detect the propagation time of each sound field excitation source at different depth layers, calculate the propagation time difference, and convert the propagation time difference into a phase control parameter; The fourth unit is configured to extract the distribution curve of the organization density of the calendered material, establish an inverse mapping relationship between the distribution curve of the organization density and the output energy of the sound field excitation source, and calculate the energy distribution coefficient of each depth layer; The fifth unit is configured to drive each sound field excitation source to exert a cooperative sound field effect on the calendered material according to the phase control parameter and the energy distribution coefficient, collect the dislocation line orientation distribution of the calendered material under the cooperative sound field effect, count the angle distribution between the dislocation line and the calendering direction, and calculate an angle concentration index; The sixth unit is configured to determine a target enhancement coefficient of the dislocation orientation based on a target calendering performance index, calculate a sound wave superposition gain according to the angle concentration index, compare the sound wave superposition gain with the target enhancement coefficient, optimize the phase difference and the energy ratio between adjacent sound field excitation sources based on the comparison result, and realize cooperative regulation of the dislocation orientation of the calendered material.

9. An electronic device, comprising: The computer program instructions are executed by the processor to realize the method of any one of claims 1-7. The computer program instructions are executed by the processor to realize the method of any one of claims 1-7. ​ ​ 10. A computer-readable storage medium having stored thereon computer program instructions, wherein, ​

Images