算例(反分析).doc

洞室开挖反分析 问题描述： 地下洞室的开挖，开挖边界与开挖影响区域取为如图所示。在开挖后测点P1（198.90,119.0）,P2(183.4,75.0),P3(214.4,103.5),开挖后通过监测发现测线1收缩0.05m，测线2收缩0.05m，测线3收缩0.08m，所用的地层参数为E＝200000kPa，0.4，kN/m3,c=12kPa，，5kPa,k=0.46。通过量测信息，反演地层的力学参数弹性模量E。 解题步骤： 定义材料 从主菜单上选择：有限元分析—材料； 单击添加，将新添加的材料改名为地层； 输入地层参数为E＝200000kPa，0.4，kN/m3,c=12kPa，，5kPa,k=0.46。 绘制几何模型 命令：rectangle↙ 第一点：110.04，19↙ 第二点：301.70，181.6↙ 命令：line↙ 第一点：183.4，75↙ 下一点：214.5，75↙ 下一点：214.4，103.5↙ 下一点：↙ 命令：line↙ 第一点：183.4，75↙ 下一点：183.4，103.5↙ 下一点：↙ 命令：point↙ 指定点的位置：198.90,119.0 指定点的位置：183.4，103.5 指定点的位置：214.4，103.5 指定点的位置：↙ 从绘制菜单上选择圆弧，3点画圆弧模式。 命令：arc 指定圆弧的起点，或者[圆心Center]:单击P3点 指定圆弧的第二点，或者[圆心Center]:单击P1点 指定圆弧的终点：单击P2点。这样图形建模完毕。 保存文件。 寻找封闭面，赋土层材料 单击命令 网格－自动寻找封闭面，完成封闭面的寻找。 选中封闭面，单击属性工具，弹出修改属性对话框。 将类型改为地层。 施加边界条件 单击工具条，施加边界条件，选中两竖直边施加水平方向上的固定约束。 选择底边施加垂直方向上的固定约束。 定义施工步与开挖过程。 单击工具条，弹出对话框定义为一个开挖步。 到第1施工步，将洞室挖去。 定义测线 定义测线前，必须先画好准备定义测线的普通直线。为此在点P1，P2，P3之间画直线。 将施工步调到施工步1 单击菜单 反分析—测线。 命令：mealine 选择作为测线的直线或圆弧线段，及作为测点的点：单击P2,P1,Line1, 删除原线段？[Yes/No]NY 在选择时也可在工程窗口中选中对应的点与直线 通过同样的步骤可以定义其他的两条测线。 完成测线定义后的工程窗口 修改测线1的属性，选中测线1，单击工具条，弹出属性对话框， 单击测量组数，详细信息，弹出如下对话框。选中施工步1添加到当前测量组。 单击测量数据，详细信息，弹出如下对话框。 在测点1处输入量测信息为收缩0.05m。 选中测线2与测线3重复上述步骤，修改测线信息。测线2收缩0.05m，测线3收缩0.08m。保存文件。 定义反分析起始步与终止步 单击 反分析――反分析起始步与终止步。 反分析控制 定义反分析目标为地层的弹性模量E。 反分析求解 单击网格－全自动生成网格。 单击反分析－反分析求解。出现反分析求解过程跟踪图， 反分析跟踪图与所选的优化分析方法密切相关。跟踪图会因所选的分析方法而不同。 查看计算结果 计算完毕后在命令窗口，可以看到反分析的结果。 原始值 反演值 反分析计算数据文件格式说明 反分析计算过程中，除生成.DAT文件外（文件格式见《同济曙光有限元正分析手册》），还将生成.MEA和.BAA文件，如果反分析优化方法不是单纯形法，还有.INP文件。以下对这些文件格式进行说明。 声明：计算文件格式仅供用户参考，由于计算文件是由同济曙光软件自动生成，因此实际使用的计算文件格式可能与以下说明的格式有所不同！ （1）.BAA文件 数据输入格式示例BACK(analysis) *READ(11,(A1,A))LL,TXT * * READ(TXT,*) K,(IDR(K,I),I=1,4),VAR0(K,1),VAR0(K,2) * * IF (IBFLAG.GT.0) THEN ！INPUT BOUNDARY LOADS PARAMETERS * READ(11,(A1,A))LL,TXT * READ(TXT,*) K,(NBCXY(K,I),I=1,5),BCXY(K,J),J=1,6) * ENDIF *----+--------------+----+--------------+----+-------+----+------+----+ 1 1 2 1 0 60000.00 0.0 2 1 3 1 1 5500.00 1500.0 E 3 1 5 1