【有限元分析】apdl建立菩萨三维模型+FLAC3d力学计算

稳定性计算采用如下过程完成，首先对三维扫描数据进行提取，建立菩萨的结构模型，然后在Ansys中通过APDL程序编程建立起三维数值模型，在FLAC3d中进行三维力学计算。

1. 计算参数及模型建立
   对拟开展保护的胁侍菩萨进行了三维扫描，其三维图如下图左。以此在Ansys中建立三维模型，模型外轮廓线为下图右。在建立三维模型过程中，为降低模型的复杂程度，对局部装饰等部位进行适当简化。
   

图 三维扫描结果及结构计算控制线图
通过Ansys编程建立模型后，进行剖分，剖分结果如下图所示。总共剖分单元数为157291个，节点数为30394个。建模程序附后。
2. 计算结果分析

图 剖分结果图
图6.12 裂隙垂直方向分布
6.3 奉国寺大佛稳定性数分析

1. 计算参数与模型建立
   利用三维数据可以比较准确的建立出大佛的外表面的三维模型，但为了降低模型复杂程度，对于部分细节，如手部手指、头部冠饰等，就行了简化处理。而内部木骨架结构主要根据调查资料而建立，由于调查资料非常有限，部分结构就行推测得到。
   表面形态主要通过对模型进行简化。
   在上述工作基础上，得到了整个三维模型。对三维数值模型进行剖分，剖分后的模型单元总数为637239，节点数162588。

图6.14 内部受力木骨架结构

图6.15 木骨架整体结构图

图6.16 模型线结构图

图6.17 模型俯视图（Ansys中）
ANSYS 模型程序
/prep7
RatioSize=1
SoilRatio=0.1
WoodRatio=0.6
SoilSize=RatioSizeSoilRatio
WoodSize=RatioSizeWoodRatio

WoodPartNum=153 !定义木构件的数量，即行数
WoodPara=8 !每个木构件的参数数量，即列数

*dim,RotYZWPart,array,WoodPartNum !定义绕x轴转动角度（YZ面转）
*dim,RotZXWPart,array,WoodPartNum !定义绕y轴转动角度（ZX面转）
*dim,XCentWPart,array,WoodPartNum !定义中心点x坐标
*dim,YCentWPart,array,WoodPartNum !定义中心点y坐标
*dim,ZCentWPart,array,WoodPartNum !定义中心点z坐标
*dim,WXWPart,array,WoodPartNum !定义截面宽度
*dim,HYWPart,array,WoodPartNum !定义截面高度
*dim,DZWPart,array,WoodPartNum !定义柱体长度
*dim,WPartVN,array,WoodPartNum
*dim,UILX,array,18
*dim,UILY,array,18
*dim,UOLX,array,18
*dim,UOLY,array,18
*dim,DILX,array,18
*dim,DILY,array,18
*dim,DOLX,array,18
*dim,DOLY,array,18
*dim,K_KUI,array,18
*dim,K_KUO,array,18
*dim,K_KDI,array,18
*dim,K_KDO,array,18
*dim,AITemp,array,18
*dim,AOTemp,array,18
*dim,ATopTemp,array,18
*dim,ABotTemp,array,18
全部内容见：
CAE有限元仿真服务