有限条法和广义变分原理在圆柱形水池分析中的应用

本方法是近似方法,适用于等截面、变截面圆柱形水池静力分析,并有很高的精度。 1.池壁的计算 假设:池壁可独立于底板、面板进行单独计算,横向荷载和结构具有轴对称性,不计竖向荷载。因而仅须考虑轴对称的径向位移。有两种计算方法。     

土质边坡极限荷载上限计算方法研究

以边坡极限分析的上限解析式为标准,在无重力情况下利用不同屈服准则的数值方法计算了相同工况下边坡的极限荷载。针对一个实际边坡工况使用解析和FLAC方法计算了边坡的极限荷载上限,并且基于强度折减法使用FLAC确定了各种荷载下边坡的安全系数。通过与解析方法计算结果的比较发现,基于Mohr-Coulomb准则的FLAC软件计算出的边坡极限荷载更加接近解析解。不同边坡顶角或不同内摩擦角的计算结果显示FLAC计算出的边坡极限荷均稍大于解析解,并且内摩擦角越大差值越大,但两者的差别最大不超过5%。     

轴对称圆薄板大挠度微分方程的数值分析方法

根据轴对称问题的特点,利用级数展开和求极限法则,证明了轴对称大挠度圆薄板在圆心处应满足的边界条件,并以圆薄板轴对称大挠度弯曲变形微分方程为基础,建立了圆心处非奇异的轴对称大挠度圆板弯曲微分方程,从而可以方便地利用现有的常微分方程数值求解方法(如变步长龙格-库塔法)对实心圆板的轴对称问题进行数值求解,又不必像摄动法那样推导复杂的公式。在数值求解轴对称圆板大挠度弯曲变形微分方程时,将非线性微分方程的求解主要归结为迭代求解圆心处三个未知边界条件的问题,即圆心处的径向膜力、圆心处的挠度、圆心处挠度的二阶导数,并提出了相应的求解方法。实例中,对于圆薄板受均布横向荷载的问题,分析了周边固支边界条件下的非线性弯曲问题,给出了中心挠度参数大范围变化时的荷载和部分边界值变化曲线,并与经典摄动解进行了对比。对比结果可见,本文方法和摄动法的解非常接近,在量纲归一化中心挠度不超过4.0时,两种方法解的相对误差均小于5.0%。另外,本文还分析了与挠度有关的液体压力作用下和集中荷载作用下周边固支圆板的非线性弯曲问题。通过算例可见:本文方法可以灵活处理不同的荷载问题;对于不同的问题,计算过程相似,不必推导复杂的计算公式,计算精度容易控制。     

用配点法计算受轴对称分布荷载的圆形大挠度板

本文以幂函数多项式作位移函数,用配点法计算了受轴对称分布荷载的圆形大挠度板。文中考虑了固定夹紧、可移夹紧、铰链支承、简单支承等四种边界条件。在计算例题中,荷载采用多项式形式、余弦函数形式的分布荷载以及均布边弯矩荷载或者由它们组合而成的荷载。通过上百个例题的计算,表明此法具有精度高、收敛快等优点。本文还将其所得的结果同用摄动法等得到的结果进行了比较。     

结构性黏土模型在FLAC中的二次开发及应用

FLAC软件广泛应用于岩石及岩土工程数值分析中,但内置材料模型库的不足影响了其应用范围及分析结果的精细化。本文基于FLAC程序提供的二次开发平台,以结构性黏土边界面模型为例,在FISH编译环境下实现了复杂本构模型的开发与应用。首先,介绍了FLAC的关键开发技术和具体实施方法;其次,给出了结构性黏土边界面模型的基本原理及应力积分方法;然后,通过FLAC单元测试与FORTRAN程序计算结果的对比,验证了FLAC二次开发在多种应力路径下的正确性;随后,通过FLAC排水平面应变剪切试验模拟剪切带的形成,测试了二次开发模型的多单元模拟能力,同时验证了其网格依赖性问题;最后,通过对土坡开挖的模拟,测试了结构性黏土模型FLAC二次开发的工程适用性。     

冲击荷载下土中浅埋箱形结构反应的模型试验研究

本文叙述了土中浅埋箱形结构模型动力试验结果。试验的目的是检验我们在文献[8、9]中,所提出的地表冲击波荷载作用下,土中箱形结构动力反应分析的两步法。试验有两种情况:一是利用室内模爆器进行的小模型(模型尺寸与实际结构之比约为1:40)试验;另一种是利用室外试验坑中进行的模拟平面波加载下的中等模型(尺寸缩比约为1:10)试验。试验的箱形结构有单孔与双孔两种。根据埋深与超压的不同,共进行了20炮次试验,试验对自由场中应力、结构周边压力、结构应变、结构运动等参数进行了测量。试验结果表明,所提出的两步分析法中关于土中结构内力分析的准静态法是可行的。     

受预张力或均布边缘径向荷载的圆薄板的大挠度问题

本文以幂函数为试函数,用配点法计算圆薄板的轴对称大挠度,圆板可先受预张力的作用。荷载为轴对称分布荷载、均布边缘力矩及径向力或它们的联合作用。计算结果表明,本法求解的收敛速度和收敛范围远远超过幂级数解法。     

有限条法和广义变分原理在圆柱形水池分析中的应用

本方法是近似方法,适用于等截面、变截面圆柱形水池静力分析,并有很高的精度。 1.池壁的计算 假设:池壁可独立于底板、面板进行单独计算,横向荷载和结构具有轴对称性,不计竖向荷载。因而仅须考虑轴对称的径向位移。有两种计算方法。 1) 有限条法 图1及图2表示池壁中面、正交曲线坐标系及单元、局部坐标系及结线位移。可由低阶平板有限条得单元径向位移     

四种非正交边界板考虑初始荷载效应的后期荷载位移近似解

基于考虑初始荷载效应情况下板的一般形式的静力平衡微分方程,运用坐标变换得到了轴对称情形,考虑初始荷载效应后圆形板的极坐标形式的静力平衡微分方程。运用Galerkin法解得了简支等边三角形板、固支椭圆板、固支圆形板和简支圆形板四种非正交边界板考虑初始荷载效应的后期荷载位移近似解。运用相关文献提出的有限元法验证了近似解的正确性。各位移近似解表达式简单、物理意义明确,清楚地反映了初始荷载及相关因素对后期荷载位移的影响。计算分析表明：初始荷载效应提高了板的弯曲刚度,减小了板的后期荷载位移;板的初始荷载效应主要受初始荷载、跨厚比及边界条件等因素的影响。     

非对称荷载作用下轴对称体应力分析的有限元半解析法

问题的提出及方法概述在工程实际中,有不少情况是属于轴对称体承受非对称荷载作用的问题.土建工程中有很多这样的事例,例如,制动力作用的路面和机场跑道,任意荷载作用的地基,承受非对称荷载的圆筒或圆柱体等.机械另件也有类似的情 ...     

Unified elastoplastic finite difference and its application

Two elastoplastic constitutive models based on the unified strength theory (UST) are established and implemented in an explicit finite difference code, fast Lagrangian analysis of continua (FLAC/FLAC3D), which includes an associated/non-associated flow rule, strain-hardening/softening, and solutions of singularities. Those two constitutive models are appropriate for metallic and strength-different (SD) materials, respectively. Two verification examples are used to compare the computation results and test data using the two-dimensional finite difference code FLAC and the finite element code ANSYS, and the two constitutive models proposed in this paper are verified. Two application examples, the large deformation of a prismatic bar and the strain-softening behavior of soft rock under a complex stress state, are analyzed using the three-dimensional code FLAC3D. The two new elastoplastic constitutive models proposed in this paper can be used in bearing capacity evaluation or stability analysis of structures built of metallic or SD materials. The effect of the intermediate principal stress on metallic or SD material structures under complex stress states, including large deformation, three-dimensional and non-association problems, can be analyzed easily using the two constitutive models proposed in this paper.     

平头弹丸正撞下钢筋混凝土靶板厚度方向的开裂

主要针对钢筋混凝土靶板在受到平头弹丸撞击下发生的厚度方向开裂的问题进行研究,并提出了一个弹丸低速撞击有限厚度板的二阶段模型。模型中第一阶段为侵彻阶段,弹丸受到混凝土介质的侵彻阻力由静阻力和速度效应引起的动阻力组成;模型中第二阶段为开裂阶段,钢筋混凝土靶板发生动态剪切破坏的最大承载力可以通过静态剪切破坏最大承载力乘以一个动态增强因子得到。该模型可以用来预测钢筋混凝土靶板发生厚度方向开裂破坏的临界能量。模型预测与实验结果吻合较好。     

Development and evaluation of an in-situ tire testing facility with variable side slip angles

An in-situ tire test rig was developed for field research on tire tractive and maneuverability performances. The Single Wheel Tester (SWT) was mounted on a tractor and a tested wheel was driven by a hydromotor, along a frame of 3 m length. In the SWT, four load cells were utilized to measure longitudinal and lateral forces, input and self-aligning torques, and two optical counters were applied to calculate forward and angular velocities. Response Surface Methodology was used to execute experimental design and to analyze the collected data. Afterwards, reduced form of a 2 Factor Interaction model was extracted to predict rolling resistance using seven factors. The test results show that increasing the normal load and side slip angle will cause an increment of rolling resistance. The incremental growth rate of the rolling resistance due to the normal load increment was observed. At higher cone index values, increasing the angular velocity reduces the rolling resistance, although at lower cone index values, the effect of angular velocity on rolling resistance is in reverse order. In addition, the increasing moisture content effect on rolling resistance at lower side slip angle values was observed.     

Aerodynamic force measurement using 3-component accelerometer force balance system in a hypersonic shock tunnel

A new three-component accelerometer force balance has been designed, calibrated and tested in hypersonic shock tunnel (HST2) of Indian Institute of Science. The newly designed balance is able to measure aerodynamic forces (within test time of one millisecond) on test models at angles of attack from 0 to 12°. Two models, a blunt cone with after body and a blunt cone with after body and frustum are used to establish the accuracy of the force balance. The tests were conducted for the above two configurations with a constant Mach number of 8 and total enthalpy of 2.0 MJ/kg. The effectiveness of the balance is demonstrated by comparing the forces and moments of measured data with AGARD models. The flow fields around the test model are simulated using a 3D axisymmetric Navier–Stokes solver and the simulated results were compared with the measured values. Measured and computed force data are matched within ±10% for two different models tested here. The accuracy of the force balance is also estimated with the Newtonian theory and the values are approximately ±10% for the axial component and ±8% for the normal and pitching moment components.      

水位变化对正倒锥体冰载荷影响的离散元分析

在海冰与锥体海洋结构的相互作用过程中,潮汐水位变化时海冰作用于锥体结构的位置改变对冰载荷具有显著影响.本文采用具有粘结破碎功能的离散元方法计算海冰与锥体作用的破坏过程.同时考虑海冰上下表面温度差异对海冰强度的影响,将离散元计算冰载荷及海冰破坏模式与渤海现场实测数据进行对比验证.离散元结果表明,海冰与正锥和倒锥碰撞时均发生弯曲破坏,且冰载荷均随水线处锥径的增大而增大.在水线处锥径相同的情况下,正锥冰载荷大于倒锥冰载荷,而正锥作用下海冰的断裂长度则较小.基于离散元计算结果和渤海现场观测资料分析了海冰与正锥、倒锥作用时冰载荷和断裂长度差异的主要原因.海冰与正倒锥交界线处作用时,一般发生弯曲破坏.当冰层中心高度与正倒锥交界线的高度相同时,海冰才会发生局部挤压破碎,但冰荷载并没有明显升高.由此可见,倒锥体结构可有效降低冰载荷从而具有较好的抗冰性能.以上研究表明离散元方法可确定海冰与锥体结构作用时的海冰破碎规律和冰载荷特性,为海洋工程结构的抗冰设计提供参考依据.      

飞机模型高速撞击钢筋混凝土载荷特性实验研究

为得到大型商用飞机撞击的冲击载荷特性及其计算方法,基于火箭橇加载试验平台,搭建了飞机模型撞击钢筋混凝土运动靶体测试系统,开展了两种不同尺寸飞机模型的撞击试验,利用高速摄影技术获得了飞机模型撞靶前的姿态、着速及飞机模型撞靶的破坏过程。采用加速度测试系统和激光干涉测速系统,分别得到撞击过程中运动靶体的加速度和速度历史,据此得到靶体受到的冲击载荷-时间曲线,二者吻合,验证了测试系统的可靠性。通过飞机模型上安装的机载存储过载测试系统,获得了撞击过程中飞机模型的负加速度-时间曲线,由此确定了修正的Riera理论模型中的静载荷项,并进一步计算得到靶体受到的冲击载荷-时间曲线,与通过测试靶体的加速度或速度得到的结果一致,验证了使用修正的Riera理论模型计算飞机模型冲击载荷的合理性及计算方法的正确性,同时确定了适合于本试验中飞机模型的修正系数α。     

Tractor tire aspect ratio effects on soil bulk density and cone index

A 580/70R38 tractor drive tire with an aspect ratio of 0.756 and a 650/75R32 tire with an aspect ratio of 0.804 were operated at two dynamic loads and two inflation pressures on a sandy loam and a clay loam with loose soil above a hardpan. Soil bulk density and cone index were measured just above the hardpan beneath the centerline and edge of the tires. The bulk densities were essentially equal for the two tires and cone indices were also essentially equal for the two tires. Soil bulk density and cone index increased with increasing dynamic load at constant inflation pressure, and with increasing inflation pressure at constant dynamic load. In comparisons of the centerline and edge locations, soil bulk density and cone index were significantly less beneath the edge than beneath the centerline of the tires. Soil compaction is not likely to be affected by the aspect ratio of radial-ply tractor drive tires when aspect ratios are between 0.75 and 0.80.     

岩石含水蠕变损伤模型的开发与应用

摘要：为研究岩石含水率对蠕变的影响,将含水开关与蠕变损伤阀值引入蠕变模型,实现了含水劣化与蠕变损伤在本构关系上的耦合,建立了岩石含水蠕变损伤模型。模型与某软岩的三轴压缩蠕变试验结果取得了较好的拟合效果,并由麦夸特法+通用全局优化（LM-UGO）算法反演出了模型参数随岩石含水率的变化规律。在C++语言和FLAC3D中Fish语言的环境下,通过推导模型本构方程的三维差分形式,并利用FLAC3D所提供的接口实现了模型的二次开发。通过对开发的岩石含水蠕变损伤模型进行算例验证和分析,验证了模型的可靠性,也得出了软岩巷道蠕变变形随围岩含水率增大而增加的结论,进一步反映出蠕变模拟计算中考虑岩石含水率的必要性。     

3D macro and micro-block models for limit analysis of out-of-plane loaded masonry walls with non-associative Coulomb friction

In this paper, a masonry system composed of a façade wall connected with two sidewalls and subjected to out-of-plane loading is investigated within the framework of three-dimensional limit analysis. Two different modeling approaches, namely macro and micro-block models, are adopted. A rigid-perfectly plastic model with dry contact interfaces governed by Coulomb failure criterion is assumed for masonry walls with regular units and staggering (non-standard limit analysis). Three classes of failure modes are investigated, involving rocking, sliding, twisting failure and combinations of them. The macro-block model is based on the assumption that the failure involves a number of cracks which separate the structure into a few macro-blocks and all the possible relative motions among micro-blocks are concentrated along the cracks. Two limiting conditions for the ultimate load factor are kinematically computed by use of minimization routines. The micro-block model is based on a concave contact formulation in which contact points are located at the corners of interfaces, allowing failure modes involving opening and sliding to be simulated. An iterative solution procedure is used to solve the non-associative friction problem, with second order cone programming (SOCP) used to allow the conic yield function to be solved directly. Both models are validated against experimental outcomes from the literature. A parametric analysis is carried out in order to highlight the influence of each geometrical and mechanical parameter on the prevalence of a mechanism over the other. The presence of an unrestrained horizontal floor system with different orientations is also analyzed.