Influence of manufacturing process-induced geometrical defects on the energy absorption capacity of polymer lattice structures

Modern additive manufacturing processes enable fabricating architected cellular materials of complex shape,which can be used for different purposes.Among them,lattice structures are increasingly used in applications requiring a compromise among lightness and suited mechanical properties,like improved energy absorption capacity and specific stiffness-to-weight and strength-to-weight ratios.A dedicated modeling strategy to assess the energy absorption capacity of lattice structures under uni-axial compression loading is presented in this work.The numerical model is developed in a non-linear framework accounting for the strain rate effect on the mechanical responses of the lattice structure.Four geometries,i.e.,cubic body centered cell,octet cell,rhombic-dodecahedron and truncated cuboctahedron 2+,are investigated.Specifically,the influence of the relative density of the representative volume element of each geometry,the strain-rate dependency of the bulk material and of the presence of the manufacturing process-induced geometrical imperfections on the energy absorption capacity of the lattice structure is investigated.The main outcome of this study points out the importance of correctly integrating geometrical imperfections into the modeling strategy when shock absorption applications are aimed for.

A simulation study of water property changes using geometrical alteration in SPC/E

We present a systematic investigation of the impact of changing the geometry structure of the SPC/E water model by performing a series of molecular dynamic simulations at 1 bar （1 bar = 105 Pa） and 298.15 K. The geometric modification includes altering the H-O-H angle range from 90° to 115° and modifying the O-H length range from 0.90 A to 1.10 A in the SPC/E model. The former is achieved by keeping the dipole moment constant by modifying the O-H length, while in the latter only the O-H length is changed. With the larger bond length and angle, we find that the liquid shows a strong quadrupole interaction and high tetrahedral structure order parameter, resulting in the enhancement of the network structure of the liquid. When the bond length or angle is reduced, the hydrogen bond lifetime and self-diffusion constant decrease due to the weakening of the intermolecular interaction. We find that modifying the water molecular bond length leading to the variation of the intermolecular interaction strength is more intensive than changing the bond angle. Through calculating the average reduced density gradient and thermal fluctuation index, it is found that the scope of vdW interaction with neighbouring water molecules is inversely proportional to the change of the bond length and angle. The effect is mainly due to a significant change of the hydrogen bond network. To study the effect of water models as a solvent whose geometry has been modified, the solutions of ions in different solvent environments are examined by introducing NaCI. During the dissolving process, NaCI ions are ideally dissolved in SPC/E water and bond with natural water more easily than with other solvent models.

Geometrically Nonlinear Flutter Analysis Based on CFD/CSD Methods and Wind Tunnel Experimental Verification

This study presents a high-speed geometrically nonlinear flutter analysis calculation method based on the highprecision computational fluid dynamics/computational structural dynamics methods.In the proposed method,the aerodynamic simulation was conducted based on computational fluid dynamics,and the structural model was established using the nonlinear finite element model and tangential stiffness matrix.First,the equilibrium position was obtained using the nonlinear static aeroelastic iteration.Second,the structural modal under a steady aerodynamic load was extracted.Finally,the generalized displacement time curve was obtained by coupling the unsteady aerodynamics and linearized structure motion equations.Moreover,if the flutter is not at a critical state,the incoming flow dynamic pressure needs to be changed,and the above steps must be repeated until the vibration amplitude are equal.Furthermore,the high-speed geometrically nonlinear flutter of the wing-body assemblymodel with a high-aspect ratio was investigated,and the correctness of the method was verified using high-speed wind tunnel experiments.The results showed that the geometric nonlinearity of the large deformation of the wing caused in-plane bending to become a key factor in flutter characteristics and significantly decreased the dynamic pressure and frequency of the nonlinear flutter compared to those of the linear flutter.

Sparsity-based efficient simulation of cluster targets electromagnetic scattering

An efficient and real-time simulation method is proposed for the dynamic electromagnetic characteristics of cluster targets to meet the requirements of engineering practical applications.First,the coordinate transformation method is used to establish a geometric model of the observation scene,which is described by the azimuth angles and elevation angles of the radar in the target reference frame and the attitude angles of the target in the radar reference frame.Then,an approach for dynamic electromagnetic scattering simulation is proposed.Finally,a fast-computing method based on sparsity in the time domain,space domain,and frequency domain is proposed.The method analyzes the sparsity-based dynamic scattering characteristic of the typical cluster targets.The error between the sparsity-based method and the benchmark is small,proving the effectiveness of the proposed method.

Effect of geometrical parameters on forming quality of high-strength TA18 titanium alloy tube in numerical control bending

The forming quality of high-strength TA18 titanium alloy tube during numerical control bending in changing bending angle β, relative bending radius R/D and tube sizes such as diameter D and wall thickness t was clarified by finite element simulation. The results show that the distribution of wall thickness change ratio Δt and cross section deformation ratio ΔD are very similar under different β; the Δt and ΔD decrease with the increase of R/D, and to obtain the qualified bent tube, the R/D must be greater than 2.0; the wall thinning ratio Δto slightly increases with larger D and t, while the wall thickening ratio Δti and ΔD increase with the larger D and smaller t; the Δto and ΔD firstly decrease and then increase, while the Δti increases, for the same D/t with the increase of D and t.

Three-dimensional spatial structure of the macro-pores and flow simulation in anthracite coal based on X-ray μ-CT scanning data

The three-dimensional(3 D) structures of pores directly affect the CH4 flow.Therefore,it is very important to analyze the3 D spatial structure of pores and to simulate the CH4 flow with the connected pores as the carrier.The result shows that the equivalent radius of pores and throats are 1-16 μm and 1.03-8.9 μm,respectively,and the throat length is 3.28-231.25 μm.The coordination number of pores concentrates around three,and the intersection point between the connectivity function and the X-axis is 3-4 μm,which indicate the macro-pores have good connectivity.During the single-channel flow,the pressure decreases along the direction of CH4 flow,and the flow velocity of CH4 decreases from the pore center to the wall.Under the dual-channel and the multi-channel flows,the pressure also decreases along the CH4 flow direction,while the velocity increases.The mean flow pressure gradually decreases with the increase of the distance from the inlet slice.The change of mean flow pressure is relatively stable in the direction horizontal to the bedding plane,while it is relatively large in the direction perpendicular to the bedding plane.The mean flow velocity in the direction horizontal to the bedding plane(Y-axis) is the largest,followed by that in the direction horizontal to the bedding plane(X-axis),and the mean flow velocity in the direction perpendicular to the bedding plane is the smallest.

Discrete element simulation of crushable rockfill materials

A discrete element method was used to study the evolution of particle crushing in a rockfill sample subjected to triaxial shear. A simple procedure was developed to generate clusters with arbitrary shapes, which resembled real rockfill particles. A theoretical method was developed to define the failure criterion for an individual particle subjected to an arbitrary set of contact forces. Then, a series of numerical tests of large-scale drained triaxial tests were conducted to simulate the behaviors of the rockfill sample. Finally, we examined the development of micro-characteristics such as particle crushing, contact characteristics, porosity, deformation, movement, and energy dissipation. The simulation results were partially compared with the laboratory experiments, and good agreement was achieved, demonstrating that the particle crushing model proposed can be used to simulate the drained triaxial test ofrockfill materials. Based on a comparison of macro behaviors of the roekfill sample and micro structures of the particles, the microscopic mechanism of the rockfill materials subjected to triaxial shear was determined qualitatively. It is shown that the crushing rate, rather than the number of crushed particles, can be used to reflect the relationship between macro- and micro-mechanical characteristics of rockfill materials. These research results further develop our understanding of the deformation mechanism of rockfill materials.

A note on the numerical simulations of flow past a wavy square-section cylinder

The flow past a square-section cylinder with a geometric disturbance is investigated by numerical simulations. The extra terms, due to the introduction of mapping transformation simulating the effect of disturbance into the transformed Navier-Stokes equations, are correctly derived, and the incorrect ones in the previous literature are pointed out and analyzed. Furthermore, the relationship between the vorticity, especially on the cylinder surface, and the disturbance is derived and explained theoretically. The computations are performed at two Reynolds numbers of 100 and 180 and three amplitudes of waviness of 0.006, 0.025 and 0.167 with another aim to explore the effects of different Reynolds numbers and disturbance on the vortex dynamics in the wake and forces on the body. Numerical results have shown that, at the mild waviness of 0.025, the Kairmain vortex shedding is suppressed completely for Re = 100, while the forced vortex dislocation is appeared in the near wake at the Reynolds number of 180. The drag reduction is up to 21.6% at Re = 100 and 25.7% at Re = 180 for the high waviness of 0.167 compared with the non-wavy cylinder. The lift and the Strouhal number varied with different Reynolds numbers and the wave steepness are also obtained.

Impacts of Disk Rock Sample Geometric Dimensions on Shear Fracture Behavior in a Punch Shear Test

Punch shear tests have been widely used to determine rock shear mechanical properties but without a standard sample geometric dimension suggestion.To investigate the impacts of sample geometric dimensions on shear behaviors in a punch shear test,simulations using Particle Flow Code were carried out.The effects of three geometric dimensions(i.e.,disk diameter,ratio of shear surface diameter to disk diameter,and ratio of disk height to shear surface diameter)were discussed.Variations of shear strength,shear stiffness,and shear dilatancy angles were studied,and the fracture processes and patterns of samples were investigated.Then,normal stress on the shear surface during test was analyzed and a suggested disk geometric dimension was given.Simulation results show that when the ratio of the shear surface diameter to the disk diameter and the ratio of disk height to the shear surface diameter is small enough,the shear strength,shear stiffness,and shear dilatancy angles are extremely sensitive to the three geometric parameters.If the ratio of surface diameter to disk diameter is too large or the ratio of disk height to surface diameter is too small,a part of the sample within the shear surface will fail due to macro tensile cracks,which is characterized by break off.Samples with a greater ratio of disk height to shear surface diameter,namely when the sample is relatively thick,crack from one end to the other while others crack from both ends towards the middle.During test,the actual normal stress on the shear surface is greater than the target value because of the extra compressive stress from the part of sample outside shear surface.

NUMERICAL SIMULATION AND EXPERIMENTS FOR PANEL FORMING BASED ON CAD/CAE SYSTEM

A comparison research about the implicit and explicit solutions of sheet forming simulation was presented. On the platform of Autoform and Ls dyna3D, a dynamic forming simulation of a sideframe, of Santana 2000 was done, and the engineering strain, the thickness distribution and the FLD between the two softwares were compared. It indicates that their results coincide with each other very much and the areas of the wrinkle and failure are the same roughly. Further, the characteristics of the two softwares in geometric model and preprocessing of the finite element were discussed and the questions which need attention provided.

NEW METHOD FOR GEOMETRIC NONLINEAR ANALYSIS OF LARGE DISPLACEMENT DRILL STRINGS

Based on the actual measured well depth, azimuth and oblique angles, a novel interpolation method to obtain the well axis is developed. The initial stress of drill string at the reference state consistent with well axis can be obtained from the curvature and the tortuosity of well axis. By using the principle of virtual work, the formula to compute the equivalent load vector of the initial stress was derived. In the derivation,the natural （curvilinear） coordinate system was adopted since both the curvature and the tortuosity were generally not zero. A set of displacement functions fully reflecting the rigid body modes was used. Some basic concepts in the finite element analysis of drill string were clarified. It is hoped that the proposed method would offer a theoretical basis for handling the geometric nonlinear problem of the drill string in a 3-D larg edisplacement wellbore.

Quasi-static Deployment Simulation for Deployable Space Truss Structures

A new method was proposed for quasi-static deployment analysis of deployable space truss structures. The structure is assumed a rigid assembly, whose constraints are classified as three categories:rigid member constraint, joint-attached kinematic constraint and boundary constraint. And their geometric constraint equations and derivative matrices are formulated. The basis of the null space and M-P inverse of the geometric constraint matrix are employed to determine the solution for quasi-static deployment analysis. The influence introduced by higher terms of constraints is evaluated subsequently. The numerical tests show that the new method is efficient.

Modeling and Simulation of Particle-Packing Structures and Their Stability Using the Distinct Element Method

A numerical method for simulating the stability of particle-packing structures is presented. The packing structures were modeled on the basis of face-centered cubic (fcc) and body-centered cubic (bcc) structures, and the stability of these structures was investigated using the distinct element method. The interaction between the particles was simplified by considering repulsive, adhesive, and damping forces, and the stability against the gravitational force was simulated. The results under a certain set of parameters showed characteristic deformation when the particles were arranged in an fcc array. Focusing on the local structure, the resulting model was divided into several domains: The bottom base, four top corners, and intermediate domains. The bottom base notably became a body-centered tetragonal (bct) structure, which corresponds to a uniaxially compressed bcc structure. Conversely, the models based on the bcc arrangement were structurally stable, as no specific deformation was observed, and a monotonously compressed bct structure was obtained. Consequently, the bcc arrangement is concluded to be more stable against uniaxial compression, such as the gravitational force, in a particle-packing system.

大气总温传感器数值计算耦合分析技术

借助自研软件开发及现有的商业软件建立阻滞室大气总温传感器模型,用UG软件构建几何模型,以商用Fluent软件模拟传感器的三维CFD数值,采用共形网格划分实现大气总温传感器的流体与固体界面网格匹配。采用Realizable k-ε湍流模型以及数值仿真批处理脚本的方法,对大气总温传感器流、固耦合仿真进行计算,提高了仿真计算效率。利用该方法对常温常压风洞对大气总温传感器的校准进行仿真分析,得出常温常压条件下马赫数大于1时,恢复修正系数基本不随马赫数变化。在常温、马赫数为0.6条件下,在静压条件±0.1atm时或总温条件为-25~45℃时,恢复修正系数变化幅度很小,数值仿真结果具有一定的准确性。

高速S-CO_(2)向心透平几何参数优化及变工况特性分析

透平是热力循环中实现热功转换的核心设备,优化透平几何参数以改善气动性能,是提升系统热效率的关键。该文以面向某船舶主机烟气余热利用的超临界二氧化碳布雷顿循环透平设计参数为基础,给出向心透平主要几何参数,采用耦合CO_(2)真实物性参数的CFX仿真方法,分析几何参数动、静叶包角以及叶轮出口几何角对透平性能的影响规律并进行参数优化,进一步研究非设计工况下的透平运行特性。结果表明:增大静叶包角会使动叶前缘出现扰动涡流并减少余速损失,增大叶轮出口几何角,动叶包角对透平效率的影响规律不变,但效率峰值点会沿动叶包角减小的方向移动;动、静叶包角和叶轮出口几何角分别为49°、21.5°、25°时透平设计点总-静等熵效率高达87.99%,相比优化前,提高0.54个百分点;透平偏离设计点±20%工况下等熵效率大于80%。结果可为超临界二氧化碳高速向心透平设计提供一定参考。

走滑断层几何结构复杂区对地震破裂传播影响的研究回顾

活动走滑断裂带上强震频发,对几何结构复杂区及其导致的分段破裂模式开展研究对于评估区域潜在地震的最大震级具有重要意义。文中从野外地质调查和数值模型研究2个方面入手,对走滑断层上地震破裂行为中的几何结构、断层分段和级联破裂等相关研究进行回顾总结。野外地质调查发现,走滑断裂带上广泛存在的断层几何结构复杂区可能会终止地震破裂的传播。前人基于地表破裂带野外调查,给出了走滑断裂带上几何结构复杂区对地震破裂方式的影响,以及这种影响程度的统计分析结果。另一方面,与断层动态破裂过程相关的数值模拟研究,不仅从动力学角度展示了断层几何结构复杂区对断层动态破裂传播的控制作用,还表明初始应力、岩石介质性质等其他因素也影响了断层动态破裂的传播。在研究程度较高的中国地震科学实验场区,基于断层探测、地球物理场密集观测和高性能计算的3D动力学模拟有望进一步深化对区域复杂断层系统中破裂行为的认识,并为判定区域最大潜在震级提供帮助。

纬编双面移圈织物多层弹簧-质点结构模型构建与实现

纬编双面移圈织物工艺结构复杂多变,现有系统在进行设计开发时难以直观地展示出织物的组织结构特点。为实现纬编双面移圈织物的组织结构模拟,本文提出了一种纬编双面移圈织物多层弹簧-质点结构模型。首先建立了成圈、集圈、移圈线圈几何结构模型,接着分析了纬编双面移圈织物结构特征,并在传统弹簧-质点模型的基础上构建了多层弹簧-质点模型,然后采用牛顿第二定律和Velocity-Verlet数值积分法求解了多层弹簧-质点模型的最终位置,最终提出了一种纬编双面移圈织物结构模拟的算法流程。本文建立的多层弹簧-质点结构模型能够反映纬编双面移圈织物的组织结构,具有较好的结构模拟效果。

锂离子电池用PET-Cu复合集流体拉伸性能研究

随着锂离子电池技术的不断发展,复合集流体因其具有显著提升电池能量密度和安全性高等优势而引起了产业界的广泛关注。而复合集流体的力学性能可靠性是保证其安全服役的前提,为此,本文针对商业用聚对苯二甲酸乙二醇酯-铜(PET-Cu)复合集流体的拉伸性能开展了系统性的研究。利用扫描电镜、激光共聚焦显微镜、X射线衍射等技术手段对复合集流体材料的表面形貌、微观结构及拉伸断裂行为进行了表征与分析。借助有限元模拟分析了不同几何形状试样在拉伸过程中的应力分布状态;通过数字图像相关技术辅助的拉伸实验研究了取样方向、应变速率及试样几何尺寸对PET-Cu复合集流体拉伸性能的影响规律,并对其应变测量进行了校正。结果表明,采用狗骨状试样并进行应变校正可以更准确地评估复合集流体的拉伸性能;PET-Cu复合集流体的拉伸性能表现出试样取样方向相关性和应变速率敏感性。此外,PET-Cu复合集流体的断裂伸长率表现出明显的试样几何尺寸效应,通过引入几何尺寸系数K可以对不同尺寸的试样拉伸性能进行合理的预测。本研究为聚合物复合集流体实际应用提供了理论依据与可靠性的论证,也为相关试验标准的建立提供了参考,有望推动高性能锂离子电池的开发。

纱罗织物三维几何模型的构建

为进一步提升纱罗组织模拟效果,提出一种纱罗织物三维几何模型。将纱组织及罗组织进行了结构分析,对绞组内的经纱及纬纱的空间几何关系进行了数学推导,并将其综合为纱罗组织的数学模型;在此基础上,提出经纱扭绞及纬纱阵列算法,计算经纬纱上每个型值点的三维坐标;最终通过样条曲线生成经纬纱的中心线轨迹,以复丝作为纱线模型,对纱罗组织进行了三维模拟。模拟效果表明:该模拟算法生成的纱罗组织结构准确,表层纱孔清晰,取得了预期效果。

极端工况下几何精确梁理论和CR理论的适用性研究

几何精确梁理论和共旋(Co-Rotational,CR)梁理论是模拟柔性梁动力响应的两种常用方法,然而,针对极端工况下这两种方法的仿真对比分析尚不充分.本文首先基于几何精确梁理论和共旋坐标理论推导了空间大转动梁的刚度矩阵、质量矩阵以及非线性求解的迭代方法,并从理论上讨论了两种方法的特点.最后编写了两种方法的有限元仿真程序,并通过具体算例对比了它们的收敛性和计算效率.仿真结果表明:在足够的单元数量下,两种方法都能提供准确的数值结果;当单元数量较少时,共旋公式的数值结果比几何精确梁理论更精确,并且在相对少量的单元情况下仍能提供足够的精度,这是由于两种方法在迭代过程中对变量更新的不同处理方式所致;在相同单元数量下,几何精确梁理论展现出更高的计算效率和更好的收敛性.本文的研究结果对于柔性梁的数值仿真具有重要的参考价值.