资源描述：
《程序使用说明及举例》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、二维弹性问题常量边界元程序使用说明及举例一、程序使用说明1本程序用来求解二维弹性(平面应变)问题，程序改编自参考资料3。2本程序的流程图如下：主程序INPUTECGHMATECCCCCCCCCCSLNPDINTERECOUTPTEC结束输入原始数据形成H和G矩阵，并根据边界条件重新排列，形成A矩阵和F向量求解方程组计算内点的位移和应力输出结果在GHMATEC中调用子程序EXTINEC和LOCINEC；EXTINEC用高斯数值积分法计算H和G矩阵副元素；LOCINEC计算H和G矩阵的对角元素。3输入数据（1）TITLE字符型的标题（2）N,L,M,(NC(k),k=1,5),
2、GE,XNUN:边界单元总数L:需计算位移与应力的内点数M:不同边界个数NC(k),k=1,5:不同边界最后结点号,不足5个以0填满GE:剪切模量XNU:泊松比（3）X(I),Y(I),I=1-N:边界单元端点I的坐标值注：外边界单元编号方向为逆时针，内边界为顺时针。（4）KODE(2*I-1),FI(2*I-1),KODE(2*I),FI(2*I),I=1,N:边界单元I的边界条件KODE(2*I-1),KODE(2*I):边界单元I的x与y方向边界条件信息，KODE值为1时表示该方向p（面力）已知，KODE值为0时表示该方向u（位移）已知。FI(2*I-1),FI(2*
3、I):边界单元I的x方向和y方向的已知面力或位移具体数值。（5）CX(I),CY(I),I=1,L:内点各点坐标4其它标识符XM,YM:每个单元的中点（结点）的坐标；DFI:方程组的右端向量，解方程后成为解向量（每个单元的未知u值或p值）；DSOL:内点位移值；SSOL:内点应力值。二、例题说明外半径为0.9内半径为0.3的厚壁圆筒，受集度集度为200的径向均布内压力。考虑对称性，取圆筒四分之一，划分为36个常量元，如下图所示，求解其应力位移。三、例题输入数据例题输入数据如下elasticproblemsusingconstantboundaryelements367136
4、0000200000.1670.00000.30000.03920.29740.07760.28980.11480.27720.15000.25980.18260.23800.21210.21210.23800.18260.25980.15000.27720.11480.28980.07760.29740.03920.30000.00000.40000.00000.50000.00000.60000.00000.70000.00000.80000.00000.90000.00000.89230.11750.86930.23290.83150.34440.77940.4500
5、0.71400.54790.63640.63640.54790.71400.45000.77940.34440.83150.23290.86930.11750.89230.00000.90000.00000.80000.00000.70000.00000.60000.00000.50000.00000.4000113.08061199.5718139.01811196.1571164.28791189.3860188.45771179.37461111.11401166.29391131.86921150.36801150.36801131.86921166.2939111
6、1.11411179.3745188.45771189.3860164.28791196.1571139.01811199.5718113.080610.000000.000010.000000.000010.000000.000010.000000.000010.000000.000010.000000.000010.000010.000010.000010.000010.000010.000010.000010.000010.000010.000010.000010.000010.000010.000010.000010.000010.000010.000010.000
7、010.000010.000010.000010.000010.000000.000010.000000.000010.000000.000010.000000.000010.000000.000010.000000.000010.00000.24750.24750.28280.28280.35360.35360.42430.42430.49500.49500.56570.56570.60100.6010四、例题输出结果例题计算输出结果如下：*************************************
