Force Equilibrium Research Articles

An Ultra‐Smooth rGO Nano‐Thin Film from a Homogeneous Thin Liquid Film Confined by a Conical Fiber Array: Toward the Highly Sensitive Pressure Sensor

AbstractReduced graphene oxide (rGO) thin films have demonstrated various advantages in flexible electronics. So far, solution processes are widely used for making rGO thin films for their mild operation conditions and low cost. However, wrinkles are frequently formed due to the capillary contraction during drying, which severely deteriorates the conductivity and the transparency of rGO thin films. Here, an ultra‐smooth rGO nano‐thin film, featuring wrinkle‐free and a rather low surface roughness of 1.38 nm is developed, which is attributable to a steady and homogeneous thin liquid film confined by a conical fiber array. The thin liquid film facilitates both in‐plane stacking of nanosheets and reducing the buried solvent under nanosheets. Notably, with the capillary flow replenishing the tri‐phase contact line evaporation, a semi‐dry film near the tri‐phase contact line is produced, that enables the force equilibrium of multi‐directional capillary forces on dispersed nanosheets for depositing a wrinkle‐free nanofilm. The as‐prepared rGO nanofilm gives a low sheet resistance of 8.3 kΩ sq−1 and a high transmittance of 91.9%, showing better performances than other reported rGO nanofilms. On this basis, it is demonstrated that a high‐sensitive pressure sensor with a detection limit as low as 0.02 Pa, highlighting the solution‐processed innovative ultra‐smooth 2D nanomaterial thin films, and devices.

Dynamics of weighted flexible ribbons in a uniform flow

This study explores the dynamics of flexible ribbons with an added weight $G$ at the tail in uniform flow, considering key parameters like inflow Reynolds number ( $Re_u$ ), mass ratio ( $M_t$ ) and aspect ratio ( ${A{\kern-4pt}R}$ ). For two-dimensional ribbons, a simplified theoretical model accurately predicts equilibrium configurations and forces. Inspired by Barois & De Langre (J. Fluid Mech., vol. 735, 2013, R2), we introduce an important control parameter ( $C_G$ ) that effectively collapses normalized forces and angle data. Vortex-induced vibration is observed, and Strouhal number ( $St$ ) scaling laws with $C_G$ are identified. In three-dimensional scenarios, the model effectively predicts lift, but its accuracy in predicting drag is limited to situations with small $Re_u$ values. The flow along the side edges mitigates pressure differences, thereby suppressing vibration and uplift, particularly noticeable in the case of narrow ribbons. This study offers new insights into the dynamics of flexible bodies in uniform flow.

Theoretical model for predicting bubble formation from submerged orifices

Bubble diameter formed at submerged orifices is crucial and challenging to predict with a reliable universal model. A theoretical model covering a wide range of pore structures, operating conditions, and liquid properties for predicting bubble diameter based on bubble force equilibrium at the bubble detachment moment was proposed and verified with good prediction accuracy. The bubble-induced liquid velocity was incorporated to improve the prediction accuracy, and its underlying mechanism on bubble formation diameter was revealed. Additionally, the influencing mechanisms of orifice size, gas inlet velocity, and principal fluid properties on bubble formation diameter were elaborated systematically based on the new mathematical model, and the abnormal phenomenon of bubble formation under high temperatures and elevated pressures is clarified. Based on the proposed high-precision mechanistic model, the theoretical foundations for designing and scaling up industrial aerated reactors and scientific boundary conditions for numerical simulation are provided.

Deformation, shape transformations, and stability of elastic rod loops within spherical confinement

Mechanical insight into the packing of slender objects within confinement is essential for understanding how polymers, filaments, or wires organize and rearrange in limited space. Here we combine theoretical modeling, numerical optimization, and experimental studies to reveal spherical packing behavior of thin elastic rod loops of homogeneous or inhomogeneous stiffness. Across varying loop lengths, a rich array of configurations including circle, saddle, figure-eight, and more intricate patterns are identified. A theoretical framework rooted in the local equilibrium of force and moment is proposed for the rod loop deformation, facilitating the determination of internal and contact forces experienced by the rods during deformation. For the confined homogeneous rod loops, their stable and metastable configurations are well described using proposed Euler rotation curves, which offer a concise and effective approach for configuration prediction. Moreover, formulated analysis on the stability and critical force for homogeneous rod loops on great circles of the spherical confinement are performed. For inhomogeneous rod loops with two segments of differing stiffness, the stiffer segment exhibits less deviation from the great circle, while the softer segment undergoes more pronounced deformation. These findings not only enhance our understanding of buckling and post-buckling phenomena but also offer insights into filament patterning within confining environments.

Analytical model for concrete crack width of steel lined reinforced concrete penstock

Cracking in steel lined reinforced concrete penstock (SLRCP) poses substantial challenges to both serviceability and safety, particularly when crack width exceeds specified limits. Despite its critical implications, accurately calculating crack width remains a persistent challenge. To address this issue, a multi-layer steel-concrete ring analytical model is derived through the concept of orthotropic and thick-walled cylinder theories. In this model, the cracked concrete is assumed to be an orthotropic material incorporating both elastic modulus reduction and tensile softening behavior based on its cracking characteristic of radial penetrating cracks uniformly distributed on the circumference. The steel liner and reinforcing steel bar are modeled as continuous steel rings using equivalent stress theory, while thick-walled cylinder theory is applied due to geometric features and bearing characteristics. Crack width is determined through equilibrium of force and compatibility of deformations at steel-concrete interfaces, utilizing these assumptions and the developed calculation model. To validate the proposed model, data from prototype observations and large-scale model tests of SLRCPs were systematically collected. The crack width results obtained from the proposed model exhibit strong agreement with the collected data, confirming its accuracy and reliability. In conclusion, the proposed analytical model provides an effective solution for calculating the crack width of SLRCP and offers valuable insights for durability evaluation.

Numerical simulation and model development of drag coefficient of bubbles in gas-liquid metal two-phase flow

The occurrence of a steam generator tube rupture (SGTR) accident in a lead-bismuth cooled fast reactor results in the formation of steam bubbles in the liquid lead-bismuth eutectic (LBE). This may degrade heat transfer and power transients in the reactor core due to the migration and accumulation of steam bubbles. To investigate the dynamics of steam bubbles flowing in liquid LBE, it is essential to develop an accurate model for the bubble drag coefficient. In this paper, a three-dimensional numerical model is first established to simulate the injection of high-pressure steam bubbles into a high-temperature LBE molten pool. The model is based on the CLSVOF method. By analyzing the trajectory, velocity, and diameter of bubbles, and combining them with the force equilibrium equation for bubbles, the values of the drag coefficient for bubbles are determined. On this basis, the suitability of current empirical drag models for bubble migration in LBE is evaluated. Finally, the optimal drag coefficient model is selected and further improved. Results reveal that the prediction error of the optimized model for the bubble drag coefficient in liquid LBE is within ±15 %.

Research on Load Estimation for Commercial Trucks Based on Air Suspension State Parameters

For commercial trucks, accurate vehicle load information is crucial to preventing overloading and uneven load distribution, thereby improving transportation efficiency and reducing transportation costs. Existing research on vehicle load estimation primarily focuses on tires, axles, frames, and leaf springs, with limited studies on load estimation for vehicles equipped with air suspension. This paper presents the design of a vehicle onboard weighing system for civil commercial trucks based on the state parameters of air springs and proposes a matching method for vehicle load estimation. The system utilizes laser distance- and gas pressure sensors to collect the state parameters of the air springs. It estimates the vehicle load and displays the load information in real time on the onboard display. The vehicle load estimation method was explored through prototype testing and theoretical derivation. A fitting model was constructed using surface fitting techniques to establish the relationship between the effective contact area and the height variation and capsule pressure of the air spring. Subsequently, a load prediction model was built based on the force equilibrium equation of the air springs. Finally, the construction and accuracy of the load prediction model for a specific type of diaphragm air spring were verified using finite element analysis. The results demonstrated that the average relative error of the load prediction model could be controlled within 1%, meeting the high accuracy requirements of the vehicle weighing system.

Effects of gravity and anisotropy on the convective instability in a nanofluid porous medium

Controlling thermal loads in heat transfer systems that produce high heat flux is a real mechanical challenge. Various cooling/heating techniques are available to achieve the reliability requirements of electronic devices. Among them, the use of nanofluids is one of the important techniques that provide long term stability, higher cooling rates, decreased pumping power needs, and cost-effectiveness. Instability in fluids typically occurs when a disturbance affects the equilibrium of forces of external, inertial, or viscous origin acting on a system. The disturbances may arise from the surroundings due to small changes in the boundary conditions, surface roughness, free stream turbulence, noise, etc. In several circumstances, these disturbances are very small and hence cannot be measured with the available instruments but can be observed after the instability sets in. In this study, we focus on examining the effect of time-periodic vertical gravity modulation on the initiation of natural convection in a porous medium saturated with a nanofluid. The porous medium is characterized by being heated at the bottom and possessing anisotropic properties in terms of permeability and thermal diffusivity. Specifically, we investigate the onset of convection in water-based nanofluids that contain metallic (Cu) and metal oxide ( A l 2 O 3 ) particles. The Khanafer-Vafai-Lightstone (KVL) nanofluid model is employed to describe the behavior of the nanofluids, and the thresholds for convective instability are determined using the Mathieu functions. In addition, we study the instabilities of synchronous and subharmonic modes and explore their transitions along with practical applications.

Eccentric compression behavior of partially CFRP wrapped seawater sea-sand concrete columns reinforced with epoxy-coated rebars

The combined use of seawater sea-sand concrete (SSC), epoxy-coated rebar (ECR) and fiber-reinforced polymer (FRP) can achieve significant economic and environmental benefits in coastal infrastructure constructions. Thus, this paper first presents an experimental and analytical study on the eccentric compression behavior of partially carbon FRP (CFRP) wrapped SSC columns reinforced with ECR. The effects of load eccentricities, thickness of CFRP strips, clear spacing ratios, and slenderness ratios were comprehensively analyzed. Test results show that as the eccentricity increased from 60 mm to 120 mm, the increment of load capacity of partially CFRP wrapped columns reduced by 27.6 % to 17.1 %. Moreover, a larger lateral deflection was detected for the specimens with a relatively high slenderness ratio, which results in an obvious second-order effect and further aggravates the nonlinearity of axial load-bending curves, especially for the specimens under a small eccentricity. In addition, an eccentricity-dependent stress-strain model for concrete was developed under FRP strip-steel hoop composite confinement, and the lateral deflection formula of eccentrically loaded columns was also derived based on force equilibrium analysis. Finally, a new theoretical formula incorporating the effect of the slenderness ratio was also presented for the determination of load and moment capacity of the partially CFRP wrapped columns based on section analysis. By comparing with the test results, the developed model exhibits satisfactory accuracy.

Axial flow impeller design and simulation for negative pressure adsorption climbing robot

Abstract Axial flow impeller is broadly applied in negative pressure adsorption climbing robots, but its design method and procedure are rarely reported. This paper used a wind turbine blade climbing NDT robot as case study, presented a comprehensive design process of the axial flow impeller, and evaluated the adsorption performance of the impeller by carrying out simulation experiments. Firstly, analysing the force equilibrium of the robot and determining minimum adsorption force required by the robot. Secondly, defining design parameters of the axial flow impeller by using fluid dynamics analysis method. Thirdly, specifying structural parameters of the impeller blades using impeller design theory and developing the CAD model of the impeller. Finally, evaluating the adsorption performance of the impeller by carrying out a series of Fluent simulation experiments. Research results showed that the total adsorption force of the impeller reached 80.30N as rotational speed is 6000rpm, which satisfied design requirements of the climbing robot. This research provides a methodological framework in designing axial flow impellers for the purpose of negative pressure adsorption.

P‐167: A simulation model for charged particle motion in INK system

This paper provides a simulation model for the electronic paper display device, the model can quickly and efficiently simulate the initial steady state of the electronic ink system, thus improving the efficiency of device development. It mainly includes the following contents: Obtaining the parameters of the model, then the built‐in electric field strength is calculated through the model parameter information and the pre‐ constructed built‐in electric field model. The motion of each charged particle is simulated to obtain the simulation relationship between the display gray scale and the applied electric field force. Finally, it make various charged particles reach the force equilibrium state, so as to realize the simulation of the voltage‐gray scale relationship of electronic ink system.

Titin mechanical unloading disrupts sarcomere tensional homeostasis triggering fast, non-canonical myocardial fibrosis

Abstract Funding Acknowledgements Type of funding sources: Public grant(s) – EU funding. Main funding source(s): We acknowledge funding from Ministerio de Ciencia e Innovación (Madrid, Spain) through grant FJC2021 047055-I and European Research Council (ERC) through grant H2020-ERC-2020-COG (ProtMechanics-Live). Background/Introduction Stresses exerted during myocardial contraction and relaxation cycles are mainly braced by mechanical proteins located in the extracellular matrix (ECM) and cardiomyocytes (CMs). Titin (TTN) is a gigantic protein that scaffolds the sarcomere (i.e. the basic unit of contraction). Since most of the CM cytoskeleton is occupied by sarcomeres, titin is an essential structural hub for the maintenance of a homeostatic level of mechanical tension. Purpose It has been shown that mechanical disruption of the ECM affects the equilibrium forces between CMs and the ECM, leading to morphological changes in sarcomeres. However, potential alteration of tensional homeostasis resulting from defective mechanics of CM proteins remains unknown due to the lack of tools for in vivo studies. Therefore, the aim of this study is to explore the in vivo response of the myocardium to the abrupt cessation of the mechanical properties of titin. Methods For this purpose, we have experimentally challenged CMs to a tensional unloading by means of a titin cleavage model (TEVs-TTN) [1], both in vitro, using neonatal CM (NMCM), and in vivo. In this model, a cassette with a tobacco etch virus protease (TEVp) recognition sequence has been included in the I-band of titin. Transducing TEVp with adeno-associated virus (myoAAV), we were able to cleave titin in living CMs, ceasing only its mechanical properties. Samples were analyzed using an array of techniques such as histology, immunofluorescence and bulk and single-nucleus RNAseq. Results Our results show that while in NMCM we obtain a complete titin cleavage, in vivo we observe mosaic expression of TEVp and no more than 30% cleavage of the protein. In both cases, TEVp expression in TEVs-TTN does not affect cell viability. However, mild myocardial titin cleavage leads to a fibrotic response in less than six days, characterized by an increase in interstitial collagen and an expansion of cardiac fibroblasts population. Transcriptional and phospo-SMAD2/3 analysis results suggest that this fibrotic response occurs without activation of the canonical TGFβ pathway. Conclusion Taken together, our results suggest that the loss of titin mechanical-structural function in cardiomyocytes derived from the cleavage of a single peptide bond elicits a fast fibrotic compensatory response in the myocardium through intercellular communication with fibroblasts.

A Generalized Limit Equilibrium-Based Platform Incorporating Simplified Bishop, Janbu and Morgenstern–Price Methods for Soil Slope Stability Problems

Limit equilibrium (LE) method is the most widely used method for slope stability analysis. Different methods based on the LE technique for the analysis of the stability of the slope have been developed. Some are based on satisfying the force equilibrium condition of the failing mass (Janbu’s method), while some focus on satisfying the moment equilibrium condition (Bishop’s method). Among these methods, the most accurate result is provided by the Morgenstern–Price method as it not only satisfies both moments as well as a force equilibrium condition but also considers the interslice shear forces (Vi) and interslice normal forces (Ei), which are neglected by most of the LE methods to avoid the condition of indeterminacy. To accommodate these forces, Morgenstern–Price (MP) gave a relation between the Vi and Ei which depends upon a scaling multiplier (λ). Thus, it becomes necessary to evaluate λ value along with the factor of safety (FS). There is barely any work discussing the detailed methodology of evaluation of λ along with FS. Method for obtaining λ along with FS have been developed and elaborated in details here. While calculating FS (MP method), evaluation of Ei is a must which is dependent upon the values of normal force at the base of each slice (Ni) and FS, which itself is dependent upon the value of Ei, making it a loop of interdependent variables. To avoid this interdependency of above stated variables, a separate formulation of Ei is given which reduces the calculations (run-time) involved. A VBA code-based platform has also been developed incorporating the generalized LE method, including Bishop’s, Janbu’s, and Morgenstern–Price methods which are represented in the form of flowcharts in this work.

Enhanced entanglement in multi-bath spin-boson models

The spin-boson model usually considers a spin coupled to a single bosonic bath. However, some physical situations require coupling of the spin to multiple environments. For example, spins interacting with phonons in three-dimensional magnetic materials. Here, we consider a spin coupled isotropically to three independent baths. We show that coupling to multiple baths can significantly increase entanglement between the spin and its environment at zero temperature. The effect of this is to reduce the spin's expectation values in the mean force equilibrium state. In contrast, the classical three-bath spin equilibrium state turns out to be entirely independent of the environmental coupling. These results reveal purely quantum effects that can arise from multi-bath couplings, with potential applications in a wide range of settings, such as magnetic materials.

Extraction of uniaxial stress–strain curve from bending test using DIC measurements

This paper presents a novel methodology for the analytical determination of the uniaxial tensile/compressive stress–strain curve derived from the four-point bending test without a pre-assumed constitutive model. While the concept of extracting the stress–strain curve from a bending test is not new, the proposed methodology differs from existing methods by recognising the force equilibrium in a deformed configuration, resulting in a temporal variation of both a bending moment and an axial force. In addition, the deformation of the beam cross-section is taken into account according to the principles of finite strain theory either via Poisson’s ratio or by assuming volume conservation. The crucial data for this methodology are obtained from an experimentally conducted bending test in which the loading force, displacement and strain fields at the top, bottom and front faces of the specimen are measured using a stereo DIC. While the method is completely analytical, it reconstructs the stress–strain curve in discrete form by solving a system of linear equations. While the existing methods for the four-point bending test allow extraction of the stress–strain curve from a single experiment only for strains up to a few per cent, our methodology allows the extraction of the stress–strain curve in tension and compression up to five times the strains obtained with methods that neglect the above features.

Effect of patient-specific scapular morphology on the glenohumeral joint force and shoulder muscle force equilibrium: a study of rotator cuff tear and osteoarthritis patients.

Introduction: Osteoarthritis (OA) and rotator cuff tear (RCT) pathologies have distinct scapular morphologies that impact disease progression. Previous studies examined the correlation between scapular morphology and glenohumeral joint biomechanics through critical shoulder angle (CSA) variations. In abduction, higher CSAs, common in RCT patients, increase vertical shear force and rotator cuff activation, while lower CSAs, common in OA patients, are associated with higher compressive force. However, the impact of the complete patient-specific scapular morphology remains unexplored due to challenges in establishing personalized models. Methods: CT data of 48 OA patients and 55 RCT patients were collected. An automated pipeline customized the AnyBody™ model with patient-specific scapular morphology and glenohumeral joint geometry. Biomechanical simulations calculated glenohumeral joint forces and instability ratios (shear-to-compressive forces). Moment arms and torques of rotator cuff and deltoid muscles were analyzed for each patient-specific geometry. Results and discussion: This study confirms the increased instability ratio on the glenohumeral joint in RCT patients during abduction (mean maximum is 32.80% higher than that in OA), while OA patients exhibit a higher vertical instability ratio in flexion (mean maximum is 24.53% higher than that in RCT) due to the increased inferior vertical shear force. This study further shows lower total joint force in OA patients than that in RCT patients (mean maximum total force for the RCT group is 11.86% greater than that for the OA group), attributed to mechanically advantageous muscle moment arms. The findings highlight the significant impact of the glenohumeral joint center positioning on muscle moment arms and the total force generated. We propose that the RCT pathomechanism is related to force magnitude, while the OA pathomechanism is associated with the shear-to-compressive loading ratio. Overall, this research contributes to the understanding of the impact of the complete 3D scapular morphology of the individual on shoulder biomechanics.

Development of Improved Confined Compression Testing Setups for Use in Stress Relaxation Testing of Viscoelastic Biomaterials.

The development of cell-based biomaterial alternatives holds significant promise in tissue engineering applications, but it requires accurate mechanical assessment. Herein, we present the development of a novel 3D-printed confined compression apparatus, fabricated using clear resin, designed to cater to the unique demands of biomaterial developers. Our objective was to enhance the precision of force measurements and improve sample visibility during compression testing. We compared the performance of our innovative 3D-printed confined compression setup to a conventional setup by performing stress relaxation testing on hydrogels with variable degrees of crosslinking. We assessed equilibrium force, aggregate modulus, and peak force. This study demonstrates that our revised setup can capture a larger range of force values while simultaneously improving accuracy. We were able to detect significant differences in force and aggregate modulus measurements of hydrogels with variable degrees of crosslinking using our revised setup, whereas these were indistinguishable with the convectional apparatus. Further, by incorporating a clear resin in the fabrication of the compression chamber, we improved sample visibility, thus enabling real-time monitoring and informed assessment of biomaterial behavior under compressive testing.

Energy‐consumption model for rotary‐wing drones

AbstractWith technological advancement, the use of drones in delivery systems has become increasingly feasible. Many companies have developed rotary‐wing drone (RWD) technologies for parcel delivery. At present, the limited endurance is the main disadvantage of RWD delivery. The energy consumption of RWDs must be carefully managed, and it is necessary to develop an effective energy‐consumption model to support RWD flight planning. Because the interaction between the forces on the RWD and its flying environment is very complex, it is challenging to estimate accurately the RWD energy consumption. This study summarizes several energy‐consumption models proposed in the literature, then we develop an RWD energy‐consumption model (called the integrated model) based on analyzing the dynamic equilibrium of forces and power consumption in flight phases (including climb, descent, hover, and horizontal flight). Computational experiments involving several commercial RWDs indicate that the integrated model is more effective than several models in the literature. In the case where an RWD completed one flight segment, on average, 87.63% of the battery capacity was consumed in the horizontal flight phase. We also analyzed the effects of the total mass and horizontal airspeed on the RWD endurance and found that a larger mass corresponded to shorter endurance, and in the experimental range of the horizontal airspeed, a higher horizontal airspeed corresponded to longer endurance. Moreover, the total mass affected the RWD endurance more significantly than the horizontal airspeed.

Disturbance range characterization and reinforcement design method of concealed bedding slope

The concealed bedding slopes exhibit sudden and extensive failure characteristics, harboring significant safety hazards. It is crucial to employ rapid and accurate methods to determine excavation disturbance ranges and reinforcement forces to prevent such slope disasters. This paper first analyzes the failure characteristics of slope excavation to reveal the mechanism of slope failure. Secondly, a novel method was proposed to calculate excavation disturbance ranges and solve for active reinforcement forces. Thirdly, the rationality of the method was validated through numerical simulation experiments. Finally, the method was applied to engineering. The results are as follows:(1) Excavation induces plastic shear failure of structural planes first, and the strength of the structural plane will be reduced, resulting in extrusion shear deformation of the rock mass at the foot of the slope under greater unidirectional load along the structural plane, forming a shear sliding surface. (2) From both strength and deformation perspectives, a method for calculating disturbance ranges in excavated slopes is proposed. This method determines key indicators for analyzing two disturbance range factors: ultimate height and disturbance thickness. The disturbance range is divided into shear plastic zones and traction sliding zones. (3) An active reinforcement design method is proposed considering the timing, scope, and force of reinforcement. The disturbance thickness is used to determine the reinforcement range, while the ultimate height is used to determine the timing of reinforcement. Active reinforcement forces are calculated based on the equilibrium of forces and deformation amounts. (4) A comparison with the Finite Element Method (FEM) shows that during excavation, the disturbance ranges and reinforcement forces calculated using this method are slightly larger than those obtained from the FEM. The calculations are fast, and as the excavation height approaches the limit slope height, the disturbance thickness and horizontal deformations closely match the FEM results, ensuring safety. It can provide monitoring, early warning and reinforcement guidance for bedding slope construction.

Modulating thermal and electrical conductivities in polymers: An approach toward extracting molecular design rules through atomistic simulations

Polymers are extensively employed in diverse industries, including electrical equipment and electronic devices. Recent technological advancements have intensified the demand for dielectric polymers with both high insulation resistance and high thermal conductivity. We employed molecular dynamics simulations to clarify the intricate relationship between molecular structures, thermal conductivity, and ionic mobility from an atomistic point of view. Examined polymers include polyethylene, polyvinyl alcohol, polyvinyl chloride, polyvinylidene fluoride, polytetrafluoroethylene, polychlorotrifluoroethylene, polyoxymethylene, and polyethylene oxide. Based on the elucidated correlations among force field parameters, we found that the parameters can be clustered into four groups: group 1 (atomic and bond parameters and force constant of the bond angle), group 2 (equilibrium angle and force constant of dihedral angle), and group 3 (side chain atom charges). Thermal conductivity showed relationships with parameters in group 1, with correlation coefficients mostly exceeding 0.7 in the absolute value. Considering the correlations between the parameters, we systematically altered the parameters within each group and computed thermal conductivity and ionic mobility. When altering the force field parameters of groups 1 and group 2, a trade-off relationship between thermal conductivity and ionic mobility becomes evident. Conversely, altering the force field parameters in group 3 increased thermal conductivity while decreasing ionic mobility, breaking the trade-off relationship. The proposed clustered-parameter variation method can predict the changes in the electrical and thermal conductivity of polymers through molecular structure modifications. The method, being a general and first-principles approach, is likely to have significant advantages in the molecular design across a diverse range of polymers.