In this work,a method is put forward to obtain the dynamic solution efficiently and accurately for a large-scale train-track-substructure(TTS)system.It is called implicit-explicit integration and multi-time-step solut...In this work,a method is put forward to obtain the dynamic solution efficiently and accurately for a large-scale train-track-substructure(TTS)system.It is called implicit-explicit integration and multi-time-step solution method(abbreviated as mI-nE-MTS method).The TTS system is divided into train-track subsystem and substruc-ture subsystem.Considering that the root cause of low effi-ciency of obtaining TTS solution lies in solving the alge-braic equation of the substructures,the high-efficient Zhai method,an explicit integration scheme,can be introduced to avoid matrix inversion process.The train-track system is solved by implicitly Park method.Moreover,it is known that the requirement of time step size differs for different sub-systems,integration methods and structural frequency response characteristics.A multi-time-step solution is pro-posed,in which time step size for the train-track subsystem and the substructure subsystem can be arbitrarily chosen once satisfying stability and precision demand,namely the time spent for m implicit integral steps is equal to n explicit integral steps,i.e.,mI=nE as mentioned above.The numeri-cal examples show the accuracy,efficiency,and engineering practicality of the proposed method.展开更多
The compact implicit integration factor (cIIF) method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF me...The compact implicit integration factor (cIIF) method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF method to some complex-valued nonlinear evolutionary equations such as the nonlinear SchrSdinger (NLS) equation and the complex Ginzburg-Landau (GL) equation. Detailed algorithm formulation and practical implementation of cIIF method are performed. The numerical results indicate that this method is very accurate and efficient.展开更多
A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first ...A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first sub-step, but evaluated and doubly used in the second sub-step. The order of two sub-steps is reversed for each transverse magnetic field component so that the cross-coupling terms are always expressed in implicit form, thus the calculation is very efficient and stable. Moreover, an improved six-point finite-difference scheme with high accuracy independent of specific structures of waveguide is also constructed to approximate the cross-coupling terms along the transverse directions. The imaginary-distance procedure is used to assess the validity and utility of the present method. The field patterns and the normalized propagation constants of the fundamental mode for a buried rectangular waveguide and a rib waveguide are presented. Solutions are in excellent agreement with the benchmark results from the modal transverse resonance method.展开更多
Integrable systems play a crucial role in physics and mathematics.In particular,the traditional(1+1)-dimensional and(2+1)-dimensional integrable systems have received significant attention due to the rarity of integra...Integrable systems play a crucial role in physics and mathematics.In particular,the traditional(1+1)-dimensional and(2+1)-dimensional integrable systems have received significant attention due to the rarity of integrable systems in higher dimensions.Recent studies have shown that abundant higher-dimensional integrable systems can be constructed from(1+1)-dimensional integrable systems by using a deformation algorithm.Here we establish a new(2+1)-dimensional Chen-Lee-Liu(C-L-L)equation using the deformation algorithm from the(1+1)-dimensional C-L-L equation.The new system is integrable with its Lax pair obtained by applying the deformation algorithm to that of the(1+1)-dimension.It is challenging to obtain the exact solutions for the new integrable system because the new system combines both the original C-L-L equation and its reciprocal transformation.The traveling wave solutions are derived in implicit function expression,and some asymmetry peakon solutions are found.展开更多
高应变率载荷作用下金属材料的变形集中于很窄的区域内,即剪切变形局部化。局部化变形带内产生严重的塑性变形,削弱材料的承载能力,甚至导致材料断裂破坏。基于有限元分析软件FEAP(Finite Element Analysis Program),采用混合有限元方法...高应变率载荷作用下金属材料的变形集中于很窄的区域内,即剪切变形局部化。局部化变形带内产生严重的塑性变形,削弱材料的承载能力,甚至导致材料断裂破坏。基于有限元分析软件FEAP(Finite Element Analysis Program),采用混合有限元方法,用Fortran语言编译适用于金属材料在高应变率下的剪切局部化问题的新单元;计算过程中采用与应变、应变率及温度相关的塑性本构关系来描述剪切带现象,同时在能量平衡方程中考虑剪切带形成过程中的热传导作用;同时考虑显式算法与隐式算法的时间离散方法,并将两种算法的结果进行对比。结果表明,虽然剪切带形成过程很短,一般为微秒量级,但剪切带形成过程中的热扩散项与塑性变形产生的热能量级相同,有效地缓解剪切带模拟的网格敏感性;对于金属材料热塑性剪切带问题,为了满足计算精度要求,显式算法需要的时间步太小,计算成本比隐式迭代高很多;而基于该单元采用隐式算法模拟热塑性剪切带问题迭代收敛稳定,计算精度高,且因为考虑了热传导作用,网格敏感性小。展开更多
基金This work was supported by the National Natural Science Foundation of China(Grant Nos.52008404,U1934217 and 11790283)Science and Technology Research and Development Program Project of China Railway Group Limited(Major Special Project,No.2020-Special-02)the National Natural Science Foundation of Hunan Province(Grant No.2021JJ30850).
文摘In this work,a method is put forward to obtain the dynamic solution efficiently and accurately for a large-scale train-track-substructure(TTS)system.It is called implicit-explicit integration and multi-time-step solution method(abbreviated as mI-nE-MTS method).The TTS system is divided into train-track subsystem and substruc-ture subsystem.Considering that the root cause of low effi-ciency of obtaining TTS solution lies in solving the alge-braic equation of the substructures,the high-efficient Zhai method,an explicit integration scheme,can be introduced to avoid matrix inversion process.The train-track system is solved by implicitly Park method.Moreover,it is known that the requirement of time step size differs for different sub-systems,integration methods and structural frequency response characteristics.A multi-time-step solution is pro-posed,in which time step size for the train-track subsystem and the substructure subsystem can be arbitrarily chosen once satisfying stability and precision demand,namely the time spent for m implicit integral steps is equal to n explicit integral steps,i.e.,mI=nE as mentioned above.The numeri-cal examples show the accuracy,efficiency,and engineering practicality of the proposed method.
文摘The compact implicit integration factor (cIIF) method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF method to some complex-valued nonlinear evolutionary equations such as the nonlinear SchrSdinger (NLS) equation and the complex Ginzburg-Landau (GL) equation. Detailed algorithm formulation and practical implementation of cIIF method are performed. The numerical results indicate that this method is very accurate and efficient.
文摘A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first sub-step, but evaluated and doubly used in the second sub-step. The order of two sub-steps is reversed for each transverse magnetic field component so that the cross-coupling terms are always expressed in implicit form, thus the calculation is very efficient and stable. Moreover, an improved six-point finite-difference scheme with high accuracy independent of specific structures of waveguide is also constructed to approximate the cross-coupling terms along the transverse directions. The imaginary-distance procedure is used to assess the validity and utility of the present method. The field patterns and the normalized propagation constants of the fundamental mode for a buried rectangular waveguide and a rib waveguide are presented. Solutions are in excellent agreement with the benchmark results from the modal transverse resonance method.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.12275144,12235007,and 11975131)K.C.Wong Magna Fund in Ningbo University。
文摘Integrable systems play a crucial role in physics and mathematics.In particular,the traditional(1+1)-dimensional and(2+1)-dimensional integrable systems have received significant attention due to the rarity of integrable systems in higher dimensions.Recent studies have shown that abundant higher-dimensional integrable systems can be constructed from(1+1)-dimensional integrable systems by using a deformation algorithm.Here we establish a new(2+1)-dimensional Chen-Lee-Liu(C-L-L)equation using the deformation algorithm from the(1+1)-dimensional C-L-L equation.The new system is integrable with its Lax pair obtained by applying the deformation algorithm to that of the(1+1)-dimension.It is challenging to obtain the exact solutions for the new integrable system because the new system combines both the original C-L-L equation and its reciprocal transformation.The traveling wave solutions are derived in implicit function expression,and some asymmetry peakon solutions are found.
文摘高应变率载荷作用下金属材料的变形集中于很窄的区域内,即剪切变形局部化。局部化变形带内产生严重的塑性变形,削弱材料的承载能力,甚至导致材料断裂破坏。基于有限元分析软件FEAP(Finite Element Analysis Program),采用混合有限元方法,用Fortran语言编译适用于金属材料在高应变率下的剪切局部化问题的新单元;计算过程中采用与应变、应变率及温度相关的塑性本构关系来描述剪切带现象,同时在能量平衡方程中考虑剪切带形成过程中的热传导作用;同时考虑显式算法与隐式算法的时间离散方法,并将两种算法的结果进行对比。结果表明,虽然剪切带形成过程很短,一般为微秒量级,但剪切带形成过程中的热扩散项与塑性变形产生的热能量级相同,有效地缓解剪切带模拟的网格敏感性;对于金属材料热塑性剪切带问题,为了满足计算精度要求,显式算法需要的时间步太小,计算成本比隐式迭代高很多;而基于该单元采用隐式算法模拟热塑性剪切带问题迭代收敛稳定,计算精度高,且因为考虑了热传导作用,网格敏感性小。