python的拱桥振型简化计算

2024-02-25 12:20

本文主要是介绍python的拱桥振型简化计算,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

from sympy import *
import numpy as np
import matplotlib.pyplot as plt
import math
#基本参数
g=9.8
l=60
f=2
Ac=8.15
Ic=3622000*(10e-8)
E=2.06*10e3
m=1435.2
#方程运算
H0=(m*g*(l)**2)/(8*f)
He=20502900
v1=1-(H0/He)
v2=1-H0/(4*He)
v3=1-H0/(9*He)
p=l**2/(8*f)
L=l*(1+(16/3)*(f/l)**2)
a21=8*E*Ac*l
a22=(m*((np.pi)**2)*(p)**2)*L
a2=(a21/a22)*10e4
c2=(2*(np.pi)/l)**2*sqrt((E*Ic*v2)/m)
c2=c2*1000#平衡系数
c1=((np.pi/l)**2)*sqrt((E*Ic*v1)/m)
c1=c1*1000#平衡系数
c3=((3*np.pi/l)**2)*sqrt((E*Ic*v1)/m)
c3=c3*1000
w=symbols('w')
G=solve([w**4-(c1**2+c3**2+(10/9)*a2)*w**2+(c1**2*c3**2+a2*(c1**2/9+c3**2))],[w])print(v1,v2,v3,a21,a22,a2)
print(c2,c1,c3)
print(G,c2)

以抛物线型拱为例,采用传统的数值计算方法过于繁琐(此题以李国豪版桥梁结构稳定与振动为例)

考虑拱桥的前几阶振型必须先考虑前几阶振型所对应的频率


(以三阶为例,一个为基频另外两个为对称频率)

对于此联立方程不只能求出振型比率

此外思考拱桥抛物线的矢高与拱桥长度的关系



采用最小2乘法拟合分析

根据上式代码计算出的结果


import matplotlib.pyplot as plt
import numpy as np
from scipy.optimize import leastsq
import math
plt.rcParams['font.sans-serif']=['SimHei']#显示中文
plt.rcParams['axes.unicode_minus']=False
X=np.array([2,2.5,3,3.5,4,4.5,5,5.5,6,6.5,7,7.5,8,8.5,9,9.5,10])
Y1=np.array([24.52,24.61,24.68,24.73,24.76,24.788,24.81,24.827,24.842,24.855,24.865,24.875,24.883,24.89,24.896,24.902,24.907])
Y2=np.array([14.75,17.9,21.05,24.17,27.23,30.214,33.086,35.81,38.352,40.662,42.7,44.44,45.877,47.035,47.955,48.683,49.26])def func1(params,x):A,B=paramsreturn A*np.sqrt(1-B*(1/x))def error1(params,x,y):return func1(params,x) - ydef  slovePara1():p0=[1,1]Para = leastsq(error1,p0,args=(X,Y1))return Paradef func2(params,x):k1,k2,k3,k4,k5,k6,k7,k8,k9,k10,k11,k12=paramsreturn 0.235702260395516*np.sqrt((k1+k2*x+k3*x**2+k4*x**3+(k5+k6*x+k7*x**2+k8*x**3+k9*x**4+k10*x**5+k11*x**6)**0.5)/(x*(k12*x**2+1)))def error2(params,x,y):return func2(params,x) - ydef slovePara2():p0=[1,1,1,1,1,1,1,1,1,1,1,1]Para = leastsq(error2,p0,args=(X,Y2))return Paradef solution():Para = slovePara1()A, B,  =Para[0]plt.scatter(X, Y1, color='green',label=u'已知样点',linewidth=0.3)x=np.linspace(2,10,50)y1=A*np.sqrt(1-B*(1/x))plt.plot(x,y1,color='red',label='w2',linewidth=2)plt.legend()Para = slovePara2()k1, k2, k3,k4,k5,k6,k7,k8,k9,k10,k11,k12  =Para[0]x=np.linspace(2,10,50)y2=0.235702260395516*np.sqrt((k1+k2*x+k3*x**2+k4*x**3+(k5+k6*x+k7*x**2+k8*x**3+k9*x**4+k10*x**5+k11*x**6)**0.5)/(x*(k12*x**2+1)))plt.plot(x,y2,color='blue',label='w1',linewidth=2)plt.legend()plt.scatter(X, Y2, color='green',linewidth=0.3)plt.title(u'w2与w1随f的变化情况')plt.show()solution()

所画图形如下



这篇关于python的拱桥振型简化计算的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!


原文地址:https://blog.csdn.net/cywtiancai/article/details/79605872
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.chinasem.cn/article/745451

相关文章

Python调用Orator ORM进行数据库操作

《Python调用OratorORM进行数据库操作》OratorORM是一个功能丰富且灵活的PythonORM库,旨在简化数据库操作,它支持多种数据库并提供了简洁且直观的API,下面我们就... 目录Orator ORM 主要特点安装使用示例总结Orator ORM 是一个功能丰富且灵活的 python O

Python使用国内镜像加速pip安装的方法讲解

《Python使用国内镜像加速pip安装的方法讲解》在Python开发中,pip是一个非常重要的工具,用于安装和管理Python的第三方库,然而,在国内使用pip安装依赖时,往往会因为网络问题而导致速... 目录一、pip 工具简介1. 什么是 pip?2. 什么是 -i 参数?二、国内镜像源的选择三、如何

python使用fastapi实现多语言国际化的操作指南

《python使用fastapi实现多语言国际化的操作指南》本文介绍了使用Python和FastAPI实现多语言国际化的操作指南,包括多语言架构技术栈、翻译管理、前端本地化、语言切换机制以及常见陷阱和... 目录多语言国际化实现指南项目多语言架构技术栈目录结构翻译工作流1. 翻译数据存储2. 翻译生成脚本

如何通过Python实现一个消息队列

《如何通过Python实现一个消息队列》这篇文章主要为大家详细介绍了如何通过Python实现一个简单的消息队列,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录如何通过 python 实现消息队列如何把 http 请求放在队列中执行1. 使用 queue.Queue 和 reque

Python如何实现PDF隐私信息检测

《Python如何实现PDF隐私信息检测》随着越来越多的个人信息以电子形式存储和传输,确保这些信息的安全至关重要,本文将介绍如何使用Python检测PDF文件中的隐私信息,需要的可以参考下... 目录项目背景技术栈代码解析功能说明运行结php果在当今,数据隐私保护变得尤为重要。随着越来越多的个人信息以电子形

使用Python快速实现链接转word文档

《使用Python快速实现链接转word文档》这篇文章主要为大家详细介绍了如何使用Python快速实现链接转word文档功能,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 演示代码展示from newspaper import Articlefrom docx import

Python Jupyter Notebook导包报错问题及解决

《PythonJupyterNotebook导包报错问题及解决》在conda环境中安装包后,JupyterNotebook导入时出现ImportError,可能是由于包版本不对应或版本太高,解决方... 目录问题解决方法重新安装Jupyter NoteBook 更改Kernel总结问题在conda上安装了

Python如何计算两个不同类型列表的相似度

《Python如何计算两个不同类型列表的相似度》在编程中,经常需要比较两个列表的相似度,尤其是当这两个列表包含不同类型的元素时,下面小编就来讲讲如何使用Python计算两个不同类型列表的相似度吧... 目录摘要引言数字类型相似度欧几里得距离曼哈顿距离字符串类型相似度Levenshtein距离Jaccard相

Python安装时常见报错以及解决方案

《Python安装时常见报错以及解决方案》:本文主要介绍在安装Python、配置环境变量、使用pip以及运行Python脚本时常见的错误及其解决方案,文中介绍的非常详细,需要的朋友可以参考下... 目录一、安装 python 时常见报错及解决方案(一)安装包下载失败(二)权限不足二、配置环境变量时常见报错及

Python中顺序结构和循环结构示例代码

《Python中顺序结构和循环结构示例代码》:本文主要介绍Python中的条件语句和循环语句,条件语句用于根据条件执行不同的代码块,循环语句用于重复执行一段代码,文章还详细说明了range函数的使... 目录一、条件语句(1)条件语句的定义(2)条件语句的语法(a)单分支 if(b)双分支 if-else(