矩阵含义
#!/usr/bin/env python# -*- coding: utf-8 -*-
# 导入必要的模块
import numpy as np
# 函数，计算弹簧单元刚度矩阵
def ele_K(K):
Ke = np.array([,[-K,K]])
return Ke
print(ele_K(4))
dimen = 1 # 一维
K1 = 4.0
K2 = 6.0
K3 = 3.0
# 节点坐标阵
node_coord = np.array()
# 单元节点阵
ele_node = np.array([,,,])
# 外力向量
F = np.array()
K = np.array()
# 总的自由度数
ndof = dimen*len(node_coord)
# 初始化总体刚度阵
K_total = np.zeros((ndof,ndof))
print(K_total)
# 计算总体刚度阵
for e in range(0,len(ele_node)):
n1 = ele_node-1
n2 = ele_node-1
K_total,)] += ele_K(K)
print(K_total)
这个代码里面，最后一句K_total,)] += ele_K(K)，里面的,,这个含义是什么？ 子列表，就像C语言里的二维数组(又称矩阵)，只是列表的线性扩展，你想想二维这个词就明白了，父列表是列，子列表是每一行
页:
[1]