s3fs 读取python python读取sas文件

转载

卫斯理 2023-11-20 00:09:49

文章标签 s3fs 读取python SAP2000 s2k Python 数据 文章分类 Python 后端开发

一、前言

当保存模型时，SAP2000 将自动将模型文本文件存成一个后缀为 .$2k 的文件，它与模型数据 .sdb 保存的同时进行保存。如果由于一些原因不能打开模型 .sdb文件时, 可以试着导入 .$2k 文件。.$2k 文本文件只是用于分析目的或者在特殊情况下复原模型，并且它不是数据库文件可替代文件。不包含所有模型信息。

从一个给定的 SAP2000 模型可以导出任意多个文件。每个文件可以包含不同的表格并且可以应用于模型不同部分。这个文件可以被其它程序使用，进行修改后再次导入到 SAP2000 程序中，或用于其它的目的。s2k 文件便是其中之一，s2k 文件即SAP2000 数据库文本 (.s2k) 文件，它是从 SAP 输出的纯文本文件。

.$2k 或 .s2k 文件是 SAP2000 文本文件，它们的格式是完全一致的，工程师可以用写字板、记事本等文字处理软件打开，对其中的内容进行修改，然后再用 SAP2000 导入，可以实现对模型的修改。

$2k 在模型数据 .sdb 保存时自动创建，s2k 需手动导出。

s3fs 读取python python读取sas文件_s3fs 读取python

s2k文件的导入

s3fs 读取python python读取sas文件_Python_02

s2k文件的创建

SAP2000版本：SAP2000 v22.0.0 。

二、读取SAP模型数据

从 .$2k 文件中可获取 SAP 模型的绝大部分信息，如点对象的标签和坐标，框架对象的起点标签、终点标签、截面指定和局部轴方向等，材料属性，截面属性等等。

本文以如下图所示的 SAP 模型为例，通过 Python 编程读取其 s2k 文件数据，并对模型数据进行处理，完成对模型基本信息的分类，生成相应的 csv 文件，以方便在其他软件中复建该模型。

s3fs 读取python python读取sas文件_s3fs 读取python_03

值得说明的是，为了适配更多版本的 s2k 文件，本文中很多 Pyhton 代码显得不是那么高明，甚至笨重啰嗦。但这些代码绝对逻辑清楚，考虑全面，通用性较强，容错能力优秀，提示信息简洁明了。如您有什么建议，还请不吝赐教，邮箱： liyang@alu.hit.edu.cn 。

本文程序以准确获取模型信息，适配尽可能多版本的 s2k 文件为主要目标，代码运行效率是否最优并不重要。速度不重要，准确最重要。

CADmep19x64.arx 此问题待解决。

#!/usr/bin/env python
# -*- coding: utf-8 -*-

'''
=============================
Author: DalNur 
Email: liyang@alu.hit.edu.cn
=============================
'''

import os
import re

s2kname = "GG-JB.s2k"  # s2k文件名称
currdir = os.getcwd()  # 当前工作目录
s2kpath = currdir + "\\" + s2kname

s2kpath = r"D:\GTC\GG-JB.s2k"  # s2k文件绝对路径

2.1 数据结构

从 s2k 文件内读取的各种模型信息，按类别存储在以下几个 Python 字典中，用于存储节点对象信息的字典 jointObjs ，用于存储框架对象信息的字典 frameObjs ，用于存储面对象信息的字典 areaObjs ，用于存储材料属性信息的字典 matProps ，用于存储框架对象截面指定信息的字典 frameSecProps ，用于存储面对象截面指定信息的字典 areaSecProps，用于存储框架对象 (均布) 荷载指定信息的字典 frameLoads，用于存储面对象 (均布) 荷载指定信息的字典 areaLoads，以及用于定位以上各种信息在 s2k 文件中的位置的字典 tables 。

字典名称	类别
tables	表格
matProps	材性
frameSecProps	框架截面
areaSecProps	面截面
jointObjs	点
frameObjs	线
areaObjs	面
frameLoads	框架荷载
areaLoads	面荷载

tables = {}  # 表格信息定位
matProps = {}  # 材料属性
frameSecProps, areaSecProps = {}, {}  # 框架截面属性/面截面属性
jointObjs, frameObjs, areaObjs = {}, {}, {}  # 点对象/框架对象/面对象
frameLoads, areaLoads = {}, {}  # 框架对象荷载/面对象荷载

表格信息定位字典的数据结构

tables = {'TABLE:  "PROGRAM CONTROL"': [4, 5], 
          'TABLE:  "ACTIVE DEGREES OF FREEDOM"': [8, 8],
          'TABLE:  "ANALYSIS OPTIONS"': [11, 11],

           . . . . . .
          
          'TABLE:  "OPTIONS - COLORS - DISPLAY"': [75786, 75794],
          'TABLE:  "OPTIONS - COLORS - OUTPUT"': [75797, 75805],
          'TABLE:  "DATABASE FORMAT TYPES"': [75808, 75808],
          'TABLE:  "PROJECT INFORMATION"': [75811, 75823]}

材料属性字典的数据结构

matProps = {'C30B': [1, 'Concrete', 'Isotropic', 2500.0, 30000000000.0, 12000000000.0, 0.25],

             . . . . . .

            'HPB300': [6, 'Rebar', 'Uniaxial', 7999.999999999999, 200000000000.0, 0, 0],

             . . . . . .

            'JTG-fpk1470': [8, 'Tendon', 'Uniaxial', 7849.999999999999, 195000000000.0, 0, 0],

             . . . . . .
            
            'Q235B': [9, 'Steel', 'Isotropic', 7799.999999999999, 206000000000.0, 79230769230.76921,0.3],

             . . . . . .

            'RIGIDBAR': [11, 'Other', 'Isotropic', 0.0, 2.8e+19, 1.4e+19, 0.0]}

框架截面字典的数据结构

frameSecProps = {
    'B1000X300X8X8XXQ80211B': [1, 'Q80211B', 'Box/Tube', 1.0, 0.3, 0.008, 0.008, 'Null', 'Null',
                               0.020544, 0.001046, 0.002451, 0.000372, 0.0, 0.016, 0.0048, 0.004903,
                               0.002478, 0.006254, 0.002659, 0.345423, 0.134507],

     . . . . . .
            
    'B600X200X10X10XXQ80211B': [29, 'Q80211B', 'Box/Tube', 0.6, 0.2, 0.01, 0.01, 'Null', 'Null',
                                0.0156, 0.000322, 0.000673, 0.000118, 0.0, 0.012, 0.004, 0.002244,
                                0.001181, 0.002862, 0.001302, 0.207754, 0.087016],

     . . . . . .
 
    'RIGIDBAR': [61, 'RIGIDBAR', 'Rectangular', 0.1, 0.1, 'Null', 'Null', 'Null', 'Null', 0.01,
                 1.4e-05, 8e-06, 8e-06, 0.0, 0.008333, 0.008333, 0.000167, 0.000167, 0.00025,
                 0.00025, 0.028868, 0.028868]}

面截面字典的数据结构

areaSecProps = {'RIG1000C35F': [1, 'C35F', 'Shell', 'Shell-Thick', 1.0],

                 . . . . . .

                'RIG200C35F': [6, 'C35F', 'Shell', 'Shell-Thick', 0.2],

                 . . . . . .

                'W500C40W': [11, 'C40W', 'Shell', 'Shell-Thick', 0.5]}

点对象字典的数据结构

jointObjs = {'2001902': [1, 'GLOBAL', 'Cartesian', 305.969, 91.824, 10.0],

              . . . . . .

             '4010545': [1682, 'GLOBAL', 'Cartesian', 93.291, 79.807, 0.928],

              . . . . . .

             '13001829': [3228, 'GLOBAL', 'Cartesian', 145.551, 55.207, 8.0]}

框架对象字典的数据结构

frameobjs = {'3000092': [1, '3000170', '3000171', 'B1200X600JXC35B', 0],

              . . . . . .
 
             '4009829': [1363, '4003302', '4005181', 'B1500X1200JXC35B', 0],

              . . . . . .

             '13005362': [4048, '13001829', '13000323', 'B800X400X25X25XXQ80211B', 0]}

面对象字典的数据结构

areaObjs = {'3000265': [1, '4', ['3000093', '3006515', '3007425', '3003723'], 'W500C40W'],

             . . . . . .
 
            '4019731': [1509, '5', ['4005397', '4005286', '4002422', '4005394', '4005340'], 'RIG120C35F'],

             . . . . . .
 
            '12000129': [2063, '4', ['12000229', '12000230', '12000161', '12000160'], 'RIG120C35F']}

框架均布荷载字典的数据结构

frameLoads = {1: ['5000001', 'DEAD', 'GLOBAL', 'Force', 'Gravity', 6000.0],

              . . . . . .
 
              555: ['8000460', 'DEAD', 'GLOBAL', 'Force', 'Gravity', 6000.0],

              . . . . . .
 
              1151: ['13001199', 'DEAD', 'GLOBAL', 'Force', 'Gravity', 6000.0]}

面均布荷载字典的数据结构

2.2 表格信息

在 s2k 文件中，模型的各种信息按类存储在不同的表格内，如轴网信息位于 TABLE: “GRID LINES” 表格内，通用材料属性信息位于 TABLE: “MATERIAL PROPERTIES 01 - GENERAL” 表格内，如下图所示。

s3fs 读取python python读取sas文件_Python_04

s2k文件

在更精确获取模型信息前，有必要知道该信息所属表格在 s2k 文件中的位置，这样有利于提高代码的运行效率，尽最大可能减少对 s2k 整个文件的遍历次数，获取 s2k 文件内各表格信息的 Python 代码如下。

tablenames, tablerows = [], []  # 表格名称及所在行号

with open(s2kpath, "r", encoding='utf-8') as f:
    s2kdata = f.readlines()  # 获取s2k文件内全部数据
    totallines = len(s2kdata)  # s2k文件总行数
    row = 0  # 行号计数器
    target_string = "TABLE"  # 目标字符串/待匹配字符串
    for line in s2kdata:
        row = row + 1

        ...
        ...
        ... 

for i in range(len(tablenames) - 1):
    tablename = tablenames[i]
    row1, row2 = tablerows[i]+1, tablerows[i + 1]-2
    tables[tablename] = [row1, row2]

row1, row2 = tablerows[-1]+1, totallines-2
tables[tablenames[-1]] = [row1, row2]
tablenames, tablerows = 0, 0  # 释放内存

print(tables)

由上述代码得到的表格信息字典 tables 如下图所示。字典 tables 的键名为每个表格的表头名称，如 TABLE: “PROGRAM CONTROL” 、TABLE: “GRID LINES” 等，与键名相对应的键值为该表格有效数据行在s2k文件中的位置信息，即始末行号，计数始于1。

如键名 TABLE: “PROGRAM CONTROL” 对应的键值为列表 [4，5]，其含义为表格 TABLE: “PROGRAM CONTROL” 的有效数据行为 s2k 文件中的第 4 行数据到第 5 行数据，通过此方式便可引用 s2k 文件中任意一表格内的数据。

s3fs 读取python python读取sas文件_Python_05

2.3 单位制

与 Abaqus、ANSYS 等有限元计算软件不同，在 SAP2000 中存在单位制，用户可在状态条最右侧的单位制下拉列表中选择模型所采用的单位制，用户切换当前单位制后，系统自动完成单位换算。

s3fs 读取python python读取sas文件_s3fs 读取python_06

单位制下拉列表

质量	长度	时间	力	压强/应力	能量	密度	弹性模量
kg	m	s	N	Pa	J	kg/m³	Pa
t	mm	s	N	MPa	N-mm	t/mm³	MPa

$s3fs 读取python python读取sas文件_数据_07$

为了规范有限元建模，同时也方便各软件间的数据交换，特此规定，无论原 SAP 模型采用何种单位制，通过 Python 处理后的模型信息均采用国际标准单位制：N-m-kg 。

通常情况下，在 SAP2000 中，工程师习惯采用 N-mm-t 单位制，而英制单位在我国不常用到，因此，本文所编写的 Python 程序，只识别 N-m-kg 和 N-mm-t 单位制，并实现 N-mm-t 向 N-m-kg 的自动换算，其他单位制，则视为非法单位，若被程序识别出，程序自动终止运行。

table = 'TABLE:  "PROGRAM CONTROL"'  # 表格名称
row1, row2 = tables[table]  # 确定有效数据行起始位置
tabledata = s2kdata[row1-1:row2]  # 表格有效数据行数据

target_strings = ["Version", "CurrUnits"]  # 待匹配字符串列表
result_strings = ["Null", "Null"]  # 匹配成功后输出的等号后的字符串/Null用于占位/防止报错

temp_line = []  # 多行数据合并为1行
for line in tabledata:
    line = line.strip()  # 去掉尾部换行符号
    line = re.split(r"  +", line)  # 2个及以上空格分割字符
    temp_line = temp_line + line

i = 0
for s in temp_line:
    target_string = target_strings[i]
    
    ...
    ...
    ...


version, currunits = result_strings
s2kname = s2kpath.split("\\")[-1]
s1 = "-->> File Name: %s." % s2kname
s2 = "-->> Total Lines: %s." % totallines
s3 = "-->> Program Name: SAP2000."
s4 = "-->> Version: %s." % version
s5 = "-->> Current Units: %s." % currunits
print("\n")
print(s1)
print(s2)
print(s3)
print(s4)
print(s5)
if currunits == '"N, mm, C"':
    sf = 0.001  # scale factor/缩放系数
    s6 = "-->> Scale Factor: %s." % sf
    print(s6)
    print("\n")
elif currunits == '"N, m, C"':
    sf = 1.0  # scale factor/缩放系数
    s6 = "-->> Scale Factor: %s." % sf
    print(s6)
    print("\n")
else:
    # pass/# 非法单位提示
    print("\n")
    s6 = "**** ERROR: Illegal units."
    s7 = 12 * " " + "Program termination."
    s8 = 12 * " " + "Please check the SAP model."
    print(s6)
    print(s7)
    print(s8)
    print("\n")
    exit()  # 非法单位程序自动终止运行

示例

运行上述代码，可获取 GG-JB.s2k 文件的基本信息及单位制信息，并确定其向国际标准单位的单位换算系数，若原 SAP 模型的单位制为 N-mm-kg，则单位换算系数 sf = 0.001，若原 SAP 模型的单位制为 N-m-kg，则单位换算系数 sf = 1.0，除此之外，其他的任何单位制都被视为非法单位，本程序不识别，程序运行至此处自动终止。

s3fs 读取python python读取sas文件_SAP2000_08

2.3 材料属性

def get_material_properties(s2kdata, tables, sf=1.0):
    """获取s2k文件内的材料属性，并将其存储到字典matProps中，返回值为材料属性字典matProps。
    :param s2kdata: s2k文件全部数据。
    :param tables: s2k文件表格信息字典。
    :param sf: 单位换算系数。
    """
    matProps = {}  # 材料属性字典
    table = 'TABLE:  "MATERIAL PROPERTIES 01 - GENERAL"'  # 表格名称
    row1, row2 = tables[table]  # 确定有效数据行起始位置
    tabledata = s2kdata[row1 - 1:row2]  # 表格有效数据行数据
    ...
    ...
    ... 
    matProps[matname] = [matnum, mattype, symtype, "Null", "Null", "Null", "Null"]
    # 材料编号、材料种类、力学分类、质量密度、弹性模量、剪切模量、泊松比。
    table = 'TABLE:  "MATERIAL PROPERTIES 02 - BASIC MECHANICAL PROPERTIES"'  # 表格名称
    row1, row2 = tables[table]  # 确定有效数据行起始位置
    tabledata = s2kdata[row1 - 1:row2]  # 表格有效数据行数据
    ...
    ...
    ... 
    matProps[matname][2:] = rho, e, g, nu
    return matProps

示例

运行上述代码，可获取 GG-JB.s2k 文件内所定义的材料属性信息，并将其存储在 matProps 字典中，各种材性的单位采用国际标准单位，即 N-m-kg。

字典 matPropss 的键名为各材料的材料名称，如 C30B、C40B 等，与键名相对应的键值为该材料的基本属性信息，依次为：材料编号、材料种类 (混凝土、钢材、钢筋等)、材料力学分类 (各向同性、单轴受力等)、质量密度、弹性模量、剪切模量、泊松比。

s3fs 读取python python读取sas文件_s3fs 读取python_09

2.4 框架截面属性

def get_frame_section_properties(s2kdata, tables, sf=1.0):
    """获取s2k文件内的框架截面属性，并将其存储到字典frameSecProps中，返回值为框架截面属性字典frameSecProps。
    :param s2kdata: s2k文件全部数据。
    :param tables: s2k文件表格信息字典。
    :param sf: 单位换算系数。
    """
    frameSecProps = {}  # 框架截面属性字典
    table = 'TABLE:  "FRAME SECTION PROPERTIES 01 - GENERAL"'  # 表格名称
    row1, row2 = tables[table]  # 确定有效数据行起始位置
    tabledata = s2kdata[row1 - 1:row2]  # 表格有效数据行数据
    secnum = 0  # 截面编号计数器
    ...
    ...
    ... 
    target_strings = ["Material", "Shape", "t3", "t2", "tf", "tw", "t2b", "tfb", "Area",
                      "TorsConst", "I33", "I22", "I23", "AS2", "AS3", "S33", "S22", "Z33",
                      "Z22", "R33", "R22"]  # 待匹配字符串
    ...
    ...
    ... 
    return frameSecProps

示例

运行上述代码，可获取 GG-JB.s2k 文件内所定义的框架截面属性信息，并将其存储在 frameSecProps 字典中，各种几何属性的单位采用国际标准单位，即 N-m-kg。

字典 frameSecProps 的键名为各框架截面的截面名称，如 B1000X500JXC35B 等，与键名相对应的键值为该框架截面的基本属性信息，依次为：框架截面编号、截面采用的材料名称、截面形状、垂直 3 轴方向上的截面尺寸、垂直 2 轴方向上的截面尺寸、上翼缘厚度、腹板厚度、下翼缘宽度、下翼缘厚度、横截面积、扭转常数、截面对局部 3 轴的惯性矩、截面对局部 2 轴的惯性矩、截面惯性积、2 方向上的剪切面积、3 方向上的剪切面积、截面对局部 3 轴的截面模量、截面对局部 2 轴的截面模量、截面对局部 3 轴的塑性截面模量、截面对局部 2 轴的塑性截面模量、截面对局部 3 轴的回转半径、截面对局部 2 轴的回转半径。

s3fs 读取python python读取sas文件_s2k_10

2.5 面截面属性

def get_area_section_properties(s2kdata, tables, sf=1.0):
    """获取s2k文件内的面截面属性，并将其存储到字典areaSecProps中，返回值为面截面属性字典areaSecProps。
    :param s2kdata: s2k文件全部数据。
    :param tables: s2k文件表格信息字典。
    :param sf: 单位换算系数。
    """
    areaSecProps = {}  # 框架截面属性字典
    table = 'TABLE:  "AREA SECTION PROPERTIES"'  # 表格名称
    row1, row2 = tables[table]  # 确定有效数据行起始位置
    tabledata = s2kdata[row1 - 1:row2]  # 表格有效数据行数据
    secnum = 0  # 截面编号计数器
    ...
    ...
    ... 
    target_strings = ["Material", "AreaType", "Type", "Thickness"]  # 待匹配字符串
    sfs = [None, None, None, sf]  # 单位换算系数
    ...
    ...
    ... 
    result = round(float(result) * sf, 6)  # 单位换算
    areaSecProps[secname][i + 1] = result
    return areaSecProps

示例

运行上述代码，可获取 GG-JB.s2k 文件内所定义的面截面属性信息，并将其存储在 areaSecProps 字典中，各种几何属性的单位采用国际标准单位，即 N-m-kg。

字典 areaSecProps 的键名为各面截面的截面名称，如 RIG1000C35F、RIG100C35F 等，与键名相对应的键值为该面截面的基本属性信息，依次为：面截面编号、面截面类型 (壳、平面、轴对称实体) 、计算类型 (厚壳、薄壳、厚板、薄板、膜等)、面截面的厚度。

s3fs 读取python python读取sas文件_SAP2000_11

2.6 节点对象

def get_joint_objects(s2kdata, tables, sf=1.0):
    """获取s2k文件内的节点对象相关信息，并将其存储到字典jointObjs中，返回值为结点对象信息字典jointObjs。
    :param s2kdata: s2k文件全部数据。
    :param tables: s2k文件表格信息字典。
    :param sf: 单位换算系数。
    """
    jointObjs = {}  # 框架截面属性字典
    table = 'TABLE:  "JOINT COORDINATES"'  # 表格名称
    row1, row2 = tables[table]  # 确定有效数据行起始位置
    tabledata = s2kdata[row1 - 1:row2]  # 表格有效数据行数据
    jointnum = 0  # 结点编号计数器
    ...
    ...
    ... 
    num = len(jointObjs)
    s = "Joints Number: %s." % num
    print(s)
    return jointObjs

示例

运行上述代码，可获取 GG-JB.s2k 文件内所定义的各节点对象信息，并将其存储在 jointObjs 字典中，各节点坐标值的单位采用国际标准单位 m 。

字典 jointObjs 的键名为各节点对象的标签名称，如 2001902、2001903 等，与键名相对应的键值为该节点对象的基本信息，依次为：结点编号、结点所采用的坐标系名称、坐标系类型 (笛卡尔、柱坐标、球坐标) 、全局坐标值 x 、全局坐标值 y 、全局坐标值 z 。

s3fs 读取python python读取sas文件_SAP2000_12

2.7 框架对象

def get_frame_objects(s2kdata, tables):
    """获取s2k文件内的框架对象相关信息，并将其存储到字典frameObjs中，返回值为框架对象信息字典frameObjs。
    :param s2kdata: s2k文件全部数据。
    :param tables: s2k文件表格信息字典。
    """
    frameObjs = {}  # 框架截面属性字典
    table = 'TABLE:  "CONNECTIVITY - FRAME"'  # 表格名称
    row1, row2 = tables[table]  # 确定有效数据行起始位置
    tabledata = s2kdata[row1 - 1:row2]  # 表格有效数据行数据
    framenum = 0  # 框架编号计数器
    ...
    ...
    ... 
    s = "Frames Number: %s." % framenum
    print(s)
    return frameObjs

示例

运行上述代码，可获取 GG-JB.s2k 文件内所定义的各框架对象信息，并将其存储在 frameObjs 字典中，框架对象信息的相关单位采用国际标准单位，即 N-m-kg 。

字典 frameObjs 的键名为各框架对象的标签名称，如 3000092、3000093 等，与键名相对应的键值为该框架对象的基本信息，依次为：框架编号 (始于1/ANSYS)、框架起点标签、框架终点标签、框架所采用的截面名称、局部轴转角。

s3fs 读取python python读取sas文件_s3fs 读取python_13

2.8 面对象

def get_area_objects(s2kdata, tables):
    """获取s2k文件内的面对象相关信息，并将其存储到字典areaObjs中，返回值为面对象信息字典areaObjs。
    :param s2kdata: s2k文件全部数据。
    :param tables: s2k文件表格信息字典。
    """
    areaObjs = {}  # 框架截面属性字典
    table = 'TABLE:  "CONNECTIVITY - AREA"'  # 表格名称
    row1, row2 = tables[table]  # 确定有效数据行起始位置
    tabledata = s2kdata[row1 - 1:row2]  # 表格有效数据行数据
    areanum = 0  # 面对象编号计数器
    for line in tabledata:
        line = line.strip()  # 去掉尾部换行符号
        line = re.split(r"  +", line)  # 2个及以上空格分割字符
        joint_labels, secname = [], "Null"  # 组成面对象的节点标签、面对象所采用的截面名称
        label = re.split(r"=+", line[0])[-1]  # 至少一个等号分割字符/获得面对象标签
        match = re.match("NumJoints", line[1])  # 匹配面对象包含的节点个数
    ...
    ...
    ... 
    return areaObjs

示例

运行上述代码，可获取 GG-JB.s2k 文件内所定义的各面对象信息，并将其存储在 areaObjs 字典中，面对象信息的相关单位采用国际标准单位，即 N-m-kg 。

字典 areaObjs 的键名为各面对象的标签名称，如 12000123、12000129 等，与键名相对应的键值为该框架对象的基本信息，依次为：面的编号 (始于1/ANSYS)、面对象所包含的节点个数、各点标签、面对象所采用的截面名称。

s3fs 读取python python读取sas文件_s3fs 读取python_14

2.9 框架均布荷载

在 SAP2000 中，对于框架填充墙或者设备荷载的处理，一般将其自重转化为线荷载指定给梁。

s3fs 读取python python读取sas文件_SAP2000_15

def get_frame_loads(s2kdata, tables, sf):
    """获取s2k文件内的框架对象荷载信息，并将其存储到字典frameLoads中，返回值为框架荷载信息字典frameLoads。
    :param s2kdata: s2k文件全部数据。
    :param tables: s2k文件表格信息字典。
    :param sf: 单位换算系数。
    """
    frameLoads = {}  # 框架对象荷载信息字典
    table = 'TABLE:  "FRAME LOADS - DISTRIBUTED"'  # 表格名称
    row1, row2 = tables[table]  # 确定有效数据行起始位置
    tabledata = s2kdata[row1 - 1:row2]  # 表格有效数据行数据
    loadnum = 0  # 荷载编号计数器
    for line in tabledata:
    ...
    ...
    ... 
    s = "Frame Loads Number: %s." % loadnum
    print(s)
    return frameLoads

示例

运行上述代码，可获取 GG-JB.s2k 文件内所定义的作用在框架对象上的均布荷载信息，并将其存储在 frameLoads 字典中，荷载信息的相关单位采用国际标准单位，即 N-m-kg 。

字典 frameLoads 的键名为框架均布荷载编号 (始于1/ANSYS)，如 1、2 等，与键名相对应的键值为作用在该框架上的均布荷载信息，依次为：框架标签、荷载模式 (恒活、活荷等)、采用的坐标系、荷载类型 (力/力矩)、荷载作用方向、作用在框架上的均布荷载数值。

因为作用在框架上的荷载不止一种，方向也不止一个，因此，很难用框架标签作为框架均布荷载字典的键名，来对作用在框架上的均布荷载进行区分。另外，由此段代码得到的是作用在框架对象全跨上的均布荷载，全跨梯形分布荷载自动进行换算，分段梯形分布荷载自动被过滤掉，不在考虑范围内。

s3fs 读取python python读取sas文件_Python_16

2.10 面均布荷载

在 SAP2000 中，面对象的分布荷载主要有两种施加方式：面均布荷载和面荷载均匀导荷到框架。后者自动将均匀面荷载转化为等效线荷载作用于面周边的框架单元，从而实现面对象的自动导荷。

s3fs 读取python python读取sas文件_SAP2000_17

面对象均布荷载的指定

s3fs 读取python python读取sas文件_SAP2000_18

def get_area_loads(s2kdata, tables, sf):
    """获取s2k文件内的面对象荷载信息，并将其存储到字典areaLoads中，返回值为面荷载信息字典areaLoads。
    :param s2kdata: s2k文件全部数据。
    :param tables: s2k文件表格信息字典。
    :param sf: 单位换算系数。
    """
    areaLoads = {}  # 框架对象荷载信息字典
    table1 = 'TABLE:  "AREA LOADS - UNIFORM"'  # 表格名称
    table2 = 'TABLE:  "AREA LOADS - UNIFORM TO FRAME"'  # 表格名称
    row1, row2 = tables[table1]  # 确定有效数据行起始位置
    tabledata1 = s2kdata[row1 - 1:row2]  # 表格有效数据行数据
    row1, row2 = tables[table2]  # 确定有效数据行起始位置
    tabledata2 = s2kdata[row1 - 1:row2]  # 表格有效数据行数据
    tabledata = tabledata1 + tabledata2  # 两个表格数据合并
    loadnum = 0  # 荷载编号计数器
    ...
    ...
    ... 
    s = "Area Loads Number: %s." % loadnum
    print(s)
    return areaLoads

2.11 边界条件

2.10 数据的精加工

2.11 创建csv文件

三、创建CAD图形

四、创建Abaqus模型

五、创建ANSYS模型

六、尾声

本文章为转载内容，我们尊重原作者对文章享有的著作权。如有内容错误或侵权问题，欢迎原作者联系我们进行内容更正或删除文章。

上一篇：python3 is needed by python3 is needed by meson

下一篇：java锁重偏向 java 重入锁

提问和评论都可以，用心的回复会被更多人看到评论

发布评论

相关文章

官方博客	全部文章	热门标签	班级博客
了解我们	网站地图	意见反馈

鸿蒙开发者社区	51CTO学堂
51CTO	软考资讯