六类有限元上机题.docVIP

注：题中 E表示材料弹性模量，μ表示泊松比，ρ表示密度。 一 静力结构分析 1 如图1所示为普通订书钉，E=2.1×105MPa，μ=0.3，横截面尺寸为宽B=0.64mm，高H=0.402mm。当订书钉被压入纸张时，约需要120N的载荷，载荷均匀地分布在订书钉上部。以下面两种情况进行有限元分析。（单位：mm） （1）订入时A、B点为铰支条件； （2）订入时A、B点为固定约束。 图1 载荷和尺寸情况 2、小型铁路桥由横截面积为3250mm2的钢制杆件组装而成。一辆火车停在桥上，其载荷施加在桥梁两侧的桁架上，单侧的桁架如图2所示，等效载荷为F1，F2，E=2.1×105MPa，μ=0.3，ρ=7.8×103kg/m3。试计算位置R处由于载荷作用而沿水平方向移动的距离以及支反力，同时，分析各个节点的位移和非单元应力。 图2 铁路桥单侧桁架及载荷情况 3如图3所示，模型参数为：E=3.0×1010Pa，A1=30m，A2=10m，B=80m，t=20m，p=2200Pa。 有关风载的确定，按照海洋井架行业标准，有以下方法： 风压（Pa）=0.6115×风速（m/s）×高度系数×形状系数 对于一般的海洋井架及建筑物，高度在30m左右，高度系数取为1.1，形状系数取为1.25，风速取为47.8m/s。换算出来后得到的风压为2200Pa。 图3 高层建筑物及其风载荷 4对于含裂纹体的结构及材料，若按照线弹性力学分析，会在裂纹的尖端处产生应力的奇异性，这时需要计算裂纹尖端处的应力强度因子（对于Ⅰ型裂纹，有K1=σ（πa）1/2），并以应力强度因子作为准则来对材料的裂纹是否扩展进行判断。图4所示为一块矩形平板，其边缘存在长为a的裂纹，板的两端承受拉应力σ0作用。利用结构上下的对成性，取矩形的一半建立有限元模型，完成看一下工作： 球裂纹的张角θ（在施加载荷前=0，θ=0） 沿直线AO，画出y方向应力σy沿x变化的曲线图。假设 σy= K1/(2πx)1/2利用回归方法估计K1。将计算结果与计算无限大平板的修正结果进行比较，其中基于无限大平板的修正关系式为K1=1.2σ0（πa）1/2） 在裂缝尖端处，进行网格的细化，重新求解（2）中的问题。 模型中的相关参数为 E=2.1×105MPa，μ=0.3，L=400mm，a=9.5mm，b=95mm，σ0=450MPa。 图4 承受拉升载荷的带裂纹平面结构 5、实体模型如图所示为一正六边形螺钉板手，在手柄端部加斜向上的面上施加扭转力100N，然后加向下的面力20N，计算出在这两种外载作用下扳手的应力分布。详细参数如下： 截面宽：10mm，正六边形边长为5.8mm，杆长：75mm，手柄长：200mm， 弯曲半径：10mm，弹性模量：2.07e11Pa 二 传热问题分析 6、 如图5所示为一桁架结构，分析下列两种情形下的节点位移和单元应力： （1）构件1、3、7和8的温度升高50°C。 （2）由于制造误差，构件9和10短了0.63mm，而6长了0.27mm，但必须进行强制装配。 桁架所用材料系数：E=200GPa，线膨胀系数α=1.25×10﹣5（/°C）。每个构件的横截面积均为100mm2。 图5 10杆件组成的桁架结构 7、受热载荷作用的后壁圆筒（图6a）的有限元建模和温度场求解。分析可知该问题属于轴对称问题，利用轴对称结构的特性，简化计算模型如图6b所示。 图6 左为a实体结构 右为b计算模型 三 瞬态动力学分析 8、如图7a所示为板-梁结构，板件上表面施加随时间变化的均布压力，计算在下列已知条件下结构的瞬态响应情况。 全部采用A3钢材料，特性：E=2×1011Pa，μ=0.3，ρ=7.8e3Kg/ m3，板壳厚度=0.02m。四条腿（梁）的几何特性：截面面积=2e-4 m2，惯性矩=2e-8 m4，宽度=0.01m，高度=0.02m。压力载荷与时间的关系曲线见下图7b所示。 图7a 质量梁-板结构及载荷示意图 压力（Pa） 10000 5000 0 1 2 4 6 时间（s） 图7b 板上压力-时间关系 四 接触分析 9、在这个实例中，我们将对一个弹簧卡子（如图8）进行接触分析，计算将卡头压进卡座和拉出卡座所需要的力。此问题属于需要输入厚度的平面应力问题，卡头和卡座的底板被认为是刚性的，因此在建模时不以考虑。 由于模型和载荷都是对称的，因此可用模型的右半部来进行计算。求解通过二个载荷步实现。 材料性质：E=2.8×103MPa (杨氏模量），μ=0.3，