柱底反力求和lisp软件_基于AutoLisp语言的RC异形柱下独立基础验算方法与流程-CSDN博客

本文链接：https://blog.csdn.net/weixin_39692761/article/details/111765754

本发明涉及独立基础设计技术领域，尤其是把cad二次开发软件技术与基础验算技术相结合，具体为一种基于autolisp语言的rc异形柱下独立基础验算方法。

背景技术：

2017年12月1日正式实施的《混凝土异形柱结构技术规程》进一步明确了异形柱的设计方法和原则，不仅增加了z形柱的相关规定，同时提出楼梯间也可以采用肢端设暗柱的异形柱，极大的扩展了异形柱的使用范围。但是，异形柱下独立基础的设计在《建筑地基基础设计规范》(以下简称《地规》)等相关基础设计规范中仍没有明确规定。

现有的结构设计软件中仅有yjk和理正工具箱可对异形柱下独立基础进行验算，其中yjk将异形柱简化成外包矩形柱进行验算，设计方法过于简略与实际情况误差较大。理正工具箱虽然在7.0以后的版本中推出了相关功能，但在验算方法和原则方面存在以下问题：(1)理正工具箱(7.0pb4sp1版)将基础受冲切验算分为“整体破坏模式”验算和“局部破坏模式”验算，每种破坏模式按冲切锥底面形状又分为“真实截面”和“外包截面”两种形式，因此一般情况下理正需要验算4种破坏情况。但是我国《地规》在验算独基时，并不存在“整体破坏模式”，该种破坏模式仅在柱子对板或筏板的冲切中需要验算；(2)“整体破坏模式”的计算方法有误，其公式仍采用《地规》第8.2.8条中独基的冲切验算公式，仅在计算冲切周长时与“局部破坏模式”有所区别，这种计算缺乏理论与试验根据，如需计算“整体破坏模式”，应按照柱对板或筏板的冲切计算方法进行计算；(3)破坏模式下的“真实截面”破坏形态不合理，其与各国规范假定的破坏形态和试验结果均不相符；(4)冲切力的计算不符合《地规》第8.2.8条的规定，仅采用冲切阴影区域内的冲切力而未采用冲切基底净反力最大值乘以冲切阴影面积，降低了基础的安全储备；(5)未验算柱边宽高比，不符合规范要求；(6)不能对z形柱下独立基础进行验算。

同时，由于柱下独立基础局部破坏的特点，手算不仅工作量大容易出错，且难以形成计算书供他人检验，这使得异形柱下独立基础的设计一直是异形柱结构设计中的痛点和难点。因此有必要发明一种适用于异形柱下独立基础冲切及剪切验算方法，其计算模型应符合实际及规范要求，计算结果应精确可靠，同时易于操作，便于设计人员掌握，并生成计算书，便于设计人员检验，从而提高设计效率。

技术实现要素：

为解决现有技术存在的上述缺陷，本发明提出了一种适用于矩形柱、十字形、l形、t形、z形柱下扩展基础冲切和剪切验算的方法，并采用dcl语言和autolisp语言编写人机对话程序、计算程序和结果输出程序，能在无规范可循的情况下实现对异形柱下独立基础的验算，通过简单的操作实现所述基础的设计，满足安全性和经济性要求。

为实现上述目的，本发明提出了一种基于autolisp语言的rc异形柱下独立基础验算方法，包括如下步骤：

s1：在autocad软件运行的情况下，打开可视化输入程序，所述可视化输入程序通过dlc语言编写，通过可视化输入程序输入所述基础的外力、混凝土强度等级、柱及基础角点坐标信息；

s2：根据所述输入的信息进行基础验算，包括如下小步：

s21：坐标变换：以步骤s1中得到的矩形基础角点为基准建立局部坐标系，以基础形心为原点，以通过形心平行于基础边的直线为轴，将基础及柱在整体坐标系中的坐标转换为局部坐标；

s22：识别柱子类型及基础边界：根据柱子对象中角点的数量及x、y坐标相同的点的个数，识别autocad图中柱子的种类；利用函数maxmin确定x、y坐标的上界和下界，为后续计算的判定提供依据；

s23：计算有效加载面：根据柱子类型和基础边界，剔除异形柱的内凹角点，将剩余角点坐标依次按顺序储存，即得到有效加载面各顶点；有效加载面定义为由异形柱各外凸顶点依次相连得到的面，其顶点数据是计算冲切锥底面的基础；

s24：验算柱边宽高比：采用异形柱外包矩形作为验算对象，若验算外包矩形各边柱边宽高比大于2.5，则程序将退出并要求设计人员重新设计基础尺寸；

s25：初步判断基底反力零应力区面积是否满足要求，当满足要求时可进一步进行计算，否则退出程序，提示设计人员增大基础尺寸；

s26：计算基底反力分布；

s27：计算零应力区面积百分比：通过所述基底反力分布，精确计算基底反力零应力区面积占基础底面积的比例，进一步验证零应力区面积是否满足要求；

s28：冲切及剪切承载力验算：以所述有效加载面的一条边为对象，计算其对应的受冲切截面面积、冲切阴影面积、冲切标记、受剪切截面面积、剪切阴影面积，从而验算基础的冲切承载力和剪切承载力，当所有边的冲切承载力和剪切承载力均满足要求时，所设计的基础才满足要求；

s3：验算结果输出：根据所述步骤s2输出计算结果，当柱边宽高比或零应力区面积百分比不满足要求时计算终止，提示设计人员修改基础尺寸，否则生成计算书。

进一步地，所述异形柱为矩形柱、l形柱、t形柱、十字形柱、z形柱中的一种；矩形柱设置有4个角点；t形设置有8个角点，其中有4个点具有相同的x或y坐标；z形设置有8个角点，其中有3个点具有相同的x或y坐标；十字形设置有12个角点。

进一步地，所述步骤s1中可视化输入程序右侧列出了适用条件，左侧提示设计人输入基础所受的外力信息；通过对话框下拉对话框选择混凝土强度等级，强度范围从c25～c60；通过对话框中“选取基础pl线

进一步地，所述步骤s21中的坐标变换通过如下公式得到：

若局部坐标系相对于整体坐标系逆时针旋转θ，局部坐标系原点在整体坐标系中的坐标为(x0，y0)时，得到局部坐标系与整体坐标系之间的关系为：

若局部坐标系相对于整体坐标系顺时针旋转θ，将-θ代入以上公式，得到局部坐标系与整体坐标系之间的关系为：

进一步地，所述步骤s23中内凹角点的特点是：角点的x坐标和y坐标既不是所有角点中坐标最大的值，也不是最小的值。

进一步地，所述步骤s25中的初步判断基底反力零应力区面积是否满足要求，根据如下公式进行判断：

式中，为x方向的偏心距比，为y方向的偏心距比；

当上述公式成立时则表示零应力区面积不超过基础底面积的15％，满足《建筑地基基础设计规范》要求，否则不满足要求。

进一步地，所述步骤s26中计算基底反力分布，所述基底反力存在以下三种情况：

情况一：当基底不存在零应力区时，基底净反力计算如下：

式中，px,y代表：基底反力设计值，n代表：柱轴力设计值；a代表：基础底面积；mx代表：基础底面绕x轴的弯矩设计值；my代表：基础底面绕y轴的弯矩设计值；ix代表：基础底面绕x轴的惯性矩；iy代表：基础底面绕y轴的惯性矩；

情况二：当基础单向受弯，基底出现零应力区时，基底反力的平均值和最大值计算如下：

式中，pave代表：基底反力平均值；pmax代表：基底反力最大值；b代表：垂直于弯矩方向的基础边长；e代表：外力的偏心距；l代表：平行于弯矩方向的基础边长；

情况三：当基础双向受弯，基底存在零应力区时，零应力区根据面积不同分为：情况a、情况b、情况c和情况d：

情况a：基础底面大面积与地基脱开，基础有三个角点应力为0，显然不能满足规范中关于基底零应力区的要求，提示设计人员重新设计基础尺寸；

情况b：基础底面大面积与地基部分脱开，基础有两个角点应力为0，未知量及基底反力最大值计算如下：

式中，l1和b1代表：与基底反力分布相关的几何参数；pmax代表：基底反力最大值；ex代表：距y轴的偏心距；ey代表：距x轴的偏心距；b代表：基础x方向宽度；l代表：基础y方向宽度；

情况c：基础底面小面积与地基部分脱开，基础有两个角点应力为0，未知量及基底反力最大值计算如下：

情况d：基础底面小面积与地基部分脱开，基础只有一个角点应力为0，采用newton法求解，计算如下：

对于非线性方程组：

令x＝(x1,x2)t，f＝(f1,f2)t，在精确解附近有近似根根据taylor公式有：

f(x)≈f(x(k))+f′(x(k))(x-x(k))

其中：

可得迭代方程：

x(k+1)＝x(k)-f′(x(k))-1f(x(k))

其中：

newton法能收敛到真实解，初值计算如下：

进一步地，所述步骤s28中冲切和剪切承载力按照《建筑地基基础设计规范gb50007-2011》第8.2.7条～第8.2.9条的规定进行验算。

进一步地，所述步骤s3中计算书为txt格式，生成在程序所在文件夹内，内容包括基础尺寸、柱子类型及尺寸、基础柱边宽高比验算过程及结果、基础的材料属性、基础外力、基础反力平均值、基础反力最大值、基础反力最小值或基础零应力区面积占基底面积的百分比，有效加载面各边对应的冲切和剪切承载力验算过程及结果以及结论。

进一步地，所述验算方法包括可视化输入模块、基础验算模块和结果输出模块，所述可视化输入模块采用dcl语言编写，所述基础验算模块和结果输出模块采用atuolisp语言编写，在autocad软件运行的情况下打开可视化输入程序，输入基础基本信息后，点选autocad软件中的图形信息，程序自动运行基础验算模块对所述基础进行验算并生成计算结果文件。

本发明的有益效果是：(1)本发明所采用的验算步骤，以及设计原则以规范为依据，简单合理清晰明了，依据此设计出的异形柱下基础既能保证安全，又经济合理；

(2)本发明可仅通过外力和基础平面尺寸，即可判断基础是否满足规范中关于零应力区面积大小的要求，为设计人员快速布置基础尺寸提供帮助；

(3)本发明可对十字型、l型、t型、z型柱下独立基础的冲切、剪切承载力进行验算，同样适用于形状类似的墙下独立基础的验算，适用范围广泛；

(4)本发明可与cad图形文件互动，读取柱及基础信息，简单方便，并拥有强大的提示功能，使得设计人员能够轻松上手，可以大大提高设计人员的工作效率；

(5)本发明提供的计算书清晰明了，可供设计人员进行校对。

附图说明

图1是本发明的原理流程框图。

图2(a)-图2(d)是异形柱的四种结构示意图。

图3是本发明与cad软件交互窗口。

图4(a)-图4(d)是四种基础双向受弯时基底零应力区的图。

图5是偏心距公式及数值分析点组成的临界曲线图。

图中：1、冲切锥底面；2、有效加载面。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

实施例1：

如图1所示，本发明所述的基于autolisp语言的rc异形柱下独立基础验算方法，本发明采用的计算步骤，参考了美国规范的相关规定，确定了异形柱冲切底面的形状，在计算过程中考察柱边宽高比保证了基底反力线性分布假定的成立，并采用数值方法精确计算基底零应力区面积，在承载力验算时，对有效加载面的每一边进行验算可保证设计的安全可靠。

本发明与一般验算方法不同的是，采用autolisp语言并通过dcl与设计人互动，可自动识别基础和柱子类型及尺寸信息，使得设计人在绘图的同时即完成计算，减小了设计人再次输入设计参数工作量。

本发明可以精确计算基底反力，当基础在双向受弯作用下基底存在零应力区时，基底反力计算较为复杂，无法得到解析解，本发明通过newton法求解非线性方程组，可准确计算基底反力，使得设计更加安全经济。

本发明提出了一种基于autolisp语言的rc异形柱下独立基础验算方法，包括如下步骤：

s1：在autocad软件运行的情况下，打开可视化输入程序，所述的可视化输入程序通过dlc语言编写，通过可视化窗口输入所述基础的外力、混凝土强度等级、柱及基础角点坐标信息；

s2：根据所述输入的信息进行基础验算，包括如下小步：

s24：验算柱边宽高比：采用异形柱外包矩形作为验算对象，若验算外包矩形各边柱边宽高比大于2.5，则程序将退出并要求设计人员重新设计基础尺寸；

s25：初步判断基底反力零应力区面积是否满足要求，当满足要求时可进一步进行计算，否则退出程序，提示设计人员增大基础尺寸；

s26：计算基底反力分布；

s27：计算零应力区面积百分比：通过所述的基底反力分布，精确计算基底反力零应力区面积占基础底面积的比例，进一步验证零应力区面积是否满足要求；

s28：冲切及剪切承载力验算：以所述的有效加载面的一条边为对象，计算其对应的受冲切截面面积、冲切阴影面积、冲切标记、受剪切截面面积、剪切阴影面积，从而验算基础的冲切承载力和剪切承载力，当所有边的冲切承载力和剪切承载力均满足要求时，所设计的基础才满足要求。

s3：验算结果输出：根据所述的步骤s2输出计算结果，当柱边宽高比或零应力区面积百分比不满足要求时计算终止，提示设计人员修改基础尺寸，否则生成计算书。

如图2(a)-图2(d)所示，所述的异形柱为矩形柱、l形柱、t形柱、十字形柱、z形柱中的一种；矩形柱设置有4个角点；t形设置有8个角点，其中有4个点具有相同的x或y坐标；z形设置有8个角点，其中有3个点具有相同的x或y坐标；十字形设置有12个角点。

如图3所示，所述步骤s1中可视化输入程序右侧列出了适用条件，左侧提示设计人输入基础所受的外力信息；通过对话框下拉对话框选择混凝土强度等级，强度范围从c25～c60；通过对话框中“选取基础pl线

所述步骤s21中的坐标变换通过如下公式得到：

若局部坐标系相对于整体坐标系逆时针旋转θ，局部坐标系原点在整体坐标系中的坐标为(x0，y0)时，得到局部坐标系与整体坐标系之间的关系为：

若局部坐标系相对于整体坐标系顺时针旋转θ，将-θ代入以上公式，得到局部坐标系与整体坐标系之间的关系为：

所述步骤s23中内凹角点的特点是：角点的x坐标和y坐标既不是所有角点中坐标最大的值，也不是最小的值。

所述步骤s25中的初步判断基底反力零应力区面积是否满足要求，根据如下公式进行判断：

式中，为x方向的偏心距比，为y方向的偏心距比。

当上述公式成立时则表示零应力区面积不超过基础底面积的15％，满足《建筑地基基础设计规范》要求，否则不满足要求。

所述步骤s26中计算基底反力分布，所述的基底反力存在以下三种情况：

情况一：当基底不存在零应力区时，基底净反力计算如下：

情况二：当基础单向受弯，基底出现零应力区时，基底反力的平均值和最大值计算如下：

情况三：当基础双向受弯，基底存在零应力区时，零应力区根据面积不同分为：情况a、情况b、情况c和情况d：

情况a：基础底面大面积与地基脱开，基础有三个角点应力为0，显然不能满足规范中关于基底零应力区的要求，提示设计人员重新设计基础尺寸；

情况b：基础底面大面积与地基部分脱开，基础有两个角点应力为0，未知量及基底反力最大值计算如下：

情况c：基础底面小面积与地基部分脱开，基础有两个角点应力为0，未知量及基底反力最大值计算如下：

情况d：基础底面小面积与地基部分脱开，基础只有一个角点应力为0，采用newton法求解，计算如下：

对于非线性方程组：

令x＝(x1,x2)t，f＝(f1,f2)t，在精确解附近有近似根根据taylor公式有：

f(x)≈f(x(k))+f′(x(k))(x-x(k))

其中：

可得迭代方程：

x(k+1)＝x(k)-f′(x(k))-1f(x(k))

其中：

newton法能收敛到真实解，初值计算如下：

所述的步骤s28中冲切和剪切承载力按照《建筑地基基础设计规范gb50007-2011》第8.2.7条～第8.2.9条的规定进行验算。

所述步骤s3中计算书为txt格式，生成在程序所在文件夹内，内容包括基础尺寸、柱子类型及尺寸、基础柱边宽高比验算过程及结果、基础的材料属性、基础外力、基础反力平均值、基础反力最大值、基础反力最小值或基础零应力区面积占基底面积的百分比，有效加载面各边对应的冲切和剪切承载力验算过程及结果以及结论。

进一步地，所述的验算方法包括可视化输入模块、基础验算模块和结果输出模块，所述的可视化输入模块采用dcl语言编写，所述的基础验算模块和结果输出模块采用autolisp语言编写，在autocad软件运行的情况下打开可视化输入程序，输入基础基本信息后，点选autocad软件中的图形信息，程序自动运行基础验算模块对所述的基础进行验算并生成计算结果文件。

实施例2：

如图1所示，本发明包括可视化输入模块、计算模块、输出模块。图1给出的流程图用于说明上述三个模块之间及子块和子块之间的逻辑关系。

可视化输入模块的功能在于，提示设计人员点选柱、基础的图形对象，通过对话框的形式，要求设计人员输入基础的高度，以及通过柱形心的轴力、弯矩和剪力的设计值，混凝土强度等级以及基础顶面的埋置深度，并将这些数据传递给计算模块。并告知设计人员本发明的适用范围以及所采取的假设。

计算模块的功能在于，将整体坐标系下柱和基础节点坐标，转换为以基础形心为原点建立的坐标系的局部坐标；根据输入模块传入的数据，验算各边柱边宽高比；计算通过基础形心的轴力、弯矩的设计值；验算零应力区是否满足《抗规》第4.2.4条；计算基础偏心情况；计算基底应力分布，给出阴影区域范围内的基底净反力最大值以及基底净反力平均值，当存在零应力区时，给出零应力区面积；验算有效加载面各边是否满足冲切承载力要求；验算有效加载面积各边是否满足剪切承载力要求；将输出的结果传递给结果输出模块。

结果输出模块的功能在于，根据计算模块传入的数据，输出txt格式的计算书。输出的计算结果包括，基础及柱子尺寸、外力设计值、钢筋及混凝土强度信息、柱子形心的局部坐标、基底反力最大值、平均值、最小值或零应力区面积比，以及各边冲切及剪切承载力的验算过程，并在计算中提示设计人员基础设计中存在的问题。

在以上三大模块中，加入了很多提示和报错功能，可以帮助设计人员了解程序所处的状态，这使得整个程序易于上手。

流程图中，列出了计算模块中主要应用到的函数，下面对这些函数的功能做一个简要的介绍。

a、坐标转换函数gtol。该函数以矩形基础为基准建立局部坐标系，以基础形心为原点，以通过形心平行于基础边的直线为轴，将基础及柱在整体坐标系中的坐标转换为局部坐标。当局部坐标系相对于整体坐标系逆时针旋转局部坐标系原点在整体坐标系中的坐标为(x0,y0)时，可以得到局部坐标系与整体坐标系之间的关系为

若局部坐标系相对于整体坐标系顺时针旋转θ，将-θ代入，可得转换关系为

b、柱子类型判别函数columntype。该函数用于判定柱子为矩形柱、l形柱、t1形柱(翼缘平行于x轴的t形柱)、t2形柱(翼缘平行于y轴的t形柱)、十字形柱、z形柱中的哪一种，为后续计算提供数据。

c、确定基础边界。用到可以确定一组点坐标的x坐标和y坐标范围的函数maxmin，确定基础范围后可以进一步为后续计算的判定提供依据。

d、剔除凹点函数elimpoint。用于剔除异形柱的凹点，求得有效加载面各边围成的多边形，以此为基础进行各边的冲切和受剪承载力计算。

e、柱边宽高比函数bhratio。用于计算与异形柱对应的矩形柱的柱边宽高比是否均满足2.5的要求。此条是基底反力成线性分布假设成立的前提，若不满足将终止程序，提示设计人员重新设计基础尺寸。判断零应力区，若零应力区的面积不满足规范要求，则结束程序，提示设计人员重新设计基础尺寸。

f、外力设计值函数exforce。用于将柱形心处的外力系转换为基础形心处的外力系，通过基础埋深与剪力的乘积考虑由剪力引起的附加弯矩。

g、基底反力函数reactforce。用于计算基底反力的分布，给出基底反力最大值、基底反力平均值、零应力区面积等三个值用于后续计算。

h、边对象函数edgeobj。以有效加载面积的一条边为对象，计算其对应的受冲切截面面积、冲切阴影面积、冲切标记(是否需要冲切验算的记号)、受剪切截面面积、剪切阴影面积等五个数据，将这些数据传递给验算冲切v1check和剪切v2check的函数，以验证该边是否满足承载力要求。

实施例3：

本发明在基础底面存在零应力区时的基底反力分布是本发明的研究重点和难点，为了求解其反力分布情况，首先应考察基础底面不存在零应力区时，以及存在零应力区时单向受弯基础的基底反力分布特点。

第一种，针对基础底面不存零应力区时的基底反力分布的情况。

此情况下，基础底面的反力分布满足线性叠加原理，基础底面某一点的反力等于轴力、x向弯矩和y向弯矩在该点引起的反力之和，以基础底面形心为原点建立的坐标系中，基础底面内任一点的应力表达式为：

其成立条件为：

即：

其中，b和l分别为对应于ex和ey的基础边长。

第二种，存在零应力区时单向受弯基础的基底反力分布的情况。

此情况下，基底反力成三角形分布，其合力的位置在距基础形心e的位置，依次建立方程可以得到：

可得：

上式中各参数的含义详见《地规》第5.2.2条。从上述公式中可以看出，零应力区的面积与外力模式相耦合，线性叠加原理不再适用，这种情况在双向受弯时变得更加复杂。

第三种，存在零应力区时双向受弯基础的基底反力分布的情况。

以下讨论基于基底反力最大值出现在第一象限，对于出现在其他象限的情况，实际应用中偏心距可取绝对值。此处只讨论下述条件成立时的情况：

当双向受弯基础的基底存在零应力区时，基底反力的分布存在三类情况。第一种情况(情况a)，基础底面大面积与地基脱开，基础有三个角点应力为0，基底反力在空间中的形状为四面体，在基础底面的投影为直角三角形，其形心在基础底面的投影坐标为(ex,ey)，与三角形的形心重合，如图4(a)所示。此时直角三角形的两直角边，应分别小于相应基础边的边长。即：

可得此情况成立的充分必要条件为：

注意到此种情况无法满足规范要求，若在程序计算中出现此情况，应提示设计人员重新设计基础尺寸。

第二类情况为基础底面与地基部分脱开，基础有两个角点应力为0，基底反力在空间中的形状为平放的三棱台，如图4(b)和图4(c)所示。以图4(c)所示情况c为例，根据三棱台的体积等于轴力，形心在基底的投影坐标为(ex,ey)，可得如下方程：

ey所在方程为l1的二次方程，通过首先对其求解，可以求解整个方程组。

同样的，对于图4(b)所示情况b，方程组为：

方程的解为

第三类情况为基础只有一个角点应力为0，基础反力在空间中的形状为切掉一角的三棱台，如图4(d)所示情况d。根据力的平衡关系及几何关系，可以得到如下方程：

上述方程组为四次二元方程，难以直接求解，通过迭代函数对其进行数值求解。

第二类情况(情况b、c)和第三类情况(情况d)的适用条件，显然当基础反力在基础底面的投影的一条直角边与基础的一条边重合时，处于第二类及第三类情况的临界状态。

情况c与情况d的临界条件为，对于情况c为b1＝b，对于情况d为b1＝0，这时两组表达式在形式上完全一致，如下式所示：

解上述方程，临界偏心距(ex0,ey0)之间有如下关系，如图5所示：

令上式可化简为：