Bi-level hybrid local search approach for three-dimensional loading problem with balancing constraints

This paper presents a bi-level hybrid local search(BHLS)algorithm for the three-dimensional loading problem with balancing constraints(3DLP-B),where several rectangular boxes with even densities but different sizes are loaded into a single cubic bin to meet the requirements of the space or capacity utilization and the balance of the center of gravity.The proposed algorithm hybridizes a novel framed-layout procedure in which the concept of the core block and its generation strategy are introduced.Once the block-loading sequence has been determined,we can load one block at a time by the designed construction heuristic.Then,the double-search is introduced;its external search procedure generates a list of compact packing patterns while its internal search procedure is used to search the core-block frames and their best distribution locations.The approach is extensively tested on weakly to strongly heterogeneous benchmark data.The results show that it has better performance in improving space utilization rate and balanced condition of the placement than existed techniques:the overall averages from 79.85%to 86.45%were obtained for the balanced cases and relatively high space-usage rate of 89.44%was achieved for the unbalanced ones.

A theory for three-dimensional response of micropolar plates

Through combined applications of the transfer-matrix method and asymptotic expansion technique,we formulate a theory to predict the three-dimensional response of micropolar plates.No ad hoc assumptions regarding through-thickness assumptions of the field variables are made,and the governing equations are two-dimensional,with the displacements and microrotations of the mid-plane as the unknowns.Once the deformation of the mid-plane is solved,a three-dimensional micropolar elastic field within the plate is generated,which is exact up to the second order except in the boundary region close to the plate edge.As an illustrative example,the bending of a clamped infinitely long plate caused by a uniformly distributed transverse force is analyzed and discussed in detail.

Oxygen tension modulates cell function in an in vitro three-dimensional glioblastoma tumor model

Hypoxia is a typical feature of the tumor microenvironment,one of the most critical factors affecting cell behavior and tumor progression.However,the lack of tumor models able to precisely emulate natural brain tumor tissue has impeded the study of the effects of hypoxia on the progression and growth of tumor cells.This study reports a three-dimensional(3D)brain tumor model obtained by encapsulating U87MG(U87)cells in a hydrogel containing type I collagen.It also documents the effect of various oxygen concentrations(1%,7%,and 21%)in the culture environment on U87 cell morphology,proliferation,viability,cell cycle,apoptosis rate,and migration.Finally,it compares two-dimensional(2D)and 3D cultures.For comparison purposes,cells cultured in flat culture dishes were used as the control(2D model).Cells cultured in the 3D model proliferated more slowly but had a higher apoptosis rate and proportion of cells in the resting phase(G0 phase)/gap I phase(G1 phase)than those cultured in the 2D model.Besides,the two models yielded significantly different cell morphologies.Finally,hypoxia(e.g.,1%O2)affected cell morphology,slowed cell growth,reduced cell viability,and increased the apoptosis rate in the 3D model.These results indicate that the constructed 3D model is effective for investigating the effects of biological and chemical factors on cell morphology and function,and can be more representative of the tumor microenvironment than 2D culture systems.The developed 3D glioblastoma tumor model is equally applicable to other studies in pharmacology and pathology.

Three-dimensional cell-based strategies for liver regeneration

Liver regeneration and the development of effective therapies for liver failure remain formidable challenges in modern medicine.In recent years,the utilization of 3D cell-based strategies has emerged as a promising approach for addressing these urgent clinical requirements.This review provides a thorough analysis of the application of 3D cell-based approaches to liver regeneration and their potential impact on patients with end-stage liver failure.Here,we discuss various 3D culture models that incorporate hepatocytes and stem cells to restore liver function and ameliorate the consequences of liver failure.Furthermore,we explored the challenges in transitioning these innovative strategies from preclinical studies to clinical applications.The collective insights presented herein highlight the significance of 3D cell-based strategies as a transformative paradigm for liver regeneration and improved patient care.

The combined application of stem cells and three-dimensional bioprinting scaffolds for the repair of spinal cord injury

Spinal cord injury is considered one of the most difficult injuries to repair and has one of the worst prognoses for injuries to the nervous system.Following surgery,the poor regenerative capacity of nerve cells and the generation of new scars can make it very difficult for the impaired nervous system to restore its neural functionality.Traditional treatments can only alleviate secondary injuries but cannot fundamentally repair the spinal cord.Consequently,there is a critical need to develop new treatments to promote functional repair after spinal cord injury.Over recent years,there have been seve ral developments in the use of stem cell therapy for the treatment of spinal cord injury.Alongside significant developments in the field of tissue engineering,three-dimensional bioprinting technology has become a hot research topic due to its ability to accurately print complex structures.This led to the loading of three-dimensional bioprinting scaffolds which provided precise cell localization.These three-dimensional bioprinting scaffolds co uld repair damaged neural circuits and had the potential to repair the damaged spinal cord.In this review,we discuss the mechanisms underlying simple stem cell therapy,the application of different types of stem cells for the treatment of spinal cord injury,and the different manufa cturing methods for three-dimensional bioprinting scaffolds.In particular,we focus on the development of three-dimensional bioprinting scaffolds for the treatment of spinal cord injury.

THREE-DIMENSIONAL ELLIPTIC CRACK UNDER IMPACT LOADING

The dynamic stress intensity factor of a three-dimensionalelliptic crack under impact loading is determined with the finiteelement method. The computation results can take into account theinfluence of time and the ratio of the wave speeds on the stressintensity factor. The present method is suitable not only forthree-dimensional dynamic crack, but also for three-dimensionaldynamic contact.

Numerical study on wave loads and motions of two ships advancing in waves by using three-dimensional translating-pulsating source

A frequency domain analysis method based on the three-dimensional translating-pulsating （3DTP） source Green function is developed to investigate wave loads and free motions of two ships advancing on parallel course in waves. Two experiments are carried out respectively to mea- sure the wave loads and the free motions for a pair of side-by- side arranged ship models advancing with an identical speed in head regular waves. For comparison, each model is also tested alone. Predictions obtained by the present solution are found in favorable agreement with the model tests and are more accurate than the traditional method based on the three dimensional pulsating （3DP） source Green function. Numer- ical resonances and peak shift can be found in the 3DP pre- dictions, which result from the wave energy trapped in the gap between two ships and the extremely inhomogeneous wave load distribution on each hull. However, they can be eliminated by 3DTP, in which the speed affects the free sur- face and most of the wave energy can be escaped from the gap. Both the experiment and the present prediction show that hydrodynamic interaction effects on wave loads and free motions are significant. The present solver may serve as a validated tool to predict wave loads and motions of two ves- sels under replenishment at sea, and may help to evaluate the hydrodynamic interaction effects on the ships safety in replenishment operation.

Study on Failure Mechanism and Bearing Capacity of Three-Dimensional Rectangular Footing Subjected to Combined Loading

This paper presents two kinematic failure mechanisms of threc-dimensional rectangular footing resting on homogeneous undrained clay foundation under uniaxial vertical loading and uniaxial moment loading. The failure mechanism under vertical loading comprises a plane strain Prandti-type mechanism over the central part of the longer side, and the size of the mechanism gradually reduces at the ends of the longer side and over the shorter side as the corner of rectangular footing is being approached where the direction of soil motion remains normal to each corresponding side respectively. The failure mechanism under moment loading comprises a plane strain scoop sliding mechanism over the central part of the longer side, and the radius of scoop sliding mechanism increases linearly at the ends of the longer side. On the basis of the kinematic failure mechanisms mentioned above, the vertical ultimate bearing capacity and the ultimate bearing capacity against moment or moment ultimate bearing capacity are obtained by use of upper bound limit analysis theory. At the same time, numerical analysis results, Skempton＇ s results and Salgado et al. ＇s results are compared with this upper bound solution. It shows that the presented failure mechanisms and plastic limit analysis predictions are validated. In order to investigate the behaviors of undrained clay foundation beneath the rectangular footing subjected to the combined loadings, numerical analysis is adopted by virtue of the general-purpose FEM software ABAQUS, where the clay is assumed to obey the Mohr-Coulomb yielding criterion. The failure envelope and the ultimate bearing capacity are achieved by the numerical analysis results with the varying aspect ratios from length L to breadth B of the rectangular footing. The failure mechanisms of rectangular footing which are subjected to the combined vertical loading V and horizontal loading H （Vertical loading V and moment loading M, and horizontal loading H and moment loading M respectively are observed in the finite element analysis. ） is explained by use of the upper bound plasticity limit analysis theory. Finally, the reason of eccentricity of failure envelope in H-M loading space is given in this study, which can not be explained by use of the traditional ＇ swipe test＇.

Development of a toroidal soft x-ray imaging system and application for investigating three-dimensional plasma on J-TEXT

A toroidal soft x-ray imaging(T-SXRI)system has been developed to investigate threedimensional(3D)plasma physics on J-TEXT.This T-SXRI system consists of three sets of SXR arrays.Two sets are newly developed and located on the vacuum chamber wall at toroidal positionsφof 126.4°and 272.6°,respectively,while one set was established previously atφ=65.50.Each set of SXR arrays consists of three arrays viewing the plasma poloidally,and hence can be used separately to obtain SXR images via the tomographic method.The sawtooth precursor oscillations are measured by T-SXRI,and the corresponding images of perturbative SXR signals are successfully reconstructed at these three toroidal positions,hence providing measurement of the 3D structure of precursor oscillations.The observed 3D structure is consistent with the helical structure of the m/n=1/1 mode.The experimental observation confirms that the T-SXRI system is able to observe 3D structures in the J-TEXT plasma.

Intraoperative application of three-dimensional printed guides in total hip arthroplasty: A systematic review

BACKGROUND Acetabular component positioning in total hip arthroplasty(THA)is of key importance to ensure satisfactory post-operative outcomes and to minimize the risk of complications.The majority of acetabular components are aligned freehand,without the use of navigation methods.Patient specific instruments(PSI)and three-dimensional(3D)printing of THA placement guides are increasingly used in primary THA to ensure optimal positioning.AIM To summarize the literature on 3D printing in THA and how they improve acetabular component alignment.METHODS PubMed was used to identify and access scientific studies reporting on different 3D printing methods used in THA.Eight studies with 236 hips in 228 patients were included.The studies could be divided into two main categories;3D printed models and 3D printed guides.RESULTS 3D printing in THA helped improve preoperative cup size planning and post-operative Harris hip scores between intervention and control groups(P=0.019,P=0.009).Otherwise,outcome measures were heterogeneous and thus difficult to compare.The overarching consensus between the studies is that the use of 3D guidance tools can assist in improving THA cup positioning and reduce the need for revision THA and the associated costs.CONCLUSION The implementation of 3D printing and PSI for primary THA can significantly improve the positioning accuracy of the acetabular cup component and reduce the number of complications caused by malpositioning.

Computer-assisted three-dimensional individualized extreme liver resection for hepatoblastoma in proximity to the major liver vasculature

BACKGROUND The management of hepatoblastoma(HB)becomes challenging when the tumor remains in close proximity to the major liver vasculature(PMV)even after a full course of neoadjuvant chemotherapy(NAC).In such cases,extreme liver resection can be considered a potential option.AIM To explore whether computer-assisted three-dimensional individualized extreme liver resection is safe and feasible for children with HB who still have PMV after a full course of NAC.METHODS We retrospectively collected data from children with HB who underwent surgical resection at our center from June 2013 to June 2023.We then analyzed the detailed clinical and three-dimensional characteristics of children with HB who still had PMV after a full course of NAC.RESULTS Sixty-seven children diagnosed with HB underwent surgical resection.The age at diagnosis was 21.4±18.8 months,and 40 boys and 27 girls were included.Fifty-nine(88.1%)patients had a single tumor,39(58.2%)of which was located in the right lobe of the liver.A total of 47 patients(70.1%)had PRE-TEXT III or IV.Thirty-nine patients(58.2%)underwent delayed resection.After a full course of NAC,16 patients still had close PMV(within 1 cm in two patients,touching in 11 patients,compressing in four patients,and showing tumor thrombus in three patients).There were 6 patients of tumors in the middle lobe of the liver,and four of those patients exhibited liver anatomy variations.These 16 children underwent extreme liver resection after comprehensive preoperative evaluation.Intraoperative procedures were performed according to the preoperative plan,and the operations were successfully performed.Currently,the 3-year event-free survival of 67 children with HB is 88%.Among the 16 children who underwent extreme liver resection,three experienced recurrence,and one died due to multiple metastases.CONCLUSION Extreme liver resection for HB that is still in close PMV after a full course of NAC is both safe and feasible.This approach not only reduces the necessity for liver transplantation but also results in a favorable prognosis.Individualized three-dimensional surgical planning is beneficial for accurate and complete resection of HB,particularly for assessing vascular involvement,remnant liver volume and anatomical variations.

Three-dimensional numerical analysis of plant-soil hydraulic interactions on pore water pressure of vegetated slope under different rainfall patterns

Understanding the pore water pressure distribution in unsaturated soil is crucial in predicting shallow landslides triggered by rainfall,mainly when dealing with different temporal patterns of rainfall intensity.However,the hydrological response of vegetated slopes,especially three-dimensional(3D)slopes covered with shrubs,under different rainfall patterns remains unclear and requires further investigation.To address this issue,this study adopts a novel 3D numerical model for simulating hydraulic interactions between the root system of the shrub and the surrounding soil.Three series of numerical parametric studies are conducted to investigate the influences of slope inclination,rainfall pattern and rainfall duration.Four rainfall patterns(advanced,bimodal,delayed,and uniform)and two rainfall durations(4-h intense and 168-h mild rainfall)are considered to study the hydrological response of the slope.The computed results show that 17%higher transpiration-induced suction is found for a steeper slope,which remains even after a short,intense rainfall with a 100-year return period.The extreme rainfalls with advanced(PA),bimodal(PB)and uniform(PU)rainfall patterns need to be considered for the short rainfall duration(4 h),while the delayed(PD)and uniform(PU)rainfall patterns are highly recommended for long rainfall durations(168 h).The presence of plants can improve slope stability markedly under extreme rainfall with a short duration(4 h).For the long duration(168 h),the benefit of the plant in preserving pore-water pressure(PWP)and slope stability may not be sufficient.

Influence of sampling on three-dimensional surface shape measurement

In order to accurately measure an object’s three-dimensional surface shape,the influence of sampling on it was studied.First,on the basis of deriving spectra expressions through the Fourier transform,the generation of CCD pixels was analyzed,and its expression was given.Then,based on the discrete expression of deformation fringes obtained after sampling,its Fourier spectrum expression was derived,resulting in an infinitely repeated"spectra island"in the frequency domain.Finally,on the basis of using a low-pass filter to remove high-order harmonic components and retaining only one fundamental frequency component,the inverse Fourier transform was used to reconstruct the signal strength.A method of reducing the sampling interval,i.e.,reducing the number of sampling points per fringe,was proposed to increase the ratio between the sampling frequency and the fundamental frequency of the grating.This was done to reconstruct the object’s surface shape more accurately under the condition of m>4.The basic principle was verified through simulation and experiment.In the simulation,the sampling intervals were 8 pixels,4 pixels,2 pixels,and 1 pixel,the maximum absolute error values obtained in the last three situations were 88.80%,38.38%,and 31.50%in the first situation,respectively,and the corresponding average absolute error values are 71.84%,43.27%,and 32.26%.It is demonstrated that the smaller the sampling interval,the better the recovery effect.Taking the same four sampling intervals in the experiment as in the simulation can also lead to the same conclusions.The simulated and experimental results show that reducing the sampling interval can improve the accuracy of object surface shape measurement and achieve better reconstruction results.

A HALF PLANE CRACK UNDER THREE-DIMENSIONAL COMBINED MODE IMPACT LOADING

The dynamic stress intensity factor history for a half plane crack in an otherwise unbounded elastic body, with the crack faces subjected to a traction distribution consisting of two pairs of suddenly-applied shear line loads is consid- ered. The analytic expression for the combined mode stress intensity factors as a function of time is obtained. The method of solution is based on the application of integral transforms and the Wiener-Hopf technique. Some features of the solutions are discussed and graphical numerical results are presented.

Three-dimensional visualization technology for guiding one-step percutaneous transhepatic cholangioscopic lithotripsy for the treatment of complex hepatolithiasis

BACKGROUND Biliary stone disease is a highly prevalent condition and a leading cause of hospitalization worldwide.Hepatolithiasis with associated strictures has high residual and recurrence rates after traditional multisession percutaneous transhepatic cholangioscopic lithotripsy(PTCSL).AIM To study one-step PTCSL using the percutaneous transhepatic one-step biliary fistulation(PTOBF)technique guided by three-dimensional(3D)visualization.METHODS This was a retrospective,single-center study analyzing,140 patients who,between October 2016 and October 2023,underwent one-step PTCSL for hepatolithiasis.The patients were divided into two groups:The 3D-PTOBF group and the PTOBF group.Stone clearance on choledochoscopy,complications,and long-term clearance and recurrence rates were assessed.RESULTS Age,total bilirubin,direct bilirubin,Child-Pugh class,and stone location were similar between the 2 groups,but there was a significant difference in bile duct strictures,with biliary strictures more common in the 3D-PTOBF group(P=0.001).The median follow-up time was 55.0(55.0,512.0)days.The immediate stone clearance ratio(88.6%vs 27.1%,P=0.000)and stricture resolution ratio(97.1%vs 78.6%,P=0.001)in the 3D-PTOBF group were significantly greater than those in the PTOBF group.Postoperative complication(8.6%vs 41.4%,P=0.000)and stone recurrence rates(7.1%vs 38.6%,P=0.000)were significantly lower in the 3D-PTOBF group.CONCLUSION Three-dimensional visualization helps make one-step PTCSL a safe,effective,and promising treatment for patients with complicated primary hepatolithiasis.The perioperative and long-term outcomes are satisfactory for patients with complicated primary hepatolithiasis.This minimally invasive method has the potential to be used as a substitute for hepatobiliary surgery.

Three-dimensional stability calculation method for high and large composite slopes formed by mining stope and inner dump in adjacent open pits

The 2D limit equilibrium method is widely used for slope stability analysis.However,with the advancement of dump engineering,composite slopes often exhibit significant 3D mechanical effects.Consequently,it is of significant importance to develop an effective 3D stability calculation method for composite slopes to enhance the design and stability control of open-pit slope engineering.Using the composite slope formed by the mining stope and inner dump in Baiyinhua No.1 and No.2 open-pit coal mine as a case study,this research investigates the failure mode of composite slopes and establishes spatial shape equations for the sliding mass.By integrating the shear resistance and sliding force of each row of microstrip columns onto the bottom surface of the strip corresponding to the main sliding surface,a novel 2D equivalent physical and mechanical parameters analysis method for the strips on the main sliding surface of 3D sliding masses is proposed.Subsequently,a comprehensive 3D stability calculation method for composite slopes is developed,and the quantitative relationship between the coordinated development distance and its 3D stability coefficients is examined.The analysis reveals that the failure mode of the composite slope is characterized by cutting-bedding sliding,with the arc serving as the side interface and the weak layer as the bottom interface,while the destabilization mechanism primarily involves shear failure.The spatial form equation of the sliding mass comprises an ellipsoid and weak plane equation.The analysis revealed that when the coordinated development distance is 1500 m,the error rate between the 3D stability calculation result and the 2D stability calculation result of the composite slope is less than 8%,thereby verifying the proposed analytical method of equivalent physical and mechanical parameters and the 3D stability calculation method for composite slopes.Furthermore,the3D stability coefficient of the composite slope exhibits an exponential correlation with the coordinated development distance,with the coefficient gradually decreasing as the coordinated development distance increases.These findings provide a theoretical guideline for designing similar slope shape parameters and conducting stability analysis.

Study of three-dimensional spatial diffuse discharge in contact electrode structure applied to air purification

In this work,based on the role of pre-ionization of the non-uniform electric field and its effect of reducing the collisional ionization coefficient,a diffuse dielectric barrier discharge plasma is formed in the open space outside the electrode structure at a lower voltage by constructing a three-dimensional non-uniform spatial electric field using a contact electrode structure.The air purification study is also carried out.Firstly,a contact electrode structure is constructed using a three-dimensional wire electrode.The distribution characteristics of the spatial electric field formed by this electrode structure are analyzed,and the effects of the non-uniform electric field and the different angles of the vertical wire on the generation of three-dimensional spatial diffuse discharge are investigated.Secondly,the copper foam contact electrode structure is constructed using copper foam material,and the effects of different mesh sizes on the electric field distribution are analyzed.The results show that as the mesh size of the copper foam becomes larger,a strong electric field region exists not only on the surface of the insulating layer,but also on the surface of the vertical wires inside the copper foam,i.e.,the strong electric field region shows a three-dimensional distribution.Besides,as the mesh size increases,the area of the vertical strong electric field also increases.However,the electric field strength on the surface of the insulating layer gradually decreases.Therefore,the appropriate mesh size can effectively increase the discharge area,which is conducive to improving the air purification efficiency.Finally,a highly permeable stacked electrode structure of multilayer wire-copper foam is designed.In combination with an ozone treatment catalyst,an air purification device is fabricated,and the air purification experiment is carried out.

Prediction of three-dimensional ocean temperature in the South China Sea based on time series gridded data and a dynamic spatiotemporal graph neural network

Ocean temperature is an important physical variable in marine ecosystems,and ocean temperature prediction is an important research objective in ocean-related fields.Currently,one of the commonly used methods for ocean temperature prediction is based on data-driven,but research on this method is mostly limited to the sea surface,with few studies on the prediction of internal ocean temperature.Existing graph neural network-based methods usually use predefined graphs or learned static graphs,which cannot capture the dynamic associations among data.In this study,we propose a novel dynamic spatiotemporal graph neural network(DSTGN)to predict threedimensional ocean temperature(3D-OT),which combines static graph learning and dynamic graph learning to automatically mine two unknown dependencies between sequences based on the original 3D-OT data without prior knowledge.Temporal and spatial dependencies in the time series were then captured using temporal and graph convolutions.We also integrated dynamic graph learning,static graph learning,graph convolution,and temporal convolution into an end-to-end framework for 3D-OT prediction using time-series grid data.In this study,we conducted prediction experiments using high-resolution 3D-OT from the Copernicus global ocean physical reanalysis,with data covering the vertical variation of temperature from the sea surface to 1000 m below the sea surface.We compared five mainstream models that are commonly used for ocean temperature prediction,and the results showed that the method achieved the best prediction results at all prediction scales.

Computed tomography three-dimensional reconstruction in the diagnosis of bleeding small intestinal polyps:A case report

BACKGROUND Computed tomography(CT)small bowel three-dimensional(3D)reconstruction is a powerful tool for the diagnosis of small bowel disease and can clearly show the intestinal lumen and wall as well as the outside structure of the wall.The horizontal axis position can show the best adjacent intestinal tube and the lesion between the intestinal tubes,while the coronal position can show the overall view of the small bowel.The ileal end of the localization of the display of excellent,and easy to quantitative measurement of the affected intestinal segments,the sagittal position for the rectum and the pre-sacral lesions show the best,for the discovery of fistulae is also helpful.Sagittal view can show rectal and presacral lesions and is useful for fistula detection.It is suitable for the assessment of inflammatory bowel disease,such as assessment of disease severity and diagnosis and differential diagnosis of the small bowel and mesenteric space-occupying lesions as well as the judgment of small bowel obstruction points.CASE SUMMARY Bleeding caused by small intestinal polyps is often difficult to diagnose in clinical practice.This study reports a 29-year-old male patient who was admitted to the hospital with black stool and abdominal pain for 3 months.Using the combination of CT-3D reconstruction and capsule endoscopy,the condition was diagnosed correctly,and the polyps were removed using single-balloon enteroscopyendoscopic retrograde cholangiopancreatography without postoperative complications.CONCLUSION The role of CT-3D in gastrointestinal diseases was confirmed.CT-3D can assist in the diagnosis and treatment of gastrointestinal diseases in combination with capsule endoscopy and small intestinal microscopy.

Research on Quantitative Identification of Three-Dimensional Connectivity of Fractured-Vuggy Reservoirs

The fractured-vuggy carbonate oil resources in the western basin of China are extremely rich.The connectivity of carbonate reservoirs is complex,and there is still a lack of clear understanding of the development and topological structure of the pore space in fractured-vuggy reservoirs.Thus,effective prediction of fractured-vuggy reservoirs is difficult.In view of this,this work employs adaptive point cloud technology to reproduce the shape and capture the characteristics of a fractured-vuggy reservoir.To identify the complex connectivity among pores,fractures,and vugs,a simplified one-dimensional connectivity model is established by using the meshless connection element method(CEM).Considering that different types of connection units have different flow characteristics,a sequential coupling calculation method that can efficiently calculate reservoir pressure and saturation is developed.By automatic history matching,the dynamic production data is fitted in real-time,and the characteristic parameters of the connection unit are inverted.Simulation results show that the three-dimensional connectivity model of the fractured-vuggy reservoir built in this work is as close as 90%of the fine grid model,while the dynamic simulation efficiency is much higher with good accuracy.