MPEMathematical Problems in Engineering1563-51471024-123XHindawi10.1155/2020/62459436245943Research ArticleApplication of Multiphysical Domain Modeling in Electrical System SimulationZhengCao1YuanjunZhou1https://orcid.org/0000-0003-2559-2914NaWang2StranoSalvatore1School of Automation Science and Electrical EngineeringBeihang UniversityBeijing 100191Chinabuaa.edu.cn2Engineering Training CenterBeihang UniversityBeijing 100191Chinabuaa.edu.cn20202892020202027122019010620202408202028920202020Copyright © 2020 Cao Zheng et al.This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The energy structure change of more electric aircraft makes the aircraft airborne system more complicated. Each subsystem realizes the transmission, interaction, and conversion of energy and information through the dynamic coupling and coordinated control of electrical energy, mechanical energy, hydraulic energy, and thermal energy. This paper applies the multiphysical domain modeling method with the parameter identification according to the original model data. Based on the power conversion relationship of electrical equipment, it defines the port with the power potential variable and flow variable and is supplemented by the information control. It can show the dynamic characteristics, power conversion, and loss characteristics of the device itself. The models can conveniently perform the large-scale system integration, which not only can build the complex electrical equipment formed by multiphysical domain models with the series connection but also can build a complex power supply system formed by multiphysical domain models with the parallel connection.

1. Introduction

The secondary energy of the more electric aircraft gradually uses electrical energy to replace hydraulic energy and pneumatic energy for the aircraft energy optimization . In the design of the aircraft electrical system, the performance of the electrical equipment in terms of energy utilization, power loss, and heat conduction needs to be taken as one of the evaluation indicators .

The electrical equipment operation of the more electric aircraft involves the transmission, interaction, and conversion of the multiphysical domain, such as electrical energy, mechanical energy, hydraulic energy, thermal energy, and information [7, 8]. Since the conversion process is dynamic, it is obviously inaccurate to discuss the energy utilization performance simply by using device efficiency, and a digital simulation method is needed [9, 10]. The energy usage of the aircraft electrical system is related to flight missions, involving energy conversion, energy transmission, and energy consumption. Therefore, it is necessary to form the actual electrical system in the structure and perform dynamic simulation under different missions .

In order to realize the simulation of the electrical system multiphysical domain power conversion, the digital model of the electrical equipment is required to be established. The multiphysical domain modeling can solve the coupling problems . The traditional multiphysical modeling methods are the physical structure modeling, the equivalent circuit modeling, and the analysis structural modeling. The structure of electrical equipment is more accurately represented to describe the characteristics of electrical equipment by the physical structure modeling. However, the problem is the operation speed. When the more complex system including nonlinear components is built, it will cause the computer calculation speed to be too slow in simulation.

The equivalent circuit modeling has the disadvantage because the structure is simplified. The more details and nonlinear characteristics cannot be accurately represented. The modeling workload by the analysis structural modeling is relatively large, and the nonlinear modeling is especially cumbersome. The input and output may be variables of different physical domains. When the electrical system integration is achieved, the system structure may be confused and inconvenient in engineering application.

The bond graph method is applied for the hydraulic system [13, 14], and the Modelica is applied for the electromechanical systems [15, 16]. It can accurately represent its power conversion and loss characteristics when expressing the device dynamic operating characteristics . In addition, because the electrical component model is used to build an electrical system, which contains a large number of nonlinear components, the transients in the simulation result in slower operation speed. Therefore, in order to achieve large-scale integrated simulation, there should be a more suitable computing speed.

This paper applies the multiphysical domain modeling method with the parameter identification according to the original model data. The analytical model based on the energy conversion process of the electrical equipment is supplemented by the information control. The model uses the energy transmission between the electrical equipment as the interface between the models, enabling the large-scale integrated electrical system simulation, and the simulation operation speed is improved.

2. Electrical Equipment Structure and Multiphysical Domain Model

Compared with the traditional aircraft, the increase in electrical equipment of more electric aircraft is the electrical drive equipment. The electromechanical actuator (EMA) [18, 19] and the electrohydraulic actuator (EHA) [20, 21] of the flight control system are the most typical multiphysical domain energy conversion devices. Because it presents a high-power short-term working mode, it has the greatest impact on the power characteristics of the power supply system. This paper analyzes the modeling method of the typical electrical equipment.

2.1. Structure of Typical Electrical Equipment

The electromechanical actuator is taken as an example to analyze the structural characteristics of the electrical equipment. The EMA is an electrical drive position servo control system for the aircraft surface. The physical model of the EMA is shown in Figure 1. It consists of the inverter, the brushless direct current motor (BLDCM), the reduction gear, and the microprocessor. The LC electromagnetic interference filter (EMI) is also on the power supply side.

Principle block diagram of the electromechanical actuator.

The electrical drive equipment shown in Figure 1 is structurally divided into two parts.

2.1.1. Power Conversion Section

From the perspective of energy conversion, the electromechanical actuator performs the conversion of the electrical energy to mechanical energy, and the conversion power is determined by the fight mission. At the same time, all parts of the equipment have the power loss. It transforms parts of the electrical energy and mechanical energy into thermal energy and forms a secondary power conversion process.

2.1.2. Information Conversion Section

The main function of the information conversion part is to complete the control during the equipment. For example, the EMA realizes the closed-loop control of the rudder surface position θa, the motor speed Ωm, and the armature current Id. In addition, it can also monitor the running status of the device and communicate with the host computer and other devices.

The physical model shown in Figure 1 is not convenient for the large-scale integration. Because of the nonlinear components, simulation transients result in the slower operation.

2.2. Model Based on Power Conversion

According to the structure of the electrical device shown in Figure 1, the model of the electrical device can be divided into two modules: the power conversion module (PCM) and the information conversion module (SCM), as shown in Figure 2:

Model structure of the electrical equipment.

The PCM represents the device power conversion process. The power conversion is mainly the conversion between two kinds of energy sources. The thermal energy generated by the power loss is the third additional energy conversion relationship for the performance analysis of the power loss and efficiency.

The SCM function of the electrical equipment is the same as that of the information conversion part in Figure 1. The equipment operational control becomes the PCM control, and monitoring and communication status can be realized at the same time.

3. Power Conversion Module of Multiphysical Domain Model3.1. Power Conversion of Multiphysical Domain

The PCM shown in Figure 2 exhibits the power conversion characteristics implemented by the electrical device, which is the power that transforms one type of physical domain power into another. For example, the generator realizes the energy conversion from mechanical energy to electrical energy, and the electrical equipment realizes the energy conversion from electrical energy to various forms of energy, including mechanical energy, hydraulic energy, pneumatic energy, and thermal energy. Therefore, the operation of the more electric aircraft electrical system is a multiphysical domain energy conversion process.

When the power conversion model shown in Figure 2 is established, the input power and output power are represented by the potential variable Y and the flow variable X, which are defined by the connection type of the multiphysical domain modeling . The power types involved in the electrical system of the more electric aircraft and the corresponding potential variables and flow variables are shown in Table 1.

Potential variables and flow variables for different power types.

 Power form Potential variable Y Flow variable X Electrical power Voltage U (V) Current I (A) Rotating machine power Speed Ω (rad/s) Torque T (N⋅m) Linear machine power Force F (N) Moving speed v (m/s) Hydraulic power Pressure p (MPa) Flow q (m3/s) Thermal power Temperature T (K) Flow (heat flow) Φ (W) Electrical power Voltage U (V) Current I (A)
3.2. Analysis Model of PCM

The PCM of the electrical equipment model can be represented as Figure 3(a). The potential variable of the input power (source power) is represented as Ys, and the flow variable is represented as Xs. The potential variable of the output power (load power) is YL, and the flow variable is XL. The power loss output is generally thermal power. The ambient temperature Tw is a potential variable, and the heat flow Φw formed by the loss is a pseudo flow variable. C is the control signal.

Input/output variables of the power transformation model. (a) Variable relationship of electrical equipment PCM. (b) Variable relationship of electromechanical actuator PCM.

According to the relationship of the PCM input/output variables shown in Figure 3(a), the analytical method can be described as(1)YLsXss=W1sW3sW2sW4sYssXLs,where the transfer functions of W1 (s), W2 (s), W3 (s), and W4 (s) reflect the dynamic characteristics between variables in the power conversion process. The power loss of the electrical device is determined by the device operating state. Therefore, the loss model can be expressed as a function of other variables, namely,(2)Φw=fwYs,Xs,YL,XL,Tw.

For the electromechanical actuator shown in Figure 1, the analytical formula of the PCM shown in Figure 3(b) can be described as follows:(3)ΩasIss=W1sW3sW2sW4sUssTLs.

The potential variable of the input power in equation (3) is the power supply voltage Us, and the flow variable is the supply current Is. The potential variable of the output power is the angular velocity Ωa of the aircraft rudder surface, and the flow variable is the resistance torque TL of the rudder surface. D is the duty cycle of the PWM signal.

3.3. System Control Model of PCM

The energy conversion relationship of the entire device is shown in equation (1). As shown in Figure 2, the energy conversion process of the electrical device needs to perform under the SCM control. In order to express the power conversion relationship under the control process, the power conversion module can be rewritten as(4)YLsXss=W1sW3sW2sW4sYssXLs=G1sG3sG2sG4syscsxLcs,where ysc and xLc are the device control and interference variables, respectively, which are generated by the SCM and have a certain functional relationship with Ys and XL. G1 (s), G2 (s), G3 (s), and G4 (s) are the transfer functions of the device control system.

The electromechanical actuator shown in Figure 1 is taken as an example and the power conversion module of the control system is(5)ΩasIss=W1sW3sW2sW4sUssTLs=G1sG3sG2sG4suscsTLcs,where the voltage applied to the motor armature winding under the SCM control is usc, and the relationship with the power source voltage Us is usc = DUs. TLc is the load torque corresponding to the electromagnetic torque. In addition to the resistance torque TL on the rudder surface, it will also include the friction torque and damping torque on the motor and reducer.

3.4. Power Loss Model of PCM

The power loss model of equation (2) is generally described according to the power consumption principle of the device. For an electrical device without energy storage characteristics, the input power Ps and the output power PL can also be approximated and described as follows:(6)Φw=PsPL=YsXsYLXL.

The electromechanical actuator shown in Figure 1 is also taken as an example. The power loss is mainly the inverter loss ∆Pec, the copper loss ∆Pcu of the motor, the iron loss ∆Pfe of the motor, the mechanical loss ∆Pm caused by the friction of the motor and the reducer, and the power loss of the electronic device which can be approximated to the motor copper loss ∆Pcu and classified into the copper loss. Therefore, the power loss can be described as(7)Φw=ΔPΣ=ΔPecΔPcu+ΔPfe+ΔPm=kecId+kcuId2+kfeΩm+kmΩm=kecId+kcuId2+kfe+kmΩm,where Id is the armature current. Due to the motor moment inertia, the electrical drive device is expressed as the device containing the energy storage component. If equation (6) is used to analyze the device loss, its dynamic characteristics will have some error.

4. Information Conversion Module of Multiphysical Domain Model

The SCM of the multiphysical domain model is used to simulate the information conversion of the electrical equipment. It can realize the operation control and the state measurement of the system, such as the functions of the controller and the sensor in Figure 1.

4.1. Control Information Conversion of SCM

For the multiphysical domain model described in Figure 2, the first function of the SCM is to simulate the electrical device control function, which forms the control signal ysc. The second function is to generate the disturbance signal xLc according to the operating state. It can be described as(8)ysct=fscYs,Xs,YL,XL,R,txLct=fLcYs,Xs,YL,XL,R,t,where R is the reference signal given by the host computer and other devices.

The electromechanical actuator shown in Figure 1 is still taken as an example. Since it is a position servo control system, the device needs to realize three closed-loop control of the rudder surface position θa, the motor speed Ωm, and the motor armature current Id. The control voltage signal should be(9)usc=DUs=WacrIdIdUs=WacrWasrΩmΩmIdUs=WacrWasrWaprθaθaΩmIdUs,where Wacr, Wasr, and Wapr are the transfer functions of current, speed, and position regulator, Id is the reference value of the armature current, Ωm is the reference value of motor speed, and θa is the reference value of the rudder surface position. The resistance torque applied to the aircraft rudder surface can be expressed by the hinge moment, and kθ is the hinge moment coefficient of the rudder surface, which is TL = kθθa. The load torque on the actuator is(10)TLc=Tf+kΩΩm+kθθa,where Tf is the friction torque of the motor and the reducer and kΩ is the speed damping coefficient of the motor.

4.2. Test of Running Status

The function of the SCM operation state measurement is to simulate the sensor function of the electrical equipment. The test signal can be used as a feedback signal of the system control and provide to the monitoring device or the host computer, which is the output information in Figure 2.

The SCM test information is defined as S, which can be expressed as(11)St=fsYs,Xs,YL,XL,R,t.

For the EMA control and disturbance signals in equations (9) and (10), the required test signals are(12)Id=ksIsΩm=kiΩaθa=0tΩadt,where Id is the motor armature current, Is is the power supply current, ks is the proportional coefficient, Ωa is the angular velocity of the aircraft rudder surface, Ωm is the motor speed, Ki is the reduction ratio of the reduction gear, and θa is the rudder surface position.

If it needs to know other signal information, such as the motor electromagnetic torque Te, it can be described by(13)Te=ktId=ktksIs,where kt is the torque factor.

5. Modeling and Verification of Electromechanical Actuator

The electromechanical actuator is used as multiphysical domain modeling object, and the consistency is verified by comparison with the original model simulation results.

5.1. Electromechanical Actuator Parameters

The multiphysical domain modeling object is the EMA model shown in Figure 1. It is composed of nonlinear characteristic components such as the inverter and the brushless DC motor, which is called the EMA original model. For example, the main parameters of an EMA are the power supply voltage 270VDC, the motor rated angular speed 600 rad/s, the maximum rudder surface torque of 11500 N·m, the electromagnetic time constant 0.000865s, and the electromechanical time constant 0.00652s.

5.2. Electromechanical Actuator Modeling

The multiphysical domain modeling is mainly to establish the PCM transfer functions of G1 (s), G2 (s), G3 (s), G4 (s) with equation (2) and the power loss with equation (4). According to the original model of the electromechanical actuator shown in Figure 1, the transfer function G1 (s)∼G4 (s) with equation (5) can be established by the principle analysis method or the signal identification method.

Since the EMA in Figure 1 contains nonlinear components, the principle analysis modeling is complex; the signal identification method is adopted. The original model shown in Figure 1 is set as the open-loop structure. A step signal is applied to the control voltage usc, and it is simulated with TL = 0 to obtain the unloaded angular velocity Ωa0 of the rudder surface and the power supply current Is0. Then, the load torque TL is set to a step signal for simulation, and the speed signal Ωa and the power source current ∆Is are obtained. The obtained signals are, respectively, identified and the PCM is described as follows:(14)G1s=Ωm0suscs=3.137.33×107s2+2.303×103s+1G2s=Is0suscs=6.69×103s+0.00058.781×107s2+2.404×103s+1G3s=ΔΩmsTLcs=1.24×103s3.727.618×107s2+2.36×103s+1G4s=ΔIssTLcs=3.688.578×107s2+2.398×103s+1.

Combining the power conversion module of equation (14) with the information conversion module of equations (9) and (10), the multiphysical domain model shown in Figure 2 can be constructed.

5.3. Simulation Verification of Electromechanical Actuator

The reference position θa of the actuator is set to 25° in the simulation. The multiphysical domain model and the original model are simultaneously simulated and compared to obtain the EMA motion speed Ωa (potential variable) and the power supply current Is (flow variable). The curves are shown in Figure 4. The dotted line is the simulation curve of the multiphysical domain model, and the solid line is the simulation curve of the original model.

Simulation results of EMA models. (a) Actuator speed. (b) Power supply current. (c) Power loss.

It can be seen from Figure 4 that the actuator speed simulation curve of the multiphysical domain model is basically consistent with the original model, while the power supply current and the power loss have some errors in the subtleties, which is due to the model linearization. The error effect can be ignored for the simulation goal proposed in this paper.

6. Integrated Method of Multiphysical Domain Model

Another advantage of the multiphysical domain model is that it facilitates large-scale integration to establish the complex electrical system. The multiphysical domain model integration is based on the energy conversion relationship, and there are two types of series integration and parallel integration.

6.1. Serial Integration of Multiphysical Domain Model

The multiphysical domain model of complex electrical equipment can be obtained by the serial integration of multiple simple multiphysical domain models.

For example, the EHA consists of the motor, the hydraulic pump, and the hydraulic actuator. If the EHA original model is available, the same EMA signal identification modeling method as equation (14) can be used. If there is no original model, it can also be modeled by the serial integration of the motor, the hydraulic pump, and the hydraulic actuator. In the integration, the consistency of the potential variable and the flow variable must be guaranteed. The structure is shown in Figure 5.

EHA model implemented by serial integration.

Since the EHA is controlled by the motor speed to realize the flow control of the hydraulic pump and the position control of the actuator, the control signal from the information conversion module is applied to the motor model. In addition, the EHA output is a linear motion; therefore, the potential variable of the output power according to Table 1 is the force Fa, and the flow variable is the moving speed va.

For example, the main parameters of the EHA are the power supply voltage 270VDC, the motor rated speed 8000r/min, the actuator rated moving speed 0.12 m/s, and the maximum force 200 kN. The EHA model is simulated with the position which is set to 80 mm. The simulation waveforms are shown in Figure 6. Among them, Figures 6(a)6(d) are the actuator output (mechanical energy), the hydraulic pump output (hydraulic energy), the motor output (mechanical energy), and the motor input (electrical energy). Figure 6(e) shows the power loss and the total loss of each EHA part.

Simulation waveforms of each EHA connection point. (a) Moving position, power, and variables of the rudder surface. (b) Hydraulic power and variables of the hydraulic pump. (c) Mechanical power and variables of the motor. (d) Electrical power and variables of the motor. (e) Power loss of EHA.

It can be seen from the simulation results that the EHA model obtained by the series integration of the multiphysical domain model can not only simulate the EHA input and output variables and power data but also observe the waveforms of each connection point in order to analyze the correctness of the model.

6.2. Parallel Integration of Multiphysical Domain Model

In the electrical system of the more electric aircraft, the power supply needs to provide electrical power for all the electrical equipment; that is, the electrical equipment forms the parallel structure on the bus bar. When all the electrical equipment is described by the multiphysical domain models, the structure is as shown in Figure 7.

Parallel integration of the power supply system.

The output power of the power supply is the electrical energy. The potential variable is the power supply voltage Us connected to the potential variable Usi of all electrical loads with the power supply side. The flow variable Is is the summation current of the electrical equipment input current, which is the power flow Isi summation of all electrical loads.(15)Us=Us1=Us2==UsmIs=Is1+Is2++Ism,

The flight control system of the aircraft is taken as an example. The flight control system consists of the left and right ailerons, left and right elevators, and the rudder. It is assumed to be driven by the EMA, established with the power supply system shown in Figure 7 and simulated.

It is assumed that the aircraft continuously completes the climb phase, the combat phase (including the yaw, roll, and dive flight attitude), and the descent phase of the flight. The control signals are obtained by the flight control system models. The corresponding reference signal θa is set for each EMA. The power characteristics of the simulation are shown in Figure 8. Figures 8(a)8(c) are the power characteristics of the elevator, the aileron, and the rudder, respectively. Figure 8(d) shows the power characteristics of the power supply. The movements of the left and right elevators and ailerons are symmetrical; therefore, the power characteristics are the same.

Flight characteristics simulation of the flight control system. (a) Elevator power characteristics. (b) Aileron power characteristics. (c) Rudder power characteristics. (d) Output power characteristics of the power supply.

It can be seen from the simulation results shown in Figure 8 that the elevator is the main drive device when the aircraft is climbing and descending. The elevator, aileron, and rudder need to be integrated operation when the aircraft is rolling and yawing, and the power conversion is complicated.

Since the electromechanical actuator belongs to the dynamic load of the short-time operation, the flight control system shows the pulse power of different sizes from the above figures. The power analysis for this type of system can only be done by the digital simulation.

7. Conclusion

The multiphysical domain model with the parameter identification according to the original model data is applied in this paper. It takes the power conversion process of the electrical equipment as the basic modeling and accurately expresses the power conversion and loss characteristics based on the dynamic characteristics of the device itself. The model based on the energy conversion relationship for system integration is to achieve the large-scale integration of the electrical system flexibly and the required digital simulation of the electrical system design.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

RoseroJ. A.OrtegaJ. A.AldabasE.RomeralL.Moving towards a more electric aircraftIEEE Aerospace and Electronic Systems Magazine20072233910.1109/maes.2007.3405002-s2.0-34047144949JonesR. I.The more electric aircraft-assessing the benefitsProceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering2002216525926910.1243/0954410023210287752-s2.0-0036926018WheelerP.Technology for the more and all electric aircraft of the futureProceedings of the International Conference on AutomaticOctober 2016Curico, ChileAveryC. R.BurrowS. G.MellorP. H.Electrical generation and distribution for the more electric aircraftProceedings of the 42th International Universalities Power Engineering ConferenceSeptember 2007Brighton, UK10071012JiY.KuhnM.Model based system level studies of more electric aircraftProceedings of Aerospace Systems and Technology ConferenceSeptember 2014Cincinnati, OH, USAShawJ. C.FletcherS. D. A.NormanP. J.GallowayS. J.More electric power system concepts for an environmentally responsible aircraftProceedings of the 47th International Universities Power Engineering ConferenceApril 2012London, UK16O’connellT.RussellG.MccarthyK.Energy management of an aircraft electrical systemProceedings of the 46th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & ExhibitJuly 2010Nashville, TN, USANairusJ.FingersR.More Electric Aircraft Power and Thermal Management, 2006Paris, FranceTechnologies for Energy Optimized Aircraft Equipment Systems ForumSchettniF.DentiE.RitoG. D.Development of a simulation platform of all electric aircraft on-board system for energy management studiesAeronautical Journal20171211239710719LawA. M.XiaoT. Y.Fan translatedW. H.Simulation Modeling and Analysis20124thBeijing, ChinaQinghua UniversityWaltersE.AmrheinM.O’connelT.INVENT modeling, simulation, analysis and optimizationProceedings of the 48th Aiaa Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace ExpositionJanuary 2010Orlando, FL, USA110ZhongJ.Coupling Design Theory and Method of Complex Electromechanical System20071stBeijing, ChinaMechanical pressGeoffR.NickK.SteveB.Bond graph modeling of a hydraulic vibration system: simulation and controlProceedings of the Integrated Modeling and Analysis in Applied Control and AutomationSeptember 2011Rome, Italy281286MutukuM.JohnK.Bond graph modeling of inter-actuator in teractions in a multi-cylinder hydraulicsystem,” World Academy of ScienceEngineering and Technology2011743241LiB.LiW.ZhangX.Modeling and simulation of aircraft power supply system based on Dymola and ModelicaProceedings of the Electrical Systems for Aircraft, Railway, Ship Propulsion and Road Vehicles & International Transportation Electrification ConferenceNovember 2017Toulouse, FranceBalsJ.JiY.KuhnM.Model based design and integration of more electric aircraft systems using modelicaProceedings of the European Power Electronics Conference and ExhibitionSeptember 2009Europe, Spain113LiuL.CaoZ.SunL. R.Energy optimization characteristic analysis of electromechanical actuator on more electric aircraftProceedings of the Communications in Computer and Information ScienceJune 2016Beijing, China4152QiaoG.LiuG.ShiZ.WangY.MaS.LimT. C.A review of electromechanical actuators for More/All Electric aircraft systemsJournal of Mechanical Engineering Science2018232224128415110.1177/09544062177498692-s2.0-85045470490RitoG. D.GalatoloR.SchettiniF.Experimental and simulation study of the dynamics of an electro-mechanical landing gear actuatorProceedings of the 30th International Council of the Aeronautical SciencesSeptember 2016Daejeon, South KoreaHuW. N.ZhouL.TianY. S.Analysis for the power loss of electrohydrostatic actuator and hydraulic actuatorProceedings of the Advanced Intelligent MechatronicsJuly 2015Busan, South Korea613616LiY.JiZ.YangL.ZhangP.XuB.ZhangJ.Thermal-fluid-structure coupling analysis for valve plate friction pair of axial piston pump in electrohydrostatic actuator (EHA) of aircraftApplied Mathematical Modelling20174783985810.1016/j.apm.2016.08.0152-s2.0-85008953604